COSMOS: A Hybrid N-Body/Hydrodynamics Code for Cosmological Problems

We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.Comment: 43 pages, 14 figures, submitted to ApJS and revised according to referee's comment

Off-Center Mergers of Clusters of Galaxies and Nonequipartition of Electrons and Ions in Intracluster Medium

We investigate the dynamical evolution of clusters of galaxies and their observational consequences during off-center mergers, explicitly considering the relaxation process between ions and electrons in intracluster medium by N-body and hydrodynamical simulations. In the contracting phase a bow shock is formed between the two subclusters. The observed temperature between two peaks in this phase depends on the viewing angle even if the geometry of the system seems to be very simple like head-on collisions. Around the most contracting epoch, when we observe merging clusters nearly along the collision axis, they look like spherical relaxed clusters with large temperature gradients. In the expanding phase, spiral bow shocks occur. As in head-on mergers, the electron temperature is significantly lower than the plasma mean one especially in the post-shock regions in the expanding phase. When the systems have relatively large angular momentum, double-peak structures in the X-ray images can survive even after the most contracting epoch. Morphological features in both X-ray images and electron temperature distribution characteristic to off-center mergers are seriously affected by the viewing angle. When the clusters are observed nearly along the collision axis, the distribution of galaxies' line-of-sight (LOS) velocities is a good indicator of mergers. In the contracting phase, an negative kurtosis and a large skewness are expected for nearly equal mass collisions and rather different mass ones, respectively. To obtain statistically significant results, about 1000 galaxies' LOS velocities are required. For nearby clusters ($z<0.05$), large redshift surveys such as 2dF will enable us to study merger dynamics.Comment: 21 pages, 7 figures. Accepted for publication in Ap

Constraining q_0 with Cluster Gas Mass Fractions: A Feasibility Study

As the largest gravitationally bound objects in the universe, clusters of galaxies may contain a fair sample of the baryonic mass fraction of the universe. Since the gas mass fraction from the hot ICM is believed to be constant in time, the value of the cosmological deceleration parameter $q_0$ can be determined by comparing the calculated gas mass fraction in nearby and distant clusters (Pen 1997). To test the potential of this method, we compare the gas fractions derived for a sample of luminous ($L_X > 10^{45}$erg s$^{-1}$), nearby clusters with those calculated for eight luminous, distant ($0.3 < z < 0.6$) clusters using ASCA and ROSAT observations. For consistency, we evaluate the gas mass fraction at a fixed physical radius of 1 $h_{50}^{-1}$ Mpc (assuming $q_0=0.0$). We find a best fit value of $q_0 = 0.07$ with -0.47 < q_0 < 0.67 at 95% confidence. We also determine the gas fraction using the method of Evrard, Metzler, & Navarro (1997) to find the total mass within $r_{500}$, the radius where the mean overdensity of matter is 500 times the critical density. In simulations, this method reduces the scatter in the determination of gravitational mass without biasing the mean. We find that it also reduces the scatter in actual observations for nearby clusters, but not as much as simulations suggest. Using this method, the best fit value is $q_0 = 0.04$ with -0.50 < q_0 < 0.64. The excellent agreement between these two methods suggests that this may be a useful technique for determining $q_0$. The constraints on $q_0$ should improve as more distant clusters are studied and precise temperature profiles are measured to large radii.Comment: 8 pages, 4 figures, uses emulateapj.sty, onecolfloat.st

The Ionized Gas Kinematics of the LMC-Type Galaxy NGC 1427A in the Fornax Cluster

NGC 1427A is a LMC-like irregular galaxy in the Fornax cluster with an extended pattern of strong star formation around one of its edges, which is probably due to some kind of interaction with the cluster environment. We present H-alpha velocities within NGC 1427A, obtained through long-slit spectroscopy at seven different positions, chosen to fall on the brightest HII regions of the galaxy. Due to its location very near the center of the cluster this object is an excellent candidate to study the effects that the cluster environment has on gas-rich galaxies embedded in it. The rotation of NGC 1427A is modeled in two different ways. The global ionized gas kinematics is reasonably well described by solid-body rotation, although on small scales it shows a chaotic behaviour. In this simple model, the collision with a smaller member of the cluster as being responsible for the peculiar morphology of NGC 1427A is very unlikely, since the only candidate intruder falls smoothly into the general velocity pattern of the main galaxy. In a more elaborate model, for which we obtain a better solution, this object does not lie in the same plane of NGC 1427A, in which case we identify it as a satellite bound to the galaxy. These results are discussed in the context of a normal irregular versus one interacting with some external agent. Based on several arguments and quantitative estimates, we argue that the passage through the hot intracluster gas of the Fornax cluster is a very likely scenario to explain the morphological properties of NGC 1427A.Comment: 31 pages, LaTeX2e, uses aas2pp4.sty and psfig.sty, including 7 Postscript figures; accepted for publication in ApJ, Vol. 530, February 200

Outer Regions of the Cluster Gaseous Atmospheres

We present a systematic study of the hot gas distribution in the outer regions of regular clusters using ROSAT PSPC data. Outside the cooling flow region, the beta-model describes the observed surface brightness closely, but not precisely. Between 0.3 and 1 virial radii, the profiles are characterized by a power law with slope, expressed in terms of the beta parameter, in the range beta=0.65 to 0.85. The values of beta in this range of radii are typically larger by ~0.05 than those derived from the global fit. There is a mild trend for the slope to increase with temperature, from ~0.68 for 3 keV clusters to ~0.8 for 10 keV clusters; however, even at high temperatures there are clusters with flat gas profiles, 0.7. Our values of beta at large radius are systematically higher, and the trend of beta with temperature is weaker than was previously found; the most likely explanation is that earlier studies were affected by an incomplete exclusion of the central cooling flow regions. For our regular clusters, the gas distribution at large radii is quite close to spherically symmetric and this is shown not to be an artifact of the sample selection. The gas density profiles are very similar when compared in the units of cluster virial radius. The radius of fixed mean gas overdensity 1000 (corresponding to the dark matter overdensity 200 for Omega=0.2) shows a tight correlation with temperature, R~T**0.5, as expected from the virial theorem for clusters with the universal gas fraction. At a given temperature, the rms scatter of the gas overdensity radius is only ~7% which translates into a 20% scatter of the gas mass fraction, including statistical scatter due to measurement uncertainties.Comment: ApJ in press, submitted 11/30/9

Radial Temperature Profiles of X-Ray--Emitting Gas Within Clusters of Galaxies

Previous analyses of ASCA data of clusters of galaxies have found conflicting results regarding the slope of the temperature profile of the hot X-ray gas within clusters, mainly because of the large, energy-dependent point spread function (PSF) of the ASCA mirrors. We present a summary of all ASCA-determined cluster temperature profiles found in the literature, and find a discrepancy in the radial temperature trend of clusters based on which PSF-correction routine is used. This uncertainty in the cluster temperature profile in turn can lead to large uncertainties in the amount of dark matter in clusters. In this study, we have used ROSAT PSPC data to obtain independent relative temperature profiles for 26 clusters, most of which have had their temperature profiles determined by ASCA. Our aim is not to measure the actual temperature values of the clusters, but to use X-ray color profiles to search for a hardening or softening of the spectra with radius for comparison to ASCA-derived profiles. The radial color profiles indicate that outside of the cooling flow region, the temperature profiles of clusters are in general constant. Within 35% of the virial radius, we find a temperature drop of 20% at 10 keV and 12% at 5 keV can be ruled out at the 99% confidence level. A subsample of non-cooling flow clusters shows that the condition of isothermality applies at very small radii too, although cooling gas complicates this determination in the cooling flow subsample. The colors predicted from the temperature profiles of a series of hydrodynamical cluster simulations match the data very well, although they cannot be used to discriminate among different cosmologies. An additional result is that the color profiles show evidence for a central peak in metallicity in low temperature clusters.Comment: 39 pages, 15 embedded Postscript figures, uses aaspp4.sty, accepted for publication in Astrophysical Journa

The Formation of a Realistic Disk Galaxy in Lambda Dominated Cosmologies

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)Comment: ApJ in press. Images and movies at http://hpcc.astro.washington.edu/faculty/fabio/galform.html Significantly expanded revised version. (9 pages vs the original 4