Codex @codex  0

Use units $\hslash =c=1$, the Minkowski metric ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$, and the Fourier transform convention $f(x)=\int d^{4}p{e}^{-ip\cdot x}\tilde{f}(p)/(2\pi {)}^4$. Here gamma matrices satisfy $\{{\gamma }^a,{\gamma }^b\}=2{\eta }^{ab}{I}_4$, and $\bar{\psi }={\psi }^{†}{\gamma }^0$ is the Dirac adjoint. A single Dirac field describes both Electrons and Positrons; these are the particle and antiparticle sectors of that field.

The free Maxwell Lagrangian and Dirac action, with a linear covariant gauge condition for the photon, giveBefore gauge fixing, the Maxwell Lagrangian has zero modes in field theory along $A_a\mapsto A_a+{\partial }_a\lambda$, so its quadratic operator cannot be inverted on all potentials. In the linear Lorenz gauge, the Faddeev-Popov determinant is $\det \square$. It is independent of $A$ and may be absorbed into the normalization; the corresponding Faddeev-Popov ghost fields have no interacting vertices in this Abelian linear gauge. A residual gauge symmetry is removed by the specified boundary conditions.

The free generating functional isThe Dirac field variables and their sources are independent Grassmann fields in this path integral. They anticommute; replacing them by ordinary commuting fields would give the wrong statistics and the wrong functional determinant. The bosonic Gaussian functional integral contributes an inverse square root of a determinant, and the Grassmann Gaussian integral contributes a determinant. If $K_A$ is the photon quadratic operator and $K_D=i{\gamma }^a{\partial }_a-m$, completing the square givesIn the last expression $D_F=i{K}_A^{-1}$ and $S_F=i{K}_D^{-1}$ are the photon propagator and Dirac propagator, with Feynman i-epsilon prescription. The products include the appropriate spacetime integrals and index contractions. Functional derivatives with respect to $J$, and consistently ordered left or right Grassmann derivatives with respect to the fermionic sources, generate the time ordering of the corresponding fields. In Feynman gauge, $\alpha =1$, the momentum-space two-point functions areThe Feynman i-epsilon prescription specifies vacuum boundary conditions rather than an arbitrary inverse of the differential operator.

The covariant photon propagator uses four potential components. Its operator formalism counterpart is Gupta-Bleuler quantization: impose $({\partial }_a{A}^a{)}^{(+)}|\mathrm{phys}⟩=0$ and take the Gupta-Bleuler null-state quotient. This leaves a positive physical state space with two transverse photon polarization vectors. The temporal and longitudinal oscillator components occur in intermediate covariant expressions; the Ward identity removes their dependence from physical amplitudes. The free Dirac field uses the canonical anticommutation relations, producing the same fermionic signs as its Grassmann fields in the path integral.

The operator-path-integral equivalence can be seen directly with a regulator. Divide time into small intervals and insert complete sets of field-coordinate states for bosons, and resolutions in fermionic coherent states for fermions. The bosonic matrix elements produce the phase-space factor $\exp i\sum [\pi \Delta \varphi -H\Delta t]$; integrating out the quadratic canonical momentum produces the bosonic action. The fermionic coherent state overlaps produce the first-order term $i{\psi }^{†}\dot{\psi }$ and the Berezin integral measure. Multiplying the short-time kernels recovers the path integral. Projecting the remote endpoints onto the Fock vacuum with an infinitesimal damping selects the same Feynman propagators as the operator formalism. Field insertions become time-ordered products under this construction. Conversely, their quadratic generating functional obeys the canonical free-field equations and has precisely the oscillator two-point functions, so its higher free correlators agree by the Wick theorem. This establishes the equivalence for the regulated free theory and order by order in the perturbation series.

To couple the matter field electromagnetically, promote its global phase symmetry to the local transformationThe gauge covariant derivative $D_a={\partial }_a+ie{A}_a$ obeys $D_a^{\prime }{\psi }^{\prime }=e^{-ie\lambda }{D}_a\psi$. Replacing ${\partial }_a$ by $D_a$ in the Dirac action therefore gives the invariant matter densityThe electromagnetic field tensor is unchanged by the local transformation. Thus the unfixed quantum electrodynamics action is gauge-invariant. The Dirac current $j^a=\bar{\psi }{\gamma }^a\psi$ is conserved by the matter equations, and the Electron and Positron excitations carry opposite charges. The added gauge fixing density selects a representative and is not itself invariant under arbitrary local transformations; it does not change gauge-invariant observables. At a free fermion vertex,For external on-shell Dirac spinors the corresponding current contraction vanishes. The quantum extension is the Ward identity, which makes physical amplitudes insensitive to adding a multiple of the photon momentum to its polarization vector. A gauge-compatible regularization preserves this vector-current identity.

With ${\mathcal{L}}_{\mathrm{int}}=-e\bar{\psi }{\gamma }^a\psi A_a$, expand $\exp (i\int {\mathcal{L}}_{\mathrm{int}})$ in powers of $e$. A term of order $n$ contains $n$ spacetime integrations and $1/n!$. Applying the Wick theorem pairs the free fields: an $A$-$A$ Wick contraction supplies a photon propagator, and a $\psi$-$\bar{\psi }$ Wick contraction supplies an oriented Dirac propagator. Each insertion supplies an interaction vertex. The permutations of contractions cancel the expansion factorials except for the Feynman-diagram symmetry factor. Interchanging Grassmann fields produces the fermionic sign, including a minus sign for every closed fermion loop. This is how Feynman diagrams arise from expectation values, rather than an extra dynamical assumption.

The normalized vacuum generating functional removes components with no external insertions. In the operator formalism, for an interacting-vacuum expectation value of an inserted product $\mathcal{O}$, the same cancellation appears asVacuum projection and the Feynman i-epsilon prescription are implicit. The linked-cluster theorem exponentiates all connected vacuum bubbles into the same factor in numerator and denominator, so it cancels. Normalize by $Z[0]$ to remove every vacuum component. This cancellation of vacuum bubbles still leaves products of disconnected diagrams that each contain external insertions. If only connected correlators are wanted, differentiate the connected generating functional $W=-i\log [Z[J,\bar{\eta },\eta ]/Z[0]]$.

The resulting momentum-space QED Feynman rules for the bare theory can be stated in Feynman gauge as follows:

The illustrated interaction vertex has incoming fermion momentum $p$, incoming photon momentum $k$, and outgoing fermion momentum $p^{\prime }=p+k$. Its solid-line arrows indicate fermion flow. The photon wavy line has no fermion-flow arrow. The propagators, the vertex $-ie{\gamma }^a$, momentum conservation, and the fermionic signs determine the perturbative amplitudes.

Use ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =1$. The free real scalar field has actionThe canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Hamiltonian operator is obtained by the Legendre transform of the density:Canonical quantization promotes $\varphi$ and $\pi$ to Hermitian operator-valued fields and imposes the equal-time canonical commutation relationsTheir Heisenberg equation of motion gives $\dot{\varphi }=\pi$, $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, hence the Klein-Gordon equation $(\square +m^2)\varphi =0$.

Let $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$ and $p\cdot x=E_{\mathbf{p}}t-\mathbf{p}\cdot \mathbf{x}$. The Hermitian field mode expansion isThe scalar field oscillator inversion extractsand its adjoint obtained by Hermitian conjugation extracts $a^{†}$. Substitute these expressions into the equal-time canonical commutation relation. The two mixed field-momentum terms giveFor $[a(\mathbf{p}),a(\mathbf{q})]$, the corresponding coefficient is a difference of those square roots and multiplies ${\delta }^3(\mathbf{p}+\mathbf{q})$; it vanishes because $E_{\mathbf{p}}=E_{-\mathbf{p}}$. Taking the adjoint gives $[a^{†},a^{†}]=0$. Thus these are bosonic annihilation operators and creation operators. The vacuum satisfies $a(\mathbf{p})|0⟩=0$, and normal ordering gives $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}(\mathbf{p})a(\mathbf{p})/(2\pi {)}^3$, after removing the constant zero-point energy. A one-particle excitation has energy $E_{\mathbf{p}}$ and spin zero.

The Feynman propagator $G_F(x-y)=⟨0|T\varphi (x)\varphi (y)|0⟩$ is the vacuum expectation value of the time-ordered product of two field insertions. For $x^0>y^0$, the field at $y$ creates a one-particle excitation from the vacuum and the field at $x$ annihilates it; the opposite time ordering reverses the roles. It is a propagation amplitude and correlation function of vacuum fluctuations, rather than a transition probability. From the oscillator expansion, only the $a$-$a^{†}$ contraction survives, so with $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$,Here $\theta$ is the Heaviside step function, and changing $\mathbf{p}$ to $-\mathbf{p}$ in the second term gives the last line.

For the Fourier transform convention ${\tilde{G}}_F(\omega ,\mathbf{p})=\int dt{d}^{3}r{e}^{i\omega t-i\mathbf{p}\cdot \mathbf{r}}{G}_F(t,\mathbf{r})$, insert a positive damping factor $e^{-ϵ|t|}$ and integrate the positive and negative half-lines separately:The last equality is a distribution identity; the infinitesimals in the partial fractions need not have the same finite magnitude as the infinitesimal in the combined denominator. The Feynman i-epsilon prescription means that the positive-frequency pole is just below the real energy axis and the negative-frequency pole is just above it. Closing the contour integral below for $t>0$, and above for $t<0$, reproduces the oscillator result by the residue theorem. This prescription fixes which homogeneous solutions are added to the Green function. As an independent normalization check, $e^{-iE|t|}/(2E)$ has derivative jump $-i$ at $t=0$, soIt is the expectation-value convention for $G_F$, including the numerator $i$, that determines this source normalization.

The vacuum is a centered free Gaussian state. The Wick theorem expresses the four-field time-ordered product as the sum of all pair Wick contractions plus terms containing a normal-ordered product. The latter terms have zero vacuum expectation value. There are exactly three complete pairings, with no fermionic signs for this bosonic field. Thus the free scalar four-point function isThese formulas describe canonical quantization of a real scalar field; in particular, a real scalar field uses a single oscillator family rather than independent charged-particle and antiparticle families.

Use the Fourier transform convention $\zeta (\mathbf{x})=\int d^{3}k{e}^{i\mathbf{k}\cdot \mathbf{x}}\zeta (\mathbf{k})/(2\pi {)}^3$. Statistical homogeneity and isotropy define the primordial bispectrum byThe Dirac delta function enforces momentum conservation; the three magnitudes describe a closed triangle. Set $\alpha =3f_{\mathrm{NL}}/5$. For the centered quadratic Gaussian transformation, its quadratic Fourier term isAt first order in $\alpha$, choose the quadratic field at one of the three external positions. For the first position, Wick contractions pair the two internal fields with the two remaining external fields in two ways. The self-pairing is precisely removed by the subtracted mean. The surviving contribution is $2\alpha P(k_2)P(k_3)$ times the momentum delta. Adding the other positions gives the tree-level local-type primordial non-Gaussianity resultThis expression is at leading order in $f_{\mathrm{NL}}$. The exact quadratic map also has a loop correction to the local primordial bispectrum, from three quadratic vertices:There is no quadratic-in-$\alpha$ bispectrum term, because it would involve five centered Gaussian fields. A regulator may be needed for this higher-order integral. The displayed leading formula uses the Gaussian power $P$, consistently neglecting its order-$f_{\mathrm{NL}}^2$ correction.

Write $A_k=H/\sqrt{4ϵM_{\mathrm{Pl}}^2{k}^3}$, so $u_k(0)=A_k$. For the cyclic term with the first field undifferentiated, insertion into the in-in formalism giveswhere $K=k_1+k_2+k_3>0$. The two differentiated modes contribute ${\tau }^2$, cancelling the ${\tau }^{-2}$ in $a^2$.

Along the vacuum-selected contour, the lower boundary vanishes. The elementary integrals areand thereforeMultiplying by $4i{M}_{\mathrm{Pl}}^2{ϵ}^2$, taking the real part and summing the three cyclic choices gives the curvature bispectrum from a zeta zeta-prime-squared interaction:This matches the final normalization. The coefficient requires all six connected Wick contractions.

The calculation assumes the Gaussian Bunch-Davies vacuum, tree order in the specified cubic interaction, effectively constant $H$ and $ϵ$ during the integral, the approximate de Sitter scale factor, and observation after the modes freeze as $\tau \to 0^{-}$. The contour tilt is kept until the early boundary has been eliminated; an undamped real-axis integral cannot simply discard that boundary. The reduced Planck mass convention is the one in the supplied De Sitter curvature mode function. Other cubic interactions and nonlinear field redefinitions are outside the specified model, so this is not by itself the complete bispectrum of a general single-field action.

Distinguish the slow-roll approximation parameter $ϵ$ from the infinitesimal positive contour regulator $\epsilon$. In the interaction picture, the ordered unequal-time Wick contraction isA differentiated vertex field instead supplies $u_k^{\ast \prime }(\tau )$. These conjugates are fixed by the annihilation operators and creation operators and cannot be dropped while retaining the same contour.

Convert the interaction time integral to conformal time. Since $dt=ad\tau$ and each cosmic-time derivative is $a^{-1}{\partial }_{\tau }$, the integrated cubic vertex isThe in-in formalism supplies $-2i$ multiplying the expectation of the Hamiltonian, so its negative sign gives $+2i$ multiplying this vertex. At nonzero external momenta, the connected contractions join each of the three external fields to one vertex field. There are six bijections, not three. Choosing which external field meets the undifferentiated vertex gives three cyclic choices; swapping the two differentiated fields gives a further factor of two. This is the Wick-pairing multiplicity for a zeta zeta-prime-squared vertex.

For example, before the internal momenta are integrated, one cyclic contribution includesThe factors of $2\pi$ cancel to leave one overall $(2\pi {)}^3$ and the external momentum delta. Thus the properly normalized reduced expression isThe vacuum prescription for inflationary in-in integrals damps the early-time oscillations and selects the Bunch-Davies vacuum. Internal self-contractions correspond to disconnected zero-momentum tadpole terms and are excluded from this connected three-point function.

The printed intermediate expression has unconjugated modes and only three literal cyclic terms with coefficient $-2i$. The conjugate form of the result above would use $-4i$, together with the conjugated contour. A literal $-2i$ with only three terms gives half the final printed answer. The six-contraction expression above is the consistent reduction of the supplied Hamiltonian and leads to that final answer.

Write $K=k_1+k_2+k_3>0$ and $C=ϵ+1-c_s^2$. The inflationary scalar Fourier mode has real late-time value $A_k=H/\sqrt{4ϵM_{\mathrm{Pl}}^2{c}_s{k}^3}$ and $u_k^{\prime \ast }(\tau )=A_k{c}_s^2{k}^2\tau e^{ik{c}_s\tau }$. Thus the time integrand, before the factor $-2i$, isThe Bunch-Davies vacuum prescription supplies damping in the past. Equivalently, use a positive regulator $d$ and then take the limit:The real part of $-2i$ times this integral is $4/(c_s^3{K}^3)$. Including the six Wick contractions yields the cubic time-derivative curvature bispectrumThe resulting primordial bispectrum is a homogeneous function of degree $-6$ in the wavenumbers, as expected for scale-independent background parameters.

The printed limit $c_s\ll 0$ is not a physical sound-speed limit for these modes. The normalization and positive-frequency choice assume $c_s>0$; a negative value is outside the positive sound-speed convention used by these modes. The likely intended limit is $0<c_s\ll 1$, which is analyzed explicitly here rather than silently substituted for the PDF.

To compare primordial non-Gaussianity, normalize by the power spectrum, not only by the raw three-point amplitude. HereAt an equilateral bispectrum configuration, the exact normalized amplitude isThus $c_s\to 1$ gives an amplitude of order $ϵ$, consistent with weak non-Gaussianity in single-field slow-roll inflation. For $0<c_s\ll 1$, this vertex gives an amplitude of order one, enhanced relative to the slow-roll limit but not diverging as $c_s^{-2}$ after power-spectrum normalization. At fixed $H,ϵ,M_{\mathrm{Pl}}$ the raw bispectrum does scale as $c_s^{-2}$, but $P_{\zeta }^2$ has the same scaling. The often-quoted larger small-sound-speed amplitudes of more general inflationary actions require their actual interaction coefficients and other vertices; they do not follow from the lone supplied vertex. These conclusions concern its tree contribution within the stated weak-perturbation approximation.

The bispectrum shape function isAt fixed perimeter, the arithmetic-geometric mean inequality bounds this by $1/27$, attained only at an equilateral bispectrum configuration. In a squeezed bispectrum configuration with $k_1=q\ll k_2\simeq k_3=k$, it behaves as $q/(8k)\to 0$. At a flattened bispectrum configuration $(k,k,2k)$ it is finite, $1/32$, with no boundary singularity. Thus the shape prefers equilateral triangles and suppresses squeezed triangles, with a broad finite flattened contribution. This is a shape preference, not a complete signal-to-noise calculation: actual observational weighting also depends on cosmological transfer functions, covariance and accessible triangles. If $C=0$, this vertex produces no bispectrum at any configuration; if $C<0$, the signed shape reverses while its absolute triangle preference is unchanged.

Use the comoving curvature perturbation mode normalization. Take $H,ϵ,c_s$ to be effectively constant at leading slow-roll approximation order, with $H,ϵ,c_s>0$, and use the Bunch-Davies vacuum selected by the early-time $i\epsilon$ contour. The regulator $\epsilon$ is distinct from the Hubble slow-roll parameter $ϵ$.

Let $C=ϵ+1-c_s^2$ and use the Fourier transform convention $\zeta (\mathbf{x})=\int d^{3}k\zeta (\mathbf{k})e^{i\mathbf{k}\cdot \mathbf{x}}/(2\pi {)}^3$. Conversion to conformal time givesThe last sign follows from substituting $a=-1/(H\tau )$, which is positive for $\tau <0$; dropping the explicit minus sign in that substitution would reverse the final primordial bispectrum.

At the order needed by the in-in formalism, the unequal-time Wick contraction isIt uses the Gaussian quantum vacuum two-point function with the displayed operator ordering; an arbitrary classical unequal-time covariance cannot replace it. For the connected three-point function, each external field contracts with one of the three fields at the vertex. There are $3!=6$ such Wick contractions, all equal. The spatial integral conserves total momentum and leaves one factor $(2\pi {)}^3{\delta }^{(3)}({\mathbf{k}}_1+{\mathbf{k}}_2+{\mathbf{k}}_3)$. The other nine pairings are disconnected tadpole terms supported on a zero external momentum; they are excluded from the connected primordial bispectrum, equivalently by subtracting the one-point function or normal ordering the vertex.

The free Gaussian odd correlator is zero. Applying the stated first-order in-in formalism therefore givesBy definition, the full connected correlator is $(2\pi {)}^3{\delta }^{(3)}({\sum }_r{\mathbf{k}}_r)B$. The assumptions are a Gaussian Bunch-Davies vacuum, the given leading cubic vertex, tree order, the late-time limit $\tau \to 0^{-}$, effectively constant background parameters, and connected nonzero external modes. No symmetry factor $1/3!$ is inserted into the given Hamiltonian.

Absorb $k_{B}T$ into the free energy, so the statistical weight is $e^{-F}$. Under the Gaussian fixed point rescaling, invariance of the gradient energy gives the engineering dimensionThe engineering dimensions of ${\varphi }^2,{\varphi }^4,{\varphi }^6$ are respectively $1,2,3$. Their coupling eigenvalues are $d-[{\varphi }^{2m}]$, hence $2,1,0$. The quadratic and quartic interactions are relevant operators; the sextic interaction is a marginal operator at this order.

The three steps of the momentum-shell renormalization group are:

Here $b=e^s>1$, so increasing $s$ means flowing toward longer distances. At leading order, the fast field has a Gaussian distribution with covarianceThis perturbative integration requires $k^2+{\mu }_0^2>0$ on the shell. By the Wick theorem, odd fast-field moments vanish and $⟨{\eta }^4⟩=3I_b^2$, $⟨{\eta }^6⟩=15I_b^3$. Therefore the first term of the cumulant expansion of a coarse-grained free energy isThis is the quartic interaction generated by a sextic interaction. Including the rescaling, the requested flow isIn spherical coordinates, $I_b=(2{\pi }^2{)}^{-1}{\int }_{\Lambda /b}^{\Lambda }{k}^2/(k^2+{\mu }_0^2)dk$. In particular, for ${\mu }_0^2=0$,The same integration shifts the quadratic coefficient by $90{\lambda }_0{I}_b^2$ before rescaling and adds an irrelevant constant to the free energy. Thus setting the bare quartic coupling to zero does not place the theory on the tricritical critical surface.

To calculate the sextic renormalization-group beta function to second order, let $V={\lambda }_0\int d^{3}x(\phi +\eta {)}^6$. The required cumulant expansion isThe subtraction removes disconnected contributions. The term $20{\lambda }_0{\phi }^3{\eta }^3$ at each vertex supplies the primitive sextic correction. The Wick contractions giveThe first term contains three propagators connecting the two vertices; its factor $6$ is the number of bijections between the three fields at each vertex. The second contains a local tadpole diagram at each vertex and one connecting propagator. The primitive contribution to the free energy is thereforeExpanding the slow fields about a common point gives a negative local sextic correction proportional towhere ${\chi }_b$ restricts a momentum to the shell. This connected Feynman diagram has two independent loop momenta. Other Wick contractions also generate lower interactions and tadpole diagrams; they must be included or subtracted consistently when fixing the tricritical critical surface.

To extract the logarithmic tricritical sextic beta function, use normal ordering to remove local tadpoles and tune the quadratic and quartic relevant operators. The massless three-dimensional propagator is $C(z)=1/(4\pi |z|)$. The same three-line contraction over short separations $a<|z|<ba$ givesIts positive coefficient and the negative sign of the second cumulant expansion term establish the sign of the flow. Matching this logarithm to a renormalized coupling, or performing consistent iterated shell integration, givesThe precise normalization of $C$ is not needed. In a sharp momentum cutoff scheme, simply forcing all three internal momenta into one infinitesimal shell does not extract this two-loop logarithm: the shell restrictions remove its linear-in-$ds$ phase space. The finite-shell cumulant and the logarithmic renormalized calculation above are distinct stages of the computation. Field normalization corrections at order ${\lambda }^2$ enter this renormalization-group beta function only at order ${\lambda }^3$.

A positive, stabilizing sextic coupling is marginally irrelevant in three dimensions on the tricritical critical surface. Reversing the convention for $s$ reverses the sign of the renormalization-group beta function. A negative sextic coupling without higher stabilizing powers makes the potential unbounded; moreover, the untuned bare choice $g=0$ generally flows away through its generated relevant operator.

Assume $m^2>0$ and take $\lambda \to 0^{+}$. The free Gaussian integral gives $Z_0=\sqrt{2\pi }/|m|$, and the normalized free Gaussian measure has covariance $C=m^{-2}$. Expanding only the interaction exponential gives the asymptotic expansionIn this zero-dimensional scalar field theory, the path integral has become an ordinary integral, but its combinatorics are precisely those of sextic scalar field theory.

The free generating functional is $⟨e^{J\varphi }{⟩}_0=e^{C{J}^2/2}$. Differentiating it proves Wick theorem: odd moments vanish, while each even moment is the sum over pairings of its factors. Each Wick contraction contributes $C$, soThus the Feynman rules are a quantum field theory propagator $C$ for each internal line and a factor $-\lambda$ for each six-valent interaction vertex. The factorials $V!$ and $(6!{)}^V$ remove vertex and half-edge labels; the remaining weight is $(-\lambda {)}^V{C}^I/S$, where $S$ is the Feynman-diagram symmetry factor, the order of the automorphism group of the resulting diagram, including its half-edge symmetries. A tadpole diagram has an additional interchange symmetry of the two ends of its loop.

There is one empty Vacuum Feynman diagram, contributing $1$. At one interaction vertex, the six half-edges form three loops. Its Feynman-diagram symmetry factor is $2^{3}3!=48$, so its contribution is $-\lambda C^3/48$.

At two interaction vertices let $r$ count the lines joining them. Each interaction vertex has $(6-r)/2$ self-loops, so $r$ can only be $0,2,4,6$. These four possibilities exhaust all Vacuum Feynman diagrams at this order. Their Feynman-diagram symmetry factors areHere the leading $2$ exchanges the two interaction vertices, $r!$ permutes the joining lines, and the remaining factors permute and reverse the self-loops. The $r=0$ graph is disconnected and must be included in $Z/Z_0$; only a connected generating functional would discard it.

Consequently,Indeed ${\sum }_r{S}_r^{-1}=77/7680$, independently agreeing with $11!!/[2(6!{)}^2]$ from Wick theorem.

This is an asymptotic expansion, rather than a convergent perturbation series: the coefficient of ${\lambda }^V$ contains $(6V-1)!!/[V!(6!{)}^V]$ and grows too rapidly for a nonzero radius of convergence. For any fixed truncation order, the Taylor expansion remainder bound for $e^{-x}$ on $x\ge 0$, followed by the finite Gaussian moment integral, proves the stated remainder estimate. For negative real $\lambda$ the original integral diverges.

Use $⟨{\phi }_G(\mathbf{k}){\phi }_G({\mathbf{k}}^{\prime })⟩=(2\pi {)}^3{\delta }^{(3)}(\mathbf{k}+{\mathbf{k}}^{\prime })P_{\phi }(k)$. For nonzero external momenta, the quadratic term in local-type primordial non-Gaussianity isInsert it in one of the three external fields. Wick theorem supplies two connected Wick contractions at each insertion, whereas the mean subtraction removes the disconnected zero-mode contribution. Thus the leading primordial bispectrum isLinear propagation by the cosmological Poisson equation gives $P_m(k)=\alpha (k{)}^2{P}_{\phi }(k)$ and multiplies the primordial bispectrum by ${\alpha }_1{\alpha }_2{\alpha }_3$. Defining the matter three-point function by its momentum-conserving Dirac delta function times $B_m$, the primordial contribution to the matter bispectrum isThe nonlinear gravitational contribution from standard perturbation theory in cosmology is excluded here.