August E. Evrard

Simulations of the formation, evolution and clustering of galaxies and quasars

The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.Numerical simulations are a primary theoretical tool to study the nonlinear gravitational growth of structure in the Universe, and to link the initial conditions of cold dark matter (CDM) cosmogonies to observations of galaxies at the present day. Without direct numerical simulation, the hierarchical build-up of structure with its threedimensional dynamics would be largely inaccessible. Since the dominant mass component, the dark matter, is assumed to consist of weakly interacting elementary particles that interact only gravitationally, such simulations use a set of discrete point particles to represent the collisionless dark matter fluid. This representation as an N-body system is obviously only a coarse approximation, and im-

Mass estimates of X-ray clusters

We use cosmological gasdynamic simulations to investigate the accuracy of galaxy cluster mass estimates based on X-ray observations. The experiments follow the formation of clusters in different cosmological models and include the effects of gravity, pressure gradients, and hydrodynamical shocks. A subset of our ensemble also allows for feedback of mass and energy from galactic winds into the intracluster medium. We find that mass estimates based on the hydrostatic, isothermal {beta}-model are remarkably accurate when evaluated at radii where the cluster mean density is between 500 and 2500 times the critical density. At lower densities, radial temperature information becomes important. In the quoted radial regime, the distribution of the estimated-to-true mass ratio, derived from 174 artificial images constructed from the simulations, is nearly unbiased and has a standard deviation of 14{percent}{endash}29{percent}. The scatter can be considerably reduced (to 8{percent}{endash}15{percent}) by using an alternative mass estimator that exploits the tightness of the mass-temperature relation found in the simulations. The improvement over {beta}-model estimates is due to the elimination of the variance contributed by the gas outer slope parameter. We discuss these findings and their implications for recent measurements of cluster baryon fractions. {copyright} {ital 1996 The American Astronomical Society.}

Properties of the Intracluster Medium in an Ensemble of Nearby Galaxy Clusters

We present a systematic analysis of the intracluster medium (ICM) in an X-ray flux limited sample of 45 galaxy clusters. Using archival ROSAT Position-Sensitive Proportional Counter (PSPC) data and published ICM temperatures, we present best-fit double and single β model profiles, and extract ICM central densities and radial distributions. We use the data and an ensemble of numerical cluster simulations to quantify sources of uncertainty for all reported parameters. We examine the ensemble properties within the context of models of structure formation and feedback from galactic winds. We present best-fit ICM mass-temperature MICM-TX relations for MICM calculated within r500 and 1 h-150 Mpc. These relations exhibit small scatter (17%), providing evidence of regularity in large, X-ray flux limited cluster ensembles. Interestingly, the slope of the MICM-TX relation (at limiting radius r500) is steeper than the self-similar expectation by 4.3 σ. We show that there is a mild dependence of ICM mass fraction fICM on TX; the clusters with ICM temperatures below 5 keV have a mean ICM mass fraction fICM=0.160±0.008, which is significantly lower than that of the hotter clusters fICM=0.212 ± 0.006 (90% confidence intervals). In apparent contradiction with previously published analyses, our large, X-ray flux limited cluster sample provides no evidence for a more extended radial ICM distribution in low-TX clusters down to the sample limit of 2.4 keV. By analyzing simulated clusters we find that density variations enhance the cluster X-ray emission and cause MICM and fICM to be overestimated by ~12%. Additionally, we use the simulations to estimate an fICM depletion factor at r500. We use the bias corrected mean fICM within the hotter cluster subsample as a lower limit on the cluster baryon fraction. In combination with nucleosynthesis constraints this measure provides a firm upper limit on the cosmological density parameter for clustered matter ΩM≤(0.36 ± 0.01) h-½50.

A MaxBCG Catalog of 13,823 Galaxy Clusters from the Sloan Digital Sky Survey

We present a catalog of galaxy clusters selected using the maxBCG red-sequence method from Sloan Digital Sky Survey photometric data. This catalog includes 13,823 clusters with velocity dispersions greater than 400 km s-1 and is the largest galaxy cluster catalog assembled to date. They are selected in an approximately volume-limited way from a 0.5 Gpc3 region covering 7500 deg2 of sky between redshifts 0.1 and 0.3. Each cluster contains between 10 and 190 E/S0 ridgeline galaxies brighter than 0.4L* within a scaled radius R200. The tight relation between ridgeline color and redshift provides an accurate photometric redshift estimate for every cluster. Photometric redshift errors are shown by comparison to spectroscopic redshifts to be small (Δ ≃ 0:01), essentially independent of redshift, and well determined throughout the redshift range. Runs of maxBCG on realistic mock catalogs suggest that the sample is more than 90% pure and more than 85% complete for clusters with masses ≥ 1 x 1014 M⊙. Spectroscopic measurements of cluster members are used to examine line-of-sight projection as a contaminant in the identification of brightest cluster galaxies and cluster member galaxies. Spectroscopic data are also used to demonstrate the correlation between optical richness and velocity dispersion. Comparison to the combined NORAS and REFLEX X-rayYselected cluster catalogs shows that X-rayYluminous clusters are found among the optically richer maxBCG clusters. This paper is the first in a series that will consider the properties of these clusters, their galaxy populations, and their implications for cosmology.

Cosmological Parameters from Observations of Galaxy Clusters

Studies of galaxy clusters have proved crucial in helping to establish the standard model of cosmology, with a Universe dominated by dark matter and dark energy. A theoretical basis that describes clusters as massive, multicomponent, quasi-equilibrium systems is growing in its capability to interpret multiwavelength observations of expanding scope and sensitivity. We review current cosmological results, including contributions to fundamental physics, obtained from observations of galaxy clusters. These results are consistent with and complementary to those from other methods. We highlight several areas of opportunity for the next few years, and emphasize the need for accurate modeling of survey selection and sources of systematic error. Capitalizing on these opportunities will require a multiwavelength approach and the application of rigorous statistical frameworks, utilizing the combined strengths of observers, simulators, and theorists.

The LX − T Relation and Intracluster Gas Fractions of X-ray Clusters

We re-examine the X-ray luminosity–temperature relation using a nearly homogeneous data set of 24 clusters selected for statistically accurate temperature measurements and absence of strong cooling flows. The data exhibit a remarkably tight power-law relation between bolometric luminosity and temperature with a slope 2.88 ± 0.15. With reasonable assumptions regarding cluster structure, we infer an upper limit on fractional variations in the intracluster gas fraction 〈(δ fgas fgas)2〉1/2≤ 15 per cent. A strictly homogeneous Ginga subset of 18 clusters places a more stringent limit of 9 per cent. Imaging data from the literature are employed to determine absolute values of fgas within spheres encompassing density contrasts δc=500 and 200 with respect to the critical density. Comparing binding mass estimates based on the virial theorem (VT) and the hydrostatic β-model (BM), we find a temperature-dependent discrepancy in fgas between the two methods caused by systematic variation of the outer slope parameter β with temperature. Mean values (for H0=50 km s−1 Mpc−1) range from f¯gas=0.10 for cool (T 4 keV) clusters using the BM at δc=200. There is evidence that cool clusters have a lower mean gas fraction than hot clusters, but it is not possible to assess the statistical significance of this effect in the present data set. The T dependence of the intracluster medium (ICM) density structure, coupled with the increase of the gas fraction with T in the VT approach, explains the steepening of the LX–T relation. The small variation about the mean gas fraction within this majority subpopulation of clusters presents an important constraint for theories of galaxy formation and supports arguments against an Einstein–de Sitter universe based on the population mean gas fraction and conventional, primordial nucleosynthesis. The apparent trend of lower gas fractions and more extended atmospheres in low-temperature systems is consistent with expectations of models incorporating the effects of galactic winds on the ICM.

The Santa Barbara Cluster Comparison Project: A Comparison of Cosmological Hydrodynamics Solutions

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including smoothed particle hydrodynamics (SPH) and grid methods with fixed, deformable, or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gasdynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be nonradiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying, although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of 2. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to within about 10%, and the ratio of specific dark matter kinetic to gas thermal energies agree to within about 5%. Various factors, including differences in the internal timing of the simulations, contribute to the spread in calculated cluster properties. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modeled.

Cosmological Constraints from the SDSS maxBCG Cluster Catalog

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness-mass relation of the clusters. Assuming a flat {Lambda}CDM cosmology, we find {sigma}{sub 8}({Omega}{sub m}/0.25){sup 0.41} = 0.832 {+-} 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness-mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find {sigma}{sub 8} = 0.807 {+-} 0.020 and {Omega}{sub m} = 0.265 {+-} 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

Galaxy clusters in Hubble volume simulations: Cosmological constraints from sky survey populations

We use gigaparticle N-body simulations to study galaxy cluster populations in Hubble volumes of ΛCDM and ΓCDM world models.

Expectations for X-ray cluster observations by the ROSAT satellite

Predictions for characteristics of the X-ray cluster population expected to be observed by the Rosat satellite are presented via theoretical models for the evolution of X-ray cluster abundances. With a square detect cell of side 4.8 arcmin employed to search for extended cluster sources, roughly 200 clusters per steradian above the all-sky survey limiting flux of 7 x 10 exp -13 ergs/s sq cm are expected. For the north ecliptic pole region of the Rosat survey, which is to cover a 10-deg radius field to a limiting flux of roughly 9 x 10 exp -14 ergs/s sq cm, about 330 cluster X-ray sources are predicted. The clusters should have a median redshift of 0.2, and 10 percent of them should have z greater than 0.4.