Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($\beta\sim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment.

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.

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 $Y_{X}$. 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 $M_{500} > 10^{14} M_{\odot} 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 \sim 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 $Y_X$ 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.

Comparisons between observed and predicted strong lensing properties of galaxy clusters have been routinely used to claim either tension or consistency with $\Lambda$CDM cosmology. However, standard approaches to such cosmological tests are unable to quantify the preference for one cosmology over another. We advocate approximating the relevant Bayes factor using a marginal likelihood that is based on the following summary statistic: the posterior probability distribution function for the parameters of the scaling relation between Einstein radii and cluster mass, $\alpha$ and $\beta$. We demonstrate, for the first time, a method of estimating the marginal likelihood using the X-ray selected $z>0.5$ MACS clusters as a case in point and employing both N-body and hydrodynamic simulations of clusters. We investigate the uncertainty in this estimate and consequential ability to compare competing cosmologies, that arises from incomplete descriptions of baryonic processes, discrepancies in cluster selection criteria, redshift distribution, and dynamical state. The relation between triaxial cluster masses at various overdensities provide a promising alternative to the strong lensing test.

The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes included in the simulation and with the thermal properties of the core. These simulations, consisting of re-simulations of 29 Lagrangian regions performed with an upgraded version of the SPH GADGET-3 code, have been run including two different sets of baryonic physics: one accounts for radiative cooling, star formation, metal enrichment and supernova (SN) feedback, and the other one further includes the effects of feedback from active galactic nuclei (AGN). In agreement with observations, we find an anti-correlation between entropy and metallicity in cluster cores, and similar radial distributions of heavy-element abundances and abundance ratios out to large cluster-centric distances (~R180). In the outskirts, namely outside of ~0.2R180, we find a remarkably homogeneous metallicity distribution, with almost flat profiles of the elements produced by either SNIa or SNII. We investigated the origin of this phenomenon and discovered that it is due to the widespread displacement of metal-rich gas by early (z>2-3) AGN powerful bursts, acting on small high-redshift haloes. Our results also indicate that the intrinsic metallicity of the hot gas for this sample is on average consistent with no evolution between z=2 and z=0, across the entire radial range.

We analyze the radial pressure profiles, the 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, AGN and/or stellar feedback. Our results are analyzed 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 center 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 $R_{200}$) in good agreement with recent observational estimates. The simulated $Y_{SZ}-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.

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 sample 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.

We present an analysis of the properties of the 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, one of which also accounts for the effect of feedback from AGN resulting from gas accretion onto super-massive 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 that are in good agreement with X-ray observations at all mass scales. However, our simulations are not able to account for the observed diversity between CC and NCC clusters: 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 NCC 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 halos into the inter-galactic 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.

Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study the ICM of inhomogeneities, focusing on the ones on the large scale that, unlike clumps, are the most difficult to identify. To this purpose we introduce the concept of residual clumpiness, C_R, that quantifies the large-scale inhomogeneity of the ICM. After showing that this quantity can be robustly defined for relaxed systems, we characterize how it varies with radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces an overestimate in the determination of the density profile from the X-ray emission, which translates into a systematic overestimate of 6 (12)% in the measurement of M_gas at R_200 for our relaxed (perturbed) cluster sample. At the same time, the increase of C_R with radius introduces also a ~2% systematic underestimate in the measurement of the hydrostatic-equilibrium mass (M_he), which adds to the previous one generating a systematic ~8.5% overestimate in f_gas in our relaxed sample. Since the residual clumpiness of the ICM is not directly observable, we study its correlation with the azimuthal scatter in the X-ray surface brightness of the halo and in the y-parameter profiles. We find that their correlation is highly significant (r_S = 0.6-0.7), allowing to define the azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as robust proxies of C_R. After providing a function that connects the two quantities, we obtain that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias in the measurement of M_gas for relaxed objects, which becomes (0+/-2)% up to 2R_200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of M_he and f_gas, although a significant halo to halo scatter remains. (abridged)

We carry out an analysis of a set of cosmological SPH hydrodynamical simulations of galaxy clusters and groups aimed at studying the total baryon budget in clusters, and how this budget is shared between the hot diffuse component and the stellar component. Using the TreePM+SPH GADGET-3 code, we carried out one set of non-radiative simulations, and two sets of simulations including radiative cooling, star formation and feedback from supernovae (SN), one of which also accounting for the effect of feedback from active galactic nuclei (AGN). The analysis is carried out with the twofold aim of studying the implication of stellar and hot gas content on the relative role played by SN and AGN feedback, and to calibrate the cluster baryon fraction and its evolution as a cosmological tool. We find that both radiative simulation sets predict a trend of stellar mass fraction with cluster mass that tends to be weaker than the observed one. However this tension depends on the particular set of observational data considered. Including the effect of AGN feedback alleviates this tension on the stellar mass and predicts values of the hot gas mass fraction and total baryon fraction to be in closer agreement with observational results. We further compute the ratio between the cluster baryon content and the cosmic baryon fraction, Y_b, as a function of cluster-centric radius and redshift. At R_500 we find for massive clusters with M_500>2\times10^{14} h^{-1} M_sun that Y_b is nearly independent of the physical processes included and characterized by a negligible redshift evolution: Y_{b,500}=0.85+/-0.03 with the error accounting for the intrinsic r.m.s. scatter within the set of simulated clusters. At smaller radii, R_2500, the typical value of Y_b slightly decreases, by an amount that depends on the physics included in the simulations, while its scatter increases by about a factor of two.

[Abridged] We present an analysis of the relation between the masses of cluster- and group-sized halos, extracted from $\Lambda$CDM cosmological N-body and hydrodynamic simulations, and their velocity dispersions, at different redshifts from $z=2$ to $z=0$. The main aim of this analysis is to understand how the implementation of baryonic physics in simulations affects such relation, i.e. to what extent the use of the velocity dispersion as a proxy for cluster mass determination is hampered by the imperfect knowledge of the baryonic physics. In our analysis we use several sets of simulations with different physics implemented. Velocity dispersions are determined using three different tracers, DM particles, subhalos, and galaxies. We confirm that DM particles trace a relation that is fully consistent with the theoretical expectations based on the virial theorem and with previous results presented in the literature. On the other hand, subhalos and galaxies trace steeper relations, and with larger values of the normalization. Such relations imply that galaxies and subhalos have a $\sim10$ per cent velocity bias relative to the DM particles, which can be either positive or negative, depending on halo mass, redshift and physics implemented in the simulation. We explain these differences as due to dynamical processes, namely dynamical friction and tidal disruption, acting on substructures and galaxies, but not on DM particles. These processes appear to be more or less effective, depending on the halo masses and the importance of baryon cooling, and may create a non-trivial dependence of the velocity bias and the $\soneD$--$\Mtwo$ relation on the tracer, the halo mass and its redshift. These results are relevant in view of the application of velocity dispersion as a proxy for cluster masses in ongoing and future large redshift surveys.

The observed abundance of giant arcs produced by galaxy cluster lenses and the measured Einstein radii have presented a source of tension for LCDM. Previous cosmological tests for high-redshift clusters (z>0.5) have suffered from small number statistics in the simulated sample and the implementation of baryonic physics is likely to affect the outcome. We analyse zoomed-in simulations of a fairly large sample of cluster-sized objects, with Mvir > 3x10^14 Msun/h, identified at z=0.25 and z=0.5, for a concordance LCDM cosmology. We start with dark matter only simulations, and then add gas hydrodynamics, with different treatments of baryonic processes: non-radiative cooling, radiative cooling with star formation and galactic winds powered by supernova explosions, and finally including the effect of AGN feedback. We find that the addition of gas in non-radiative simulations does not change the strong lensing predictions significantly, but gas cooling and star formation together significantly increase the number of expected giant arcs and the Einstein radii, particularly for lower redshift clusters and lower source redshifts. Further inclusion of AGN feedback reduces the predicted strong lensing efficiencies such that the lensing probability distributions becomes closer to those obtained for simulations including only dark matter. Our results indicate that the inclusion of baryonic physics in simulations will not solve the arc-statistics problem at low redshifts, when the physical processes included provide a realistic description of cooling in the central regions of galaxy clusters. [Abridged]

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

We introduce a generalized scaling law, M_tot = 10^K A^a B^b, to look for the minimum scatter in reconstructing the total mass of hydrodynamically simulated X-ray galaxy clusters, given gas mass M_gas, luminosity L and temperature T. We find a locus in the plane of the logarithmic slopes $a$ and $b$ of the scaling relations where the scatter in mass is minimized. This locus corresponds to b_M = -3/2 a_M +3/2 and b_L = -2 a_L +3/2 for A=M_gas and L, respectively, and B=T. Along these axes, all the known scaling relations can be identified (at different levels of scatter), plus a new one defined as M_tot ~ (LT)^(1/2). Simple formula to evaluate the expected evolution with redshift in the self-similar scenario are provided. In this scenario, no evolution of the scaling relations is predicted for the cases (b_M=0, a_M=1) and (b_L=7/2, a_L=-1), respectively. Once the single quantities are normalized to the average values of the sample under considerations, the normalizations K corresponding to the region with minimum scatter are very close to zero. The combination of these relations allows to reduce the number of free parameters of the fitting function that relates X-ray observables to the total mass and includes the self-similar redshift evolution.

We present a detailed study of scaling relations between total cluster mass and three mass proxies based on X-ray observables: temperature of the intra-cluster medium, gas mass and the product of the two, Y_X. Our analysis is based on two sets of high-resolution hydrodynamical simulations performed with the TreePM-SPH GADGET code. The first set includes about 140 clusters with masses above 5x10^13 M_sun/h (30 having mass above 10^15 M_sun/h), that have been simulated with (i) non-radiative physics and including (ii) cooling, star formation, chemical enrichment and the effect of supernova feedback triggering galactic ejecta. This large statistics is used to quantify the robustness of the scaling relations, to determine their redshift evolution and to calibrate their intrinsic scatter and its distribution. We use a smaller set of clusters including 18 halos with masses above 5x10^13 M_sun/h to test the robustness of mass proxies against changing the physical processes included in simulations (thermal conduction, artificial viscosity, cooling and star formation, galactic winds and AGN feedback). We find the M-Y_X scaling relation to be the least sensitive one to variations of the ICM physics, with its slope and redshift evolution close to the self-similar model predictions. The distribution of the scatter around the best-fitting relations is close to a log-normal one. M_gas has the smallest scatter in mass, with values of sigma_lnM = 0.04-0.06, depending on the physics included in the simulation, and with a mild dependence on redshift. The M-T relation is the one with the largest scatter, with sigma_lnM > 0.1 at z=0, increasing to > 0.15 at z=1. The intrinsic scatter in the M-Y_X relation is slightly larger than in the M-M_gas relation. These results confirm that both Y_X and M_gas mass proxies are well suited for cosmological applications of future large X-ray surveys. [abridged]

We present a study of the effect of AGN feedback on metal enrichment and thermal properties of the intracluster medium (ICM) in hydrodynamical simulations. The cosmological simulations are performed for a set of clusters using a version of the TreePM-SPH Gadget code that follows chemo-dynamical evolution by accounting for metal enrichment by different stellar populations. Besides runs not including any efficient form of energy feedback, we carry out simulations including: (i) kinetic feedback in the form of galactic winds triggered by supernova explosions; (ii) AGN feedback from gas accretion onto super-massive black holes (BHs); (iii) AGN feedback in which a 'radio mode' is included. We find that AGN feedback is able to quench star formation in the brightest cluster galaxies at z<4 and provides correct temperature profiles in the central regions of galaxy groups. However, its effect is not sufficient to create cool cores in massive clusters. AGN feedback creates a widespread enrichment in the outskirts of clusters, thanks to its efficiency in displacing enriched gas from galactic halos to the inter-galactic medium at relatively high redshift. Iron abundance profiles are in better agreement with observations, with a more pristine enrichment of the ICM around and beyond the cluster virial regions. From the pattern of the relative abundances of Silicon and Iron, we conclude that a significant fraction of ICM enrichment in simulations is contributed by a diffuse population of intra-cluster stars. Our simulations also predict that profiles of Z_Si/Z_Fe abundance ratio do not increase at least out to 0.5 R_vir. Our results clearly show that different sources of energy feedback leave distinct imprints in the enrichment pattern of the ICM, that are more evident when looking at cluster external regions.

We present a comparison between simulation results and X-ray observational data on the evolution of the metallicity of the intra-cluster medium (ICM). The simulations of galaxy clusters were performed with the Tree-SPH Gadget2 code that includes a detailed model of chemical evolution, by assuming three different shapes for the stellar initial mass function (IMF), namely the Salpeter (1955), Kroupa (2001) and Arimoto-Yoshii (1987) IMF. Our simulations predict significant radial gradients of the Iron abundance, which extend over the whole cluster virialized region. At larger radii, we do not detect any flattening of the metallicity profiles. As for the evolution of the ICM metal (Iron) abundance out to z=1, we find that it is determined by the combined action of (i) the sinking of already enriched gas, (ii) the ongoing metal production in galaxies and (iii) the locking of ICM metals in newborn stars. As a result, rather than suppressing the metallicity evolution, stopping star formation at z=1 has the effect of producing an even too fast evolution of the emission-weighted ICM metallicity with too high values at low redshift. Finally, we compare simulations with the observed rate of type-Ia supernovae per unit B-band luminosity (SnU_B). We find that our simulated clusters do not reproduce the decreasing trend of SnU_B at low redshift, unless star formation is truncated at z=1.

The study of the metal enrichment of the intra-cluster and inter-galactic media (ICM and IGM) represents a direct means to reconstruct the past history of star formation, the role of feedback processes and the gas-dynamical processes which determine the evolution of the cosmic baryons. In this paper we review the approaches that have been followed so far to model the enrichment of the ICM in a cosmological context. While our presentation will be focused on the role played by hydrodynamical simulations, we will also discuss other approaches based on semi-analytical models of galaxy formation, also critically discussing pros and cons of the different methods. We will first review the concept of the model of chemical evolution to be implemented in any chemo-dynamical description. We will emphasise how the predictions of this model critically depend on the choice of the stellar initial mass function, on the stellar life-times and on the stellar yields. We will then overview the comparisons presented so far between X-ray observations of the ICM enrichment and model predictions. We will show how the most recent chemo-dynamical models are able to capture the basic features of the observed metal content of the ICM and its evolution. We will conclude by highlighting the open questions in this study and the direction of improvements for cosmological chemo-dynamical models of the next generation.