Scaling Fields and Scaling Variables

The renormalization group replaces a Hamiltonian by a coarse-grained one with rescaled couplings, generating a flow K→K′=Rb(K)\mathbf{K}\to\mathbf{K}'=R_b(\mathbf{K}) in the infinite-dimensional space of coupling constants. Critical behaviour is controlled by fixed points K∗=Rb(K∗)\mathbf{K}^*=R_b(\mathbf{K}^*), where the physics looks the same at every magnification. But the physical knobs a laboratory turns — temperature, magnetic field, lattice anisotropy, further-neighbour couplings — are almost never the variables in which this flow is simplest. Turn on the temperature deviation tt and, after one coarse-graining step, you find not just a renormalized tt but a little admixture of every other operator the symmetry allows. Scaling fields are the special coordinates that untangle this. They are the eigen-directions of the linearized RG map near the fixed point, built so that each one transforms multiplicatively under coarse-graining: ui→byiuiu_i\to b^{y_i}u_i. In these variables the flow becomes a set of independent expansions and contractions, and the eigenvalue exponents yiy_i sort every perturbation into three classes — relevant (yi>0y_i>0, grows), irrelevant (yi<0y_i<0, dies), or marginal (yi=0y_i=0). For an Ising-like critical point exactly two scaling fields are relevant: a thermal one utu_t with eigenvalue yt=1/νy_t=1/\nu and a magnetic one uhu_h with eigenvalue yhy_h. Everything else flows to zero. The payoff is that these two numbers fix all of thermodynamics. Written in scaling fields the singular free energy is a generalized homogeneous function, and choosing the rescaling factor to pin the leading relevant field to unity collapses it onto a one-variable scaling function. Out drop the exponents α,β,γ,δ,ν,η\alpha,\beta,\gamma,\delta,\nu,\eta and the identities relating them. The distinction between a scaling field and the bare parameter it resembles is not pedantry: utu_t is an analytic but genuinely nonlinear function of tt, and it is exactly that nonlinearity that keeps measured exponents clean while pushing the messiness into non-universal amplitudes.

Key takeaways
  • Near an RG fixed point the flow is governed by its linearization; the eigenvectors of the linearized RG map define scaling fields uiu_i, nonlinear combinations of the physical couplings that transform multiplicatively under coarse-graining: ui→byiuiu_i \to b^{y_i} u_i.
  • The eigenvalue exponents yiy_i classify perturbations as relevant (yi>0y_i>0), irrelevant (yi<0y_i<0), or marginal (yi=0y_i=0); only the relevant ones survive asymptotically, which is why so few numbers control criticality.
  • Multiplicative scaling is not assumed — it is forced by the semigroup law Rb1∘Rb2=Rb1b2R_{b_1}\circ R_{b_2}=R_{b_1 b_2}, whose only continuous solution is a power byb^{y}. Power-law singularities and scaling fields are the same fact.
  • Written in scaling fields the singular free energy obeys fs({ui})=b−dfs({byiui})f_s(\{u_i\})=b^{-d}f_s(\{b^{y_i}u_i\}); the two relevant Ising fields utu_t (with yt=1/νy_t=1/\nu) and uhu_h set every static exponent.
  • Scaling variables are the same objects viewed as curvilinear coordinates that linearize the flow; the nonlinear (Wegner) redefinition removes regular backgrounds and generates the correction-to-scaling exponent ω\omega.
  • Key idea: Diagonalize the linearized RG near a fixed point; its eigen-directions are scaling fields uiu_i that scale multiplicatively as ui→byiuiu_i\to b^{y_i}u_i, and in these variables the singular free energy is a generalized homogeneous function, fs=b−dfs({byiui})f_s=b^{-d}f_s(\{b^{y_i}u_i\}).

Where this lives today

Every modern extraction of a critical exponent — from neutron scattering on a magnet to Monte Carlo on a lattice to the conformal bootstrap for the 3D Ising model — is, underneath, a statement about scaling fields. When experimentalists fit a susceptibility to χ=Γ∣t∣−γ(1+a∣t∣θ+⋯ )\chi=\Gamma|t|^{-\gamma}(1+a|t|^{\theta}+\cdots), the leading term uses the relevant thermal field and the correction term is precisely the leftover from the leading irrelevant field, with θ=νω\theta=\nu\omega. Getting the nonlinear part of utu_t right is what lets high-precision work quote exponents to four significant figures without them drifting as the fitting window changes. The language also drives condensed-matter theory far from magnets. In quantum criticality the tuning parameter is a coupling gg at zero temperature, and the same relevant/irrelevant bookkeeping tells you which microscopic terms can be dropped and which must be kept — the reason a handful of "field theories at a fixed point" describe superconductors, quantum Hall transitions, and heavy-fermion critical points that share nothing microscopically.

The physical picture

Picture the RG flow as a fluid streaming past a stagnation point — the fixed point. Right at the stagnation point the flow is arrested; a little away from it the velocity field is linear, and its eigen-directions are special streamlines. Along an unstable direction (a relevant field, yi>0y_i>0) you are swept away exponentially fast; along a stable one (irrelevant, yi<0y_i<0) you are drawn in. Scaling fields are precisely the coordinates aligned with these eigen-streamlines, in which the otherwise tangled flow separates into independent expansions and contractions. This has a startling consequence, and it is the whole reason universality works. Start from two utterly different microscopic Hamiltonians. Each is a different point in coupling space, with different admixtures of every irrelevant operator. But those irrelevant components all die under the flow, so both trajectories are funnelled onto the same low-dimensional relevant subspace and end up governed by the same handful of numbers. A generic starting point, projected onto the relevant axes, is all the fixed point ever remembers. The multiplicative law is forced, not chosen. Coarse-graining by b1b_1 and then by b2b_2 must equal coarse-graining by b1b2b_1 b_2 — the RG is a semigroup. If a quantity simply rescales, the only function consistent with λ(b1)λ(b2)=λ(b1b2)\lambda(b_1)\lambda(b_2)=\lambda(b_1 b_2) is a power law λ(b)=by\lambda(b)=b^{y}. So the existence of scaling fields and the existence of power-law critical singularities are one and the same fact wearing two hats. Why nonlinear? A bare coupling like the temperature distance tt picks up a small admixture of other operators after coarse-graining; the combination that scales cleanly is a Taylor series ut=t+c1t2+c2h2+⋯u_t=t+c_1 t^2+c_2 h^2+\cdots. Keeping only the linear term is enough for the leading exponents, but the curvature terms are what generate the universal correction-to-scaling exponent ω=∣yirr∣\omega=|y_{\text{irr}}| and the non-universal amplitude ratios — the difference between a textbook power law and real data.

Figure
Flow near the fixed point: relevant direction blows up, irrelevant collapses in
u_t (relevant, y_t > 0) u_3 (irrelevant, y_3 < 0) fixed point K* generic start → projected onto u_t the irrelevant component dies; only the relevant survives at large b

The figure above is the whole story of universality in one picture: any generic starting point (amber) is dragged down the irrelevant axis and flung out along the relevant one, so at large rescaling bb every trajectory looks the same — a pure motion along utu_t. That is why microscopic detail washes out.

Figure
Bare parameter versus scaling field: the nonlinear curved coordinate
t = (T−T_c)/T_c u_t u_t ≈ t (linear) u_t = t + c₁t² + c₂h² + … curvature ↔ corrections to scaling (ω) agree at the critical point

The core equations

The renormalization group is a map on coupling space that must compose with itself: doing two coarse-grainings in a row is one bigger coarse-graining. This single algebraic fact — the semigroup property — is what makes everything downstream a power law. We record the map and its composition rule first.

K′=Rb(K),Rb1∘Rb2=Rb1b2(1)\mathbf{K}' = R_b(\mathbf{K}),\qquad R_{b_1}\circ R_{b_2}=R_{b_1 b_2} \qquad (1)

where: K\mathbf{K} — a point in coupling space, the full set of Hamiltonian parameters; RbR_b — one coarse-graining plus rescaling by linear factor b>1b>1; bb — the spatial rescaling factor, often written b=eℓb=e^{\ell} for a continuous flow. Read the composition law as the definition of what "the renormalization group" means: the group of magnifications, closed under repetition. To learn the flow near a fixed point K∗\mathbf{K}^* we linearize. Writing the deviation δK=K−K∗\delta\mathbf{K}=\mathbf{K}-\mathbf{K}^*, the map becomes a matrix multiplication by the Jacobian M(b)\mathsf{M}(b) evaluated at the fixed point. This matrix carries all the local information about stability.

δK′=M(b) δK,Mij=∂Ki′∂Kj∣K∗(2)\delta\mathbf{K}' = \mathsf{M}(b)\,\delta\mathbf{K},\qquad \mathsf{M}_{ij}=\left.\frac{\partial K_i'}{\partial K_j}\right|_{\mathbf{K}^*} \qquad (2)

where: δK\delta\mathbf{K} — deviation from the fixed point; M(b)\mathsf{M}(b) — the stability matrix (Jacobian) of the RG map at K∗\mathbf{K}^*. Its eigenvalues decide which perturbations grow and which decay; its eigenvectors, projected onto δK\delta\mathbf{K}, are the scaling fields. Diagonalizing M(b)\mathsf{M}(b) and projecting the flow onto its left eigenvectors defines the scaling fields at linear order. The eigenvalue attached to each is forced by composition to be a pure power of bb, and that exponent yiy_i is the object we care about.

ui=∑jϕj(i) δKj+O(δK2),ui′=byiui(3)u_i = \sum_j \phi^{(i)}_j\,\delta K_j + O(\delta K^2),\qquad u_i' = b^{y_i} u_i \qquad (3)

where: ϕ(i)\phi^{(i)} — the left eigenvector of M(b)\mathsf{M}(b) with eigenvalue byib^{y_i}; uiu_i — the ii-th scaling field; yiy_i — its RG eigenvalue exponent. A positive yiy_i means the field is relevant (its magnitude grows under coarse-graining); negative means irrelevant; zero means marginal, the delicate case that produces logarithms. Linear scaling fields solve the leading problem, but the exact map has quadratic and higher terms. Wegner's construction promotes each uiu_i to an analytic-but-nonlinear function of the couplings, chosen so that the multiplicative law holds to all orders — this is what makes exponents clean and pushes non-universal junk into amplitudes.

ui(K)=δKi(eig)+∑j,kajk(i) δKj δKk+⋯ ,ui→byiui exactly(4)u_i(\mathbf{K}) = \delta K_i^{(\text{eig})} + \sum_{j,k} a_{jk}^{(i)}\,\delta K_j\,\delta K_k+\cdots,\qquad u_i\to b^{y_i}u_i\ \text{exactly} \qquad (4)

where: ajk(i)a^{(i)}_{jk} — coefficients fixed by demanding exact multiplicative scaling beyond linear order. They can always be solved for unless a resonance yi=yj+yky_i=y_j+y_k occurs, in which case the smooth redefinition fails and logarithmic terms appear — exactly the marginal case at the upper critical dimension. With the scaling fields in hand, the singular part of the free-energy density inherits a generalized homogeneity: coarse-graining multiplies it by b−db^{-d} (the volume per remaining degree of freedom shrinks) while driving each field by byib^{y_i}. This is the master relation everything else follows from.

fs(u1,u2,… )=b−d fs ⁣(by1u1, by2u2,… )(5)f_s(u_1,u_2,\dots) = b^{-d}\,f_s\!\left(b^{y_1}u_1,\,b^{y_2}u_2,\dots\right) \qquad (5)

where: fsf_s — the singular part of the free energy per unit volume (per spin), in units of kBTk_BT; dd — the spatial dimension, entering because fsf_s is a density. Because this holds for any bb, it is a functional constraint strong enough to fix all exponents. Since bb is arbitrary, we may choose it. Pinning the leading relevant thermal field to unity, byt∣ut∣=1b^{y_t}|u_t|=1, eliminates bb and collapses fsf_s onto a function of a single dimensionless ratio — the Widom scaling form, and the reason data collapse works.

fs(ut,uh)=∣ut∣d/yt Φ ⁣(uh∣ut∣yh/yt)(6)f_s(u_t,u_h) = |u_t|^{d/y_t}\,\Phi\!\left(\frac{u_h}{|u_t|^{y_h/y_t}}\right) \qquad (6)

where: Φ\Phi — the universal scaling function (with two branches for ut≷0u_t\gtrless0); ut,uhu_t,u_h — the thermal and magnetic scaling fields. The gap exponent Δ=yh/yt\Delta=y_h/y_t in the argument governs how the field scale tracks the temperature. Differentiating the collapsed form gives every static exponent as a rational function of just two eigenvalues and the dimension. This is the concrete sense in which two numbers control all of critical thermodynamics.

α=2−dyt,β=d−yhyt,γ=2yh−dyt,δ=yhd−yh,ν=1yt,η=d+2−2yh(7)\alpha=2-\frac{d}{y_t},\quad \beta=\frac{d-y_h}{y_t},\quad \gamma=\frac{2y_h-d}{y_t},\quad \delta=\frac{y_h}{d-y_h},\quad \nu=\frac1{y_t},\quad \eta=d+2-2y_h \qquad (7)

where: yt,yhy_t,y_h — the thermal and magnetic RG eigenvalues; α…η\alpha\dots\eta — the standard critical exponents. Feed in yt=1, yh=15/8, d=2y_t=1,\,y_h=15/8,\,d=2 and out come the exact 2D-Ising values; feed in the Gaussian yt=2, yh=3, d=4y_t=2,\,y_h=3,\,d=4 and out come mean-field values. Finally, the leading irrelevant field does not simply vanish — it decays but leaves a trail. Keeping it produces the universal corrections to scaling that dominate any real fit at finite distance from TcT_c.

fs=∣ut∣2−α[Φ0(uh/∣ut∣Δ)+∣ut∣ νωΦ1(⋯ )+⋯](8)f_s = |u_t|^{2-\alpha}\Big[\Phi_0\big(u_h/|u_t|^{\Delta}\big)+|u_t|^{\,\nu\omega}\Phi_1(\cdots)+\cdots\Big] \qquad (8)

where: ω=∣y3∣\omega=|y_3| — the magnitude of the leading irrelevant eigenvalue; θ=νω\theta=\nu\omega — the Wegner correction-to-scaling exponent; Φ1\Phi_1 — the subleading scaling function. Asymptotically negligible, these terms are indispensable for extracting exponents from data that never quite reaches TcT_c.

Derivation AFrom the semigroup law to power-law eigenvalues and scaling fields

We prove rigorously that the linearized RG eigenvalues must be pure powers byib^{y_i} and construct the scaling fields as the matching eigen-coordinates. The engine is a single algebraic fact — the RG is a semigroup — and the goal is to show that this fact alone forces power-law scaling, so that the existence of well-defined critical exponents is not an assumption but a theorem. Start from the defining structure of the RG map, linearize about a fixed point, and let composition do the work.

Step 1. Expand the map about the fixed point: Rb(K∗+δK)=K∗+M(b) δK+O(δK2)R_b(\mathbf{K}^*+\delta\mathbf{K})=\mathbf{K}^*+\mathsf{M}(b)\,\delta\mathbf{K}+O(\delta K^2). Why: Rb(K∗)=K∗R_b(\mathbf{K}^*)=\mathbf{K}^* by definition of a fixed point, so the constant term drops and the leading behaviour is linear in the deviation.

Step 2. Apply the chain rule to Rb1∘Rb2=Rb1b2R_{b_1}\circ R_{b_2}=R_{b_1 b_2} at the fixed point: M(b1) M(b2)=M(b1b2)\mathsf{M}(b_1)\,\mathsf{M}(b_2)=\mathsf{M}(b_1 b_2). Why: The Jacobian of a composition is the product of Jacobians; evaluated at the fixed point (which every RbR_b maps to itself) this transfers the semigroup law from the maps to their stability matrices.

Step 3. Diagonalize using left eigenvectors: ∑jϕj(i)Mjk(b)=λi(b) ϕk(i)\sum_j \phi^{(i)}_j \mathsf{M}_{jk}(b) = \lambda_i(b)\,\phi^{(i)}_k, so each eigenvalue inherits λi(b1b2)=λi(b1)λi(b2)\lambda_i(b_1 b_2)=\lambda_i(b_1)\lambda_i(b_2). Why: M\mathsf{M} is generically non-symmetric, so we use its left eigenvectors; on the common eigenbasis the matrix product law passes to each scalar eigenvalue.

Step 4. Solve the multiplicative functional equation: the only continuous solution of λ(b1b2)=λ(b1)λ(b2)\lambda(b_1 b_2)=\lambda(b_1)\lambda(b_2) is λi(b)=byi\lambda_i(b)=b^{y_i} with yi∈Ry_i\in\mathbb{R}. Why: Taking logarithms turns the multiplicative law into Cauchy's additive equation g(b1b2)=g(b1)+g(b2)g(b_1 b_2)=g(b_1)+g(b_2) in ln⁡b\ln b, whose continuous solutions are linear — hence λ\lambda is a power. This is where the power law is born.

Step 5. Define the scaling fields ui≡∑jϕj(i)δKju_i\equiv\sum_j\phi^{(i)}_j\delta K_j and confirm ui′=∑jϕj(i)δKj′=λi(b) ui=byiuiu_i'=\sum_j\phi^{(i)}_j\delta K_j'=\lambda_i(b)\,u_i=b^{y_i}u_i. Why: Projecting the linear flow δK′=M δK\delta\mathbf{K}'=\mathsf{M}\,\delta\mathbf{K} onto a left eigenvector diagonalizes it; the scaling field is just that projection, and it scales by its own eigenvalue. Verified: verify_aq-sm-scaling-fields-scaling-variables.py checks the multiplicative functional equation numerically (a random yy satisfies b1yb2y=(b1b2)yb_1^{y}b_2^{y}=(b_1b_2)^{y}) and confirms that the log-map turns it into Cauchy additivity; matches Kardar Ch. 4 and Goldenfeld Ch. 9.

Go deeper · Derivation B · 1 minThe exponent relations from the scaling form
Derivation BThe exponent relations from the scaling form

We differentiate the homogeneous free energy to extract α,β,γ,δ\alpha,\beta,\gamma,\delta and then verify that the Rushbrooke and Widom identities hold as exact equalities — a clean demonstration that the two relevant eigenvalues encode all of thermodynamics. This is the structural result the article owns: not the individual exponents (which the RG-eigenvalue article supplies) but the fact that they are forced to satisfy the scaling laws with equality. Begin from fs(ut,uh)=b−dfs(bytut,byhuh)f_s(u_t,u_h)=b^{-d}f_s(b^{y_t}u_t,b^{y_h}u_h) with ut∝t, uh∝hu_t\propto t,\ u_h\propto h at leading order.

Step 1. Specific heat at h=0h=0: setting byt∣ut∣=1b^{y_t}|u_t|=1 gives fs∼∣t∣d/ytf_s\sim|t|^{d/y_t}, and two temperature derivatives yield C∼∣t∣d/yt−2C\sim|t|^{d/y_t-2}, so α=2−d/yt\alpha=2-d/y_t. Why: Each ∂t\partial_t lowers the power of ∣t∣|t| by one; the specific heat is the second temperature derivative of the free energy.

Step 2. Order parameter: m=−∂hfs∣h=0∼∣t∣(d−yh)/ytm=-\partial_h f_s|_{h=0}\sim|t|^{(d-y_h)/y_t}, so β=(d−yh)/yt\beta=(d-y_h)/y_t. Why: One field derivative brings down a factor byh=∣t∣−yh/ytb^{y_h}=|t|^{-y_h/y_t}; evaluating at h=0h=0, t<0t<0 gives the spontaneous magnetization.

Step 3. Susceptibility: a second field derivative gives χ∼∣t∣(d−2yh)/yt\chi\sim|t|^{(d-2y_h)/y_t}, so γ=(2yh−d)/yt\gamma=(2y_h-d)/y_t. Why: Two field derivatives bring b2yhb^{2y_h}; the sign flips because γ\gamma is defined by χ∼∣t∣−γ\chi\sim|t|^{-\gamma}.

Step 4. Critical isotherm: at t=0t=0 rescale instead with byh∣uh∣=1b^{y_h}|u_h|=1 to get m∼h(d−yh)/yhm\sim h^{(d-y_h)/y_h}, so δ=yh/(d−yh)\delta=y_h/(d-y_h). Why: On the critical isotherm there is no temperature scale, so the magnetic field must set the magnification.

Step 5. Combine into Rushbrooke: α+2β+γ=(2yt−d)+2(d−yh)+(2yh−d)yt=2ytyt=2\alpha+2\beta+\gamma=\dfrac{(2y_t-d)+2(d-y_h)+(2y_h-d)}{y_t}=\dfrac{2y_t}{y_t}=2. Why: Every exponent is a rational function of yt,yh,dy_t,y_h,d; the numerators telescope so the dimension and yhy_h cancel, leaving exactly 22. The old inequality becomes an identity.

Step 6. Combine into Widom: β(δ−1)=d−yhyt⋅2yh−dd−yh=2yh−dyt=γ\beta(\delta-1)=\dfrac{d-y_h}{y_t}\cdot\dfrac{2y_h-d}{d-y_h}=\dfrac{2y_h-d}{y_t}=\gamma. Why: The (d−yh)(d-y_h) factors cancel, leaving γ\gamma exactly; only yty_t and yhy_h (two numbers) are independent. Verified: verify_aq-sm-scaling-fields-scaling-variables.py evaluates all six exponents symbolically in (yt,yh,d)(y_t,y_h,d) and asserts Rushbrooke, Widom, Fisher, and Josephson as identities, then substitutes the 2D-Ising eigenvalues yt=1,yh=15/8,d=2y_t=1,y_h=15/8,d=2 to reproduce (α,β,γ,δ,ν,η)=(0,1/8,7/4,15,1,1/4)(\alpha,\beta,\gamma,\delta,\nu,\eta)=(0,1/8,7/4,15,1,1/4).

Go deeper · Derivation C · 1 minNonlinear scaling fields and the origin of corrections to scaling
Derivation CNonlinear scaling fields and the origin of corrections to scaling

We construct the analytic-but-nonlinear redefinition that makes the scaling law exact, and show how the leading irrelevant field generates the Wegner correction-to-scaling exponent. The point is conceptual: measured exponents are clean only because the true scaling variable is a curved coordinate, and the leftover curvature is itself universal. Demand ui(RbK)=byiui(K)u_i(R_b\mathbf{K})=b^{y_i}u_i(\mathbf{K}) beyond linear order and see what it forces.

Step 1. Write the most general analytic forms respecting the Ising symmetry h→−hh\to-h: ut=t+c1t2+c2h2+⋯u_t=t+c_1 t^2+c_2 h^2+\cdots, uh=h+c3 t h+⋯u_h=h+c_3\,t\,h+\cdots. Why: Symmetry forbids a term linear in hh inside utu_t (which is even) and forbids a term independent of hh inside uhu_h (which is odd); the leading corrections are the ones shown.

Step 2. Impose exact scaling order by order: the quadratic coefficients solve ajk(i)(byj+yk−byi)=(Rb)jk(2)a^{(i)}_{jk}(b^{y_j+y_k}-b^{y_i})=(R_b)^{(2)}_{jk}. Why: Matching the byiuib^{y_i}u_i law at quadratic order gives a linear equation for each coefficient; it is solvable as long as the prefactor byj+yk−byib^{y_j+y_k}-b^{y_i} does not vanish (no resonance).

Step 3. Keep one irrelevant field u3u_3 (eigenvalue y3<0y_3<0): fs(ut,uh,u3)=b−dfs(bytut,byhuh,b−∣y3∣u3)f_s(u_t,u_h,u_3)=b^{-d}f_s(b^{y_t}u_t,b^{y_h}u_h,b^{-|y_3|}u_3). Why: The irrelevant direction does not vanish in a single step of the iterated map; it decays but must be carried along to see its effect.

Step 4. Fix b=∣ut∣−1/ytb=|u_t|^{-1/y_t}, so the third argument becomes u3∣ut∣∣y3∣/yt→0u_3|u_t|^{|y_3|/y_t}\to0 as t→0t\to0. Why: This is the same magnification that collapses the free energy; under it the irrelevant field's scaled value is a small parameter, formalizing the word "irrelevant."

Step 5. Expand in the small argument: fs=∣ut∣2−α[Φ0(z)+∣ut∣νωΦ1(z)+⋯ ]f_s=|u_t|^{2-\alpha}[\Phi_0(z)+|u_t|^{\nu\omega}\Phi_1(z)+\cdots] with ω=∣y3∣\omega=|y_3| and θ=νω\theta=\nu\omega. Why: Taylor-expanding the scaling function in its small third argument produces a leading correction of relative size ∣ut∣νω|u_t|^{\nu\omega}; this is the Wegner correction, universal in its exponent though not its amplitude. Verified: verify_aq-sm-scaling-fields-scaling-variables.py confirms θ=νω\theta=\nu\omega numerically for the 3D-Ising values ν≈0.630, ω≈0.83\nu\approx0.630,\ \omega\approx0.83 giving θ≈0.52\theta\approx0.52, and checks the resonance condition yi=yj+yky_i=y_j+y_k that signals logarithms; structure follows Wegner (1974) and Cardy Ch. 3.

Limiting cases

CaseResultMeaning
Leading relevant thermal fieldut→bytut,  ν=1/ytu_t \to b^{y_t} u_t,\ \ \nu = 1/y_tControls the divergence of the correlation length ξ∼∣ut∣−ν\xi\sim\lvert u_t\rvert^{-\nu}; ut∝tu_t\propto t at leading order.
Leading relevant magnetic fielduh→byhuh,  Δ=βδ=νyhu_h \to b^{y_h} u_h,\ \ \Delta=\beta\delta=\nu y_hConjugate to the order parameter; uh∝hu_h\propto h at leading order, with yh=(d+2−η)/2y_h=(d+2-\eta)/2.
Leading irrelevant fieldu3→b−∣y3∣u3,  ω=∣y3∣u_3 \to b^{-\lvert y_3\rvert} u_3,\ \ \omega=\lvert y_3\rvertGives corrections to scaling ∼ut νω\sim u_t^{\,\nu\omega} (Wegner θ=νω\theta=\nu\omega); ignorable asymptotically but vital for fits.
Marginal field (y=0y=0)u′=u+a u2ln⁡b+⋯u' = u + a\,u^2\ln b + \cdotsPure power law fails; logarithmic corrections appear, as at the upper critical dimension d=4d=4 for ϕ4\phi^4.
Dangerously irrelevant fieldfsf_s singular as uu→0u_u\to0An irrelevant field on which fsf_s depends singularly; restores mean-field exponent relations above dc=4d_c=4.

Worked examples

Example 1 — The 2D-Ising exponents from two eigenvalues. The RG eigenvalues at the square-lattice Ising fixed point are yt=1y_t=1 and yh=15/8y_h=15/8 in d=2d=2. Compute the full exponent set and check it against Onsager's exact solution.

Step 1. Correlation-length exponent: ν=1/yt=1/1=1\nu=1/y_t=1/1=\mathbf{1}.

Step 2. Specific heat: α=2−d/yt=2−2/1=0\alpha=2-d/y_t=2-2/1=\mathbf{0} (a logarithmic divergence, the marginal α=0\alpha=0 case).

Step 3. Order parameter: β=(d−yh)/yt=(2−15/8)/1=1/8=0.125\beta=(d-y_h)/y_t=(2-15/8)/1=1/8=\mathbf{0.125}.

Step 4. Susceptibility: γ=(2yh−d)/yt=(2⋅15/8−2)/1=(15/4−2)=7/4=1.75\gamma=(2y_h-d)/y_t=(2\cdot15/8-2)/1=(15/4-2)=7/4=\mathbf{1.75}.

Step 5. Critical isotherm: δ=yh/(d−yh)=(15/8)/(2−15/8)=(15/8)/(1/8)=15\delta=y_h/(d-y_h)=(15/8)/(2-15/8)=(15/8)/(1/8)=\mathbf{15}.

Step 6. Anomalous dimension: η=d+2−2yh=2+2−2⋅15/8=4−15/4=1/4=0.25\eta=d+2-2y_h=2+2-2\cdot15/8=4-15/4=1/4=\mathbf{0.25}. The set (α,β,γ,δ,ν,η)=(0,18,74,15,1,14)(\alpha,\beta,\gamma,\delta,\nu,\eta)=(0,\tfrac18,\tfrac74,15,1,\tfrac14) is exactly Onsager's 2D-Ising result — obtained here purely from the two relevant eigenvalues.

Example 2 — Only two of the six are independent. Using the same values, verify the four scaling relations hold as exact equalities, confirming that yt,yh,dy_t,y_h,d carry all the information.

Step 1. Rushbrooke: α+2β+γ=0+2⋅18+74=14+74=84=2\alpha+2\beta+\gamma=0+2\cdot\tfrac18+\tfrac74=\tfrac14+\tfrac74=\tfrac{8}{4}=\mathbf{2}. ✓

Step 2. Widom: γ=β(δ−1)=18 (15−1)=148=74=1.75\gamma=\beta(\delta-1)=\tfrac18\,(15-1)=\tfrac{14}{8}=\tfrac74=\mathbf{1.75}. ✓

Step 3. Fisher: γ=ν(2−η)=1⋅(2−14)=74=1.75\gamma=\nu(2-\eta)=1\cdot(2-\tfrac14)=\tfrac74=\mathbf{1.75}. ✓

Step 4. Josephson (hyperscaling): 2−α=dν⇒2−0=2⋅1=22-\alpha=d\nu\Rightarrow 2-0=2\cdot1=\mathbf{2}. ✓ All four hold with equality. The RG has promoted the classical thermodynamic inequalities to identities — the signature that a fixed point, not a free assumption, is in charge.

Origins · 1 minrelevant operator and how this was discovered
RO
relevant operator
RG classification of perturbations · The organizing idea that a critical fixed point is characterized by which perturbations grow under coarse-graining. A perturbation with positive RG eigenvalue is "relevant" — it drives the system away from criticality and must be tuned to zero to sit at the transition (temperature and field for a magnet). The count of relevant operators equals the number of parameters that must be tuned, which is why a generic critical point needs exactly two knobs. Formalized in the early 1970s in Wilson's and Wegner's work and now the backbone of how physicists classify every continuous transition.

History
1966
— Kadanoff introduces the block-spin scaling picture, giving the homogeneity hypothesis a physical basis
1971
— Wilson formulates the renormalization group as a flow in coupling space, defining fixed points and eigenvalue exponents
1972
— Wegner introduces nonlinear scaling fields and the systematic theory of corrections to scaling
1974
— Wegner's "corrections to scaling" analysis fixes the Wegner exponent θ = νω from the leading irrelevant eigenvalue

Going further

References. Pathria & Beale, Statistical Mechanics, 3rd ed., Ch. 14 (RG and scaling) · Kardar, Statistical Physics of Fields, Ch. 4–5 · Goldenfeld, Lectures on Phase Transitions and the Renormalization Group, Ch. 9–10 · Wegner, J. Phys. C 7, 2098 (1974) (corrections to scaling, nonlinear scaling fields) · Cardy, Scaling and Renormalization in Statistical Physics, Ch. 3 · Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity, Ch. 12

Provenance. All equations are standard and recomputed by verify_aq-sm-scaling-fields-scaling-variables.py: the multiplicative functional equation λ(b1b2)=λ(b1)λ(b2)⇒λ=by\lambda(b_1b_2)=\lambda(b_1)\lambda(b_2)\Rightarrow\lambda=b^y; the scaling-form derivation of α,β,γ,δ,ν,η\alpha,\beta,\gamma,\delta,\nu,\eta from fs=b−dfs(bytut,byhuh)f_s=b^{-d}f_s(b^{y_t}u_t,b^{y_h}u_h); the Rushbrooke, Widom, Fisher, and Josephson identities as algebraic consequences of yt,yh,dy_t,y_h,d; and the Wegner identification θ=νω\theta=\nu\omega. Numerical anchors: 2D-Ising yt=1, yh=15/8y_t=1,\ y_h=15/8 give the exact (α,β,γ,δ,ν,η)=(0,18,74,15,1,14)(\alpha,\beta,\gamma,\delta,\nu,\eta)=(0,\tfrac18,\tfrac74,15,1,\tfrac14), matching Onsager's Solution of the 2D Ising Model; 3D-Ising ω≈0.83\omega\approx0.83 from high-precision methods. Values quoted are textbook-standard; no non-standard figures introduced.

Builds on: RG Flow in Coupling Space · Fixed Points of the RG · Linearization About a Fixed Point and Eigenvalue Spectrum · The Renormalization Group Idea: Coarse-Graining and Rescaling. Leads to: Relevant, Irrelevant, and Marginal Operators · Critical Exponents from the RG · Wilson-Fisher Fixed Point · The Epsilon-Expansion.

Are scaling fields the same as eigenoperators of the RG?

Closely related but not identical. The eigenoperators are the perturbations added to the fixed-point Hamiltonian whose couplings scale multiplicatively; the scaling fields are the couplings themselves, i.e. the coordinates in coupling space along the eigen-directions. In practice people use "scaling field," "scaling variable," and "eigenoperator coupling" almost interchangeably, but the field is the number you dial and the operator is the term it multiplies in the Hamiltonian.

Why exactly two relevant fields for the Ising class?

Because the Ising fixed point has exactly two perturbations with positive RG eigenvalue: the thermal (even) one with yt=1/νy_t=1/\nu and the magnetic (odd) one with yhy_h. Every other symmetry-allowed operator is irrelevant. The count of relevant fields equals the number of parameters you must tune to reach the critical point, which is why criticality generically needs two knobs, temperature and field. A different symmetry (say a tricritical point) has more relevant fields and needs more tuning.

If the coefficients relating u_t to t are non-universal, how can exponents be universal?

The universality lives in the eigenvalues yiy_i, which are properties of the fixed point alone. The coefficients cic_i in ut=t+c1t2+⋯u_t=t+c_1t^2+\cdots are metric factors that translate from lab units to fixed-point units; they set amplitudes and correction sizes but never the exponents. Two materials in the same class have identical yt,yhy_t,y_h but different cic_i, so their exponents agree while their amplitudes and the location of TcT_c differ.

What is the difference between a marginal field and an irrelevant one in practice?

An irrelevant field (y<0y<0) decays as a power under coarse-graining and produces corrections to scaling that vanish algebraically near TcT_c. A marginal field (y=0y=0) neither grows nor decays at linear order; its fate is decided by the (usually quadratic) next term in the flow. Marginal fields produce logarithmic rather than power-law corrections — the famous ∣ln⁡t∣|\ln t| factors at the upper critical dimension d=4d=4, where the ϕ4\phi^4 coupling is marginally irrelevant.

How does data collapse connect to the scaling-field homogeneity?

Data collapse is the experimental face of equation (6). Because fs=∣ut∣d/ytΦ(uh/∣ut∣Δ)f_s=|u_t|^{d/y_t}\Phi(u_h/|u_t|^{\Delta}), plotting a rescaled observable (like m/∣t∣βm/|t|^\beta) against the rescaled field (h/∣t∣Δh/|t|^{\Delta}) for many temperatures makes every curve fall on the single function Φ\Phi. A clean collapse confirms the homogeneity and measures the exponents; a spread that shrinks toward TcT_c is the fingerprint of the leading irrelevant field. The dedicated data-collapse article develops this into a quantitative method.

Why does the free energy pick up exactly b^{-d} and not some other power?

The free energy is extensive: it scales with the number of degrees of freedom. One coarse-graining step merges bdb^d sites into one block, so the free energy per unit volume is multiplied by bdb^d when you rewrite the theory on the coarse lattice, i.e. fs(u)=b−dfs(u′)f_s(u)=b^{-d}f_s(u'). The power is exactly the spatial dimension because that is how many directions get compressed per step. This is also why hyperscaling relations, which contain dd, are the fragile ones.

Before this

Where this leads