Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGeneral framework for cosmological dark matter bounds using $N$-body simulations
88
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a general framework for obtaining robust bounds on the nature of dark matter using cosmological $N$-body simulations and Lyman-alpha forest data. We construct an emulator of hydrodynamical simulations, which is a flexible, accurate and computationally-efficient model for predicting the response of the Lyman-alpha forest flux power spectrum to different dark matter models, the state of the intergalactic medium (IGM) and the primordial power spectrum. The emulator combines a flexible parameterization for the small-scale suppression in the matter power spectrum arising in non-cold dark matter models, with an improved IGM model. We then demonstrate how to optimize the emulator for the case of ultra-light axion dark matter, presenting tests of convergence. We also carry out cross-validation tests of the accuracy of flux power spectrum prediction. This framework can be optimized for the analysis of many other dark matter candidates, e.g., warm or interacting dark matter. Our work demonstrates that a combination of an optimized emulator and cosmological effective theories, where many models are described by a single set of equations, is a powerful approach for robust and computationally-efficient inference from the cosmic large-scale structure.
rate research
Read More
We demonstrate that testing for self-similarity in scale-free simulations provides an excellent tool to quantify the resolution at small scales of cosmological N-body simulations. Analysing two-point correlation functions measured in simulations using ABACUS, we show how observed deviations from self-similarity reveal the range of time and distance scales in which convergence is obtained. While the well-converged scales show accuracy below 1 percent, our results show that, with a small force softening length, the spatial resolution is essentially determined by the mass resolution. At later times the lower cut-off scale on convergence evolves in comoving units as $a^{-1/2}$ ($a$ being the scale factor), consistent with a hypothesis that it is set by two-body collisionality. A corollary of our results is that N-body simulations, particularly at high red-shift, contain a significant spatial range in which clustering appears converged with respect to the time-stepping and force softening but has not actually converged to the physical continuum result. The method developed can be applied to determine the resolution of any clustering statistic and extended to infer resolution limits for non-scale-free simulations.
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.
We revisit constraints on annihilating dark matter based on the global 21cm signature observed by EDGES. For this purpose, we used the numerical data of the latest N-body simulation performed by state-of-art standard in order to estimate the boost factor at high redshifts ($z$ = 10 - 100), which enhances the annihilation of dark matter in course of structure formation. By taking into account to what fraction injected energy from dark matter annihilation contributes to ionization, excitation and heating of intergalactic medium during dark ages, we estimated how large the global 21cm absorption can be. By assuming the thermal freezeout scenario, we find that $m_{rm DM} < 15$ GeV and $m_{rm DM} < 3$ GeV have been excluded at 95$%$ C.L. for the annihilation modes into $bbar{b}$ and $e^+ e^-$, respectively.
Gravitational softening length is one of the key parameters to properly set up a cosmological $N$-body simulation. In this paper, we perform a large suit of high-resolution $N$-body simulations to revise the optimal softening scheme proposed by Power et al. (P03). Our finding is that P03 optimal scheme works well but is over conservative. Using smaller softening lengths than that of P03 can achieve higher spatial resolution and numerically convergent results on both circular velocity and density profiles. However using an over small softening length overpredicts matter density at the inner most region of dark matter haloes. We empirically explore a better optimal softening scheme based on P03 form and find that a small modification works well. This work will be useful for setting up cosmological simulations.
We present a new method for generating initial conditions for numerical cosmological simulations in which massive neutrinos are treated as an extra set of N-body (collisionless) particles. It allows us to accurately follow the density field for both Cold Dark Matter (CDM) and neutrinos at both high and low redshifts. At high redshifts, the new method is able to reduce the shot noise in the neutrino power spectrum by a factor of more than $10^7$ compared to previous methods, where the power spectrum was dominated by shot noise at all scales. We find that our new approach also helps to reduce the noise on the total matter power spectrum on large scales, whereas on small scales the results agree with previous simulations. Our new method also allows for a systematic study of clustering of the low velocity tail of the distribution function of neutrinos. This method also allows for the study of the evolution of the overall velocity distribution as a function of the environment determined by the CDM field.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
