The Renormalization Group for the Working Mathematician

In which throwing away detail, on purpose and repeatedly, turns out to be the whole idea.

In about 1977 a post-graduate mentor was asked by a nine-year-old what the cutting-edge fields in advanced mathematics were. He answered: renormalization group theory — and it will probably win a Nobel prize. Kenneth Wilson collected it in 1982. The method collected a second one in 2004. This page is what he was pointing at.

1. The question that needed answering

Water boils. A magnet loses its magnetism at the Curie point. A random network of pipes suddenly percolates. A forest fire either dies out or spreads forever. These systems share no ingredients: one is molecules with dipole interactions, one is electron spins, one is pure graph combinatorics.

Near their transitions, they share numbers.

Measure how the order parameter grows below the transition and you get a power law with an exponent β. Measure how the susceptibility diverges and you get γ. Measure how far correlations reach and you get ν. Wildly different microscopic systems return the same values, often to several decimal places. Liquid–gas and uniaxial magnets: same numbers. This is universalityDifferent microscopic systems exhibiting identical critical exponents and scaling functions. Not a coincidence and not an analogy — a consequence of what survives repeated coarse-graining., and before 1971 nobody could explain it.

microscopic models repeated coarse-graining one fixed point uniaxial magnet liquid–gas binary alloy lattice gas K* the fixed point ν = 0.62997 β = 0.32642 γ = 1.23708 3D Ising class:
Figure 1 — the universality funnel. Four microscopic models that share no physical ingredients flow, under repeated coarse-graining, to the same fixed point — and inherit its exponents. The 3D Ising values shown are from the conformal bootstrap.5 The whole of this page is an explanation of that arrow.

2. Why criticality is the special case

The key object is the correlation lengthThe distance over which fluctuations remain statistically related. Below it, the system looks correlated; beyond it, independent. Written ξ (xi). ξ: the distance over which one part of the system still “knows about” another. Away from a transition it is a few lattice spacings, and the system looks like independent patches. As you tune towards the critical point it diverges:

ξ ∼ |t|−ν,    t = (T − Tc) / Tc

At t = 0 there is no characteristic length left. Fluctuations exist at every scale simultaneously, nested in one another, and a photograph of the system zoomed by any factor is statistically the same photograph. The system is scale invariant.

T >> Tₜ ξ ≈ a few sites T → Tₜ ξ growing T = Tₜ ξ = ∞ — every scale at once
Figure 2 — the correlation length diverges. Away from criticality the system is a mosaic of independent patches. At the critical point the patches are of every size at once and no scale is preferred. Scale invariance is what makes the next step legal.
The condition, stated once and enforced throughout

Everything below requires ξ → ∞. A first-order transition — ice to water, most condensation, and quite possibly many transitions that look critical from a distance — has a finite correlation length, latent heat, and coexisting phases. There is no scale invariance, so there is no fixed point, no universality and no exponents. Ask whether a transition is continuous before reaching for any of this.

3. Kadanoff’s idea: throw away detail on purpose

Leo Kadanoff, 1966. If the system looks the same at every scale, then describe it at a coarser scale and see what happens to the description.

Take a lattice of Ising spins. Group them into blocks of side b. Replace each block by a single spin, by majority rule. Now shrink all lengths by b so the new lattice has the same spacing as the old one. You have a system of the same type — spins on a lattice — but with different interaction strengths.

stage 1 — the fine lattice a = lattice spacing  ·  ξ = 27a stage 2 — block by majority each 3×3 block → one spin by majority rule detail below scale 3a is integrated out — gone for good stage 3 — rescale lengths by b same lattice spacing a, new couplings K′, and ξ′ = ξ/b = 9a
Figure 3 — one step of the renormalization group. Block, average, rescale. The system that comes out is the same kind of system with different coupling constants, and its correlation length has been divided by b. Repeating this is the entire method. Notice what is lost: everything below the block size, permanently. That loss is not a defect — it is the mechanism.
Definition — the RG transformation

Let K = (K1, K2, …) be the full set of coupling constants describing a system — nearest-neighbour coupling, next-nearest, four-spin terms, everything the symmetry allows. A coarse-graining by factor b defines a map

K′ = Rb(K)

The map is not invertible — information was discarded — so Rb forms a semigroup, not a group: Rb ∘ Rb′ = Rbb′. The name “renormalization group” is, strictly, a misnomer that stuck.

4. Fixed points, and what flows where

Iterating Rb traces a trajectory through the space of all possible couplings. The interesting points are the ones that do not move:

Definition — fixed point
K* = Rb(K*)

A system sitting at K* is unchanged by coarse-graining: it looks identical at every scale. Since each step divides the correlation length by b, a fixed point must have ξ* = ξ*/b, so ξ* is either 0 or ∞. The ξ = 0 fixed points are the boring ones — perfect order and perfect disorder. The ξ = ∞ fixed point is the critical point.

K₁ K₂ critical surface ξ = ∞ on this line T = 0 ordered ξ = 0, stable T = ∞ disordered ξ = 0, stable K* critical, ξ = ∞ UNSTABLE — one way out
Figure 4 — the flow. Almost every starting system flows to a boring fixed point: perfectly ordered or perfectly disordered. The critical fixed point sits on a separatrix — the critical surface — and is unstable. That instability is why you have to tune a system to criticality: temperature must be set to a particular value to land on the surface. Codimension is the number of knobs you must set. One relevant direction, one knob, codimension 1.

5. Linearize, and universality falls out

Near the fixed point write K = K* + δK and linearize:

δK′ = M δK,   Mij = ∂K′i/∂Kj at K*

Diagonalize M. Its eigenvalues are conventionally written λi = byi, and the eigenvectors are the scaling directions. Three cases, and everything hangs on them:

CaseNameBehaviour under coarse-grainingConsequence
yi > 0relevantgrows; the flow runs away from the fixed pointmust be tuned to zero to reach criticality — one knob each
yi < 0irrelevantshrinks; the flow forgets itthis is universality — the microscopic details live here
yi = 0marginalundecided at linear orderneeds higher order; can give logarithms
K* relevant — y > 0 flow runs AWAY irrelevant — y < 0 flow runs IN Three systems with completely different microscopic couplings collapse onto the same relevant axis — after which they are indistinguishable, and share every exponent.
Figure 5 — why universality is not a coincidence. Irrelevant directions contract under the flow. Two systems differing only in irrelevant couplings become the same system after enough steps. What survives is only the relevant directions, and there are very few of them.
A universality class is fixed by a short list: spatial dimensionality, the symmetry and number of components of the order parameter, and the range of the interactions. Everything else is irrelevant — technically, not rhetorically.

Exponents are eigenvalues

For a system with one thermal relevant direction yt and one magnetic one yh, in dimension d, all the familiar exponents are arithmetic on those two numbers:

ExponentIn terms of eigenvaluesWhat it measures
ν1 / ytdivergence of the correlation length
α2 − d/ytspecific-heat singularity
β(d − yh) / ytgrowth of the order parameter
γ(2yh − d) / ytdivergence of susceptibility
δyh / (d − yh)response at the critical isotherm
ηd + 2 − 2yhanomalous dimension of the correlator

Two independent numbers generate six exponents, which is why they obey identities — Rushbrooke α + 2β + γ = 2, Widom γ = β(δ−1), Fisher γ = ν(2−η), and hyperscaling 2 − α = dν. These were empirical puzzles before the RG and are corollaries after it. That is the single most convincing thing the theory did.

Hyperscaling is the one that breaks

The first three relations follow from the two-eigenvalue structure alone. 2 − α = dν additionally assumes the free energy scales purely as ξ−d, and it fails above the upper critical dimension, in the presence of dangerous irrelevant variables, and in random-field systems. A measured hyperscaling violation is therefore informative, not necessarily an error — but it does rule out the simple picture.

6. Worked example A — the 1D Ising chain, exactly

This one you can do by hand in ten minutes, and it is honest: it produces a real answer and the answer is “no phase transition”.

Take a chain with coupling K = J/kBT. Sum out every second spin — decimationA real-space RG in which a subset of degrees of freedom is summed over exactly, leaving an effective model for the rest. Exact in 1D; approximate in higher dimensions., with b = 2. For the middle spin σ2 between neighbours σ1, σ3:

∑σ2=±1 eK(σ1σ2 + σ2σ3) = 2 cosh(K(σ1+σ3))

Demand that this equal eg + K′σ1σ3. Taking the two cases σ1σ3 = ±1 and dividing:

The 1D Ising recursion
K′ = ½ ln cosh(2K)

Fixed points: K* = 0 (infinite temperature) and K* = ∞ (zero temperature). Nothing in between.

Check the stability. For small K, cosh(2K) ≈ 1 + 2K² so K′ ≈ K² < K — the flow moves towards zero. For large K, cosh(2K) ≈ e2K/2 so K′ ≈ K − ½ln2 — it decreases by a constant each step and eventually falls to zero too.

K* = 0 T = ∞, stable K* = ∞ T = 0, unstable every finite temperature flows to disorder ∴ no ordered phase at any T > 0 in one dimension
Figure 6 — the 1D Ising flow. Both fixed points are at ξ = 0 or unreachable; there is no critical fixed point at finite coupling, so there is no phase transition. This reproduces a known exact result — and it shows the RG returning a negative answer cleanly, which is how you know the method has content.

7. Worked example B — Wilson, Fisher, and dimension as a dial

The 1966 picture was physically compelling and computationally stuck. Wilson’s move (1971) was to coarse-grain in momentum space instead: integrate out field modes in a shell between Λ/b and Λ, then rescale. Applied to the φ4 theory this gives flow equations for the mass and quartic coupling that can actually be solved — perturbatively, in an unexpected small parameter.

Wilson and Fisher, 1972: treat the dimension itself as continuous and expand around four. Write d = 4 − ε.

d > 4 Gaussian STABLE u → 0 : mean-field exponents are exact d < 4 Gaussian UNSTABLE Wilson–Fisher STABLE, u* ~ ε a new fixed point takes over and sets the exponents Upper critical dimension dₜ = 4 for this class. Above it, fluctuations are irrelevant and mean-field theory is exact. Below it, fluctuations are relevant and mean-field is wrong. The exchange of stability at d = 4 is the whole content of the epsilon-expansion.
Figure 7 — exchange of stability at four dimensions. The Gaussian fixed point hands control to the Wilson–Fisher fixed point as dimension drops below four. Because the new fixed point sits at u* ∼ ε, everything about it can be computed as a series in ε.

For the one-component (Ising) case the results begin

ν = ½ + ε/12 + O(ε²),     η = ε²/54 + O(ε³)
Read the small print on that number

In three dimensions ε = 1, which is not small. First order gives ν ≈ 0.583 against a true value of 0.62997 — about 7% off. The series is asymptotic, not convergent; adding terms naively eventually makes it worse, and the good modern values come from resummation or from the conformal bootstrap, not from truncating this. The ε-expansion is a proof of mechanism, not a calculator.

Exponents across a few classes

Universality classβγνNotes
Mean field / Landau1/211/2exact for d > 4
2D Ising (Onsager, exact)1/87/41α = 0 (logarithmic)
3D Ising (bootstrap)0.326421.237080.62997the most precisely known non-trivial class5
2D percolation (exact)5/3643/184/3geometric, no Hamiltonian needed
Directed percolation (1+1)0.2765—ν⊥ 1.0970the generic absorbing-state class6

8. The different renormalization groups

“The RG” names a strategy, not a single algorithm. The strategy is: define a coarse-graining, follow what it does to the description, find the fixed points. A dozen quite different machines implement it.

FlavourCoarse-graining is…Best atOrigin
Real-space / block-spingrouping sites, majority ruleintuition; exact in 1D by decimationKadanoff 1966
Momentum-shell (Wilsonian)integrating out a shell of high-k modesthe ε-expansion; critical exponentsWilson 19711
Field-theoreticchanging the renormalization scale μ; beta functionsparticle physics; running couplingsGell-Mann & Low 1954; Callan, Symanzik 1970
Functional / exact (FRG)a flowing effective action Γknon-perturbative problems; no small parameter neededWetterich 1993
Numerical RGiterative diagonalization of a discretized bathimpurity problems — it solved the Kondo problemWilson 1975
DMRGkeeping the highest-weight eigenstates of the reduced density matrix1D quantum chains, to machine precisionWhite 19922
Tensor-network / MERA / TNRremoving short-range entanglement, then coarse-grainingquantum states; flows in wavefunction spaceVidal 2007; Evenbly & Vidal 2015
Conformal bootstrapnot an RG at allthe most precise exponents we have; uses only symmetry and consistencyRattazzi et al. 2008

The last row is there deliberately. The bootstrap gets better numbers than any RG by never coarse-graining anything — it imposes conformal symmetry and crossing consistency and squeezes the allowed region until one point is left. Worth knowing that the RG is the explanatory framework and no longer the most precise computational one.

Beta functions, and the second Nobel

In the field-theoretic version one tracks how a coupling depends on the energy scale at which you measure it:

β(g) = μ ∂g/∂μ

The sign is everything. In QED β > 0: the coupling grows with energy. In QCD the non-abelian gluon self-interaction flips it — β < 0 — so the coupling shrinks at high energy and quarks behave as if free when probed hard. That is asymptotic freedom, and it is why the RG collected a second Nobel prize in 2004.3

energy scale μ → g QED   β > 0 coupling grows — Landau pole QCD   β < 0 asymptotic freedom A zero of β(g) is a fixed point. A negative slope through it means the coupling is driven there at high energy.
Figure 8 — running couplings. Same mathematics as the critical-phenomena flow, different application: here the flow parameter is energy rather than block size. Fixed points of β are scale-invariant theories in both readings.

9. Where the renormalization group does not reach

A tutorial that only sells the method is not useful. Four boundaries, all sharp.

First-order transitions

No diverging correlation length, no scale invariance, no fixed point. There is nothing for the RG to compute. Before applying any of this, establish that the transition is continuous.

Systems with no local order parameter

The Landau–Ginzburg–Wilson framework begins by writing a free energy as a functional of a local order parameter. Topological orderA gapped quantum phase distinguished by no local measurement — no symmetry broken, no local order parameter. Classified instead by anyon content and modular data. is defined by not having one. Three independent sources state the consequence in print:

Tensor-network RG does flow in wavefunction space and does find fixed points — but there is no scaling theory around them, so there are no exponents. This is an open problem, precisely stated, not an oversight.

Combinatorial and computational structures

An autocatalytic set is a set-theoretic object: reflexively autocatalytic and food-generated, or not. There is no scale parameter, no coarse-graining that preserves closure in a controlled way, and hence no RG. Computation is worse: the natural invariant — Turing degree — is degenerate, since almost everything non-trivial is universal. The nearest thing to an exception is coarse-graining of cellular automata.10

Statics and dynamics are different classifications

This one causes real confusion, so it gets its own box.

β and z are indexed by different things

Static exponents are fixed by dimensionality, order-parameter symmetry and interaction range. The dynamic exponent z — relating relaxation time to length, τ ∼ ξz — is fixed by something else entirely: conservation laws and mode couplings. Hohenberg and Halperin classified dynamics into separate models A, B, C… for exactly this reason.4 A non-conserved order parameter gives z ≈ 2.17; a conserved one gives z = 4 — with identical static exponents.

So two systems sharing β while differing in z is not an anomaly and not a partial failure. It is what the theory predicts. Their static classification agrees and their dynamic classification does not. Reporting that as a puzzle means the two classifications have been silently merged.

10. Is the RG the mechanism of domain convergence?

Partly, and the boundary is worth being exact about — because “different domains share numbers” is precisely the observation the RG was built to explain, and precisely the observation that is easiest to over-claim.

What the RG genuinely delivers. A mechanism — not an analogy — by which systems with unrelated microphysics end up with identical critical behaviour. The convergence is real, it is derivable, and it comes with an index: the short list of dimensionality, symmetry and range. Two systems converge if and only if those agree.

What it does not deliver. Convergence of anything else. Sharing a universality class means sharing exponents and scaling functions — it does not mean sharing a mechanism, a structure, or a vocabulary. Two systems in the 3D Ising class have nothing else in common, and that is the point of the theorem rather than a limitation of it.

The renormalization group explains why unrelated systems share numbers.
It says nothing whatever about why they should share structure.

So: as a candidate mechanism for domain convergence, the RG is the strongest one available, and it comes with a built-in test. If two domains are claimed to converge, the RG says to check whether their dimensionality, order-parameter symmetry and interaction range agree — and if they do not, the shared numbers need a different explanation. It also says which quantities should transfer (static exponents) and which have no business transferring (dynamic ones, absent matching conservation laws). That asymmetry is a prediction, and it is checkable.

11. What to remember

  1. Criticality means the correlation length diverges and no scale is special.
  2. Coarse-grain, rescale, repeat. The map on couplings is the renormalization group — a semigroup, because information is destroyed at each step.
  3. Fixed points are scale-invariant theories. Critical points are the unstable ones.
  4. Linearize about the fixed point. Relevant directions must be tuned; irrelevant ones are forgotten.
  5. Universality is the forgetting. The class is indexed by dimension, symmetry, and range — nothing else.
  6. Exponents are eigenvalues, which is why they satisfy identities.
  7. Static and dynamic exponents belong to different classifications. Do not merge them.
  8. No diverging correlation length, no local order parameter, or no scale parameter at all — then no RG.

Coda.

Around 1977, a nine-year-old on a college campus asked his mentor what was cutting-edge in advanced mathematics. The answer was renormalization group theory, with the prediction that it would win a Nobel prize.

Kenneth Wilson received the 1982 Nobel Prize in Physics “for his theory for critical phenomena in connection with phase transitions”. The prediction was right within five years. It was right a second time in 2004, when Gross, Politzer and Wilczek received the prize for asymptotic freedom — a statement about the sign of a beta function, and therefore also a renormalization-group result.

Written for Bruce Stephenson, and for Bruce A Metzger, who was right.

  1. Wilson, K. G. (1971). Renormalization group and critical phenomena I & II. Phys. Rev. B 4:3174 and 3184. doi:10.1103/PhysRevB.4.3174 ↑
  2. White, S. R. (1992). Density matrix formulation for quantum renormalization groups. Phys. Rev. Lett. 69:2863. doi:10.1103/PhysRevLett.69.2863 ↑
  3. Gross, D. J. & Wilczek, F. (1973). Phys. Rev. Lett. 30:1343; Politzer, H. D. (1973). Phys. Rev. Lett. 30:1346. Nobel Prize in Physics 2004. ↑
  4. Hohenberg, P. C. & Halperin, B. I. (1977). Theory of dynamic critical phenomena. Rev. Mod. Phys. 49:435. doi:10.1103/RevModPhys.49.435 ↑
  5. Kos, F., Poland, D., Simmons-Duffin, D. & Vichi, A. (2016). Precision islands in the Ising and O(N) models. JHEP 08:036. doi:10.1007/JHEP08(2016)036 ↑
  6. Hinrichsen, H. (2000). Non-equilibrium critical phenomena and phase transitions into absorbing states. Adv. Phys. 49:815. doi:10.1080/00018730050198152 ↑
  7. Burnell, F. J. (2018). Anyon condensation and its applications. Annu. Rev. Condens. Matter Phys. 9:307. doi:10.1146/annurev-conmatphys-033117-054154 ↑
  8. Eliëns, I. S., Romers, J. C. & Bais, F. A. (2014). Diagrammatics for Bose condensation in anyon theories. Phys. Rev. B 90:195130. doi:10.1103/PhysRevB.90.195130 ↑
  9. Kong, L. (2014). Anyon condensation and tensor categories. Nucl. Phys. B 886:436. doi:10.1016/j.nuclphysb.2014.07.003 ↑
  10. Israeli, N. & Goldenfeld, N. (2004). Computational irreducibility and the predictability of complex physical systems. Phys. Rev. Lett. 92:074105. doi:10.1103/PhysRevLett.92.074105 ↑
← Web Annex