No Arabic abstract
Accurate and precise measurements of masses of galaxy clusters are key to derive robust constraints on cosmological parameters. Rising evidence from observations, however, confirms that X-ray masses, obtained under the assumption of hydrostatic equilibrium, might be underestimated, as previously predicted by cosmological simulations. We analyse more than 300 simulated massive clusters, from `The Three Hundred Project, and investigate the connection between mass bias and several diagnostics extracted from synthetic X-ray images of these simulated clusters. We find that the azimuthal scatter measured in 12 sectors of the X-ray flux maps is a statistically significant indication of the presence of an intrinsic (i.e. 3D) clumpy gas distribution. We verify that a robust correction to the hydrostatic mass bias can be inferred when estimates of the gas inhomogeneity from X-ray maps (such as the azimuthal scatter or the gas ellipticity) are combined with the asymptotic external slope of the gas density or pressure profiles, which can be respectively derived from X-ray and millimetric (Sunyaev-Zeldovich effect) observations. We also obtain that mass measurements based on either gas density and temperature or gas density and pressure result in similar distributions of the mass bias. In both cases, we provide corrections that help reduce both the dispersion and skewness of the mass bias distribution. These are effective even when irregular clusters are included leading to interesting implications for the modelling and correction of hydrostatic mass bias in cosmological analyses of current and future X-ray and SZ cluster surveys.
Galaxy clusters have a triaxial matter distribution. The weak-lensing signal, an important part in cosmological studies, measures the projected mass of all matter along the line-of-sight, and therefore changes with the orientation of the cluster. Studies suggest that the shape of the brightest cluster galaxy (BCG) in the centre of the cluster traces the underlying halo shape, enabling a method to account for projection effects. We use 324 simulated clusters at four redshifts between 0.1 and 0.6 from `The Three Hundred Project to quantify correlations between the orientation and shape of the BCG and the halo. We find that haloes and their embedded BCGs are aligned, with an average $sim$20 degree angle between their major axes. The bias in weak lensing cluster mass estimates correlates with the orientation of both the halo and the BCG. Mimicking observations, we compute the projected shape of the BCG, as a measure of the BCG orientation, and find that it is most strongly correlated to the weak-lensing mass for relaxed clusters. We also test a 2-dimensional cluster relaxation proxy measured from BCG mass isocontours. The concentration of stellar mass in the projected BCG core compared to the total stellar mass provides an alternative proxy for the BCG orientation. We find that the concentration does not correlate to the weak-lensing mass bias, but does correlate with the true halo mass. These results indicate that the BCG shape and orientation for large samples of relaxed clusters can provide information to improve weak-lensing mass estimates.
Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the textsc{Mock-X} analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-$y$ and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias $bapprox0.13$ with cluster mass, independent of hydrodynamical method or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to $b=0.3$ for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias.
We carry out a systematic study of the recently discovered fundamental plane of galaxy clusters (CFP) using a sample of ~250 simulated clusters from the 300th project, focusing on the stability of the plane against different temperature definitions and its dependence on the dynamical relaxation state of clusters. The CFP is characterised in the form of $T propto M_s^alpha r_s^beta$, defined with the gas temperature ($T$) and the characteristic halo scale radius and mass ($r_s$ and $M_s$) assuming an NFW halo description. We explore two definitions of weighted temperatures, namely mass-weighted and spectroscopic-like temperatures, in three radial ranges: [0.1, 1.0]$r_{200}$, [0.15,1.0]$r_{500}$, and [50,500]$h^{-1}$ kpc. We find that 300th clusters at $z=0$ lie on a thin plane whose parameters ($alpha, beta$) and dispersion (0.015--0.030 dex) depend on the gas temperature definition. The CFP for mass-weighted temperatures is closer to the virial equilibrium expectation ($alpha=1, beta=-1$) with a smaller dispersion. When gas temperatures are measured inside 500$h^{-1}$ kpc, which is close to the median value of $r_s$, the resulting CFP deviates the most from the virial expectation and shifts towards the similarity solution for a secondary infall model ($alpha=1.5, beta=-2$). Independently of the temperature definition, we find that clusters at $z=1$ form a CFP similar to the virial expectation. At all epochs, the CFP remains well defined throughout the evolution of the cluster population. The CFP of relaxed clusters is always close to the virial expectation, with a milder evolution than for the unrelaxed case. We find that only systems formed over the last 4 Gyr have a CFP that is closer to the self-similar solution. All these findings are compatible with the CFP obtained for a CLASH subsample excluding the hottest clusters with $T_X>12$ keV.
Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.
We analyse the gas content evolution of infalling haloes in cluster environments from THE THREE HUNDRED project, a collection of 324 numerically modelled galaxy clusters. The haloes in our sample were selected within $5R_{200}$ of the main cluster halo at $z=0$ and have total halo mass $M_{200}geq10^{11} h^{-1} M_{odot}$. We track their main progenitors and study their gas evolution since their crossing into the infall region, which we define as $1-4R_{200}$. Studying the radial trends of our populations using both the full phase space information and a line-of-sight projection, we confirm the Arthur et al. (2019) result and identify a characteristic radius around $1.7R_{200}$ in 3D and at $R_{200}$ in projection at which infalling haloes lose nearly all of the gas prior their infall. Splitting the trends by subhalo status we show that subhaloes residing in group-mass and low-mass host haloes in the infall region follow similar radial gas-loss trends as their hosts, whereas subhaloes of cluster-mass host haloes are stripped of their gas much further out. Our results show that infalling objects suffer significant gaseous disruption that correlates with time-since-infall, cluster-centric distance and host mass, and that the gaseous disruption they experience is a combination of subhalo pre-processing and object gas depletion at a radius which behaves like an accretion shock.