Search

The last few years have seen the development of a promising theoretical framework for statistics of the cosmic large-scale structure -- the theory of large deviations (LDT) for modelling weak-lensing one-point statistics in the mildly non-linear regime. The goal of this series of papers is to make the leap and lay out the steps to perform an actual data analysis with this theoretical tool. Building upon the LDT framework, in this work (Paper I) we demonstrate how to accurately model the Probability Distribution Function (PDF) of a reconstructed Kaiser-Squires convergence field under a realistic mask, that of the third data release of the Dark Energy Survey (DES). We also present how weak lensing systematics and higher-order lensing corrections due to intrinsic alignments, shear biases, photo-$z$ errors and baryonic feedback can be incorporated in the modelling of the reconstructed convergence PDF. In an upcoming work (Paper II) we will then demonstrate the robustness of our modelling through simulated likelihood analyses, the final step required before applying our method to actual data., Comment: 37 pages, 9 figures and 2 tables. Submitted to JCAP, comments welcome

Detecting primordial non-Gaussianity on mildly non-linear scales requires precise modelling of late-time structure formation. Accurately predicting the impact of non-linear gravitational collapse, non-linear tracer bias and baryonic physics on the variance and higher order moments of the cosmic density field is challenging, as they strongly depend on the tails of the probability distribution function (PDF) of density fluctuations. A way around this problem is to directly analyse the bulk of the PDF instead. For this purpose we devise a new method to predict the impact of general non-Gaussian initial conditions on the late-time density PDF. With this formalism we show that - even when marginalizing over potential ignorance of the amplitude and slope of the non-linear power spectrum - an analysis of the PDF at mildly non-linear densities can measure the amplitude of different primordial bispectrum shapes to an accuracy of $\Delta f_{\mathrm{NL}}^{\mathrm{loc}} = \pm 3.1\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{equi}} = \pm 10.0\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{ortho}} = \pm 17.0\ $. This assumes a joint analysis of the PDF on smoothing scales of $15$Mpc/$h$ and $30$Mpc/$h$ in a survey volume of $V=100(\mathrm{Gpc}/h)^3$ at $z=1$, analysing only densities of $\delta(15\mathrm{Mpc}/h) \in [-0.4, 0.5]$ ($\approx 87\%$ of probability) and $\delta(30\mathrm{Mpc}/h) \in [-0.3, 0.4]$ ($\approx 95\%$ of probability). Note that a formalism closely related to ours was already successfully applied to observational data \citep{Gruen2018, Friedrich2018}, demonstrating that the methodology developed here can indeed be carried over to real data analysis.

The cosmic large-scale structure of our Universe is comprised of baryons and cold dark matter (CDM). Yet it is customary to treat these two components as a combined single-matter fluid with vanishing pressure, which is justified only for sufficiently large scales and late times. Here we go beyond the single-fluid approximation and develop the perturbation theory for two gravitationally coupled fluids while still assuming vanishing pressure. We mostly focus on perturbative expansions in powers of $D$ (or $D_+$), the linear structure growth of matter in a $��$CDM Universe with cosmological constant $��$. We derive in particular (1) explicit recursion relations for the two fluid densities, (2) complementary all-order results in the Lagrangian-coordinates approach, as well as (3) the associated component wavefunctions in a semi-classical approach to cosmic large-scale structure. In our companion paper (Hahn et al. 2020) we apply these new theoretical results to generate novel higher-order initial conditions for cosmological hydrodynamical simulations., 15 pages + appendices, one figure, MNRAS in press

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

We calculate the dark matter halo correlation function in redshift space using the Gaussian streaming model (GSM). To determine the scale dependent functions entering the streaming model we use local Lagrangian bias together with Convolution Lagrangian perturbation theory (CLPT) which constitutes an approximation to the Post-Zel'dovich approximation. On the basis of N-body simulations we demonstrate that a smoothing of the initial conditions with the Lagrangian radius improves the Zel'dovich approximation and its ability to predict the displacement field of proto-halos. Based on this observation we implement a "truncated" CLPT by smoothing the initial power spectrum and investigate the dependence of the streaming model ingredients on the smoothing scale. We find that the real space correlation functions of halos and their mean pairwise velocity are optimised if the coarse graining scale is chosen to be 1 Mpc/h at z=0, while the pairwise velocity dispersion is optimised if the smoothing scale is chosen to be the Lagrangian size of the halo. We compare theoretical results for the halo correlation function in redshift space to measurements within the Horizon Run 2 N-body simulation halo catalog. We find that this simple two-filter smoothing procedure in the spirit of the truncated Zel'dovich approximation significantly improves the GSM+CLPT prediction of the redshift space halo correlation function over the whole mass range from large galaxy to galaxy cluster-sized halos., Comment: 24 pages, 18 figures

We investigate large-scale structure formation of collisionless dark matter in the phase space description based on the Vlasov equation whose nonlinearity is induced solely by gravitational interaction according to the Poisson equation. Determining the time-evolution of density and peculiar velocity demands solving the full Vlasov hierarchy for the moments of the phase space distribution function. In the presence of long-range interaction no consistent truncation of the hierarchy is known apart from the pressure- less fluid (dust) model which is incapable of describing virialization due to the occurrence of shell-crossing singularities and the inability to generate vorticity and higher cumulants like velocity dispersion. Our goal is to find a simple ansatz for the phase space distribution function that approximates the full Vlasov distribution function without pathologies in a controlled way and therefore can serve as theoretical N-body double and as a replacement for the dust model. We argue that the coarse-grained Wigner probability distribution obtained from a wave function fulfilling the Schr��dinger-Poisson equation (SPE) is the sought-after function. We show that its evolution equation approximates the Vlasov equation and therefore also the dust equations before shell-crossing, but cures the shell-crossing singularities and is able to describe regions of multi-streaming and virialization. This feature was already employed in simulations by Widrow & Kaiser (1993). The coarse-grained Wigner ansatz allows to calculate all higher moments from density and velocity analytically, thereby incorporating nonzero higher cumulants in a self-consistent manner. On this basis we are able to show that the Schr��dinger method (ScM) automatically closes the corresponding hierarchy such that it suffices to solve the SPE in order to directly determine density and velocity and all higher cumulants., 19 p. 9 Figs, corresponds to published version

We investigate large scale structure formation of collisionless dark matter in the phase space description based on the Vlasov-Poisson equation. We present the Schr��dinger method, originally proposed by Widrow and Kaiser (1993) as numerical technique based on the Schr��dinger Poisson equation, as an analytical tool which is superior to the common standard pressureless fluid model. Whereas the dust model fails and develops singularities at shell crossing the Schr��dinger method encompasses multi-streaming and even virialization., Proceedings of IAU Symposium 308 "The Zeldovich Universe: Genesis and Growth of the Cosmic Web", 23-28 June 2014, Tallinn, Estonia

Pinning down the total neutrino mass and the dark energy equation of state is a key aim for upcoming galaxy surveys. Weak lensing is a unique probe of the total matter distribution whose non-Gaussian statistics can be quantified by the one-point probability distribution function (PDF) of the lensing convergence. We calculate the convergence PDF on mildly non-linear scales from first principles using large-deviation statistics, accounting for dark energy and the total neutrino mass. For the first time, we comprehensively validate the cosmology-dependence of the convergence PDF model against large suites of simulated lensing maps, demonstrating its percent-level precision and accuracy. We show that fast simulation codes can provide highly accurate covariance matrices, which can be combined with the theoretical PDF model to perform forecasts and eliminate the need for relying on expensive N-body simulations. Our theoretical model allows us to perform the first forecast for the convergence PDF that varies the full set of $\Lambda$CDM parameters. Our Fisher forecasts establish that the constraining power of the convergence PDF compares favourably to the two-point correlation function for a Euclid-like survey area at a single source redshift. When combined with a CMB prior from Planck, the PDF constrains both the neutrino mass $M_\nu$ and the dark energy equation of state $w_0$ more strongly than the two-point correlation function., Comment: Accepted by MNRAS

Since dark matter almost exclusively interacts gravitationally, the phase-space dynamics is described by the Vlasov-Poisson equation. A key characteristic is its infinite cumulant hierarchy, a tower of coupled evolution equations for the cumulants of the phase-space distribution. While on large scales the matter distribution is well described as a fluid and the hierarchy can be truncated, smaller scales are in the multi-stream regime in which all higher-order cumulants are sourced through nonlinear gravitational collapse. This regime is crucial for the formation of bound structures and the emergence of characteristic properties such as their density profiles. We present a novel closure strategy for the cumulant hierarchy that is inspired by finitely generated cumulants and hence beyond truncation. This constitutes a constructive approach for reducing nonlinear phase-space dynamics of Vlasov-Poisson to a closed system of equations in position space. Using this idea, we derive Schr\"odinger-Poisson as approximate quantal method for solving classical dynamics of Vlasov-Poisson with cold initial conditions. Our deduction complements the common reverse inference of the Schr\"odinger-Vlasov relation using a semi-classical limit of quantum mechanics and provides a clearer picture of the correspondence between classical and quantum dynamics. Our framework outlines an essential first step towards constructing approximate methods for Vlasov-like systems in cosmology and plasma physics with different initial conditions and potentials., Comment: 17 pages, 1 figure

Large deviation statistics is implemented to predict the statistics of cosmic densities in cylinders applicable to photometric surveys. It yields few percent accurate analytical predictions for the one-point probability distribution function (PDF) of densities in concentric or compensated cylinders; and also captures the density-dependence of their angular clustering (cylinder bias). All predictions are found to be in excellent agreement with the cosmological simulation Horizon Run 4 in the quasi-linear regime where standard perturbation theory normally breaks down. These results are combined with a simple local bias model that relates dark matter and tracer densities in cylinders and validated on simulated halo catalogues. This formalism can be used to probe cosmology with existing and upcoming photometric surveys like DES, Euclid or WFIRST containing billions of galaxies., Comment: 14 pages, 16 figures, comments welcome