Search

Context. To assume hydrostatic equilibrium between the intracluster medium and the gravitational potential of galaxy clusters is an extensively used method to investigate their total masses. Aims. We want to test hydrostatic masses obtained with an observational code in the context of the Spectrum-Roentgen-Gamma/eROSITA survey. Methods. We used the hydrostatic modeling code MBProj2 to fit surface-brightness profiles to simulated clusters with idealized properties as well as to a sample of 93 clusters taken from the Magneticum Pathfinder simulations. We investigated the latter under the assumption of idealized observational conditions and also for realistic eROSITA data quality. The comparison of the fitted cumulative total mass profiles and the true mass profiles provided by the simulations allows us to gain knowledge both about the validity of hydrostatic equilibrium in each cluster and the reliability of our approach. Furthermore, we used the true profiles for gas density and pressure to compute hydrostatic mass profiles based on theory for every cluster. Results. For an idealized cluster that was simulated to fulfill perfect hydrostatic equilibrium, we find that the cumulative total mass at the true r500 and r200 can be reproduced with deviations of less than 7%. For the clusters from the Magneticum Pathfinder simulations under idealized observational conditions, the median values of the fitted cumulative total masses at the true r500 and r200 are in agreement with our expectations, taking into account the hydrostatic mass bias. Nevertheless, we find a tendency towards steeper cumulative total mass profiles in the outskirts than expected. For realistic eROSITA data quality, this steepness problem intensifies for clusters with high redshifts and leads to excessive cumulative total masses at r200. For the hydrostatic masses based on the true profiles known from the simulations, we find good agreement with our expectations concerning the hydrostatic mass.

Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

We report results from deep Suzaku and mostly snapshot Chandra observations of four nearby galaxy groups: MKW4, Antlia, RXJ1159+5531, and ESO3060170. Their peak temperatures vary over 2-3 keV, making them the smallest systems with gas properties constrained to their viral radii. The average Fe abundance in the outskirts (R $>$ 0.25R$_{200}$) of their intragroup medium (IGrM) is $Z_{\rm Fe}=0.309\pm0.018$ $Z_\odot$ with $\chi^2$ = 14 for 12 degrees of freedom, which is remarkably uniform and strikingly similar to that of massive galaxy clusters, and is fully consistent with the numerical predictions from the IllustrisTNG cosmological simulation. Our results support an early-enrichment scenario among galactic systems over an order of magnitude in mass, even before their formation. When integrated out to R$_{200}$, we start to see a tension between the measured Fe content in ICM and what is expected from supernovae yields. We further constrain their O, Mg, Si, S, and Ni abundances. The abundance ratios of those elements relative to Fe are consistent with the predictions (if available) from IllustrisTNG. Their Type Ia supernovae fraction varies between 14%-21%. A pure core collapsed supernovae enrichment at group outskirts can be ruled out. Their cumulative iron-mass-to-light ratios within R$_{200}$ are half that of the Perseus cluster, which may imply that galaxy groups do not retain all of their enriched gas due to their shallower gravitational potential wells, or that groups and clusters may have different star formation histories., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Major mergers between galaxy clusters can produce large turbulent and bulk flow velocities in the intra-cluster medium and thus imprint diagnostic features in X-ray spectral emission lines from heavy ions. As demonstrated by Hitomi in observations of the Perseus cluster, measurements of gas velocities in clusters from high-resolution X-ray spectra will be achievable with upcoming X-ray calorimeters like those on board XRISM, Athena, or a Lynx like mission. We investigate this possibility for interesting locations across a major cluster merger from a hydrodynamical simulation, via X-ray synthetic spectra with a few eV energy resolution. We observe the system from directions perpendicular to the plane of the merger and along the merger axis. In these extreme geometrical configurations, we find clear non-Gaussian shapes of the iron He-like K_alpha line at 6.7keV. The velocity dispersion predicted from the simulations can be retrieved for the brightest 100ks pointings with XRISM Resolve, despite some discrepancy related to the complex non-Gaussian line shapes. Measurements in faint regions require however high S/N and the larger collecting area of the Athena X-IFU calorimeter is thus needed. With the latter, we also investigate the gas temperature and velocity gradient across the merger bow shock edge, from 20"-wide annuli extracted from a single 1Ms X-IFU pointing. We find best-fit temperature and velocity dispersion values that are consistent with predictions from the simulations within 1-sigma, but the uncertainties on the inferred velocity dispersion are too large to place any stringent constraints on the shallow gradient downstream of the shock. We also present simulated images of the thermal and kinetic Sunyaev-Zeldovich effects, using the above viewing configurations, and compare the results at angular resolutions appropriate for future observatories such as CMB-S4 and AtLAST., Comment: 22 pages, 18 figures, re-submitted after minor revision following referee report

Galaxy groups host the majority of matter and more than half of all the galaxies in the Universe. Their hot ($10^7$ K), X-ray emitting intra-group medium (IGrM) reveals emission lines typical of many elements synthesized by stars and supernovae. Because their gravitational potentials are shallower than those of rich galaxy clusters, groups are ideal targets for studying, through X-ray observations, feedback effects, which leave important marks on their gas and metal contents. Here, we review the history and present status of the chemical abundances in the IGrM probed by X-ray spectroscopy. We discuss the limitations of our current knowledge, in particular due to uncertainties in the modeling of the Fe-L shell by plasma codes, and coverage of the volume beyond the central region. We further summarize the constraints on the abundance pattern at the group mass scale and the insight it provides to the history of chemical enrichment. Parallel to the observational efforts, we review the progress made by both cosmological hydrodynamical simulations and controlled high-resolution 3D simulations to reproduce the radial distribution of metals in the IGrM, the dependence on system mass from group to cluster scales, and the role of AGN and SN feedback in producing the observed phenomenology. Finally, we highlight future prospects in this field, where progress will be driven both by a much richer sample of X-ray emitting groups identified with eROSITA, and by a revolution in the study of X-ray spectra expected from micro-calorimeters onboard XRISM and ATHENA., Comment: 43 pages, 15 figures. Accepted for publication in Universe. This review article is part of the special issue "The Physical Properties of the Groups of Galaxies", edited by L. Lovisari and S. Ettori. Published in MDPI - Universe https://www.mdpi.com/journal/universe/special_issues/PPGG"

T. H. Reiprich et al., arXiv:2012.08491v1, [Context]: Inferences about dark matter, dark energy, and the missing baryons all depend on the accuracy of our model of large-scale structure evolution. In particular, with cosmological simulations in our model of the Universe, we trace the growth of structure, and visualize the build-up of bigger structures from smaller ones and of gaseous filaments connecting galaxy clusters., [Aims]: Here we aim to reveal the complexity of the large-scale structure assembly process in great detail and on scales from tens of kiloparsecs up to more than 10 Mpc with new sensitive large-scale observations from the latest generation of instruments. We also aim to compare our findings with expectations from our cosmological model., [Methods]: We used dedicated SRG/eROSITA performance verification (PV) X-ray, ASKAP/EMU Early Science radio, and DECam optical observations of a ~15 deg2 region around the nearby interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process, the morphological complexity in the inner parts of the clusters, and the (re-)acceleration of plasma. We also used complementary Sunyaev-Zeldovich (SZ) effect data from the Planck survey and custom-made Galactic total (neutral plus molecular) hydrogen column density maps based on the HI4PI and IRAS surveys. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite., [Results]: We trace the irregular morphology of warm and hot gas of the main clusters from their centers out to well beyond their characteristic radii, r200. Between the two main cluster systems, we observe an emission bridge on large scale and with good spatial resolution. This bridge includes a known galaxy group but this can only partially explain the emission. Most gas in the bridge appears hot, but thanks to eROSITA’s unique soft response and large field of view, we discover some tantalizing hints for warm, truly primordial filamentary gas connecting the clusters. Several matter clumps physically surrounding the system are detected. For the “Northern Clump,” we provide evidence that it is falling towards A3391 from the X-ray hot gas morphology and radio lobe structure of its central AGN. Moreover, the shapes of these X-ray and radio structures appear to be formed by gas well beyond the virial radius, r100, of A3391, thereby providing an indirect way of probing the gas in this elusive environment. Many of the extended sources in the field detected by eROSITA are also known clusters or new clusters in the background, including a known SZ cluster at redshift z = 1. We find roughly an order of magnitude more cluster candidates than the SPT and ACT surveys together in the same area. We discover an emission filament north of the virial radius of A3391 connecting to the Northern Clump. Furthermore, the absorption-corrected eROSITA surface brightness map shows that this emission filament extends south of A3395 and beyond an extended X-ray-emitting object (the “Little Southern Clump”) towards another galaxy cluster, all at the same redshift. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA PV observation field. The Northern and Southern Filament are each detected at >4σ. The Planck SZ map additionally appears to support the presence of both new filaments. Furthermore, the DECam galaxy density map shows galaxy overdensities in the same regions. Overall, the new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution; the similarities of features found in a similar system in the Magneticum simulation are striking. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure., N.O. acknowledges support by a Fund for the Promotion of Joint International Research (JSPS KAKENHI Grant Number 16KK0101). K.D. acknowledges support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 3907833. V.B. acknowledges support by the DFG project nr. 415510302 (BI 1699/2-1). Partial support for LR comes from US National Science Foundation grant AST17-14205 to the University of Minnesota.

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., 23 pages, 23 figures; submitted to A&A

We investigated the origin and gas properties of a simulated galaxy cluster pair, to connect simulation predictions to the SRG/eROSITA X-ray observations of the Abell 3391/95 field. The simulated system has been extracted from the (352 cMpc/h)^3 volume of the Magneticum Pathfinder cosmological simulations at z=0.07. We tracked back in time the main progenitors of the pair clusters and surrounding groups to study the assembly history of the system and its evolution. Similarly to the observed A3391/95 system, the simulated pair is embedded in a complex network of gas filaments, with structures aligned over more than 20 projected Mpc and the whole region collapsing towards the central overdense node. The spheres of influence (3*R200) of the two main clusters already overlap at z=0.07, but their virial boundaries are still physically separated. The diffuse gas located in the interconnecting bridge closely reflects the WHIM, with typical temperature of ~1 keV and overdensity $\sim 100$, with respect to the mean baryon density of the Universe, and lower enrichment level compared to the ICM in clusters. We find that most of the bridge gas collapsed from directions roughly orthogonal to the intra-cluster gas accretion directions, and its origin is mostly unrelated to the two cluster progenitors. We find clear signatures in the surrounding groups of infall motion towards the pair, such as significant radial velocities and slowdown of gas compared to dark matter. These findings further support the picture of the Northern Clump (MCXC J0621.7-5242) cluster infalling along a cosmic gas filament towards Abell 3391, possibly merging with it. We conclude that, in such a configuration, the pair clusters of the A3391/95-like system are in a pre-merger phase, and did not interact yet. The diffuse gas in the interconnecting bridge is mostly warm filament gas, rather than tidally-stripped cluster gas., 19 pages, 15 figures in the main text; 2 pages, 2 figures as appendices. Improved version after referee report. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

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., Comment: 21 pages, 17 figure, recommended for publication in A&A

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., 24 pages, 11 figures, 5 tables, replaced to match accepted version

This article describes a data centre hosting a web portal for accessing and sharing the output of large, cosmological, hydro-dynamical simulations with a broad scientific community. It also allows users to receive related scientific data products by directly processing the raw simulation data on a remote computing cluster. The data centre has a multi-layer structure: a web portal, a job control layer, a computing cluster and a HPC storage system. The outer layer enables users to choose an object from the simulations. Objects can be selected by visually inspecting 2D maps of the simulation data, by performing highly compounded and elaborated queries or graphically by plotting arbitrary combinations of properties. The user can run analysis tools on a chosen object. These services allow users to run analysis tools on the raw simulation data. The job control layer is responsible for handling and performing the analysis jobs, which are executed on a computing cluster. The innermost layer is formed by a HPC storage system which hosts the large, raw simulation data. The following services are available for the users: (I) ClusterInspect visualizes properties of member galaxies of a selected galaxy cluster; (II) SimCut returns the raw data of a sub-volume around a selected object from a simulation, containing all the original, hydro-dynamical quantities; (III) Smac creates idealized 2D maps of various, physical quantities and observables of a selected object; (IV) Phox generates virtual X-ray observations with specifications of various current and upcoming instruments.

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., Comment: 16 pages, 17 figures, 3 tables

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermodynamics, kinematics, and chemical composition of the circumgalactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying at least 1600 square degrees over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view ($\sim1$ deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5 arcsec in the soft X-ray band, and a low and stable instrumental background ensuring high sensitivity to faint, extended emission., Comment: White paper submitted in response to ESA's Voyage 2050 Call. Accepted for publication in Experimental Astronomy

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data., 11 pages, 6 figures, MNRAS accepted

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), 14 pages, 7 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics. v3: added missing references

In this study, we present a modelling of the X-ray emission from the simulated SMBHs within the cosmological hydrodynamical Magneticum Pathfinder Simulation, in order to study the statistical properties of the resulting X-ray Active Galactic Nuclei (AGN) population and their expected contribution to the X-ray flux from galaxy clusters. The simulations reproduce the evolution of the observed unabsorbed AGN bolometric luminosity functions (LFs) up to redshift z~2, consistently with previous works. Furthermore, we study the evolution of the LFs in the soft (SXR) and hard (HXR) X-ray bands by means of synthetic X-ray data generated with the PHOX simulator, that includes an observationally-motivated modelling of an instrinsic absorption component, mimicking the torus around the AGN. The reconstructed SXR and HXR AGN LFs present a remarkable agreement with observational data up to z~2 when an additional obscuration fraction for Compton-thick AGN is assumed, although a discrepancy still exists for the SXR LF at z=2.3. With this approach, we also generate full eROSITA mock observations to predict the level of contamination due to AGN of the intra-cluster medium (ICM) X-ray emission, which can affect cluster detection especially at high redshifts. We find that, at z~1-1.5, for 20-40% of the clusters with M500>3e13 Msun/h, the AGN counts in the observed SXR band exceed by more than a factor of 2 the counts from the whole ICM., in press, MNRAS Main Journal

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 (, 22 pages, 15 figures

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., Comment: 11 pages, 7 figures, A&A in press. This is one of three companion papers (including Ghirardini et al. and Ettori et al.) on the X-COP sample. X-COP data products are available for download on https://www.astro.unige.ch/xcop

The distribution of chemical elements in the hot intracluster medium (ICM) retains valuable information about the enrichment and star formation histories of galaxy clusters, and on the feedback and dynamical processes driving the evolution of the cosmic baryons. In the present study we review the progresses made so far in the modelling of the ICM chemical enrichment in a cosmological context, focusing in particular on cosmological hydrodynamical simulations. We will review the key aspects of embedding chemical evolution models into hydrodynamical simulations, with special attention to the crucial assumptions on the initial stellar mass function, stellar lifetimes and metal yields, and to the numerical limitations of the modelling. At a second stage, we will overview the main simulation results obtained in the last decades and compare them to X-ray observations of the ICM enrichment patterns. In particular, we will discuss how state-of-the-art simulations are able to reproduce the observed radial distribution of metals in the ICM, from the core to the outskirts, the chemical diversity depending on cluster thermo-dynamical properties, the evolution of ICM metallicity and its dependency on the system mass from group to cluster scales. Finally, we will discuss the limitations still present in modern cosmological, chemical, hydrodynamical simulations and the perspectives for improving the theoretical modelling of the ICM enrichment in galaxy clusters in the future., Comment: 38 pages. Review paper. Accepted for publication on Space Science Reviews. This is the companion review of "Enrichment of the hot intracluster medium: observations"

The X-ray Integral Field Unit (X-IFU) is the cryogenic imaging spectrometer onboard the ESA L2 mission Athena. With its array of almost 3840 superconducting Transition Edge Sensors micro-calorimeters, the X-IFU will provide spatially resolved (5" over the field of view) high-resolution spectroscopy (2.5 eV FWHM up to 7 keV) in the 0.2-12 keV energy band. These transformational capabilities will allow the X-IFU to probe the Hot and Energetic Universe, and notably measure the physical properties of large-scale structures with unprecedented accuracy. Starting from numerically-simulated massive (1014 Mo) galaxy clusters at different steps of their evolution, we investigate the capabilities of the X-IFU in recovering chemical abundances, redshift and gas temperature spatial distributions across time, making use of full field-of-view End-To-End simulations of X-IFU observations. This work serve as feasibility study for the Chemical Enrichment of the Universe science objective. We show that using 100 ks observations, the X-IFU will provide an unprecedented spatially-accurate knowledge of the physics of the ICM (abundances, temperature, bulk-motion). We also demonstrate that challenges related to the data analysis of extended sources with very high-resolution spectrometers (e.g. binning, line of sight mixing, particle background) need to be thoroughly addressed to maximise the science of the instrument.

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, MNRAS, accepted

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., 18 pages; submitted to MNRAS

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., Comment: 22 pages, 14 figures, 2 tables. Submitted to MNRAS

Why is our universe as we observe it? Will it be the same forever? Understanding the nature of the main constituents of the universe is crucial to obtain a precise description of the way in which it reached its present state. Nowadays, many independent observations support a picture in which the matter content of the universe is shared between an ordinary and observable baryonic component ( ~ 5?%) and an invisible dark matter ( ~ 23?%). The remaining ~ 72?% of the universe content is in the form of a completely mysterious dark energy field. This composition emphasizes that, while ~ 95?% of our universe represents a major uncertainty for us, even the minor contribution from normal and, apparently, known matter entails important challenges for cosmologists.

We study the thermal structure of the intra-cluster medium (ICM) in a set of cosmological hydrodynamical cluster simulations performed with an SPH numerical scheme employing an artificial conductivity (AC) term. We explore the effects of this term on the ICM temperature and entropy profiles, thermal distribution, velocity field and expected X-ray emission. We find that in adiabatic runs the artificial conductivity favours (i) the formation of an entropy core, raising and flattening the central entropy profiles, in better agreement with findings from Eulerian codes; and (ii) a systematic reduction of the cold gas component. In fact, the cluster large-scale structure and dynamical state are preserved across different runs, but the improved gas mixing enabled by the AC term strongly increases the stripping rate of gas from the cold clumps moving through the ICM. This in turn reduces the production of turbulence generated by the instabilities which develop because of the interaction between clumps and ambient ICM. We then find that turbulent motions, enhanced by the time-dependent artificial viscosity scheme we use, are rather damped by the AC term. The ICM synthetic X-ray emission substantially mirrors the changes in its thermo-dynamical structure, stressing the robustness of the AC impact. All these effects are softened by the introduction of radiative cooling but still present, especially a partial suppression of cold gas. Therefore, not only the physics accounted for, but also the numerical approach itself can have an impact in shaping the ICM thermo-dynamical structure and ultimately in the use of SPH cluster simulations for cosmological studies., Comment: 23 pages, 17 figures; MNRAS, in press

We study the baryonic, chemical and dynamical properties of a significantly large sample of early proto-galaxies in the first 500 Myr of the Universe (redshift z>9), obtained from high-resolution numerical, N-body, hydrodynamical, chemistry simulations including atomic and molecular networks, gas cooling, star formation, stellar evolution and metal spreading for population III and population II-I regimes according to proper stellar yields and lifetimes. We find that first star formation events take place in halos with dark-matter mass M_{DM} > 2\times 10^6 M_{sun}. Early star forming objects have: molecular fractions from x_{mol} < 10^{-4} in quiescent structures up to x_{mol} > 0.1 in active regions; star formation rates SFR ~10^{-8}-10^{-3}M_{sun}/yr; and metallicities in the range ~10^{-8}-10^{-2}Z_{sun}. Roughly ~10% of high-z haloes host population II-I star formation and dominate the cosmic SFR density. They usually are bursty objects with mean specific SFR around ~10 Gyr^{-1} at z~9 and increasing with redshift up to ~10^2 Gyr^{-1}. Stellar feedback effects alter the baryonic content of the haloes and locally affect their chemical and thermo-dynamical properties, as reflected by the broadening of various physical relations. The establishment of gaseous rotationally-supported cores is quite uncommon, weakly related to the large-scale dark-matter behaviour and evolving in an intermittent fashion. The colder, molecular-rich phase tends to maintain any established rotational motion longer with respect to the hotter, metal-rich component, which is very sensitive to environmental processes. While the fraction of haloes featuring a significant amount of co-rotating, molecular-rich gas increases with cosmic time (from a few percent at z~20 up to ~5-15% at z~9), the chaotic nature of metal-enriched material does not lead to particular trends., Comment: 20 pages, 19 figures; MNRAS, in press

Non-thermal motions in the intra-cluster medium (ICM) are believed to play a non-negligible role in the pressure support to the total gravitating mass of galaxy clusters. Future high-precision X-ray missions will allow us to directly detect the signature of these motions from highly resolved spectra of the ICM. With this perspective, we will present here a study on a set of clusters extracted from a cosmological hydrodynamical simulation, for which we simulated X-ray observations and characterized them via the non-thermal velocity amplitude of the ICM. We have applied the X-ray virtual telescope PHOX to generate mock observations of the simulated clusters and estimated in an observational fashion global X-ray properties. We included ICM non-thermal motions as additional parameter to further identify non-regular clusters. In particular, the influence of internal gas motions on the location of clusters in the LX-T plane was investigated and was found that objects with the largest amount of non-thermal velocity have on average higher X-ray luminosities. (© 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

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., 14 pages, 13 figures, 1 table; submitted

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., Comment: 19 pages, 12 figures. Accepted for publication in MNRAS

Gas motions in the hot intracluster medium of galaxy clusters have an important effect on the mass determination of the clusters through X-ray observations. The corresponding dynamical pressure has to be accounted for in addition to the hydrostatic pressure support to achieve a precise mass measurement. An analysis of the velocity structure of the ICM for simulated cluster-size haloes, especially focusing on rotational patterns, has been performed, demonstrating them to be an intermittent phenomenon, strongly related to the internal dynamics of substructures. We find that the expected build-up of rotation due to mass assembly gets easily destroyed by passages of gas-rich substructures close to the central region. Though, if a typical rotation pattern is established, the corresponding mass contribution is estimated to be up to ~17% of the total mass in the innermost region, and one has to account for it. Extending the analysis to a larger sample of simulated haloes we statistically observe that (i) the distribution of the rotational component of the gas velocity in the innermost region has typical values of ~200-300 km/s; (ii) except for few outliers, there is no monotonic increase of the rotational velocity with decreasing redshift, as we would expect from approaching a relaxed configuration. Therefore, the hypothesis that the build-up of rotation is strongly influenced by internal dynamics is confirmed, and minor events like gas-rich substructures passing close to the equatorial plane can easily destroy any ordered rotational pattern.

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., Comment: 6 pages, 4 figures, accepted in ApJL, v2 contains some modifications on the text (results unchanged)

Revealing hidden patterns in astronomical data is often the path to fundamental scientific breakthroughs; meanwhile the complexity of scientific inquiry increases as more subtle relationships are sought. Contemporary data analysis problems often elude the capabilities of classical statistical techniques, suggesting the use of cutting edge statistical methods. In this light, astronomers have overlooked a whole family of statistical techniques for exploratory data analysis and robust regression, the so-called Generalized Linear Models (GLMs). In this paper -- the first in a series aimed at illustrating the power of these methods in astronomical applications -- we elucidate the potential of a particular class of GLMs for handling binary/binomial data, the so-called logit and probit regression techniques, from both a maximum likelihood and a Bayesian perspective. As a case in point, we present the use of these GLMs to explore the conditions of star formation activity and metal enrichment in primordial minihaloes from cosmological hydro-simulations including detailed chemistry, gas physics, and stellar feedback. We predict that for a dark mini-halo with metallicity $\approx 1.3 \times 10^{-4} Z_{\bigodot}$, an increase of $1.2 \times 10^{-2}$ in the gas molecular fraction, increases the probability of star formation occurrence by a factor of 75%. Finally, we highlight the use of receiver operating characteristic curves as a diagnostic for binary classifiers, and ultimately we use these to demonstrate the competitive predictive performance of GLMs against the popular technique of artificial neural networks., Comment: 20 pages, 10 figures, 3 tables, accepted for publication in Astronomy and Computing

X-ray astronomy is an important tool in the astrophysicist's toolkit to investigate high-energy astrophysical phenomena. Theoretical numerical simulations of astrophysical sources are fully three-dimensional representations of physical quantities such as density, temperature, and pressure, whereas astronomical observations are two-dimensional projections of the emission generated via mechanisms dependent on these quantities. To bridge the gap between simulations and observations, algorithms for generating synthetic observations of simulated data have been developed. We present an implementation of such an algorithm in the yt analysis software package. We describe the underlying model for generating the X-ray photons, the important role that yt and other Python packages play in its implementation, and present a detailed workable example of the creation of simulated X-ray observations., Fixed typos from first version. 8 pages, 5 figures. Appears in the Proceedings of the 13th Python in Science Conference (SciPy 2014) at http://conference.scipy.org/proceedings/scipy2014/zuhone.html with video of the talk at http://youtu.be/fUMq6rmNshc

We present the X-ray properties and scaling relations of a large sample of clusters extracted from the Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC) dataset. We focus on a sub-sample of 179 clusters at redshift z~0.11, with 3.2e14M_sun/h, Comment: 17 pages, 14 figures; MNRAS, in press

We present a novel approach, based on robust principal components analysis (RPCA) and maximal information coefficient (MIC), to study the redshift dependence of halo baryonic properties. Our data are composed of a set of different physical quantities for primordial minihaloes: dark-matter mass ($M_{\mathrm{dm}}$), gas mass ($M_{\mathrm{gas}}$), stellar mass ($M_{\mathrm{star}}$), molecular fraction (${\mathrm{x_{mol}}}$), metallicity ($Z$), star formation rate (SFR) and temperature. We find that $M_{\mathrm{dm}}$ and $M_{\mathrm{gas}}$ are dominant factors for variance, particularly at high redshift. Nonetheless, with the emergence of the first stars and subsequent feedback mechanisms, ${\mathrm{x_{mol}}}$, SFR and $Z$ start to have a more dominant role. Standard PCA gives three principal components (PCs) capable to explain more than 97 per cent of the data variance at any redshift (two PCs usually accounting for no less than 92 per cent), whilst the first PC from the RPCA analysis explains no less than 84 per cent of the total variance in the entire redshift range (with two PCs explaining $ \gtrsim 95$ per cent anytime). Our analysis also suggests that all the gaseous properties have a stronger correlation with $M_{\mathrm{gas}}$ than with $M_{\mathrm{dm}}$, while $M_{\mathrm{gas}}$ has a deeper correlation with ${\mathrm{x_{mol}}}$ than with $Z$ or SFR. This indicates the crucial role of gas molecular content to initiate star formation and consequent metal pollution from Population III and Population II/I regimes in primordial galaxies. Finally, a comparison between MIC and Spearman correlation coefficient shows that the former is a more reliable indicator when halo properties are weakly correlated., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

We present a novel, virtual X-ray observatory designed to obtain synthetic observations from hydro-numerical simulations, named PHOX. In particular, we provide a description of the code constituting the photon simulator and of the new approach implemented. We apply PHOX to simulated galaxy clusters in order to demonstrate its capabilities. In fact, X-ray observations of clusters of galaxies continue to provide us with an increasingly detailed picture of their structure and of the underlying physical phenomena governing the gaseous component, which dominates their baryonic content. Therefore, it is fundamental to find the most direct and faithful way to compare such observational data with hydrodynamical simulations of cluster-like objects, which can currently include various complex physical processes. Here, we present and analyse synthetic Suzaku observations of two cluster-size haloes obtained by processing with PHOX the hydrodynamical simulation of the large-scale, filament-like region in which they reside. Taking advantage of the simulated data, we test the results inferred from the X-ray analysis of the mock observations against the underlying, known solution. Remarkably, we are able to recover the theoretical temperature distribution of the two haloes by means of the multi-temperature fitting of the synthetic spectra. Moreover, the shapes of the reconstructed distributions allow us to trace the different thermal structure that distinguishes the dynamical state of the two haloes.

Non-thermal motions in the intra-cluster medium (ICM) are believed to play a non-negligible role in the pressure support to the total gravitating mass of galaxy clusters. Future X-ray missions, such as ASTRO-H and ATHENA, will eventually allow us to directly detect the signature of these motions from high-resolution spectra of the ICM. In this paper, we present a study on a set of clusters extracted from a cosmological hydrodynamical simulation, devoted to explore the role of non-thermal velocity amplitude in characterising the cluster state and the relation between observed X-ray properties. In order to reach this goal, we apply the X-ray virtual telescope PHOX to generate synthetic observations of the simulated clusters with both Chandra and ATHENA, the latter used as an example for the performance of very high-resolution X-ray telescopes. From Chandra spectra we extract global properties, e.g. luminosity and temperature, and from ATHENA spectra we estimate the gas velocity dispersion along the line of sight from the broadening of heavy-ion emission lines (e.g. Fe). We further extend the analysis to the relation between non-thermal velocity dispersion of the gas and the L_X-T scaling law for the simulated clusters. Interestingly, we find a clear dependence of slope and scatter on the selection criterion for the clusters, based on the level of significance of non-thermal motions. Namely, the scatter in the relation is significantly reduced by the exclusion of the clusters, for which we estimate the highest turbulent velocities. Such velocity diagnostics appears therefore as a promising independent way to identify disturbed clusters, in addition to the commonly used morphological inspection., Comment: 17 pages, 15 figures, abstract abridged. Accepted for publication in MNRAS

We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, warm (10^4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked via the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new IFU measurements in Perseus cluster. The ensemble multiphase gas distributions are characterized by substantial spectral line broadening (100-200 km/s) with mild line shift. On the other hand, pencil-beam detections sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km/s, with increased scatter due to the turbulence intermittency. We present new ensemble sigma_v of the warm Halpha+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). We show the physically motivated criterion C = t_cool/t_eddy ~ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters/groups. The ensemble method can be applied to many available datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions., 18 pages, 5 figures, 2 tables; ApJ in press - typos corrected, formatting improved (new emulateapj mod included here)

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., Comment: 13 pages, 10 figures, Accepted by MNRAS 2018 June 17