Codex @codex  0

Use a selected translationally invariant pure thermodynamic phase and define the connected lattice correlation functionAway from criticality its long-distance tail has an exponential factor $e^{-r/\xi }$, possibly multiplied by a power. Equivalently ${\xi }^{-1}=-{\lim }_{r\to \infty }{r}^{-1}\log |G(r)|$ when this limit exists; different directions can have different lengths on an anisotropic lattice. At a thermodynamic critical point the correlation length diverges and the tail becomes algebraic. Differentiate $M=N^{-1}⟨{\sum }_n{\sigma }_n⟩$ with respect to the physical field appearing in the energy Statistical Hamiltonian. With ${\beta }_{\mathrm{th}}=1/(k_{B}T)$, the correlation-function susceptibility sum rule isIn physical continuum coordinates the sum becomes $a^{-D}\int d^{D}rG(r)$. Differentiating with respect to the dimensionless source ${\beta }_{\mathrm{th}}h$ instead removes the explicit inverse-temperature factor. Below the transition, take a one-sided pure-phase limit rather than a macroscopic symmetric mixture.

A normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )$ assigns a probability or a delta constraint to a coarse configuration for each microscopic one, with ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$. Define the blocked Boltzmann weight bySumming over ${\sigma }^{\prime }$ proves exact equality of the partition functions. Repeated real-space renormalization group steps compose the kernels and integrate out shorter-scale variables. To retain equality one must keep all generated interactions and the field-independent constant; a finite truncation is an approximation. Long-distance observables are mapped through the corresponding coarse observables, not identified blindly with microscopic spins. In physical units,A subsequent coordinate rescaling restores the reference lattice spacing while expressing couplings in cutoff units.

Let $F(u,C)=-({\beta }_{\mathrm{th}}N{)}^{-1}\log Z=C+f(u)$ be free energy per microscopic site. Invariance of $Z$ givesThe constant transforms as $C_1=b^D{C}_0+c(u_0)$. Therefore, with $g(u)=b^{-D}c(u)$, the identity-operator contribution to renormalization-group free energy givesThe second equality follows by substitution and induction. The inhomogeneous term records eliminated-mode entropy, determinants, and normalization factors, all multiplying the identity operator; it is not an interaction with the retained spins. Strictly, this $f$ can contain regular background terms. After an analytic subtraction of an inhomogeneous renormalization recursion, the genuinely singular part satisfies the homogeneous scaling equation, except for possible logarithmic resonances.

A renormalization-group fixed point satisfies $R_b(u_{\ast })=u_{\ast }$. Linearize in scaling coordinates $s_i$: $s_i^{\prime }=b^{{\lambda }_i}{s}_i$. A relevant operator has ${\lambda }_i>0$, an irrelevant operator has ${\lambda }_i<0$, and a marginal operator has ${\lambda }_i=0$ and needs nonlinear analysis. These ${\lambda }_i$ are logarithmic scaling eigenvalues; the discrete Jacobian eigenvalues are $b^{{\lambda }_i}$, not ${\lambda }_i$. The critical surface is the stable manifold obtained by tuning all relevant coordinates to zero. A repulsive renormalization-group trajectory leaves the fixed point along a relevant direction as the coarse-graining length grows. Stable and unstable manifolds become curved away from linear order.

The left section has $h=0$, one irrelevant coordinate $w$ and a relevant thermal coordinate $t$. Its critical surface is $t=0$ in that section. With two relevant variables, the full critical surface has codimension two: the second section shows repulsion in the $(t,h)$ plane at $w=0$.

Assume an isolated fixed point with two relevant scaling fields, analytic changes from physical controls to $(t,h)$, short-range isotropic scaling, no dangerously irrelevant variable in the thermodynamic scaling function, and no marginal or resonant logarithm. Let ${\lambda }_t,{\lambda }_h>0$ denote their scaling eigenvalues. The homogeneous singular free-energy density obeysChoose a stopping scale $L=b^p\simeq |t{|}^{-1/{\lambda }_t}$ so the renormalized thermal field is order one. Irrelevant fields vanish at this scale. With $s=D/{\lambda }_t$ and $\Delta ={\lambda }_h/{\lambda }_t$,Continuous scale interpolation absorbs the immaterial integer stopping-scale convention. Nonuniversal thermal, field and free-energy metric factors can be included explicitly; the normalized functions depend only on the universality class. The $+$ and $-$ labels refer to $t>0$ and $t<0$, which flow to different sides of the critical surface.

For clarity about the additive terms, let $a(t,h)$ solve $a(t,h)=b^{-D}a(b^{{\lambda }_t}t,b^{{\lambda }_h}h)+g(t,h)$. If $g=\sum g_{mn}{t}^m{h}^n$, its analytic coefficients areprovided the denominators do not vanish. Subtracting $a$ makes the recursion homogeneous. A vanishing denominator is a resonance and generically produces logarithms, which the simple power-law hypothesis excludes. In a convention where the regular field-dependent background has been subtracted, write ${\Phi }_{\pm }=f_{\pm }+C_{\pm }$ and $I_{\pm }(t)=C_0+a(t,0)$. This yields the requested formThe $C_{\pm }$ contain branch-dependent matching contributions; only $f_{\pm }+C_{\pm }$ is physically determined. The $I_{\pm }$ are nonsingular functions analytic through the thermodynamic critical point and represent the same regular Taylor background on the two sides, though either branch may be used to describe them. For the full free energy at nonzero field, the regular background is generally $I(t,h)$, not just $I(t)$. Adding an analytic identity-operator contribution proportional to $h^2$ provides a direct counterexample to an exact unrestricted $I(t)$ formula; it changes no singular exponents. The printed expression is justified for the leading singular scaling form with that analytic field background removed. Confluent corrections from irrelevant fields are also omitted in this leading form.

Differentiate the singular scaling form, with one-sided derivatives in the ordered phase. It givesThe length transforms as $\xi (t,h)=L\xi (L^{{\lambda }_t}t,L^{{\lambda }_h}h)$, so $\nu =1/{\lambda }_t$. Thus the Rushbrooke scaling relation and hyperscaling relation follow by direct substitution:The same reasoning also gives $\gamma =\beta (\delta -1)$. Hyperscaling relation relies on the assumptions above and can fail above an upper critical dimension because of a stabilizing dangerously irrelevant coupling.

For the universal specific-heat amplitude ratio, write $F_s=A_f|a_{t}t{|}^s{\Phi }_{\pm }(0)$ at $h=0$. Two temperature derivatives multiply both branches by the same dimensional and thermal metric factors. Away from logarithmic resonances and exceptional zero amplitudes, $A_{+}/A_{-}={\Phi }_{+}(0)/{\Phi }_{-}(0)$ with the same heat-capacity sign convention on both sides. HenceThe functions and their matched branch constants are tied to the fixed-point theory; the constants are not independent arbitrary experimental parameters. If $\alpha <0$, this means the amplitudes of the subleading singular contribution after subtracting the analytic specific-heat background. A finite background ratio itself need not be universal.

Take a positive gradient coefficient ${\kappa }^{-1}$, a positive squared mass, a finite-volume regulator and an ultraviolet cutoff. Use the dimensionless statistical action convention $Z=\int \mathcal{D}\varphi e^{-H}$: inverse temperature is absorbed into the coefficients. If $H$ denotes dimensional energy instead, the exponent below gains a factor ${\beta }_{\mathrm{th}}$ and the covariance gains ${\beta }_{\mathrm{th}}^{-1}$. These conventions cannot be mixed.

Choose the Fourier transform conventionIntegration of the quadratic gradient and source terms givesComplete the square with ${\tilde{\varphi }}_{\mathrm{cl}}(p)=\tilde{\Delta }(p{)}^{-1}\tilde{J}(p)$. Translation of the regulated real Gaussian variables leaves their measure unchanged, so the Gaussian integral givesHere $Z_G=Z_0[0]$ is the source-free Gaussian partition function. For $n$ independent regulated real coordinates and a positive quadratic matrix $A$, it is $(2\pi {)}^{n/2}(\det A{)}^{-1/2}$ under Lebesgue normalization. The continuum determinant needs the stated regulator and a measure convention. Reality pairs the nonzero Fourier modes, so they must not be counted as two independent real fields.

A connected correlation function subtracts the product of the means. The logarithm $W[J]=\log Z_0[J]$ generates its cumulants: differentiating $W$ twice with respect to the source, using the measure $d^{D}p/(2\pi {)}^D$ in the functional derivative convention, yieldsThe source changes the mean but not the covariance; the stated zero-source limit therefore follows immediately. Momentum conservation uses $p+q$, appropriate for a real field without complex conjugation on either displayed factor.

For the continuum massive Gaussian field correlation tail, set $\mu =\sqrt{\kappa }m>0$ and $r=|x|$. Fourier inversion givesThe large-argument Modified Bessel function of the second kind has $K_{\rho }(z)\sim \sqrt{\pi /(2z)}{e}^{-z}$. Hence $G_0(r)$ is proportional to $r^{-(D-1)/2}{e}^{-\mu r}$ at large distance andAt zero mass the exponential length diverges. For $D>2$ the critical continuum tail is algebraic, $G_0(r)\propto r^{2-D}$; in lower dimensions infrared regulation needs additional care.

For a uniform source $h$, perform a momentum-shell renormalization group step: separate slow modes $|p|<\Lambda /b$ and fast modes $\Lambda /b<|p|<\Lambda$. In a Gaussian field theory they decouple; the uniform source couples only to the retained zero mode. Integrating fast modes contributes only a field-independent determinant to the free energy. Restore the cutoff by $x^{\prime }=x/b$, $p^{\prime }=bp$ andThe factors from the measure, gradient and two fields cancel in the gradient term. The mass term gains $b^2$, and the linear source gains $b^{(D+2)/2}$. ThereforeHere $m$ is the positive mass; for signed thermal coordinates it is the squared-mass relation that applies. The Gaussian fixed point has thermal scaling eigenvalue ${\lambda }_t=2$, source eigenvalue ${\lambda }_h=(D+2)/2$ and zero anomalous dimension.

Under the formal homogeneous Gaussian scaling hypothesis, substitution of those eigenvalues givesThe last two values are independently visible in $\xi \propto t^{-1/2}$ and $\chi ={\tilde{G}}_0(0)=1/m^2\propto t^{-1}$, agreeing with Landau-Ginzburg theory. The specific-heat power follows directly from the determinant: after two thermal derivatives its singular contribution is proportional toThe change of variables $p=\sqrt{\kappa t}q$ proves this power because the rescaled integral converges at large $q$ exactly when $D<4$. The two-dimensional free energy itself has a $t\log t$ resonance, while its second derivative still has the indicated $t^{-1}$ power.

The printed magnetization exponent is a formal Gaussian scaling index, not a proved ordered-phase law of the purely quadratic model. The ordered-phase obstruction in a Gaussian scalar model is explicit: the uniform field has energy density $m^2{M}^2/2-hM$. At $h=0,m^2<0$ it tends to $-\infty$ as $|M|\to \infty$, and the partition function diverges. For $m^2>0$, the mean is $M=h/m^2$ and tends to zero as $h\to 0$, so there is no spontaneously ordered branch. Thus no literal $M\sim |t{|}^{\beta }$ for $t<0$ can be established from this Statistical Hamiltonian alone. The value $(D-2)/4$ follows from the scaling dimension of the field if one assumes such a branch; for $D\le 2$ its nonpositive value further warns against that interpretation. Adding a stabilizing quartic interaction defines an ordered phase, but below four dimensions that interaction is relevant and generally leads to the Wilson-Fisher fixed point, not these Gaussian thermodynamic powers.

For the Gaussian specific-heat infrared threshold, at $D=4$ the same determinant integral grows like $\log ({\Lambda }^2/t)$, so $\alpha =0$ as a power exponent, with a logarithmic divergence. For $D>4$, it has a finite cutoff-dependent limit; the leading heat capacity is a regular constant and has no divergent power, again described by $\alpha =0$ in the requested convention. After subtraction of this background the Gaussian singular term can still have the negative power index $2-D/2$; it is not correct to erase that distinction. In a stable quartic Landau-Ginzburg theory above four dimensions, the quartic coefficient is a dangerously irrelevant coupling: it is needed for the ordered phase and contributes $f_s\propto -t^2/u$, giving the mean-field specific-heat jump and $\alpha =0$ while invalidating naive hyperscaling relation. At four dimensions, marginal interactions can add further logarithmic corrections. The agreement with Landau theory is in the leading heat-capacity power; it does not assert the absence of logarithms or supply a missing stable Gaussian ordered phase.

Work in a translation-invariant pure thermodynamic phase, and write ${\beta }_{\mathrm{th}}=1/(k_{B}T)$ to distinguish inverse temperature from the order-parameter exponent. The connected correlation function isAway from criticality in the massive scalar order channel, it decays exponentially, possibly multiplied by an algebraic factor. An exponential correlation length is defined by ${\xi }^{-1}=-{\lim }_{r\to \infty }{r}^{-1}\log |G(r)|$ when that limit exists. At a continuous transition $\xi$ diverges and the critical decay is algebraic. In an ordered symmetric mixture, a nondecaying contribution can remain; selecting a pure branch prevents confusing it with connected critical fluctuations.

Let $M={\sum }_{\mathbf{n}}{\sigma }_{\mathbf{n}}$ and $m=⟨M⟩/N$. Differentiating the partition function with the energy term $-hM$ gives the correlation-function susceptibility sum ruleIf the source is instead the dimensionless ${\beta }_{\mathrm{th}}h$, the explicit inverse-temperature factor is absent. In continuum physical coordinates the lattice sum is approximately $a^{-D}\int d^{D}rG(\mathbf{r})$.

Choose a normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )\ge 0$ such that ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$, replacing the sum by an integral for continuous variables. A deterministic example is a product of delta functions imposing that each coarse variable equals the average of spins in a block of $b^D$ sites, with its field normalization included consistently. Define the effective Boltzmann weight byKernel normalization makes the blocked partition function exactly equal to the original. Its coarse-grained variables retain the long-distance observables through their defining relation to the original fields. An exact RG step generally generates all symmetry-allowed operators; keeping only a small coupling set is an approximation, not an exact closure assumption. The geometrical coarse lattice hasso its physical volume is unchanged. A subsequent coordinate rescaling may restore the numerical lattice spacing to its initial value, which is the equivalent rescaled-coordinate convention.

Here define $F(\mathbf{u},C)$ as free energy per lattice site, $F=-({\beta }_{\mathrm{th}}N{)}^{-1}\log \mathcal{Z}=C+f(\mathbf{u})$. Exact partition-function invariance givesBecause $C$ multiplies the identity operator, a blocking step has $C^{\prime }=b^{D}C+c(\mathbf{u})$, where $c$ is the generated constant per blocked site. With $g(\mathbf{u})=b^{-D}c(\mathbf{u})$ this impliesThe identity-operator contribution to renormalization-group free energy includes the eliminated modes' entropy, connected vacuum terms and field-measure normalization. It cannot be discarded when calculating an absolute free energy, even though it cancels out of normalized correlation functions.

There is a useful convention behind the word “singular” in this inhomogeneous equation. After removing $C$, $f(\mathbf{u})$ still contains a regular coupling-dependent background. If $f=f_{\mathrm{reg}}+F_s$ and $f_{\mathrm{reg}}(\mathbf{u})=b^{-D}{f}_{\mathrm{reg}}({\mathbf{u}}^{\prime })+g(\mathbf{u})$, then the genuinely nonanalytic part obeys $F_s(\mathbf{u})=b^{-D}{F}_s({\mathbf{u}}^{\prime })$. Thus the printed inhomogeneous representative and homogeneous singular scaling below are consistent after background subtraction. At resonances or marginal points this subtraction may leave additive or multiplicative logarithms; a strictly homogeneous pure-power form is then qualified accordingly.

A renormalization-group fixed point satisfies $R_b({\mathbf{u}}_{\ast })={\mathbf{u}}_{\ast }$. Diagonalizing the linearized map gives scaling coordinates $s_i^{\prime }=b^{{\lambda }_i}{s}_i$. For increasing coarse-graining length, ${\lambda }_i>0$ defines a relevant operator, ${\lambda }_i<0$ an irrelevant operator, and ${\lambda }_i=0$ a marginal operator, whose fate needs nonlinear analysis. Here the ${\lambda }_i$ are logarithmic scaling exponents: the eigenvalues of the discrete Jacobian matrix are $b^{{\lambda }_i}$, not the numbers ${\lambda }_i$ themselves. The critical surface is the stable manifold flowing into the critical fixed point after every relevant scaling field is tuned to zero. A repulsive renormalization-group trajectory leaves the fixed point along relevant directions and flows toward a different long-distance phase.

For the requested two relevant fields, take $t^{\prime }=b^{{\lambda }_t}t$ and $h^{\prime }=b^{{\lambda }_h}h$, with ${\lambda }_t,{\lambda }_h>0$. Ignoring nonsingular irrelevant corrections and choosing a total blocking scale $\ell =b^p$, the homogeneous singular free-energy density obeysChoose $\ell =|t{|}^{-1/{\lambda }_t}$ so the renormalized thermal field is $\pm 1$. This provesMetric factors and regular analytic redefinitions of the scaling fields only change amplitudes. Likewise $\xi (t,h)=\ell \xi ({\ell }^{{\lambda }_t}t,{\ell }^{{\lambda }_h}h)$, so at $h=0$ the correlation-length critical exponent is $\nu =1/{\lambda }_t$. Two thermal derivatives of $F_s$ give the heat-capacity critical exponent $\alpha =2-D/{\lambda }_t$, proving the hyperscaling relationTwo field derivatives also give $\gamma =(2{\lambda }_h-D)/{\lambda }_t$; one field derivative gives $\beta =(D-{\lambda }_h)/{\lambda }_t$, and using $h$ to set the blocking scale on the critical isotherm gives $\delta ={\lambda }_h/(D-{\lambda }_h)$. These statements assume that no dangerously irrelevant coupling makes the scaling function singular as it is removed. In particular the naive hyperscaling form need not describe the quartic ordered phase above four dimensions, where its stabilizing quartic coupling is dangerously irrelevant. That restriction reconciles this derivation with the mean-field exponents and Question 3.

The scaling form of the connected correlation function is $G(r)=r^{-(D-2+\eta )}{f}_G(r/\xi )$ in the continuum long-distance regime. For $r\ll \xi$, $f_G$ tends to a finite nonzero constant. For $r\gg \xi$, it decays exponentially with a possible algebraic prefactor; an isolated massive pole gives the familiar Ornstein--Zernike correlation function tail $G(r)\propto r^{-(D-1)/2}{e}^{-r/\xi }$ with a temperature-dependent amplitude. Exponential decay is the essential feature; the same critical algebraic power need not remain the exact large-distance prefactor.

Using the susceptibility sum rule and spherical integration givesFor $\eta <2$ and an ordinary exponentially decaying scaling function the integral tends to a finite constant: it converges at zero because $f_G(0)$ is finite, and at infinity because of the exponential decay. Microscopic distances add a regular background. Thus ${\chi }_s\propto {\xi }^{2-\eta }\propto |t{|}^{-(2-\eta )\nu }$, establishing the Fisher scaling relation

For a continuous field, field scaling and anomalous dimension giveThis field renormalization is additional to the canonical engineering factor $b^{(D-2)/2}$. It makes $G(r)=b^{-2x_{\sigma }}{G}^{\prime }(r/b)$. Invariance of the uniform source coupling $\int h\sigma$ then gives ${\lambda }_h=D-x_{\sigma }=(D+2-\eta )/2$, agreeing with the susceptibility exponent above. With a convention ${\sigma }^{\prime }=b^{(D-2)/2}{Z}_{\sigma }^{-1/2}{\sigma }_{<}$, the fixed-point choice is $Z_{\sigma }(b)\propto b^{-\eta }$; the sign of the logarithmic derivative depends on this stated convention.

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

If the singular free-energy density is set by one correlated region, then $f_s\sim {\xi }^{-d}$. With reduced temperature $t$, the definitions $\xi \sim |t{|}^{-\nu }$ and $f_s\sim |t{|}^{2-\alpha }$ therefore giveThis is the hyperscaling relation; it assumes that no dangerously irrelevant coupling introduces an additional scale.