Nonlinear Regime Of Structure Formation Research Articles

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

ABSTRACT We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A bacusS ummitN-body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

Sensitivity of Cosmological Parameter Estimation to Nonlinear Prescription from Cosmic Shear

Several ongoing and upcoming large-scale structure surveys aim to explore the nonlinear regime of structure formation with high precision. Making reliable cosmological inferences from these observations necessitates precise theoretical modeling of the mildly nonlinear regime. In this work we explore how the choice of nonlinear prescription would impact parameter estimation from cosmic shear measurements for a Euclid-like survey. Specifically, we employ two different nonlinear prescriptions of halofit and the Effective Field Theory of the Large Scale Structure and compare their measurements for the three different cosmological scenarios of ΛCDM, wCDM, and (w 0, w a ) CDM. We also investigate the impact of different nonlinear cutoff schemes on parameter estimation. We find that the predicted errors on most parameters shrink considerably as smaller scales are included in the analysis, with the amount depending on the nonlinear prescription and the cutoff scheme used. We use predictions from the halofit model to analyze the mock data from DarkSky N-body simulations and quantify the parameter bias introduced in the measurements due to the choice of nonlinear prescription. We observe that σ 8 and n s have the largest measurement bias induced by inaccuracies of the halofit model.

The e-MANTIS emulator: fast predictions of the non-linear matter power spectrum in f(R)CDM cosmology

ABSTRACT In order to probe modifications of gravity at cosmological scales, one needs accurate theoretical predictions. N-body simulations are required to explore the non-linear regime of structure formation but are very time consuming. In this work, we release a new public emulator, dubbed e-mantis, that performs an accurate and fast interpolation between the predictions of f(R) modified gravity cosmological simulations, run with ecosmog. We sample a wide 3D parameter space given by the current background scalar field value $10^{-7}\lt \left|f_{R_0}\right|\lt 10^{-4}$, matter density 0.24 < Ωm < 0.39, and primordial power spectrum normalization 0.6 < σ8 < 1.0, with 110 points sampled from a Latin hypercube. For each model we perform pairs of f(R)CDM and ΛCDM simulations covering an effective volume of $\left(560 \, h^{-1}\, \mathrm{Mpc}\right)^3$ with a mass resolution of ∼2 × 1010h−1M⊙. We build an emulator for the matter power spectrum boost B(k) = Pf(R)(k)/PΛCDM(k) using a Gaussian process regression method. The boost is mostly independent of h, ns, and Ωb, which reduces the dimensionality of the relevant cosmological parameter space. Additionally, it is more robust against statistical and systematic errors than the raw power spectrum, thus strongly reducing our computational needs. According to our dedicated study of numerical systematics, the resulting emulator has an estimated maximum error of 3 per cent across the whole cosmological parameter space, for scales $0.03 \ h\, \mathrm{Mpc}^{-1} \lt k \lt 7 \ h\, \mathrm{Mpc}^{-1}$, and redshifts 0 < z < 2, while in most cases the accuracy is better than 1 per cent. Such an emulator could be used to constrain f(R) gravity with weak lensing analyses.

The Mira–Titan Universe – IV. High-precision power spectrum emulation

ABSTRACT Modern cosmological surveys are delivering data sets characterized by unprecedented quality and statistical completeness; this trend is expected to continue in the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the non-linear regime of structure formation where cosmological simulations are the primary means of obtaining the required information. The computational cost of sufficiently resolved large-volume simulations makes it prohibitive to run very large ensembles. Nevertheless, precision emulators built on a tractable number of high-quality simulations can be used to build very fast prediction schemes to enable a variety of cosmological inference studies. We have recently introduced the Mira–Titan Universe simulation suite designed to construct emulators for a range of cosmological probes. This gravity-only set of simulations covers the standard six cosmological parameters {ωm, ωb, σ8, h, ns, w0} and, in addition, includes massive neutrinos and a dynamical dark energy equation of state {ων, wa}. In this paper, we present the final emulator for the matter power spectrum based on 111 cosmological simulations, each covering a (2.1 Gpc)3 volume and evolving 32003 particles. An additional set of 1776 lower resolution simulations and TimeRG perturbation theory results for the power spectrum are used to cover scales straddling the linear to mildly non-linear regimes (maximum wavenumber k = 5 Mpc−1). The emulator provides predictions at the 2–3 per cent level of accuracy over a wide range of cosmological parameters and is publicly released as part of this paper.

Magrathea-Pathfinder: a 3D adaptive-mesh code for geodesic ray tracing in N-body simulations

We introduce MAGRATHEA-PATHFINDER, a relativistic ray-tracing framework that can reconstruct the past light cone of observers in cosmological simulations. The code directly computes the 3D trajectory of light rays through the null geodesic equations, with the weak-field limit as its only approximation. This approach offers high levels of versatility while removing the need for many of the standard ray-tracing approximations such as plane-parallel, Born, or multiple-lens. Moreover, the use of adaptive integration steps and interpolation strategies based on adaptive-mesh refinement grids allows MAGRATHEA-PATHFINDER to accurately account for the nonlinear regime of structure formation and fully take advantage of the small-scale gravitational clustering. To handle very large N-body simulations, the framework has been designed as a high-performance computing post-processing tool relying on a hybrid paralleliza-tion that combines MPI tasks with C++11 std::threads. In this paper, we describe how realistic cosmological observables can be computed from numerical simulation using ray-tracing techniques. We discuss in particular the production of simulated catalogues and sky maps that account for all the observational effects considering first-order metric perturbations (such as peculiar velocities, gravitational potential, integrated Sachs-Wolfe, time-delay, and gravitational lensing). We perform convergence tests of our gravitational lensing algorithms and conduct performance benchmarks of the null geodesic integration procedures. MAGRATHEA-PATHFINDER introduces sophisticated ray-tracing tools to make the link between the space of N-body simulations and light-cone observables. This should provide new ways of exploring existing cosmological probes and building new ones beyond standard assumptions in order to prepare for the next generation of large-scale structure surveys.

Theoretical and numerical perspectives on cosmic distance averages

The interpretation of cosmological observations relies on a notion of an average Universe, which is usually considered as the homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker (FLRW) model. However, inhomogeneities may statistically bias the observational averages with respect to FLRW, notably for distance measurements, due to a number of effects such as gravitational lensing and redshift perturbations. In this article, we review the main known theoretical results on average distance measures in cosmology, based on second-order perturbation theory, and we fill in some of their gaps. We then comprehensively test these theoretical predictions against ray tracing in a high-resolution dark-matter N-body simulation. This method allows us to describe the effect of small-scale inhomogeneities deep into the non-linear regime of structure formation on light propagation up to z = 10. We find that numerical results are in remarkably good agreement with theoretical predictions in the limit of super-sample variance. No unexpectedly large bias originates from very small scales, whose effect is fully encoded in the non-linear power spectrum. Specifically, the directional average of the inverse amplification and the source-averaged amplification are compatible with unity; the change in area of surfaces of constant cosmic time is compatible with zero; the biases on other distance measures, which can reach slightly less than 1% at high redshift, are well understood. As a side product, we also confront the predictions of the recent finite-beam formalism with numerical data and find excellent agreement.

Modelling the matter bispectrum at small scales in modified gravity

Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. {We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to kmax = 4 h/Mpc, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, f(R) gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as Euclid. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.}

The galaxy–halo connection in modified gravity cosmologies: environment dependence of galaxy luminosity function

Abstract We investigate the dependence of the galaxy–halo connection and galaxy density field in modified gravity models using the N-body simulations for f(R) and nDGP models at z = 0. Because of the screening mechanisms employed by these models, chameleon and Vainshtein, haloes are clustered differently in the non-linear regime of structure formation. We quantify their deviations in the galaxy density field from the standard Λ cold dark matter (ΛCDM) model under different environments. We populate galaxies in haloes via the (sub)halo abundance matching. Our main results are as follows: (1) The galaxy–halo connection strongly depends on the gravity model; a maximum variation of ${\sim }40{{\ \rm per\ cent}}$ is observed between halo occupational distribution (HOD) parameters; (2) f(R) gravity models predict an excess of galaxies in low-density environments of ${\sim }10{{\ \rm per\ cent}}$ but predict a deficit of ${\sim }10{{\ \rm per\ cent}}$ at high-density environments for |fR0| = 10−4 and 10−6 while |fR0| = 10−5 predicts more high-density structures; nDGP models are consistent with ΛCDM; (3) different gravity models predict different dependences of the galaxy luminosity function (GLF) with the environment, especially in void-like regions we find differences around ${\sim }10{{\ \rm per\ cent}}$ for the f(R) models while nDPG models remain closer to ΛCDM for low-luminosity galaxies but there is a deficit of ${\sim }11{{\ \rm per\ cent}}$ for high-luminosity galaxies in all environments. We conclude that the dependence of the GLF with environment might provide a test to distinguish between gravity models and their screening mechanisms from the ΛCDM. We provide HOD parameters for the gravity models analysed in this paper.

Measuring the growth of structure with intensity mapping surveys

Line intensity mapping offers a new avenue for constraining cosmological parameters in the high redshift Universe. However measurements of the growth of structure, a sensitive probe of gravity, are affected by a well known degeneracy with astrophysical parameters, encoded in the mean brightness temperature of the specific line. In this work we show how to break this degeneracy, to a level that could allow constraints of the amplitude of cosmological fluctuations at the percent level, using information in the mildly non-linear regime of structure formation as described by Lagrangian Perturbation Theory. We focus on the 21-cm line with forecasts for HIRAX and the proposed Stage II experiment as illustrations.

Large-scale structure probes of modified gravity

Observations of the evolution of large-scale structures in the Universe provides unique tools to confront Einstein’s theory of General Relativity on cosmological scales. We review weak gravitational lensing and galaxy clustering studies, discussing how these can be used in combination in order to constrain a range of different modified gravity theories. We argue that in order to maximise the future information gain from these probes, theoretical effort will be required in order to model the impact of beyond-Einstein gravity in the nonlinear regime of structure formation.

Complete super-sample lensing covariance in the response approach

We derive the complete super-sample covariance (SSC) of the matter and weak lensing convergence power spectra using the power spectrum response formalism to accurately describe the coupling of super- to sub-survey modes. The SSC term is completely characterized by the survey window function, the nonlinear matter power spectrum and the full first-order nonlinear power spectrum response function, which describes the response to super-survey density and tidal field perturbations. Generalized separate universe simulations can efficiently measure these responses in the nonlinear regime of structure formation, which is necessary for lensing applications. We derive the lensing SSC formulae for two cases: one under the Limber and flat-sky approximations, and a more general one that goes beyond the Limber approximation in the super-survey mode and is valid for curved sky applications. Quantitatively, we find that for sky fractions fsky ≈ 0.3 and a single source redshift at zS=1, the use of the flat-sky and Limber approximation underestimates the total SSC contribution by ≈ 10%. The contribution from super-survey tidal fields to the lensing SSC, which has not been included in cosmological analyses so far, is shown to represent about 5% of the total lensing covariance on multipoles ℓ1,ℓ2 ≳ 300. The SSC is the dominant off-diagonal contribution to the total lensing covariance, making it appropriate to include these tidal terms and beyond flat-sky/Limber corrections in cosmic shear analyses.

Probing screened modified gravity with nonlinear structure formation

We review the effects of modified gravity theories, in particular the symmetron and [Formula: see text] gravity, on the nonlinear regime of structure formation. In particular, we investigate the velocity dispersion of galaxy clusters as a function of the halo masses, how the matter power spectra depend on the coupling, range, and screening scale of the fifth force, and on possible ways of detecting violations of the equivalence principle using the mass inferred via lensing methods versus the mass inferred via dynamical methods. Furthermore, we show how one could use different voids statistics as one of the most promising probes of modified gravity.

Beyond Kaiser bias: mildly non-linear two-point statistics of densities in distant spheres

Simple parameter-free analytic bias functions for the two-point correlation of densities in spheres at large separation are presented. These bias functions generalize the so-called Kaiser bias to the mildly non-linear regime for arbitrary density contrasts. The derivation is carried out in the context of large deviation statistics while relying on the spherical collapse model. A logarithmic transformation provides a saddle approximation which is valid for the whole range of densities and shown to be accurate against the 30 Gpc cube state-of-the-art Horizon Run 4 simulation. Special configurations of two concentric spheres that allow to identify peaks are employed to obtain the conditional bias and a proxy to BBKS extrema correlation functions. These analytic bias functions should be used jointly with extended perturbation theory to predict two-point clustering statistics as they capture the non-linear regime of structure formation at the percent level down to scales of about 10 Mpc/h at redshift 0. Conversely, the joint statistics also provide us with optimal dark matter two-point correlation estimates which can be applied either universally to all spheres or to a restricted set of biased (over- or underdense) pairs. Based on a simple fiducial survey, this estimator is shown to perform five times better than usual two-point function estimators. Extracting more information from correlations of different types of objects should prove essential in the context of upcoming surveys like Euclid, DESI, PFS or LSST.

Extensive investigation of the generalized dark matter model

The Cold Dark Matter (CDM) model, wherein the dark matter is treated as a pressureless perfect fluid, provides a good fit to galactic and cosmological data. With the advent of precision cosmology, it should be asked whether this simplest model needs to be extended, and whether doing so could improve our understanding of the properties of dark matter. One established parameterisation for generalising the CDM fluid is the Generalised Dark Matter (GDM) model, in which dark matter is an imperfect fluid with pressure and shear viscosity that fulfill certain closure equations. We investigate these closure equations and the three new parametric functions they contain: the background equation of state w, the speed of sound c_s^2 and the viscosity c_{vis}^2. Taking these functions to be constant parameters, we analyse an exact solution of the perturbed Einstein equations in a GDM-dominated universe and discuss the main effects of the three parameters on the Cosmic Microwave Background (CMB). Our analysis suggests that the CMB alone is not able to distinguish between the GDM sound speed and viscosity parameters, but that other observables, such as the matter power spectrum, are required to break this degeneracy. We also consider other descriptions of imperfect fluids that have a non-perturbative definition and relate these to the GDM model. In particular, we consider scalar fields, an effective field theory (EFT) of fluids, an EFT of Large Scale Structure, non-equilibrium thermodynamics and tightly-coupled fluids. These descriptions could be used to extend the GDM model into the nonlinear regime of structure formation, which is necessary if the wealth of data available on those scales is to be employed in constraining the model. We also derive the initial conditions for adiabatic and isocurvature perturbations and provide the result in a form ready for implementation in Einstein-Boltzmann solvers.

Nonlinear structure formation in gravity theories beyond general relativity

We investigate the effects of modified gravity theories, in particular, the symmetron and f(R) gravity, on the nonlinear regime of structure formation. In particular, we investigate the velocity dispersion of galaxy clusters as a function of the halo masses, how the matter power spectra depend on the coupling, range and screening scale of the fifth force, and on possible ways of detecting violations of the equivalence principle using the mass inferred via lensing methods versus the mass inferred via dynamical methods.

Weak lensing of large scale structure in the presence of screening

A number of alternatives to general relativity exhibit gravitational screening in the non-linear regime of structure formation. We describe a set of algorithms that can produce weak lensing maps of large scale structure in such theories and can be used to generate mock surveys for cosmological analysis. By analysing a few basic statistics we indicate how these alternatives can be distinguished from general relativity with future weak lensing surveys.

The anisotropic line correlation function as a probe of anisotropies in galaxy surveys

We propose an anisotropic generalisation of the line correlation function (ALCF) to separate and quantify phase information in the large-scale structure of galaxies. The line correlation function probes the strictly non-linear regime of structure formation and since phase information drops out of the power spectrum, the line correlation function provides a complementary tool to commonly used techniques based on two-point statistics. Furthermore, it is independent of linear bias as well as the Gaussian variance on the modulus of the density field and thus may also prove to be advantageous compared to the bispectrum or similar higher-order statistics for certain cases. For future applications it is vital, though, to be able to account for observational effects that cause anisotropies in the distribution of galaxies. Based on a number of numerical studies, we find that our ALCF is well suited to accomplish this task and we demonstrate how the Alcock-Paczynski effect and kinematical redshift-space distortions can in principle be measured via the ALCF.

The information content of galaxy surveys

A large fraction of this thesis is dedicated to the study of the information content of random fields with heavy tails, in particular the lognormal field, a model for the matter density fluctuation field. It is well known that in the nonlinear regime of structure formation, the matter fluctuation field develops such large tails. It has also been suggested that fields with large tails are not necessarily well described by the hierarchy of $N$-point functions. In this thesis, we are able to make this last statement precise and with the help of the lognormal model to quantify precisely its implications for inference on cosmological parameters : we find as our main result that only a tiny fraction of the total Fisher information of the field is still contained in the hierarchy of $N$-point moments in the nonlinear regime, rendering parameter inference from such moments very inefficient. We show that the hierarchy fails to capture the information that is contained in the underdense regions, which at the same time are found to be the most rich in information. We find further our results to be very consistent with numerical analysis using $N$-body simulations. We also discuss these issues with the help of explicit families of fields with the same hierarchy of $N$-point moments defined in this work. A similar analysis is then applied to the convergence field, the weighted projection of the matter density fluctuation field along the line of sight, with similar conclusions. We also show how simple mappings can correct for this inadequacy, consistently with previous findings in the literature (Abridged) .

Statistics of Caustics in Large-Scale Structure Formation

AbstractThe cosmic web is a complex spatial pattern of walls, filaments, cluster nodes and underdense void regions. It emerged through gravitational amplification from the Gaussian primordial density field. Here we infer analytical expressions for the spatial statistics of caustics in the evolving large-scale mass distribution. In our analysis, following the quasi-linear Zel'dovich formalism and confined to the 1D and 2D situation, we compute number density and correlation properties of caustics in cosmic density fields that evolve from Gaussian primordial conditions. The analysis can be straightforwardly extended to the 3D situation. We moreover, are currently extending the approach to the non-linear regime of structure formation by including higher order Lagrangian approximations and Lagrangian effective field theory.

THE COYOTE UNIVERSE EXTENDED: PRECISION EMULATION OF THE MATTER POWER SPECTRUM

Modern sky surveys are returning precision measurements of cosmological statistics such as weak lensing shear correlations, the distribution of galaxies, and cluster abundance. To fully exploit these observations, theorists must provide predictions that are at least as accurate as the measurements, as well as robust estimates of systematic errors that are inherent to the modeling process. In the nonlinear regime of structure formation, this challenge can only be overcome by developing a large-scale, multi-physics simulation capability covering a range of cosmological models and astrophysical processes. As a first step to achieving this goal, we have recently developed a prediction scheme for the matter power spectrum (a so-called emulator), accurate at the 1% level out to k~1/Mpc and z=1 for wCDM cosmologies based on a set of high-accuracy N-body simulations. It is highly desirable to increase the range in both redshift and wavenumber and to extend the reach in cosmological parameter space. To make progress in this direction, while minimizing computational cost, we present a strategy that maximally re-uses the original simulations. We demonstrate improvement over the original spatial dynamic range by an order of magnitude, reaching k~10 h/Mpc, a four-fold increase in redshift coverage, to z=4, and now include the Hubble parameter as a new independent variable. To further the range in k and z, a new set of nested simulations run at modest cost is added to the original set. The extension in h is performed by including perturbation theory results within a multi-scale procedure for building the emulator. This economical methodology still gives excellent error control, ~5% near the edges of the domain of applicability of the emulator. A public domain code for the new emulator is released as part of the work presented in this paper.