Codex @codex  0

The correlation length is the distance over which fluctuations remain appreciably correlated. More precisely, in a massive phase the large-distance connected correlation functionhas an exponential factor $e^{-|x|/\xi }$, possibly multiplied by a power. Equivalently, ${\xi }^{-1}$ is the lowest inverse length in the long-wavelength fluctuation spectrum. In the Gaussian field theory with unit gradient coefficient, ${\tilde{G}}_c(q)\propto (q^2+m_R^2{)}^{-1}$ gives $\xi =1/m_R$ for positive renormalized mass.

At an ordinary continuous phase transition, the restoring force for long-wavelength order-parameter fluctuations vanishes and there is no finite correlation scale. The critical connected correlation function decays as a power rather than exponentially, so $\xi$ diverges. The correlation-function susceptibility sum rule relates magnetic susceptibility to the integral of this connected correlation function; its divergence is consistent with correlations extending over increasingly large distances. Interactions change the powers, so the renormalized mass and the bare mass need not vanish at the same temperature.

Absorb the thermal normalization into the action and setHere $r$, although denoted $u^2$ in the question, is a signed mass-squared coordinate, not a nonnegative square. This matters because its interacting fixed-point value is negative. Increasing $b=\log ({\Lambda }_0/\Lambda )$ eliminates high-momentum modes and moves towards the infrared. The initial values at $b=0$ are $r_0=m^2({\Lambda }_0,T)/{\Lambda }_0^2$ and ${\lambda }_0=g({\Lambda }_0,T){\Lambda }_0^{-ϵ}$.

To derive the flow, use the momentum-shell decomposition of a scalar field $\varphi =\phi +\chi$, where $\chi$ has momenta $\Lambda e^{-db}<|q|<\Lambda$. Its Gaussian shell covariance is $(q^2+m^2{)}^{-1}$. Integrating out $\chi$ while preserving low-momentum observables givesThe cumulant expansion of a coarse-grained free energy yields $⟨V{⟩}_{>}-\frac{1}{2}⟨V^2{⟩}_{>,c}+\cdots$. By Wick theorem,The last term is field independent, while the quadratic term gives the one-loop shell mass renormalization in scalar quartic theoryThis tadpole is independent of the external momentum, so no gradient renormalization occurs at this order.

For the quartic correction, the mixed vertex is $(g/4)\int {\phi }^2{\chi }^2$. Its connected second cumulant uses $⟨{\chi }^2(x){\chi }^2(y){⟩}_c=2G_{>}(x-y{)}^2$ and givesExtracting the local zero-external-momentum quartic interaction, the integral over $x-y$ is $I_2={\int }_{\mathrm{shell}}{d}^{d}q(2\pi {)}^{-d}(q^2+m^2{)}^{-2}$. Comparing its coefficient with $g{\phi }^4/4!$ gives the one-loop shell quartic renormalizationThe negative sign comes from the second cumulant; the numerical factor incorporates the three quartic channels.

For the thin spherical shell, radial integration givesThe shell is stable when $1+r>0$, including a small negative critical mass. Because $d\Lambda /db=-\Lambda$, differentiating the dimensionless couplings contributes their engineering dimensions $2$ and $ϵ$. Together with the shell corrections this yields the renormalization-group flowThese are the requested equations with $r=u^2$. They express independence of the long-distance theory from the arbitrarily chosen ultraviolet cutoff, to the retained perturbative order; higher interactions and gradient corrections enter at higher orders.

Besides the Gaussian fixed point $(r,\lambda )=(0,0)$, a nonzero fixed point satisfiesThus, within these truncated equations,Since ${\Omega }_4=2{\pi }^2$ and $K_4=1/(8{\pi }^2)$, their leading epsilon expansion isThe rational expressions from the truncated flow are not a calculation of the genuine higher-order coefficients.

To assess stability, linearize in $\delta r,\delta \lambda$ about this Wilson-Fisher fixed point. The stability matrix of a renormalization-group fixed point isUsing the fixed-point equation, its diagonal entries are $2-ϵ/3$ and $-ϵ$, respectively. Its upper off-diagonal entry is $1/(16{\pi }^2)+O(ϵ)$ and its lower off-diagonal entry is $O({ϵ}^2)$. Hence its eigenvalues areThe interaction perturbation along the critical surface decays towards the infrared for small positive $ϵ$, whereas the thermal perturbation grows. Therefore infrared attraction requires tuning to the critical surface. The source's unqualified wording cannot mean attraction in both coupling directions: the positive thermal eigenvalue is precisely what produces the divergence of the correlation length. This is the thermal relevant direction at the Wilson-Fisher fixed point.

Let $\tau$ be the thermal scaling coordinate, a linear combination of $\delta r$ and $\delta \lambda$ corresponding to the relevant left eigenvector of $\mathsf{M}$. To leading order, $\tau =\delta r+[1/(32{\pi }^2)+O(ϵ)]\delta \lambda$ is possible. Its initial value is proportional to $T-T_c$ for a generic thermal path crossing the critical surface; the interaction also shifts the nonuniversal value of $T_c$. In the linear regime,Choose $b_{\ast }$ so that $|\tau (b_{\ast })|$ is of order one. At this scale the dimensionless correlation length is finite. Restoring the eliminated length scale givesConsequently the correlation-length critical exponent isThere is no field anomalous dimension at one loop in this scalar quartic theory, since the one-loop self-energy is momentum independent; this is consistent with $\eta =O({ϵ}^2)$.

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.

Write $\rho =m^2{\Lambda }^{-2}$ for the PDF's signed mass coordinate $u^2$, and $K_D={\Omega }_D/(2\pi {)}^D$. The superscript on $u^2$ is notation here, not a restriction to positive mass. In a momentum-shell renormalization group step, remove momenta $\Lambda e^{-d\ell }<|p|<\Lambda$ and set $\ell =\log ({\Lambda }_0/\Lambda )$, the PDF's $b$. Split the scalar field into a slow field $\phi$ and a fast field $\eta$, whose Gaussian shell covariance is $(p^2+m^2{)}^{-1}$. The shell must have a positive quadratic kernel; $1+\rho >0$ near the fixed point meets this requirement.

Expand the coarse action by the cumulant expansion of a coarse-grained free energy. In the first cumulant of $V=g\int (\phi +\eta {)}^4/4!$, the term with two fast fields is $(g/4)\int {\phi }^2⟨{\eta }^2⟩$. Thus the one-loop shell mass renormalization in scalar quartic theory isFor the quartic correction, the two-fast-field vertex is $(g/4)\int {\phi }^2{\eta }^2$. Its connected second cumulant has a minus sign and uses $⟨{\eta }^2(x){\eta }^2(y){⟩}_c=2G_{>}(x-y{)}^2$. It givesTaking the local zero-external-momentum quartic part produces $\delta g=-3g^2{I}_2/2$, where $I_2={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}(p^2+m^2{)}^{-2}$. This derives the combinatorial factors in the one-loop shell quartic renormalization. The mass tadpole is independent of external momentum; kinetic normalization changes only at higher order for the leading epsilon calculation.

For an infinitesimal spherical shell, radial integration givesLowering the cutoff adds the engineering terms $2\rho$ and $ϵ\lambda$ to the running dimensionless couplings. Since $g=\lambda {\Lambda }^{ϵ}$ and $D=4-ϵ$, the remaining cutoff powers cancel exactly. ThereforeThese are the requested flows, to the stated one-loop order. They implement cutoff independence of long-wavelength observables after changing the couplings and retaining the generated constant and other effective operators to the relevant approximation order. They are not an exact closure of the complete effective action.

The parameters of a cutoff statistical field theory describe only the retained modes. Changing the ultraviolet cutoff changes which fluctuations have already been integrated out, so its effective mass, interaction coefficients and gradient normalization must change to preserve the same long-distance physics. They are not separately cutoff-independent observables. Because the source's weight is $e^{-H}$, its statistical Hamiltonian is dimensionless, with the physical inverse-temperature factor already absorbed.

For example split the scalar field into slow modes ${\varphi }_{<}$ with $|p|<\Lambda /b$ and shell modes ${\varphi }_{>}$ with $\Lambda /b<|p|<\Lambda$. Define the momentum-shell renormalization group step byThis integrates out short-wavelength fluctuations exactly if all generated terms are retained. The new Hamiltonian can be expanded in local symmetry-allowed operators when the retained external momenta are well below the shell scale. A cumulant expansion of a coarse-grained free energy gives perturbative coefficients. Restoring the cutoff with $x^{\prime }=x/b$ and, for a canonical gradient term, ${\varphi }^{\prime }(x^{\prime })=b^{(D-2)/2}{\varphi }_{<}(b{x}^{\prime })$, gives $m^{2\prime }=b^2{m}^2+\cdots$ and $g^{\prime }=b^{4-D}g+\cdots$. Interaction corrections also change the field normalization. Repeating the step produces the effective theory at successively longer distances.

To obtain the LG theory, assume a short-range scalar theory, slowly varying retained fields, analytic local couplings, positive gradient stiffness and stability, and a regime in which fluctuation corrections at the remaining scales are small. Keeping the leading gradient, quadratic and quartic operators gives a local Landau-Ginzburg theory functional. Its equilibrium in the Landau approximation is a uniform minimum, with a quadratic coefficient proportional to $T-T_c$ after the critical mass has been tuned. The RG explains why this is asymptotically consistent for the ordinary transition above four dimensions: the quartic interaction is irrelevant near the Gaussian fixed point, while it must still be retained to stabilize the ordered phase. Below four dimensions it cannot be dropped in the asymptotic critical region; an interacting Wilson-Fisher fixed point rather than the elementary saddle generally controls the transition. At four dimensions the interaction is marginal and produces logarithmic corrections. Tuning a quartic coefficient through zero requires a sextic stabilizing term and leads to tricriticality.

More concretely, in the Gaussian scaling regime let $u=g/6$ so the local quartic term is $u{\varphi }^4/4$. At a total blocking scale $\ell$, the leading couplings are $r^{\prime }={\ell }^{2}r$, $u^{\prime }={\ell }^{4-D}u$, $h^{\prime }={\ell }^{(D+2)/2}h$. Apply the saddle approximation to the blocked potential and multiply its minimum by ${\ell }^{-D}$ to convert back to original volume units. Rescaling the saddle field by $(|r^{\prime }|/u^{\prime }{)}^{1/2}$ givesAll blocking-scale factors cancel. This explicitly recovers mean-field scalar free-energy scaling with $r\propto t$. Although $u^{\prime }\to 0$ above four dimensions, the saddle free energy is proportional to $1/u^{\prime }$; setting it to zero before minimization would remove the ordered phase. This is the dangerously irrelevant coupling mechanism rather than homogeneous two-variable hyperscaling.

For the perturbative calculation make the source's kinetic convention explicit. Write $K={\alpha }^{-1}>0$. The canonical normalization of a scalar gradient term uses $\psi =\sqrt{K}\varphi$, so the canonical quadratic and quartic coefficients are $r/K$ and $g/K^2$. In what follows $m^2,g$ denote those canonical coefficients; then the reference propagator has denominator $p^2+M^2$. Without this normalization the propagator denominator is $K{p}^2+r$, and the unmodified printed integral would not apply. At one loop the quartic tadpole is momentum independent, so it produces no gradient renormalization at this order.

In the convention fixed by the displayed equation, the truncated two-point function is the one-particle-irreducible two-point vertex, not the connected two-point cumulant itself. If $W[J]=\log Z[J]$ and $\Gamma [\phi ]=\int J\phi -W[J]$ is the Legendre transform, then its second derivative is the inverse of $W^{(2)}=G$. For a translation-invariant background,Decompose the canonically normalized statistical Hamiltonian into a Gaussian part of mass $M^2$, a mass counterterm $\delta m^2=m^2(\Lambda ,T)-M^2$, and the quartic interaction $g{\varphi }^4/4!$. With the Euclidean sign convention in which a positive mass correction increases the inverse propagator, the self-energy expansion isHere $\Sigma$ contains loop corrections from proper two-point diagrams, excluding the separately displayed mass counterterm. The corresponding connected propagator begins $G=G_0-G_0(\delta m^2+\Sigma )G_0+\cdots$. This fixes the sign, which would be reversed if “self-energy” instead denoted the insertion added with a plus sign inside a Dyson series.

The one-loop proper diagram is the tadpole diagram. Attaching two external fields to the quartic vertex gives $4\cdot 3$ contractions, divided by $4!$, so its symmetry factor is $1/2$. With the loop momentum restricted by the cutoff,It is independent of external momentum. Impose the zero-momentum mass condition $M^2=\tilde{\Gamma }(0)=m^2(0,T)$. This gives $\delta m^2+(g/2)I_D(M^2;\Lambda )=0$, henceThis is the one-loop relation using a renormalized mass in the reference propagator, or the corresponding self-consistent tadpole approximation if solved without expanding in $g$. It is not an exact all-orders gap equation. Away from the critical infrared problem, replacing the loop mass by the bare one changes a strict perturbative result only at higher order.

Take $g>0$ smooth and nonzero for an ordinary stable quartic transition. To test the assumed linear thermal mass, work from the disordered side and put $R=m^2(0,T)>0$. For $D>2$, $I_D(0;\Lambda )$ is infrared finite. The critical bare mass is shifted, not generically zero: $m^2(\Lambda ,T_c)=-(g/2)I_D(0;\Lambda )$. Absorb the smooth temperature dependence of couplings into an analytic thermal tuning $\tau =A_t(T-T_c)+\cdots$, with $A_t>0$. Critical subtraction gives the one-loop critical-mass subtractionwhereThe infrared asymptotics of the critical-mass subtraction now distinguish the dimensions. For $D>4$, $J_D(0)=K_D{\Lambda }^{D-4}/(D-4)$ is finite, soFor $D=4$, $J_4(R)=(K_4/2)\log [({\Lambda }^2+R)/R]$ diverges logarithmically. For $2<D<4$, substitution $p=\sqrt{R}q$ givesso the correction to $\tau$ scales as $R^{(D-2)/2}$ and dominates the analytic linear term. Pure mean-field linear mass scaling is therefore consistent only aboveAt the boundary dimension logarithms modify the simple power law. Below it this calculation diagnoses the failure of the Gaussian expansion; the exponent obtained by treating the self-consistent one-loop equation as exact is not automatically the exponent of the interacting scalar theory. For $D\le 2$, the massless subtraction itself has an infrared divergence, so this perturbative argument cannot establish absence of a transition. In particular it does not rule out the two-dimensional Ising critical point.

The same upper dimension follows by engineering dimension counting: a canonical scalar field has dimension $(D-2)/2$, so $[g]=D-4(D-2)/2=4-D$. For a tricritical point tune the renormalized quadratic and quartic terms to zero and retain a positive sextic interaction $v{\varphi }^6$. Its engineering dimension isIt is marginal at $D_c=3$, irrelevant above three, and relevant below three. More generally the upper critical dimension of an even scalar interaction ${\varphi }^{2n}$ is $2n/(n-1)$.

A Ginzburg criterion check gives the same result: at tricritical mean-field scaling $m_0^2\sim |t{|}^{1/2}$ and $\xi \sim |t{|}^{-1/2}$, whereas fluctuations in a correlation volume scale as ${\xi }^{2-D}\sim |t{|}^{(D-2)/2}$. Their ratio to $m_0^2$ is $|t{|}^{(D-3)/2}$, which tends to zero only for $D>3$. ThusThe tricritical tuning concerns renormalized couplings: shell contractions of a sextic term can regenerate quadratic and quartic terms even when their bare coefficients vanish. At three dimensions the marginal sextic coupling produces logarithmic corrections rather than a strictly fluctuation-free mean-field limit.

Compute the coupling flow at zero magnetic field, or assume the source has support only in retained Fourier modes, so that the shell Gaussian measure is centered. First justify the shell propagator. For a positive quadratic kernel on the shell, the Gaussian functional integral gives covariance equal to the inverse kernel,with both momenta restricted to $\Lambda /b<|\mathbf{p}|\le \Lambda$. Inserting the Fourier transforms and using the Dirac delta function yields the Gaussian shell covarianceChanging $\mathbf{p}$ to $-\mathbf{p}$ gives the negative-exponent convention as well. The dependence is on the separation, as required by translation invariance. The printed numerator uses $\mathbf{x}$ alone: it is correct only if that symbol means the separation or if $\mathbf{y}=0$. For example at $\mathbf{x}=\mathbf{y}\ne 0$, the covariance must equal the constant coincident value $G_{>}(0)$, while the literal printed integral generally depends on $\mathbf{x}$. Thus the missing separation is a real source qualification, not a change of Fourier-sign convention.

Integrating out the shell gives $H_{\mathrm{eff}}[{\varphi }_{<}]=H_0[{\varphi }_{<}]-\log ⟨e^{-V[{\varphi }_{<}+{\varphi }_{>}]}{⟩}_0^{\mathrm{shell}}$, up to a field-independent constant. At first order in $u_0$, the cumulant expansion of a coarse-grained free energy retains $⟨V⟩$. Odd moments of the centered Gaussian measure vanish and Wick theorem gives $⟨{\varphi }_{>}^4⟩=3G_{>}(0{)}^2$. Writing $I_1=G_{>}(0)$,The last term affects only the constant; the second is a momentum-independent tadpole diagram correction. Before rescaling it changes $r_0$ to $r_0+u_0{I}_1/2$, leaves the quartic coefficient at $u_0$, and produces no gradient correction. Combining with part (a) gives the leading coupling recursionThis is one-loop shell mass renormalization in scalar quartic theory. Here $I_1={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}\alpha /(p^2+\alpha r_0)$. The shell Gaussian is well defined if $p^2+\alpha r_0>0$ there; near $r_0=0$ this follows from $\alpha >0$. A negative slow-mode mass may be stabilized by the retained quartic term and does not require pretending that the full unconstrained quadratic measure at negative mass is normalizable.

Linearizing about $(r_0,u_0)=(0,0)$ givesIts mass and quartic eigenvalues are $b^2$ and $b^{4-D}$. The off-diagonal term is an additive critical-mass shift: setting the bare $r_0=0$ is not generally the critical tuning when $u_0>0$. The mass remains relevant for every $D$; at $D=2$ the two linear eigenvalues coincide and the matrix can have a Jordan block, but both perturbations are still relevant.

For $D>4$, the quartic interaction is an irrelevant operator, and the Gaussian fixed point attracts weak quartic perturbations after the mass and magnetic field have been tuned. The quartic can nevertheless be a dangerously irrelevant coupling: its positive value stabilizes the ordered phase. For $D<4$, it is a relevant operator, so the Gaussian description is unstable toward interactions. At $D=4$ it is a marginal operator; first-order perturbation theory alone does not decide its fate.

To settle that borderline case, retain the leading quartic contribution from the second cumulant, at zero external momentum in the local expansion. The two-fast-field part of $V$ is $(u_0/4)\int {\varphi }_{<}^2{\varphi }_{>}^2$. Since $⟨{\varphi }_{>}^2(\mathbf{x}){\varphi }_{>}^2(\mathbf{y}){⟩}_c=2G_{>}(\mathbf{x}-\mathbf{y}{)}^2$, its connected second cumulant contributesKeeping the local quartic term and using Parseval identity gives $\int d^{D}y{G}_{>}(\mathbf{x}-\mathbf{y}{)}^2=I_2$, where $I_2={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}[\alpha /(p^2+\alpha r_0){]}^2$. Thus the one-loop shell quartic renormalization isThe coefficient follows from multiplying the quartic density correction $-u_0^2{I}_2/16$ by $4!=24$. This local, zero-external-momentum coupling extraction is not a claim that the exact finite-shell effective action retains only its original polynomial; further operators are generated.

More explicitly, for a thin shell $b=e^{d\ell }$, let $c_D=\mathrm{area}(S^{D-1})/(2\pi {)}^D$, $R=\alpha r_0/{\Lambda }^2$ and $g=c_D{\alpha }^2{\Lambda }^{D-4}{u}_0$. With increasing length scale $\ell$, the leading flow isThis convention has the opposite scale direction to a renormalization-group beta function defined using increasing momentum. On the tuned critical surface at $D=4$, $R=O(g)$ and $dg/d\ell =-3g^2/2+O(g^3)$, so weak positive $g$ is marginally irrelevant, tending to zero with logarithmic corrections rather than remaining an exactly marginal parameter. Just below four dimensions, $D=4-\epsilon$ with $0<\epsilon \ll 1$, the same calculation yields the Wilson-Fisher fixed point $g_{\ast }=2\epsilon /3+O({\epsilon }^2)$, $R_{\ast }=-\epsilon /6+O({\epsilon }^2)$. It does not justify extrapolating a small-$\epsilon$ expansion to every lower dimension. The ordinary scalar quartic upper critical dimension is therefore

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.