-
Machine Learning methods to estimate observational properties of galaxy clusters in large volume cosmological N-body simulations
Authors:
Daniel de Andres,
Gustavo Yepes,
Federico Sembolini,
Gonzalo Martínez-Muñoz,
Weiguang Cui,
Francisco Robledo,
Chia-Hsun Chuang,
Elena Rasia
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 November, 2022; v1 submitted 22 April, 2022;
originally announced April 2022.
-
Exploring the hydrostatic mass bias in MUSIC clusters: application to the NIKA2 mock sample
Authors:
Giulia Gianfagna,
Marco De Petris,
Gustavo Yepes,
Federico De Luca,
Federico Sembolini,
Weiguang Cui,
Veronica Biffi,
Florian Kéruzoré,
Juan Macías-Pérez,
Frédéric Mayet,
Laurence Perotto,
Elena Rasia,
Florian Ruppin
Abstract:
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…
▽ More
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$σ$. 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.
△ Less
Submitted 1 February, 2021; v1 submitted 7 October, 2020;
originally announced October 2020.
-
Confirmation of NIKA2 investigation of the Sunyaev-Zel'dovich effect by using synthetic clusters of galaxies
Authors:
M. De Petris,
F. Ruppin,
F. Sembolini,
R. Adam,
A. S. Baldi,
G. Cialone,
B. Comis,
F. De Luca,
G. Gianfagna,
F. Kéruzoré,
J. Macìas-Pérez,
F. Mayet,
L. Perotto,
G. Yepes
Abstract:
The NIKA2 Sunyaev-Zel'dovich Large Program (SZLP) is focused on mapping the thermal SZ signal of a representative sample of selected Planck and ACT clusters spanning the redshift range 0.5<$z$<0.9. Hydrodynamical N-body simulations prove to be a powerful tool to endorse NIKA2 capabilities for estimating the impact of IntraCluster Medium (ICM) disturbances when recovering the pressure radial profil…
▽ More
The NIKA2 Sunyaev-Zel'dovich Large Program (SZLP) is focused on mapping the thermal SZ signal of a representative sample of selected Planck and ACT clusters spanning the redshift range 0.5<$z$<0.9. Hydrodynamical N-body simulations prove to be a powerful tool to endorse NIKA2 capabilities for estimating the impact of IntraCluster Medium (ICM) disturbances when recovering the pressure radial profiles. For this goal we employ a subsample of objects, carefully extracted from the catalog "Marenostrum MUltidark SImulations of galaxy Clusters" (MUSIC), spanning equivalent redshift and mass ranges as the LPSZ. The joint analysis of real observations of the tSZ with NIKA2 and Planck enables to validate the NIKA2 pipeline and to estimate the ICM pressure profiles. Moreover, the possibility to identify "a priori" the dynamical state of the selected synthetic clusters allows us to verify the impact on the recovered ICM profile shapes and their scatters. Morphological analysis of maps of the Compton parameter seems to be a way to observationally segregate the sample based on the dynamical state in relaxed and disturbed synthetic clusters.
△ Less
Submitted 8 November, 2019;
originally announced November 2019.
-
Impact of ICM disturbances on the mean pressure profile of galaxy clusters: a prospective study of the NIKA2 SZ large program with MUSIC synthetic clusters
Authors:
F. Ruppin,
F. Sembolini,
M. De Petris,
R. Adam,
G. Cialone,
J. F. Macías-Pérez,
F. Mayet,
L. Perotto,
G. Yepes
Abstract:
$Context…
▽ More
$Context$. The mean pressure profile of the galaxy cluster population plays an essential role in cosmological analyses. An accurate characterization of the shape, intrinsic scatter, and redshift evolution of this profile is necessary to estimate some of the biases and systematic effects that currently prevent cosmological analyses based on thermal Sunyaev-Zel'dovich (tSZ) surveys from obtaining precise and unbiased cosmological constraints. This is one of the main goals of the ongoing NIKA2 SZ large program, which aims at mapping the tSZ signal of a representative cluster sample selected from the $Planck$ and ACT catalogs and spans a redshift range $0.5 < z < 0.9$. $Aims$. To estimate the impact of intracluster medium (ICM) disturbances that can be detected by NIKA2 on the mean pressure profile of galaxy clusters, we realized a study based on a synthetic cluster sample that is similar to that of the NIKA2 SZ large program. $Methods$. To reach this goal we employed the hydrodynamical simulation $Marenostrum~MUltidark~SImulations~of~galaxy~Clusters$ (MUSIC). We simulated realistic NIKA2 and $Planck$ tSZ observations, which were jointly analyzed to estimate the ICM pressure profile of each cluster. A comparison of the deprojected profiles with the true radial profiles directly extracted from the MUSIC simulation allowed us to validate the NIKA2 tSZ pipeline and to study the impact of ICM disturbances on the characterization of the ICM pressure distribution even at high redshift. After normalizing each profile by the integrated quantities estimated under the hydrostatic equilibrium hypothesis, we evaluated the mean pressure profile of the twin sample and show that it is compatible with that extracted directly from the MUSIC simulation in the scale range that can be recovered by NIKA2. We studied the impact of cluster dynamical state on both its shape and associated scatter. (abridged)
△ Less
Submitted 17 October, 2019; v1 submitted 14 January, 2019;
originally announced January 2019.
-
Kinetic Sunyaev--Zel'dovich effect in rotating galaxy clusters from MUSIC simulations
Authors:
Anna Silvia Baldi,
Marco De Petris,
Federico Sembolini,
Gustavo Yepes,
Weiguang Cui,
Luca Lamagna
Abstract:
The masses of galaxy clusters are a key tool to constrain cosmology through the physics of large-scale structure formation and accretion. Mass estimates based on X-ray and Sunyaev--Zel'dovich measurements have been found to be affected by the contribution of non-thermal pressure components, due e.g. to kinetic gas energy. The characterization of possible ordered motions (e.g. rotation) of the intr…
▽ More
The masses of galaxy clusters are a key tool to constrain cosmology through the physics of large-scale structure formation and accretion. Mass estimates based on X-ray and Sunyaev--Zel'dovich measurements have been found to be affected by the contribution of non-thermal pressure components, due e.g. to kinetic gas energy. The characterization of possible ordered motions (e.g. rotation) of the intra-cluster medium could be important to recover cluster masses accurately. We update the study of gas rotation in clusters through the maps of the kinetic Sunyaev--Zel'dovich effect, using a large sample of massive synthetic galaxy clusters ($ M_{vir} > 5\times 10^{14} h^{-1}$M$_\odot$ at $z~=~0 $) from MUSIC high-resolution simulations. We select few relaxed objects showing peculiar rotational features, as outlined in a companion work. To verify whether it is possible to reconstruct the expected radial profile of the rotational velocity, we fit the maps to a theoretical model accounting for a specific rotational law, referred as the vp2b model. We find that our procedure allows to recover the parameters describing the gas rotational velocity profile within two standard deviations, both with and without accounting for the bulk velocity of the cluster. The amplitude of the temperature distortion produced by the rotation is consistent with theoretical estimates found in the literature, and it is of the order of 23 per cent of the maximum signal produced by the cluster bulk motion. We also recover the bulk velocity projected on the line of sight consistently with the simulation true value.
△ Less
Submitted 18 May, 2018;
originally announced May 2018.
-
Morphological estimators on Sunyaev--Zel'dovich maps of MUSIC clusters of galaxies
Authors:
Giammarco Cialone,
Marco De Petris,
Federico Sembolini,
Gustavo Yepes,
Anna Silvia Baldi,
Elena Rasia
Abstract:
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. S…
▽ More
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.
△ Less
Submitted 27 March, 2018; v1 submitted 10 August, 2017;
originally announced August 2017.
-
nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region
Authors:
Jake Arthur,
Frazer R. Pearce,
Meghan E. Gray,
Pascal J. Elahi,
Alexander Knebe,
Alexander M. Beck,
Weiguang Cui,
Daniel Cunnama,
Romeel Davé,
Sean February,
Shuiyao Huang,
Neal Katz,
Scott T. Kay,
Ian G. McCarthy,
Giuseppe Murante,
Valentin Perret,
Chris Power,
Ewald Puchwein,
Alexandro Saro,
Federico Sembolini,
Romain Teyssier,
Gustavo Yepes
Abstract:
We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a…
▽ More
We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid 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\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing 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 subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.
△ Less
Submitted 23 September, 2016;
originally announced September 2016.
-
On the coherent rotation of diffuse matter in numerical simulations of galaxy clusters
Authors:
Anna Silvia Baldi,
Marco De Petris,
Federico Sembolini,
Gustavo Yepes,
Luca Lamagna,
Elena Rasia
Abstract:
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…
▽ More
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.
△ Less
Submitted 7 November, 2016; v1 submitted 7 June, 2016;
originally announced June 2016.
-
nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties
Authors:
Weiguang Cui,
Chris Power,
Alexander Knebe,
Scott T. Kay,
Federico Sembolini,
Pascal J. Elahi,
Gustavo Yepes,
Frazer Pearce,
Daniel Cunnama,
Alexander M. Beck,
Claudio Dalla Vecchia,
Romeel Davé,
Sean February,
Shuiyao Huang,
Alex Hobbs,
Neal Katz,
Ian G. McCarthy,
Giuseppe Murante,
Valentin Perret,
Ewald Puchwein,
Justin I. Read,
Alexandro Saro,
Romain Teyssier,
Robert J. Thacker
Abstract:
Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a…
▽ More
Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.
△ Less
Submitted 22 February, 2016;
originally announced February 2016.
-
nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes
Authors:
Pascal J. Elahi,
Alexander Knebe,
Frazer R. Pearce,
Chris Power,
Gustavo Yepes,
Weiguang Cui,
Daniel Cunnama,
Scott T. Kay,
Federico Sembolini,
Alexander M. Beck,
Romeel Davé,
Sean February,
Shuiyao Huang,
Neal Katz,
Ian G. McCarthy,
Giuseppe Murante,
Valentin Perret,
Ewald Puchwein,
Alexandro Saro,
Romain Teyssier
Abstract:
We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) 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/gala…
▽ More
We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) 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 full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.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 $\approx20\%$. 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}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN 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 in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.
△ Less
Submitted 10 February, 2016; v1 submitted 25 November, 2015;
originally announced November 2015.
-
nIFTy galaxy cluster simulations II: radiative models
Authors:
Federico Sembolini,
Pascal Jahan Elahi,
Frazer R. Pearce,
Chris Power,
Alexander Knebe,
Scott T. Kay,
Weiguang Cui,
Gustavo Yepes,
Alexander M. Beck,
Stefano Borgani,
Daniel Cunnama,
Romeel Davé,
Sean February,
Shuiyao Huang,
Neal Katz,
Ian G. McCarthy,
Giuseppe Murante,
Richard D. A. Newton,
Valentin Perret,
Alexandro Saro,
Joop Schaye,
Romain Teyssier
Abstract:
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 $Λ$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 an…
▽ More
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 $Λ$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.
△ Less
Submitted 11 November, 2015;
originally announced November 2015.
-
nIFTy galaxy cluster simulations I: dark matter & non-radiative models
Authors:
Federico Sembolini,
Gustavo Yepes,
Frazer R. Pearce,
Alexander Knebe,
Scott T. Kay,
Chris Power,
Weiguang Cui,
Alexander M. Beck,
Stefano Borgani,
Claudio Dalla Vecchia,
Romeel Davé,
Pascal Jahan Elahi,
Sean February,
Shuiyao Huang,
Alex Hobbs,
Neal Katz,
Erwin Lau,
Ian G. McCarthy,
Giuseppe Murante,
Daisuke Nagai,
Kaylea Nelson,
Richard D. A. Newton,
Ewald Puchwein,
Justin I. Read,
Alexandro Saro
, et al. (2 additional authors not shown)
Abstract:
We have simulated the formation of a galaxy cluster in a $Λ$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 advanc…
▽ More
We have simulated the formation of a galaxy cluster in a $Λ$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.
△ Less
Submitted 20 March, 2015;
originally announced March 2015.
-
The MUSIC of CLASH: predictions on the concentration-mass relation
Authors:
M. Meneghetti,
E. Rasia,
J. Vega,
J. Merten,
M. Postman,
G. Yepes,
F. Sembolini,
M. Donahue,
S. Ettori,
K. Umetsu,
I. Balestra,
M. Bartelmann,
N. Benitez,
A. Biviano,
R. Bouwens,
L. Bradley,
T. Broadhurst,
D. Coe,
N. Czakon,
M. De Petris,
H. Ford,
C. Giocoli,
S. Gottloeber,
C. Grillo,
L. Infante
, et al. (20 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 April, 2014; v1 submitted 4 April, 2014;
originally announced April 2014.
-
The MUSIC of Galaxy Clusters II: X-ray global properties and scaling relations
Authors:
Veronica Biffi,
Federico Sembolini,
Marco De Petris,
Riccardo Valdarnini,
Gustavo Yepes,
Stefan Gottlöber
Abstract:
We present the X-ray properties and scaling relations of a large sample of clusters extracted from the Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC) dataset. We focus on a sub-sample of 179 clusters at redshift z~0.11, with 3.2e14M_sun/h<M_vir<2e15Msun/h, complete in mass. We employed the X-ray photon simulator PHOX to obtain synthetic Chandra Observations and derive observable-like…
▽ More
We present the X-ray properties and scaling relations of a large sample of clusters extracted from the Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC) dataset. We focus on a sub-sample of 179 clusters at redshift z~0.11, with 3.2e14M_sun/h<M_vir<2e15Msun/h, 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 under-estimate 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_g500) mass; integrated Compton parameter (Y_SZ) of the Sunyaev-Zel'dovich (SZ) thermal effect; Y_X=M_g500 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%. The study of scaling relations among X-ray, intrinsic and SZ properties indicates that 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.
△ Less
Submitted 13 January, 2014;
originally announced January 2014.
-
The MUSIC of Galaxy Clusters III: Properties, evolution and Y-M scaling relation of protoclusters of galaxies
Authors:
Federico Sembolini,
Marco De Petris,
Gustavo Yepes,
Emma Foschi,
Luca Lamagna,
Stefan Gottlöber
Abstract:
In this work we study the properties of protoclusters of galaxies by employing the MUSIC set of hydrodynamical simulations, featuring a mass-limited sample of 282 resimulated clusters with available merger trees up to high redshift, and we trace the cluster formation back to $z$ = 1.5, 2.3 and 4. We study the features and redshift evolution of the mass and the spatial distribution for all the clus…
▽ More
In this work we study the properties of protoclusters of galaxies by employing the MUSIC set of hydrodynamical simulations, featuring a mass-limited sample of 282 resimulated clusters with available merger trees up to high redshift, and we trace the cluster formation back to $z$ = 1.5, 2.3 and 4. We study the features and redshift evolution of the mass and the spatial distribution for all the cluster progenitors and for the protoclusters, which we define as the most massive progenitors of the clusters identified at $z$ = 0. A natural extension to redshifts larger than 1 is applied to the estimate of the baryon content also in terms of gas and stars budgets: no remarkable variations with redshift are discovered. Furthermore, motivated by the proven potential of Sunyaev-Zel'dovich surveys to blindly search for faint distant objects, we focus on the scaling relation between total object mass and integrated Compton $y$-parameter, and we check for the possibility to extend the mass-observable paradigm to the protocluster regime, far beyond the redshift of 1, to account for the properties of the simulated objects. We find that the slope of this scaling law is steeper than what expected for a self-similarity assumption among these objects, and it increases with redshift mainly for the synthetic clusters where radiative processes, such as radiative cooling, heating processes of the gas due to UV background, star formation and supernovae feedback, are included. We use three different criteria to account for the dynamical state of the protoclusters, and find no significant dependence of the scaling parameters from the level of relaxation. Based on this, we exclude that the dynamical state is the cause of the observed deviations from self-similarity.
△ Less
Submitted 20 September, 2013;
originally announced September 2013.
-
The MUSIC of Galaxy Clusters I: Baryon properties and Scaling Relations of the thermal Sunyaev-Zel'dovich Effect
Authors:
Federico Sembolini,
Gustavo Yepes,
Marco De Petris,
Stefan Gottloeber,
Luca Lamagna,
Barbara Comis
Abstract:
We introduce the Marenostrum-MultiDark SImulations of galaxy Clusters (MUSIC) Dataset, one of the largest sample of hydrodynamically simulated galaxy clusters with more than 500 clusters and 2000 groups. The objects have been selected from two large N-body simulations and have been resimulated at high resolution using SPH together with relevant physical processes (cooling, UV photoionization, star…
▽ More
We introduce the Marenostrum-MultiDark SImulations of galaxy Clusters (MUSIC) Dataset, one of the largest sample of hydrodynamically simulated galaxy clusters with more than 500 clusters and 2000 groups. The objects have been selected from two large N-body simulations and have been resimulated at high resolution using SPH together with relevant physical processes (cooling, UV photoionization, star formation and different feedback processes). We focus on the analysis of the baryon content (gas and star) of clusters in the MUSIC dataset both as a function of aperture radius and redshift. The results from our simulations are compared with the most recent observational estimates of the gas fraction in galaxy clusters at different overdensity radii. When the effects of cooling and stellar feedbacks are included, the MUSIC clusters show a good agreement with the most recent observed gas fractions quoted in the literature. A clear dependence of the gas fractions with the total cluster mass is also evident. The impact of the aperture radius choice, when comparing integrated quantities at different redshifts, is tested: the standard definition of radius at a fixed overdensity with respect to critical density is compared with a definition based on the redshift dependent overdensity with respect to background density. We also present a detailed analysis of the scaling relations of the thermal SZ (Sunyaev Zel'dovich) Effect derived from MUSIC clusters. The integrated SZ brightness, Y, is related to the cluster total mass, M, as well as, the M-Y counterpart, more suitable for observational applications. Both laws are consistent with predictions from the self-similar model, showing a very low scatter. The effects of the gas fraction on the Y-M scaling and the presence of a possible redshift dependence on the Y-M scaling relation are also explored.
△ Less
Submitted 13 November, 2012; v1 submitted 18 July, 2012;
originally announced July 2012.
-
Baryon and Sunyaev-Zel'dovich effect properties of MareNostrum and MultiDark simulated clusters (MUSIC)
Authors:
Federico Sembolini,
Gustavo Yepes,
Marco De Petris
Abstract:
We report the first results of the MUSIC project. It consists of two data sets of resimulates clusters extracted from two large dark matter only simulations: Marenostrum Universe and Multidark. In total, the MUSIC contains more than 400 clusters resimulated with high resolution both with radiative and non-radiative physics included. Here we present the first results on the properties of the baryon…
▽ More
We report the first results of the MUSIC project. It consists of two data sets of resimulates clusters extracted from two large dark matter only simulations: Marenostrum Universe and Multidark. In total, the MUSIC contains more than 400 clusters resimulated with high resolution both with radiative and non-radiative physics included. Here we present the first results on the properties of the baryon content and the Sunyaev Zeldovich scaling relations.
△ Less
Submitted 16 November, 2011;
originally announced November 2011.