Critical Propagators and the Upper Critical Dimension

The previous page derived the continuum scalar field from the Ising model by an exact Hubbard–Stratonovich transformation. The Gaussian part of that field theory already knows a lot: it knows which momentum mode becomes soft, how the correlation length diverges, and why the two-point function near criticality has the universal form of a massive free scalar propagator.

This page asks the next question. Once the continuum action contains an interaction such as $\lambda \phi^4$, when is the Gaussian propagator trustworthy? The answer is not simply “when $\lambda$ is small.” Near a critical point the correlation length $\xi$ becomes large, so fluctuations are sampled over larger and larger regions. The true expansion parameter is a scale-dependent one,

$$g(\xi)\sim \lambda \xi^{4-D}.$$

Thus the Ising $\phi^4$ interaction becomes less important at long distances for $D>4$, marginal for $D=4$, and more important for $D<4$. The number

$$D_c=4$$

is the upper critical dimension. For the ordinary short-range scalar transition, bulk mean-field exponents hold above it. Below it, the quartic interaction destabilizes the Gaussian fixed point. Near four dimensions the Wilson–Fisher fixed point provides a controlled replacement. Power counting by itself does not establish that a finite-temperature transition exists in every dimension.

Required background. Hubbard–Stratonovich transformation and the continuum field supplies the lattice quadratic kernel, the auxiliary-field propagator, and the continuum $\phi^4$ action used below.

Critical modes of a lattice kernel

Propagator normalization

Throughout this page $D$ denotes the Euclidean spatial dimension of the statistical system. We use

$$\int_k \equiv \int {d^D k\over (2\pi)^D},$$

usually with an ultraviolet cutoff $\Lambda\sim a^{-1}$ inherited from the lattice spacing $a$.

The continuum Gaussian action is normalized as

$$S_0[\phi]= {1\over2}\int_k (k^2+m^2)\phi(k)\phi(-k),$$

Here and below a momentum-space propagator is the coefficient after removing the momentum-conserving delta function:

$$\langle\phi(k)\phi(k')\rangle_0 =(2\pi)^D\delta^{(D)}(k+k')G_0(k), \qquad G_0(k)={1\over k^2+m^2}.$$

In the interacting symmetric phase, use the bare mass $m_0$ in the reference propagator and define the self-energy by

$$G(k)^{-1}=G_{0,\mathrm{bare}}(k)^{-1}+\Sigma(k) =k^2+m_0^2+\Sigma(k).$$

Some QFT books define the self-energy with the opposite sign. Nothing physical depends on this convention; what matters is the exact inverse two-point vertex $\Gamma^{(2)}(k)=G(k)^{-1}$.

For an interacting theory, $\Gamma^{(2)}(0)=\chi^{-1}$ is exactly the inverse susceptibility, not automatically the square of the exponential correlation-length mass. If

$$\Gamma^{(2)}(k)=\chi^{-1}+Z_2 k^2+O(k^4),$$

then the second-moment correlation length is $\xi_2^2=Z_2\chi$. After rescaling the field so that $Z_2=1$, it is convenient to write $m^2=\chi^{-1}=\xi_2^{-2}$, with $\chi$ now referring to that normalized field. The exponential correlation length has the same critical exponent in the scaling regime considered here, but its amplitude need not equal that of $\xi_2$, even in the scaling limit.

Start with a translation-invariant ferromagnetic interaction. In the spin language one may write schematically

$$H=-{1\over2}\sum_{x,y}J(x-y)\sigma_x\sigma_y,$$

with real even couplings $J(r)=J(-r)$, nonnegative off-site couplings, and finite range or exponential decay. These hypotheses make $V(k)$ analytic and ensure $V(k)\leq V(0)$. Assume also a nondegenerate isotropic quadratic maximum; an anisotropic positive stiffness would give a tensor of correlation lengths. Its Fourier transform is

$$V(k)=\sum_r J(r)e^{ik\cdot r}.$$

The Hubbard–Stratonovich Gaussian propagator from the previous page has the form

$$G_{\mathrm{lat}}(k)\propto {1\over 1-\beta V(k)},$$

up to a smooth nonzero normalization factor. The important information is the zero of the inverse propagator. In the Gaussian or random-phase approximation, an instability occurs when some eigenvalue of the interaction kernel reaches

$$1-\beta_c^{\mathrm G} V(k_\star)=0.$$

For a ferromagnet, the largest value of $V(k)$ is at $k_\star=0$, so

$$\boxed{\beta_c^{\mathrm G} V(0)=1.}$$

If the maximum is instead at a nonzero wavevector $Q$, the soft field is not the uniform magnetization. The continuum expansion must be made around $Q$, which describes antiferromagnetic or modulated order. This is the same mathematical mechanism with a different ordering wavevector.

Under the stated analyticity and isotropy assumptions,

$$V(k)=V(0)-\rho k^2+O(k^4), \qquad \rho>0.$$

Define the reduced distance from the Gaussian critical point by

$$\tau=1-\beta V(0).$$

On the disordered side, $\tau>0$. Then

$$1-\beta V(k) = \tau+\beta\rho k^2+O(k^4).$$

Write the smooth numerator to leading order as $A>0$ and set $c=\beta\rho>0$. Dividing the whole denominator by $c$ gives the Ornstein–Zernike form

$$\boxed{ G_{\mathrm{lat}}(k)\simeq {A\over\tau+c k^2} ={Z\over k^2+m^2}, \qquad Z={A\over c},\quad m^2={\tau\over c}. }$$

Thus $m^2\propto T-T_c^{\mathrm G}$ on the disordered side. For this quadratic propagator, $m$ is the inverse second-moment length and also the continuum pole mass:

$$\xi_2=m^{-1}.$$

Therefore the Gaussian theory predicts

$$\xi\sim |T-T_c^{\mathrm G}|^{-1/2}.$$

This is the mean-field value $\nu_{\mathrm{MF}}=1/2$.

The figure separates the stiffness $c$, which sets the length, from the field amplitude $A$.

The soft mode lies at the maximum of the lattice kernel. On the disordered side of an isotropic ferromagnetic Gaussian instability, $c=\beta\rho>0$, $m^2=\tau/c$, and $Z=A/c$. The quadratic propagator has $\xi_2=1/m$; a different ordering wavevector requires expansion around that wavevector. The kernel curve is schematic, not a quantitative lattice spectrum.

For the nearest-neighbor hypercubic Ising model,

$$V(k)=2J\sum_{\mu=1}^D\cos(k_\mu a).$$

At small $k$,

$$V(k)=2DJ-Ja^2k^2+O(k^4a^4).$$

Thus

$$1-\beta V(k) = (1-2D\beta J)+\beta J a^2 k^2+\cdots.$$

The Gaussian critical value is

$$\beta_c^{\mathrm{G}}J={1\over 2D},$$

and

$$G_\sigma^{(0)}(k)\simeq {1\over (1-2D\beta J)+\beta Ja^2 k^2}.$$

This last expression is the dimensionless spin susceptibility of the Gaussian approximation, before canonical field normalization. The exact critical point is shifted by fluctuations in dimensions where the interaction is important. The Gaussian calculation nevertheless identifies the soft mode and the analytic structure around it. The superscript on $\beta_c^{\mathrm G}$ is essential; $1/(2D)$ is not the exact nearest-neighbor Ising critical coupling.

Position-space propagator and the correlation length

For the canonically normalized field ($Z=1$), the continuum Gaussian propagator at a noncoincident point is

$$G_0(x)=\int_k {e^{ik\cdot x}\over k^2+m^2}.$$

The formulas in this section use the continuum Green distribution, with the ultraviolet regulator removed at $x\ne0$. An abrupt finite momentum cutoff adds oscillatory long-distance tails, so one must not extract an exponential length from that cutoff artifact. Loop integrals later retain their stated cutoff. The distinction between the pole length, the second-moment length and sharp-cutoff ringing is explicit in Wilson and Kogut 1974, § 3, pp. 98–99.

At criticality, $m=0$, dimensional analysis alone gives

$$G_0(x)\propto {1\over |x|^{D-2}},$$

for $D>2$. With the normalization above, the exact coefficient is

$$\boxed{ G_0(x) = {\Gamma\!\left({D\over2}-1\right)\over 4\pi^{D/2}} {1\over |x|^{D-2}}, \qquad m=0,\quad D>2. }$$

Away from criticality,

$$G_0(x) = {1\over (2\pi)^{D/2}} \left({m\over |x|}\right)^{D/2-1} K_{D/2-1}(m|x|),$$

where $K_\nu$ is a modified Bessel function and $m>0$. This formula has two important limits. For $D>2$ and distances much shorter than $\xi=1/m$,

$$|x|\ll \xi, \qquad G_0(x)\sim {1\over |x|^{D-2}},$$

so the system looks critical. In $D=2$ the corresponding behavior is logarithmic:

$$G_m(R)={1\over2\pi}K_0(mR) \simeq-{\log(mR/2)+\gamma_E\over2\pi}, \qquad mR\ll1.$$

The massless field has no finite translation-invariant covariance of this form in two dimensions, but the difference $G_m(R)-G_m(R_0)$ tends to $-\log(R/R_0)/(2\pi)$ as $m\to0$. For distances much longer than $\xi$,

$$|x|\gg \xi, \qquad G_0(x)\sim {e^{-|x|/\xi}\over |x|^{(D-1)/2}},$$

so correlations are exponentially suppressed.

The susceptibility is the zero-momentum two-point function,

$$\chi=G_0(k=0)={1\over m^2}.$$

Since $m^2\propto T-T_c^{\mathrm G}$, Gaussian theory gives

$$\chi\sim |T-T_c^{\mathrm G}|^{-1}.$$

Thus $\gamma_{\mathrm{MF}}=1$. At this level the anomalous dimension is also zero, because at criticality

$$G_0(x)\sim {1\over |x|^{D-2}} = {1\over |x|^{D-2+\eta}} \quad\Longrightarrow\quad \eta_{\mathrm{MF}}=0.$$

The central issue is whether interactions preserve these exponents.

Dyson equation and critical mass tuning

The continuum action near the Ising critical point is

$$S[\phi] = \int d^D x\, \left[ {1\over2}(\partial\phi)^2 +{1\over2}m_0^2\phi^2 +{\lambda\over4!}\phi^4 +\cdots \right].$$

The dots include higher even powers and higher derivatives. For the present page the leading interaction is the quartic term. Perturbatively, the exact two-point function is organized by one-particle-irreducible self-energy insertions:

$$G(k)={1\over k^2+m_0^2+\Sigma(k)}.$$

The physical susceptibility is not fixed by the bare parameter $m_0$. At zero momentum the exact inverse propagator gives

$$\boxed{ \chi^{-1}\equiv \Gamma^{(2)}(0)=G(0)^{-1} =m_0^2+\Sigma(0). }$$

The critical point is the value of the microscopic parameters for which $\chi$ diverges,

$$\chi^{-1}=0.$$

Thus criticality is a tuning condition. The bare mass must cancel the fluctuation correction:

$$m_{0,c}^2+\Sigma_c(0)=0.$$

Near the critical point it is often cleaner to subtract the critical value:

$$\chi^{-1} = m_0^2-m_{0,c}^2 + \left[\Sigma(0;m)-\Sigma_c(0)\right].$$

This formula is the statistical-mechanics version of mass renormalization. The lattice cutoff regulates ultraviolet divergences; it does not remove critical infrared singularities. The separation between the physical susceptibility scale and the microscopic bare mass is still essential. When the coefficient of $k^2$ has been normalized to one, we may abbreviate $\chi^{-1}$ as $m^2$ and identify it with $\xi_2^{-2}$, as in the convention note.

Self-energy insertions shift the inverse propagator. The critical point is defined by $\Gamma^{(2)}(0)=0$, not by the naive vanishing of the bare mass. With $G^{-1}=G_0^{-1}+\Sigma$, expanding $G$ produces alternating signs. In $\phi^4$ theory the one-loop tadpole shifts $\Gamma^{(2)}(0)$; momentum-dependent corrections begin at two loops.

Using a propagator with the renormalized scaling mass $m=\xi_2^{-1}$—equivalently, reorganizing perturbation theory around the physical quadratic term—the one-loop tadpole gives

$$\boxed{ \Sigma_1={\lambda\over2}\int^\Lambda {d^D q\over (2\pi)^D} {1\over q^2+m^2}. }$$

The factor $1/2$ follows from the conventional $\lambda\phi^4/4!$ normalization. Other normalizations of the quartic term move this numerical factor but not the scaling.

Write

$$I_D(m)=\int^\Lambda {d^D q\over (2\pi)^D}{1\over q^2+m^2} = {S_{D-1}\over(2\pi)^D} \int_0^\Lambda dq\,{q^{D-1}\over q^2+m^2},$$

where

$$S_{D-1}={2\pi^{D/2}\over \Gamma(D/2)}$$

is the area of the unit $(D-1)$-sphere. The ultraviolet behavior is

$$I_D(0)\sim {S_{D-1}\over(2\pi)^D} {\Lambda^{D-2}\over D-2}, \qquad D>2.$$

This cutoff-dependent constant shifts the critical temperature. It is not universal. The universal question is how the remaining long-distance part behaves as $m\to0$.

For $2<D<4$,

$$I_D(m)-I_D(0)\sim -C_D m^{D-2}, \qquad C_D>0.$$

For $D=4$,

$$I_4(m)-I_4(0) \sim -{m^2\over 8\pi^2}\log{\Lambda\over m} +\text{analytic terms}.$$

For $D>4$, the leading small-$m$ correction is analytic in $m^2$ once the cutoff is kept fixed:

$$I_D(m)-I_D(0)\sim -A_D \Lambda^{D-4}m^2+\cdots.$$

For noninteger $D>4$, a subleading nonanalytic term proportional to $m^{D-2}$ can also occur; it does not overturn the leading comparison with $m^2$. At even dimensions further logarithms appear in subleading orders.

The borderline behavior at $D=4$ is the first warning that four dimensions are special. The next section gives a cleaner and more general explanation.

Power counting at the Gaussian fixed point

At criticality the Gaussian action is

$$S_0={1\over2}\int d^D x\,(\partial\phi)^2.$$

Under the coordinate dilation

$$x'=b x,$$

the derivative scales as

$$\partial\mapsto b^{-1}\partial,$$

and the measure scales as

$$d^D x\mapsto b^D d^D x.$$

The transformed field that leaves $S_0$ invariant is

$$\phi'(x')=b^{-(D-2)/2}\phi(x).$$

Equivalently, the engineering dimension of the scalar field is

$$\boxed{ [\phi]={D-2\over2}. }$$

Now consider the quartic interaction

$$S_{\mathrm{int}}={\lambda\over4!}\int d^D x\,\phi^4.$$

The operator $\phi^4$ has dimension

$$[\phi^4]=2(D-2).$$

Since the action is dimensionless,

$$[\lambda]-D+2(D-2)=0.$$

Therefore

$$\boxed{ [\lambda]=4-D. }$$

This single equation encodes the upper critical dimension:

$$\begin{array}{ccl} D>4 &:& \lambda \text{ has negative dimension, so it is irrelevant at long distance},\\ D=4 &:& \lambda \text{ is classically marginal},\\ D<4 &:& \lambda \text{ has positive dimension, so it is relevant at long distance}. \end{array}$$

At a length scale $R$, the associated dimensionless coupling is

$$\boxed{ g(R)\sim \lambda R^{4-D}. }$$

Thus $g(R)\to0$ as $R\to\infty$ for $D>4$. This is the precise sense in which mean-field theory becomes exact at the longest distances above four dimensions. For $D<4$, $g(R)$ grows with $R$, so the Gaussian approximation eventually breaks down no matter how small the microscopic coupling was.

At $D=4$, power counting alone cannot decide the fate of $\lambda$. Quantum or statistical fluctuations generate logarithms, and the next page will turn those logarithms into a renormalization-group flow.

Position-space Ginzburg test

There is a useful way to see the same result in a correlation function. Initially keep a UV regulator of length $a$, a positive mass as an IR regulator, and, if a finite box is used, its length $L$. Work at separations $a\ll R\ll\xi,L$. A finite periodic box alone would leave the massless zero mode unregulated. For $D>2$, the critical Gaussian power law in that range is

$$G_0(R)\sim {1\over R^{D-2}}.$$

The composite operator $\phi^2$ is the continuum representative of the local energy-density perturbation. After subtracting its expectation value, its Gaussian connected two-point function behaves as

$$C_0(R) \equiv \langle \phi^2(x_1)\phi^2(x_2)\rangle_{0,c} =2G_0(R)^2 \sim {1\over R^{2D-4}}, \qquad R=|x_1-x_2|.$$

The subscript $c$ is essential: the raw expectation also contains $\langle\phi^2(x_1)\rangle\langle\phi^2(x_2)\rangle$. At a regulator where coincident fields exist, Wick’s theorem gives that disconnected term plus the two cross contractions displayed above.

Consider the first-order graph in which all four fields at one interaction vertex contract with the two composite insertions. Expanding $e^{-S_{\mathrm{int}}}$ gives a minus sign. There are $\binom42 2!2!=24$ contractions, cancelling $4!$, so this particular graph is

$$\delta C^{\mathrm{ins}}(R) =-\lambda \int d^D z\, G_0(x_1-z)^2G_0(x_2-z)^2.$$

The propagators and integral here first retain the regulators. Tadpole graphs, the adjustment of the mass and renormalization of the composite operator are separate terms in the full connected correlator.

To isolate a contribution at the observation scale, restrict the vertex to

$$\mathcal A_R=\left\{z: {R\over3}<|z-x_i|<2R\ \text{for }i=1,2\right\}.$$

This region stays away from both endpoints and from infinity; its shape is fixed under $z=x_1+Ry$. Its volume scales as

$$\int_{\mathcal A_R}d^D z\sim R^D,$$

and each of the four propagators is of order $R^{-(D-2)}$. The magnitude of the restricted graph therefore scales as

$$|\delta C_R^{\mathrm{ins}}| \sim \lambda {R^D\over R^{4D-8}} = {\lambda\over R^{3D-8}}.$$

Here $\lambda\geq0$, as appropriate to the stable quartic potential. Its relative magnitude is

$$\boxed{ {|\delta C_R^{\mathrm{ins}}|\over C_0(R)} \sim \lambda R^{4-D}. }$$

This is the same dimensionless coupling $g(R)$ found by power counting. It is a scale-local estimate, not a claim that the full massless insertion is finite in every dimension.

Indeed, removing both regulators in that full graph imposes two different tests. Near either endpoint, a radial distance $s$ gives an integral $\int_0 ds\,s^{3-D}$, finite only for $D<4$. Far from both endpoints, the radial integral is $\int^\infty ds\,s^{7-3D}$, finite only for $D>8/3$. Thus the unregulated graph converges precisely in

$${8\over3}<D<4.$$

At the endpoints the corresponding divergence is logarithmic. UV counterterms do not remove the far-distance divergence at $D\leq8/3$; keep the mass or volume regulator there. The engineering argument for $D_c=4$ remains valid outside this convergence window, whereas this unregulated massless graph does not.

A finite three-dimensional check

For $D=3$ and $R>0$, $G_0(R)=1/(4\pi R)$ and the full insertion is convergent. A Feynman parameter evaluates the required convolution:

$$\begin{aligned} \int {d^3z\over |z|^2|z-x|^2} &=\int_0^1 du\int {d^3w\over[w^2+u(1-u)R^2]^2}\\ &={\pi^2\over R}\int_0^1 {du\over\sqrt{u(1-u)}} ={\pi^3\over R},\qquad R=|x|. \end{aligned}$$

The shift is $w=z-(1-u)x$. Hence this graph, with the stated $\lambda/4!$ normalization, gives

$$C_0(R)={1\over8\pi^2R^2},\qquad \delta C^{\mathrm{ins}}(R)=-{\lambda\over256\pi R},\qquad {\delta C^{\mathrm{ins}}(R)\over C_0(R)}=-{\pi\lambda R\over32}.$$

This checks the sign, numerical combinatorics and the power $\lambda R$ for this graph. It does not replace the other first-order contributions to the full renormalized composite correlator. The figure summarizes the regional estimate and its separate convergence qualification.

For $D>2$ in the Gaussian scaling range, the connected free composite correlator is $C_0=2G_0^2$. Restricting the quartic vertex to $\mathcal A_R$ gives $|\delta C_R^{\mathrm{ins}}|/C_0\sim\lambda R^{4-D}$. The full massless graph instead requires $8/3<D<4$; outside that window retain the appropriate regulators and renormalization terms. Lines and vertex positions are schematic.

The same estimate is often called the Ginzburg criterion. Evaluate the dimensionless coupling at the correlation length:

$$g(\xi)\sim \lambda \xi^{4-D}.$$

As $T\to T_c$, $\xi\to\infty$. Therefore:

• for $D>4$, $g(\xi)\to0$, so the Gaussian or mean-field approximation becomes asymptotically reliable;
• for $D=4$, $g(\xi)$ is marginal and receives logarithmic corrections;
• for $D<4$, $g(\xi)\to\infty$, so critical fluctuations invalidate mean-field exponents.

These statements concern the growth of the Gaussian expansion parameter where a continuous transition exists. They do not derive the interacting exponents or establish a transition in every dimension. For small positive $\epsilon$ in $D=4-\epsilon$, the Wilson–Fisher fixed point is a controlled replacement; see Wilson and Fisher 1972, pp. 241–243.

The random-walk meaning of 4 equals 2 plus 2

There is a beautiful interpretation of the number four. The Gaussian propagator has the Schwinger representation

$${1\over k^2+m^2} = \int_0^\infty dL\,e^{-L(k^2+m^2)}.$$

Fourier transforming gives

$$G_0(x) = \int_0^\infty dL\, {1\over(4\pi L)^{D/2}} \exp\left(-{x^2\over4L}-m^2L\right).$$

The parameter $L$ behaves like the length, or proper time, of a Brownian path. At criticality $m=0$, the exponential factor

$$\exp\left(-{R^2\over4L}\right)$$

says that paths contributing to a displacement $R$ typically have

$$L\sim R^2.$$

Thus a free critical propagator is geometrically like a random walk of fractal dimension $2$.

The quartic interaction couples two such walks when they meet. The elementary intersection-dimension estimate says that two sets of fractal dimension $2$ have a positive-dimensional intersection when

$$2+2>D.$$

The borderline is

$$D=2+2=4.$$

This is the geometric version of the upper critical dimension. Above four dimensions, intersections of two long independent walks are sufficiently sparse that the local interaction becomes irrelevant; below four dimensions they proliferate and reshape the long-distance theory. At the borderline, the estimate is marginal and logarithms decide the result. This is an interpretation of the power counting, not a substitute for the renormalization-group argument.

This random-walk viewpoint will reappear later in the course when field-theoretic propagators are related to sums over paths, and eventually when random surfaces enter the discussion.

Mean field above four dimensions and its caveat

For $D>D_c$, the long-distance critical two-point function is Gaussian:

$$G(k)\sim {1\over k^2+m^2}.$$

The critical exponents take their mean-field values:

$$\nu={1\over2}, \qquad \eta=0, \qquad \gamma=1.$$

However, one should be careful with the phrase “the interaction is irrelevant.” The quartic coupling is irrelevant at the Gaussian fixed point in the renormalization-group sense, but it is still needed to stabilize the ordered phase. If $m^2<0$, the potential

$${1\over2}m^2\phi^2$$

is unbounded unless the quartic term is kept. Thus $\lambda$ is dangerously irrelevant above four dimensions: it does not change the leading long-distance two-point critical behavior, but it enters thermodynamic quantities such as the magnetization amplitude in the ordered phase.

This subtlety is one reason critical phenomena are more delicate than a first reading of power counting suggests. Power counting identifies which fixed point is stable; it does not automatically tell us which amplitudes and scaling relations remain ordinary. In particular, bulk mean-field exponents hold for $D>4$, but hyperscaling and naive finite-size scaling fail because they are sensitive to the dangerously irrelevant coupling.

Example: the one-loop mass shift

Let us make the mass-renormalization statement concrete. In the symmetric phase, use

$$S[\phi]= \int d^D x\, \left[ {1\over2}(\partial\phi)^2 +{1\over2}m_0^2\phi^2 +{\lambda\over4!}\phi^4 \right].$$

At one loop,

$$m^2=m_0^2+{\lambda\over2}I_D(m)+O(\lambda^2).$$

Using $m$ rather than $m_0$ in the internal line is a convenient self-consistent notation; replacing it by the zeroth-order mass changes the result only beyond the displayed order away from the critical infrared singularity.

The critical bare mass is obtained by setting $m=0$:

$$m_{0,c}^2=-{\lambda\over2}I_D(0)+O(\lambda^2).$$

Therefore

$$m^2 = (m_0^2-m_{0,c}^2) + {\lambda\over2}\left[I_D(m)-I_D(0)\right] +O(\lambda^2).$$

The first term is the experimentally or microscopically tunable distance from the critical point. In a thermal system it is proportional to $T-T_c$ after nonuniversal normalization. The second term is the fluctuation correction.

For $D>4$, the bracket is proportional to $m^2$ times a cutoff-dependent coefficient, so it can be absorbed into a finite renormalization of the slope relating $m^2$ to $T-T_c$. Mean-field scaling survives.

For $D=4$, the bracket contains

$$-m^2\log{\Lambda\over m},$$

so logarithmic corrections appear.

For $2<D<4$, the bracket contains a nonanalytic contribution

$$-m^{D-2}.$$

This term becomes more singular than $m^2$ as $m\to0$. At and below $D=2$, even $I_D(0)$ has an infrared divergence, so the subtraction must be reformulated. In either case ordinary perturbation theory around the Gaussian point has lost control in the infrared. Near four dimensions, the $\epsilon$ expansion repairs this by moving the expansion point from the Gaussian fixed point to a nearby interacting fixed point.

Summary

The lattice interaction kernel determines which mode becomes critical. For a ferromagnet,

$$\beta_c^{\mathrm G} V(0)=1, \qquad V(k)=V(0)-\rho k^2+\cdots,$$

and the Gaussian propagator near the critical point takes the universal form

$$G_0(k)\simeq {Z\over k^2+m^2}, \qquad m=\xi^{-1}, \qquad m^2\propto |T-T_c^{\mathrm G}|.$$

At criticality for $D>2$,

$$G_0(x)\sim {1\over |x|^{D-2}},$$

so Gaussian theory predicts $\nu=1/2$, $\eta=0$, and $\gamma=1$.

Interactions shift the critical point through the self-energy:

$$G(k)^{-1}=k^2+m_0^2+\Sigma(k), \qquad \chi^{-1}=G(0)^{-1}.$$

After normalizing the $k^2$ coefficient, $\chi^{-1}=m^2=\xi_2^{-2}$.

The one-loop tadpole gives

$$\Sigma_1={\lambda\over2}\int^\Lambda {d^D q\over(2\pi)^D}{1\over q^2+m^2},$$

so the bare mass must be tuned against fluctuation corrections to reach $m=0$.

The quartic coupling has engineering dimension

$$[\lambda]=4-D.$$

Equivalently, its dimensionless strength at length scale $R$ is

$$g(R)\sim \lambda R^{4-D}.$$

Thus

$$D_c=4$$

is the upper critical dimension of the Ising $\phi^4$ theory. Above four dimensions the Gaussian fixed point controls the critical exponents. Below four dimensions the quartic interaction grows at long distances, and the correct critical theory is interacting.

Common pitfalls

Calling the Gaussian instability exact. The condition $\beta_c^{\mathrm G} V(0)=1$ is not generally the exact critical temperature. It is the point at which the quadratic approximation becomes unstable; fluctuations shift it.

Equating every definition of mass. Exactly, $G(0)^{-1}=\chi^{-1}$. It equals $\xi_2^{-2}$ only after normalizing the $k^2$ coefficient, while the exponential correlation length comes from the nearest complex-momentum singularity. These definitions share a critical exponent but can differ by finite factors.

Discarding the tadpole as a divergent nuisance. In a lattice statistical model it is finite, but it still changes the relation between microscopic temperature and physical correlation length. Renormalization is already present conceptually before taking a continuum cutoff to infinity.

Mistaking power counting for the full RG. Power counting identifies $D=4$ as the marginal dimension. At exactly four dimensions logarithms decide the flow; below four dimensions the Wilson–Fisher fixed point replaces the Gaussian fixed point.

Dropping the quartic term above four dimensions. It is irrelevant for leading bulk critical exponents but still stabilizes the ordered phase and affects amplitudes, hyperscaling, and finite-size scaling. This is the standard dangerously irrelevant caveat.

Exercises

Exercise 1

Let

$$G_{\mathrm{lat}}(k)={A\over 1-\beta V(k)}$$

and suppose

$$V(k)=V(0)-\rho k^2+O(k^4), \qquad \rho>0.$$

Assuming $\beta_c^{\mathrm G}V(0)=1$, show that near $\beta_c^{\mathrm G}$ the propagator can be written in the form

$$G_{\mathrm{lat}}(k)\simeq {Z\over k^2+m^2}.$$

Find $m^2$ in terms of $\beta$, $V(0)$, and $\rho$ up to an overall normalization convention.

Solution

Using the small-$k$ expansion,

$$1-\beta V(k) = 1-\beta V(0)+\beta\rho k^2+O(k^4).$$

Define

$$\tau=1-\beta V(0).$$

Then

$$G_{\mathrm{lat}}(k)\simeq {A\over \tau+\beta\rho k^2}.$$

Factor out $\beta\rho$ from the denominator:

$$G_{\mathrm{lat}}(k) \simeq {A\over \beta\rho} {1\over k^2+\tau/(\beta\rho)}.$$

Thus

$$Z={A\over\beta\rho}, \qquad m^2={\tau\over\beta\rho} ={1-\beta V(0)\over \beta\rho}.$$

Since $\beta_c^{\mathrm G}V(0)=1$, near the Gaussian transition

$$m^2\propto \beta_c^{\mathrm G}-\beta$$

on the disordered side. In terms of temperature this is proportional to $T-T_c^{\mathrm G}$, up to a positive nonuniversal constant.

Exercise 2

Use dimensional analysis to show that

$$G_0(x)=\int {d^D k\over(2\pi)^D}{e^{ik\cdot x}\over k^2}$$

scales as $|x|^{2-D}$ for $D>2$ and $x\ne0$, with the Fourier integral understood as the continuum Green distribution. Explain why one cannot set $D=2$ in the resulting coefficient to obtain a finite constant covariance.

Solution

Let $r=|x|$ and rescale the integration variable by

$$q=kr.$$

Then

$$d^D k={d^D q\over r^D}, \qquad k^2={q^2\over r^2}, \qquad e^{ik\cdot x}=e^{iq\cdot \hat x}.$$

Therefore

$$G_0(x) = \int {d^D q\over(2\pi)^D} {r^2\over r^D} {e^{iq\cdot \hat x}\over q^2} = r^{2-D} \int {d^D q\over(2\pi)^D}{e^{iq\cdot \hat x}\over q^2}.$$

For $D>2$ the infrared integral is locally integrable. At $x\ne0$, removing the UV regulator in the distributional sense leaves a dimensionless coefficient independent of $r$. Thus

$$G_0(x)\propto {1\over r^{D-2}}.$$

For $D>2$ the standard normalization gives

$$G_0(x)= {\Gamma(D/2-1)\over4\pi^{D/2}} {1\over r^{D-2}}.$$

At $D=2$ the coefficient has a pole and the massless infrared integral diverges. With $m>0$, the small-$m$ result is logarithmic; the difference between two nonzero separations tends to $-\log(r/r_0)/(2\pi)$. A scale-independent finite constant would miss that dependence.

Exercise 3

For

$$I_D(m)=\int^\Lambda {d^D q\over(2\pi)^D}{1\over q^2+m^2},$$

show by power counting that $I_D(0)$ has a UV divergence proportional to $\Lambda^{D-2}$ for $D>2$. What happens at $D=2$?

Solution

Using spherical coordinates,

$$I_D(0) = {S_{D-1}\over(2\pi)^D} \int_0^\Lambda dq\,{q^{D-1}\over q^2} = {S_{D-1}\over(2\pi)^D} \int_0^\Lambda dq\,q^{D-3}.$$

For $D>2$,

$$\int_0^\Lambda dq\,q^{D-3} = {\Lambda^{D-2}\over D-2},$$

up to the infrared lower limit. Therefore

$$I_D(0)\sim {S_{D-1}\over(2\pi)^D} {\Lambda^{D-2}\over D-2}.$$

At $D=2$, the integral becomes

$$\int_0^\Lambda {dq\over q},$$

which is logarithmic. In the massless theory this logarithm is also infrared divergent. With nonzero $m$, the denominator $q^2+m^2$ cuts off the infrared region and produces a logarithm of $\Lambda/m$.

Exercise 4

Find the engineering dimensions of $\phi$ and $\lambda$ in

$$S=\int d^D x\, \left[ {1\over2}(\partial\phi)^2 +{\lambda\over4!}\phi^4 \right].$$

Use them to identify the upper critical dimension.

Solution

The action is dimensionless. Since $[d^D x]=-D$ and $[\partial]=1$, the kinetic term gives

$$-D+2+2[\phi]=0.$$

Hence

$$[\phi]={D-2\over2}.$$

The quartic interaction gives

$$-D+[\lambda]+4[\phi]=0.$$

Substituting the field dimension,

$$[\lambda]=D-4[\phi] = D-2(D-2) = 4-D.$$

The coupling is dimensionless when

$$D=4.$$

Therefore the upper critical dimension of the $\phi^4$ interaction is

$$D_c=4.$$

For $D>4$, $[\lambda]<0$ and the interaction is irrelevant at the Gaussian fixed point. For $D<4$, $[\lambda]>0$ and it is relevant.

Exercise 5

In the Gaussian scaling range $a\ll R\ll\xi,L$ for $D>2$, let

$$G_0(R)\sim {1\over R^{D-2}}.$$

Define the connected composite correlator

$$C(R)=\langle \phi^2(x_1)\phi^2(x_2)\rangle_c, \qquad R=|x_1-x_2|.$$

For $\lambda\phi^4/4!$, estimate the graph in which all four vertex fields contract with the endpoints, restricting the vertex to $\mathcal A_R$ defined above. Show that its relative magnitude scales as $\lambda R^{4-D}$. Give the convergence window if one instead integrates this graph over all space with the massless continuum propagator.

Solution

The Gaussian correlator is

$$C_0(R)=2G_0(R)^2 \sim {1\over R^{2D-4}}.$$

The specified restricted graph is

$$\delta C_R^{\mathrm{ins}} =-\lambda \int_{\mathcal A_R} d^D z\, G_0(x_1-z)^2G_0(x_2-z)^2.$$

For a scaling estimate at separation $R$, the integration region has volume $R^D$, and each propagator contributes a factor $R^{-(D-2)}$. Since there are four propagators,

$$|\delta C_R^{\mathrm{ins}}| \sim \lambda R^D \left(R^{-(D-2)}\right)^4 = \lambda R^{D-4D+8} = \lambda R^{-3D+8}.$$

Therefore

$${|\delta C_R^{\mathrm{ins}}|\over C_0(R)} \sim {\lambda R^{-3D+8}\over R^{-2D+4}} = \lambda R^{4-D}.$$

This relative magnitude decreases with $R$ for $D>4$, grows for $D<4$, and is marginal by power counting for $D=4$. Without the restriction, endpoint convergence requires $\int_0 ds\,s^{3-D}$ to converge, while convergence at infinity requires $\int^\infty ds\,s^{7-3D}$. Thus the full massless graph is finite only for $8/3<D<4$. Subtracting the disconnected expectation does not cure these separate UV or IR singularities.

Exercise 6

Use the Schwinger representation

$${1\over k^2+m^2}=\int_0^\infty dL\,e^{-L(k^2+m^2)}$$

to show that a critical Gaussian propagator describes paths with typical length $L\sim R^2$ for endpoint separation $R$.

Solution

Fourier transforming the Schwinger representation gives

$$G_0(x) = \int_0^\infty dL \int {d^D k\over(2\pi)^D} e^{ik\cdot x}e^{-Lk^2}e^{-Lm^2}.$$

The Gaussian momentum integral is

$$\int {d^D k\over(2\pi)^D} e^{ik\cdot x}e^{-Lk^2} = {1\over(4\pi L)^{D/2}} \exp\left(-{x^2\over4L}\right).$$

Thus

$$G_0(x) = \int_0^\infty dL\, {1\over(4\pi L)^{D/2}} \exp\left(-{R^2\over4L}-m^2L\right).$$

At criticality $m=0$. The exponential factor suppresses $L\ll R^2$, because then $R^2/(4L)$ is large. The power $L^{-D/2}$ suppresses very large $L$ in dimensions where the integral is convergent. The natural scaling variable is

$$u={L\over R^2}.$$

Therefore the dominant path length scales as

$$L\sim R^2.$$

This is the Brownian scaling relation: the path has fractal dimension $2$.

References

• Wilson, Kenneth G., and Michael E. Fisher. “Critical Exponents in 3.99 Dimensions.” Physical Review Letters 28, no. 4 (1972): 240–243. DOI.
• Wilson, Kenneth G., and John Kogut. “The Renormalization Group and the $\epsilon$ Expansion.” Physics Reports 12, no. 2 (1974): 75–199. DOI.

Further reading

• Cardy, John. Scaling and Renormalization in Statistical Physics. Cambridge Lecture Notes in Physics 5. Cambridge: Cambridge University Press, 1996. Chapters 2–4.
• Ginzburg, Vitaly L. “Some Remarks on Phase Transitions of the Second Kind and the Microscopic Theory of Ferroelectric Materials.” Russian original, Fizika Tverdogo Tela 2, no. 9 (1960): 2031–2043. Author’s institutional bibliography, item 374.
• Goldenfeld, Nigel. Lectures on Phase Transitions and the Renormalization Group. Frontiers in Physics 85. Reading, MA: Addison-Wesley, 1992. Chapters 4–6.
• Polyakov, Alexander M. Gauge Fields and Strings. Contemporary Concepts in Physics 3. Chur: Harwood Academic Publishers, 1987. Chapters 1–3 and 9.
• Zinn-Justin, Jean. Quantum Field Theory and Critical Phenomena. 5th ed. Oxford: Oxford University Press, 2021. Chapters 14–17.