We measure and quantify properties of galactic outflows and diffuse gas at z ≥ 1 in cosmological hydrodynamical simulations. Our novel subresolution model, Multi-Phase Particle Integrator (MUPPI), ...implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at z = 2 by analysing central galaxies having a stellar mass higher than 109 M⊙. The outflow velocity and mass outflow rate (
$\dot{M}_{\rm out}$
) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor (η) is between 0.2 and 10. The comparison effective model generates a constant outflow velocity, and a negative correlation of η with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than 80 per cent over z = 1–5. The outflow velocity correlation with SFR becomes flatter at z = 1, and η displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and effective models produce significant outflows at ∼1/10 of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the effective model. Our MUPPI model, using local properties of gas in the subresolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties, and consistent with observations.
New challenges in Astronomy and Astrophysics (AA) are urging the need for many exceptionally computationally intensive simulations. “Exascale” (and beyond) computational facilities are mandatory to ...address the size of theoretical problems and data coming from the new generation of observational facilities in AA. Currently, the High-Performance Computing (HPC) sector is undergoing a profound phase of innovation, in which the primary challenge to the achievement of the “Exascale” is the power consumption. The goal of this work is to give some insights about performance and energy footprint of contemporary architectures for a real astrophysical application in an HPC context. We use a state-of-the-art N-body application that we re-engineered and optimized to exploit the heterogeneous underlying hardware fully. We quantitatively evaluate the impact of computation on energy consumption when running on four different platforms. Two of them represent the current HPC systems (Intel-based and equipped with NVIDIA GPUs), one is a micro-cluster based on ARM-MPSoC, and one is a “prototype towards Exascale” equipped with ARM-MPSoCs tightly coupled with FPGAs. We investigate the behavior of the different devices where the high-end GPUs excel in terms of time-to-solution while MPSoC-FPGA systems outperform GPUs in power consumption. Our experience reveals that considering FPGAs for computationally intensive application seems very promising, as their performance is improving to meet the requirements of scientific applications. This work can be a reference for future platform development for astrophysics applications where computationally intensive calculations are required.
ABSTRACT
We employ a set of Magneticum cosmological hydrodynamic simulations that span over 15 different cosmologies, and extract masses and concentrations of all well-resolved haloes between z = ...0 and 1 for critical overdensities $\Delta _\textrm {vir}, \Delta _{200c}, \Delta _{500c}, \Delta _{2500c}$ and mean overdensity Δ200m. We provide the first mass–concentration (Mc) relation and sparsity relation (i.e. MΔ1 − MΔ2 mass conversion) of hydrodynamic simulations that is modelled by mass, redshift, and cosmological parameters Ωm, Ωb, σ8, h0 as a tool for observational studies. We also quantify the impact that the Mc relation scatter and the assumption of Navarro–Frank–White (NFW) density profiles have on the uncertainty of the sparsity relation. We find that converting masses with the aid of an Mc relation carries an additional fractional scatter ($\approx 4{{\ \rm per\ cent}}$) originated from deviations from the assumed NFW density profile. For this reason, we provide a direct mass–mass conversion relation fit that depends on redshift and cosmological parameters. We release the package hydro_mc, a python tool that perform all kind of conversions presented in this paper.
ABSTRACT Galaxy cluster masses derived from observations of weak lensing suffer from a number of biases affecting the accuracy of mass-observable relations calibrated from such observations. In ...particular, the choice of the cluster centre plays a prominent role in biasing inferred masses. In the past, empirical miscentring distributions have been used to address this issue. Using hydrodynamic simulations, we aim to test the accuracy of weak lensing mass bias predictions based on such miscentring distributions by comparing the results to mass biases computed directly using intracluster medium (ICM)-based centres from the same simulation. We construct models for fitting masses to both centred and miscentred Navarro–Frenk–White profiles of reduced shear, and model the resulting distributions of mass bias with normal and lognormal distributions. We find that the standard approach of using miscentring distributions leads to an overestimation of cluster masses at levels of between 2 per cent and 6 per cent when compared to the analysis in which actual simulated ICM centres are used, even when the underlying miscentring distributions match in terms of the miscentring amplitude. While we find that neither lognormal nor normal distributions are generally reliable for accurately modelling the shapes of the mass bias distributions, both models can serve as reasonable approximations in practice.
ABSTRACT
Dark matter (DM) self-interactions have been proposed to solve problems on small length scales within the standard cold DM cosmology. Here, we investigate the effects of DM self-interactions ...in merging systems of galaxies and galaxy clusters with equal and unequal mass ratios. We perform N-body DM-only simulations of idealized setups to study the effects of DM self-interactions that are elastic and velocity-independent. We go beyond the commonly adopted assumption of large-angle (rare) DM scatterings, paying attention to the impact of small-angle (frequent) scatterings on astrophysical observables and related quantities. Specifically, we focus on DM-galaxy offsets, galaxy–galaxy distances, halo shapes, morphology, and the phase–space distribution. Moreover, we compare two methods to identify peaks: one based on the gravitational potential and one based on isodensity contours. We find that the results are sensitive to the peak finding method, which poses a challenge for the analysis of merging systems in simulations and observations, especially for minor mergers. Large DM-galaxy offsets can occur in minor mergers, especially with frequent self-interactions. The subhalo tends to dissolve quickly for these cases. While clusters in late merger phases lead to potentially large differences between rare and frequent scatterings, we believe that these differences are non-trivial to extract from observations. We therefore study the galaxy/star populations which remain distinct even after the DM haloes have coalesced. We find that these collisionless tracers behave differently for rare and frequent scatterings, potentially giving a handle to learn about the micro-physics of DM.
ABSTRACT
Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent ...DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self-scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function, and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall, we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide – for the first time – upper limits on the cross-section for frequent self-interactions.
ABSTRACT
Self-interacting dark matter (SIDM) models have the potential to solve the small-scale problems that arise in the cold dark matter paradigm. Simulations are a powerful tool for studying SIDM ...in the context of astrophysics, but it is numerically challenging to study differential cross-sections that favour small-angle scattering (as in light-mediator models). Here, we present a novel approach to model frequent scattering based on an effective drag force, which we have implemented into the N-body code gadget-3. In a range of test problems, we demonstrate that our implementation accurately models frequent scattering. Our implementation can be used to study differences between SIDM models that predict rare and frequent scattering. We simulate core formation in isolated dark matter haloes, as well as major mergers of galaxy clusters and find that SIDM models with rare and frequent interactions make different predictions. In particular, frequent interactions are able to produce larger offsets between the distribution of galaxies and dark matter in equal-mass mergers.
ABSTRACT
Dark matter self-interactions may have the capability to solve or at least mitigate small-scale problems of the cosmological standard model, Lambda cold dark matter. There are a variety of ...self-interacting dark matter models that lead to distinguishable astrophysical predictions and hence varying success in explaining observations. Studies of dark matter (DM) density cores on various mass scales suggest a velocity-dependent scattering cross-section. In this work, we investigate how a velocity dependence alters the evolution of the DM distribution for frequent DM scatterings and compare to the velocity-independent case. We demonstrate that these cases are qualitatively different using a test problem. Moreover, we study the evolution of the density profile of idealized DM haloes and find that a velocity dependence can lead to larger core sizes and different time-scales of core formation and core collapse. In cosmological simulations, we investigate the effect of velocity-dependent self-interaction on haloes and satellites in the mass range of ≈1011–$10^{14} \, \mathrm{M_\odot }$. We study the abundance of satellites, density, and shape profiles and try to infer qualitative differences between velocity-dependent and velocity-independent scatterings as well as between frequent and rare self-interactions. We find that a strongly velocity-dependent cross-section can significantly amplify the diversity of rotation curves, independent of the angular dependence of the differential cross-section. We further find that the abundance of satellites in general depends on both the velocity dependence and the scattering angle, although the latter is less important for strongly velocity-dependent cross-sections.
The Three Hundred: Msub–Vcirc relation Srivastava, Atulit; Cui, Weiguang; Meneghetti, Massimo ...
Monthly Notices of the Royal Astronomical Society,
03/2024, Letnik:
528, Številka:
3
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
In this study, we investigate a recent finding based on strong lensing observations, which suggests that the sub-haloes observed in clusters exhibit greater compactness compared to those ...predicted by ΛCDM simulations. To address this discrepancy, we compare the cumulative sub-halo mass function and the Msub–Vcirc relation between observed clusters and 324 simulated clusters from $\rm \small {The\,Three\,\,Hundred}$ project, focusing on the hydrodynamic resimulations using $\rm \small {Gadget-X}$ and $\rm \small {Gizmo-Simba}$ baryonic models. The cumulative sub-halo mass function of $\rm \small {Gizmo-Simba}$ simulated clusters aligns with observations, while $\rm \small {Gadget-X}$ simulations exhibit discrepancies in the lower sub-halo mass range, possibly due to its strong supernova feedback. Both $\rm \small {Gadget-X}$ and $\rm \small {Gizmo-Simba}$ simulations demonstrate a redshift evolution of the sub-halo mass function and the Vcirc function, with slightly fewer sub-haloes observed at lower redshifts. Neither the $\rm \small {Gadget-X}$ nor $\rm \small {Gizmo-Simba}$ (albeit a little closer) simulated clusters’ predictions for the Msub–Vcirc relation align with the observational result. Further investigations on the correlation between sub-halo/halo properties and the discrepancy in the Msub–Vcirc relation reveal that the sub-halo’s half mass radius and galaxy stellar age, the baryon fraction, and sub-halo distance from the cluster’s centre, as well as the halo relaxation state, play important roles on reproducing this relation. Nonetheless, challenges persist in accurately reproducing the observed Msub–Vcirc relationship within our current hydrodynamic cluster simulation that adheres to the standard ΛCDM cosmology. These challenges may stem from shortcomings in our baryon modelling, numerical intricacies within the simulation, or even potential limitations of the ΛCDM framework.