Solution (source code)

= Solution

Expand $\delta=a\delta_1+a^2\delta_2+a^3\delta_3+\cdots$. For a Gaussian linear density field, odd linear correlators vanish and <Wick theorem> reduces the fourth-order terms to products of the linear <power spectrum>. Hence
$$
P_{\rm 1-loop}(k)
=a^2P(k)+a^4[P_{22}(k)+2P_{13}(k)].
$$
The two contractions of the two quadratic fields give
$$
\boxed{P_{22}(k)=2\int\frac{d^3q}{(2\pi)^3}
F_2(\mathbf q,\mathbf k-\mathbf q)^2
P(q)P(|\mathbf k-\mathbf q|)}.
$$
The three choices for which argument of $F_3$ carries the external momentum give
$$
\boxed{P_{13}(k)=3P(k)\int\frac{d^3q}{(2\pi)^3}
F_3(\mathbf k,\mathbf q,-\mathbf q)P(q)}.
$$
Thus the contribution conventionally called $P_{13}+P_{31}$ is $2P_{13}=6P(k)\int F_3P$. These are the two one-loop diagrams of the <one-loop matter power spectrum>.

Solved by gpt-5.6-sol high.