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

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

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

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

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 – V. Investigation of the cluster infall region

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$ cold dark matter galaxy cluster studied by Elahi et al. at $z$ = 0. The 1.1 × 10$^{15}$ $h^{−1}$ M$_\odot$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and sub-grid schemes. All models completed a dark-matter-only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_\text{DM}$ = 9.01 × 10$^8$ $h^{−1}$ M$_\odot$ and $m_\text{gas}$ = 1.9 × 10$^8$ $h^{−1}$ M$_\odot$, respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60 per cent of haloes in the infall region with mass ~10$^{12.5}$–10$^{14}$ $h^{−1}$ M$_\odot$, including 2–3 group-sized haloes (>10$^{13}$ $h^{−1}$ M$_\odot$). However, we find that only ~10 per cent of objects in the infall region are sub-haloes residing in haloes, which may suggest that there is not much ongoing pre-processing occurring in the infall region at $z$ = 0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in sub-grid schemes and calibration procedures that each model uses. Models that do not include active galactic nucleus feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.The authors would like the acknowledge the Centre for High Performance Computing in Rosebank, Cape Town, for financial support and for hosting the ‘Comparison Cape Town’ workshop in 2016, July. The authors would further like to acknowledge the support of the International Centre for Radio Astronomy Research (ICRAR) node at the University of Western Australia (UWA) in hosting the precursor workshop ‘Perth Simulated Cluster Comparison’ in 2015, March; the financial support of the UWA Research Collaboration Award 2014 and 2015 schemes; the financial support of the ARC Centre of Excellence for All Sky Astrophysics (CAASTRO) CE110001020 and ARC Discovery Projects DP130100117 and DP140100198. We would also like to thank the Instituto de Fisica Teorica (IFT-UAM/CSIC in Madrid) for its support, via the Centro de Excelencia Severo Ochoa Program under Grant No. SEV- 2012-0249, during the three-week workshop ‘nIFTy Cosmology’ in 2014, where the foundation for the whole comparison project was established. JA acknowledges support from a post-graduate award from STFC. PJE is supported by the SSimPL programme and the Sydney Institute for Astronomy (SIfA) and Australian Research Council (ARC) grants DP130100117 and DP140100198. AK is supported by the Ministerio de Econom´ıa y Competitividad (MINECO) in Spain through grant AYA2012-31101 as well as the ConsoliderIngenio 2010 Programme of the Spanish Ministerio de Ciencia e Innovacion (MICINN) under grant MultiDark CSD2009-00064. ´ He also acknowledges support from the ARC grant DP140100198. He further thanks Noonday Underground for surface noise. STK acknowledges support from STFC through grant ST/L000768/1. CP acknowledges the support of the ARC through Future Fellowship FT130100041 and Discovery Project DP140100198. WC and CP acknowledge the support of ARC DP130100117. GY and FS acknowledge support from MINECO (Spain) through the grant AYA 2012-31101. GY thanks also the Red Espanola de Supercomputa- ˜ cion for granting the computing time in the Marenostrum Supercomputer at BSC, where all the MUSIC simulations have been performed. AMB is supported by the DFG Research Unit 1254 ‘Magnetisation of interstellar and intergalactic media’ and by the DFG Cluster of Excellence ‘Universe’. GM acknowledge support from the PRIN-MIUR 2012 Grant ‘The Evolution of Cosmic Baryons’ funded by the Italian Minister of University and Research, by the PRIN-INAF 2012 Grant ‘Multi-scale Simulations of Cosmic Structures’, by the INFN INDARK Grant and by the ‘Consorzio per la Fisica di Trieste’. IGM acknowledges support from an STFC Advanced Fellowship. EP acknowledges support by the ERC grant ‘The Emergence of Structure During the Epoch of Reionization’

nIFTy galaxy cluster simulations II: radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction

nIFTy galaxy cluster simulations I: dark matter & non-radiative models

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods