Search

Cosmological perturbation theory is known to converge poorly for predicting the spherical collapse and void evolution of collisionless matter. Using the exact parametric solution as a testing ground, we develop two asymptotic methods in spherical symmetry that resolve the gravitational evolution to much higher accuracy than Lagrangian perturbation theory (LPT), which is the current gold standard in the literature. One of the methods selects a stable fixed-point solution of the renormalization-group flow equation, thereby predicting already at the leading order the critical exponent of the phase transition to collapsed structures. The other method completes the truncated LPT series far into the UV regime, by adding a non-analytic term that captures the critical nature of the gravitational collapse. We find that the UV method most accurately resolves the evolution of the nonlinear density as well as its one-point probability distribution function. Similarly accurate predictions are achieved with the renormalization-group method, especially when paired with Pad\'e approximants. Further, our results yield new, very accurate, formulas to relate linear and nonlinear density contrasts. Finally, we chart possible ways on how to adapt our methods to the case of cosmological random field initial conditions., Comment: 11 pages + appendices and references, 10 figures, v3: updated Figs. 1,2,3,5 now showing sign changes in the particle positions after collapse; matches published version

The cosmological fluid equations describe the early gravitational dynamics of cold dark matter (CDM), exposed to a uniform component of dark energy, the cosmological constant $\Lambda$. Perturbative predictions for the fluid equations typically assume that the impact of $\Lambda$ on CDM can be encapsulated by a refined growing mode $D$ of linear density fluctuations. Here we solve, to arbitrary high perturbative orders, the nonlinear fluid equations with an {\it Ansatz} for the fluid variables in increasing powers of $D$. We show that $\Lambda$ begins to populate the solutions starting at the fifth order in this strict $D$-expansion. By applying suitable resummation techniques, we recast these solutions to a standard perturbative series where not $D$, but essentially the initial gravitational potential serves as the bookkeeping parameter within the expansion. Then, by using the refined growth functions at second and third order in standard perturbation theory, we determine the matter power spectrum to one-loop accuracy as well as the leading-order contribution to the matter bispectrum. We find that employing our refined growth functions impacts the total power- and bispectra at a precision that is below one percent at late times. However, for the power spectrum, we find a characteristic scale-dependent suppression that is fairly similar to what is observed in massive neutrino cosmologies. Therefore, we recommend employing our refined growth functions in order to reduce theoretical uncertainties for analysing data in related pipelines., Comment: v1: 10 pages, v2: 11 pages, figure 3 shows also results for z=1, added references, MNRAS accepted

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 ΛCDM Universe, with the inclusion of a cosmological constant Λ > 0.

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., Comment: v2, matches published version at JCAP

Inaccuracies in the initial conditions for cosmological N-body simulations could easily be the largest source of systematic error in predicting the non-linear large-scale structure. From the theory side, initial conditions are usually provided by using low-order truncations of the displacement field from Lagrangian perturbation theory, with the first and second-order approximations being the most common ones. Here we investigate the improvement brought by using initial conditions based on third-order Lagrangian perturbation theory (3LPT). We show that with 3LPT, truncation errors are vastly suppressed, thereby opening the portal to initializing simulations accurately as late as z=12 (for the resolution we consider). We analyse the competing effects of perturbative truncation and particle discreteness on various summary statistics. Discreteness errors are essentially decaying modes and thus get strongly amplified for earlier initialization times. We show that late starting times with 3LPT provide the most accurate configuration, which we find to coincide with the continuum fluid limit within 1 per cent for the power- and bispectrum at z=0 up to the particle Nyquist wave number of our simulations (k ~ 3h/Mpc). In conclusion, to suppress non-fluid artefacts, we recommend initializing simulations as late as possible with 3LPT. We make our 3LPT initial condition generator publicly available., Comment: 22 pages, 13 figures, minor revisions (matches published version), code available from https://bitbucket.org/ohahn/monofonic

Perturbation theory is an indispensable tool for studying the cosmic large-scale structure, and establishing its limits is therefore of utmost importance. One crucial limitation of perturbation theory is shell-crossing, which is the instance when cold-dark-matter trajectories intersect for the first time. We investigate Lagrangian perturbation theory (LPT) at very high orders in the vicinity of the first shell-crossing for random initial data in a realistic three-dimensional Universe. For this, we have numerically implemented the all-order recursion relations for the matter trajectories, from which the convergence of the LPT series at shell-crossing is established. Convergence studies performed at large orders reveal the nature of the convergence-limiting singularities. These singularities are not the well-known density singularities at shell-crossing but occur at later times when LPT already ceased to provide physically meaningful results.

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., Comment: 28 pages, 9 figures

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., Comment: 18 pages, added figure 6 that displays the generation of vorticity, matches published version at PRD

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., Comment: minors revisions, published at JCAP

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., 22 pages, 3 figures; accepted for publication at JCAP

The present is a companion paper to "A contemporary look at Hermann Hankel's 1861 pioneering work on Lagrangian fluid dynamics" by Frisch, Grimberg and Villone (2017). Here we present the English translation of the 1861 prize manuscript from G\"ottingen University "Zur allgemeinen Theorie der Bewegung der Fl\"ussigkeiten" (On the general theory of the motion of the fluids) of Hermann Hankel (1839-1873), which was originally submitted in Latin and then translated into German by the Author for publication. We also provide the English translation of two important reports on the manuscript, one written by Bernhard Riemann and the other by Wilhelm Eduard Weber, during the assessment process for the prize. Finally we give a short biography of Hermann Hankel with his complete bibliography., Comment: 44 pages, see the companion paper by Frisch, Grimberg and Villone (2017), v2: minor revisions including change of title, accepted for publication in EPJH

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$., v1: 12 pages; v2: 14 pages, title changed, added references and a new section about theoretical limitations. v2 matches published version in PRD

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., Comment: v2: 5 pages, 4 figures. Published at MNRAS Letters

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, 20 pages, 8 figures

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., Comment: 9 pages; received 2017 May 24, and accepted 2017 June 21 at MNRAS

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., 10 pages + appendix, 4 figures - v2: typos corrected; corresponds to published version

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}$., Comment: 25 pages, 11 figures; accepted for publication in JCAP

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., Comment: 17 pages, v3: minor revisions together with some new equations, matches published version

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., v2: 15 pages, 6 figures, revised version to match published version; v3: minor fixes regarding scaling of the relative energy density; v4: references added

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., Comment: 22 pages, 6 figures, accepted JCAP version

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., V2: Second order density added and other expressions simplified. References updated. 15 pages

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., Comment: 16 pages, 4 figures, published in MNRAS, this paper introduces a convergent formulation of Lagrangian perturbation theory for LCDM

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., 5 pages, 3 figures, revised text and figures, matches published version in PRD

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., 44 pages, 14 figures; included scalar Boltzmann hierarchy in the massless case and plots of integral kernels; accepted by JCAP

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., 7 pages, 2 figures, invited contribution to the Classical and Quantum Gravity focus issue on "Relativistic Effects in Cosmology", edited by Kazuya Koyama; v2 small modifications, matches published version

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., 11 pages, 3 figures, version published in PRD

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., Comment: 19 pages, two figures, improved discussion, matches published version

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., Comment: Title adjusted, added a table for clarity, matches published version

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., Comment: 14 pages, results unchanged, extended discussion, matches published version

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., Comment: 5 pages, corrected a typo, accepted for publication in MNRAS Letters

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., 12 pages, clarified the notation in the appendix

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., 33 pages, 14 of them appendix, 4 figures, matches published version in JCAP

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., Comment: 18 pages, 4 figures; v2: references added, matches published version

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., Comment: 31 pages, matches published version in JCAP, added an extra section which discusses and motivates the choice of initial conditions, extended the title for the sake of precision

We present a novel approach to generate higher-order initial conditions (ICs) for cosmological simulations that take into account the distinct evolution of baryons and dark matter. We focus on the numerical implementation and the validation of its performance, based on both collisionless N-body simulations and full hydrodynamic Eulerian and Lagrangian simulations. We improve in various ways over previous approaches that were limited to first-order Lagrangian perturbation theory (LPT). Specifically, we (1) generalize nth-order LPT to multi-fluid systems, allowing 2LPT or 3LPT ICs for two-fluid simulations, (2) employ a novel propagator perturbation theory to set up ICs for Eulerian codes that are fully consistent with 1LPT or 2LPT, (3) demonstrate that our ICs resolve previous problems of two-fluid simulations by using variations in particle masses that eliminate spurious deviations from expected perturbative results, (4) show that the improvements achieved by going to higher-order PT are comparable to those seen for single-fluid ICs, and (5) demonstrate the excellent (i.e., few per cent level) agreement between Eulerian and Lagrangian simulations, once high-quality initial conditions are used. The rigorous development of the underlying perturbation theory is presented in a companion paper. All presented algorithms are implemented in the Monofonic Music-2 package that we make publicly available., Comment: 20 pages, 14 figures, MNRAS in press, revised version discusses how relative velocity mode can be included to leading order, code available from https://bitbucket.org/ohahn/monofonic