Kinetic Field Theory (KFT) is a statistical field theory for an ensemble of point-like classical particles in or out of equilibrium. We review its application to cosmological structure formation. Beginning with the construction of the generating functional of the theory, we describe in detail how the theory needs to be adapted to reflect the expanding spatial background and the homogeneous and isotropic, correlated initial conditions for cosmic structures. Based on the generating functional, we develop three main approaches to non-linear, late-time cosmic structures, which rest either on the Taylor expansion of an interaction operator, suitable averaging procedures for the interaction term, or a resummation of perturbation terms. We show how an analytic, parameter-free equation for the non-linear cosmic power spectrum can be derived. We explain how the theory can be used to derive the density profile of gravitationally bound structures and use it to derive power spectra of cosmic velocity densities. We further clarify how KFT relates to the BBGKY hierarchy. We then proceed to apply kinetic field theory to fluids, introduce a reformulation of KFT in terms of macroscopic quantities which leads to a resummation scheme, and use this to describe mixtures of gas and dark matter. We discuss how KFT can be applied to study cosmic structure formation with modified theories of gravity. As an example for an application to a non-cosmological particle ensemble, we show results on the spatial correlation function of cold Rydberg atoms derived from KFT.

In earlier work, we have developed a nonequilibrium statistical field theory description of cosmic structure formation, dubbed Kinetic Field Theory (KFT), which is based on the Hamiltonian phase-space dynamics of classical particles and thus remains valid beyond shell-crossing. Here, we present an exact reformulation of the KFT framework that allows to resum an infinite subset of terms appearing in the original perturbative expansion of KFT. We develop the general formalism of this resummed KFT, including a diagrammatic language for the resummed perturbation theory, and compute the lowest-order results for the power spectra of the dark matter density contrast and momentum density. This allows us to derive analytically how the linear growth of the largest structures emerges from Newtonian particle dynamics alone, which, to our knowledge, is the first time this has been achieved.

We derive a non-perturbative, closed, analytic equation for the non-linear power spectrum of cosmic density fluctuations. This result is based on our kinetic field theory (KFT) of cosmic structure formation and on evaluating particle interactions in the Born approximation. At redshift zero, relative deviations between our analytic result and typical numerical results are ~ 15 % on average up to wave numbers of k <= 10 h/Mpc. The theory underlying our analytic equation is fully specified once the statistical properties of the initial density and momentum fluctuations are set. It has no further adjustable parameters. Apart from this equation, our main result is that the characteristic non-linear deformations of the power spectrum at late cosmic times are determined by the initial momentum correlations and a partial compensation between diffusion damping and particle interactions.

In recent work we developed a description of cosmic large scale structure formation in terms of non-equilibrium ensembles of classical particles, with time evolution obtained in the framework of a statistical field theory. In these works, the initial Gaussian correlations between particles have so far been treated perturbatively or restricted to pure momentum correlations. Here we treat the correlations between all phase-space coordinates exactly by adopting a diagrammatic language for the different forms of correlations, directly inspired by the Mayer cluster expansion. We will demonstrate that explicit expressions for phase-space density cumulants of arbitrary $n$-point order, which fully capture the non-linear coupling of free streaming kinematics due to initial correlations, can be obtained from a simple set of Feynman rules. These cumulants will be the foundation for further investigations of interacting perturbation theory.

In earlier work, we have developed a Kinetic Field Theory (KFT) for cosmological structure formation and showed that the non-linear density-fluctuation power spectrum known from numerical simulations can be reproduced quite well even if particle interactions are taken into account to first order only. Besides approximating gravitational interactions, we had to truncate the initial correlation hierarchy of particle momenta at the second order. Here, we substantially simplify KFT. We show that its central object, the free generating functional, can be factorized, taking the full hierarchy of momentum correlations into account. The factors appearing in the generating functional, which we identify as non-linearly evolved density-fluctuation power spectra, have a universal form and can thus be tabulated for fast access in perturbation schemes. In this paper, we focus on a complete evaluation of the free generating functional of KFT, not including particle interactions yet. This implies that the non-linearly evolved power spectra contain a damping term which reflects that structures are being wiped out at late times by free streaming. Once particle interactions will be taken into account, they will compensate this damping. If we suppress this damping in a way suggested by the fluctuation-dissipation relations of KFT, our results show that the complete hierarchy of initial momentum correlations is responsible for a large part of the characteristic non-linear deformation and the mode transport in the density-fluctuation power spectrum. Without any adjustable parameters, KFT accurately reproduces the scale at which non-linear evolution sets in. Finally, we further develop perturbation theory based on the factorization of the generating functional and propose a diagrammatic scheme for the perturbation terms.

Building upon the recent pioneering work by Mazenko and by Das and Mazenko, we develop a microscopic, non-equilibrium, statistical field theory for initially correlated canonical ensembles of classical microscopic particles obeying Hamiltonian dynamics. Our primary target is cosmic structure formation, where initial Gaussian correlations in phase space are believed to be set by inflation. We give an exact expression for the generating functional of this theory and work out suitable approximations. We specify the initial correlations by a power spectrum and derive general expressions for the correlators of the density and the response field. We derive simple closed expressions for the lowest-order contributions to the non-linear cosmological power spectrum, valid for arbitrary wave numbers. We further calculate the bispectrum expected in this theory within these approximations and the power spectrum of cosmic density fluctuations to first order in the gravitational interaction, using a recent improvement of the Zel'dovich approximation. We show that, with a modification motivated by the adhesion approximation, the non-linear growth of the density power spectrum found in numerical simulations of cosmic structure evolution is reproduced well to redshift zero and for arbitrary wave numbers even within first-order perturbation theory. Our results present the first fully analytic calculation of the non-linear power spectrum of cosmic structures.

Recently Mazenko and Das and Mazenko introduced a non-equilibrium field theoretical approach to describe the statistical properties of a classical particle ensemble starting from the microscopic equations of motion of each individual particle. We use this theory to investigate the transition from those microscopic degrees of freedom to the evolution equations of the macroscopic observables of the ensemble. For the free theory, we recover the continuity and Jeans equations of a collisionless gas. For a theory containing two-particle interactions in a canonical perturbation series, we find the macroscopic evolution equations to be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy (BBGKY hierarchy) with a truncation criterion depending on the order in perturbation theory. This establishes a direct link between the classical and the field-theoretical approaches to kinetic theory that might serve as a starting point to investigate kinetic theory beyond the classical limits.

It was recently shown by Bartelmann et al. how correlated initial conditions can be introduced into the statistical field theory for classical particles pioneered by Das and Mazenko. In this paper we extend this development from the canonical to the grand canonical ensemble for a system satisfying statistical homogeneity and isotropy. We do this by translating the probability distribution for the initial phase space coordinates of the particles into an easy diagrammatic representation and then using a variant of the Mayer cluster expansion to sum over particle numbers. The grand canonical generating functional is then used in a structured approach to the derivation of the non-interacting cumulants of the two core collective fields, the density $\rho$ and the response field $B$. As a side-product we find several theorems pertaining to these cumulants which will be useful when investigating the interacting regime of the theory in future work.

We calculate the power spectrum of density fluctuations in the statistical non-equilibrium field theory for classical, microscopic degrees of freedom to first order in the interaction potential. We specialise our result to cosmology by choosing appropriate initial conditions and propagators and show that the non-linear growth of the density power spectrum found in numerical simulations of cosmic structure evolution is reproduced well to redshift zero and for arbitrary wave numbers. The main difference of our approach to ordinary cosmological perturbation theory is that we do not perturb a dynamical equation for the density contrast. Rather, we transport the initial phase-space distribution of a canonical particle ensemble forward in time and extract any collective information from it at the time needed. Since even small perturbations of particle trajectories can lead to large fluctuations in density, our approach allows to reach high density contrast already at first order in the perturbations of the particle trajectories. We argue why the expected asymptotic behaviour of the non-linear power spectrum at large wave numbers can be reproduced in our approach at any order of the perturbation series.

We use the non-equlibrium statistical field theory for classical particles, recently developed by Mazenko and Das and Mazenko, together with the free generating functional we have previously derived for point sets initially correlated in phase space, to calculate the time evolution of power spectra in the free theory, i.e. neglecting particle interactions. We provide expressions taking linear and quadratic momentum correlations into account. Up to this point, the expressions are general with respect to the free propagator of the microscopic degrees of freedom. We then specialise the propagator to that expected for particles in cosmology treated within the Zel'dovich approximation and show that, to linear order in the momentum correlations, the linear growth of the cosmological power spectrum is reproduced. Quadratic momentum correlations return a first contribution to the non-linear evolution of the power spectrum, for which we derive a simple closed expression valid for arbitrary wave numbers. This expression is a convolution of the initial density power spectrum with itself, multiplied by a mode-coupling kernel. We also derive the bispectrum expected in this theory within these approximations and show that its connected part reproduces almost, but not quite, the bispectrum expected in Eulerian perturbation theory of the density contrast.