Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAccelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory
69
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The linear matter power spectrum is an essential ingredient in all theoretical models for interpreting large-scale-structure observables. Although Boltzmann codes such as CLASS or CAMB are very efficient at computing the linear spectrum, the analysis of data usually requires $10^4$-$10^6$ evaluations, which means this task can be the most computationally expensive aspect of data analysis. Here, we address this problem by building a neural network emulator that provides the linear theory matter power spectrum in about one millisecond with 0.3% accuracy over $10^{-3} le k [h,{rm Mpc}^{-1}] < 30$. We train this emulator with more than 150,000 measurements, spanning a broad cosmological parameter space that includes massive neutrinos and dynamical dark energy. We show that the parameter range and accuracy of our emulator is enough to get unbiased cosmological constraints in the analysis of a Euclid-like weak lensing survey. Complementing this emulator, we train 15 other emulators for the cross-spectra of various linear fields in Eulerian space, as predicted by 2nd-order Lagrangian Perturbation theory, which can be used to accelerate perturbative bias descriptions of galaxy clustering. Our emulators are especially designed to be used in combination with emulators for the nonlinear matter power spectrum and for baryonic effects, all of which are publicly available at http://www.dipc.org/bacco.
rate research
Read More
The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional Fast Fourier Transforms. This is exact and a few orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. The method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.
We present a simple numerical scheme for perturbation theory (PT) calculations of large-scale structure. Solving the evolution equations for perturbations numerically, we construct the PT kernels as building blocks of statistical calculations, from which the power spectrum and/or correlation function can be systematically computed. The scheme is especially applicable to the generalized structure formation including modified gravity, in which the analytic construction of PT kernels is intractable. As an illustration, we show several examples for power spectrum calculations in $f(R)$ gravity and $Lambda$CDM models.
Perturbation theory (PT) calculation of large-scale structure has been used to interpret the observed non-linear statistics of large-scale structure at the quasi-linear regime. In particular, the so-called standard perturbation theory (SPT) provides a basis for the analytical computation of the higher-order quantities of large-scale structure. Here, we present a novel, grid-based algorithm for the SPT calculation, hence named GridSPT, to generate the higher-order density and velocity fields from a given linear power spectrum. Taking advantage of the Fast Fourier Transform, the GridSPT quickly generates the nonlinear density fields at each order, from which we calculate the statistical quantities such as non-linear power spectrum and bispectrum. Comparing the density fields (to fifth order) from GridSPT with those from the full N-body simulations in the configuration space, we find that GridSPT accurately reproduces the N-body result on large scales. The agreement worsens with smaller smoothing radius, particularly for the under-dense regions where we find that 2LPT (second-order Lagrangian perturbation theory) algorithm performs well.
The form of the primordial power spectrum (PPS) of cosmological scalar (matter density) perturbations is not yet constrained satisfactorily in spite of the tremendous amount of information from the Cosmic Microwave Background (CMB) data. While a smooth power-law-like form of the PPS is consistent with the CMB data, some PPS with small non-smooth features at large scales can also fit the CMB temperature and polarization data with similar statistical evidence. Future CMB surveys cannot help distinguish all such models due to the cosmic variance at large angular scales. In this paper, we study how well we can differentiate be- tween such featured forms of the PPS not otherwise distinguishable using CMB data. We ran 15 N-body DESI-like simulations of these models to explore this approach. Showing that statistics such as the halo mass function and the two-point correlation function are not able to distinguish these models in a DESI-like survey, we advocate to avoid reducing the dimensionality of the problem by demonstrating that the use of a simple three-dimensional count-in-cell density field can be much more effective for the purpose of model distinction.
We consider the growth of primordial dark matter halos seeded by three crossed initial sine waves of various amplitudes. Using a Lagrangian treatment of cosmological gravitational dynamics, we examine the convergence properties of a high-order perturbative expansion in the vicinity of shell-crossing, by comparing the analytical results with state-of-the-art high resolution Vlasov-Poisson simulations. Based on a quantitative exploration of parameter space, we study explicitly for the first time the convergence speed of the perturbative series, and find, in agreement with intuition, that it slows down when going from quasi one-dimensional initial conditions (one sine wave dominating) to quasi triaxial symmetry (three sine waves with same amplitude). In most cases, the system structure at collapse time is, as expected, very similar to what is obtained with simple one-dimensional dynamics, except in the quasi-triaxial regime, where the phase-space sheet presents a velocity spike. In all cases, the perturbative series exhibits a generic convergence behavior as fast as an exponential of a power-law of the order of the expansion, allowing one to numerically extrapolate it to infinite order. The results of such an extrapolation agree remarkably well with the simulations, even at shell-crossing.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
