Codex @codex  0

For the complex scalar field, the canonical momenta are $\pi ={\dot{\varphi }}^{†}$ and ${\pi }^{†}=\dot{\varphi }$. Split $H=H_0+H_I$. The interaction picture is obtained from the Schrödinger picture by $|{\Psi }_I(t)⟩=e^{i{H}_{0}t}|{\Psi }_S(t)⟩$ and $O_I(t)=e^{i{H}_{0}t}{O}_S{e}^{-i{H}_{0}t}$, choosing a common time origin. The fields therefore obey the free Klein-Gordon equation, while the states evolve through $H_I(t)$. The canonical quantization of a complex scalar field givesThe adjoint expansion contains $a^{†}$ and $b$. Particle and antiparticle modes are independent; their nonzero commutators areAll cross-species commutators, and all pairs of creation operators or of annihilation operators, vanish. These formulas reproduce $[\varphi (\mathbf{x}),\pi (\mathbf{y})]=[{\varphi }^{†}(\mathbf{x}),{\pi }^{†}(\mathbf{y})]=i{\delta }^{(3)}(\mathbf{x}-\mathbf{y})$ and the remaining vanishing equal-time commutators.

There are no derivative interactions, so the interaction Hamiltonian is $H_I(t)=-\int d^{3}x{\mathcal{L}}_I(x)$. Differentiating the interaction picture state gives $i{\partial }_t|{\Psi }_I⟩=H_I(t)|{\Psi }_I⟩$. Its time-evolution operator satisfiesIteration gives ordered integration regions $t_1>\cdots >t_n$. Each is $1/n!$ of the full region with a time-ordered product, yielding the Dyson seriesWith the usual asymptotic scattering prescription, taking $t_0\to -\infty$ and $t\to +\infty$ provesTime ordering is essential because different interaction Hamiltonians need not commute.

At first order, use the connected tree scattering amplitude, equivalently the normal-ordered quartic vertex with the vacuum and one-particle mass contributions removed. This qualification matters if the local product is interpreted literally without normal ordering: single-vertex tadpole diagrams and vacuum diagrams are not the requested connected two-particle amplitude. The complex scalar quartic contact vertex is obtained fromIn a fully connected matrix element, an incoming particle and outgoing antiparticle attach to the two $\varphi$ factors; an outgoing particle and incoming antiparticle attach to the two ${\varphi }^{†}$ factors. Each pair has two assignments, giving $2!2!=4$. The external relativistic normalization of a one-particle state cancels each mode measure: for example $[\varphi (x),a^{†}(p_1)]=e^{-i{p}_{1}x}$ and $[b(q_2),\varphi (x)]=e^{i{q}_{2}x}$. Consequently the connected integrand is $4e^{-i(p_1+q_1-p_2-q_2)\cdot x}$, andThus the invariant connected amplitude isThe Dirac delta function expresses four-momentum conservation.

For the identity contribution, commute $a(p_2)$ through $a^{†}(p_1)$ and independently $b(q_2)$ through $b^{†}(q_1)$; all terms leaving a creation operator against the left Fock vacuum or an annihilation operator against the right Fock vacuum vanish. HenceThere is no exchange term between the two distinct particle and antiparticle species.

All three fields have mass $m$. In the interaction picture the quantum fields evolve with the free Hamiltonian operator, so their equations areThe interactions evolve the states instead. For $\chi \in \{\varphi ,\psi ,\Phi \}$ the appropriate mode expansion of a free field isIndependent species have commuting creation and annihilation operators. The full set of canonical commutation relations isThese follow from the separate equal-time canonical commutation relations for the three real scalar fields. For comparison, in the Heisenberg picture the full Euler-Lagrange field equations contain interactions: $(\square +m^2){\varphi }_H=-g{\varphi }_H{\Phi }_H$, $(\square +m^2){\psi }_H=-g{\psi }_H{\Phi }_H$, and $(\square +m^2){\Phi }_H=-\frac{g}{2}({\varphi }_H^2+{\psi }_H^2)$. These are distinct from the free equations used in the interaction picture.

There are no derivative couplings, so the interaction Hamiltonian is $H_I(t)=-\int d^{3}x{\mathcal{L}}_I(x)$. The interaction picture evolution operator obeysIterating its integral equation gives the Dyson seriesThe ordered integration region is equivalent to the unrestricted region with time ordering and the factor $1/n!$. Taking the scattering limit $t_0\to -\infty$, $t\to \infty$, with the usual adiabatic prescription for perturbation theory, therefore yields the S-matrixAll fields inside this expression are interaction picture fields.

Use relativistic normalization of a one-particle state and take $|i⟩=a_{\varphi }^{†}(\mathbf{p})a_{\psi }^{†}(\mathbf{q})|0⟩$, $|f⟩=a_{\varphi }^{†}({\mathbf{p}}^{\prime })a_{\psi }^{†}({\mathbf{q}}^{\prime })|0⟩$. Define the connected scattering amplitude by $⟨f|S|i{⟩}_{\mathrm{conn}}=i(2\pi {)}^4{\delta }^4(p^{\prime }+q^{\prime }-p-q)\mathcal{M}$. The identity contribution describes no scattering; the term linear in $g$ cannot connect these four external particles. At quadratic order, the connected term in the Dyson series isBy Wick theorem, a connected Wick contraction must attach the two external $\varphi$ legs to ${\varphi }_I^2$ at one vertex, the two external $\psi$ legs to ${\psi }_I^2$ at the other, and contract the remaining two ${\Phi }_I$ fields. The two assignments of species to $x,y$ cancel the Dyson factor $1/2!$. At either vertex, the $2!$ assignments of its identical fields to the incoming and outgoing particle cancel the interaction factor $1/2$. Equivalently, the coefficient $-g^2/8$ is multiplied by $2\times 2\times 2$, giving $-g^2$. Therefore, writing $D_F(x-y)=⟨0|T{\Phi }_I(x){\Phi }_I(y)|0⟩$,Insert the Feynman propagator $D_F(x-y)=\int \frac{d^{4}r}{(2\pi {)}^4}\frac{i{e}^{-ir\cdot (x-y)}}{r^2-m^2+i0}$. The two vertex integrals impose $r=p^{\prime }-p$ and $p^{\prime }+q^{\prime }=p+q$. With the Mandelstam invariant $t=(p-p^{\prime }{)}^2$, this givesThis is distinguishable scalar scattering by a shared mediator with mediator mass $m$. In physical elastic kinematics $t\le 0$, so the denominator has no pole for $m>0$. There is only the displayed $t$-channel tree-level Feynman diagram: an $s$- or $u$-channel would require an absent mixed $\varphi \psi \Phi$ vertex. Each vertex contributes $-ig$ and the exchanged $\Phi$ contributes $i/(t-m^2+i0)$, agreeing with the direct Wick contraction calculation.

The interaction contains no time derivatives, so $\pi =\dot{\varphi }$. The Legendre transform in mechanics givesThusHere $H_0=\int {\mathcal{H}}_0$ is the free Hamiltonian operator of mass $m$, and $H_I=\int {\mathcal{H}}_I$ is the interaction Hamiltonian.

In the interaction picture, operators evolve with $H_0$ and the states with $H_I(t)=e^{i{H}_{0}t}{H}_I{e}^{-i{H}_{0}t}$. With $U_I(t_0,t_0)=I$, the state evolution satisfies $i{\partial }_t{U}_I(t,t_0)=H_I(t)U_I(t,t_0)$. Integrating once and iterating gives the Dyson seriesEquivalently it is the time ordering of the exponential. Taking the infinite-time scattering limit, with adiabatic switching understood, gives the scattering matrixThe field in the interaction is the free interaction-picture field, which is essential for applying the free-field contractions.

Use relativistic normalization of a one-particle state, $|\mathbf{p}⟩=\sqrt{2E_{\mathbf{p}}}{a}^{†}(\mathbf{p})|0⟩$, so each external scalar contraction contributes its plane-wave phase without an extra $1/\sqrt{2E}$ factor. For the connected on-shell two-to-two scattering amplitude, the first Dyson series term is $-i\lambda \int d^{4}x{\varphi }_I(x{)}^4/4!$. The Wick theorem gives $4!$ ways to attach the two incoming and two outgoing external particles to the four fields, canceling the denominator. ConsequentlywhereThe delta function expresses four-momentum conservation. The connected scattering convention removes spectator, vacuum and external self-energy contributions; equivalently at this order one can normal-order the interaction and keep the physical external mass fixed. These contributions should not be mistaken for additional connected four-point diagrams.

Use ${\eta }_{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =c=1$. In the interaction picture, fields evolve with the free Hamiltonian, while interactions enter the state evolution through the Dyson series. Write the Klein-Gordon field as $\varphi ={\varphi }^{(+)}+{\varphi }^{(-)}$, with the positive-frequency part of a quantum field containing annihilation operators and ${\varphi }^{(-)}=({\varphi }^{(+)}{)}^{†}$ containing creation operators. The mode normalization in the PDF corresponds to $[a_{\mathbf{p}},a_{\mathbf{q}}^{†}]=(2\pi {)}^3{\delta }^3(\mathbf{p}-\mathbf{q})$ and $E_{\mathbf{p}}=\sqrt{|\mathbf{p}{|}^2+m^2}$.

The time-ordered product places fields with the latest time on the left. For two bosonic insertions,The normal-ordered product $:\varphi (x_1)\cdots \varphi (x_n):$ instead moves every creation operator to the left of every annihilation operator, without adding the commutators generated by that rearrangement. Its vacuum expectation value is zero whenever a nonempty set of free-field factors remains.

The mode commutator gives the positive-frequency Wightman functionDefine the Feynman propagator here to be the two-point function itself, rather than that function divided by $i$. ThenEquivalently,The positive-energy pole is below and the negative-energy pole above the real axis. Closing the energy contour integral below for positive time separation and above for negative time separation gives the two terms above. Thus this is the scalar Feynman propagator pole prescription, and $(\square +m^2){\Delta }_F=-i{\delta }^4$. In particular, $T\{\varphi (x)\varphi (y)\}=:\varphi (x)\varphi (y):+{\Delta }_F(x-y)$: the difference between time ordering and normal ordering is the Wick contraction.

For $n$ free bosonic fields, the Wick theorem stateswhere $P$ runs over all sets of disjoint pairs of indices, including the empty set, and $V(P)$ is the set of paired indices. The empty normal-ordered product equals one. This is an operator identity for free fields, understood with smearing or a regulator when products of operator-valued distributions require it.

Here is an inductive proof of the scalar Wick theorem. Sort the arguments chronologically so that $x_1^0$ is the latest. For any remaining set $I$, commuting the annihilation part of ${\varphi }_1$ through the creation parts in a normal-ordered product givesThe creation part of ${\varphi }_1$ is already in the correct position. Every commutator on the right is a scalar, so no further operator rearrangement is needed for it. Because $x_1$ is latest, that commutator is ${\Delta }_F(x_1-x_j)$. Assume the Wick theorem for the other $n-1$ fields and multiply its expansion by ${\varphi }_1$. The first term of the displayed identity gives every pairing in which index 1 is unpaired; the sum gives every pairing in which 1 is paired with exactly one $j$. Each final partial pairing has exactly one such construction. The one-field identity supplies the base case, proving the theorem by mathematical induction. Applying time ordering restores arbitrary argument order.

Taking the vacuum expectation value retains only complete pairings. Odd free correlators vanish; for example the free scalar four-point function is ${\Delta }_{12}{\Delta }_{34}+{\Delta }_{13}{\Delta }_{24}+{\Delta }_{14}{\Delta }_{23}$. These facts explain why propagators are the edges of Feynman diagrams.

For the interacting quartic scalar field theory, the interaction Hamiltonian density is $\lambda {\varphi }^4/4!$. With a single time ordering acting on external insertions and interaction vertices, the unnormalized numerator expands asApply the proved Wick theorem term by term. Each complete pairing connects external insertions and four-legged vertices by Feynman propagators. There are $4!$ assignments of four distinct incident legs at a vertex, canceling its factorial. Permuting identical integrated vertices cancels the corresponding part of $v!$. Permutations still leaving the labeled graph unchanged give its Feynman-diagram symmetry factor. For example,whereas the two-point tadpole diagram has coefficient $-i\lambda /2$, because choosing the two vertex legs contracted with the labeled external fields gives $4\cdot 3$ possibilities and the other two contract with each other.

Thus the position-space rules attach ${\Delta }_F$ to every edge, $-i\lambda \int d^{4}z$ to every four-legged interaction vertex, and divide by the graph's Feynman-diagram symmetry factor. Fourier transformation gives the momentum-space Feynman rules:Enforce momentum conservation at each vertex, retaining the overall four-momentum delta function. External lines of these correlation functions remain unamputated. The unnormalized expression includes Vacuum Feynman diagrams; dividing by $N_0=⟨0|S|0⟩$ cancels those disconnected vacuum factors. They should not be discarded without that normalization.

There is a time-ordering convention in the printed correlation-function notation worth making explicit. The usual interacting correlator uses the single global time ordering in $N_n$ above. If the displayed $T\{{\varphi }_1\cdots {\varphi }_n\}S$ is read literally with $S$ outside that ordering, the external block and the vertex block are separately time ordered: their cross-contractions are $D^{+}(x_i-z)$ because each external field stands to the left of every vertex field. The same pairing construction then uses ${\Delta }_F$ for lines within either block and $D^{+}$ for lines between blocks; in momentum space the latter is $2\pi \theta (p^0)\delta (p^2-m^2)$. This literal product is different from the conventional globally time-ordered interacting correlator, and does not have the ordinary external Feynman propagators derived above.

Let $X={\partial }^{-2}\zeta$, with an inverse spatial Laplacian defined by a boundary or zero-mode prescription. It depends on $\zeta$, not on its time derivative. Differentiation gives the canonical momentumCount $\zeta ,\pi ,X$ as first order. Inverting perturbatively givesFor the Legendre transform in mechanics, write $v=\dot{\zeta }=\pi +\Delta$ with $\Delta$ second order. The quadratic kinetic contribution is $\pi v-v^2/2={\pi }^2/2-{\Delta }^2/2$, whose correction is fourth order; in the cubic terms one may consequently replace $v$ by $\pi$. The Hamiltonian density is thereforeThe cubic interaction Hamiltonian density isIn the free interaction picture, ${\pi }_I={\dot{\zeta }}_I$. Thus here ${\mathcal{H}}_{\mathrm{int}}=-{\mathcal{L}}_3$ with the time derivatives interpreted as free fields. This equality follows from the cubic Legendre transform for scalar derivative interactions; at higher orders additional terms can appear.

Use the paper's mode normalization, effectively $M_p=1$, and set $C=ϵ+1-c_s^2$, $K=k_1+k_2+k_3$, $s_2=k_1{k}_2+k_1{k}_3+k_2{k}_3$, $s_3=k_1{k}_2{k}_3$. Write $⟨{\zeta }_1{\zeta }_2{\zeta }_3{⟩}_c=(2\pi {)}^3{\delta }^3(\sum {\mathbf{k}}_i)B$. There are two source consistency issues: the supplied interaction picture expansion makes an external-before-vertex Wick contraction contain $u_k({\tau }^{\prime })$, with positive exponential $e^{ik{c}_s{\tau }^{\prime }}$; it converges at $-\infty (1-i0)$, whereas the PDF prints $-\infty (1+i\epsilon )$. Also the stated positive interaction has the opposite overall bispectrum sign from the printed target. The following calculation identifies both issues directly.

For the stated mode expansion, the ordered contraction isA spatial derivative at the vertex supplies $i{p}_j$; the two differentiated legs consequently give $-{\mathbf{k}}_i\cdot {\mathbf{k}}_j$. There are two Wick contractions for each choice of the undifferentiated external leg. The extra $a({\tau }^{\prime })$ in the stipulated in-in formalism integral makes the vertex coefficient $a^2ϵC/c_s^2$, using $a^2=1/(H^2{\tau }^{\prime 2})$. This is consistent if the named $H_{\mathrm{int}}$ is the cosmic-time generator expressed as a function of $\tau$, since $dt=ad\tau$; if it were already the conformal-time generator, that extra $a$ would have to be removed.

For the spatial-gradient cubic curvature bispectrum, define the late-time scalar gradient vertex integralWith $\omega =c_{s}K$, the first two terms of the product obey $({\tau }^{\prime -2}-i\omega {\tau }^{\prime -1})e^{i\omega {\tau }^{\prime }}=-{\partial }_{{\tau }^{\prime }}(e^{i\omega {\tau }^{\prime }}/{\tau }^{\prime })$. The other two use ${\int }_{-\infty }^0{e}^{i\omega {\tau }^{\prime }}d{\tau }^{\prime }=-i/\omega$ and ${\int }_{-\infty }^0{\tau }^{\prime }{e}^{i\omega {\tau }^{\prime }}d{\tau }^{\prime }=1/{\omega }^2$, with the same vacuum prescription. The endpoint divergence $-1/{\tau }_f$ is real, whileIn particular the finite $-K$ term comes from the endpoint of the total derivative and must not be discarded along with the real divergence.

Let $S={\mathbf{k}}_1\cdot {\mathbf{k}}_2+{\mathbf{k}}_1\cdot {\mathbf{k}}_3+{\mathbf{k}}_2\cdot {\mathbf{k}}_3$. The ordered vertex expectation, after stripping the momentum delta, is $Z=-2S{H}^{4}CI/(64{ϵ}^2{c}_s^5{k}_1^3{k}_2^3{k}_3^3)$. Since $\mathrm{Re}(-2iZ)=2\mathrm{Im}Z$, the literal positive Hamiltonian givesThe printed positive-sign target is obtained by replacing the displayed interaction Hamiltonian by its negative, with the convergent contour above. This is also the expected sign if a positive cubic $\zeta (\partial \zeta {)}^2$ coefficient was intended as a Lagrangian term, by the interaction Hamiltonian relation proved in question 3(a). Under that explicit repair,Restoring the overall $(2\pi {)}^3{\delta }^3(\sum {\mathbf{k}}_i)$ reproduces all three cyclic terms in the requested expression. An equilateral bispectrum configuration is a concrete sign check: for $k_i=k$, $S=-3k^2/2$, $Q=-17k/9$, so $SQ=17k^3/6>0$. The literal interaction gives $B=-17H^{4}C/(96{ϵ}^2{c}_s^4{k}^6)$, while the printed target is positive. This cannot be repaired by merely changing the overall definition of $\zeta$, because its field expansion and interaction were specified together. The contour typo and Hamiltonian sign are therefore documented source repairs, rather than silently altered contractions.

For a centred Gaussian random field, Wick theorem expresses each even correlation function as a sum over all pairings of two-point functions; odd moments vanish. For example,For six fields there are $5\cdot 3\cdot 1=15$ pairings. In a tree-level connected three-point calculation with one cubic vertex, each of the three external fields must contract with a different vertex field. This gives $3!=6$ connected Wick contractions. Pairings connecting two external fields instead are disconnected tadpole contributions, omitted from the primordial bispectrum.

Put $\mathcal{F}=(1-c_s^2)\mathcal{A}/c_s^2$. Using $dt=ad\tau$ and $\dot{\zeta }={\zeta }^{\prime }/a$, the integrated interaction Hamiltonian becomesThe unequal-time Wick contraction required by the displayed operator order isThe conjugate is essential: the equal-time power spectrum alone does not specify the erroneous unstarred replacement printed in the target expression. The in-in bispectrum conjugation rule fixes the vacuum convergence and final sign.

Fourier transforming the vertex gives $(2\pi {)}^3{\delta }^{(3)}({\mathbf{p}}_1+{\mathbf{p}}_2+{\mathbf{p}}_3)$ with three measures $d^3{p}_j/(2\pi {)}^3$. The six pairings supply the product of three external-vertex delta functions, each carrying $(2\pi {)}^3$. Integrating all three internal momenta leaves one overall $(2\pi {)}^3{\delta }^{(3)}({\mathbf{k}}_1+{\mathbf{k}}_2+{\mathbf{k}}_3)$. Thus the normalized connected contraction formula isThis spells out what the printed permutations and delta functions must represent: all six permutations, full Fourier normalization, and conjugated vertex modes. Three cyclic permutations alone would miss a factor of two.

At cubic order, the expansion of $P(X,\varphi )$ produces the schematic operatorsThe derivative self-interactions involving $P_{XX}$ and $P_{XXX}$ are generally largest when $c_s\ll 1$, while coefficients containing explicit $\varphi$ derivatives of a slowly varying $P$ are commonly slow-variation suppressed. The former therefore tend to dominate primordial non-Gaussianity.

The terms contributing to $\phi {\phi }^{\prime 2}$ areThusEquivalently, if the complete time-dependent coefficient is called $\lambda$, then $\lambda (\tau )=a^2{\lambda }_0$.

Treat ${\lambda }_0$ as constant at leading slow variation. To cubic order the corresponding interaction Hamiltonian isFor $q_i=c_s{k}_i$, the mode function above satisfiesWriting $K=k_1+k_2+k_3$, the needed regulated integral isThere are three choices for the undifferentiated leg and two contractions interchanging the differentiated legs. The stated in-in formalism formula therefore giveswithThe powers of $c_s$ cancel for this vertex with the normalization specified in part (i). Reversing the convention for the sign of $H_I$ reverses the displayed overall sign but not the momentum shape of the primordial bispectrum.