No Arabic abstract
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 argue that the global signal of neutral hydrogen 21cm line can be a powerful probe of primordial power spectrum on small scales. Since the amplitude of small scale primordial fluctuations is important to determine the early structure formation and the timing when the sources of Lyman ${alpha}$ photons are produced, they in turn affect the neutral hydrogen 21cm line signal. We show that the information of the position of the absorption trough can severely constrain the small scale amplitude of primordial fluctuations once astrophysical parameters relevant to the 21cm line signal are fixed. We also discuss how the uncertainties of astrophysical parameters affect the constraints.
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.
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.
Several satellite missions have uncovered a series of potential anomalies in the fluctuation spectrum of the cosmic microwave background temperature, including: (1) an unexpectedly low level of correlation at large angles, manifested via the angular correlation function, C(theta); and (2) missing power in the low multipole moments of the angular power spectrum, C_ell. Their origin is still debated, however, due to a persistent lack of clarity concerning the seeding of quantum fluctuations in the early Universe. A likely explanation for the first of these appears to be a cutoff, k_min=(3.14 +/- 0.36) x 10^{-4} Mpc^{-1}, in the primordial power spectrum, P(k). Our goal in this paper is twofold: (1) we examine whether the same k_min can also self-consistently explain the missing power at large angles, and (2) we confirm that the of this cutoff in P(k) does not adversely affect the remarkable consistency between the prediction of Planck-LCDM and the Planck measurements at ell > 30. We use the publicly available code CAMB to calculate the angular power spectrum, based on a line-of-sight approach. The code is modified slightly to include the additional parameter (i.e., k_min) characterizing the primordial power spectrum. In addition to this cutoff, the code optimizes all of the usual standard-model parameters. In fitting the angular power spectrum, we find an optimized cutoff, k_min = 2.04^{+1.4}_{-0.79} x 10^{-4} Mpc^{-1}, when using the whole range of ells, and k_min=3.3^{+1.7}_{-1.3} x 10^{-4} Mpc^{-1}, when fitting only the range ell < 30, where the Sachs-Wolfe effect is dominant. These are fully consistent with the value inferred from C(theta), suggesting that both of these large-angle anomalies may be due to the same truncation in P(k).
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.