Rotation in galaxy clusters from MUSIC simulations with the kinetic Sunyaev-Zel'dovich effect

We propose in this work its application for the detection of possible coherent rotational motions in the hot intra-cluster medium. We select a sample of massive, relaxed and rotating galaxy clusters from Marenostrum-mUltidark SImulations of galaxy Clusters (MUSIC), and we produce mock maps of the temperature distortion produced by the kinetic Sunyaev-Zel'dovich effect by exploring six different lines of sight, in the best observational condition. These maps are compared with the expected signal computed from a suitable theoretical model in two cases: (i) focusing only on the contribution from the rotation, and (ii) accounting also for the cluster bulk motion. We find that the parameters of the model assumed for the radial profile of the rotational velocity, averaged over the considered lines of sight, are in agreement within two standard deviations at most with independent estimates from the simulation data, without being significantly affected by the presence of the cluster bulk term. The amplitude of the rotational signal is, on average, of the order of 23 per cent of the total signal accounting also for the cluster bulk motion, and its values are consistent with the literature. The projected bulk velocity of the cluster is also recovered at the different lines of sight, with values in agreement with the simulation dataASB acknowledges funding from Sapienza Università di Roma - Progetti per Avvio alla Ricerca Anno 2017, prot. AR11715C82402BC

On the coherent rotation of diffuse matter in numerical simulations of galaxy clusters

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

Morphological estimators on Sunyaev-Zel'dovich maps of MUSIC clusters of galaxies

The determination of the morphology of galaxy clusters has important repercussions for cosmological and astrophysical studies of them. In this paper, we address the morphological characterization of synthetic maps of the Sunyaev-Zel'dovich (SZ) effect for a sample of 258 massive clusters (Mvir> 5×1014h-1M⊙at z=0), extracted from theMUSIC hydrodynamical simulations. Specifically, we use five known morphological parameters (which are already used in X-ray) and two newly introduced ones, and we combine them 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 red shifts between 0.43 and 0.82. For each parameter, we test its stability and efficiency in discriminating the true cluster dynamical state, measured by theoretical indicators. The combined parameter is more efficient at discriminating between relaxed and disturbed clusters. This parameter had a mild correlation with the hydrostatic mass (~0.3) and a strong correlation (~0.8) with the offset between the SZ centroid and the cluster centre of mass. The latter quantity is, thus, the most accessible and efficient indicator of the dynamical state for SZ studiesThis work has been partially supported by funding from Sapienza University of Rome - Progetti di Ricerca Anno 2015 prot. C26A15LXNR. GY and FS acknowledge financial support from MINECO/FEDER under research grant AYA2015-63810-P. ER acknowledge financial contribution from the agreement ASI-INAF n 2017-14-H.

Rotation in galaxy clusters from MUSIC simulations with the kinetic Sunyaev-Zel'dovich effect

The Sunyaev-Zel'dovich effect in galaxy clusters is a unique probe for studying astrophysics and cosmology. We propose in this work its application for the detection of possible coherent rotational motions in the hot intra-cluster medium. We select a sample of massive, relaxed and rotating galaxy clusters from Marenostrum-mUltidark SImulations of galaxy Clusters (MUSIC), and we produce mock maps of the temperature distortion produced by the kinetic Sunyaev-Zel'dovich effect by exploring six different lines of sight, in the best observational condition. These maps are compared with the expected signal computed from a suitable theoretical model in two cases: (i) focusing only on the contribution from the rotation, and (ii) accounting also for the cluster bulk motion. We find that the parameters of the model assumed for the radial profile of the rotational velocity, averaged over the considered lines of sight, are in agreement within two standard deviations at most with independent estimates from the simulation data, without being significantly affected by the presence of the cluster bulk term. The amplitude of the rotational signal is, on average, of the order of 23 per cent of the total signal accounting also for the cluster bulk motion, and its values are consistent with the literature. The projected bulk velocity of the cluster is also recovered at the different lines of sight, with values in agreement with the simulation data

Machine learning methods to estimate observational properties of galaxy clusters in large volume cosmological N-body simulations

This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record Monthly Notices of the Royal Astronomical Society 518.1 (2023): 111-129 is available online at: https://academic.oup.com/mnras/article-abstract/518/1/111/6795309?redirectedFrom=fulltext#no-access-messageIn 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 (DM)-only cluster-size haloes. The training set is built from the three hundred project that consists of a series of zoomed hydrodynamical simulations of cluster-size regions extracted from the 1 Gpc volume MultiDark DM-only simulation (MDPL2). We use as target variables a set of baryonic properties for the intracluster gas and stars derived from the hydrodynamical simulations and correlate them with the properties of the DM haloes from the MDPL2 N-body simulation. The different ML models are trained from this data base and subsequently used to infer the same baryonic properties for the whole range of cluster-size haloes identified in the MDPL2. We also test the robustness of the predictions of the models against mass resolution of the DM haloes and conclude that their inferred baryonic properties are rather insensitive to their DM properties that 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 haloes 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 availableThe authors thank the anonymous referee for his/her invaluable comments and suggestions, without which this work would be in complete. D.d.A., W.C. and G.Y. would like to thank Ministerio de Ciencia e Innovación for financial support under project grant PID2021-122603NB-C21. WC is supported by the STFC AGP Grant ST/V000594/1 and the Atracción de Talento Contract no. 2020-T1/TIC-19882 granted by the Comunidad de Madrid in Spain. He also thanks the Ministerio de Ciencia e Innovación (Spain) for financial support under Project grant PID2021-122603NB C21. He further acknowledges the science research grants from the China Manned Space Project with NO. CMS-CSST-2021-A01 and CMS-CSST-2021-B01. G.M. acknowledges financial support from PID2019-106827GB-I00/AEI / 10.13039/501100011033 The CosmoSim database used in this paper is a service by the Leibniz Institute for Astrophysics Potsdam (AIP). The MultiDark database was developed in cooperation with the Spanish MultiDark Con solider Project CSD2009-00064. The authors acknowledge The Red Española de Supercomputación for granting computing time for running the hydrodynamical simulations of The300 galaxy cluster project in the Marenostrum supercomputer at the Barcelona Super-computing Cente

The MUSIC of Galaxy Clusters II: X-ray global properties and scaling relations

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) data set. We focus on a sub-sample of 179 clusters at redshift z similar to 0.11, with 3.2 x 10(14) h(-1) M-circle dot < M-vir < 2 x 10(15) h(-1) M-circle dot, complete in mass. We employed the X-ray photon simulator PHOX to obtain synthetic Chandra observations and derive observable-like global properties of the intracluster medium (ICM), as X-ray temperature (T-X) and luminosity (L-X). T-X is found to slightly underestimate the true mass-weighted temperature, although tracing fairly well the cluster total mass. We also study the effects of T-X on scaling relations with cluster intrinsic properties: total (M-500 and gas M-g,M-500 mass; integrated Compton parameter (Y-SZ) of the Sunyaev-Zel'dovich (SZ) thermal effect; Y-X = M-g,M-500 T-X. We confirm that Y-X is a very good mass proxy, with a scatter on M-500-Y-X and Y-SZ-Y-X lower than 5 per cent. The study of scaling relations among X-ray, intrinsic and SZ properties indicates that simulated MUSIC clusters reasonably resemble the self-similar prediction, especially for correlations involving T-X. The observational approach also allows for a more direct comparison with real clusters, from which we find deviations mainly due to the physical description of the ICM, affecting T-X and, particularly, L-X

Impact of baryons on the cluster mass function and cosmological parameter determination

Recent results by the Planck collaboration have shown that cosmological parameters derived from the cosmic microwave background anisotropies and cluster number counts are in tension, with the latter preferring lower values of the matter density parameter, $\Omega_\mathrm{m}$, and power spectrum amplitude, $\sigma_8$. Motivated by this, we investigate the extent to which the tension may be ameliorated once the effect of baryonic depletion on the cluster mass function is taken into account. We use the large-volume Millennium Gas simulations in our study, including one where the gas is pre-heated at high redshift and one where the gas is heated by stars and active galactic nuclei (in the latter, the self-gravity of the baryons and radiative cooling are omitted). In both cases, the cluster baryon fractions are in reasonably good agreement with the data at low redshift, showing significant depletion of baryons with respect to the cosmic mean. As a result, it is found that the cluster abundance in these simulations is around 15 per cent lower than the commonly-adopted fit to dark matter simulations by Tinker et al (2008) for the mass range $10^{14}-10^{14.5}h^{-1} \mathrm{M}_\odot$. Ignoring this effect produces a significant artificial shift in cosmological parameters which can be expressed as $\Delta[\sigma_8(\Omega_\mathrm{m}/0.27)^{0.38}]\simeq -0.03$ at $z=0.17$ (the median redshift of the $\mathit{Planck}$ cluster sample) for the feedback model. While this shift is not sufficient to fully explain the $\mathit{Planck}$ discrepancy, it is clear that such an effect cannot be ignored in future precision measurements of cosmological parameters with clusters. Finally, we outline a simple, model-independent procedure that attempts to correct for the effect of baryonic depletion and show that it works if the baryon-dark matter back-reaction is negligible.Comment: 10 pages, 5 figures, Accepted by MNRA

The MUSIC of CLASH: predictions on the concentration-mass relation

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

nIFTY galaxy cluster simulations - III. The similarity and diversity of galaxies and subhaloes

We examine subhaloes and galaxies residing in a simulated $\Lambda$ cold dark matter galaxy cluster ($M^{crit} _{200}$ = 1.1 × 10$^{15}$ $h^{−1}$ M$_\odot$) produced by hydrodynamical codes ranging from classic smooth particle hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity-only, non-radiative hydrodynamics and $\textit{full feedback physics}$ runs by looking at the overall subhalo/galaxy distribution and on an individual object basis. We find that the subhalo population is reproduced to within $\lesssim$10 per cent for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim$0.1 dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by ≈20 per cent. The galaxy populations also show striking code-to-code variations. Although the Tully–Fisher relation is similar in almost all codes, the number of galaxies with 10$^9$ $h^{−1}$ M$_\odot$ $\lesssim$ $M_∗$ $\lesssim$ 10$^{12}$ $h^{−1}$ M$_\odot$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of ~0.5 dex in stellar mass. Moreover, systematic differences exist, with some codes producing galaxies 70 per cent smaller than others. The diversity partially arises from the inclusion/absence of active galactic nucleus feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies $\textit{in all environments}$ remains a challenge. We argue that even basic galaxy properties, such as stellar mass to halo mass, should be treated with errors bars of ~0.2–0.4 dex

nIFTY galaxy cluster simulations - III. The similarity and diversity of galaxies and subhaloes

We examine subhaloes and galaxies residing in a simulated $\Lambda$ cold dark matter galaxy cluster ($M^{crit} _{200}$ = 1.1 × 10$^{15}$ $h^{−1}$ M$_\odot$) produced by hydrodynamical codes ranging from classic smooth particle hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity-only, non-radiative hydrodynamics and $\textit{full feedback physics}$ runs by looking at the overall subhalo/galaxy distribution and on an individual object basis. We find that the subhalo population is reproduced to within $\lesssim$10 per cent for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim$0.1 dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by ≈20 per cent. The galaxy populations also show striking code-to-code variations. Although the Tully–Fisher relation is similar in almost all codes, the number of galaxies with 10$^9$ $h^{−1}$ M$_\odot$ $\lesssim$ $M_∗$ $\lesssim$ 10$^{12}$ $h^{−1}$ M$_\odot$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of ~0.5 dex in stellar mass. Moreover, systematic differences exist, with some codes producing galaxies 70 per cent smaller than others. The diversity partially arises from the inclusion/absence of active galactic nucleus feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies $\textit{in all environments}$ remains a challenge. We argue that even basic galaxy properties, such as stellar mass to halo mass, should be treated with errors bars of ~0.2–0.4 dex