<figure><img src="image_1.png" /></figure>

HAL
open science

# One-point statistics of the Lagrangian displacement field; Addendum to: one-point remapping of Lagrangian perturbation theory in the mildly non-linear regime of cosmic structure formation

Florent Leclercq, Jens Jasche, Benjamin Wandelt

<figure><img src="image_2.png" /></figure>

## ► To cite this version:

Florent Leclercq, Jens Jasche, Benjamin Wandelt. One-point statistics of the Lagrangian displacement field; Addendum to: one-point remapping of Lagrangian perturbation theory in the mildly non-linear regime of cosmic structure formation. Journal of Cosmology and Astroparticle Physics, 2015, 2015 (04), pp.026-026. 10.1088/1475-7516/2015/04/A01. hal-05347508

HAL Id: hal-05347508

https://hal.science/hal-05347508v1

Submitted on 5 Nov 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_4.png" /></figure>

OPEN ACCESS

Addendum: one-point remapping of Lagrangian
perturbation theory in the mildly non-linear regime
of cosmic structure formation

To cite this article: Florent Leclercq et al JCAP04(2015)026

View the [article online](#) for updates and enhancements.

You may also like

- The dynamics of neutron star crusts:  
  *Lagrangian perturbation theory for a relativistic superfluid-elastic system*  
  N Andersson, B Haskell, G L Comer et al.

- A Lagrangian perturbation theory in the presence of massive neutrinos
  Alejandro Aviles and Arka Banerjee

- *Lagrangian Perturbation Theory for Biased Tracers: Significance of the Number Conservation*
  Peter Espenshade and Jaiyul Yoo
Addendum: one-point remapping of
Lagrangian perturbation theory in the
mildly non-linear regime of cosmic
structure formation

Florent Leclercq$^{a,b,c}$ Jens Jasche$^d$ and Benjamin Wandelt$^{a,b,e,f}$

E-mail: florent.leclercq@polytechnique.org, jasche@iap.fr, wandelt@iap.fr

Received March 28, 2015

Accepted March 30, 2015

Published April 20, 2015

**Addendum to:** JCAP11(2013)048

**ArXiv ePrint:** 1305.4642

<figure><img src="image_5.png" /></figure>

<footnote><sup>a</sup>Institut d'Astrophysique de Paris (IAP),
UMR 7095, CNRS — UPMC Université Paris 6, Sorbonne Universités,
98bis boulevard Arago, F-75014 Paris, France</footnote>

<footnote><sup>b</sup>Institut Lagrange de Paris (ILP), Sorbonne Universités,
98bis boulevard Arago, F-75014 Paris, France</footnote>

<footnote><sup>c</sup>École polytechnique ParisTech, Route de Saclay,
F-91128 Palaiseau, France</footnote>

<footnote><sup>d</sup>Excellence Cluster Universe, Technische Universität München,
Boltzmannstrasse 2, 85748 Garching, Germany</footnote>

<footnote><sup>e</sup>Department of Physics, University of Illinois at Urbana-Champaign,
Urbana, IL 61801, U.S.A.</footnote>

<footnote><sup>f</sup>Department of Astronomy, University of Illinois at Urbana-Champaign,
Urbana, IL 61801, U.S.A.</footnote>
**One-point statistics of the Lagrangian displacement field.** The remapping procedure, described in section 2 of the main paper, relies on the Eulerian density contrast. As noted by previous authors [see in particular 1], in the Lagrangian representation of the LSS, it is natural to use the divergence of the displacement field $\psi$ instead of the Eulerian density contrast $\delta$. This addendum provides additional comments on the one-point statistics of $\psi$ and comparatively analyzes key features of $\psi$ and $\delta$.

In the Lagrangian frame, the quantity of interest is not the position, but the displacement field $\Psi(\mathbf{q})$ which maps the initial comoving particle position $\mathbf{q}$ to its final comoving Eulerian position $\mathbf{x}$ (see e.g. [2] or [3] for overviews),

$$ \mathbf{x} = \mathbf{q} + \Psi(\mathbf{q}). \qquad (1) $$

It is important to note that, though $\Psi(\mathbf{q})$ is a priori a full three-dimensional vector field, it is curl-free up to second order in Lagrangian perturbation theory (appendix D in [4] or [3] for a review). We did not consider contributions beyond 2LPT. After publication of the main paper, Chan [5] analyzed the non-linear evolution of $\Psi$, splitting it into its scalar and vector parts (the so-called “Helmholtz decomposition”). Looking at two-point statistics, he found that shell-crossing leads to a suppression of small-scale power in the scalar part, and, subdominantly, to the generation of a vector contribution.

Let $\psi(\mathbf{q}) = \nabla_{\mathbf{q}} \cdot \Psi(\mathbf{q})$ denote the divergence of the displacement field, where $\nabla_{\mathbf{q}}$ is the divergence operator in Lagrangian coordinates. $\psi$ quantifies the angle-averaged spatial-stretching of the Lagrangian dark matter “sheet” in comoving coordinates [1]. Let $\mathcal{P}_{\psi,\text{LPT}}$ and $\mathcal{P}_{\psi,\text{Nbody}}$ be the one-point probability distribution functions (PDF) for the divergence of the displacement field in the LPT and in the full N-body fields, respectively. We denote by $\mathcal{P}_{\delta}$ the corresponding PDFs for the Eulerian density contrast.

In figure 1, we show the PDFs of $\psi$ for the ZA, 2LPT and full N-body gravity. The most important feature of $\psi$ is that, whatever the model for structure formation, the PDF exhibits reduced non-Gaussianity compared to the PDF for the density contrast $\delta$ (see the upper panel of figure 7 in the main paper for comparison). The main reason is that $\mathcal{P}_{\delta}$, unlike $\mathcal{P}_{\psi}$, is tied down to zero at $\delta = -1$. It is highly non-Gaussian in the final conditions, both in N-body simulations and in approximations to the true dynamics. For a quantitative analysis, we looked at the first and second-order non-Gaussianity statistics: the skewness $\gamma_1$ and the excess kurtosis $\gamma_2$,

$$ \gamma_1 = \frac{\mu_3}{\sigma^3} \quad \text{and} \quad \gamma_2 = \frac{\mu_4}{\sigma^4} - 3, \qquad (2) $$

where $\mu_n$ is the n-th moment about the mean and $\sigma$ is the standard deviation. We estimated $\gamma_1$ and $\gamma_2$ at redshift zero in our simulations, in the one-point statistics of the density contrast $\delta$ and of the divergence of the displacement field $\psi$. The results are shown in table 1. In all cases, we found that both $\gamma_1$ and $\gamma_2$ are much smaller when measured from $\mathcal{P}_{\psi}$ instead of $\mathcal{P}_{\delta}$.

At linear order in Lagrangian perturbation theory (the Zel'dovich approximation), the divergence of the displacement field is proportional to the density contrast in the initial conditions, $\delta(\mathbf{q})$, scaling with the negative growth factor, $-D_1(\tau)$:

$$ \psi^{(1)}(\mathbf{q}, \tau) = \nabla_{\mathbf{q}} \cdot \Psi^{(1)}(\mathbf{q}, \tau) = -D_1(\tau)\delta(\mathbf{q}). \qquad (3) $$

Since we take Gaussian initial conditions, the PDF for $\psi$ is Gaussian at any time with the ZA. In full gravity, non-linear evolution slightly breaks Gaussianity. $\mathcal{P}_{\psi,\text{Nbody}}$ is slightly skewed towards negative values while its mode gets shifted around $\psi \approx 1$. Taking into account non-local effects, 2LPT tries to get closer to the shape observed in N-body simulations, but the correct skewness is overshot and the PDF is exceedingly peaked.
<figure><img src="image_6.png" /><figcaption>**Figure 1.** Redshift-zero probability distribution function for the divergence of the displacement field $\psi$, computed from eight 1024 Mpc/h-box simulations of 512³ particles. A quantitative analysis of the deviation from Gaussianity of these PDFs is given in table 1. The particle distribution is determined using: a full N-body simulation (purple curve), the Zel'dovich approximation (ZA, light red curve) and second-order Lagrangian perturbation theory (2LPT, light blue curve). The vertical line at $\psi = -3$ represents the collapse barrier about which $\psi$ values bob around after gravitational collapse. A bump at this value is visible with full gravity, but LPT is unable to reproduce this feature. This regime corresponds to virialized, overdense clusters.</figcaption></figure>

<table><caption>**Table 1.** Non-Gaussianity parameters (the skewness $\gamma_1$ and the excess kurtosis $\gamma_2$) of the redshift-zero probability distribution functions $P_\delta$ and $P_\psi$ of the density contrast $\delta$ and the divergence of the displacement field $\psi$, respectively. The confidence intervals given correspond to the 1-$\sigma$ standard deviations among eight realizations. In all cases, $\gamma_1$ and $\gamma_2$ are reduced when measured from $\psi$ instead of $\delta$.</caption><thead><tr><th>Model</th><th>P<sub>δ</sub></th><th>P<sub>ψ</sub></th></tr></thead><tbody><tr><td colspan="3" style="text-align: center;">Skewness γ<sub>1</sub></td></tr><tr><td>ZA</td><td>2.36 ± 0.01</td><td>-0.0067 ± 0.0001</td></tr><tr><td>2LPT</td><td>2.83 ± 0.01</td><td>-1.5750 ± 0.0002</td></tr><tr><td>N-body</td><td>5.14 ± 0.05</td><td>-0.4274 ± 0.0001</td></tr><tr><td colspan="3" style="text-align: center;">Excess kurtosis γ<sub>2</sub></td></tr><tr><td>ZA</td><td>9.95 ± 0.09</td><td>-2.2154 × 10<sup>-6</sup> ± 0.0003</td></tr><tr><td>2LPT</td><td>13.91 ± 0.15</td><td>3.544 ± 0.0011</td></tr><tr><td>N-body</td><td>62.60 ± 2.75</td><td>-0.2778 ± 0.0004</td></tr></tbody></table>

Figure 2 shows a slice of the divergence of the displacement field, measured at redshift zero for particles occupying a flat 512²-pixel Lagrangian sheet from one of our simulations. For comparison, see also the figures in [1, 6, 7]. We used the color scheme of the latter paper, suggesting a topographical analogy when working in Lagrangian coordinates. As structures take shape, $\psi$ departs from its initial value; it takes positive values in underdensities and negative values in overdensities. The shape of voids (the “mountains”) is found to be reasonably similar in LPT and in the N-body simulation. For this reason, the influence of late-time non-linear effects in voids is milder as compared to overdense structures, which makes them easier to relate to the initial conditions. However, in overdense regions where $\psi$ decreases, it is not allowed to take arbitrary values: where gravitational collapse occurs,
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Slices of the divergence of the displacement field, $\psi$, on a Lagrangian sheet of 512$^2$ particles from a 512$^3$-particle simulation of box size 1024 Mpc/h, run to redshift zero. For clarity we show only a 200 Mpc/h region. Each pixel corresponds to a particle. The particle distribution is determined using respectively a full $N$-body simulation, the Zel'dovich approximation (ZA) and second-order Lagrangian perturbation theory (2LPT). In the upper left panel, the density contrast $\delta$ in the $N$-body simulation is shown, after binning on a 512$^3$-voxel grid. To guide the eye, some clusters and voids are identified by yellow and purple dots, respectively. The “lakes”, Lagrangian regions that have collapsed to form halos, are only visible in the $N$-body simulation, while the “mountains”, Lagrangian regions corresponding to cosmic voids, are well reproduced by LPT.</figcaption></figure>

"lakes" form and $\psi$ gets stuck around a collapse barrier, $\psi \approx -3$. As expected, these "lakes", corresponding to virialized clusters, can only be found in $N$-body simulations, since LPT fails to accurately describe the highly non-linear physics involved. A small bump at $\psi = -3$ is visible in $P_{\psi,Nbody}$ (see figure 1). We checked that this bump gets more visible in higher mass-resolution simulations (200 Mpc/h box for 256$^3$ particles), where matter is more clustered. This means that part of the information about gravitational clustering can be found in the one-point statistics of $\psi$. Of course, the complete description of halos requires to precisely account for the shape of the "lakes", which can only be done via higher-order correlation functions. More generally, it is possible to use Lagrangian information in order to classify structures of the cosmic web. In particular, DIVA [8] uses the shear of the displacement field and ORIGAMI [9] the number of phase-space folds. As pointed out by Falck & Neyrinck [10], while these techniques cannot be straightforwardly used for the analysis of galaxy surveys, where we lack Lagrangian information, recently proposed techniques for physical inference of the initial conditions [11, 12] should allow their use with observational data.
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Left panel. Two-dimensional histograms comparing particle densities evolved with full $N$-body dynamics (the $x$-axis) to densities in the LPT-evolved particle distributions (the $y$-axis). The red lines show the ideal $y = x$ locus. A turn-up at low densities is visible with 2LPT, meaning that some overdense regions are predicted where there should be deep voids. Right panel. Same plot for the divergence of the displacement field $\psi$. Negative $\psi$ corresponds to overdensities and positive $\psi$ correspond to underdensities. The dotted blue line shows the collapse barrier at $\psi = -3$ where particle get clustered in full gravity. The scatter is bigger with $\psi$ than with $\delta$, in particular in overdensities, since with LPT, particles do not cluster. The turn-up at low densities with 2LPT, observed with the density contrast, is also visible with the divergence of the displacement field.</figcaption></figure>

Figure 3 shows two-dimensional histograms comparing $N$-body simulations to the LPT realizations for the density contrast $\delta$ and the divergence of the displacement field $\psi$. At this point, it is useful to note that a good mapping exists in the case where the relation shown is monotonic and the scatter is narrow. As pointed out by Neyrinck [1, 13], matter in the substructure of 2LPT-voids has incorrect statistical properties: there are overdense particles in the low density region of the 2LPT $\delta$-scatter plot. This degeneracy is also visible in the $\psi > 0$ region of the 2LPT $\psi$-scatter plot. On average, the scatter is bigger with $\psi$ than with $\delta$, in particular in overdensities ($\psi < 0$), since with LPT, particles do not cluster: $\psi$ takes any value between 2 and $-3$ where it should remain around $-3$.

Summing up our discussions in this addendum, we analyzed the relative merits of the Lagrangian divergence of the displacement field $\psi$, and the Eulerian density contrast $\delta$. The important differences are the following:

1. $\Psi$ being irrotational up to order two, its divergence $\psi$ contains nearly all information on the displacement field in one dimension, instead of three. The collapse barrier at $\psi = -3$ is visible in $P_{\psi}$ for $N$-body simulations but not for LPT. A part of the information about non-linear gravitational clustering is therefore encoded in the one-point statistics of $\psi$.
2. $\psi$ exhibits much fewer gravitationally-induced non-Gaussian features than $\delta$ in the final conditions (figure 1 and table 1).

3. However, the values of $\psi$ are more scattered than the values of $\delta$ with respect to the true dynamics (figure 3), meaning that an unambiguous mapping is more difficult.

## Acknowledgments

FL acknowledges funding from an AMX grant (École polytechnique ParisTech). JJ is partially supported by a Feodor Lynen Fellowship by the Alexander von Humboldt foundation. BW acknowledges funding from an ANR Chaire d'Excellence (ANR-10-CEXC-004-01) and the UPMC Chaire Internationale in Theoretical Cosmology. This work has been done within the Labex Institut Lagrange de Paris (reference ANR-10-LABX-63) part of the Idex SUPER, and received financial state aid managed by the Agence Nationale de la Recherche, as part of the programme Investissements d'avenir under the reference ANR-11-IDEX-0004-02. This research was supported by the DFG cluster of excellence "Origin and Structure of the Universe".

## References

[1] M.C. Neyrinck, *Quantifying distortions of the Lagrangian dark-matter mesh in cosmology*, Mon. Not. Roy. Astron. Soc. **428** (2013) 141 [arXiv:1204.1326] [INSPIRE].

[2] F.R. Bouchet, S. Colombi, E. Hivon and R. Juszkiewicz, *Perturbative Lagrangian approach to gravitational instability*, Astron. Astrophys. **296** (1995) 575 [astro-ph/9406013] [INSPIRE].

[3] F. Bernardeau, S. Colombi, E. Gaztanaga and R. Scoccimarro, *Large scale structure of the universe and cosmological perturbation theory*, Phys. Rept. **367** (2002) 1 [astro-ph/0112551] [INSPIRE].

[4] F. Bernardeau, *The nonlinear evolution of rare events*, Astrophys. J. **427** (1994) 51 [astro-ph/9311066] [INSPIRE].

[5] K.C. Chan, *Helmholtz decomposition of the lagrangian displacement*, Phys. Rev. D **89** (2014) 083515 [arXiv:1309.2243] [INSPIRE].

[6] R. Mohayaee, H. Mathis, S. Colombi and J. Silk, *Reconstruction of primordial density fields*, Mon. Not. Roy. Astron. Soc. **365** (2006) 939 [astro-ph/0501217] [INSPIRE].

[7] S. Pueblas and R. Scoccimarro, *Generation of vorticity and velocity dispersion by orbit crossing*, Phys. Rev. D **80** (2009) 043504 [arXiv:0809.4606] [INSPIRE].

[8] G. Lavaux and B.D. Wandelt, *Precision cosmology with voids: definition, methods, dynamics*, Mon. Not. Roy. Astron. Soc. **403** (2010) 1392 [arXiv:0906.4101] [INSPIRE].

[9] B.L. Falck, M.C. Neyrinck and A.S. Szalay, *ORIGAMI: delineating halos using phase-space folds*, Astrophys. J. **754** (2012) 126 [arXiv:1201.2353] [INSPIRE].

[10] B. Falck and M.C. Neyrinck, *The persistent percolation of single-stream voids*, arXiv:1410.4751 [INSPIRE].

[11] J. Jasche and B.D. Wandelt, *Bayesian physical reconstruction of initial conditions from large scale structure surveys*, Mon. Not. Roy. Astron. Soc. **432** (2013) 894 [arXiv:1203.3639] [INSPIRE].

[12] J. Jasche, F. Leclercq and B.D. Wandelt, *Past and present cosmic structure in the SDSS DR7 main sample*, JCAP **01** (2015) 036 [arXiv:1409.6308] [INSPIRE].

[13] V. Sahni and S. Shandarin, *Accuracy of Lagrangian approximations in voids*, Mon. Not. Roy. Astron. Soc. **282** (1996) 641.