No Arabic abstract
We review recent progress in the description of the formation and evolution of galaxy clusters in a cosmological context by using numerical simulations. We focus our presentation on the comparison between simulated and observed X-ray properties, while we will also discuss numerical predictions on properties of the galaxy population in clusters. Many of the salient observed properties of clusters, such as X-ray scaling relations, radial profiles of entropy and density of the intracluster gas, and radial distribution of galaxies are reproduced quite well. In particular, the outer regions of cluster at radii beyond about 10 per cent of the virial radius are quite regular and exhibit scaling with mass remarkably close to that expected in the simplest case in which only the action of gravity determines the evolution of the intra-cluster gas. However, simulations generally fail at reproducing the observed cool-core structure of clusters: simulated clusters generally exhibit a significant excess of gas cooling in their central regions, which causes an overestimate of the star formation and incorrect temperature and entropy profiles. The total baryon fraction in clusters is below the mean universal value, by an amount which depends on the cluster-centric distance and the physics included in the simulations, with interesting tensions between observed stellar and gas fractions in clusters and predictions of simulations. Besides their important implications for the cosmological application of clusters, these puzzles also point towards the important role played by additional physical processes, beyond those already included in the simulations. We review the role played by these processes, along with the difficulty for their implementation, and discuss the outlook for the future progress in numerical modeling of clusters.
Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $sigma_mathrm{v}$ with the density $rho$ can be used to define the pseudo-entropy $S(r)=sigma_mathrm{v}^2/rho^{,2/3}$, whose profile is well-described by a simple power-law $Spropto,r^{,alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profiles as traced by substructures, while DM particles profiles are characterized by an almost universal behavior, thus suggesting that the level of pseudo-entropy could represent a potential low-scatter mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$times10^{11}$ M$_odot$ to 2$times10^{12}$ M$_odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.
Studies of cluster mass and velocity anisotropy profiles are useful tests of dark matter models, and of the assembly history of clusters of galaxies. These studies might be affected by unknown systematics caused by projection effects. We aim at testing observational methods for the determination of mass and velocity anisotropy profiles of clusters of galaxies. Particularly, we focus on the MAMPOSSt technique (Mamon et al. 2013). We use results from two semi-analytic models of galaxy formation coupled with high-resolution N-body cosmological simulations, the catalog of De Lucia & Blaizot (2007) and the FIRE catalog based on the new GAlaxy Evolution and Assembly model. We test the reliability of the Jeans equation in recovering the true mass profile when full projected phase-space information is available. We examine the reliability of the MAMPOSSt method in estimating the true mass and velocity anisotropy profiles of the simulated halos when only projected phase-space information is available, as in observations. The spherical Jeans equation provides a reliable tool for the determination of cluster mass profiles, also for subsamples of tracers separated by galaxy color. Results are equally good for prolate and oblate clusters. Using only projected phase-space information, MAMPOSSt provides estimates of the mass profile with a standard deviation of 35-69 %, and a negative bias of 7-17 %, nearly independent of radius, and that we attribute to the presence of interlopers in the projected samples. The bias changes sign, that is, the mass is over-estimated, for prolate clusters with their major axis aligned along the line-of-sight. MAMPOSSt measures the velocity anisotropy profiles accurately in the inner cluster regions, with a slight overestimate in the outer regions, both for the whole sample of observationally-identified cluster members and separately for red and blue galaxies.
We present results from a new set of 30 cosmological simulations of galaxy clusters, including the effects of radiative cooling, star formation, supernova feedback, black hole growth and AGN feedback. We first demonstrate that our AGN model is capable of reproducing the observed cluster pressure profile at redshift, z~0, once the AGN heating temperature of the targeted particles is made to scale with the final virial temperature of the halo. This allows the ejected gas to reach larger radii in higher-mass clusters than would be possible had a fixed heating temperature been used. Such a model also successfully reduces the star formation rate in brightest cluster galaxies and broadly reproduces a number of other observational properties at low redshift, including baryon, gas and star fractions; entropy profiles outside the core; and the X-ray luminosity-mass relation. Our results are consistent with the notion that the excess entropy is generated via selective removal of the densest material through radiative cooling; supernova and AGN feedback largely serve as regulation mechanisms, moving heated gas out of galaxies and away from cluster cores. However, our simulations fail to address a number of serious issues; for example, they are incapable of reproducing the shape and diversity of the observed entropy profiles within the core region. We also show that the stellar and black hole masses are sensitive to numerical resolution, particularly the gravitational softening length; a smaller value leads to more efficient black hole growth at early times and a smaller central galaxy.
The uniformity of the intra-cluster medium (ICM) enrichment level in the outskirts of nearby galaxy clusters suggests that chemical elements were deposited and widely spread into the intergalactic medium before the cluster formation. This observational evidence is supported by numerical findings from cosmological hydrodynamical simulations, as presented in Biffi et al. (2017), including the effect of thermal feedback from active galactic nuclei. Here, we further investigate this picture, by tracing back in time the spatial origin and metallicity evolution of the gas residing at z=0 in the outskirts of simulated galaxy clusters. In these regions, we find a large distribution of iron abundances, including a component of highly-enriched gas, already present at z=2. At z>1, the gas in the present-day outskirts was distributed over tens of virial radii from the the main cluster and had been already enriched within high-redshift haloes. At z=2, about 40% of the most Fe-rich gas at z=0 was not residing in any halo more massive than 1e11 Msun/h in the region and yet its average iron abundance was already 0.4, w.r.t. the solar value by Anders & Grevesse (1989). This confirms that the in situ enrichment of the ICM in the outskirts of present-day clusters does not play a significant role, and its uniform metal abundance is rather the consequence of the accretion of both low-metallicity and pre-enriched (at z>2) gas, from the diffuse component and through merging substructures. These findings do not depend on the mass of the cluster nor on its core properties.