ABSTRACT In this paper, we investigate the level of hydrostatic equilibrium (HE) in the intracluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical ...simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which are stellar and active galactic nucleus feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces via the gas accelerations generated, we effectively examine the level of HE in every object of the sample and its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10%-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces the deviation from HE very well, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions ( ), relative to thermal ones. As a general result, the HE-deviation and mass bias, at a given distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.
ABSTRACT We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core ...(CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and active galactic nucleus (AGN) feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of CC systems, to nearly flat core isentropic profiles, characteristic of NCC systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and observations. Furthermore, we also find that simulated CC clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic CC structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.
We present an analysis of the properties of the intracluster medium (ICM) in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the treepm+sph
...gadget-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae (SNe), one of which also accounts for the effect of feedback from active galactic nuclei (AGN) resulting from gas accretion on to supermassive black holes. These simulations are analysed with the aim of studying the relative role played by SN and AGN feedback on the general properties of the diffuse hot baryons in galaxy clusters and groups: scaling relations, temperature, entropy and pressure radial profiles, and ICM chemical enrichment. We find that simulations including AGN feedback produce scaling relations between X-ray observable quantities that are in good agreement with observations at all mass scales. Observed pressure profiles are also shown to be quite well reproduced in our radiative simulations, especially when AGN feedback is included. However, our simulations are not able to account for the observed diversity between cool-core and non-cool-core clusters, as revealed by X-ray observations: unlike for observations, we find that temperature and entropy profiles of relaxed and unrelaxed clusters are quite similar and resemble more the observed behaviour of non-cool-core clusters. As for the pattern of metal enrichment, we find that an enhanced level of iron abundance is produced by AGN feedback with respect to the case of purely SN feedback. As a result, while simulations including AGN produce values of iron abundance in groups in agreement with observations, they over-enrich the ICM in massive clusters. The efficiency of AGN feedback in displacing enriched gas from haloes into the intergalactic medium at high redshift also creates a widespread enrichment in the outskirts of clusters and produces profiles of iron abundance whose slope is in better agreement with observations. By analysing the pattern of the relative abundances of silicon and iron and the fraction of metals in the stellar phase, our results clearly show that different sources of energy feedback leave different imprints in the enrichment pattern of the hot ICM and stars. Our results confirm that including AGN feedback goes in the right direction of reconciling simulation predictions and observations for several observational ICM properties. Still a number of important discrepancies highlight that the model still needs to be improved to produce the correct interplay between cooling and feedback in central cluster regions.
Abstract
We analyse the radial pressure profiles, the intracluster medium (ICM) clumping factor and the Sunyaev–Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups ...identified in a set of hydrodynamical simulations based on an updated version of the treepm–SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among other processes, active galactic nucleus (AGN) and/or stellar feedback. Our results are analysed as a function of redshift, ICM physics, cluster mass and cluster cool-coreness or dynamical state. In general, the mean pressure profiles obtained for our sample of groups and clusters show a good agreement with X-ray and SZ observations. Simulated cool-core (CC) and non-cool-core (NCC) clusters also show a good match with real data. We obtain in all cases a small (if any) redshift evolution of the pressure profiles of massive clusters, at least back to z = 1. We find that the clumpiness of gas density and pressure increases with the distance from the cluster centre and with the dynamical activity. The inclusion of AGN feedback in our simulations generates values for the gas clumping ($\sqrt{C}_{\rho }\sim 1.2$ at R200) in good agreement with recent observational estimates. The simulated YSZ–M scaling relations are in good accordance with several observed samples, especially for massive clusters. As for the scatter of these relations, we obtain a clear dependence on the cluster dynamical state, whereas this distinction is not so evident when looking at the subsamples of CC and NCC clusters.
Abstract
We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, ...X-ray temperature, and YX. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics gadget-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with M500 > 1014 M⊙E(z)−1, mimicking the typical selection of Sunyaev–Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at z ∼ 2. The results of the analysis show a general agreement with observations. The values of the slope of the mass–gas mass and mass–temperature relations at z = 2 are 10 per cent lower with respect to z = 0 due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range z = 0–1, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and YX is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM.
Aims.
We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range 0 <
z
< 4, comparing them to recent observational ...studies; we also investigated the effect of varying the parameters of the star formation model on galaxy properties such as SFR, star-formation efficiency, and gas fraction.
Methods.
We analyse a set of zoom-in cosmological hydrodynamical simulations centred on 12 clusters. The simulations are carried out with the GADGET-3 Tree-PM smoothed-particle hydro-dynamics code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback.
Results.
Simulations do not reproduce the high values of SFR observed within protocluster cores, where the values of SFR are underpredicted by a factor ≳4 both at
z
∼ 2 and
z
∼ 4. The difference arises as simulations are unable to reproduce the observed starburst population and is greater at
z
∼ 2 because simulations underpredict the normalisation of the main sequence (MS) of star forming galaxies (i.e. the correlation between stellar mass and SFR) by a factor of ∼3. As the low normalisation of the MS seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to overly underpredicted gas fractions or overly low star formation efficiencies. Our results are stable against varying several parameters of the star formation subgrid model and do not depend on the details of AGN feedback.
Conclusions.
The subgrid model for star formation, introduced to reproduce the self-regulated evolution of quiescent galaxies, is not suitable to describe violent events like high-redshift starbursts. We find that this conclusion holds, independently of the parameter choice for the star formation and AGN models. The increasing number of multi-wavelength high-redshift observations will help to improve the current star formation model, which is needed to fully recover the observed star formation history of galaxy clusters.
Context. The correlations between the properties of the brightest cluster galaxy (BCG) and the mass of its central super-massive black hole (SMBH) have been extensively studied from a theoretical and ...observational angle. More recently, relations connecting the SMBH mass and global properties of the hosting cluster, such as temperature and mass, were observed. Aims. We investigate the correlation between SMBH mass and cluster mass and temperature, their establishment and evolution. We compare their scatter to that of the classical MBH − MBCG relation. Moreover, we study how gas accretion and BH-BH mergers contribute to SMBH growth across cosmic time. Methods. We employed 135 groups and clusters with a mass range 1.4 × 1013 M⊙ − 2.5 × 1015 M⊙ extracted from a set of 29 zoom-in cosmological hydro-dynamical simulations where the baryonic physics is treated with various sub-grid models, including feedback by active galactic nuclei. Results. In our simulations we find that MBH correlates well with M500 and T500, with the scatter around these relations compatible within 2σ with the scatter around MBH − MBCG at z = 0. The MBH − M500 relation evolves with time, becoming shallower at lower redshift as a direct consequence of hierarchical structure formation. On average, in our simulations the contribution of gas accretion to the total SMBH mass dominates for the majority of the cosmic time (z > 0.4), while in the last 2 Gyr the BH-BH mergers become a larger contributor. During this last process, substructures hosting SMBHs are disrupted in the merger process with the BCG and the unbound stars enrich the diffuse stellar component rather than increase BCG mass. Conclusions. From the results obtained in our simulations with simple sub-grid models we conclude that the scatter around the MBH − T500 relation is comparable to the scatter around the MBH − MBCG relation and that, given the observational difficulties related to the estimation of the BCG mass, clusters temperature and mass can be a useful proxy for the SMBHs mass, especially at high redshift.
ABSTRACT
Nowadays, Machine Learning techniques offer fast and efficient solutions for classification problems that would require intensive computational resources via traditional methods. We examine ...the use of a supervised Random Forest to classify stars in simulated galaxy clusters after subtracting the member galaxies. These dynamically different components are interpreted as the individual properties of the stars in the Brightest Cluster Galaxy (BCG) and IntraCluster Light (ICL). We employ matched stellar catalogues (built from the different dynamical properties of BCG and ICL) of 29 simulated clusters from the DIANOGA set to train and test the classifier. The input features are cluster mass, normalized particle cluster-centric distance, and rest-frame velocity. The model is found to correctly identify most of the stars, while the larger errors are exhibited at the BCG outskirts, where the differences between the physical properties of the two components are less obvious. We investigate the robustness of the classifier to numerical resolution, redshift dependence (up to z = 1), and included astrophysical models. We claim that our classifier provides consistent results in simulations for z < 1, at different resolution levels and with significantly different subgrid models. The phase-space structure is examined to assess whether the general properties of the stellar components are recovered: (i) the transition radius between BCG-dominated and ICL-dominated region is identified at 0.04 R200; (ii) the BCG outskirts (>0.1 R200) is significantly affected by uncertainties in the classification process. In conclusion, this work suggests the importance of employing Machine Learning to speed up a computationally expensive classification in simulations.
ABSTRACT
Contradictory results have been reported on the time evolution of the alignment between clusters and their brightest cluster galaxy (BCG). We study this topic by analysing cosmological ...hydrosimulations of 24 massive clusters with $M_{200}|_{z=0} \gtrsim 10^{15}\, \rm {\, M_{\odot }}$, plus 5 less massive with $1 \times 10^{14} \lesssim M_{200}|_{z=0} \lesssim 7 \times 10^{14}\, \rm {\, M_{\odot }}$, which have already proven to produce realistic BCG masses. We compute the BCG alignment with both the distribution of cluster galaxies and the dark matter (DM) halo. At redshift z = 0, the major axes of the simulated BCGs and their host cluster galaxy distributions are aligned on average within 20°. The BCG alignment with the DM halo is even tighter. The alignment persists up to z ≲ 2 with no evident evolution. This result continues, although with a weaker signal, when considering the projected alignment. The cluster alignment with the surrounding distribution of matter (3R200) is already in place at z ∼ 4 with a typical angle of 35°, before the BCG–cluster alignment develops. The BCG turns out to be also aligned with the same matter distribution, albeit always to a lesser extent. These results taken together might imply that the BCG–cluster alignment occurs in an outside–in fashion. Depending on their frequency and geometry, mergers can promote, destroy or weaken the alignments. Clusters that do not experience recent major mergers are typically more relaxed and aligned with their BCG. In turn, accretions closer to the cluster elongation axis tend to improve the alignment as opposed to accretions closer to the cluster minor axis.