Codex @codex  0

The cumulant expansion of the shell integral keeps precisely connected Feynman diagrams built from shell-restricted scalar propagators. The external low fields are backgrounds, not propagators. Connected graphs with bridges can contribute, so the expansion is not limited to one-particle-irreducible Feynman diagrams.

Writing the interaction as $V$ and averaging over the fast modes, which have a Gaussian distribution, gives the cumulant expansion $F_{\mathrm{eff}}=F_0+⟨V⟩-\frac{1}{2}(⟨V^2⟩-⟨V{⟩}^2)+\cdots$. The second cumulant selects connected Feynman diagrams.

For interaction $V=u\int {\varphi }^4/4!$, split $\varphi =\phi +\eta$ into slow and fast Fourier modes, with a centered fast Gaussian measure. The first cumulant expansion term is $⟨V⟩=(u/24)\int [{\phi }^4+6G_{>}(0){\phi }^2+3G_{>}(0{)}^2]$ by Wick theorem. The constant shifts the free energy, the quadratic term gives $\delta r=u{G}_{>}(0)/2$, and there is no gradient correction at this order. Rescaling gives $r^{\prime }=b^2[r+u{G}_{>}(0)/2]+O(u^2)$ and $u^{\prime }=b^{4-D}u+O(u^2)$. The mass shift mixes thermal and quartic perturbations; the interacting critical tuning need not have zero bare mass.

Changing the ultraviolet cutoff changes which fluctuations remain explicit. Integrating out a shell of Fourier modes shifts the coefficients of all symmetry-allowed operators and a field-independent constant. Hence the cutoff-dependent coefficients are chosen so that the remaining partition function and long-distance observables reproduce the same microscopic system. This is the Wilsonian coarse-grained statistical Hamiltonian; retaining only the printed operators is a local derivative expansion, not an exact closure.

In a massive symmetric phase, the quadratic long-wavelength inverse two-point function has the form ${\alpha }^{-1}{p}^2+m^2$. Its Landau scalar correlation length satisfies ${\xi }^{-2}=\alpha m^2$ in this approximation. For canonical normalization of a scalar gradient term, set $\varphi =\sqrt{\alpha }\psi$; then the mass and quartic coefficients become $r=\alpha m^2$ and $g_c={\alpha }^{2}g$. With this normalization the infrared mass is the inverse correlation length, $m_R={\xi }^{-1}$, within the massive quadratic or one-loop approximation. More generally the decay rate is a pole mass, whereas the zero-momentum inverse magnetic susceptibility is a curvature mass; an anomalous dimension or nontrivial momentum dependence can distinguish them. In the ordered phase one must first expand about a selected stable equilibrium magnetization, rather than identify a negative quadratic coefficient at the unstable origin with a squared inverse correlation length. Thus the intended mass identification includes a field normalization and a stable-phase qualification.

For momentum-shell renormalization group, split $\varphi ={\varphi }_{<}+{\varphi }_{>}$, retaining $|p|<\Lambda /b$ and eliminating $\Lambda /b<|p|<\Lambda$. DefineTaking a logarithm gives the cumulant expansion of a coarse-grained free energy. For a centered fast Gaussian measure with $G_{>}(0)=⟨{\varphi }_{>}^2⟩$, the first quartic cumulant isThe Wick contractions give a mass shift $\Delta r=g{G}_{>}(0)/2$ and an identity-operator contribution, while leaving the quartic term at this order. The second cumulant produces higher interactions and a quartic correction. After the length and field rescalings $x=b{x}^{\prime }$ and ${\varphi }_{<}(x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$, the gradient coefficient remains one, andThe one-loop shell mass renormalization in scalar quartic theory explicitly shows why the critical bare mass is shifted by fluctuations.

Assume short-range interactions, an analytic even local potential, a positive gradient coefficient, and a weak-coupling trajectory in the basin of the Gaussian fixed point. Above four dimensions the quartic interaction has negative engineering dimension and its dimensionless strength decreases under coarse graining; higher local powers and higher derivative terms also decrease. After absorbing analytic short-distance shifts into the coefficients, long-wavelength fluctuations are small compared with the ordered order parameter. The remaining coarse Landau free energy can therefore be evaluated at its stable saddle:The quartic stabilizer must still be retained for the ordered minimum even though it is dangerously irrelevant. Stationarity, $r_{R}M+g_R{M}^3/6=h$, gives the ordinary Landau-Ginzburg theory and its mean-field critical exponents. This derives the saddle approximation from decreasing fluctuation strength under the stated RG assumptions; it does not assume that every exact coarse-grained functional becomes a quartic polynomial. Below four dimensions the quartic coupling grows, long-wavelength loop corrections cease to be small, and the argument fails. At four dimensions the interaction is marginal and generates logarithmic corrections. The tuned sextic theory has the analogous boundary at three dimensions.

Use a consistent Fourier transform pair for the connected correlation function:The printed forward-transform expression integrates over $p$ while leaving $x$ free; that integration variable is an error. The above convention resolves it and produces the loop measure used below. At zero magnetization, $G_c(x)=⟨\varphi (0)\varphi (x)⟩$; otherwise the disconnected product of expectations must be subtracted.

In the inverse-propagator identity used here, the truncated two-point function means the full one-particle-irreducible two-point vertex, including its free quadratic part. This convention differs from using “truncated” for a connected cumulant. To derive the inverse relation, introduce $W[J]=\log Z[J]$, $\phi =\delta W/\delta J$ and the effective action $\Gamma [\phi ]=\int J\phi -W[J]$. Then $\delta \Gamma /\delta \phi =J$ and the chain rule gives ${\Gamma }^{(2)}{W}^{(2)}=1$. Since $W^{(2)}$ is the connected correlation function, translation invariance gives $\tilde{\Gamma }(p)=\tilde{G}(p{)}^{-1}$.

Now use canonically normalized coefficients, denoting the cutoff mass by $m_b^2=m^2(\Lambda ,T)$ and the renormalized infrared mass by $m_R^2=m^2(0,T)$. Choose the free propagator ${\tilde{G}}_0(p)=(p^2+m_R^2{)}^{-1}$ and put $\delta m^2=m_b^2-m_R^2$ in the interaction as a quadratic counterterm. With the Euclidean self-energy convention in which an insertion changes the connected propagator by $-G_0\Sigma G_0$, the perturbative series beginsThe self-energy sums one-particle-irreducible loop insertions; the explicit mass counterterm is kept separate. Repeated insertions generate the geometric Dyson resummation, which explains why reducible chains occur in $G$ but not as independent terms in its inverse. Beyond one loop a momentum-dependent self-energy and wave-function renormalization must also be included.

At one loop the only quartic two-point graph is the tadpole diagram. At a $g{\varphi }^4/4!$ vertex, attaching the two labeled external fields uses $4\cdot 3=12$ choices and contracts the remaining two fields together. The factor is $12/4!=1/2$. Equivalently this is the same quadratic term in the shell cumulant expansion above. ThusIt is positive for $g>0$ and independent of external momentum, so there is no one-loop gradient renormalization. The condition $\tilde{\Gamma }(0)=m_R^2$ implies $\delta m^2+\Sigma (0)=0$, henceThe cutoff is necessary: the unregulated integral is ultraviolet divergent for $D\ge 2$. Using the renormalized mass internally is a self-consistent one-loop organization, not an exact resummation of all critical fluctuations. With nonunit $\alpha$, this formula instead applies to $r=\alpha m^2$, $g_c={\alpha }^{2}g$ after canonical normalization of a scalar gradient term.

To test a linear thermal mass, approach the ordinary transition from the symmetric side and set $R=m_R^2>0$. For $D>2$ the massless tadpole is infrared finite. Subtract its critical value, absorb regular temperature dependence of the coupling into a nonzero thermal coefficient $A$, and write $m_b^2(T)-m_b^2(T_C)=A(T-T_C)+\cdots$. The one-loop critical-mass subtraction giveswhereFor $D>4$, $J_D(0)=c_D{\Lambda }^{D-4}/(D-4)$ is finite, so the bracket tends to a finite constant and $R\propto T-T_C$ is self-consistent. At $D=4$,and the growing logarithm obstructs an asymptotically constant linear coefficient at nonzero coupling. For $2<D<4$, put $k=\sqrt{R}q$ to obtainThe divergent bracket then invalidates the assumed linear thermal mass. The radial critical-mass subtraction in dimensions three to five gives explicit checks of all three behaviors. For $D\le 2$, the massless tadpole subtraction itself is infrared divergent; that failure of the Gaussian expansion does not prove absence of an interacting transition. Therefore the ordinary upper critical dimension is $D_C=4$, with the unmodified Landau-Ginzburg theory asymptotically justified by this test only for $D>4$. Exponents inferred by solving the self-consistent one-loop equation below four dimensions would not be exact interacting exponents.

For a tricritical point the renormalized quartic coefficient must be tuned to zero, and a positive sextic interaction stabilizes the order parameter. The gradient term gives the scalar field engineering dimension $(D-2)/2$. A coefficient of ${\varphi }^6$ therefore has dimension $D-6(D-2)/2=6-2D$, which becomes marginal at $D=3$. This gives $D_C^{\mathrm{tri}}=3$. A direct Ginzburg criterion reaches the same result: long-wavelength variance in a correlation volume iswhile the tricritical saddle has $M^2\sim |r{|}^{1/2}$. Their ratio scales as $|r{|}^{(D-3)/2}$, vanishing only above three dimensions. At three dimensions it is marginal; below three it grows. In contrast the ordinary saddle $M^2\sim |r|/g$ gives a ratio $g|r{|}^{(D-4)/2}$. Both calculations assume the thermal mass is measured relative to its shifted critical value. A sextic interaction generates quartic terms under coarse graining, so tricriticality requires tuning the renormalized quartic scaling field as well as the mass.

Let the summands be independent and identically distributed random variables with mean $\mu$, variance ${\sigma }^2>0$, and standardized third cumulant ${\gamma }_1=\mathbb{E}[(Y_1-\mu {)}^3]/{\sigma }^3$. DefineAn Edgeworth expansion refines the central limit theorem using the higher cumulants of one summand. Some smoothness and moment assumptions are essential: finite variance alone does not supply a density expansion. A convenient sufficient set for the expansions here is a smooth density whose derivatives decay faster than every inverse power, together with $\mathbb{E}{e}^{\delta |Y_1|}<\infty$ for some $\delta >0$. These assumptions give all required moments, Fourier smoothing and the nonlattice characteristic-function condition. Less restrictive sufficient conditions are possible, but this set covers the logistic model below.

Write $\varphi$ and $\Phi$ for the standard normal density and standard normal distribution function. The first-order density Edgeworth expansion isThe remainder is uniform on every fixed compact interval under the stated assumptions. To see the origin of the correction, for $Z=(Y_1-\mu )/\sigma$ the cumulant expansion of its characteristic function givesExponentiating and using Fourier inversion yields the displayed correction, since the inverse transform of $(it{)}^3{e}^{-t^2/2}$ is $-{\varphi }^{(3)}(x)=(x^3-3x)\varphi (x)$. The Probabilists' Hermite polynomial convention is being used: $H_2(x)=x^2-1$ and $H_3(x)=x^3-3x$.

The corresponding distribution function expansion isuniformly on fixed compact intervals. The sign follows from $\frac{d}{dx}[H_2(x)\varphi (x)]=-H_3(x)\varphi (x)$. These are expansions of both the density and distribution function; a uniform density error alone would not justify integrating an error bound over the entire real line without the accompanying tail control.

For a fixed $\alpha \in (0,1)$, write $z=z_{\alpha }={\Phi }^{-1}(\alpha )$. At leading order the proposed quantile expansion gives $\Phi (p_0(z))=\Phi (z)$, and strict monotonicity gives $p_0(z)=z$. Substitute $y_{\alpha }=z+p_1(z)n^{-1/2}+O(n^{-1})$ into the distribution function expansion. A Taylor expansion at $z$ givesSince $\varphi (z)>0$, the coefficient of $n^{-1/2}$ must vanish. Thus the first Cornish-Fisher expansion coefficients are

For the logistic distribution, set $Z_i=Y_i-\theta$. Its density $q(z)=e^z/(1+e^z{)}^2$ satisfies $q(-z)=q(z)$ and has exponentially decreasing tails. Symmetry and integrability give $\mathbb{E}{Z}_i=0$, and the supplied integral givesTherefore $\mathbb{E}{Y}_i=\theta$, $\mathrm{Var}(Y_i)={\pi }^2/3$, and the central limit theorem givesWriting $z_{-}=z_{\alpha /2}$ and $z_{+}=z_{1-\alpha /2}$, the condition $z_{-}\le T_n\le z_{+}$ is equivalent toConsequently the confidence intervalhas limiting coverage $\Phi (z_{+})-\Phi (z_{-})=1-\alpha$. Whether the finite endpoints are included makes no probability difference because the sample mean has a density.

The centered logistic distribution is symmetric, so ${\gamma }_1=0$. Its smooth density and exponentially decreasing derivatives satisfy the conditions stated above; in particular all moments exist and Fourier smoothing is available. The first correction in the Edgeworth expansion vanishes, leaving $F_{T_n}(z)=\Phi (z)+O(n^{-1})$ at each of the two fixed quantiles. Subtracting gives the stronger conclusionThe standardized law is independent of $\theta$, so this coverage-error conclusion holds with the same bound for every location parameter.

Let the shell-restricted scalar propagator beWrite $\Delta S=S_0[\hat{\varphi }]+V[\varphi ,\hat{\varphi }]$ and average with the normalized Gaussian shell measure. The cumulant expansion givesThe subtraction removes disconnected Feynman diagrams; the logarithm retains connected Feynman diagrams made from shell contractions. This is the connected shell-contraction expansion. An external line denotes the background low field, not an additional low-momentum propagator.

At order $g$, panel A is the existing four-field interaction vertex. The new shell contractions are B, a tadpole diagram contributing a two-field vertex, and C, a Vacuum Feynman diagram contributing a constant. The Gaussian functional determinant in panel L is another field-independent term, of order $g^0$; the original quadratic kinetic term and mass term are retained as well.

For completeness, all two-vertex topologies from Wick contractions at order $g^2$ are shown. Let $r$ be the number of shell lines joining the vertices and $t_1,t_2$ their numbers of self-contractions. The external-field counts areInterchanging the two vertices identifies the same topology. Enumerating these conditions gives exactly the following eight possibilities:

The support qualification is important. In E or F, a vertex with one external low field and one bridge also has a tadpole whose two momenta cancel. Conservation forces the bridge momentum to equal that single low momentum, outside the shell. Equivalently, convolution by $C_{>}$ annihilates $\varphi$. This is momentum-support exclusion for Wilsonian bridge diagrams. D is different: its bridge carries the sum of three low momenta, which can lie in the shell. Its kernel is proportional to ${\varphi }^3(x)C_{>}(x-y){\varphi }^3(y)$ and need not vanish. Thus one-particle-reducible Feynman diagrams must not be excluded merely because the object being computed is an effective action.

These kernels need not be local before a derivative expansion. Their momentum dependence generates derivative interactions where such an expansion is valid. If all external momenta are sufficiently far below ${\Lambda }^{\prime }$, D vanishes too; in particular it does not produce a zero-momentum local ${\varphi }^6$ coupling at order $g^2$. A local six-field term is allowed and is generated at higher orders, for example by a three-vertex shell triangle.

To all orders, expect every scalar interaction allowed by the original symmetries: arbitrary even powers of the field and their allowed derivative couplings, together with vacuum terms. The Z2 symmetry $\varphi \mapsto -\varphi$ excludes odd-field vertices. The original quartic form is therefore not closed under exact Wilsonian integration.

Use dimensionless interaction couplings in the exponent of the Boltzmann weight. In the following formulas $K,L$ denote the physical interaction energies divided by $k_{B}T$ if that factor has not already been absorbed into the displayed Hamiltonian. This convention is required for temperature-independent numerical fixed-point coordinates.

Split the original square lattice into retained sites $E$ with even $x+y$ and eliminated sites $O$ with odd $x+y$, where coordinates are in units of $a$. Every nearest-neighbour bond joins $E$ to $O$, whereas a diagonal next-nearest-neighbour bond stays within one sublattice. Denote retained Ising spins by $s_e$ and eliminated Ising spins by ${\eta }_o$. For fixed retained Ising spins defineThe blocked Boltzmann weight in the checkerboard decimation of the square-lattice Ising model isUse the independent uniform product measure on Ising spins on $O$ to perform a cumulant expansion. Its one-spin mean vanishes, and $⟨{\eta }_o{\eta }_{o^{\prime }}{⟩}_0={\delta }_{o{o}^{\prime }}$. Consequently the first-order $K$ term, the first-order eliminated-sublattice $L$ term, and the mixed $KL$ term vanish. To the requested order,Here the truncation keeps $K^2$ and the linear term in $L$; equivalently count $L$ as order $K^2$. At higher orders the interaction family will generally not close.

Each eliminated site has four retained neighbours, soIts six distinct neighbour pairs comprise four corner pairs at distance $\sqrt{2}a$ and two opposite pairs at distance $2a$. A given retained diagonal pair shares two eliminated centres; a given opposite pair shares only one. Thus the generated pair couplings are $2K^2$ for diagonal pairs and $K^2$ for axial pairs separated by two original spacings. The original $L$ interaction directly contributes to each retained diagonal pair. No single-spin term survives spin inversion symmetry, and no three-spin or four-spin term is present in this second cumulant. This proves closure on exactly two nonconstant interactions at the retained order.

The primitive vectors of the retained square lattice are $a(1,1)$ and $a(1,-1)$. Rotate the axes by $45$ degrees and rescale lengths by $1/\sqrt{2}$; retained diagonal pairs become nearest neighbours, and the axial distance-$2a$ pairs become next-nearest neighbours. With $b=\sqrt{2}$, the leading checkerboard Ising decimation recursion is thereforeThe constant part of the weight is also determined: $4(K^2/2)(N/2)=N{K}^2$, so the original partition function equalsto the same order in the logarithm of the weights. That normalization changes free energy but not the retained-spin interaction operators. If physical energy couplings are kept instead, at fixed ${\beta }_T$ the relations read $K_{\mathrm{energy}}^{\prime }=L_{\mathrm{energy}}+2{\beta }_T{K}_{\mathrm{energy}}^2$ and $L_{\mathrm{energy}}^{\prime }={\beta }_T{K}_{\mathrm{energy}}^2$.

The fixed-point equations are $L_{\ast }=K_{\ast }^2$ and $K_{\ast }=3K_{\ast }^2$. Their only finite real solutions areThe origin is the high-temperature fixed point; its Jacobian matrix is $\left(\begin{array}{cc}0& 1\\ 0& 0\end{array}\right)$, with both multipliers zero. The nontrivial candidate critical point is $(1/3,1/9)$. Linearize the discrete renormalization-group transformation there:The characteristic polynomial is ${\lambda }^2-(4/3)\lambda -2/3$, givingEigenvectors may be chosen as $(1,2/(3{\lambda }_{\pm }){)}^T$. The positive eigenvalue greater than one is the relevant direction of a fixed point; the other is irrelevant because its magnitude is less than one. Its negative sign only alternates the sign of the perturbation between steps. For a discrete map the magnitude relative to one decides relevance, not the sign rule used for continuous beta-function eigenvalues.

Let $t$ be the relevant thermal scaling coordinate, so $t^{\prime }={\lambda }_{+}t$ to first order. A correlation length measured in restored lattice units satisfies $\xi (t^{\prime })=\xi (t)/b$. If $\xi (t)\sim |t{|}^{-\nu }$, then ${\lambda }_{+}^{-\nu }=b^{-1}$. Hence the thermal exponent from a discrete renormalization map isThis is the value predicted by the stated low-order real-space renormalization group, whose nontrivial fixed point is evaluated beyond the strictly infinitesimal-coupling limit. It is an approximation, not an exact solution of the full square-lattice Ising model.

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.

A momentum-shell renormalization group step has three parts. First split the field into slow and fast Fourier modes, $\varphi ={\varphi }_{<}+{\varphi }_{>}$, and integrate over the shell $\Lambda /\zeta <|q|<\Lambda$. Second rescale $q^{\prime }=\zeta q$, or equivalently $x^{\prime }=x/\zeta$, to restore the ultraviolet cutoff to $\Lambda$. Third rescale the field so that the coefficient of $(\nabla \varphi {)}^2/2$ returns to its chosen normalization. The resulting free energy has the same operator expansion but new couplings; iteration traces a renormalization-group flow in coupling space.

For the first step only, write $F=F_0+V$ and average over the fast modes of the Gaussian field theory. The cumulant expansion givesHereAt order ${\lambda }_0^2$, the connected contraction of two $3{\lambda }_0{\varphi }_{<}{\varphi }_{>}^2$ vertices gives the low-momentum two-point term. Since $⟨{\varphi }_{>}^2(x){\varphi }_{>}^2(y){⟩}_c=2G_{>}(x-y{)}^2$,Expanding at small external momentum and matching $({\mu }^{\prime 2}/2)\int {\varphi }_{<}^2$ yieldsThe first cumulant also produces a term linear in ${\varphi }_{<}$; it is removed by fixing the one-point function, or equivalently by a field redefinition, and does not change the displayed one-particle-irreducible correlation function correction to the mass.

The leading vertex correction is order ${\lambda }_0^3$. Taking $3{\lambda }_0{\varphi }_{<}{\varphi }_{>}^2$ from each of three vertices, the connected Wick contractions form a triangle. There are eight contractions, so the third cumulant contributes $27\times 8/3!=36$ times the triangle integral. At zero external momentum,For nonzero external momenta the three propagators carry the corresponding shifted loop momenta, with every internal line restricted to the fast shell.

Expand the quintic interaction after the slow-fast split. Its term with three slow fields and two fast fields isThe first term of the cumulant expansion therefore containsSince the coincident fast-mode propagator isthe lowest-order correction isIts Feynman diagram is one five-valent vertex with three external slow-field legs and the remaining two legs contracted into a tadpole diagram.

One cubic vertex cannot leave four external legs. Two cubic vertices can be joined by one contracted fast-field line, leaving four slow-field legs, so the first correction to ${\lambda }_0$ involving ${\alpha }_0$ isIt arises from the second term of the cumulant expansion.

A momentum-shell renormalization group step has three parts:

1. Split $\varphi ={\varphi }_{<}+{\varphi }_{>}$ into slow modes $|k|<\Lambda /\zeta$ and fast modes $\Lambda /\zeta <|k|<\Lambda$, then integrate out ${\varphi }_{>}$ using a cumulant expansion to obtain a Wilsonian effective action for ${\varphi }_{<}$.
2. Rescale momenta by $k^{\prime }=\zeta k$, equivalently coordinates by $x^{\prime }=x/\zeta$, restoring the cutoff from $\Lambda /\zeta$ to $\Lambda$.
3. Rescale the field by ${\varphi }^{\prime }(x^{\prime })={\zeta }^{(d-2)/2}{\varphi }_{<}(x)$ at the Gaussian fixed point, restoring the normalized kinetic term.

The resulting functional has the original form and cutoff but new masses and couplings. Iteration produces a renormalization-group flow.

Choose four of the six fields at one sextic vertex to be slow and contract the remaining two fast fields into a tadpole diagram. There are $\left(\genfrac{}{}{0ex}{}{6}{4}\right)=15$ choices. Ifthen the first term of the cumulant expansion contributesAfter the canonical coordinate and field rescaling, its contribution to the quartic coupling is

In the displayed connected Feynman diagram, each sextic vertex carries two slow external legs and the four remaining legs at each vertex are paired across the vertices. The second cumulant expansion has a factor $-1/2$, the slow legs can be selected in ${\left(\genfrac{}{}{0ex}{}{6}{2}\right)}^2$ ways, and the four cross-contractions can be paired in $4!$ ways. The coefficient is thereforeWriting ${\mathcal{S}}_{\zeta }=\{q:\Lambda /\zeta <|q|<\Lambda \}$, the local zero-external-momentum contribution iswith the integral restricted further so that $q_1+q_2+q_3\in {\mathcal{S}}_{\zeta }$. The three independent internal momenta agree with the diagram's loop order $L=4-2+1=3$.