For plane-parallel pressureless matter in an Einstein-de Sitter universe, before shell crossing, the symmetrized growing-mode standard perturbation theory density kernel is
In particular . An unsymmetrized representation of the one-dimensional beta mode-coupling kernel is ; its symmetric part is the expression above.
Expand , where the standard perturbation theory density kernel convolves linear fields. Through fourth order in the Gaussian linear field, the only power-spectrum diagrams are the tree contraction , the loop joining two vertices , and the two orderings that join an vertex to a linear leg, . Their expressions are
Figure 1.
One-loop matter-power-spectrum diagrams and schematic present-day contributions
. The three required contraction topologies are P11, P22, and P13 plus P31. At low wavenumber the one-loop sum approaches the linear spectrum; the loop terms become appreciable around 0.1 inverse megaparsecs times h, and fixed-order perturbation theory eventually fails. Details.
For symmetrized standard perturbation theory density kernels, the two one-loop contractions are
and
The factors and count the relevant Wick contractions. Since the cross-correlation occurs in both orders, the one-loop matter power spectrum is
In standard perturbation theory in cosmology, is an -leg vertex carrying the symmetrized standard perturbation theory density kernel . For Gaussian random fields, a connected four-point diagram with loops contains linear power spectra, hence linear fields.
At tree level the required perturbative-order partitions are
The first is a star with one vertex and three linear external legs. The second has two vertices joined by one internal power-spectrum line, with one linear external leg attached to each vertex. At one loop all five partitions of eight into four positive parts are required:
These are conventionally denoted , , , , and ; all inequivalent placements of external legs and all Wick contractions are understood. This list is the complete set of one-loop matter trispectrum topologies.
Let and . With , the star contribution is
where . For the exchange contribution, sum over the six unordered choices of linear external legs and let be the complementary pair:
Thus the connected tree trispectrum is .
The pressureless single-stream equations cease to be a complete one-loop prediction because loop momenta probe short nonlinear scales and produce ultraviolet-sensitive terms. A consistent result must include the effective field theory of large-scale structure: effective-stress counterterms, including their insertions into the tree topologies, and stochastic terms at the order allowed by mass conservation and momentum conservation. Their coefficients absorb short-scale dependence and must be fitted or matched. If the observable is a biased tracer rather than matter itself, the corresponding renormalized bias expansion is also required.
The ultraviolet softness of the second-order density kernel requires
for fixed as . In this limit the two hard vectors are opposite,
The constant part of the stated kernel is therefore
It must vanish, giving
This recovers the standard standard perturbation theory density kernel .
Use the Fourier convention with . Neglect anisotropic stress and assume the pressureless velocity is irrotational, so
Fourier transforming in the nonlinear continuity equation gives
where the alpha mode-coupling kernel and its symmetrization are
Taking the divergence of the Euler equation gives
with the symmetric beta mode-coupling kernel
In the Einstein-de Sitter universe, and . At linear order the ansatz gives . At second order, the continuity and Euler equations become
where each integral includes the momentum-conserving measure above. Eliminating yields
with the standard perturbation theory density kernel
For Gaussian initial conditions and the Einstein-de Sitter second-order standard perturbation theory density kernel, the unsmoothed density contrast obeys
There are three choices of the second-order field and two connected Wick contractions; the angular average of is . Thus the reduced skewness is . Smoothing changes this result, and the unsmoothed calculation requires a finite variance.