A new truncation scheme based on the cumulant expansion of the one-particle phase-space distribution function for dark matter particles is developed. Extending the method of moments in relativistic kinetic theory, we derive evolution equations which supplement the covariant conservation of the energy-momentum tensor and particle number current. Truncating the cumulant expansion we obtain a closed, covariant and hyperbolic system of equations which can be used to model the evolution of a general relativistic non-ideal fluid. As a working example we consider a Friedmann-Lema\^itre-Robertson-Walker cosmology with dynamic pressure and solve for the time evolution of the effective equation of state parameter.

Dark matter evolution during the process of cosmological structure formation can be described in terms of a one-particle irreducible effective action at a characteristic scale $k_m$ and a loop expansion below this scale, based on the effective propagators and vertices. We calculate the form of the effective vertices and compute the bispectrum of density perturbations within a one-loop approximation. We find that the effective vertices play a subdominant role as compared to the effective viscosity and sound velocity that modify the (inverse) propagators. For the bispectrum we reproduce the results of standard perturbation theory in the range where it is applicable, and find a slightly improved agreement with N-body simulations at larger wavenumbers.

Cosmological perturbation theory for the late Universe dominated by dark matter is extended beyond the perfect fluid approximation by taking the dark matter velocity dispersion tensor as an additional field into account. A proper tensor decomposition of the latter leads to two additional scalar fields, as well as a vector and a tensor field. Most importantly, the trace of the velocity dispersion tensor can have a spatially homogeneous and isotropic expectation value. While it decays at early times, we show that a back-reaction effect quadratic in perturbations makes it grow strongly at late times. We compare sterile neutrinos as a candidate for comparatively warm dark matter to weakly interacting massive particles as a rather cold dark matter candidate and show that the late time growth of velocity dispersion is stronger for the latter. Another feature of a non-vanishing velocity dispersion expectation value is that it destroys the apparent self-consistency of the single-stream approximation and allows thereby to treat times and scales beyond shell-crossing.

Treating dark matter at large scales as an effectively viscous fluid provides an improved framework for the calculation of the density and velocity power spectra compared to the standard assumption of an ideal pressureless fluid. We discuss how this framework can be made concrete through an appropriate coarse-graining procedure. We also review results that demonstrate that it improves the convergence of cosmological perturbation theory.

Following an approach of Matarrese and Pietroni, we derive the functional renormalization group (RG) flow of the effective action of cosmological large-scale structures. Perturbative solutions of this RG flow equation are shown to be consistent with standard cosmological perturbation theory. Non-perturbative approximate solutions can be obtained by truncating the a priori infinite set of possible effective actions to a finite subspace. Using for the truncated effective action a form dictated by dissipative fluid dynamics, we derive RG flow equations for the scale dependence of the effective viscosity and sound velocity of non-interacting dark matter, and we solve them numerically. Physically, the effective viscosity and sound velocity account for the interactions of long-wavelength fluctuations with the spectrum of smaller-scale perturbations. We find that the RG flow exhibits an attractor behaviour in the IR that significantly reduces the dependence of the effective viscosity and sound velocity on the input values at the UV scale. This allows for a self-contained computation of matter and velocity power spectra for which the sensitivity to UV modes is under control.

Cosmological perturbations of sufficiently long wavelength admit a fluid dynamic description. We consider modes with wavevectors below a scale $k_m$ for which the dynamics is only mildly non-linear. The leading effect of modes above that scale can be accounted for by effective non-equilibrium viscosity and pressure terms. For mildly non-linear scales, these mainly arise from momentum transport within the ideal and cold but inhomogeneous fluid, while momentum transport due to more microscopic degrees of freedom is suppressed. As a consequence, concrete expressions with no free parameters, except the matching scale $k_m$, can be derived from matching evolution equations to standard cosmological perturbation theory. Two-loop calculations of the matter power spectrum in the viscous theory lead to excellent agreement with $N$-body simulations up to scales $k=0.2 \, h/$Mpc. The convergence properties in the ultraviolet are better than for standard perturbation theory and the results are robust with respect to variations of the matching scale.

The dissipation of energy from local velocity perturbations in the cosmological fluid affects the time evolution of spatially averaged fluid dynamic fields and the cosmological solution of Einstein's field equations. We show how this backreaction effect depends on shear and bulk viscosity and other material properties of the dark sector, as well as the spectrum of perturbations. If sufficiently large, this effect could account for the acceleration of the cosmological expansion.