Momentum-Shell Renormalization Group

Real-space block spins are intuitive, but they carry a nagging ambiguity: what exactly counts as a "short-distance degree of freedom" when interactions spread across a lattice? Kenneth Wilson's answer was to move to Fourier space, where the question has a crisp answer — a short-distance degree of freedom is a field mode ϕ(q)\phi(\mathbf{q}) with large wavevector. Define the theory with an ultraviolet cutoff Λ\Lambda (the smoothed inverse lattice spacing), and coarse-graining becomes something you can do exactly for the free theory and systematically for the interacting one: integrate out the modes in a thin shell Λ/b<∣q∣<Λ\Lambda/b<|\mathbf{q}|<\Lambda near the cutoff. The procedure is three clean steps. First coarse-grain: split the field into slow modes (∣q∣<Λ/b|\mathbf{q}|<\Lambda/b) and fast modes (in the shell), and integrate out the fast modes. Second rescale momenta, q′=b q\mathbf{q}'=b\,\mathbf{q}, to restore the cutoff to Λ\Lambda. Third rescale the field, ϕ′=ϕ/ζ\phi'=\phi/\zeta, to bring the gradient term back to canonical form. The composition is a transformation on the space of couplings, and taking the shell infinitesimally thin (b=eℓb=e^{\ell}, ℓ→0\ell\to0) turns the discrete map into a continuous flow. This article owns the foundations of that machinery: the exact Gaussian shell integration, the engineering dimensions that fall out of pure rescaling, and the tree-level result that the mass term is always relevant (yr=2y_r=2) while the quartic coupling has yu=4−dy_u=4-d — which changes sign at the upper critical dimension dc=4d_c=4. That single sign change is the whole story of why mean-field theory works above four dimensions and fails below. The full nonperturbative payoff — the Wilson–Fisher fixed point and the ϵ=4−d\epsilon=4-d expansion of the exponents — is developed in its own article; here we lay the scaffold on which it stands, and quote its leading result by name.

Key takeaways
  • Wilson's momentum-shell RG coarse-grains in Fourier space: integrate out field modes in a thin shell Λ/b<∣q∣<Λ\Lambda/b<|\mathbf{q}|<\Lambda, then rescale momenta and fields to restore the original cutoff.
  • Demanding the gradient term stay canonical fixes the field rescaling ζ=b(d+2−η)/2\zeta=b^{(d+2-\eta)/2} (with η=0\eta=0 at the Gaussian level); loop corrections to the q2q^2 coefficient define the anomalous dimension η\eta.
  • Pure rescaling gives the tree-level recursions r′=b2rr'=b^2 r and u′=b4−duu'=b^{4-d}u, so yr=2y_r=2 (mass always relevant) and yu=4−dy_u=4-d.
  • yuy_u changes sign at d=4d=4: irrelevant above (Gaussian fixed point, mean-field exponents), relevant below — this is the upper critical dimension dc=4d_c=4.
  • Taking the shell infinitesimal turns the discrete map into continuous beta-function flow dgi/dℓ=βi({g})dg_i/d\ell=\beta_i(\{g\}); the interacting Wilson–Fisher fixed point and ϵ\epsilon-expansion build on this foundation.
  • Key idea: integrate out a thin momentum shell, then rescale momenta and fields to restore the theory; the resulting flow of couplings has fixed points and linearized eigenvalues that give universal critical exponents, and the tree-level power counting alone fixes dc=4d_c=4.
The core insight

Integrate out a thin momentum shell near the cutoff, then rescale momenta and fields to restore the original theory: the resulting transformation of the couplings, taken in the infinitesimal-shell limit, is a continuous flow whose fixed points and linearized eigenvalues give universal critical exponents.

Where this lives today

The momentum-shell RG is the analytic workhorse of critical phenomena. Applied to the single-component ϕ4\phi^4 theory it gives the ϵ\epsilon-expansion of the 3D Ising exponents; extended to an nn-component order parameter it covers the XY (n=2n=2), Heisenberg (n=3n=3), and polymer (n→0n\to0) universality classes from one machine, with the Wilson–Fisher coupling scaling as u∗∝ϵ/(n+8)u^*\propto\epsilon/(n+8). Its modern descendant is the functional (exact) renormalization group of Wetterich and Polchinski, which replaces the sharp shell with a smooth momentum regulator and flows the entire effective action rather than a few couplings — the tool of choice for strongly correlated electrons, the quark–gluon plasma, and even quantum gravity's asymptotic-safety program. The same tree-level power-counting that fixes dc=4d_c=4 here is what particle physicists call "dimensional analysis of operators," and it decides which interactions matter at low energy in every effective field theory, from the Standard Model to the theory of sound in a fluid. Wilson's thin shell turned out to be one of the most portable ideas in physics.

The physical picture

Picture the field's information as spread across momentum scales. A microscope with finite resolution simply cannot see modes near the cutoff; the RG asks what effective theory describes the modes it can see. Integrating the shell is exactly "averaging out" the finest unresolved wiggles — but done in Fourier space, where "fine" means "large ∣q∣|\mathbf{q}|." This is the same coarse-grain-and-rescale program as the renormalization group idea and its real-space block-spin cousin, now with an unambiguous notion of "short distance."

Figure
Momentum-shell coarse-graining: integrate out the fast modes in the shell Λ/b < |q| < Λ, keep the slow modes, then rescale q → bq
Fourier space (radial): modes filled up to cutoff Λ Λ Λ/b slow φ< shell: fast φ> (integrated out) rescale q → bq Cutoff restored to Λ; theory comparable again Λ slow φ′

Rescaling is what makes the procedure a true flow rather than a one-shot blur. After removing the shell, the slow theory has a smaller cutoff Λ/b\Lambda/b; magnifying lengths (equivalently, shrinking momenta back to Λ\Lambda) lets us compare the new theory with the old on the same footing. A fixed point is the self-similar theory that looks identical at every magnification — scale invariance at criticality, expressed in momentum space. The field rescaling factor ζ\zeta is the subtle ingredient, and it is worth pinning down carefully because a common convention slip lives here. The Gaussian part of the action dictates its value: demanding that the coefficient of the q2q^2 gradient term stay fixed after coarse-graining and momentum rescaling forces a specific ζ\zeta. When loop corrections shift the propagator's momentum dependence, absorbing that shift into ζ\zeta introduces a small correction — the anomalous dimension η\eta — measuring how the field's scaling deviates from naive dimensional analysis.

Figure
Tree-level eigenvalues: the mass r is always relevant (y_r = 2); the quartic u has y_u = 4 − d, changing sign at d = 4
d y 0 y_r = 2 (mass always relevant) y_u = 4 − d d_c = 4 d < 4: u relevant d > 4: u irrelevant 2 3 4 5

The last figure is the punchline of the whole method's foundation. Two straight lines: the mass eigenvalue yr=2y_r=2 sits flat and positive (the mass is always relevant), while the quartic eigenvalue yu=4−dy_u=4-d slopes downward and crosses zero exactly at d=4d=4. Above four dimensions the quartic coupling is irrelevant and mean-field theory is exact; below, it is relevant and drives the system to a new fixed point. That crossing point is the upper critical dimension. Relevant couplings (like the reduced temperature) grow under the flow and drive the system away from criticality; irrelevant ones decay, which is why universality holds — the number of relevant directions equals the number of parameters you must tune, typically two for the Ising class.

The core equations

We work with the Landau–Ginzburg–Wilson action for a scalar order-parameter field, defined with an ultraviolet momentum cutoff Λ\Lambda.

Z=∫Dϕ  e−H[ϕ],H=∫∣q∣<Λddq(2π)d 12(r+q2)∣ϕ(q)∣2+u ⁣∫(ϕ2)2(1)Z=\int\mathcal{D}\phi\;e^{-\mathcal{H}[\phi]},\qquad \mathcal{H}=\int_{|\mathbf{q}|<\Lambda}\frac{d^dq}{(2\pi)^d}\,\tfrac12\big(r+q^2\big)|\phi(\mathbf{q})|^2 + u\!\int(\phi^2)^2 \qquad (1)

where: ϕ(q)\phi(\mathbf{q}) — Fourier mode of the order-parameter field; Λ\Lambda — ultraviolet cutoff (inverse microscopic length); rr — quadratic coupling, ∝(T−Tc)\propto(T-T_c) at mean-field level; q2q^2 — gradient (stiffness) coefficient, normalized to one; uu — quartic self-coupling. The Gaussian part 12(r+q2)∣ϕ∣2\tfrac12(r+q^2)|\phi|^2 is the free theory; uu is the interaction that makes criticality nontrivial. The first step splits the field by wavevector into slow and fast modes.

ϕ(q)=ϕ<(q) θ ⁣(Λb−∣q∣)+ϕ>(q) θ ⁣(∣q∣−Λb)θ ⁣(Λ−∣q∣)(2)\phi(\mathbf{q})=\phi^{<}(\mathbf{q})\,\theta\!\big(\tfrac{\Lambda}{b}-|\mathbf{q}|\big)+\phi^{>}(\mathbf{q})\,\theta\!\big(|\mathbf{q}|-\tfrac{\Lambda}{b}\big)\theta\!\big(\Lambda-|\mathbf{q}|\big) \qquad (2)

where: ϕ<\phi^{<} — slow modes retained (∣q∣<Λ/b|\mathbf{q}|<\Lambda/b); ϕ>\phi^{>} — fast modes in the shell Λ/b<∣q∣<Λ\Lambda/b<|\mathbf{q}|<\Lambda; b>1b>1 — RG scale factor; θ\theta — Heaviside step enforcing the shell boundaries. Coarse-graining is then the exact functional integral over the fast modes,

e−H′[ϕ<]=∫Dϕ>  e−H[ϕ<+ϕ>],(3)e^{-\mathcal{H}'[\phi^{<}]}=\int\mathcal{D}\phi^{>}\;e^{-\mathcal{H}[\phi^{<}+\phi^{>}]}, \qquad (3)

producing the effective action H′\mathcal{H}' for the slow modes. Now the rescaling. Restore the cutoff with q′=b q\mathbf{q}'=b\,\mathbf{q} and rescale the field by a factor ζ\zeta chosen to keep the gradient term canonical.

q′=b q,ϕ′(q′)=ζ−1ϕ<(q),ζ=b(d+2−η)/2(4)\mathbf{q}'=b\,\mathbf{q},\qquad \phi'(\mathbf{q}')=\zeta^{-1}\phi^{<}(\mathbf{q}),\qquad \zeta=b^{(d+2-\eta)/2} \qquad (4)

where: q′\mathbf{q}' — rescaled momentum restoring the cutoff to Λ\Lambda; ζ\zeta — field rescaling factor; η\eta — anomalous dimension. At the Gaussian level η=0\eta=0 and ζ=b(d+2)/2\zeta=b^{(d+2)/2}; the derivation below shows this follows from demanding gradient-term invariance. (A frequent slip: the momentum-space field-rescale factor is ζ=b(d+2−η)/2\zeta=b^{(d+2-\eta)/2}, not b(d−2+η)/2b^{(d-2+\eta)/2} — the latter is the real-space scaling dimension Δϕ=(d−2+η)/2\Delta_\phi=(d-2+\eta)/2 of the field, a related but distinct object. Confusing the two flips the tree-level exponents.) Pure rescaling of each term of the action gives its tree-level recursion — the engineering dimensions.

r′=b2 r,u′=b4−d u(Gaussian / tree level)(5)r'=b^{2}\,r,\qquad u'=b^{4-d}\,u\qquad(\text{Gaussian / tree level}) \qquad (5)

where: r′r' — rescaled quadratic coupling; u′u' — rescaled quartic coupling. So yr=2y_r=2 (mass always relevant) and yu=4−dy_u=4-d, which is positive for d<4d<4, zero at d=4d=4, negative for d>4d>4. Taking the shell infinitesimally thin turns the discrete map into continuous flow.

dgidℓ=βi({g}),b=eℓ(6)\frac{dg_i}{d\ell}=\beta_i(\{g\}),\qquad b=e^{\ell} \qquad (6)

where: ℓ=ln⁡b\ell=\ln b — logarithmic RG "time"; βi\beta_i — the beta function of coupling gig_i. Fixed points (βi=0\beta_i=0), the eigenvalues of the linearized flow, and the relevant/irrelevant classification then follow by ordinary-differential-equation analysis. The eigenvalues map to scaling exponents through

λi=byi,yi=ln⁡λiln⁡b,(7)\lambda_i=b^{y_i},\qquad y_i=\frac{\ln\lambda_i}{\ln b}, \qquad (7)

interpreted by linearization about a fixed point and the operator classification, with the exponents themselves owned by critical exponents from the RG.

The core derivations

Derivation AThe Gaussian fixed point and the engineering dimensions

The free theory is where momentum-shell RG is exact, and everything universal about the tree level lives here. The goal is to show that integrating out the shell in the Gaussian theory is trivial, that the rescaling alone determines the field factor ζ\zeta and the recursions r′=b2rr'=b^2 r, u′=b4−duu'=b^{4-d}u, and hence that the Gaussian fixed point sits at the origin with yt=2y_t=2, yu=4−dy_u=4-d. The strategy is pure power-counting: track how the measure, the momentum, and the fields each scale under q′=bq\mathbf{q}'=b\mathbf{q}, ϕ′=ζ−1ϕ<\phi'=\zeta^{-1}\phi^<, and demand each term of the action return to canonical form. The one thing to get right is the direction of every power of bb — a single sign error here propagates into wrong exponents.

Step 1. Shell integration is trivial. With u=0u=0 the action is a decoupled Gaussian, so ∫Dϕ>e−H0[ϕ>]\int\mathcal{D}\phi^{>}e^{-\mathcal{H}_0[\phi^{>}]} is a field-independent constant; H′[ϕ<]\mathcal{H}'[\phi^{<}] is just H0\mathcal{H}_0 restricted to ∣q∣<Λ/b|\mathbf{q}|<\Lambda/b. No new couplings are generated.

Step 2. Fix ζ\zeta from the gradient term. Under q′=bq\mathbf{q}'=b\mathbf{q} the measure gives ddq=b−dddq′d^dq=b^{-d}d^dq' and q2=b−2q′2q^2=b^{-2}q'^2; the field gives ϕ<=ζϕ′\phi^<=\zeta\phi'. Then

12 ⁣∫ddq q2∣ϕ<∣2=12 b−d−2ζ2 ⁣∫ddq′ q′2∣ϕ′∣2.\tfrac12\!\int d^dq\,q^2|\phi^<|^2=\tfrac12\,b^{-d-2}\zeta^2\!\int d^dq'\,q'^2|\phi'|^2.

Canonical form (coefficient back to 12\tfrac12) requires b−d−2ζ2=1b^{-d-2}\zeta^2=1, i.e. ζ=b(d+2)/2\zeta=b^{(d+2)/2} (so η=0\eta=0 at the Gaussian fixed point). Why this and not b(d−2)/2b^{(d-2)/2}: the b−db^{-d} from the measure and b−2b^{-2} from q2q^2 must be undone by ζ2\zeta^2; solving gives ζ2=bd+2\zeta^2=b^{d+2}. The exponent (d−2)/2(d-2)/2 is the real-space scaling dimension of the field, obtained by a different bookkeeping.

Step 3. Recursion for rr. The mass term 12r∫ddq∣ϕ∣2\tfrac12 r\int d^dq|\phi|^2 carries one measure and two fields:

r b−dζ2=r b−dbd+2=r b2 ⇒ r′=b2r,yr=2.r\,b^{-d}\zeta^2=r\,b^{-d}b^{d+2}=r\,b^{2}\ \Rightarrow\ r'=b^2 r,\quad y_r=2.

Step 4. Recursion for uu. The quartic u∫(ddq)3(ϕ)4u\int(d^dq)^3(\phi)^4 carries three measures and four fields:

u b−3dζ4=u b−3db2d+4=u b4−d ⇒ u′=b4−du,yu=4−d.u\,b^{-3d}\zeta^4=u\,b^{-3d}b^{2d+4}=u\,b^{4-d}\ \Rightarrow\ u'=b^{4-d}u,\quad y_u=4-d.

Step 5. Fixed point and eigenvalues. The only fixed point is r∗=u∗=0r^*=u^*=0 (the Gaussian fixed point). The map is already diagonal, so the eigenvalues are b2b^2 and b4−db^{4-d}: yt=2y_t=2 gives ν=1/yt=1/2\nu=1/y_t=1/2 (the mean-field value), and uu is irrelevant precisely for d>4d>4, marking the upper critical dimension dc=4d_c=4. Verified: `verify_momentum-shell-rg.py` solves b−d−2ζ2=1b^{-d-2}\zeta^2=1 for ζ=b(d+2)/2\zeta=b^{(d+2)/2}, derives r′=b2rr'=b^2 r and u′=b4−duu'=b^{4-d}u by power-counting, and confirms the sign change of yuy_u at d=4d=4 (checks 1–5).

Go deeper · Derivation B · 2 minThe Gaussian shell integral and the onset of interactions
Derivation BThe Gaussian shell integral and the onset of interactions

Having fixed the tree level, we set up — but do not fully execute — the interacting calculation, showing concretely where the Wilson–Fisher fixed point comes from while leaving its full construction to the dedicated article. The goal is to write down the Gaussian shell propagator, confirm the shell integral is finite and positive, and exhibit the structure of the one-loop corrections to rr and uu that generate a nontrivial fixed point for d<4d<4. The strategy: expand e−Hinte^{-\mathcal{H}_{\rm int}} in the quartic coupling and average over the fast modes with the Gaussian weight, whose propagator in the shell is 1/(r+q2)1/(r+q^2). We stop at identifying the fixed-point condition; the loop integrals and the ϵ\epsilon-expansion of the exponents are owned by the Wilson–Fisher article.

Step 1. The shell propagator. For fast modes in the shell, the Gaussian two-point function is G>(q)=1/(r+q2)G^{>}(\mathbf{q})=1/(r+q^2), and the one-loop corrections come from contracting fast-mode pairs in the interaction.

Step 2. The shell integral is finite. The relevant integral ∫Λ/b<∣q∣<Λddq(2π)d 1r+q2\int_{\Lambda/b<|\mathbf{q}|<\Lambda}\frac{d^dq}{(2\pi)^d}\,\frac{1}{r+q^2} runs over a compact shell of a smooth positive integrand, so it is finite and positive; for a thin shell it is proportional to the shell width dℓ=ln⁡bd\ell=\ln b.

Step 3. Structure of the corrections. The tadpole (one contracted pair in (ϕ2)2(\phi^2)^2) shifts rr, and the bubble (two four-point contractions) shifts uu:

drdℓ=2r+Au1+r,dudℓ=(4−d) u−Bu2(1+r)2,\frac{dr}{d\ell}=2r+\frac{Au}{1+r},\qquad \frac{du}{d\ell}=(4-d)\,u-\frac{Bu^2}{(1+r)^2},

with A,B>0A,B>0 collecting combinatorial and angular (unit-sphere) factors. The tree-level pieces 2r2r and (4−d)u(4-d)u are exactly the engineering dimensions from the previous derivation; the loop terms are the new physics.

Step 4. Where the new fixed point comes from. Setting du/dℓ=0du/d\ell=0 gives either u∗=0u^*=0 (Gaussian) or the nonzero root

u∗=(4−d)(1+r∗)2B=ϵB+O(ϵ2),ϵ=4−d,u^*=\frac{(4-d)(1+r^*)^2}{B}=\frac{\epsilon}{B}+O(\epsilon^2),\qquad \epsilon=4-d,

which exists only for ϵ>0\epsilon>0 (d<4d<4). This is the Wilson–Fisher fixed point; its full construction, stability analysis, and the ϵ\epsilon-expansion of the exponents are developed there. We have only located the condition that produces it. Verified: `verify_momentum-shell-rg.py` confirms the Gaussian shell integral is finite and positive (d=3d=3 shell), and that du/dℓ=0du/d\ell=0 yields the nonzero root u∗=ϵ(1+r∗)2/Bu^*=\epsilon(1+r^*)^2/B (checks 6–7); the quoted ν=1/2+ϵ/12\nu=1/2+\epsilon/12 is verified as the Wilson–Fisher result (check 8), attributed by name.

Limiting cases

The behaviour of the shell RG depends sharply on dimension and on the chosen scale factor. The exponents in the last column are owned by the critical-exponents and Wilson–Fisher articles; here they interpret the flow this method sets up.

CaseResultMeaning
d>4d>4 (above upper critical dimension)yu=4−d<0y_u=4-d<0The quartic coupling is irrelevant; the flow reaches the Gaussian fixed point and mean-field exponents (ν=1/2\nu=1/2, η=0\eta=0) hold exactly.
d=4d=4 (marginal)du/dℓ∝−u2du/d\ell\propto -u^2uu is marginally irrelevant: it decays logarithmically, giving mean-field exponents with logarithmic corrections.
d<4d<4 (Wilson–Fisher regime)u∗∼ϵ=4−du^*\sim\epsilon=4-dA new stable fixed point at O(ϵ)O(\epsilon) controls criticality; nontrivial exponents organized in powers of ϵ\epsilon (owned by the WF article).
Infinitesimal shell b=eℓ, ℓ→0b=e^\ell,\ \ell\to0dgi/dℓ=βi({g})dg_i/d\ell=\beta_i(\{g\})Discrete recursion becomes autonomous ODE flow — the natural, standard form of the momentum-shell RG.

Worked examples

Example 1 — The engineering dimensions, from power-counting alone. We compute the tree-level eigenvalues explicitly, the calculation that fixes the upper critical dimension.

Step 1. From the Gaussian gradient term, ζ=b(d+2)/2\zeta=b^{(d+2)/2} (Step 2 of the first derivation). Take a concrete d=3d=3: ζ=b5/2\zeta=b^{5/2}.

Step 2. Mass recursion: r′=r b−dζ2=r b−3b5=r b2r'=r\,b^{-d}\zeta^2=r\,b^{-3}b^{5}=r\,b^{2}, so yr=2y_r=2, independent of dd.

Step 3. Quartic recursion: u′=u b−3dζ4=u b−9b10=u b1u'=u\,b^{-3d}\zeta^4=u\,b^{-9}b^{10}=u\,b^{1}, so yu=1y_u=1 at d=3d=3. Since yu=4−dy_u=4-d, this is positive — the quartic is relevant in three dimensions, so mean-field theory fails and a Wilson–Fisher fixed point exists.

Step 4. Repeat at d=5d=5: ζ=b7/2\zeta=b^{7/2}, u′=u b−15b14=u b−1u'=u\,b^{-15}b^{14}=u\,b^{-1}, so yu=−1<0y_u=-1<0 — irrelevant, mean-field exact. The sign flip between d=3d=3 and d=5d=5 brackets the upper critical dimension at dc=4d_c=4, where yu=0y_u=0.

Example 2 — The Gaussian fixed point gives the mean-field correlation-length exponent. We check that the tree-level thermal eigenvalue reproduces the known mean-field ν\nu.

Step 1. At the Gaussian fixed point the thermal eigenvalue is yt=yr=2y_t=y_r=2 (the mass is the thermal scaling field, r∝T−Tcr\propto T-T_c).

Step 2. The correlation-length exponent is ν=1/yt=1/2\nu=1/y_t=1/2.

Step 3. Cross-check against mean-field theory: Landau theory gives ξ∼∣t∣−1/2\xi\sim|t|^{-1/2}, i.e. ν=1/2\nu=1/2. They agree — as they must, since above dc=4d_c=4 the Gaussian fixed point is the exact critical theory. ν=1/2\nu=1/2 (mean-field), exact for d>4d>4.

Step 4. For d<4d<4 the flow leaves the Gaussian fixed point along the relevant uu direction and settles at the Wilson–Fisher fixed point, where yt=2−ϵ/3+O(ϵ2)y_t=2-\epsilon/3+O(\epsilon^2) and hence ν=1/2+ϵ/12+⋯\nu=1/2+\epsilon/12+\cdots — the standard one-loop result, quoted here by name from the Wilson–Fisher article and verified numerically.

Origins · 1 minwilsonian renormalization and how this was discovered
WR
wilsonian renormalization
Wilson's momentum-space formulation of the RG (1971) · Wilsonian renormalization is the viewpoint, introduced by Kenneth Wilson around 1971, that a physical theory is defined with a finite cutoff and that lowering the cutoff by integrating out high-momentum modes generates an effective theory with running couplings. In this picture renormalizability is not a formal requirement but a statement about which couplings are relevant near a fixed point, and the same machinery that tames ultraviolet divergences in quantum field theory computes critical exponents in statistical mechanics. The momentum-shell implementation — integrate a thin shell, rescale, repeat — is the concrete realization that turned this viewpoint into a calculational tool and earned Wilson the 1982 Nobel Prize.

History
1954
— Gell-Mann and Low introduce renormalization-group functional equations in quantum electrodynamics (Physical Review 95, 1300), the field-theory origin of the idea
1966
— Kadanoff's block-spin scaling supplies the real-space intuition Wilson would formalize
1971
— Wilson formulates the momentum-space renormalization group and the notion of integrating out modes shell by shell (Physical Review B 4)
1972
— Wilson and Fisher introduce the ε = 4 − d expansion, computing critical exponents in 3.99 dimensions (Physical Review Letters 28, 240)
1982
— Wilson receives the Nobel Prize in Physics for the theory of critical phenomena in connection with phase transitions

Going further

References. Kardar, Statistical Physics of Fields (2007), Ch. 5 · Pathria & Beale, Statistical Mechanics (3rd ed.), Ch. 14 · Goldenfeld, Lectures on Phase Transitions and the RG (1992), Ch. 9 · Wilson & Kogut, Phys. Rep. 12 (1974) 75 · Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity (2nd ed.), Ch. 12.

Provenance. Equations and structure follow Kardar, Statistical Physics of Fields Ch. 5, and Goldenfeld Ch. 9. A verifier should recompute: (1) the engineering dimensions yr=2y_r=2, yu=4−dy_u=4-d from the rescaling power-counting; (2) the field rescaling ζ=b(d+2−η)/2\zeta=b^{(d+2-\eta)/2} (momentum space) from gradient-term invariance; (3) the finiteness of the Gaussian shell integral; (4) the fixed-point condition u∗=ϵ(1+r∗)2/Bu^*=\epsilon(1+r^*)^2/B and the quoted one-loop ν=1/2+ϵ/12\nu=1/2+\epsilon/12 (owned by the Wilson–Fisher article). Standard reference numbers: mean-field ν=1/2\nu=1/2; 3D Ising ν≈0.630\nu\approx0.630, η≈0.036\eta\approx0.036, ω≈0.8\omega\approx0.8 (conformal bootstrap / Monte Carlo).

Note on a corrected base error. The source draft's physical-picture prose wrote the field-rescale factor as ζ=b(d−2+η)/2\zeta=b^{(d-2+\eta)/2}; the correct momentum-space factor is ζ=b(d+2−η)/2\zeta=b^{(d+2-\eta)/2} (the former is the real-space scaling dimension Δϕ\Delta_\phi). The tree-level recursions r′=b2rr'=b^2 r, u′=b4−duu'=b^{4-d}u follow only from the correct factor; this has been fixed throughout and flagged as a teaching point.

Builds on: The Renormalization Group Idea: Coarse-Graining and Rescaling · Block-Spin (Real-Space) Renormalization Group · RG Flow in Coupling Space. Leads to: Wilson-Fisher Fixed Point · Gaussian Fixed Point and Upper Critical Dimension · The Epsilon-Expansion · Critical Exponents from the RG.

Why work in momentum space instead of real space like the block-spin method?

Because "short-distance degree of freedom" is ambiguous on a lattice — interactions spread, and it is not obvious which spins to average — but in Fourier space it is crisp: a short-distance mode is one with large wavevector ∣q∣|\mathbf{q}|. Integrating out a thin momentum shell near the cutoff is therefore unambiguous, and for the Gaussian theory it is exact. The trade-off is that you need a continuum field-theory description (the Landau–Ginzburg–Wilson action) rather than a lattice model, and you gain a systematic small parameter, ϵ=4−d\epsilon=4-d, that the real-space method lacks. The two approaches give the same fixed points and exponents.

What exactly is the anomalous dimension η\eta, and why is it zero at the Gaussian fixed point?

η\eta measures how the field's scaling deviates from naive dimensional analysis. At the Gaussian (free) fixed point the field rescaling ζ=b(d+2)/2\zeta=b^{(d+2)/2} keeps the q2q^2 gradient term exactly canonical, so there is no deviation and η=0\eta=0. Once interactions are switched on, loop corrections modify the momentum dependence of the propagator — the coefficient of q2q^2 shifts — and to keep the gradient term canonical you must adjust ζ\zeta to b(d+2−η)/2b^{(d+2-\eta)/2}. That required adjustment is η\eta. It is small: for the 3D Ising class η≈0.036\eta\approx0.036, so the field scales almost, but not quite, by its engineering dimension.

Why does the quartic coupling being irrelevant above four dimensions mean mean-field theory is exact there?

Mean-field (Landau) theory ignores fluctuations, keeping only the Gaussian part of the action plus the uniform value of the quartic term. The RG makes precise when that is justified: if the quartic coupling uu is irrelevant (yu=4−d<0y_u=4-d<0, i.e. d>4d>4), it flows to zero under coarse-graining, so at long distances the theory is Gaussian and mean-field exponents hold exactly. For d<4d<4, uu is relevant, grows under the flow, and the system is driven to the interacting Wilson–Fisher fixed point where fluctuations shift the exponents. This is the Ginzburg criterion made rigorous.

Where does the Wilson–Fisher fixed point come from, and why isn't it derived in full here?

It comes from the competition in the quartic flow du/dℓ=(4−d)u−Bu2/(1+r)2du/d\ell=(4-d)u-Bu^2/(1+r)^2: the tree-level term (4−d)u(4-d)u pushes uu up (for d<4d<4) while the one-loop bubble term −Bu2-Bu^2 pushes it down, and they balance at u∗∝ϵ=4−du^*\propto\epsilon=4-d. This article sets up that condition to show where the fixed point originates, but its full construction — the loop integrals, the stability analysis, and the ϵ\epsilon-expansion of every exponent — is owned by the dedicated Wilson–Fisher fixed-point article to keep the scope clean. Here we establish only the Gaussian foundation and the tree-level power-counting the WF calculation is built on.

Why is r∗r^* negative at the interacting fixed point?

The tadpole correction +Au/(1+r)+Au/(1+r) to dr/dℓdr/d\ell is positive when u∗>0u^*>0, so for the flow to sit still (dr/dℓ=0dr/d\ell=0) the tree-level term 2r2r must be negative — forcing r∗<0r^*<0 to leading order. Physically, the interaction shifts the critical temperature: the bare mass parameter rr at which the system is actually critical lies below its mean-field value of zero, because fluctuations help disorder the system, so you must start colder (more negative rr) to reach criticality. The precise value is part of the Wilson–Fisher calculation.

Is the sharp momentum-shell cutoff a problem?

At one loop it is fine and physically transparent. Beyond one loop a sharp cutoff generates non-analytic, scheme-dependent contributions, and the shell procedure does not commute cleanly with the infinitesimal-shell limit for every observable. Practitioners therefore switch to smooth cutoffs (as in the functional RG) or to dimensional regularization for higher-order calculations. None of this affects the universal quantities — the exponents are scheme-independent — but it does affect intermediate, non-universal numbers, so the sharp cutoff is best regarded as a pedagogical and one-loop tool rather than a high-precision engine.

Before this

Where this leads