Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCosmological Perturbation Theory as a Tool for Estimating Box-Scale Effects in N-body Simulations
185
0
0.0
(
0
)
Authors
Chris Orban
Ask ChatGPT about the research
No Arabic abstract
In performing cosmological N-body simulations, it is widely appreciated that the growth of structure on the largest scales within a simulation box will be inhibited by the finite size of the simulation volume. Following ideas set forth in Seto (1999), this paper shows that standard (a.k.a. 1-loop) cosmological perturbation theory (SPT) can be used to predict, in an approximate way, the deleterious effect of the box scale on the power spectrum of density fluctuations in simulation volumes. Alternatively, this approach can be used to quickly estimate post facto the effect of the box scale on power spectrum results from existing simulations. In this way SPT can help determine whether larger box sizes or other more-sophisticated methods are needed to achieve a particular level of precision for a given application (e.g. simulations to measure the non-linear evolution of baryon acoustic oscillations). I focus on SPT in this note and show that its predictions differ only by about a factor of two or less from the measured suppression inferred from both powerlaw and $Lambda$CDM $N$-body simulations. It should be possible to improve the accuracy of these predictions through using more-sophisticated perturbation theory models. An appendix compares power spectrum measurements from the powerlaw simulations at outputs where box-scale effects are minimal to perturbation theory models and previously-published fitting functions. These power spectrum measurements are included with this paper to aid efforts to develop new perturbation theory models.
rate research
Read More
We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the synchronous-comoving gauge, allows us to go beyond the Newtonian predictions and to calculate relativistic corrections to it. We use an Einstein--de Sitter model, including both growing and decaying modes in our solutions. The impact of our results should be assessed through the implementation of the featured displacement in cosmological $N$-body simulations.
Cosmology is entering an era of percent level precision due to current large observational surveys. This precision in observation is now demanding more accuracy from numerical methods and cosmological simulations. In this paper, we study the accuracy of $N$-body numerical simulations and their dependence on changes in the initial conditions and in the simulation algorithms. For this purpose, we use a series of cosmological $N$-body simulations with varying initial conditions. We test the influence of the initial conditions, namely the pre-initial configuration (preIC), the order of the Lagrangian perturbation theory (LPT), and the initial redshift, on the statistics associated with the large scale structures of the universe such as the halo mass function, the density power spectrum, and the maximal extent of the large scale structures. We find that glass or grid pre-initial conditions give similar results at $zlesssim 2$. However, the initial excess of power in the glass initial conditions yields a subtle difference in the power spectra and the mass function at high redshifts. The LPT order used to generate the ICs of the simulations is found to play a crucial role. First-order LPT (1LPT) simulations underestimate the number of massive haloes with respect to second-order (2LPT) ones, typically by 2% at $10^{14} h^{-1} M_odot$ for an initial redshift of 23, and the small-scale power with an underestimation of 6% near the Nyquist frequency for $z_mathrm{ini} = 23$. Moreover, at higher redshifts, the high-mass end of the mass function is significantly underestimated in 1LPT simulations. On the other hand, when the LPT order is fixed, the starting redshift has a systematic impact on the low-mass end of the halo mass function.
We investigate the energy release due to the large-scale structure formation and the subsequent transfer of energy from larger to smaller scales. We calculate the power spectra for the large-scale velocity field and show that the coupling of modes results in a transfer of power predominately from larger to smaller scales. We use the concept of cumulative energy for calculating which energy amount is deposited into the small scales during the cosmological structure evolution. To estimate the contribution due to the gravitational interaction only we perform our investigations by means of dark matter simulations. The global mean of the energy transfer increases with redshift $sim (z+1)^{3}$; this can be traced back to the similar evolution of the merging rates of dark matter halos. The global mean energy transfer can be decomposed into its local contributions, which allows to determine the energy injection per mass into a local volume. The obtained energy injection rates are at least comparable with other energy sources driving the interstellar turbulence as, e.g. by the supernova kinetic feedback. On that basis we make the crude assumption that processes causing this energy transfer from large to small scales, e.g. the merging of halos, may contribute substantially to drive the ISM turbulence which may eventually result in star formation on much smaller scales. We propose that the ratio of the local energy injection rate to the energy already stored within small-scale motions is a rough measure for the probability of the local star formation efficiency applicable within cosmological large-scale n-body 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.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
