Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $\sigma_\mathrm{v}$ with the density $\rho$ can be used to define the pseudo-entropy $S(r)=\sigma_\mathrm{v}^2/\rho^{2/3}$, whose profile is well-described by a simple power-law $S\propto r^{\alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profile as traced by substructures, while pseudo-entropy profiles as traced by DM particles are characterized by an almost universal behavior and can be used as a low-scatter ($\sim6-7\%$) mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.

We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

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. We analyze a set of zoom-in cosmological hydrodynamical simulations centred on twelve clusters. The simulations are carried out with the GADGET-3 TreePM/SPH code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback. Simulations do not reproduce the high values of SFR observed within protoclusters cores, where the values of SFR are underpredicted by a factor $\gtrsim 4$ both at $z\sim2$ and $z\sim 4$. The difference arises as simulations are unable to reproduce the observed starburst population and is worsened at $z\sim 2$ because simulations underpredict the normalization of the main sequence of star forming galaxies (i.e., the correlation between stellar mass and SFR) by a factor of $\sim 3$. As the low normalization of the main sequence seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to a too low predicted gas fractions or too 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 the AGN feedback.

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 analyzing cosmological hydro-simulations of 24 massive clusters with $M_{200}|_{z=0} \gtrsim 10^{15}\, M_\odot$, plus 5 less massive with $1 \times 10^{14} \lesssim M_{200}|_{z=0} \lesssim 7 \times 10^{14}\, 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$^\circ$. The BCG alignment with the DM halo is even tighter. The alignment persists up to $z\lesssim2$ 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 ($3R_{200}$) is already in place at $z\sim4$ with a typical angle of $35^\circ$, 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.

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.

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.

Recently, relations connecting the SMBH mass of central galaxies and global properties of the hosting cluster, such as temperature and mass, were observed. 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 $M_{\rm BH}-M_{\rm BCG}$ relation. We study how gas accretion and BH-BH mergers contribute to SMBH growth across cosmic time. We employed 135 groups and clusters with a mass range $1.4\times 10^{13}M_{\odot}-2.5\times 10^{15} M_{\odot}$ 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 AGN. In our simulations we find that $M_{\rm BH}$ correlates well with $M_{500}$ and $T_{500}$, with the scatter around these relations compatible within $2\sigma$ with the scatter around $M_{\rm BH}-M_{\rm BCG}$ at $z=0$. The $M_{\rm BH}-M_{500}$ 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. From the results obtained in our simulations with simple sub-grid models we conclude that the scatter around the $M_{\rm BH}-T_{500}$ relation is comparable to the scatter around the $M_{\rm BH}-M_{\rm BCG}$ 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.

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-\sigma_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

The possibility of determining the value of the Hubble constant using observations of galaxy clusters in X-ray and microwave wavelengths through the Sunyaev Zel\'dovich (SZ) effect has long been known. Previous measurements have been plagued by relatively large errors in the observational data and severe biases induced, for example, by cluster triaxiality and clumpiness. The advent of \textit{Planck} allows us to map the Compton parameter y, that is, the amplitude of the SZ effect, with unprecedented accuracy at intermediate cluster-centric radii, which in turn allows performing a detailed spatially resolved comparison with X-ray measurements. Given such higher quality observational data, we developed a Bayesian approach that combines informed priors on the physics of the intracluster medium obtained from hydrodynamical simulations of massive clusters with measurement uncertainties. We apply our method to a sample of 61 galaxy clusters with redshifts up to z < 0.5 observed with Planck and XMM-Newton observations and find H_0=67 \pm 3 km s^{-1} Mpc^{-1}.

We use TheThreeHundred project, a suite of 324 resimulated massive galaxy clusters embedded in a broad range of environments, to investigate (i) how the gas content of surrounding haloes correlates with phase-space position at $z=0$, and (ii) to investigate the role that ram pressure plays in this correlation. By stacking all 324 normalised phase-space planes containing 169287 haloes and subhaloes, we show that the halo gas content is tightly correlated with phase-space position. At $\sim\,1.5-2\,\text{R}_{\text{200}}$ of the cluster dark matter halo, we find an extremely steep decline in the halo gas content of infalling haloes and subhaloes irrespective of cluster mass, possibly indicating the presence of an accretion shock. We also find that subhaloes are particularly gas-poor, even in the cluster outskirts, which could indicate active regions of ongoing pre-processing. By modelling the instantaneous ram pressure experienced by each halo and subhalo at $z=0$, we show that the ram pressure intensity is also well correlated with phase-space position, which is again irrespective of cluster mass. In fact, we show that regions in the phase-space plane with high differential velocity between a halo or subhalo and its local gas environment, are almost mutually exclusive with high halo gas content regions. This suggests a causal link between the gas content of objects and the instantaneous ram pressure they experience, where the dominant factor is the differential velocity.

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.

According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_\Delta\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses.

The hot plasma in galaxy clusters is expected to be heated to high temperatures through shocks and adiabatic compression. The thermodynamical properties of the gas encode information on the processes leading to the thermalization of the gas in the cluster's potential well as well as non-gravitational processes such as gas cooling, AGN feedback and kinetic energy. In this work we present the radial profiles of the thermodynamic properties of the intracluster medium (ICM) out to the virial radius for a sample of 12 galaxy clusters selected from the Planck all-sky survey. We determine the universal profiles of gas density, temperature, pressure, and entropy over more than two decades in radius. We exploit jointly X-ray information from XMM and Sunyaev-Zel'dovich constraints from Planck to recover thermodynamic properties out to 2 R500. We provide average functional forms for the radial dependence of the main quantities and quantify the slope and intrinsic scatter of the population as a function of radius. We find that gas density and pressure profiles steepen steadily with radius, in excellent agreement with previous observational results. Entropy profiles beyond R500 closely follow the predictions for the gravitational collapse of structures. The scatter in all thermodynamical quantities reaches a minimum in the range [0.2-0.8] R500 and increases outwards. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region R > 0.3 R500) follow remarkably well the predictions of simple gravitational collapse and require little non-gravitational corrections.

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons.

Answers to the metal production of the Universe can be found in galaxy clusters, notably within their Intra-Cluster Medium (ICM). The X-ray Integral Field Unit (X-IFU) on board the next-generation European X-ray observatory Athena (2030s) will provide the necessary leap forward in spatially-resolved spectroscopy required to disentangle the intricate mechanisms responsible for this chemical enrichment. In this paper, we investigate the future capabilities of the X-IFU in probing the hot gas within galaxy clusters. From a test sample of four clusters extracted from cosmological hydrodynamical simulations, we present comprehensive synthetic observations of these clusters at different redshifts (up to z = 2) and within the scaled radius R500 performed using the instrument simulator SIXTE. Through 100 ks exposures, we demonstrate that the X-IFU will provide spatially-resolved mapping of the ICM physical properties with little to no biases (<5%) and well within statistical uncertainties. The detailed study of abundance profiles and abundance ratios within R500 also highlights the power of the X-IFU in providing constraints on the various enrichment models. From synthetic observations out to z = 2, we also quantify its ability to track the chemical elements across cosmic time with excellent accuracy, and thereby to investigate the evolution of metal production mechanisms as well as the link to the stellar initial mass-function. Our study demonstrates the unprecedented capabilities of the X-IFU in unveiling the properties of the ICM but also stresses the data analysis challenges faced by future high-resolution X-ray missions such as Athena.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

Observations and numerical simulations have shown that the relation between the mass scaled with the critical density of the universe and the X-ray temperature of galaxy clusters is approximately represented by $M_\Delta \propto T_X^{3/2}$ (e.g. $\Delta=500$). This relation is often interpreted as evidence that clusters are in virial equilibrium. However, the recently discovered fundamental plane (FP) of clusters indicates that the temperature of clusters primarily depends on a combination of the characteristic mass $M_s$ and radius $r_s$ of the Navarro-Frenk-White profile rather than $M_\Delta$. Moreover, the angle of the FP revealed that clusters are not in virial equilibrium because of continuous mass accretion from the surrounding matter. By considering both the FP and the mass dependence of the cluster concentration parameter, we show that this paradox can be solved and the relation $M_\Delta \propto T_X^{3/2}$ actually reflects the central structure of clusters. We also find that the intrinsic scatter in the halo concentration-mass relation can largely account for the spread of clusters on the FP. We also show that X-ray data alone form the FP and the angle and the position are consistent with those of the FP constructed from gravitational lensing data. We demonstrate that a possible shift between the two FPs can be used to calibrate cluster masses obtained via X-ray observations.

The X-ray Integral Field Unit (X-IFU) that will be on board the Athena telescope will provide an unprecedented view of the intracluster medium (ICM) kinematics through the observation of gas velocity, $v$, and velocity dispersion, $w$, via centroid-shift and broadening of emission lines, respectively. The improvement of data quality and quantity requires an assessment of the systematics associated with this new data analysis, namely biases, statistical and systematic errors, and possible correlations between the different measured quantities. We have developed an end-to-end X-IFU simulator that mimics a full X-ray spectral fitting analysis on a set of mock event lists, obtained using SIXTE. We have applied it to three hydrodynamical simulations of a Coma-like cluster that include the injection of turbulence. This allowed us to assess the ability of X-IFU to map five physical quantities in the cluster core: emission measure, temperature, metal abundance, velocity and velocity dispersion. Finally, starting from our measurements maps, we computed the 2D structure function (SF) of emission measure fluctuations, $v$ and $w$ and compared them with those derived directly from the simulations. All quantities match with the input projected values without bias; the systematic errors were below 5%, except for velocity dispersion whose error reaches about 15%. Moreover, all measurements prove to be statistically independent, indicating the robustness of the fitting method. Most importantly, we recover the slope of the SFs in the inertial regime with excellent accuracy, but we observe a systematic excess in the normalization of both SF$_v$ and SF$_w$ ascribed to the simplistic assumption of uniform and (bi-)Gaussian measurement errors. Our work highlights the excellent capabilities of Athena X-IFU in probing the thermodynamic and kinematic properties of the ICM. (abridged)

We analyze the stellar growth of Brightest Cluster Galaxies (BCGs) produced by cosmological zoom-in hydrodynamical simulations of the formation of massive galaxy clusters. The evolution of the stellar mass content is studied considering different apertures, and tracking backwards either the main progenitor of the $z=0$ BCG or that of the cluster hosting the BCG at $z=0$. Both methods lead to similar results up to $z \simeq 1.5$. The simulated BCGs masses at $z=0$ are in agreement with recent observations. In the redshift interval from $z=1$ to $z=0$ we find growth factors 1.3, 1.6 and 3.6 for stellar masses within 30kpc, 50kpc and 10% of $R_{500}$ respectively. The first two factors, and in general the mass evolution in this redshift range, are in agreement with most recent observations. The last larger factor is similar to the growth factor obtained by a semi-analytical model (SAM). Half of the star particles that end up in the inner 50 kpc was typically formed by redshift $\sim$ 3.7, while the assembly of half of the BCGs stellar mass occurs on average at lower redshifts $\sim 1.5$. This assembly redshift correlates with the mass attained by the cluster at high $z \gtrsim 1.3$, due to the broader range of the progenitor clusters at high-$z$. The assembly redshift of BCGs decreases with increasing apertures. Our results are compatible with the {\it inside-out} scenario. Simulated BCGs could lack intense enough star formation (SF) at high redshift, while possibly exhibit an excess of residual SF at low redshift.

In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane.

The determination of the morphology of galaxy clusters has important repercussion on their cosmological and astrophysical studies. In this paper we address the morphological characterisation of synthetic maps of the Sunyaev--Zel'dovich (SZ) effect produced for a sample of 258 massive clusters ($M_{vir}>5\times10^{14}h^{-1}$M$_\odot$ at $z=0$), extracted from the MUSIC hydrodynamical simulations. Specifically, we apply five known morphological parameters, already used in X-ray, two newly introduced ones, and we combine them together in a single parameter. We analyse two sets of simulations obtained with different prescriptions of the gas physics (non radiative and with cooling, star formation and stellar feedback) at four redshifts between 0.43 and 0.82. For each parameter we test its stability and efficiency to discriminate the true cluster dynamical state, measured by theoretical indicators. The combined parameter discriminates more efficiently relaxed and disturbed clusters. This parameter had a mild correlation with the hydrostatic mass ($\sim 0.3$) and a strong correlation ($\sim 0.8$) with the offset between the SZ centroid and the cluster centre of mass. The latter quantity results as the most accessible and efficient indicator of the dynamical state for SZ studies.

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.

We jointly analyze Bolocam Sunyaev-Zeldovich (SZ) effect and Chandra X-ray data for a set of 45 clusters to derive gas density and temperature profiles without using spectroscopic information. The sample spans the mass and redshift range $3 \times 10^{14} M_{\odot} \le M_{500} \le 25 \times 10^{14} M_{\odot}$ and $0.15\le z \le 0.89$. We define cool-core (CC) and non-cool core (NCC) subsamples based on the central X-ray luminosity, and 17/45 clusters are classified as CC. In general, the profiles derived from our analysis are found to be in good agreement with previous analyses, and profile constraints beyond $r_{500}$ are obtained for 34/45 clusters. In approximately 30% of the CC clusters our analysis shows a central temperature drop with a statistical significance of $>3\sigma$; this modest detection fraction is due mainly to a combination of coarse angular resolution and modest S/N in the SZ data. Most clusters are consistent with an isothermal profile at the largest radii near $r_{500}$, although 9/45 show a significant temperature decrease with increasing radius. The sample mean density profile is in good agreement with previous studies, and shows a minimum intrinsic scatter of approximately 10% near $0.5 \times r_{500}$. The sample mean temperature profile is consistent with isothermal, and has an intrinsic scatter of approximately 50% independent of radius. This scatter is significantly higher compared to earlier X-ray-only studies, which find intrinsic scatters near 10%, likely due to a combination of unaccounted for non-idealities in the SZ noise, projection effects, and sample selection.

We derive and compare the fractions of cool-core clusters in the {\em Planck} Early Sunyaev-Zel'dovich sample of 164 clusters with $z \leq 0.35$ and in a flux-limited X-ray sample of 100 clusters with $z \leq 0.30$, using {\em Chandra} observations. We use four metrics to identify cool-core clusters: 1) the concentration parameter: the ratio of the integrated emissivity profile within 0.15 $r_{500}$ to that within $r_{500}$, and 2) the ratio of the integrated emissivity profile within 40 kpc to that within 400 kpc, 3) the cuspiness of the gas density profile: the negative of the logarithmic derivative of the gas density with respect to the radius, measured at 0.04 $r_{500}$, and 4) the central gas density, measured at 0.01 $r_{500}$. We find that the sample of X-ray selected clusters, as characterized by each of these metrics, contains a significantly larger fraction of cool-core clusters compared to the sample of SZ selected clusters (44$\pm$7\% vs. 28$\pm$4\% using the concentration parameter in the 0.15--1.0 $r_{500}$ range, 61$\pm$8\% vs. 36$\pm$5\% using the concentration parameter in the 40--400 kpc range, 64$\pm$8\% vs. 38$\pm$5\% using the cuspiness, and 53$\pm$7\% vs. 39$\pm$5\% using the central gas density). Qualitatively, cool-core clusters are more X-ray luminous at fixed mass. Hence, our X-ray flux-limited sample, compared to the approximately mass-limited SZ sample, is over-represented with cool-core clusters. We describe a simple quantitative model that uses the excess luminosity of cool-core clusters compared to non-cool-core clusters at fixed mass to successfully predict the observed fraction of cool-core clusters in X-ray selected samples.

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.

We present a study on the coherent rotation of the intracluster medium and dark matter components of simulated galaxy clusters extracted from a volume-limited sample of the MUSIC project. The set is re-simulated with three different recipes for the gas physics: $(i)$ non-radiative, $(ii)$ radiative without AGN feedback, and $(iii)$ radiative with AGN feedback. Our analysis is based on the 146 most massive clusters identified as relaxed, 57 per cent of the total sample. We classify these objects as rotating and non-rotating according to the gas spin parameter, a quantity that can be related to cluster observations. We find that 4 per cent of the relaxed sample is rotating according to our criterion. By looking at the radial profiles of their specific angular momentum vector, we find that the solid body model is not a suitable description of rotational motions. The radial profiles of the velocity of the dark matter show a prevalence of the random velocity dispersion. Instead, the intracluster medium profiles are characterized by a comparable contribution from the tangential velocity and the dispersion. In general, the dark matter component dominates the dynamics of the clusters, as suggested by the correlation between its angular momentum and the gas one, and by the lack of relevant differences among the three sets of simulations.

In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster 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 stellar and AGN 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, 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 very well the deviation from HE, 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 (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.

The Sunyaev-Zeldovich (SZ) effect is a promising tool for detecting the presence of hot gas out to the galaxy cluster peripheries. We developed a spectral imaging algorithm dedicated to the SZ observations of nearby galaxy clusters with Planck, with the aim of revealing gas density anisotropies related to the filamentary accretion of materials, or pressure discontinuities induced by the propagation of shock fronts. To optimize an unavoidable trade-off between angular resolution and precision of the SZ flux measurements, the algorithm performs a multiscale analysis of the SZ maps as well as of other extended components, such as the cosmic microwave background (CMB) anisotropies and the Galactic thermal dust. The demixing of the SZ signal is tackled through kernel weighted likelihood maximizations. The CMB anisotropies are further analyzed through a wavelet analysis, while the Galactic foregrounds and SZ maps are analyzed via a curvelet analysis that best preserves their anisotropic details. The algorithm performance has been tested against mock observations of galaxy clusters obtained by simulating the Planck High Frequency Instrument and by pointing a few characteristic positions in the sky. These tests suggest that Planck should easily allow us to detect filaments in the cluster peripheries and detect large-scale shocks in colliding galaxy clusters that feature favorable geometry.

Morphology is often used to infer the state of relaxation of galaxy clusters. The regularity, symmetry, and degree to which a cluster is centrally concentrated inform quantitative measures of cluster morphology. The Cluster Lensing and Supernova survey with Hubble Space Telescope (CLASH) used weak and strong lensing to measure the distribution of matter within a sample of 25 clusters, 20 of which were deemed to be relaxed based on their X-ray morphology and alignment of the X-ray emission with the BCG. Towards a quantitative characterization of this important sample of clusters, we present uniformly estimated X-ray morphological statistics for all 25 CLASH clusters. We compare X-ray morphologies of CLASH clusters with those identically measured for a large sample of simulated clusters from the MUSIC-2 simulations, selected by mass. We confirm a threshold in X-ray surface brightness concentration of C>0.4 for cool-core clusters, where C is the ratio of X-ray emission inside 100 kpc/h70 compared to inside 500 kpc/h70. We report and compare morphologies of these clusters inferred from Sunyaev-Zeldovich Effect (SZE) maps of the hot gas and in from projected mass maps based on strong and weak lensing. We find a strong agreement in alignments of the orientation of major axes for the lensing, X-ray, and SZE maps of nearly all of the CLASH clusters at radii of 500 kpc (approximately 0.5R500 for these clusters). We also find a striking alignment of clusters shapes at the 500 kpc scale, as measured with X-ray, SZE, and lensing, with that of the near-infrared stellar light at 10 kpc scales for the 20 "relaxed" clusters. This strong alignment indicates a powerful coupling between the cluster- and galaxy-scale galaxy formation processes.

By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $\Omega_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $\Omega_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.

We present an initial study of the mass and evolutionary state of a massive and distant cluster, RCS2 J232727.6-020437. This cluster, at z=0.6986, is the richest cluster discovered in the RCS2 project. The mass measurements presented in this paper are derived from all possible mass proxies: X-ray measurements, weak-lensing shear, strong lensing, Sunyaev Zel'dovich effect decrement, the velocity distribution of cluster member galaxies, and galaxy richness. While each of these observables probe the mass of the cluster at a different radius, they all indicate that RCS2 J232727.6-020437 is among the most massive clusters at this redshift, with an estimated mass of M_200 ~3 x10^15 h^-1 Msun. In this paper, we demonstrate that the various observables are all reasonably consistent with each other to within their uncertainties. RCS2 J232727.6-020437 appears to be well relaxed -- with circular and concentric X-ray isophotes, with a cool core, and no indication of significant substructure in extensive galaxy velocity data.

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 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 on 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 cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core 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 cool-core 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.

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle-hydrodynamics (SPH) and adaptive-mesh-refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in no-radiative simulations, while it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with the observed one in the central regions of clusters. Around R_500 the temperature inhomogeneities of the SPH simulations can generate twice the typical hydrostatic-equilibrium mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

[Abridged] In the past 15 years, the concentration-mass relation has been investigated diffusely in theoretical studies. On the other hand, only recently has this relation been derived from X-ray observations. When that happened, the results caused a certain level of concern: the X-ray normalizations and slopes were found significantly dissimilar from those predicted by theory. We analyzed 52 objects, simulated each time with different physical recipes for the baryonic component, as well as 60 synthetic X-ray images, to determine if these discrepancies are real or artificial. In particular, we investigate how the simulated concentration-mass relation depends (1) on the radial range used to derive the concentration, (2) on the presence of baryons in the simulations, and on the prescription used to reproduce the gas. Finally, we evaluate (3) how the results differ when adopting an X-ray approach for the analysis and (4) how the selection functions based on X-ray luminosity can impact the results. All effects studied go in the direction of alleviating the discrepancy between observations and simulations, although with different significance: while the fitting radial range and the baryonic component play only a minor role, the X-ray approach and selection function have profound repercussion on the resulting concentration-mass relation.

Major astrophysical questions related to the formation and evolution of structures, and more specifically of galaxy groups and clusters, will still be open in the coming decade and beyond: what is the interplay of galaxy, supermassive black hole, and intergalactic gas evolution in the most massive objects in the Universe - galaxy groups and clusters? What are the processes driving the evolution of chemical enrichment of the hot diffuse gas in large-scale structures? How and when did the first galaxy groups in the Universe, massive enough to bind more than 10^7 K gas, form? Focussing on the period when groups and clusters assembled (0.5<z<2.5), we show that, due to the continuum and line emission of this hot intergalactic gas at X-ray wavelengths, Athena+, combining high sensitivity with excellent spectral and spatial resolution, will deliver breakthrough observations in view of the aforementioned issues. Indeed, the physical and chemical properties of the hot intra-cluster gas, and their evolution across time, are a key to understand the co-evolution of galaxy and supermassive black hole within their environments.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

The concentration-mass relations proposed by Prada et al. (2012) and by Duffy et al. (2008) on the scales of galaxy clusters show some of the largest discrepancies among all the works present in literature. This is surprising because they are both derived from the analysis of dark-matter halos forming in LCDM simulations with similar set-ups. With the help of analytic models and numerical simulations we investigate the origin of this discrepancy focusing on the procedures used to derive the concentrations (circular velocity ratios versus density profile fitting) and on the selection criteria used to bin the halos (binning in maximum circular velocity versus binning in mass). We find that both steps of the analysis have a large impact on the resulting c-M relation. In particular, we show that the two c-M relations can be entirely reconciled, if we account for these methodological differences. Our analysis demonstrates that the concentration estimates are sensitive to the largely different radial scales probed by a particular measurement method. This implies that concentrations derived with different techniques (both in observations and in simulations) must be compared over the same radial range.

The classification of galaxy clusters according to their X-ray appearance is a powerful tool to discriminate between regular clusters (associated to relaxed objects) and disturbed ones (linked to dynamically active systems). The compilation of the two subsamples is a necessary step both for cosmological studies - oriented towards spherical and virialized systems- and for astrophysical investigations - focused on phenomena typically present in highly disturbed clusters such as turbulence, particle re-acceleration, magneto-astrophysics . In this paper, we review several morphological parameters: asymmetry and fluctuation of the X-ray surface brightness, hardness ratios, X-ray surface-brightness concentration, centroid shift, and third-order power ratio. We test them against 60 Chandra-like images obtained from hydrodynamical simulations through the X-ray Map Simulator and visually classified as regular and disturbed. The best performances are registered when the parameters are computed using the largest possible region (either within R500 or 1000 kpc). The best indicators are the third-order power ratio, the asymmetry parameter, and the X-ray-surface-brightness concentration. All their combinations offer an efficient way to distinguish between the two morphological classes achieving values of purity extremely close to 1. A new parameter, M, is defined. It combines the strengths of the aforementioned indicators and, therefore, resulted to be the most effective parameter analyzed.

Precision cosmology studies based on wide-field surveys of galaxy clusters benefit from constraints on intrinsic scatter in mass-observable relationships. In principle, two-parameter models combining direct measurements of galaxy cluster structural variation with mass proxies such as X-ray luminosity and temperature can be used to constrain scatter in the relationship between the mass proxy and the cluster's halo mass and to measure the redshift evolution of that scatter. One candidate for quantifying cluster substructure is the ICM temperature inhomogeneity inferred from X-ray spectral properties, an example of which is T_HBR, the ratio of hardband to broadband spectral-fit temperatures. In this paper we test the effectiveness of T_HBR as an indicator of scatter in the mass-temperature relation using 118 galaxy clusters simulated with radiative cooling and feedback. We find that, while T_HBR is correlated with clusters' departures \delta lnT_X from the mean M-T_X relation, the effect is modest.

[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]

Clusters of galaxies are known to be dynamically active systems, yet X-ray studies of the low redshift population exhibit tight scaling laws. In this work, we extend previous studies of this apparent paradox using numerical simulations of two extreme merger cases, one is a high Mach number (above 2.5) satellite merger similar to the "bullet cluster" and the other a merger of nearly equal mass progenitors. Creating X-ray images densely sampled in time, we construct TX, Mgas, and YX measures within R500 and compare to the calibrations of Kravtsov et al. (2006). We find that these extreme merger cases respect the scaling relations, for both intrinsic measures and for measures derived from appropriately masked, synthetic Chandra X-ray images. The masking procedure plays a critical role in the X-ray temperature calculation while it is irrelevant in the X-ray gas mass derivation. Mis-centering up to 100 kpc does not influence the result. The observationally determined radius R500 might conduce to systematic shifts in Mgas, and YX which increase the total mass scatter.

We present the baryon fractions of 2MASS groups and clusters as a function of cluster richness using total and gas masses measured from stacked ROSAT X-ray data and stellar masses estimated from the infrared galaxy catalogs. We detect X-ray emission even in the outskirts of clusters, beyond r_200 for richness classes with X-ray temperatures above 1 keV. This enables us to more accurately determine the total gas mass in these groups and clusters. We find that the optically selected groups and clusters have flatter temperature profiles and higher stellar-to-gas mass ratios than the individually studied, X-ray bright clusters. We also find that the stellar mass in poor groups with temperatures below 1 keV is comparable to the gas mass in these systems. Combining these results with individual measurements for clusters, groups, and galaxies from the literature, we find a break in the baryon fraction at ~1 keV. Above this temperature, the baryon fraction scales with temperature as f_b \propto T^0.20\pm0.03. We see significantly smaller baryon fractions below this temperature, and the baryon fraction of poor groups joins smoothly onto that of systems with still shallower potential wells such as normal and dwarf galaxies where the baryon fraction scales with the inferred velocity dispersion as f_b \propto \sigma^1.6. The small scatter in the baryon fraction at any given potential well depth favors a universal baryon loss mechanism and a preheating model for the baryon loss. The scatter is, however, larger for less massive systems. Finally, we note that although the broken power-law relation can be inferred from data points in the literature alone, the consistency between the baryon fractions for poor groups and massive galaxies inspires us to fit the two categories of objects (galaxies and clusters) with one relation.

We aim at investigating potential biases in lensing and X-ray methods to measure the cluster mass profiles. We do so by performing realistic simulations of lensing and X-ray observations that are subsequently analyzed using observational techniques. The resulting mass estimates are compared among them and with the input models. Three clusters obtained from state-of-the-art hydrodynamical simulations, each of which has been projected along three independent lines-of-sight, are used for this analysis. We find that strong lensing models can be trusted over a limited region around the cluster core. Extrapolating the strong lensing mass models to outside the Einstein ring can lead to significant biases in the mass estimates, if the BCG is not modeled properly for example. Weak lensing mass measurements can be largely affected by substructures, depending on the method implemented to convert the shear into a mass estimate. Using non-parametric methods which combine weak and strong lensing data, the projected masses within R200 can be constrained with a precision of ~10%. De-projection of lensing masses increases the scatter around the true masses by more than a factor of two due to cluster triaxiality. X-ray mass measurements have much smaller scatter (about a factor of two smaller than the lensing masses) but they are generally biased low by 5-20%. This bias is ascribable to bulk motions in the gas of our simulated clusters. Using the lensing and the X-ray masses as proxies for the true and the hydrostatic equilibrium masses of the simulated clusters and averaging over the cluster sample we are able to measure the lack of hydrostatic equilibrium in the systems we have investigated.

The joint likelihood of observable cluster signals reflects the astrophysical evolution of the coupled baryonic and dark matter components in massive halos, and its knowledge will enhance cosmological parameter constraints in the coming era of large, multi-wavelength cluster surveys. We present a computational study of intrinsic covariance in cluster properties using halo populations derived from Millennium Gas Simulations (MGS). The MGS are re-simulations of the original 500 Mpc/h Millennium Simulation performed with gas dynamics under two different physical treatments: shock heating driven by gravity only (GO) and a second treatment with cooling and preheating (PH). We examine relationships among structural properties and observable X-ray and Sunyaev-Zel'dovich (SZ) signals for samples of thousands of halos with M_200 > 5 \times 10^{13} Msun/h and z < 2. While the X-ray scaling behavior of PH model halos at low-redshift offers a good match to local clusters, the model exhibits non-standard features testable with larger surveys, including weakly running slopes in hot gas observable--mass relations and ~10% departures from self-similar redshift evolution for 10^14 Msun/h halos at redshift z ~ 1. We find that the form of the joint likelihood of signal pairs is generally well-described by a multivariate, log-normal distribution, especially in the PH case which exhibits less halo substructure than the GO model. At fixed mass and epoch, joint deviations of signal pairs display mainly positive correlations, especially the thermal SZ effect paired with either hot gas fraction (r=0.88/0.69 for PH/GO at z=0) or X-ray temperature (r=0.62/0.83). We discuss halo mass selection by signal pairs, and find a minimum mass scatter of 4% in the \PH model by combining thermal SZ and gas fraction measurements.

The hot, diffuse gas that fills the largest overdense structures in the Universe -- clusters of galaxies and a web of giant filaments connecting them -- provides us with tools to address a wide array of fundamental astrophysical and cosmological questions via observations in the X-ray band. Clusters are sensitive cosmological probes. To utilize their full potential for precision cosmology in the following decades, we must precisely understand their physics -- from their cool cores stirred by jets produced by the central supermassive black hole (itself fed by inflow of intracluster gas), to their outskirts, where the infall of intergalactic medium (IGM) drives shocks and accelerates cosmic rays. Beyond the cluster confines lies the virtually unexplored warm IGM, believed to contain most of the baryonic matter in the present-day Universe. As a depository of all the matter ever ejected from galaxies, it carries unique information on the history of energy and metal production in the Universe. Currently planned major observatories, such as Astro-H and IXO, will make deep inroads into these areas, but to see the most interesting parts of the picture will require an almost science-fiction-grade facility with tens of m^2 of effective area, subarcsecond angular resolution, a matching imaging calorimeter and a super high-dispersion spectrograph, such as Generation-X.

Anisotropic gravitational radiation from a coalescing black hole binary is known to impart recoil velocities of up to ~1000 km/s to the remnant black hole. In this context, we study the motion of a recoiling black hole inside a galaxy modelled as an Hernquist sphere, and the signature that the hole imprints on the hot gas, using N-body/SPH simulations. Ejection of the black hole results in a sudden expansion of the gas ending with the formation of a gaseous core, similarly to what is seen for the stars. A cometary tail of particles bound to the black hole is initially released along its trail. As the black hole moves on a return orbit, a nearly spherical swarm of hot gaseous particles forms at every apocentre: this feature can live up to ~ 100 Myr. If the recoil velocity exceeds the sound speed initially, the black hole shocks the gas in the form of a Mach cone in density near each super-sonic pericentric passage. We find that the X-ray fingerprint of a recoiling black hole can be detected in Chandra X-ray maps out to a distance of Virgo. For exceptionally massive black holes the Mach cone and the wakes could be observed out to a few hundred of Mpc. Detection of the Mach cone is found to become of twofold importance: i) as a probe of high-velocity recoils and ii) as an assessment of the scatter of the mass-sigma relation at large black hole masses.

(Abridged)The X-ray measurements of the ICM metallicity are becoming more frequent due to the availability of powerful X-ray telescope with excellent spatial and spectral resolutions. The information which can be extracted from the measurements of the alpha-elements, like Oxygen, Magnesium and Silicon with respect to the Iron abundance is extremely important to better understand the stellar formation and its evolutionary history. In this paper we investigate possible source of bias connected to the plasma physics when recovering metal abundances from X-ray spectra. To do this we analyze 6 simulated galaxy clusters processed through the new version of our X-ray MAp Simulator, which allows to create mock XMM-Newton EPIC MOS1 and MOS2 observations. By comparing the spectroscopic results to the input values we find that: i) Fe is recovered with high accuracy for both hot (T>3 keV) and cold (T<2 keV) systems; at intermediate temperatures, however, we find a systematic overestimate which depends on the number counts; ii) O is well recovered in cold clusters, while in hot systems its measure may overestimate by a factor up to 2-3; iii) Being a weak line, the measurement of Mg is always difficult; despite of this, for cold systems (T<2 keV) we do not find any systematic behavior, while for very hot systems (T>5 keV) the spectroscopic measurement may be strongly overestimated up to a factor of 4; iv) Si is well recovered for all the clusters in our sample. We investigate in detail the nature of the systematic effects and biases found. We conclude that they are mainly connected with the multi-temperature nature of the projected observed spectra and to the intrinsic limitation of the XMM-Newton EPIC spectral resolution that does not always allow to disentangle among the emission lines produced by different elements.

We present a precise estimate of the bulk virial scaling relation of halos formed via hierarchical clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to cosmological parameters, the presence of a trace, dissipationless gas component, and numerical resolution down to a limit of ~1000 particles. The dark matter velocity dispersion scales with total mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026) log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly log-normal, with scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher dispersions containing 10% of the population that are merger transients. We combine this relation with the halo mass function in LCDM models, and show that a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69, favored by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with the current generation of direct simulations of cluster formation. A higher normalization, S_8 = 0.80, alleviates this tension and implies that the hot gas fraction within r_{500} is (0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent Sunyaev-Zel'dovich observations.

We explore how the behavior of galaxy cluster scaling relations are affected by flux-limited selection biases and intrinsic covariance among observable properties. Our models presume log-normal covariance between luminosity (L) and temperature (T) at fixed mass (M), centered on evolving, power-law mean relations as a function of host halo mass. Selection can mimic evolution; the \lm and \lt relations from shallow X-ray flux-limited samples will deviate from mass-limited expectations at nearly all scales while the relations from deep surveys ($10^{-14} \cgsflux$) become complete, and therefore unbiased, at masses above $\sims 2 \times 10^{14} \hinv \msol$. We derive expressions for low-order moments of the luminosity distribution at fixed temperature, and show that the slope and scatter of the \lt relation observed in flux-limited samples is sensitive to the assumed \lt correlation coefficient. In addition, \lt covariance affects the redshift behavior of halo counts and mean luminosity in a manner that is nearly degenerate with intrinsic population evolution.

The origin of the recently reported systematic bias in the spectroscopic temperature of galaxy clusters is investigated using cosmological hydrodynamical simulations. We find that the local inhomogeneities of the gas temperature and density, after corrected for the global radial profiles, have nearly a universal distribution that resembles the log-normal function. Based on this log-normal approximation for the fluctuations in the intra-cluster medium, we develop an analytical model that explains the bias in the spectroscopic temperature discovered recently. We conclude that the multi-phase nature of the intra-cluster medium not only from the radial profiles but also from the local inhomogeneities plays an essential role in producing the systematic bias.

Using the results of a high-resolution, cosmological hydrodynamical re-simulation of a supercluster-like region we investigate the physical properties of the gas located along the filaments and bridges which constitute the so-called cosmic web. First we analyze the main characteristics of the density, temperature and velocity fields, which have quite different distributions, reflecting the complex dynamics of the structure formation process. Then we quantify the signals which originate from the matter in the filaments by considering different observables. Inside the cosmic web, we find that the halo density is about 4-6 times larger than in the neighbouring region; the bremsstrahlung X-ray surface brightness reaches at most 10e-16 erg/s/cm^2/armin^2; the Compton-y parameter due to the thermal Sunyaev-Zel'dovich effect is about 10e-6; the reduced shear produced by the weak lensing effect is ~ 0.01-0.02. These results confirm the difficulty of an observational detection of the cosmic web. Finally we find that projection effects of the filamentary network can affect the estimates of the properties of single clusters, increasing their X-ray luminosity by less than 10% and their central Compton-y parameter by up to 30% cent.

We examine the systematics affecting the X-ray mass estimators applied to a set of five simulated galaxy clusters. They have been processed through the X-ray Map Simulator, X-MAS, to provide Chandra-like long exposures that are analyzed to reconstruct the gas temperature, density, and mass profiles used as input. We find that at R_2500 the mass profile obtained via a direct application of the hydrostatic equilibrium equation is consistent within 1 sigma with the actual mass; although we notice this estimator shows high statistical errors due to high level of Chandra background. Instead, the poorness of the beta-model in describing the gas density profile makes the evaluated masses to be underestimated by \sim 40 per cent with respect to the true mass, both with an isothermal and a polytropic temperature profile. We also test ways to recover the mass by adopting an analytic mass model, such as those proposed by Navarro et al. (1997) and Rasia et al. (2004), and fitting the temperature profile expected from the hydrostatic equilibrium equation to the observed one. These methods and the one of the hydrostatic equilibrium equation provide a more robust mass estimation than the ones based on the beta-model. In the present work the main limitation for a precise mass reconstruction is to ascribe to the relatively high level of the background chosen to reproduce the Chandra one. After artificially reducing it by a factor of 100, we find that the estimated mass significantly underestimates the true mass profiles. This is manly due (i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to the bias towards lower values of the X-ray temperature measurements because of the complex thermal structure of the emitting plasma.

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. Here we show that, if the intra-cluster medium is thermally complex, T_ew fails at reproducing T_spec. We propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we also find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and beta-model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X--ray temperature function could bias low sigma_8 by 10-20 per cent. This may alleviate the tension between the value of sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is thermally complex, T_ew fails at reproducing T_spec, and propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T_sl relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and \beta--model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X-ray temperature function could bias low \sigma_8 by 10-20 per cent. This may alleviate the tension between the value of \sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

Theoretical studies of the physical processes in clusters of galaxies are mainly based on the results of numerical simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We show that the projected spectroscopic temperature of clusters obtained from X-ray observations is always lower than the emission-weighed temperature. This bias is related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. In real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps. We propose a new formula, the spectroscopic-like temperature function that better approximates the spectroscopic temperature, making simulations more directly comparable to observations

Theoretical studies of the physical processes guiding the formation and evolution of galaxies and galaxy clusters in the X-ray are mainly based on the results of numerical hydrodynamical N-body simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We demonstrate that the projected spectroscopic temperature of thermally complex clusters obtained from X-ray observations is always lower than the emission-weighed temperature, which is widely used in the analysis of numerical simulations. We show that this temperature bias is mainly related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has important implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. Because of this bias, in real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps, and may even not be detected. This may explain why, although numerical simulations predict that shock fronts are a quite common feature in clusters of galaxies, to date there are very few observations of objects in which they are clearly seen. To fix this problem we propose a new formula, the spectroscopic-like temperature function, and show that, for temperature larger than 3 keV, it approximates the spectroscopic temperature better than few per cent, making simulations more directly comparable to observations.

Using high-resolution non-radiative hydrodynamic simulations of galaxy clusters we obtain simple analytic formulae for DM and gas distribution in the spherical approximation. We derive fits for the DM density, velocity dispersion and velocity anisotropy. We use these models to test the dynamical equilibrium hypothesis through the Jeans equation: we find that this is satisfied to good accuracy. This result show that our fits constitute a self-consistent dynamical model. We then extend our analysis to the gas, studying its density, temperature and velocity structure, with no further hypothesis on the dynamical status or state equation. Gas and DM show self-similar density profiles down to 0.06 R_vir, while at smaller radii the gas produces a central core. Gas temperatures are almost isothermal out to 0.2 R_vir, then steeply decrease reaching at R_vir a value a factor of 2 lower. We find that the gas is not at rest inside R_vir: velocity dispersions are increasing functions of the radius, motions are isotropic to slightly tangential, and contribute to the total pressure support. We test our model using a generalization of the hydrostatic equilibrium equation, where the gas motion is properly taken into account. We find that our fits provide an accurate description of the system: the gas is in equilibrium and is a good tracer the overall cluster potential if all terms are taken into account, while simpler assumptions, as the beta-model, cause systematic mass underestimates. We also find that, if gas velocities are neglected, then a simple isothermal model fares better at large radii than a non-isothermal one. (abridged)

We study the formation of clusters of galaxies using high-resolution hydrodynamic cosmological simulations that include the effect of thermal conduction with an effective isotropic conductivity of 1/3 the classical Spitzer value. We find that, both for a hot ($T_{\rm ew}\simeq 12$ keV) and several cold ($T_{\rm ew}\simeq 2$ keV) galaxy clusters, the baryonic fraction converted into stars does not change significantly when thermal conduction is included. However, the temperature profiles are modified, particularly in our simulated hot system, where an extended isothermal core is readily formed. As a consequence of heat flowing from the inner regions of the cluster both to its outer parts and into its innermost resolved regions, the entropy profile is altered as well. This effect is almost negligible for the cold cluster, as expected based on the strong temperature dependence of the conductivity. Our results demonstrate that while thermal conduction can have a significant influence on the properties of the intra--cluster medium of rich galaxy clusters, it appears unlikely to provide by itself a solution for the overcooling problem in clusters, or to explain the current discrepancies between the observed and simulated properties of the intra--cluster medium.

We study the impact of merger events on the strong lensing properties of galaxy clusters. Previous lensing simulations were not able to resolve dynamical time scales of cluster lenses, which arise on time scales which are of order a Gyr. In this case study, we first describe qualitatively with an analytic model how some of the lensing properties of clusters are expected to change during merging events. We then analyse a numerically simulated lens model for the variation in its efficiency for producing both tangential and radial arcs while a massive substructure falls onto the main cluster body. We find that: (1) during the merger, the shape of the critical lines and caustics changes substantially; (2) the lensing cross sections for long and thin arcs can grow by one order of magnitude and reach their maxima when the extent of the critical curves is largest; (3) the cross section for radial arcs also grows, but the cluster can efficiently produce this kind of arcs only while the merging substructure crosses the main cluster centre; (4) while the arc cross sections pass through their maxima as the merger proceeds, the cluster's X-ray emission increases by a factor of $\sim5$. Thus, we conclude that accounting for these dynamical processes is very important for arc statistics studies. In particular, they may provide a possible explanation for the arc statistics problem.

Numerical hydro-N-body simulations are very important tools for making theoretical predictions for the formation of galaxy clusters. They show that the atmospheres of clusters of galaxies have quite complex angular and thermal structures. The full understanding of the physical processes behind these features can be only achieved by direct comparison of observations to hydro-N-body simulations. Although simple in principle, these comparisons are not always trivial. In fact, real data are convolved with the instrument response which may substantially influence the apparent properties of the studied features. To overcome this problem we build the software package X-MAS devoted to simulate X-ray observations of galaxy clusters obtained from hydro-N-body simulations. In this paper we present how this software package works and discuss its application to the simulation of Chandra ACIS-S3 observations. We compare some of the main physical properties of the input data to the ones derived from simulated observations after performing a standard imaging and spectral analysis. We show that, if the thermal structure of the cluster along a particular line of sight is quite complex, the projected spectroscopic temperature obtained from the observation is significantly lower than the emission-weighed value inferred directly from hydrodynamical simulation. This implies that much attention must be paid in the theoretical interpretation of observational temperatures.

Arc statistics is known to be a powerful cosmological tool. Numerical lensing simulations show that orders of magnitude differences in the number of {\em giant} arcs on the whole sky are expected in different cosmological models. In this paper, we discuss the analytic and numerical methods in arc statistics and show that analytic models fail to reproduce the efficiency for strong lensing of more realistic numerical cluster models. Then, we discuss two recent extensions of the lensing simulations, i.e. the effects of cD galaxies in the lensing clusters and the impact of cluster mergers on arc statistics. We show that cD galaxies can increase the lensing cross sections for long and thin arcs by perhaps up to $\sim50%$, while major mergers can change the cluster efficiency for producing such arcs by up to one order of magnitude.