Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new usernIFTy galaxy cluster simulations I: dark matter & non-radiative models
206
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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
rate research
Read 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 $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.
In this work we calculate the Sunyaev-Zeldovich (SZ) effect due to the $e^+e^-$ from dark matter (DM) annihilation in galaxy clusters. Two candidates of DM particle, (1) the weakly-interacting massive particle (WIMP) and (2) the light dark matter (LDM) are investigated. For each case, we also consider several DM profiles with and without central cusp. We generally find smaller signals than previously reported. Moreover, the diffusion of electrons and positrons in the galaxy clusters, which was generally thought to be negligible, is considered and found to have significant effect on the central electron/positron distribution for DM profile with large spatial gradient. We find that the SZ effect from WIMP is almost always non-observable, even for the highly cuspy DM profile, and using the next generation SZ interferometer such as ALMA. Although the signal of the LDM is much larger than that of the WIMP, the final SZ effect is still very small due to the smoothing effect of diffusion. Only for the configuration with large central cusp and extremely small diffusion effect, the LDM induced SZ effect might have a bit chance of being detected.
Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxys cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the clusters outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a clusters outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history.
We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1times10^{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.01times 10^8h^{-1}text{M}_{odot}$ and $m_{text{gas}}=1.9times 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.
Abridged: We simulate a massive galaxy cluster in a LCDM Universe using three different approaches to solving the equations of non-radiative hydrodynamics: `classic Smoothed Particle Hydrodynamics (SPH); a novel SPH with a higher order dissipation switch (SPHS); and adaptive mesh refinement (AMR). We find that SPHS and AMR are in excellent agreement, with both forming a well-defined entropy core that rapidly converges with increasing mass and force resolution. By contrast, SPH exhibits rather different behaviour. At low redshift, entropy decreases systematically with decreasing cluster-centric radius, converging on ever lower central values with increasing resolution. At higher redshift, SPH is in better agreement with SPHS and AMR but shows much poorer numerical convergence. We trace these discrepancies to artificial surface tension in SPH at phase boundaries. At early times, the passage of massive substructures close to the cluster centre stirs and shocks gas to build an entropy core. At later times, artificial surface tension causes low entropy gas to sink artificially to the centre of the cluster. We use SPHS to study the contribution of numerical versus physical dissipation on the entropy core, and argue that numerical dissipation is required to ensure single-valued fluid quantities in converging flows. However, provided this dissipation occurs only at the resolution limit, and provided that it does not propagate errors to larger scales, its effect is benign. There is no requirement to build `sub-grid models of unresolved turbulence for galaxy cluster simulations. We conclude that entropy cores in non-radiative simulations of galaxy clusters are physical, resulting from entropy generation in shocked gas during cluster assembly, putting to rest the long-standing puzzle of cluster entropy cores in AMR simulations versus their apparent absence in classic SPH simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
