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.

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.

We discuss opportunities that may arise from subjecting high-multiplicity events in relativistic heavy ion collisions to an analysis similar to the one used in cosmology for the study of fluctuations of the Cosmic Microwave Background (CMB). To this end, we discuss examples of how pertinent features of heavy ion collisions including global characteristics, signatures of collective flow and event-wise fluctuations are visually represented in a Mollweide projection commonly used in CMB analysis, and how they are statistically analyzed in an expansion over spherical harmonic functions. If applied to the characterization of purely azimuthal dependent phenomena such as collective flow, the expansion coefficients of spherical harmonics are seen to contain redundancies compared to the set of harmonic flow coefficients commonly used in heavy ion collisions. Our exploratory study indicates, however, that these redundancies may offer novel opportunities for a detailed characterization of those event-wise fluctuations that remain after subtraction of the dominant collective flow signatures. By construction, the proposed approach allows also for the characterization of more complex collective phenomena like higher-order flow and other sources of fluctuations, and it may be extended to the characterization of phenomena of non-collective origin such as jets.