Accurate cosmological simulations that include the effect of non-linear matter clustering as well as of massive neutrinos are essential for measuring the neutrino mass scale from upcoming galaxy surveys. Typically, Newtonian simulations are employed and the neutrino distribution is sampled with a large number of particles in order to beat down the shot noise. Here we perform very efficient simulations with light massive neutrinos which require virtually no extra cost over matter-only simulations, and furthermore do not require tracer particles for the neutrinos. Instead, we use a weak-field dictionary based on the recently developed Newtonion motion approach, where Newtonian simulations for matter are paired with a linear relativistic Boltzmann code to allow for an absorption of the neutrino evolution into a time-dependent coordinate transformation. For this, only minimal modifications in existing N-body codes are required, which we have explicitly implemented in $\textit{gevolution}$ and $\textit{gadget-2}$. Our fast method determines the non-linear matter power spectrum to permille-level precision when compared against state-of-the-art simulations that have been performed for $0.1\,{\rm eV} \leq \sum m_\nu \leq 0.3\,$eV.

Gravitational collapse of cold dark matter leads to infinite-density caustics that seed the primordial dark-matter halos in the large-scale structure. The development of these caustics begins, generically, as an almost one-dimensional phenomenon with the formation of pancakes. Focusing on the one-dimensional case, we identify a landscape of so far unknown singularities in the particle acceleration that emerge after the first crossing of particle trajectories. We complement our fully analytical studies by high-resolution N-body simulations and find outstanding agreement, particularly shortly after the first crossing. We develop the methods in 1D but outline briefly the necessary steps for the 3D case.

We chart a path toward solving for the nonlinear gravitational dynamics of cold dark matter by relying on a semiclassical description using the propagator. The evolution of the propagator is given by a Schr\"odinger equation, where the small parameter $\hbar$ acts as a softening scale that regulates singularities at shell-crossing. The leading-order propagator, called free propagator, is the semiclassical equivalent of the Zel'dovich approximation (ZA), that describes inertial particle motion along straight trajectories. At next-to-leading order, we solve for the propagator perturbatively and obtain, in the classical limit the displacement field from second-order Lagrangian perturbation theory (LPT). The associated velocity naturally includes an additional term that would be considered as third order in LPT. We show that this term is actually needed to preserve the underlying Hamiltonian structure, and ignoring it could lead to the spurious excitation of vorticity in certain implementations of second-order LPT. We show that for sufficiently small $\hbar$ the corresponding propagator solutions closely resemble LPT, with the additions that spurious vorticity is avoided and the dynamics at shell-crossing is regularised. Our analytical results possess a symplectic structure that allows us to advance numerical schemes for the large-scale structure. For times shortly after shell-crossing, we explore the generation of vorticity, which in our method does not involve any explicit multi-stream averaging, but instead arises naturally as a conserved topological charge.

We report the findings of new exact analytical solutions to the cosmological fluid equations, namely for the case where the initial conditions are perturbatively close to a spherical top-hat profile. To do so we enable a fluid description in a Lagrangian-coordinates approach, and prove the convergence of the Taylor-series representation of the Lagrangian displacement field until the time of collapse ("shell-crossing"). This allows the determination of the time for quasi-spherical collapse, which is shown to happen generically earlier than in the spherical case. For pedagogical reasons, calculations are first given for a spatially flat universe that is only filled with a non-relativistic component of cold dark matter (CDM). Then, the methodology is updated to a $\Lambda$CDM Universe, with the inclusion of a cosmological constant $\Lambda>0$.

Observables of cosmic structures are usually not the underlying matter field but biased tracers of matter, such as galaxies or halos. We show how the bias found in Newtonian N-body simulations can be interpreted in terms of the weak-field limit of General Relativity (GR). For this we employ standard Newtonian simulations of cold dark matter and incorporate GR/radiation via a weak-field dictionary that we have recently developed. We find that even when a simple local biasing scheme is employed in the Newtonian simulation, the relativistic bias becomes inherently scale-dependent due to the presence of radiation and GR corrections. This scale-dependence could be in principle observed on large scales in upcoming surveys. As a working example, we apply our methodology to Newtonian simulations for the spherical collapse and recover permille-level agreement between the approaches for extracting the relativistic bias on all considered scales.

We show how Newtonian cosmological simulations can be employed to investigate the non-linear evolution of two particle species in a relativistic context. We discuss the application for massive neutrinos and other multi-species systems such as Cold Dark Matter (CDM) plus baryons or Warm Dark Matter (WDM). We propose a method that allows us to perform simulations including massive neutrinos and general relativistic effects at almost the same computational cost as ordinary CDM only N-body simulations, employing tailor-made initial conditions and a dictionary for the interpretation of the simulation output.

We study the relativistic dynamics of a pressure-less and irrotational fluid of dark matter (CDM) with a cosmological constant ($\Lambda$), up to second order in cosmological perturbation theory. In our analysis we also account for vector and tensor perturbations and include primordial non-Gaussianity. We consider three gauges: the synchronous-comoving gauge, the Poisson gauge and the total matter gauge, where the first is the unique relativistic Lagrangian frame of reference, and the latters are convenient gauge choices for Eulerian frames. Our starting point is the metric and fluid variables in the Poisson gauge up to second order. We then perform the gauge transformations to the synchronous-comoving gauge and subsequently to the total matter gauge. Our expressions for the metrics, densities, velocities, and the gauge generators are novel and coincide with known results in the limit of a vanishing cosmological constant.

We show how standard Newtonian N-body simulations can be interpreted in terms of the weak-field limit of general relativity by employing the recently developed Newtonian motion gauge. Our framework allows the inclusion of radiation perturbations and the non-linear evolution of matter. We show how to construct the weak-field metric by combining Newtonian simulations with results from Einstein-Boltzmann codes. We discuss observational effects on weak lensing and ray tracing, identifying important relativistic corrections.

We consider the impact of neutrino self-interactions described by an effective four-fermion coupling on cosmological observations. Implementing the exact Boltzmann hierarchy for interacting neutrinos first derived in [arxiv:1409.1577] into the Boltzmann solver CLASS, we perform a detailed numerical analysis of the effects of the interaction on the cosmic microwave background (CMB) anisotropies, and compare our results with known approximations in the literature. While we find good agreement between our exact approach and the relaxation time approximation used in some recent studies, the popular $\left( c_{\text{eff}}^2,c_{\text{vis}}^2 \right)$-parameterisation fails to reproduce the correct scale dependence of the CMB temperature power spectrum. We then proceed to derive constraints on the effective coupling constant $G_{\text{eff}}$ using currently available cosmological data via an MCMC analysis. Interestingly, our results reveal a bimodal posterior distribution, where one mode represents the standard $\Lambda$CDM limit with $G_{\rm eff} \lesssim 10^8 \, G_{\rm F}$, and the other a scenario in which neutrinos self-interact with an effective coupling constant $G_{\rm eff} \simeq 3 \times 10^9 \, G_{\rm F}$.

Blow-up of solutions for the cosmological fluid equations, often dubbed shell-crossing or orbit crossing, denotes the breakdown of the single-stream regime of the cold-dark-matter fluid. At this instant, the velocity becomes multi-valued and the density singular. Shell-crossing is well understood in one dimension (1D), but not in higher dimensions. This paper is about quasi-one-dimensional (Q1D) flow that depends on all three coordinates but differs only slightly from a strictly 1D flow, thereby allowing a perturbative treatment of shell-crossing using the Euler--Poisson equations written in Lagrangian coordinates. The signature of shell-crossing is then just the vanishing of the Jacobian of the Lagrangian map, a regular perturbation problem. In essence the problem of the first shell-crossing, which is highly singular in Eulerian coordinates, has been desingularized by switching to Lagrangian coordinates, and can then be handled by perturbation theory. Here, all-order recursion relations are obtained for the time-Taylor coefficients of the displacement field, and it is shown that the Taylor series has an infinite radius of convergence. This allows the determination of the time and location of the first shell-crossing, which is generically shown to be taking place earlier than for the unperturbed 1D flow. The time variable used for these statements is not the cosmic time $t$ but the linear growth time $\tau \sim t^{2/3}$. For simplicity, calculations are restricted to an Einstein--de Sitter universe in the Newtonian approximation, and tailored initial data are used. However it is straightforward to relax these limitations, if needed.

Newtonian N-body simulations have been employed successfully over the past decades for the simulation of the cosmological large-scale structure. Such simulations usually ignore radiation perturbations (photons and massless neutrinos) and the impact of general relativity (GR) beyond the background expansion. This approximation can be relaxed and we discuss three different approaches that are accurate to leading order in GR. For simulations that start at redshift less than about 100 we find that the presence of early radiation typically leads to percent-level effects on the numerical power spectra at large scales. Our numerical results agree across the three methods, and we conclude that all of the three methods are suitable for simulations in a standard cosmology. Two of the methods modify the N-body evolution directly, while the third method can be applied as a post-processing prescription.

Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling method can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for LambdaCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple "N-body gauge" for z<50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable "forwards approach" for such cases.

What are the fundamental limitations of reconstructing the properties of dark energy, given cosmological observations in the quasilinear regime in a range of redshifts, to be as precise as required? The aim of this paper is to address this question by constructing model-independent observables, whilst completely ignoring practical problems of real-world observations. Non-Gaussianities already present in the initial conditions are not directly accessible from observations, because of a perfect degeneracy with the non-Gaussianities arising from the nonlinear matter evolution in generalized dark energy models. By imposing a specific set of evolution equations that should cover a range of dark energy cosmologies, we however find a constraint equation for the linear structure growth rate $f_1$ expressed in terms of model-independent observables. Entire classes of dark energy models which do not satisfy this constraint equation could be ruled out, and for models satisfying it we could reconstruct e.g. the nonlocal bias parameters $b_1$ and $b_2$.

Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian approximation, thus discarding general relativistic effects such as the coupling between matter and radiation ($\equiv$ photons and neutrinos). In this paper we investigate novel hybrid simulations which incorporate interactions between radiation and matter to the leading order in General Relativity, whilst evolving the matter dynamics in full non-linearity according to Newtonian theory. Our hybrid simulations come with a relativistic space-time and make it possible to investigate structure formation in a unified framework. In the present work we focus on simulations initialized at $z=99$, and show that the extracted matter power spectrum receives up to $3\%$ corrections on very large scales through radiation. Our numerical findings compare favourably with linear analytical results from Fidler et al. (2016), from which we deduce that there cannot be any significant non-linear mode-coupling induced through linear radiation corrections.

We study the nonlinear gravitational dynamics of a universe filled with a pressureless fluid and a cosmological constant $\Lambda$ in the context of Newtonian gravity, and in the relativistic post-Friedmann approach proposed in paper I [I. Milillo et al., Phys. Rev. D 92, 023519 (2015).]. The post-Friedmann approximation scheme is based on the $1/c$ expansion of the space-time metric and the energy-momentum tensor, and includes nonlinear Newtonian cosmology. Here we establish the nonlinear post-Friedmann framework in the Lagrangian-coordinates approach for structure formation. For this we first identify a Lagrangian gauge which is suitable for incorporating nonzero vorticity. We analyze our results in two limits: at the leading order we recover the fully nonlinear Newtonian cosmological equations in the Lagrangian formulation, and we provide a space-time metric consistent from the perspective of general relativity. We then linearize our expressions and recover the relativistic results at first order in cosmological perturbation theory. Therefore, the introduced approximation scheme provides a unified treatment for the two leading-order regimes, from the small scales described by Newtonian gravity to the large linear scale, where first-order relativistic cosmological perturbation theory gives a very good description of structure formation.

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

Certain inflationary models as well as realisations of phase transitions in the early Universe predict the formation of primordial black holes. For most mass ranges, the fraction of matter in the form of primordial black holes is limited by many different observations on various scales. Primordial black holes are assumed to be formed when overdensities that cross the horizon have Schwarzschild radii larger than the horizon. Traditionally it was therefore assumed that primordial black-hole masses were equal to the horizon mass at their time of formation. However, detailed calculations of their collapse show that primordial black holes formed at each point in time should rather form a spectrum of different masses, obeying critical scaling. Though this has been known for more than fifteen years, the effect of this scaling behaviour is largely ignored when considering predictions for primordial black hole mass spectra. In this paper we consider the critical collapse scaling for a variety of models which produce primordial black holes, and find that it generally leads to a shift, broadening and an overall decrease of the mass contained in primordial black holes. This effect is model and parameter dependent and cannot be contained by a constant rescaling of the spectrum; it can become important and should be taken into account when comparing to observational constraints.

We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the synchronous-comoving gauge, allows us to go beyond the Newtonian predictions and to calculate relativistic corrections to it. We use an Einstein--de Sitter model, including both growing and decaying modes in our solutions. The impact of our results should be assessed through the implementation of the featured displacement in cosmological $N$-body simulations.

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

It is shown here that in a flat, cold dark matter (CDM) dominated Universe with positive cosmological constant ($\Lambda$), modelled in terms of a Newtonian and collisionless fluid, particle trajectories are analytical in time (representable by a convergent Taylor series) until at least a finite time after decoupling. The time variable used for this statement is the cosmic scale factor, i.e., the "$a$-time", and not the cosmic time. For this, a Lagrangian-coordinates formulation of the Euler-Poisson equations is employed, originally used by Cauchy for 3-D incompressible flow. Temporal analyticity for $\Lambda$CDM is found to be a consequence of novel explicit all-order recursion relations for the $a$-time Taylor coefficients of the Lagrangian displacement field, from which we derive the convergence of the $a$-time Taylor series. A lower bound for the $a$-time where analyticity is guaranteed and shell-crossing is ruled out is obtained, whose value depends only on $\Lambda$ and on the initial spatial smoothness of the density field. The largest time interval is achieved when $\Lambda$ vanishes, i.e., for an Einstein-de Sitter universe. Analyticity holds also if, instead of the $a$-time, one uses the linear structure growth $D$-time, but no simple recursion relations are then obtained. The analyticity result also holds when a curvature term is included in the Friedmann equation for the background, but inclusion of a radiation term arising from the primordial era spoils analyticity.

Starting from the collisional Boltzmann equation, we derive for the first time and from first principles the Boltzmann hierarchy for neutrinos including interactions with a scalar particle. Such interactions appear, for example, in majoron-like models of neutrino mass generation. We study two limits of the scalar mass: (i) An extremely massive scalar whose only role is to mediate an effective 4-fermion neutrino-neutrino interaction, and (ii) a massless scalar that can be produced in abundance and thus demands its own Boltzmann hierarchy. In contrast to, e.g., the first-order Boltzmann hierarchy for Thomson-scattering photons, our interacting neutrino/scalar Boltzmann hierarchies contain additional momentum-dependent collision terms arising from a non-negligible energy transfer in the neutrino-neutrino and neutrino-scalar interactions. This necessitates that we track each momentum mode of the phase space distributions individually, even if the particles were massless. Comparing our hierarchy with the commonly used $(c_{\rm eff}^2,c_{\rm vis}^2)$-parameterisation, we find no formal correspondence between the two approaches, which raises the question of whether the latter parameterisation even has an interpretation in terms of particle scattering. Lastly, although we have invoked majoron-like models as a motivation for our study, our treatment is in fact generally applicable to all scenarios in which the neutrino and/or other ultrarelativistic fermions interact with scalar particles.

We investigate the purely spatial Lagrangian coordinate transformation from the Lagrangian to the basic Eulerian frame. We demonstrate three techniques for extracting the relativistic displacement field from a given solution in the Lagrangian frame. These techniques are (a) from defining a local set of Eulerian coordinates embedded into the Lagrangian frame; (b) from performing a specific gauge transformation; and (c) from a fully non-perturbative approach based on the ADM split. The latter approach shows that this decomposition is not tied to a specific perturbative formulation for the solution of the Einstein equations. Rather, it can be defined at the level of the non-perturbative coordinate change from the Lagrangian to the Eulerian description. Studying such different techniques is useful because it allows us to compare and develop further the various approximation techniques available in the Lagrangian formulation. We find that one has to solve the gravitational wave equation in the relativistic analysis, otherwise the corresponding Newtonian limit will necessarily contain spurious non-propagating tensor artefacts at second order in the Eulerian frame. We also derive the magnetic part of the Weyl tensor in the Lagrangian frame, and find that it is not only excited by gravitational waves but also by tensor perturbations which are induced through the non-linear frame-dragging. We apply our findings to calculate for the first time the relativistic displacement field, up to second order, for a $\Lambda$CDM Universe in the presence of a local primordial non-Gaussian component. Finally, we also comment on recent claims about whether mass conservation in the Lagrangian frame is violated.

We discuss the relation between the output of Newtonian N-body simulations on scales that approach or exceed the particle horizon to the description of General Relativity. At leading order, the Zeldovich approximation is correct on large scales, coinciding with the General Relativistic result. At second order in the initial metric potential, the trajectories of particles deviate from the second order Newtonian result and hence the validity of 2LPT initial conditions should be reassessed when used in very large simulations. We also advocate using the expression for the synchronous gauge density as a well behaved measure of density fluctuations on such scales.

We investigate gravitational analogue models to describe slowly rotating objects (e.g., dark-matter halos, or boson stars) in terms of Bose-Einstein condensates, trapped in their own gravitational potentials. We begin with a modified Gross-Pitaevskii equation, and show that the resulting background equations of motion are stable, as long as the rotational component is treated as a small perturbation. The dynamics of the fluctuations of the velocity potential are effectively governed by the Klein-Gordon equation of a "Eulerian metric," where we derive the latter by the use of a relativistic Lagrangian extrapolation. Superradiant scattering on such objects is studied. We derive conditions for its occurence and estimate its strength. Our investigations might give an observational handle to phenomenologically constrain Bose-Einstein condensates.

The physical interpretation of cold dark matter perturbations is clarified by associating Bertschinger's Poisson gauge with a Eulerian/observer's frame of reference. We obtain such an association by using a Lagrangian approach to relativistic cosmological structure formation. Explicitly, we begin with the second-order solution of the Einstein equations in a synchronous/comoving coordinate system---which defines the Lagrangian frame, and transform it to a Poissonian coordinate system. The generating vector of this coordinate/gauge transformation is found to be the relativistic displacement field. The metric perturbations in the Poissonian coordinate system contain known results from standard/Eulerian Newtonian perturbation theory, but contain also purely relativistic corrections. On sub-horizon scales these relativistic corrections are dominated by the Newtonian bulk part. These corrections however set up non-linear constraints for the density and for the velocity which become important on scales close to the horizon. Furthermore, we report the occurence of a transverse component in the displacement field, and find that it induces a non-linear frame dragging as seen in the observer's frame, which is sub-dominant at late-times and sub-horizon scales. Finally, we find two other gauges which can be associated with a Eulerian frame. We argue that the Poisson gauge is to be preferred because it comes with the simplest physical interpretation.

We derive a recursion relation in the framework of Lagrangian perturbation theory, appropriate for studying the inhomogeneities of the large scale structure of the universe. We use the fact that the perturbative expansion of the matter density contrast is in one-to-one correspondence with standard perturbation theory (SPT) at any order. This correspondence has been recently shown to be valid up to fourth order for a non-relativistic, irrotational and dust-like component. Assuming it to be valid at arbitrary (higher) order, we express the Lagrangian displacement field in terms of the perturbative kernels of SPT, which are itself given by their own and well-known recursion relation. We argue that the Lagrangian solution always contains more non-linear information in comparison with the SPT solution, (mainly) if the non-perturbative density contrast is restored after the displacement field is obtained.

We describe the irrotational dust component of the universe in terms of a relativistic gradient expansion and transform the resulting synchronous metric to a Newtonian coordinate system. The two metrics are connected via a space-like displacement field and a time-like perturbation, providing a relativistic generalization of the transformation from Lagrangian to Eulerian coordinates. The relativistic part of the displacement field generates already at initial time a non-local density perturbation at second order. This is a purely relativistic effect since it originates from space-time mixing. We give two options, the passive and the active approach, on how to include the relativistic corrections for example in N-body simulations. In the passive approach we treat the corrections as a non-Gaussian modification of the initial Gaussian field (primordial non-Gaussianity could be incorporated as well). The induced non-Gaussianity depends on scale and the redshift at which initial conditions are set, with f_NL ~ few for small enough scales and redshifts. In the active approach we show how to use the relativistic trajectory to obtain the initial displacement and velocity of particles for N-body simulations without modifying the initial Gaussian field.

We show how the Zel'dovich approximation and the second order displacement field of Lagrangian perturbation theory can be obtained from a general relativistic gradient expansion in \Lambda{}CDM cosmology. The displacement field arises as a result of a second order non-local coordinate transformation which brings the synchronous/comoving metric into a Newtonian form. We find that, with a small modification, the Zel'dovich approximation holds even on scales comparable to the horizon. The corresponding density perturbation is not related to the Newtonian potential via the usual Poisson equation but via a modified Helmholtz equation. This is a consequence of causality not present in the Newtonian theory. The second order displacement field receives relativistic corrections that are subdominant on short scales but are comparable to the second order Newtonian result on scales approaching the horizon. The corrections are easy to include when setting up initial conditions in large N-body simulations.

This is part two in a series of papers in which we investigate an approach based on Lagrangian perturbation theory (LPT) to study the non-linear evolution of the large-scale structure distribution in the universe. Firstly, we compute the matter bispectrum in real space using LPT up one-loop order, for both Gaussian and non-Gaussian initial conditions. In the initial position limit, we find that the one-loop bispectrum computed in this manner is identical to its counterpart obtained from standard Eulerian perturbation theory (SPT). Furthermore, the LPT formalism allows for a simple reorganisation of the perturbative series corresponding to the resummation of an infinite series of perturbations in SPT. Applying this method, we find a resummed one-loop bispectrum that compares favourably with results from N-body simulations. We generalise the resummation method also to the computation of the redshift-space bispectrum up to one loop.

We investigate the Lagrangian perturbation theory of a homogeneous and isotropic universe in the non-relativistic limit, and derive the solutions up to the fourth order. These solutions are needed for example for the next-to-leading order correction of the (resummed) Lagrangian matter bispectrum, which we study in an accompanying paper. We focus on flat cosmologies with a vanishing cosmological constant, and provide an in-depth description of two complementary approaches used in the current literature. Both approaches are solved with two different sets of initial conditions---both appropriate for modelling the large-scale structure. Afterwards we consider only the fastest growing mode solution, which is not affected by either of these choices of initial conditions. Under the reasonable approximation that the linear density contrast is evaluated at the initial Lagrangian position of the fluid particle, we obtain the nth-order displacement field in the so-called initial position limit: the nth order displacement field consists of 3(n-1) integrals over n linear density contrasts, and obeys self-similarity. Then, we find exact relations between the series in Lagrangian and Eulerian perturbation theory, leading to identical predictions for the density contrast and the peculiar-velocity divergence up to the fourth order.

The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape information currently contained in the SDSS DR7 halo power spectrum (HPS) is by far superseded by geometric information derived from the baryonic features. An immediate corollary is that contrary to popular beliefs, the upper limit on the neutrino mass m_\nu presently derived from LSS combined with cosmic microwave background (CMB) data does not in fact arise from the possible small-scale power suppression due to neutrino free-streaming, if we limit the model framework to minimal LambdaCDM+m_\nu. However, in more complicated models, such as those extended with extra light degrees of freedom and a dark energy equation of state parameter w differing from -1, shape information becomes crucial for the resolution of parameter degeneracies. This conclusion will remain true even when data from the Planck surveyor become available. In the course of our analysis, we introduce a new dewiggling procedure that allows us to extend consistently the use of the SDSS HPS to models with an arbitrary sound horizon at decoupling. All the cases considered here are compatible with the conservative 95%-bounds \sum m_\nu < 1.16 eV, N_eff = 4.8 \pm 2.0.