ABSTRACT
We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By ...comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher mass and dynamically perturbed systems. Similar trends are also seen for the azimuthal temperature scatter and non-thermal pressure fraction, both of which correlate with density fluctuations, with these values also increasing with redshift. However, in agreement with recent work, we find only a weak correlation between the clumping, or its proxies, and the hydrostatic mass bias. To reduce the effect of clumping in the projected profiles, we compute the azimuthal median following recent observational studies, and find this reduces the scatter in the bias. We also attempt to correct the cluster masses by using a non-thermal pressure term and find overcorrected mass estimates (1 − b = 0.86 to 1 − b = 1.15) from 3D gas profiles but improved mass estimates (1 − b = 0.75 to 1 − b = 0.85) from projected gas profiles, with the caveat of systematically increased scatter. We conclude that the cluster-averaged mass bias is minimized from applying a non-thermal pressure correction (1 − b = 0.85) with more modest reductions from selecting clusters that have low clumping (1 − b = 0.79) or are dynamically relaxed (1 − b = 0.80). However, the latter selection is most effective at minimizing the scatter for individual objects. Such results can be tested with next-generation X-ray missions equipped with high-resolution spectrometers such as Athena.
We use the BAHAMAS (BAryons and HAloes of MAssive Systems) and MACSIS (MAssive ClusterS and Intercluster Structures) hydrodynamic simulations to quantify the impact of baryons on the mass ...distribution and dynamics of massive galaxy clusters, as well as the bias in X-ray and weak lensing mass estimates. These simulations use the subgrid physics models calibrated in the BAHAMAS project, which include feedback from both supernovae and active galactic nuclei. They form a cluster population covering almost two orders of magnitude in mass, with more than 3500 clusters with masses greater than 10 super( 14) M... at z = 0. We start by characterizing the clusters in terms of their spin, shape and density profile, before considering the bias in both weak lensing and hydrostatic mass estimates. Whilst including baryonic effects leads to more spherical, centrally concentrated clusters, the median weak lensing mass bias is unaffected by the presence of baryons. In both the dark matter only and hydrodynamic simulations, the weak lensing measurements underestimate cluster masses by ...10 per cent for clusters with M200 = 10 super( 15) M... and this bias tends to zero at higher masses. We also consider the hydrostatic bias when using both the true density and temperature profiles, and those derived from X-ray spectroscopy. When using spectroscopic temperatures and densities, the hydrostatic bias decreases as a function of mass, leading to a bias of ...40 per cent for clusters with M500 greater than or equal to 10 super( 15) M... This is due to the presence of cooler gas in the cluster outskirts. Using mass weighted temperatures and the true density profile reduces this bias to 5-15 per cent. (ProQuest: ... denotes formulae/symbols omitted.)
Abstract
We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, ...a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range 14.0 ≤ log10(M200 c/M⊙) ≤ 15.4, with their sophisticated subgrid physics modelling and high numerical resolution (subkpc gravitational softening), making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5–10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased (<5 per cent) estimator of the velocity dispersion. This result holds out to z = 1.5 and is relatively insensitive to the choice of cluster aperture, varying by less than 5 per cent between r500 c and r200 m. We show that the velocity bias is a function of the time spent by a galaxy inside the cluster environment. Selecting galaxies by their total mass results in a larger bias because a larger fraction of objects have only recently entered the cluster and these have a velocity bias above unity. Galaxies that entered more than 4 Gyr ago become progressively colder with time, as expected from dynamical friction. We conclude that velocity bias should not be a major issue when estimating cluster masses from kinematic methods.
We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching ...a wide range of observed properties. MACSIS extends the recent BAryons and HAloes of MAssive Systems simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of the effect of selecting only the hottest systems. We study the observable-mass scaling relations and the X-ray luminosity-temperature relation over the complete observed cluster mass range. As expected, we find that the slope of these scaling relations and the evolution of their normalization with redshift depart significantly from the self-similar predictions. However, for a sample of hot clusters with core-excised temperatures k sub( B)T greater than or equal to 5 keV, the normalization and the slope of the observable-mass relations and their evolution are significantly closer to self-similar. The exception is the temperature-mass relation, for which the increased importance of non-thermal pressure support and biased X-ray temperatures leads to a greater departure from self-similarity in the hottest systems. As a consequence, these also affect the slope and evolution of the normalization in the luminosity-temperature relation. The median hot gas profiles show good agreement with observational data at z = 0 and z = 1, with their evolution again departing significantly from the self-similar prediction. However, selecting a hot sample of clusters yields profiles that evolve significantly closer to the self-similar prediction. In conclusion, our results show that understanding the selection function is vital for robust calibration of cluster properties with mass and redshift.
We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the Wilkinson Microwave Anisotropy Probe year 5 ...(WMAP5) cosmology. The median relation between concentration and mass is adequately described by a power law for halo masses in the range 1011–1015 h−1 M⊙ and redshifts z < 2, regardless of whether the halo density profiles are fitted using Navarro, Frenk & White or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of σ8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 1011 h−1 M⊙ to 16 per cent at 1014 h−1 M⊙. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters.
Abstract
We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M
200c = 1014–1015.4 M⊙) with baryon particle masses of ∼106 ...M⊙. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the ‘Cluster–EAGLE’ project. They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses M
star > 1010 M⊙, this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (M
star ≲ 1010 M⊙) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (>5 r
200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ∼0.3 dex even well beyond r
200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated.
ABSTRACT
The Sunyaev–Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain ${\gt }10\, 000$ simulated galaxy groups and ...clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (tSZ) effect. The third measure determines the lowest order relativistic correction to the tSZ signal, which is seeing increased observational relevance. Our procedure allows us to accurately model the relativistic SZ (rSZ) contribution and we show that a ${\gtrsim}10\!-\!40{{\ \rm per\ cent}}$ underestimation of this rSZ cluster temperature is expected when applying standard X-ray relations. The correction also exhibits significant mass and redshift evolution, as we demonstrate here. We present the mass dependence of each temperature measure alongside their profiles and a short analysis of the temperature dispersion as derived from the aforementioned simulations. We also discuss a new relation connecting the temperature and Compton-y parameter, which can be directly used for rSZ modelling. Simple fits to the obtained scaling relations and profiles are provided. These should be useful for future studies of the rSZ effect and its relevance to cluster cosmology.
ABSTRACT
The kinetic Sunyaev–Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster ...motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulations. To utilize the rotational kSZ as a cosmological probe, simulations offer some of the most comprehensive data sets that can inform the modelling of this signal. In this work, we use the MACSIS data set to investigate the rotational kSZ effect in massive clusters specifically. Based on these models, we test stacking approaches and estimate the amplitude of the combined signal with varying mass, dynamical state, redshift, and map-alignment geometry. We find that the dark matter, galaxy and gas spins are generally misaligned, an effect that can cause a suboptimal estimation of the rotational kSZ effect when based on galaxy motions. Furthermore, we provide halo-spin–mass scaling relations that can be used to build a statistical model of the rotational kSZ. The rotational kSZ contribution, which is largest in massive unrelaxed clusters (≳100 $\mu$K), could be relevant to studies of higher order CMB temperature signals, such as the moving lens effect. The limited mass range of the MACSIS sample strongly motivates an extended investigation of the rotational kSZ effect in large-volume simulations to refine the modelling, particularly towards lower mass and higher redshift, and provide forecasts for upcoming cosmological CMB experiments (e.g. Simons Observatory, SKA-2) and X-ray observations (e.g. Athena/X-IFU).
The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, ...we draw on a large suite of high-resolution cosmological hydrodynamical simulations to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The effects of baryons on the dark matter distribution are typically not well described by adiabatic contraction models. In the inner 10 per cent of the virial radius the models are only successful if we allow their parameters to vary with baryonic physics, halo mass and redshift, thereby removing all predictive power. On larger scales the profiles from dark matter only simulations consistently provide better fits than adiabatic contraction models, even when we allow the parameters of the latter models to vary. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anticorrelation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that active galactic nuclei (AGN) feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disc galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons. We therefore conclude that the disagreement between the concentrations inferred from observations of groups of galaxies and predictions from simulations that was identified by Duffy et al. is not alleviated by the inclusion of baryons.