In the hierarchical scenario of structure formation, galaxy clusters are the ultimate virialised products in mass and time. Hot baryons in the intracluster medium (ICM) and cold baryons in galaxies inhabit a dark matter dominated halo. Internal proce
sses, accretion, and mergers can perturb the equilibrium, which is established only at later times. However, the cosmic time when thermalisation is effective is still to be assessed. Here we show that massive clusters in the observed universe attained an advanced thermal equilibrium $sim~1.8~text{Gyr}$ ago, at redshift $z =0.14pm0.06$, when the universe was $11.7pm0.7~text{Gyr}$ old. Hot gas is mostly thermalised after the time when cosmic densities of matter and dark energy match. We find in a statistically nearly complete and homogeneous sample of 120 clusters from the {it Planck} Early Sunyaev-Zeldovich (ESZ) sample that the kinetic energy traced by the galaxy velocity dispersion is a faithful probe of the gravitational energy since a look back time of at least $sim5.4~text{Gyr}$, whereas the efficiency of hot gas in converting kinetic to thermal energy, as measured through X-ray observations in the core-excised area within $r_{500}$, steadily increases with time. The evolution is detected at the $sim 98$ per cent probability level. Our results demonstrate that halo mass accretion history plays a larger role for cluster thermal equilibrium than radiative physics. The evolution of hot gas is strictly connected to the cosmic structure formation.
Recent observations reveal that, at a given stellar mass, blue galaxies tend to live in haloes with lower mass while red galaxies live in more massive host haloes. The physical driver behind this is still unclear because theoretical models predict th
at, at the same halo mass, galaxies with high stellar masses tend to live in early-formed haloes which naively leads to an opposite trend. Here, we show that the {sc Simba} simulation quantitatively reproduces the colour bimodality in SHMR and reveals an inverse relationship between halo formation time and galaxy transition time. It suggests that the origin of this bimodality is rooted in the intrinsic variations of the cold gas content due to halo assembly bias. {sc Simba}s SHMR bimodality quantitatively relies on two aspects of its input physics: (1) Jet-mode AGN feedback, which quenches galaxies and sets the qualitative trend; and (2) X-ray AGN feedback, which fully quenches galaxies and yields better agreement with observations. The interplay between the growth of cold gas and the AGN quenching in {sc Simba} results in the observed SHMR bimodality.
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 distributi
ons 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.
Recent numerical studies of the dark matter density profiles of massive galaxy clusters ($M_{rm halo} > 10^{15}$M$_{odot}$) show that their median radial mass density profile remains unchanged up to $z > 1$, displaying a highly self-similar evolution
. We verify this by using the data set of the THE THREE HUNDRED project, i.e. 324 cluster-sized haloes as found in full physics hydrodynamical simulations. We track the progenitors of the mass-complete sample of clusters at $z=0$, and find that their median shape is already in place by $z=2.5$. However, selecting a dynamically relaxed subsample ($sim16$ per cent of the clusters), we observe a shift of the scale radius $r_s$ towards larger values at earlier times. Classifying the whole sample by formation time, this evolution is understood as a result of a two-phase halo mass accretion process. Early-forming clusters -- identified as relaxed today -- have already entered their slow accretion phase, hence their mass growth occurs mostly at the outskirts. Late-forming clusters -- which are still unrelaxed today -- are in their fast accretion phase, thus the central region of the clusters is still growing. We conclude that the density profile of galaxy clusters shows a profound self-similarity out to redshifts $zsim2.5$. This result holds for both gas and total density profiles when including baryonic physics, as reported here for two rather distinct sub-grid models.
Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical samp
le of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.
[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of
the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component.
