No Arabic abstract
Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3{sigma} constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6{sigma} constraint within the same conditions.
We derive a non-perturbative equation for the large scale structure power spectrum of long-wavelength modes. Thereby, we use an operator product expansion together with relations between the three-point function and power spectrum in the soft limit. The resulting equation encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field theory (EFT) of large scale structures receive most of their contribution from modes close to the non-linear scale. Consequently, the UV dependence found in explicit computations within standard perturbation theory stems mostly from counter-term(s). We confront a simplified version of our non-perturbative equation against existent numerical simulations, and find good agreement within the expected uncertainties. Our approach can in principle be used to precisely infer the relevance of the leading order EFT coefficient(s) using small volume simulations in an `anisotropic separate universe framework. Our results suggest that the importance of these coefficient(s) is a $sim 10 %$ effect, and plausibly smaller.
We show how the non-linearity of general relativity generates a characteristic non-Gaussian signal in cosmological large-scale structure that we calculate at all perturbative orders in a large scale limit. Newtonian gravity and general relativity provide complementary theoretical frameworks for modelling large-scale structure in $Lambda$CDM cosmology; a relativistic approach is essential to determine initial conditions which can then be used in Newtonian simulations studying the non-linear evolution of the matter density. Most inflationary models in the very early universe predict an almost Gaussian distribution for the primordial metric perturbation, $zeta$. However, we argue that it is the Ricci curvature of comoving-orthogonal spatial hypersurfaces, $R$, that drives structure formation at large scales. We show how the non-linear relation between the spatial curvature, $R$, and the metric perturbation, $zeta$, translates into a specific non-Gaussian contribution to the initial comoving matter density that we calculate for the simple case of an initially Gaussian $zeta$. Our analysis shows the non-linear signature of Einsteins gravity in large-scale structure.
In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.
Modeling the large-scale structure of the universe on nonlinear scales has the potential to substantially increase the science return of upcoming surveys by increasing the number of modes available for model comparisons. One way to achieve this is to model nonlinear scales perturbatively. Unfortunately, this involves high-dimensional loop integrals that are cumbersome to evaluate. Trying to simplify this, we show how two-loop (next-to-next-to-leading order) corrections to the density power spectrum can be reduced to low-dimensional, radial integrals. Many of those can be evaluated with a one-dimensional Fast Fourier Transform, which is significantly faster than the five-dimensional Monte-Carlo integrals that are needed otherwise. The general idea of this FFT-PT method is to switch between Fourier and position space to avoid convolutions and integrate over orientations, leaving only radial integrals. This reformulation is independent of the underlying shape of the initial linear density power spectrum and should easily accommodate features such as those from baryonic acoustic oscillations. We also discuss how to account for halo bias and redshift space distortions.
(ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB temperature anisotropy bispectrum significantly tighten the observational limits, but they are still far from the PNG level predicted by the simplest models of inflation. Probing levels below CMB sensitivities will require other methods, such as searching for the statistical imprint of PNG on galactic halo clustering. During the epoch of reionization (EoR), the first stars and galaxies released radiation into the intergalactic medium (IGM) that created ionized patches whose large-scale geometry and evolution reflected the underlying abundance and large-scale clustering of the star-forming galaxies. This statistical connection between ionized patches in the IGM and galactic halos suggests that observing reionization may be another way to constrain PNG. We employ the linear perturbation theory of reionization and semi-analytic models based on the excursion-set formalism to model the effects of PNG on the EoR. We quantify the effects of PNG on the large-scale structure of reionization by deriving the ionized density bias, i.e. ratio of ionized atomic to total matter overdensities in Fourier space, at small wavenumber. Just as previous studies found that PNG creates a scale-dependent signature in the halo bias, so, too, we find a scale-dependent signature in the ionized density bias. Our results, which differ significantly from previous attempts in the literature to characterize this PNG signature, will be applied elsewhere to predict its observable consequences, e.g. in the cosmic 21cm background.