We derive the Edgeworth streaming model (ESM) for the redshift space correlation function starting from an arbitrary distribution function for biased tracers of dark matter by considering its two-point statistics and show that it reduces to the Gaussian streaming model (GSM) when neglecting non-Gaussianities. We test the accuracy of the GSM and ESM independent of perturbation theory using the Horizon Run 2 N-body halo catalog. While the monopole of the redshift space halo correlation function is well described by the GSM, higher multipoles improve upon including the leading order non-Gaussian correction in the ESM: the GSM quadrupole breaks down on scales below 30 Mpc/h whereas the ESM stays accurate to 2% within statistical errors down to 10 Mpc/h. To predict the scale dependent functions entering the streaming model we employ Convolution Lagrangian perturbation theory (CLPT) based on the dust model and local Lagrangian bias. Since dark matter halos carry an intrinsic length scale given by their Lagrangian radius, we extend CLPT to the coarse-grained dust model and consider two different smoothing approaches operating in Eulerian and Lagrangian space, respectively. The coarse-graining in Eulerian space features modified fluid dynamics different from dust while the coarse-graining in Lagrangian space is performed in the initial conditions with subsequent single streaming dust dynamics, implemented by smoothing the initial power spectrum in the spirit of the truncated Zel'dovich approximation. Finally, we compare the predictions of the different coarse-grained models for the streaming model ingredients to N-body measurements and comment on the proper choice of both the tracer distribution function and the smoothing scale. Since the perturbative methods we considered are not yet accurate enough on small scales, the GSM is sufficient when applied to perturbation theory.

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

We investigate the relation between the standard Newtonian equations for a pressureless fluid (dust) and the Einstein equations in a double expansion in small scales and small metric perturbations. We find that parts of the Einstein equations can be rewritten as a closed system of two coupled differential equations for the scalar and transverse vector metric perturbations in Poisson gauge. It is then shown that this system is equivalent to the Newtonian system of continuity and Euler equations. Brustein and Riotto (2011) conjectured the equivalence of these systems in the special case where vector perturbations were neglected. We show that this approach does not lead to the Euler equation but to a physically different one with large deviations already in the 1-loop power spectrum. We show that it is also possible to consistently set to zero the vector perturbations which strongly constrains the allowed initial conditions, in particular excluding Gaussian ones such that inclusion of vector perturbations is inevitable in the cosmological context. In addition we derive nonlinear equations for the gravitational slip and tensor perturbations, thereby extending Newtonian gravity of a dust fluid to account for nonlinear light propagation effects and dust-induced gravitational waves.

Contrary to recent claims in the literature, Newtonian N-body simulations of collisionless Dark Matter in a LambdaCDM background are compatible with general relativity and are not invalidated by general relativistic effects at the linear level. This verdict is based on four facts. (1) The system of linearized Einstein equations and conservation laws is well-posed in the gauge invariant formulation and physically meaningful. (2) Comparing general relativity with its Newtonian approximation at a given order in perturbation theory is only meaningful at the level of observables. (3) The dynamics of observables describing a dust fluid in general relativity and its Newtonian approximation agree at the linear level. Any disagreement for observables on the lightcone are well-known, of which the most dominant is gravitational lensing. (4) Large fluctuations in the Hubble parameter contribute significantly only to gravitational lensing effects. Therefore, these fluctuations are not in conflict with Newtonian N-body simulations beyond what has already been carefully taken into account using ray tracing technology.