Search

The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a programme to study a minimally biased sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. Accurate and precise mass measurements are required to exploit CHEX-MATE as an astrophysical laboratory and a calibration sample for cosmological probes in the era of large surveys. We measured masses based on the galaxy dynamics, which are highly complementary to weak-lensing or X-ray estimates. We analysed the sample with a uniform pipeline that is stable both for poorly sampled or rich clusters - using spectroscopic redshifts from public (NED, SDSS, and DESI) or private archives - and dedicated observational programmes. We modelled the halo mass density and the anisotropy profile. Membership is confirmed with a cleaning procedure in phase space. We derived masses from measured velocity dispersions under the assumed model. We measured dynamical masses for 101 CHEX-MATE clusters with at least ten confirmed members within the virial radius r_200c. Estimated redshifts and velocity dispersions agree with literature values when available. Validation with weak-lensing masses shows agreement within 8+-16(stat.)+-5(sys.)%, and confirms dynamical masses as an unbiased proxy. Comparison with {\it Planck} masses shows them to be biased low by 34+-3(stat.)+-5(sys.)%. A follow-up spectroscopic campaign is underway to cover the full CHEX-MATE sample., Comment: 18 pages, 12 figures. In press on A&A

The observation of a massive galaxy with an extremely low dark matter content (i.e. NGC 1277) has posed questions about how such objects form and evolve in a hierarchical universe. We here report on the finding of several massive, dark matter-deficient galaxies in a set of 324 galaxy clusters theoretically modelled by means of full-physics hydrodynamical simulations. We first focus on two example galaxies selected amongst the most massive and dark matter-deficient ones. By tracing the evolution of these galaxies, we find that their lack of dark matter is a result of multiple pericentre passages. While orbiting their host halo, tidal interactions gradually strip away dark matter while preserving the stellar component. A statistical analysis of all massive satellite galaxies in the simulated clusters shows that the stellar-to-total mass ratio today is strongly influenced by the number of orbits and the distance at pericentres. Galaxies with more orbits and closer pericentres are more dark matter-deficient. Additionally, we find that massive, dark matter-deficient galaxies at the present day are either the remnants of very massive galaxies at infall or former central galaxies of infalling groups. We conclude that such massive yet dark matter-deficient galaxies exist and are natural by-products of typical cluster galaxy evolution, with no specific requirement for an exotic formation scenario., Comment: 12 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

Hydrodynamical cosmological simulations are ideal laboratories where the evolution of the cosmic web can be studied. This allows for easier insight into the nature of the filaments. We investigate how the intrinsic properties of filaments are evolving in areas extracted from a larger cosmological simulation. We aim to identify significant trends in the properties of Warm-Hot Intergalactic Medium (WHIM) and suggest possible explanations. To study the filaments and their contents, we select a subset of regions from the Dianoga simulation. We analysed these regions that were simulated with different baryon physics, namely with and without the AGN feedback. We construct the cosmic web using the Sub-space Constrained Mean Shift (SCMS) algorithm and the Sequential Chain Algorithm for Resolving Filaments (SCARF). We examined the basic physical properties of filaments (length, shape, mass, radius) and analysed different gas phases (hot, WHIM and colder gas components) within those structures. The evolution of the global filament properties and the properties of the gas phases were studied in the redshift range $0 < z < 1.48$. Within our simulations, the detected filaments have, on average, lengths below $9$ Mpc. The filaments' shape correlates with their length; the longer they are, the more likely they are curved. We find that the scaling relation between mass $M$ and length $L$ of the filaments is well described by the power law $M \propto L^{1.7}$. The radial density profile is widening with redshift, meaning that the radius of the filaments is getting larger over time. The fraction of gas mass in the WHIM phase does not depend on the model and is rising towards lower redshifts. However, the included baryon physics has a strong impact on the metallicity of gas in filaments, indicating that the AGN feedback impacts the metal content already at redshifts of $z \sim 2$., Comment: 16 pages, 13 figures

The intra-cluster medium is prone to turbulent motion that will contribute to the non-thermal heating of the gas, complicating the use of galaxy clusters as cosmological probes. Indirect approaches can estimate the intensity and structure of turbulent motions by studying the associated fluctuations in gas density and X-ray surface brightness. In this work, we want to constrain the gas density fluctuations at work in the CHEX-MATE sample to obtain a detailed view of their properties in a large population of clusters. We use a simulation-based approach to constrain the parameters of the power spectrum of density fluctuations, assuming a Kolmogorov-like spectrum and including the sample variance, further providing an approximate likelihood for each cluster. This method requires clusters to be not too disturbed, as fluctuations can originate from dynamic processes such as merging. Accordingly, we remove the less relaxed clusters (centroid shift $w>0.02$) from our sample, resulting in a sample of 64 clusters. We define different subsets of CHEX-MATE to determine properties of density fluctuations as a function of dynamical state, mass and redshift, and investigate the correlation with the presence or not of a radio halo. We found a positive correlation between the dynamical state and density fluctuation variance, a non-trivial behaviour with mass and no specific trend with redshift or the presence/absence of a radio halo. The injection scale is mostly constrained by the core region. The slope in the inertial range is consistent with Kolmogorov theory. When interpreted as originating from turbulent motion, the density fluctuations in $R_{500}$ yield an average Mach number of $M_{3D}\simeq 0.4\pm 0.2$, an associated non-thermal pressure support of $ P_{turb}/P_{tot}\simeq (9\pm 6) \%$ or a hydrostatic mass bias $b_{turb}\simeq 0.09\pm 0.06$, in line with what is expected from the literature., Comment: Accepted for publication in A&A, abstract slightly abridged for Arxiv

We characterise the intra-cluster light (ICL) in ensembles of full-physics cluster simulations from The Three Hundred project, a suite of 324 hydrodynamical resimulations of cluster-sized halos. We identify the ICL as those stellar particles bound to the potential of the cluster itself, but not to any of its substructures, and separate the brightest cluster galaxy (BCG) by means of a fixed 50 kpc aperture. We find the total BCG+ICL mass to be in agreement with state-of-the-art observations of galaxy clusters. The ICL mass fraction of our clusters is between 30 and 50 per cent of the total stellar mass within $R_{500}$, while the BCG represents around 10 percent. We further find no trend of the ICL fraction with cluster halo mass, at least not in the range $[0.2,3]\cdot10^{15}h^{-1}M_\odot$ considered here. For the dynamical state, characterised both by theoretical estimators and by the recent merging history of the cluster, there is a clear correlation, such that more relaxed clusters and those that have undergone fewer recent mergers have a higher ICL fraction. Finally, we investigate the possibility of using the ICL to explore the dark matter (DM) component of galaxy clusters. We compute the volumetric density profile for the DM and ICL components and show that, up to $R_{500}$, the ratio between the two can be described by a power law. Working with the velocity dispersion profiles instead, we show that the ratio can be fit by a straight line. Providing the parameters of these fits, we show how the ICL can be used to infer DM properties., Comment: 18 pages, 14 figures. Accepted for publication in Astronomy & Astrophysics

Synthesized in the cores of stars and supernovae, most metals disperse over cosmic scales and are ultimately deposited well outside the gravitational potential of their host galaxies. Since their presence is well visible through their X-ray emission lines in the hot gas pervading galaxy clusters, measuring metal abundances in the intracluster medium (ICM) offers us a unique view of chemical enrichment of the Universe as a whole. Despite extraordinary progress in the field thanks to four decades of X-ray spectroscopy using CCD (and gratings) instruments, understanding the precise stellar origins of the bulk of metals, and when the latter were mixed on Mpc scales, requires an X-ray mission capable of spatial, non-dispersive high resolution spectroscopy covering at least the soft X-ray band over a large field of view. In this White Paper, we demonstrate how the Line Emission Mapper (LEM) probe mission concept will revolutionize our current picture of the ICM enrichment. Specifically, we show that LEM will be able to (i) spatially map the distribution of ten key chemical elements out to the virial radius of a nearby relaxed cluster and (ii) measure metal abundances in serendipitously discovered high-redshift protoclusters. Altogether, these key observables will allow us to constrain the chemical history of the largest gravitationally bound structures of the Universe. They will also solve key questions such as the universality of the initial mass function (IMF) and the initial metallicity of the stellar populations producing these metals, as well as the relative contribution of asymptotic giant branch (AGB) stars, core-collapse, and Type Ia supernovae to enrich the cosmic web over Mpc scales. Concrete observing strategies are also briefly discussed., Comment: 19 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity

High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field-of-view, and effective area are essential to separate cluster emissions from the Galactic foreground and efficiently map the cluster outskirts. Several mission concepts that meet these criteria have been proposed recently, e.g., LEM, HUBS, and SuperDIOS. This theoretical study explores what information on ICM physics could be recovered with such missions and the associated challenges. We emphasize the need for a comprehensive comparison between simulations and observations to interpret the high-resolution spectroscopic observations correctly. Using Line Emission Mapper (LEM) characteristics as an example, we demonstrate that it enables the use of soft X-ray emission lines (e.g., O VII/VIII and Fe-L complex) from the cluster outskirts to measure the thermodynamic, chemical, and kinematic properties of the gas up to $r_{200}$ and beyond. By generating mock observations with full backgrounds, analysing their images/spectra with observational approaches, and comparing the recovered characteristics with true ones from simulations, we develop six key science drivers for future missions, including the exploration of multiphase gas in galaxy clusters (e.g., temperature fluctuations, phase-space distributions), metallicity, ICM gas bulk motions and turbulence power spectra, ICM-cosmic filament interactions, and advances for cluster cosmology., Comment: 24 pages, 26 figures, submitted to MNRAS. Comments are welcome

Previous studies have revealed that the estimated probability of galaxy-galaxy strong lensing in observed galaxy clusters exceeds the expectations from the $\Lambda$ Cold Dark Matter cosmological model by one order of magnitude. We aim to understand the origin of this excess by analyzing a larger set of simulated galaxy clusters and investigating how the theoretical expectations vary under different adopted prescriptions and numerical implementations of star formation and feedback in simulations. We perform a ray-tracing analysis of 324 galaxy clusters from the Three Hundred project, comparing the Gadget-X and Gizmo-Simba runs. These simulations, which start from the same initial conditions, are performed with different implementations of hydrodynamics and galaxy formation models tailored to match different observational properties of the Intra-Cluster-Medium and cluster galaxies. We find that galaxies in the Gizmo-Simba simulations develop denser stellar cores than their Gadget-X counterparts. Consequently, their probability for galaxy-galaxy strong lensing is higher by a factor of $\sim 3$. This increment is still insufficient to fill the gap with observations, as a discrepancy by a factor $\sim 4$ still persists. In addition, we find that several simulated galaxies have Einstein radii that are too large compared to observations. We conclude that a persistent excess of galaxy-galaxy strong lensing exists in observed galaxy clusters. The origin of this discrepancy with theoretical predictions is still unexplained in the framework of the cosmological hydrodynamical simulations. This might signal a hitherto unknown issue with either the simulation methods or our assumptions regarding the standard cosmological model., Comment: 6 pages, 2 figures, accepted for publication on A&A Letters. See companion paper Srivastava et al. (2023)

The distribution of baryons provides a significant way to understand the formation of galaxy clusters by revealing the details of its internal structure and changes over time. In this paper, we present theoretical studies on the scaled profiles of physical properties associated with the baryonic components, including gas density, temperature, metallicity, pressure and entropy as well as stellar mass, metallicity and satellite galaxy number density in galaxy clusters from $z=4$ to $z=0$ by tracking their progenitors. These mass-complete simulated galaxy clusters are coming from THE THREE HUNDRED with two runs: GIZMO-SIMBA and Gadget-X. Through comparisons between the two simulations, and with observed profiles which are generally available at low redshift, we find that (1) the agreements between the two runs and observations are mostly at outer radii $r \gtrsim 0.3r_{500}$, in line with the self-similarity assumption. While Gadget-X shows better agreements with the observed gas profiles in the central regions compared to GIZMO-SIMBA; (2) the evolution trends are generally consistent between the two simulations with slightly better consistency at outer radii. In detail, the gas density profile shows less discrepancy than the temperature and entropy profiles at high redshift. The differences in the cluster centre and gas properties imply different behaviours of the AGN models between Gadget-X and GIZMO-SIMBA, with the latter, maybe too strong for this cluster simulation. The high-redshift difference may be caused by the star formation and feedback models or hydrodynamics treatment, which requires observation constraints and understanding., Comment: 20 pages, 20 figures, accepted in MNRAS

The Sunyaev Zel'dovich (SZ) effect is sensitive to the pressure of ionized gas inside galaxy clusters, which is in turn controlled largely by the gravitational potential of the cluster. Changing the concentration parameter describing the cluster mass distribution impacts the gravitational potential and thus the cluster SZ signal, with implications for cosmological and other analyses of SZ-selected clusters. We investigate the concentration-SZ relation in theory and simulations. We find that the impact of concentration on the inner SZ profile ($R \lesssim 0.75 R_{200c}$) can be captured with standard polytropic gas models. However, we find that such models do a poor job of reproducing the outer SZ profiles ($R \gtrsim 0.75 R_{200c}$) and the relation between the integrated SZ signal, $Y$, and concentration. This disagreement results from a sharp truncation of the gas pressure profile near the splashback radius, likely caused by virial shocks. We develop a simple description of the truncation that leads to a good match with the simulated SZ profiles out to several $R_{200c}$ for clusters of varying mass and concentration, and that also accurately predicts the concentration-$Y$ relationship. Finally, we determine how inference of the linear bias parameter and splashback radius for SZ-selected clusters can be biased by ignoring the concentration dependence of the SZ signal, finding that bias to the former is essentially negligible, while bias to the latter can be as much as 2\%., Comment: 16 pages, 11 figures; replaced to match version accepted by MNRAS. Revised version uses an updated model for nonthermal pressure support, but conclusions are not impacted

Galaxy clusters are the most massive gravitationally bound structures in the Universe, comprising thousands of galaxies and pervaded by a diffuse, hot ``intracluster medium'' (ICM) that dominates the baryonic content of these systems. The formation and evolution of the ICM across cosmic time is thought to be driven by the continuous accretion of matter from the large-scale filamentary surroundings and dramatic merger events with other clusters or groups. Until now, however, direct observations of the intracluster gas have been limited only to mature clusters in the latter three-quarters of the history of the Universe, and we have been lacking a direct view of the hot, thermalized cluster atmosphere at the epoch when the first massive clusters formed. Here we report the detection (about $6\sigma$) of the thermal Sunyaev-Zeldovich (SZ) effect in the direction of a protocluster. In fact, the SZ signal reveals the ICM thermal energy in a way that is insensitive to cosmological dimming, making it ideal for tracing the thermal history of cosmic structures. This result indicates the presence of a nascent ICM within the Spiderweb protocluster at redshift $z=2.156$, around 10 billion years ago. The amplitude and morphology of the detected signal show that the SZ effect from the protocluster is lower than expected from dynamical considerations and comparable with that of lower-redshift group-scale systems, consistent with expectations for a dynamically active progenitor of a local galaxy cluster., Comment: 18 pages (main text + methods), 10 figures, 2 tables; published online in Nature on March 29th, 2023

The relationship linking a galaxy cluster's total mass with the concentration of its mass profile and its redshift is a fundamental prediction of the Cold Dark Matter (CDM) paradigm of cosmic structure formation. However, confronting those predictions with observations is complicated by the fact that simulated clusters are not representative of observed samples where detailed mass profile constraints are possible. In this work, we calculate the Symmetry-Peakiness-Alignment (SPA) morphology metrics for maps of X-ray emissivity from THE THREE HUNDRED project hydrodynamical simulations of galaxy clusters at four redshifts, and thereby select a sample of morphologically relaxed, simulated clusters, using observational criteria. These clusters have on average earlier formation times than the full sample, confirming that they are both morphologically and dynamically more relaxed than typical. We constrain the concentration-mass-redshift relation of both the relaxed and complete sample of simulated clusters, assuming power-law dependences on mass ($\kappa_m$) and $1+z$ ($\kappa_\zeta$), finding $\kappa_m = -0.12 \pm 0.07$ and $\kappa_\zeta = -0.27 \pm 0.19$ for the relaxed subsample. From an equivalently selected sample of massive, relaxed clusters observed with ${\it Chandra}$, we find $\kappa_m = -0.12 \pm 0.08$ and $\kappa_\zeta = -0.48 \pm 0.19$, in good agreement with the simulation predictions. The simulated and observed samples also agree well on the average concentration at a pivot mass and redshift providing further validation of the $\Lambda$CDM paradigm in the properties of the largest gravitationally collapsed structures observed. This also represents the first clear detection of decreasing concentration with redshift, a longstanding prediction of simulations, in data., Comment: 10 pages, 10 figures. Accepted by MNRAS

The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole feedback and flows of baryonic matter into and out of galaxies. These processes are best studied in X-rays, and emission-line mapping is the pressing need in this area. LEM will use a large microcalorimeter array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft X-ray line emission from objects that constitute galactic ecosystems. These include supernova remnants, star-forming regions, superbubbles, galactic outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs in other galaxies), the Circumgalactic Medium in the Milky Way and other galaxies, and the Intergalactic Medium at the outskirts and beyond the confines of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV band will make it possible to disentangle the faintest emission lines in those objects from the bright Milky Way foreground, providing groundbreaking measurements of the physics of these plasmas, from temperatures, densities, chemical composition to gas dynamics. While LEM's main focus is on galaxy formation, it will provide transformative capability for all classes of astrophysical objects, from the Earth's magnetosphere, planets and comets to the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling gas in galaxy clusters. In addition to pointed observations, LEM will perform a shallow all-sky survey that will dramatically expand the discovery space., Comment: 18 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. v2: All-sky survey figure expanded, references fixed. v3: Added energy resolution measurements for prototype detector array. v4: Author list and reference fixes

We use a set of about 300 simulated clusters from The Three Hundred Project to calculate their hydrostatic masses and evaluate the associated bias by comparing them with the true cluster mass. Over a redshift range from 0.07 to 1.3, we study the dependence of the hydrostatic bias on redshift, concentration, mass growth, dynamical state, mass, and halo shapes. We find almost no correlation between the bias and any of these parameters. However, there is a clear evidence that the scatter of the mass-bias distribution is larger for low-concentrated objects, high mass growth, and more generically for disturbed systems. Moreover, we carefully study the evolution of the bias of twelve clusters throughout a major-merger event. We find that the hydrostatic-mass bias follows a particular evolution track along the merger process: to an initial significant increase of the bias recorded at the begin of merger, a constant plateaus follows until the end of merge, when there is a dramatic decrease in the bias before the cluster finally become relaxed again. This large variation of the bias is in agreement with the large scatter of the hydrostatic bias for dynamical disturbed clusters. These objects should be avoided in cosmological studies because their exact relaxation phase is difficult to predict, hence their mass bias cannot be trivially accounted for., Comment: 11 pages, 8 figures. Accepted for publication in MNRAS

The reservoir of molecular gas (H$_{\rm 2}$) represents the fuel for the star formation (SF) of a galaxy. Connecting the star formation rate (SFR) to the available H$_{\rm 2}$ is key to accurately model SF in cosmological simulations of galaxy formation. We investigate how modifying the underlying modelling of H$_{\rm 2}$ and the description of stellar feedback in low-metallicity environments (LMF, i.e. low-metallicity stellar feedback) in cosmological, zoomed-in simulations of a Milky Way-size halo influences the formation history of the forming, spiral galaxy and its final properties. We exploit two different models to compute the molecular fraction of cold gas (f$_{\rm H_{\rm 2}}$): $i)$ the theoretical model by Krumholz et al. (2009b) and $ii)$ the phenomenological prescription by Blitz & Rosolowsky (2006). We find that the model adopted to estimate f$_{\rm H_{\rm 2}}$ plays a key role in determining final properties and in shaping the morphology of the galaxy. The clumpier interstellar medium (ISM) and the more complex H$_{\rm 2}$ distribution that the Krumholz et al. (2009b) model predicts result in better agreement with observations of nearby disc galaxies. This shows how crucial it is to link the SFR to the physical properties of the star-forming, molecular ISM. The additional source of energy that LMF supplies in a metal-poor ISM is key in controlling SF at high redshift and in regulating the reservoir of SF across cosmic time. Not only is LMF able to regulate cooling properties of the ISM, but it also reduces the stellar mass of the galaxy bulge. These findings can foster the improvement of the numerical modelling of SF in cosmological simulations., Comment: accepted for publication in MNRAS

We have combined X-ray observations from Chandra with Sunyaev-Zel'dovich (SZ) effect data from Planck and Bolocam to measure intra-cluster medium pressure profiles from 0.03R$_{500}$ $\le$ R $\le$ 5R$_{500}$ for a sample of 21 low-$z$ galaxy clusters with a median redshift $\langle z \rangle = 0.08$ and a median mass $\langle \textrm{M}_{500} \rangle = 6.1 \times 10^{14}$ M$_{\odot}$ and a sample of 19 mid-$z$ galaxy clusters with $\langle z \rangle = 0.50$ and $\langle \textrm{M}_{500} \rangle = 10.6 \times 10^{14}$ M$_{\odot}$. The mean scaled pressure in the low-$z$ sample is lower at small radii and higher at large radii, a trend that is accurately reproduced in similarly selected samples from The300 simulations. This difference appears to be primarily due to dynamical state at small radii, evolution at intermediate radii, and a combination of evolution and mass dependence at large radii. Furthermore, the overall flattening of the mean scaled pressure profile in the low-$z$ sample compared to the mid-$z$ sample is consistent with expectations due to differences in mass accretion rate and the fractional impact of feedback mechanisms. In agreement with previous studies, the fractional scatter about the mean scaled pressure profile reaches a minimum of $\simeq 20$ per cent near 0.5R$_{500}$. This scatter is consistent between the low-$z$ and mid-$z$ samples at all radii, suggesting it is not strongly impacted by sample selection, and this general behavior is reproduced in The300 simulations. Finally, analytic functions that approximately describe the mass and redshift trends in mean pressure profile shape are provided., Comment: Accepted for publication in ApJ. Updated to closely match the published version

In this work, we performed an analysis of the X-ray morphology of the 118 CHEX-MATE (Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation) galaxy clusters, with the aim to provide a classification of their dynamical state. To investigate the link between the X-ray appearance and the dynamical state, we considered four morphological parameters: the surface brightness concentration, the centroid shift, and the second- and third-order power ratios. These indicators result to be: strongly correlated with each other, powerful in identifying the disturbed and relaxed population, characterised by a unimodal distribution and not strongly influenced by systematic uncertainties. In order to obtain a continuous classification of the CHEX-MATE objects, we combined these four parameters in a single quantity, M, which represents the grade of relaxation of a system. On the basis of the M value, we identified the most extreme systems of the sample, finding 15 very relaxed and 27 very disturbed galaxy clusters. From a comparison with previous analysis on X-ray selected samples, we confirmed that the Sunyaev-Zeldovich (SZ) clusters tend to be more disturbed. Finally, by applying our analysis on a simulated sample, we found a general agreement between the observed and simulated results, with the only exception of the concentration. This latter behaviour, is partially related to the presence of particles with high smoothed-particle hydrodynamics density in the central regions of the simulated clusters due to the action of the idealised isotropic thermal Active Galactic Nuclei (AGN) feedback., Comment: 24 pages, 16 figures, 10 tables. Accepted for publication in A&A

In this paper we study the applicability of a set of supervised machine learning (ML) models specifically trained to infer observed related properties of the baryonic component (stars and gas) from a set of features of dark matter only cluster-size halos. The training set is built from THE THREE HUNDRED project which consists of a series of zoomed hydrodynamical simulations of cluster-size regions extracted from the 1 Gpc volume MultiDark dark-matter only simulation (MDPL2). We use as target variables a set of baryonic properties for the intra cluster gas and stars derived from the hydrodynamical simulations and correlate them with the properties of the dark matter halos from the MDPL2 N-body simulation. The different ML models are trained from this database and subsequently used to infer the same baryonic properties for the whole range of cluster-size halos identified in the MDPL2. We also test the robustness of the predictions of the models against mass resolution of the dark matter halos and conclude that their inferred baryonic properties are rather insensitive to their DM properties which are resolved with almost an order of magnitude smaller number of particles. We conclude that the ML models presented in this paper can be used as an accurate and computationally efficient tool for populating cluster-size halos with observational related baryonic properties in large volume N-body simulations making them more valuable for comparison with full sky galaxy cluster surveys at different wavelengths. We make the best ML trained model publicly available., Comment: 20 pages, 10 figures, published in MNRAS

In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm., Comment: 11 pages, 10 figures, abstract is redacted to fit arXiv character count limit

Meneghetti et al. (2020) recently reported an excess of galaxy-galaxy strong lensing (GGSL) in galaxy clusters compared to expectations from the LCDM cosmological model. Theoretical estimates of the GGSL probability are based on the analysis of numerical hydrodynamical simulations in the LCDM cosmology. We quantify the impact of the numerical resolution and AGN feedback scheme adopted in cosmological simulations on the predicted GGSL probability and determine if varying these simulation properties can alleviate the gap with observations. We repeat the analysis of Meneghetti et al. (2020) on cluster-size halos simulated with different mass and force resolutions and implementing several independent AGN feedback schemes. We find that improving the mass resolution by a factor of ten and twenty-five, while using the same galaxy formation model that includes AGN feedback, does not affect the GGSL probability. We find similar results regarding the choice of gravitational softening. On the contrary, adopting an AGN feedback scheme that is less efficient at suppressing gas cooling and star formation leads to an increase in the GGSL probability by a factor between three and six. However, we notice that such simulations form overly massive subhalos whose contribution to the lensing cross-section would be significant while their Einstein radii are too large to be consistent with the observations. The primary contributors to the observed GGSL cross-sections are subhalos with smaller masses, that are compact enough to become critical for lensing. The population with these required characteristics appears to be absent in simulations., Comment: 13 pages, 11 figures. Submitted for publication on Astronomy and Astrophysics

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments., Comment: 20 pages, 12 figures, matched to published version in MNRAS

Several methods are used to evaluate, from observational data, the dynamical state of galaxy clusters. Among them, the morphological analysis of cluster images is well suited for this purpose. We report a new approach to the morphology, which consists in analytically modelling the images with a set of orthogonal functions, the Zernike polynomials (ZPs). We validated the method on mock high-resolution Compton parameter maps of synthetic galaxy clusters from THE THREE HUNDRED project. To classify the maps for their morphology we defined a single parameter, $\mathcal{C}$, by combining the contribution of some ZPs in the modelling. We verify that $\mathcal{C}$ is linearly correlated with a combination of common morphological parameters and also with a proper 3D dynamical-state indicator available for the synthetic clusters we used. We also show the early results of the Zernike modelling applied on Compton parameter maps of local clusters ($z < 0.1$) observed by the $\textit{Planck}$ satellite. At last, we report the preliminary results of this kind of morphological analysis on mock X-ray maps of THE THREE HUNDRED clusters., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

Hydrostatic equilibrium (HE) is often used in observations to estimate galaxy clusters masses. We use a set of almost 300 simulated clusters from The Three Hundred Project, to estimate the cluster HE mass and the bias deriving from it. We study the dependence of the bias on several dynamical state indicators across a redshift range from 0.07 to 1.3, finding no dependence between them. Moreover, we focus our attention on the evolution of the HE bias during the merger phase, where the bias even reaches negative values due to an overestimation of the mass with HE., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available., Comment: [v1]: 8 figures, 14 pages in main text + 5 figures, 3 pages in Appendix. [v2]: Minor text edits. Accepted in MNRAS

Galaxy clusters have a triaxial matter distribution. The weak-lensing signal, an important part in cosmological studies, measures the projected mass of all matter along the line-of-sight, and therefore changes with the orientation of the cluster. Studies suggest that the shape of the brightest cluster galaxy (BCG) in the centre of the cluster traces the underlying halo shape, enabling a method to account for projection effects. We use 324 simulated clusters at four redshifts between 0.1 and 0.6 from `The Three Hundred Project' to quantify correlations between the orientation and shape of the BCG and the halo. We find that haloes and their embedded BCGs are aligned, with an average $\sim$20 degree angle between their major axes. The bias in weak lensing cluster mass estimates correlates with the orientation of both the halo and the BCG. Mimicking observations, we compute the projected shape of the BCG, as a measure of the BCG orientation, and find that it is most strongly correlated to the weak-lensing mass for relaxed clusters. We also test a 2-dimensional cluster relaxation proxy measured from BCG mass isocontours. The concentration of stellar mass in the projected BCG core compared to the total stellar mass provides an alternative proxy for the BCG orientation. We find that the concentration does not correlate to the weak-lensing mass bias, but does correlate with the true halo mass. These results indicate that the BCG shape and orientation for large samples of relaxed clusters can provide information to improve weak-lensing mass estimates., Comment: 15 pages, 9 figures, Figure 7 is key plot, Updated to match version accepted by MNRAS

We perform a systematic study of the recently discovered fundamental plane of galaxy clusters (CFP) using ~250 simulated clusters from The Three Hundred project, focusing on the stability of the plane against different temperature definitions and its dependence on the dynamical relaxation state of clusters. The CFP is characterised by $T\propto M_\mathrm{s}^\alpha\,r_\mathrm{s}^\beta$, defined with the gas temperature ($T$) and the characteristic halo scale radius and mass ($r_\mathrm{s}$ and $M_\mathrm{s}$) assuming a Navarro-Frenk-White halo description. We explore two definitions of weighted temperatures, namely mass-weighted and spectroscopic-like temperatures, in three radial ranges. The Three Hundred project clusters at $z=0$ lie on a thin plane whose parameters ($\alpha,\beta$) and dispersion ($0.015-0.030$ dex) depend on the gas temperature definition. The CFP for mass-weighted temperatures is closer to the virial equilibrium expectation ($\alpha=1,\beta=-1$) with a smaller dispersion. For gas temperatures measured within $500h^{-1}$ kpc, the resulting CFP deviates the most from the virial expectation and shifts towards the similarity solution for a secondary infall model ($\alpha=1.5,\beta=-2$). Independently of the temperature definition, we find that clusters at $z=1$ and relaxed clusters form a CFP similar to the virial expectation, unlike disturbed clusters exhibiting stronger evolution. Only systems formed over the last 4 Gyr present a CFP that is closer to the similarity solution. All these findings are compatible with the CFP obtained for a CLASH subsample excluding the hottest clusters with $T_\mathrm{X}>12$ keV., Comment: 19 pages, 10 figures, 9 tables. Accepted for publication in MNRAS

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

Inferring line-of-sight distances from redshifts in and around galaxy clusters is complicated by peculiar velocities, a phenomenon known as the "Fingers of God" (FoG). This presents a significant challenge for finding filaments in large observational data sets as these artificial elongations can be wrongly identified as cosmic web filaments by extraction algorithms. Upcoming targeted wide-field spectroscopic surveys of galaxy clusters and their infall regions such as the WEAVE Wide-Field Cluster Survey motivate our investigation of the impact of FoG on finding filaments connected to clusters. Using zoom-in resimulations of 324 massive galaxy clusters and their outskirts from The ThreeHundred project, we test methods typically applied to large-scale spectroscopic data sets. This paper describes our investigation of whether a statistical compression of the FoG of cluster centres and galaxy groups can lead to correct filament extractions in the cluster outskirts. We find that within 5 R200 (~15 Mpc/h) statistically correcting for FoG elongations of virialized regions does not achieve reliable filament networks compared to reference filament networks based on true positions. This is due to the complex flowing motions of galaxies towards filaments in addition to the cluster infall, which overwhelm the signal of the filaments relative to the volume we probe. While information from spectroscopic redshifts is still important to isolate the cluster regions, and thereby reduce background and foreground interlopers, we expect future spectroscopic surveys of galaxy cluster outskirts to rely on 2D positions of galaxies to extract cosmic filaments., Comment: 12 pages, 9 figures

We report the non-thermal pressure fraction (Pnt/Ptot) obtained from a three-dimensional triaxial analysis of 16 galaxy clusters in the CLASH sample using gravitational lensing (GL) data primarily from Subaru and HST, X-ray spectroscopic imaging from Chandra, and Sunyaev-Zel'dovich effect (SZE) data from Planck and Bolocam. Our results span the approximate radial range 0.015-0.4R200m (35-1000 kpc). At cluster-centric radii smaller than 0.1R200m the ensemble average Pnt/Ptot is consistent with zero with an upper limit of nine per cent, indicating that heating from active galactic nuclei and other relevant processes does not produce significant deviations from hydrostatic equilibrium (HSE). The ensemble average Pnt/Ptot increases outside of this radius to approximately 20 per cent at 0.4R200m, as expected from simulations, due to newly accreted material thermalizing via a series of shocks. Also in agreement with simulations, we find significant cluster-to-cluster variation in Pnt/Ptot and little difference in the ensemble average Pnt/Ptot based on dynamical state. We conclude that on average, even for diverse samples, HSE-derived masses in the very central regions of galaxy clusters require only modest corrections due to non-thermal motions., Comment: published in MNRAS

Galaxy cluster outskirts are described by complex velocity fields induced by diffuse material collapsing towards filaments, gas and galaxies falling into clusters, and gas shock processes triggered by substructures. A simple scenario that describes the large-scale tidal fields of the cosmic web is not able to fully account for this variety, nor for the differences between gas and collisionless dark matter. We have studied the filamentary structure in zoom-in resimulations centred on 324 clusters from The ThreeHundred project, focusing on differences between dark and baryonic matter. This paper describes the properties of filaments around clusters out to five $R_{200}$, based on the diffuse filament medium where haloes had been removed. For this, we stack the remaining particles of all simulated volumes to calculate the average profiles of dark matter and gas filaments. We find that filaments increase their thickness closer to nodes and detect signatures of gas turbulence at a distance of $\sim 2 h^{-1}\rm{Mpc}$ from the cluster. These are absent in dark matter. Both gas and dark matter collapse towards filament spines at a rate of $\sim 200 h^{-1} \rm{km ~ s^{-1}} $. We see that gas preferentially enters the cluster as part of filaments, and leaves the cluster centre outside filaments. We further see evidence for an accretion shock just outside the cluster. For dark matter, this preference is less obvious. We argue that this difference is related to the turbulent environment. This indicates that filaments act as highways to fuel the inner regions of clusters with gas and galaxies., Comment: 14 pages, 12 figures, accepted for publication in MNRAS

We study the connection between morphology and dynamical state of the simulated galaxy clusters in $z\in[0,1.031]$ from THE THREE HUNDRED Project. We quantify cluster dynamical state using a combination of dynamical indicators from theoretical measures and compare this combined parameter, $\chi$, with the results from morphological classifications. The dynamical state of the cluster sample shows a continuous distribution from dynamically relaxed, more abundant at lower redshift, to hybrid and disturbed. The dynamical state presents a clear dependence on the radius, with internal regions more relaxed than outskirts. The morphology from multi-wavelength mock observation of clusters in X-ray, optical, and Sunyaev-Zel'dovich (SZ) effect images, is quantified by $M$ -- a combination of six parameters for X-ray and SZ maps and the offsets between the optical position of the Brightest Central Galaxy (BCG) and the X-ray/SZ centroids. All the morphological parameters are highly correlated with each other, while they show a moderately strong correlation with the dynamical $\chi$ parameter. The X-ray or SZ peaks are less affected by the dynamical state than centroids, which results in reliable tracers of the cluster density peak. The principal source of contamination in the relaxed cluster fraction, inferred from morphological parameters, is due to dynamically hybrid clusters. Compared to individual parameters, which consider only one aspect of cluster property (e.g. only clumping or asymmetry), the combined morphological and dynamical parameters ($M$ and $\chi$) collect more information and provide a single and more accurate estimation of the cluster dynamical state., Comment: 18 pages, 12 figures. Accepted for publication in MNRAS

We explore the persistence of the alignment of brightest cluster galaxies (BCGs) with their local environment. We find that a significant fraction of BCGs do not coincide with the centroid of the X-ray gas distribution and/or show peculiar velocities (they are not at rest with respect to the cluster mean). Despite this, we find that BCGs are generally aligned with the cluster mass distribution even when they have significant offsets from the X-ray centre and significant peculiar velocities. The large offsets are not consistent with simple theoretical models. To account for these observations BCGs must undergo mergers preferentially along their major axis, the main infall direction. Such BCGs may be oscillating within the cluster potential after having been displaced by mergers or collisions, or the dark matter halo itself may not yet be relaxed., Comment: In press, MNRAS

Clusters of galaxies are useful tools to constrain cosmological parameters, only if their masses can be correctly inferred from observations. In particular, X-ray and Sunyaev-Zeldovich (SZ) effect observations can be used to derive masses within the framework of the hydrostatic equilibrium. Therefore, it is crucial to have a good control of the possible mass biases that can be introduced when this hypothesis is not valid. In this work, we analyzed a set of 260 synthetic clusters from the MUSIC simulation project, at redshifts $0 \leq z \leq 0.82$. We estimate the hydrostatic mass of the MUSIC clusters from X-ray only (temperature and density) and from X-ray and SZ (density and pressure). Then, we compare them with the true 3D dynamical mass. The biases are of the order of 20%. We find that using the temperature instead of the pressure leads to a smaller bias, although the two values are compatible within 1$\sigma$. Non-thermal contributions to the total pressure support, arising from bulk motion and turbulence of the gas, are also computed and show that they are sufficient to account for this bias. We also present a study of the correlation between the mass bias and the dynamical state of the clusters. A clear correlation is shown between the relaxation state of the clusters and the bias factor. We applied the same analysis on a subsample of 32 objects, already selected for supporting the NIKA2 SZ Large Program., Comment: 19 pages, 16 figures. Accepted for publication in MNRAS

The knowledge of the dynamical state of galaxy clusters allows to alleviate systematics when observational data from these objects are applied in cosmological studies. Evidence of correlation between the state and the morphology of the clusters is well studied. The morphology can be inferred by images of the surface brightness in the X-ray band and of the thermal component of the Sunyaev-Zel'dovich (tSZ) effect in the millimetre range. For this purpose, we apply, for the first time, the Zernike polynomial decomposition, a common analytic approach mostly used in adaptive optics to recover aberrated radiation wavefronts at the telescopes pupil plane. With this novel way we expect to correctly infer the morphology of clusters and so possibly, their dynamical state. To verify the reliability of this new approach we use more than 300 synthetic clusters selected in THE THREE HUNDRED project at different redshifts ranging from 0 up to 1.03. Mock maps of the tSZ, quantified with the Compton parameter, $y$-maps, are modelled with Zernike polynomials inside $R_{500}$, the cluster reference radius. We verify that it is possible to discriminate the morphology of each cluster by estimating the contribution of the different polynomials to the fit of the map. The results of this new method are correlated with those of a previous analysis made on the same catalogue, using two parameters that combine either morphological or dynamical-state probes. We underline that instrumental angular resolution of the maps has an impact mainly when we extend this approach to high-redshift clusters., Comment: 16 pages, 10 figures, 6 tables (including appendix A and B). Accepted for publication in MNRAS

Cold dark matter (CDM) constitutes most of the matter in the Universe. The interplay between dark and luminous matter in dense cosmic environments like galaxy clusters is studied theoretically using cosmological simulations. Observed gravitational lensing is used to test and characterize the properties of substructures - the small-scale distribution of dark matter - in clusters. An apt metric, the probability of strong lensing events produced by dark matter substructure, is devised and computed for 11 galaxy clusters. We report that observed cluster substructures are more efficient lenses than predicted by CDM simulations, by more than an order of magnitude. We suggest that hitherto undiagnosed systematic issues with simulations or incorrect assumptions about the properties of dark matter could explain our results., Comment: 56 pages, 15 figures. Published in Science

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

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., Comment: 20 pages, 18 figures

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., Comment: 10 pages, 8 figures. MNRAS accepted. Comments welcome

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

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., Comment: 17 pages, 12 figures, accepted to MNRAS

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., Comment: Invited review article, to be published in "From Dark Haloes to Visible Galaxies", special issue of Galaxies

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}., Comment: 15 pages, 9 figures, accepted for publication by A&A

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., Comment: 20 pages, 8 figures, 7 tables. MNRAS published version

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., Comment: Published on ApJ. Matched to published version

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., Comment: MNRAS, accepted

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., Comment: MNRAS, accepted

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., Comment: Replaced with a revised version to match the ApJ accepted version

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., Comment: 42 pages, 52 figures

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., Comment: Updated to match the MNRAS published version. 15 pages, 9 figures, 9 tables

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., Comment: 18 pages, 6 figures, 5 tables. Accepted for publication in ApJ