Skip to content

Dimensional Transmutation and Mass Gaps

The previous page ended with the one-loop asymptotically free flow

β(g)=−b016π2g3+O(g5),b0>0,\beta(g)=-{b_0\over16\pi^2}g^3+O(g^5), \qquad b_0>0,

and therefore

1g2(μ)=1g2(μ0)+b08π2log⁡μμ0.{1\over g^2(\mu)}={1\over g^2(\mu_0)}+{b_0\over8\pi^2}\log{\mu\over\mu_0}.

This formula describes dimensional transmutation: instead of specifying a dimensionless coupling at an arbitrary subtraction scale, one may specify an RG-invariant scale such as ΛYM\Lambda_{\mathrm{YM}} or ΛQCD\Lambda_{\mathrm{QCD}}. The numerical value of this scale depends on the renormalization scheme; a measurable mass requires additional dynamical information. This lesson derives the scale, explains its continuum-limit meaning, and compares it with the exponential scale in a specified weak-coupling pairing model.

The same idea also explains why perturbation theory can be excellent at short distances while giving no direct perturbative expansion for the masses of bound states. If the generated scale is

Λ∼μexp⁡[−8π2b0g2(μ)],\Lambda\sim \mu\exp\left[-{8\pi^2\over b_0g^2(\mu)}\right],

then it is non-analytic at g=0g=0. No finite Taylor series in gg can see it. This is the first warning that mass gaps in asymptotically free theories are nonperturbative, even though their scale is already visible in perturbative RG.

Required background. QED and Yang–Mills Beta Functions fixes the beta-function normalization, derives asymptotic freedom, and supplies the one-loop running used throughout this page.

Scale variables. Momentum scales such as μ\mu and pp denote positive Euclidean magnitudes. On this page

β(g)=μdgdμ,β(g)=−b016π2g3+O(g5)\beta(g)=\mu {dg\over d\mu}, \qquad \beta(g)=-{b_0\over16\pi^2}g^3+O(g^5)

for an asymptotically free coupling. The UV cutoff is written as Λ0=a−1\Lambda_0=a^{-1}, where aa may be a lattice spacing or any short-distance regulator length. When the precise one-loop coefficient is unimportant, we write positive constants as cc.

Start from the one-loop equation

ddlog⁡μ1g2(μ)=b08π2.{d\over d\log\mu}{1\over g^2(\mu)}={b_0\over8\pi^2}.

Integrating gives

1g2(μ)=b08π2log⁡μΛ,{1\over g^2(\mu)}={b_0\over8\pi^2}\log{\mu\over\Lambda},

where the integration constant has been written as a scale Λ\Lambda. Solving for g2(μ)g^2(\mu) gives

g2(μ)=8π2b0log⁡(μ/Λ)(μ≫Λ).\boxed{ g^2(\mu)={8\pi^2\over b_0\log(\mu/\Lambda)} } \qquad (\mu\gg\Lambda).

Equivalently,

Λ=μexp⁡[−8π2b0g2(μ)].\boxed{ \Lambda=\mu\exp\left[-{8\pi^2\over b_0g^2(\mu)}\right]. }

This is the RG-invariant scale. At one loop it is exactly independent of μ\mu. Beyond one loop its definition is convention-dependent by a multiplicative constant, but the statement that the theory generates a scale is not convention-dependent. Changing the renormalization scheme changes the numerical value assigned to Λ\Lambda; it does not change dimensionless physical predictions such as mass ratios.

The all-orders construction makes this precise. In a specified renormalization scheme S\mathcal S, define formally

ΛS=μexp⁡[−∫g(μ)dg′βS(g′)]\boxed{ \Lambda_{\mathcal S} =\mu\exp\left[-\int^{g(\mu)}{dg'\over\beta_{\mathcal S}(g')}\right] }

with the finite part of the integration constant included in the definition of the scheme. Differentiating shows dΛS/dlog⁡μ=0d\Lambda_{\mathcal S}/d\log\mu=0. If

β(g)=−b0g316π2−b1g5(16π2)2+⋯ ,\beta(g)=-{b_0g^3\over16\pi^2} -{b_1g^5\over(16\pi^2)^2}+\cdots,

then at weak coupling

ΛS=μexp⁡[−8π2b0g2(μ)]×(b0g2(μ)16π2)−b1/(2b02)[1+O(g2)].\begin{aligned} \Lambda_{\mathcal S}={}& \mu\exp\left[-{8\pi^2\over b_0g^2(\mu)}\right]\\ &\times \left({b_0g^2(\mu)\over16\pi^2}\right)^{-b_1/(2b_0^2)} \left[1+O(g^2)\right]. \end{aligned}

A redefinition g2↦g2+c1g4+⋯g^2\mapsto g^2+c_1g^4+\cdots changes the numerical value called ΛS\Lambda_{\mathcal S} by a finite multiplicative factor. It cannot change dimensionless physical predictions such as mass ratios Mi/MjM_i/M_j, nor can it remove the essential singularity e−const/g2e^{-\mathrm{const}/g^2}.

The one-loop asymptotically free running coupling defines an RG-invariant scale

The one-loop asymptotically free flow trades a dimensionless coupling for a dimensionful scale. Its straight inverse-coupling trajectory extrapolates to zero near μ=Λ\mu=\Lambda; the extrapolation is not a controlled infrared prediction. With a UV cutoff a−1a^{-1}, the RG scale behaves as aΛ∼e−const/g02a\Lambda\sim e^{-\mathrm{const}/g_0^2}. Original schematic of the one-loop approximation.

The name “dimensional transmutation” is literal. For this Euclidean cutoff illustration, write the pure Yang–Mills action in terms of a rescaled connection A=g0A\mathcal A=g_0A and curvature Fμνa=∂μAνa−∂νAμa+fabcAμbAνc\mathcal F^a_{\mu\nu}=\partial_\mu\mathcal A^a_\nu-\partial_\nu\mathcal A^a_\mu+f^{abc}\mathcal A^b_\mu\mathcal A^c_\nu. With the topological angle fixed to zero, the action is

SE=14g02∫d4x FμνaFμνaS_E={1\over4g_0^2}\int d^4x\,\mathcal F^a_{\mu\nu}\mathcal F^a_{\mu\nu}

at the cutoff scale Λ0=a−1\Lambda_0=a^{-1}. Classically, g0g_0 is dimensionless and there is no mass parameter. Quantum mechanically, the pair (a,g0)(a,g_0) determines

ΛYM∼1aexp⁡[−8π2b0g02]\boxed{ \Lambda_{\mathrm{YM}}\sim {1\over a}\exp\left[-{8\pi^2\over b_0g_0^2}\right] }

up to scheme-dependent multiplicative factors and higher-loop corrections. To keep ΛYM\Lambda_{\mathrm{YM}} finite while taking a→0a\to0, one must tune

g02(a)≃8π2b0log⁡(1/aΛYM).g_0^2(a)\simeq {8\pi^2\over b_0\log(1/a\Lambda_{\mathrm{YM}})}.

Thus the continuum limit is not obtained by keeping the bare coupling fixed. It is obtained by approaching the UV fixed point g0=0g_0=0 along a trajectory that keeps the physical scale fixed.

The same result can be written in the language of a running coupling measured at momentum pp below the cutoff. At one loop,

1g2(p)=1g02+b08π2log⁡pΛ0,Λ0=1a.{1\over g^2(p)}={1\over g_0^2}+{b_0\over8\pi^2}\log{p\over\Lambda_0}, \qquad \Lambda_0={1\over a}.

Since p<Λ0p<\Lambda_0, the logarithm is negative. As pp is lowered, 1/g2(p)1/g^2(p) decreases and the coupling grows. The formal scale at which the one-loop denominator vanishes is

Λ=Λ0exp⁡[−8π2b0g02].\Lambda=\Lambda_0\exp\left[-{8\pi^2\over b_0g_0^2}\right].

The zero of this denominator is a one-loop extrapolation. It cannot establish that the exact coupling diverges. If the flow continues toward strong coupling, perturbation theory fails before that formal pole. More generally, asymptotic freedom alone does not determine the infrared endpoint: higher-order terms can instead produce an infrared fixed point. For example, the beta function β(g)=−bg3(1−g2/g∗2)\beta(g)=-b g^3(1-g^2/g_*^2) with b>0b>0 and small g∗>0g_*>0 is asymptotically free but approaches g∗g_* from below as the scale is lowered. Its leading term alone misses that endpoint.

For p≫Λp\gg\Lambda,

g2(p)=8π2b0log⁡(p/Λ)\boxed{ g^2(p)={8\pi^2\over b_0\log(p/\Lambda)} }

is small when the logarithm is sufficiently large. Near the apparent one-loop pole, this formula supplies no controlled value for the coupling. For pure Yang–Mills theory the expected infrared regime is strongly coupled; establishing its spectrum requires nonperturbative physics.

This is why Yang–Mills theory and QCD can be weakly coupled at short distance and strongly coupled at long distance. Short-distance gauge-invariant observables admit a weakly coupled partonic description. Confinement, bound states, and a mass gap concern the infrared and require nonperturbative input.

It is useful to contrast this with QED. For massless charged matter, the perturbative beta function has the opposite sign,

β(e)=+b16π2e3+O(e5),b>0.\beta(e)=+{b\over16\pi^2}e^3+O(e^5), \qquad b>0.

In momentum variables, the charge grows toward the ultraviolet. In coordinate-space variables, with RR a distance scale and μ∼1/R\mu\sim1/R,

de2(R)dlog⁡R=−b8π2e4(R)+O(e6).{d e^2(R)\over d\log R}=-{b\over8\pi^2}e^4(R)+O(e^6).

Therefore

e2(R)=e2(a)1+b8π2e2(a)log⁡(R/a)e^2(R)={e^2(a)\over 1+{b\over8\pi^2}e^2(a)\log(R/a)}

at large RR in the massless theory. The potential between static charges has the schematic form

V(R)∼e2(R)R.V(R)\sim {e^2(R)\over R}.

If the beta function has no nonzero fixed point, then

e2(R)∼1log⁡(R/a)(R→∞),e^2(R)\sim {1\over \log(R/a)} \qquad (R\to\infty),

so the long-distance charge is screened. In the opposite direction, the one-loop formula develops a UV Landau pole. In the absence of a separate non-Gaussian ultraviolet fixed point, perturbative QED therefore does not furnish the asymptotically free continuum limit available in Yang–Mills theory.

This comparison is conceptually important. Both theories have logarithmic running. But in QED, the massless photon remains in the spectrum and the perturbative flow does not by itself produce a confinement scale. In pure Yang–Mills theory, the coupling grows in the infrared, and the natural expectation is a spectrum of gauge-invariant massive states with masses of order ΛYM\Lambda_{\mathrm{YM}}.

A relativistic quantum theory has a mass gap if the vacuum is separated from the rest of the physical spectrum by a positive energy. After choosing a vacuum sector and subtracting its energy, the infinite-volume definition is

Mgap=inf⁡(spec⁡H∖{0})>0.\boxed{ M_{\mathrm{gap}} =\inf\bigl(\operatorname{spec}H\setminus\{0\}\bigr)>0. }

This formulation does not assume that the first excitation is a discrete finite-volume level. For a Lorentz-invariant theory, MgapM_{\mathrm{gap}} is the smallest invariant mass in the physical Hilbert space. In Euclidean correlation functions, the same statement appears as exponential decay. If O\mathcal O is a local gauge-invariant operator with overlap with the lightest state in its channel, then at large Euclidean separation r=∣x∣r=|x|,

⟨O(x)O(0)⟩c∼e−MO∣x∣∣x∣α\langle \mathcal O(x)\mathcal O(0)\rangle_c \sim {e^{-M_{\mathcal O}|x|}\over |x|^\alpha}

for some power α\alpha. The smallest such MOM_{\mathcal O} over all nontrivial gauge-invariant channels is the mass gap.

Equivalently, when the asymptotic coefficient is nonzero,

MO=−lim⁡r→∞1rlog⁡∣⟨O(r)O(0)⟩c∣.M_{\mathcal O} =-\lim_{r\to\infty}{1\over r} \log\left|\langle\mathcal O(r)\mathcal O(0)\rangle_c\right|.

A spectrum with a vacuum, a mass gap, massive one-particle states, and a multiparticle continuum

A mass gap is a positive separation between the vacuum and the lightest physical excitation. In pure Yang–Mills theory the expected particles are gauge-invariant glueballs, with masses Mn=cnΛYMM_n=c_n\Lambda_{\mathrm{YM}}.

For pure Yang–Mills theory, the expected physical particles are glueballs: gauge-invariant excitations created, for example, by local operators such as tr⁡FμνFμν\operatorname{tr}F_{\mu\nu}F_{\mu\nu}. Their masses should have the form

Mn=cnΛYM,M_n=c_n\Lambda_{\mathrm{YM}},

where the cnc_n are dimensionless numbers. Perturbative RG determines the existence and scaling of ΛYM\Lambda_{\mathrm{YM}}, but it does not determine the constants cnc_n. Those constants are genuinely nonperturbative.

There is a useful distinction here:

RG perturbation theory: sets the scale Λ,\text{RG perturbation theory: }\quad \text{sets the scale }\Lambda,

but

strong dynamics: decides the spectrum in units of Λ.\text{strong dynamics: }\quad \text{decides the spectrum in units of }\Lambda.

This distinction keeps the logic honest. Asymptotic freedom strongly suggests where the mass gap should live, but the existence of a mass gap in four-dimensional pure Yang–Mills theory is a nonperturbative statement. The RG argument supplies the only possible scale; it does not by itself prove which gauge-invariant states exist or that the lightest one has strictly positive mass.

In practice one diagnoses a gap by looking for exponential decay in Euclidean correlators of gauge-invariant operators, not by looking for a pole in a gauge-dependent gluon propagator. The latter can be useful in a fixed gauge, but it is not itself the definition of a physical particle. A mass gap and confinement are also logically distinct properties: pure Yang–Mills theory is expected to have both, but the definition of either one does not imply the other in every quantum field theory.

Consider a dimensionless physical amplitude A\mathcal A depending on a characteristic Euclidean momentum pp and a coupling g(μ)g(\mu) defined at a subtraction scale μ\mu. RG invariance says schematically

(μ∂∂μ+β(g)∂∂g)A(pμ,g(μ))=0.\left(\mu{\partial\over\partial\mu}+\beta(g){\partial\over\partial g}\right)\mathcal A\left({p\over\mu},g(\mu)\right)=0.

For a single-scale observable, the solution may be written in terms of the running coupling at the physical momentum,

A=F(g(p)).\mathcal A=F(g(p)).

For p≫Λp\gg\Lambda, this becomes a perturbative expansion in

g2(p)=8π2b0log⁡(p/Λ).g^2(p)={8\pi^2\over b_0\log(p/\Lambda)}.

But the same RG equation also says that, after the arbitrary subtraction scale is eliminated, the only dimensionless argument left is

pΛ.{p\over\Lambda}.

Thus

A(p)=F(pΛ)\boxed{ \mathcal A(p)=\mathcal F\left({p\over\Lambda}\right) }

for a dimensionless observable in a massless one-coupling theory. More generally, if an observable O\mathcal O has mass dimension dOd_{\mathcal O}, then

O(pi)=ΛdO F(piΛ).\boxed{ \mathcal O(p_i)=\Lambda^{d_{\mathcal O}}\, \mathcal F\left({p_i\over\Lambda}\right). }

This is the practical content of dimensional transmutation. At finite cutoff, the same dimensionless amplitude can be parametrized as

A=A(g02,ap),\mathcal A=\mathcal A(g_0^2,ap),

but after taking the continuum limit along a trajectory of fixed Λ\Lambda,

A⟶F(pΛ).\mathcal A\longrightarrow \mathcal F\left({p\over\Lambda}\right).

Examples fix the dimensions. A glueball mass has dO=1d_{\mathcal O}=1, so M=cΛM=c\Lambda. A string tension has dO=2d_{\mathcal O}=2, so σ=cσΛ2\sigma=c_\sigma\Lambda^2. A dimensionless short-distance scattering amplitude depends on p/Λp/\Lambda and becomes expandable in 1/log⁡(p/Λ)1/\log(p/\Lambda) at large momentum.

For QCD with quark masses, additional dimensionless ratios remain:

A=F(pΛQCD,m^1ΛQCD,m^2ΛQCD,…).\mathcal A=\mathcal F\left({p\over\Lambda_{\mathrm{QCD}}},{\widehat m_1\over\Lambda_{\mathrm{QCD}}},{\widehat m_2\over\Lambda_{\mathrm{QCD}}},\ldots\right).

Here m^i\widehat m_i are RG-invariant quark-mass parameters with fixed normalization conventions. Equivalently one can use running masses in a specified scheme, retaining their scale dependence together with that of the coupling. Pure Yang–Mills has no such mass ratios. At fixed gauge group and topological angle, its dimensionless mass ratios are predictions of the theory.

The continuum limit of an asymptotically free gauge theory is structurally similar to the continuum limit of a statistical system at a second-order critical point. In both cases a microscopic length aa is sent to zero while a physical correlation length is held fixed.

In a critical spin system, use the dimensionless reduced temperature t=(T−Tc)/Tct=(T-T_c)/T_c with Tc>0T_c>0. On a specified side of the transition,

ξlat∼Aξ∣t∣−ν,\xi_{\mathrm{lat}}\sim A_\xi |t|^{-\nu},

where ξlat\xi_{\mathrm{lat}} is measured in lattice units and Aξ>0A_\xi>0 is a dimensionless, nonuniversal amplitude. The physical correlation length is

ξphys=aξlat,\xi_{\mathrm{phys}}=a\xi_{\mathrm{lat}},

so the physical mass scale is

M=1ξphys∼1aAξ∣t∣ν.M={1\over \xi_{\mathrm{phys}}} \sim {1\over aA_\xi}|t|^\nu.

To keep MM finite as a→0a\to0, one tunes T→TcT\to T_c.

In asymptotically free Yang–Mills theory, the weak-coupling UV fixed point plays the role of the critical point. The physical scale is

M∼ΛYM∼1aexp⁡[−8π2b0g02].M\sim \Lambda_{\mathrm{YM}} \sim {1\over a}\exp\left[-{8\pi^2\over b_0g_0^2}\right].

To keep MM finite as a→0a\to0, one tunes g02→0g_0^2\to0.

Compare the exponential bare-coupling tuning with the power-law thermal tuning in the two panels below. Their shared continuum-limit condition does not make the two mechanisms equivalent.

Yang–Mills exponential coupling tuning and critical-system power-law temperature tuning both keep a physical mass fixed as lattice spacing vanishes; the Yang–Mills mass requires separate nonperturbative input.

Send the microscopic spacing aa to zero while tuning a bare parameter so that a physical mass MM remains finite. The thermal variable is the reduced temperature t=(T−Tc)/Tct=(T-T_c)/T_c, and AξA_\xi is the dimensionless correlation-length amplitude. The Yang–Mills relation assumes a nonzero mass proportional to its RG scale; one-loop running does not prove that mass gap. Original schematic of exponential versus power-law tuning.

The analogy is not perfect. The Ising temperature perturbation is relevant at the critical point and produces power-law scaling. The Yang–Mills coupling is marginally relevant in the infrared, so the scaling is exponential in 1/g021/g_0^2. But the logic is the same: a continuum QFT is obtained by making the correlation length large compared with the cutoff.

Essential singularities and invisible perturbative masses

Section titled “Essential singularities and invisible perturbative masses”

The exponential dependence

Ma∼exp⁡[−8π2b0g02]M a\sim \exp\left[-{8\pi^2\over b_0g_0^2}\right]

has an essential singularity at g0=0g_0=0. For every nonnegative integer nn,

lim⁡g0→0exp⁡[−const/g02]g0n=0.\lim_{g_0\to0} {\exp[-\mathrm{const}/g_0^2]\over g_0^n}=0.

Therefore a mass gap of this type is invisible in ordinary perturbation theory. Any finite perturbative calculation around g0=0g_0=0 gives a power series in g0g_0 and logarithms; it cannot produce e−1/g02e^{-1/g_0^2}.

Removing the regulator and probing shorter physical distances are different operations. At fixed p/Λp/\Lambda, tuning the bare coupling gives

1g2(p)=b08π2log⁡Λ0Λ+b08π2log⁡pΛ0=b08π2log⁡pΛ.\begin{aligned} {1\over g^2(p)} &={b_0\over8\pi^2}\log{\Lambda_0\over\Lambda} +{b_0\over8\pi^2}\log{p\over\Lambda_0}\\ &={b_0\over8\pi^2}\log{p\over\Lambda}. \end{aligned}

The cutoff Λ0=a−1\Lambda_0=a^{-1} cancels: g0g_0 tends to zero as a→0a\to0, while the renormalized coupling g(p)g(p) stays fixed to this order. Regulator removal suppresses cutoff effects; increasing the physical ratio p/Λp/\Lambda improves the weak-coupling expansion. For a theory with MM of order Λ\Lambda, perturbation theory controls sufficiently large

p≫M,p\gg M,

while the mass gap concerns

p∼M.p\sim M.

The RG flow connects the two regimes, but it does not make the infrared perturbative.

A similar exponential scale appears in the Cooper channel. Consider a balanced three-dimensional two-component Fermi liquid with coherent, degenerate partners (k,↑)(\mathbf k,\uparrow) and (−k,↓)(-\mathbf k,\downarrow). Let ξ\xi be energy relative to the Fermi surface, and retain an instantaneous separable ss-wave attraction V<0V<0 in the sharp energy shell ∣ξ∣<ΛF|\xi|<\Lambda_F. Write ν=Npair(0)\nu=N_{\mathrm{pair}}(0) for the constant density of states of one species, per volume and per energy. Thus VV has units of energy times volume, while

gC=νVg_C=\nu V

is dimensionless and negative. The subscript distinguishes this channel coupling from the gauge coupling above.

The zero-temperature frequency integral in the normal-state pair bubble gives ∫dω/[2π(ω2+ξ2)]=1/(2∣ξ∣)\int d\omega/[2\pi(\omega^2+\xi^2)]=1/(2|\xi|). Including both sides of the Fermi surface therefore gives

Πpp(E)=ν∫E<∣ξ∣<ΛFdξ2∣ξ∣=νlog⁡ΛFE.\Pi_{pp}(E)=\nu\int_{E<|\xi|<\Lambda_F}{d\xi\over2|\xi|} =\nu\log{\Lambda_F\over E}.

Summing the leading ladders gives V(E)=V0/[1+V0Πpp(E)]V(E)=V_0/[1+V_0\Pi_{pp}(E)]. With ℓ=log⁡(ΛF/E)\ell=\log(\Lambda_F/E),

dVdℓ=−νV2,dgCdℓ=−gC2,gC(ℓ)=gC01+gC0ℓ,E∗=ΛFe−1/∣gC0∣.\begin{aligned} {dV\over d\ell}&=-\nu V^2, &{dg_C\over d\ell}&=-g_C^2,\\ g_C(\ell)&={g_{C0}\over1+g_{C0}\ell}, &E_*&=\Lambda_F e^{-1/|g_{C0}|}. \end{aligned}

Here E∗E_* is the formal normal-state pole scale; weak-coupling control is lost before the pole. It does not by itself calculate a gap or a transition temperature. The angular-channel flow has the same normalized form in Shankar 1993, v2, §VI.C, pp. 101–102, Eqs. (383)–(384), PDF; his channel variable has a different phase-space normalization.

To calculate a gap, add the homogeneous zero-temperature BCS mean-field approximation for this reduced-shell model, assuming ∣gC0∣≪1|g_{C0}|\ll1, ΛF≪EF\Lambda_F\ll E_F, and no competing instability or pair-breaking mismatch. The pairing field Δ>0\Delta>0 mixes particle and hole states through

hξ=(ξΔΔ−ξ),det⁡(E−hξ)=E2−ξ2−Δ2.h_\xi=\begin{pmatrix}\xi&\Delta\\\Delta&-\xi\end{pmatrix}, \qquad \det(E-h_\xi)=E^2-\xi^2-\Delta^2.

Diagonalization gives uξvξ=Δ/(2ξ2+Δ2)u_\xi v_\xi=\Delta/(2\sqrt{\xi^2+\Delta^2}). In the corresponding paired ground state, ⟨c−k↓ck↑⟩=−uξvξ\langle c_{-\mathbf k\downarrow}c_{\mathbf k\uparrow}\rangle=-u_\xi v_\xi. The self-consistency condition Δ=V0∫d3k/(2π)3 ⟨c−k↓ck↑⟩\Delta=V_0\int d^3k/(2\pi)^3\,\langle c_{-\mathbf k\downarrow}c_{\mathbf k\uparrow}\rangle therefore closes the BCS gap equation:

1=∣V0∣ν∫−ΛFΛFdξ2ξ2+Δ2=∣gC0∣arsinh⁡ΛFΔ,Δ=ΛFsinh⁡(1/∣gC0∣)=2ΛFe−1/∣gC0∣[1+O ⁣(e−2/∣gC0∣)].\begin{aligned} 1&=|V_0|\nu\int_{-\Lambda_F}^{\Lambda_F} {d\xi\over2\sqrt{\xi^2+\Delta^2}} =|g_{C0}|\operatorname{arsinh}{\Lambda_F\over\Delta},\\ \Delta&={\Lambda_F\over\sinh(1/|g_{C0}|)} =2\Lambda_F e^{-1/|g_{C0}|} \left[1+O\!\left(e^{-2/|g_{C0}|}\right)\right]. \end{aligned}

This is the constant-matrix-element construction of Bardeen, Cooper, and Schrieffer 1957, §II, pp. 1181–1183, Eqs. (2.34)–(2.40) and (2.50), with their positive attraction magnitude replaced by ∣V0∣|V_0| and their ϵ0\epsilon_0 by Δ\Delta. Their product of the one-spin density of states and interaction matrix element becomes ν∣V0∣\nu|V_0| in the per-volume convention here.

The exact saddle has minimum positive quasiparticle energy Δ\Delta when the dispersion crosses ξ=0\xi=0. The negative BdG branch is its particle–hole partner; the separation of the branches there is 2Δ2\Delta. In this model Δ/E∗→2\Delta/E_*\to2 at weak coupling. That prefactor comes from the specified gap integral, not from the normal-state leading logarithm. Follow the conditional step between the two calculations in the figure.

The normal-state Cooper coupling defines a breakdown scale E star; only after imposing the stated BCS saddle assumptions does the gap equation yield Delta and particle–hole branches with minimum positive energy Delta.

The upper normal-state calculation uses gC=νVg_C=\nu V and defines E∗E_*. The dashed arrow requires the stated constant-DOS, sharp-shell, instantaneous BCS mean-field model at zero temperature; its gap obeys Δ/E∗→2\Delta/E_*\to2 at weak coupling. Below, solid curves show the exact normalized BdG branches ±1+(ξ/Δ)2\pm\sqrt{1+(\xi/\Delta)^2} and dashed lines the normal crossing. The positive minimum is Δ\Delta; the negative branch is a particle–hole partner. Original schematic of the conditional calculation with exact dimensionless dispersion curves, not material data.

The analogy has limits. A Fermi surface is not Lorentz invariant. A nonzero quasiparticle gap is not a proof of phase coherence or a gap to every many-body excitation: a neutral superfluid can retain a gapless phase mode. Mismatch, temperature, or competing channels can also stop the normal-state flow. The shared conclusion is that a logarithm can generate an exponentially small scale; identifying the resulting state requires additional dynamics.

What is and is not proved by the one-loop flow

Section titled “What is and is not proved by the one-loop flow”

The leading beta function establishes weak ultraviolet running near the Gaussian fixed point when b0>0b_0>0. Its formal infrared pole indicates the limit of the one-loop extrapolation; it neither excludes an infrared fixed point nor proves confinement, a discrete glueball spectrum, or a positive mass gap.

For pure Yang–Mills theory, the physical expectation is

Mgap∼ΛYM,M_{\mathrm{gap}}\sim \Lambda_{\mathrm{YM}},

and lattice gauge theory strongly supports this picture. But from the viewpoint of continuum perturbation theory, the statement remains nonperturbative. This is exactly why dimensional transmutation is so important: it shows how the scale can exist without pretending that weak-coupling diagrams calculate the infrared spectrum.

A good slogan is:

Asymptotic freedom explains why the scale is generated;\text{Asymptotic freedom explains why the scale is generated;} strong dynamics decides what happens at that scale.\text{strong dynamics decides what happens at that scale.}

An asymptotically free theory with

β(g)=−b016π2g3,b0>0,\beta(g)=-{b_0\over16\pi^2}g^3, \qquad b_0>0,

has the one-loop running

1g2(μ)=b08π2log⁡μΛ.{1\over g^2(\mu)}={b_0\over8\pi^2}\log{\mu\over\Lambda}.

The integration constant may be written as the RG-invariant scale

Λ=μexp⁡[−8π2b0g2(μ)].\Lambda=\mu\exp\left[-{8\pi^2\over b_0g^2(\mu)}\right].

This is dimensional transmutation: the dimensionless coupling is traded for a dimensionful scale.

With a cutoff aa, the continuum limit of pure Yang–Mills theory is obtained by taking

a→0,g02→0,1aexp⁡[−8π2b0g02]=fixed.a\to0, \qquad g_0^2\to0, \qquad {1\over a}\exp\left[-{8\pi^2\over b_0g_0^2}\right] =\text{fixed}.

This resembles the critical continuum limit of a statistical system, where T→TcT\to T_c while aξlata\xi_{\mathrm{lat}} is held fixed.

A mass gap means a positive separation between the vacuum and the lightest physical state. In pure Yang–Mills theory, the expected spectrum consists of gauge-invariant massive states with

Mn=cnΛYM.M_n=c_n\Lambda_{\mathrm{YM}}.

The scale ΛYM\Lambda_{\mathrm{YM}} is visible from perturbative RG, but the constants cnc_n and the existence of the gap are nonperturbative.

Claiming that perturbation theory computes the gap. Perturbation theory determines the high-momentum running and identifies the RG-invariant scale. The gap itself is infrared and nonperturbative.

Treating a one-loop pole as an exact singularity. The pole is a warning that the chosen perturbative variables have failed, not a controlled statement about the exact coupling.

Confusing scheme dependence with physical arbitrariness. A coupling redefinition changes the numerical value assigned to ΛS\Lambda_{\mathcal S}, but physical masses and mass ratios do not depend on that convention. One fixes a scheme for Λ\Lambda and quotes nonperturbative quantities in that scheme or as dimensionless ratios.

Holding the bare coupling fixed while removing the cutoff. Fixed g0g_0 leaves the correlation length finite in cutoff units. The continuum limit requires g0→0g_0\to0 so that aΛ→0a\Lambda\to0 while Λ\Lambda in physical units remains fixed.

Equating a mass gap with confinement. They are distinct infrared properties. Pure Yang–Mills theory is expected to possess both, but one-loop asymptotic freedom proves neither.

Assuming every marginal interaction generates a gap. Some marginal interactions are exactly marginal or marginally irrelevant. The sign of the beta function and the interaction channel matter.

Exercise 1: Verify one-loop RG invariance of the transmutation scale

Section titled “Exercise 1: Verify one-loop RG invariance of the transmutation scale”

Let

β(g)=−b016π2g3,b0>0.\beta(g)=-{b_0\over16\pi^2}g^3, \qquad b_0>0.

Show that

Λ=μexp⁡[−8π2b0g2(μ)]\Lambda=\mu\exp\left[-{8\pi^2\over b_0g^2(\mu)}\right]

is independent of μ\mu at one loop.

Solution

First compute

ddlog⁡μ1g2=−2g3dgdlog⁡μ=−2g3β(g).{d\over d\log\mu}{1\over g^2} =-{2\over g^3}{dg\over d\log\mu} =-{2\over g^3}\beta(g).

Using the beta function,

ddlog⁡μ1g2=−2g3[−b016π2g3]=b08π2.{d\over d\log\mu}{1\over g^2} =-{2\over g^3}\left[-{b_0\over16\pi^2}g^3\right] ={b_0\over8\pi^2}.

Now differentiate

log⁡Λ=log⁡μ−8π2b0g2(μ).\log\Lambda=\log\mu-{8\pi^2\over b_0g^2(\mu)}.

This gives

dlog⁡Λdlog⁡μ=1−8π2b0ddlog⁡μ1g2=1−8π2b0b08π2=0.{d\log\Lambda\over d\log\mu} =1-{8\pi^2\over b_0}{d\over d\log\mu}{1\over g^2} =1-{8\pi^2\over b_0}{b_0\over8\pi^2}=0.

Thus Λ\Lambda is RG invariant at one loop.

Exercise 2: Tune the bare coupling toward the continuum limit

Section titled “Exercise 2: Tune the bare coupling toward the continuum limit”

With cutoff Λ0=a−1\Lambda_0=a^{-1}, suppose

Λ=Λ0exp⁡[−8π2b0g02].\Lambda=\Lambda_0\exp\left[-{8\pi^2\over b_0g_0^2}\right].

Solve for g02g_0^2 in terms of aa and fixed Λ\Lambda. What happens as a→0a\to0? Compare g(p)g(p) at the two cutoffs Λ0=100Λ\Lambda_0=100\Lambda and 1000Λ1000\Lambda, keeping p=10Λp=10\Lambda fixed, within the one-loop model.

Solution

Since Λ0=a−1\Lambda_0=a^{-1},

Λ=1aexp⁡[−8π2b0g02].\Lambda={1\over a}\exp\left[-{8\pi^2\over b_0g_0^2}\right].

Multiply by aa and take the logarithm:

log⁡(aΛ)=−8π2b0g02.\log(a\Lambda)=-{8\pi^2\over b_0g_0^2}.

Equivalently,

log⁡1aΛ=8π2b0g02.\log{1\over a\Lambda}={8\pi^2\over b_0g_0^2}.

Therefore

g02(a)=8π2b0log⁡(1/aΛ).\boxed{ g_0^2(a)={8\pi^2\over b_0\log(1/a\Lambda)}. }

As a→0a\to0 at fixed Λ\Lambda, the logarithm diverges and hence

g02(a)→0.g_0^2(a)\to0.

The continuum limit approaches the weak-coupling UV fixed point. In the two-cutoff comparison, set C=b0/(8π2)C=b_0/(8\pi^2). Then g0−2g_0^{-2} changes from Clog⁡100C\log100 to Clog⁡1000C\log1000, but

g−2(10Λ)=Clog⁡100+Clog⁡(10/100)=Clog⁡1000+Clog⁡(10/1000)=Clog⁡10.\begin{aligned} g^{-2}(10\Lambda) &=C\log100+C\log(10/100)\\ &=C\log1000+C\log(10/1000) =C\log10. \end{aligned}

Thus the renormalized coupling at this fixed physical momentum does not decrease when the cutoff is increased. The comparison tests cutoff cancellation within the one-loop model, not the numerical accuracy of that model at p/Λ=10p/\Lambda=10.

Exercise 3: Derive the scaling form after dimensional transmutation

Section titled “Exercise 3: Derive the scaling form after dimensional transmutation”

Assume a dimensionless physical observable satisfies the RG equation

(μ∂∂μ+β(g)∂∂g)A(pμ,g(μ))=0.\left(\mu{\partial\over\partial\mu}+\beta(g){\partial\over\partial g}\right) \mathcal A\left({p\over\mu},g(\mu)\right)=0.

Explain why, in a massless one-coupling asymptotically free theory, it can be written as

A=F(pΛ).\mathcal A=\mathcal F\left({p\over\Lambda}\right).
Solution

The observable is dimensionless, so before solving the RG equation it can only depend on the dimensionless ratio p/μp/\mu and on the running coupling g(μ)g(\mu).

The RG equation says that changing μ\mu while moving g(μ)g(\mu) along the RG trajectory leaves A\mathcal A unchanged. Therefore the only invariant data on the RG trajectory can appear. For a one-coupling asymptotically free theory, the invariant data are encoded in

Λ=μexp⁡[−8π2b0g2(μ)]\Lambda=\mu\exp\left[-{8\pi^2\over b_0g^2(\mu)}\right]

at one loop, with higher-loop refinements changing only the precise scheme definition of Λ\Lambda.

Since pp is the only external scale and Λ\Lambda is the only intrinsic scale, dimensional analysis then gives

A=F(pΛ).\boxed{ \mathcal A=\mathcal F\left({p\over\Lambda}\right). }

At p≫Λp\gg\Lambda, the same statement can be rewritten as a perturbative expansion in g(p)g(p).

Exercise 4: Solve the coordinate-space QED flow

Section titled “Exercise 4: Solve the coordinate-space QED flow”

In coordinate-space QED, suppose

de2(R)dlog⁡R=−Be4(R),B>0.{d e^2(R)\over d\log R}=-B e^4(R), \qquad B>0.

Solve for e2(R)e^2(R) with boundary condition e(a)=e0e(a)=e_0. What is the large-RR behavior?

Solution

Write

ddlog⁡R1e2(R)=−1e4(R)de2(R)dlog⁡R=B.{d\over d\log R}{1\over e^2(R)} =-{1\over e^4(R)}{d e^2(R)\over d\log R} =B.

Integrating from aa to RR gives

1e2(R)=1e02+Blog⁡Ra.{1\over e^2(R)}={1\over e_0^2}+B\log{R\over a}.

Thus

e2(R)=e021+Be02log⁡(R/a).\boxed{ e^2(R)={e_0^2\over1+Be_0^2\log(R/a)}. }

For R≫aR\gg a,

e2(R)∼1Blog⁡(R/a).e^2(R)\sim {1\over B\log(R/a)}.

The charge is screened at long distance. This running alone does not generate a mass gap for the photon.

Exercise 5: Tune a statistical system to its critical limit

Section titled “Exercise 5: Tune a statistical system to its critical limit”

Near an Ising critical point, use t=(T−Tc)/Tct=(T-T_c)/T_c and a fixed dimensionless amplitude Aξ>0A_\xi>0:

ξlat∼Aξ∣t∣−ν.\xi_{\mathrm{lat}}\sim A_\xi |t|^{-\nu}.

The physical correlation length is ξphys=aξlat\xi_{\mathrm{phys}}=a\xi_{\mathrm{lat}} and the physical mass is M=1/ξphysM=1/\xi_{\mathrm{phys}}. Determine how tt must scale with aa if MM is held fixed while a→0a\to0. Translate back to the dimensionful temperature difference.

Solution

The physical mass is

M=1aξlat∼1aAξ∣t∣ν.M={1\over a\xi_{\mathrm{lat}}} \sim {1\over aA_\xi}|t|^\nu.

Holding MM fixed gives

∣t∣ν∼aAξM.|t|^\nu\sim aA_\xi M.

Therefore

∣t∣∼(aAξM)1/ν,∣T−Tc∣∼Tc(aAξM)1/ν.\boxed{ |t|\sim (aA_\xi M)^{1/\nu}, \qquad |T-T_c|\sim T_c(aA_\xi M)^{1/\nu}. }

As a→0a\to0, the temperature must be tuned to TcT_c. This is the statistical-mechanics analog of tuning g0→0g_0\to0 in an asymptotically free gauge theory while keeping the physical mass scale fixed.

Exercise 6: Distinguish the Cooper pole from the BCS gap

Section titled “Exercise 6: Distinguish the Cooper pole from the BCS gap”

For the dimensionful channel interaction VV, take a constant one-species density of states ν\nu and define gC=νVg_C=\nu V. The leading flow is

dgCdℓ=−gC2,gC0<0.{dg_C\over d\ell}=-g_C^2, \qquad g_{C0}<0.

Find the formal pole scale E∗E_*. Under the reduced-shell BCS assumptions in the text, solve the gap integral and compare Δ/E∗\Delta/E_* for gC0=−0.20g_{C0}=-0.20. Would replacing ν\nu by the spin-summed density of states while leaving V0V_0 unchanged preserve the prediction?

Solution

Integrate

dgCgC2=−dℓ.{dg_C\over g_C^2}=-d\ell.

This gives

−1gC(ℓ)=−1gC0−ℓ,-{1\over g_C(\ell)}=-{1\over g_{C0}}-\ell,

or

1gC(ℓ)=1gC0+ℓ.{1\over g_C(\ell)}={1\over g_{C0}}+\ell.

Since gC0<0g_{C0}<0, the denominator reaches zero at

ℓ∗=1∣gC0∣.\ell_*={1\over|g_{C0}|}.

If the running energy scale is E=ΛFe−ℓE=\Lambda_F e^{-\ell}, the formal normal-state pole is

E∗=ΛFe−ℓ∗=ΛFexp⁡[−1∣gC0∣].\boxed{ E_*=\Lambda_F e^{-\ell_*} =\Lambda_F\exp\left[-{1\over|g_{C0}|}\right]. }

The BCS integral is a separate calculation. With ξ=Δsinh⁡u\xi=\Delta\sinh u, it gives

1∣gC0∣=∫0ΛFdξξ2+Δ2=arsinh⁡ΛFΔ.{1\over|g_{C0}|} =\int_0^{\Lambda_F}{d\xi\over\sqrt{\xi^2+\Delta^2}} =\operatorname{arsinh}{\Lambda_F\over\Delta}.

Consequently,

ΔE∗=21−e−2/∣gC0∣,E∗ΛF=e−5≃0.00673795,ΔΛF=1sinh⁡5≃0.01347651,ΔE∗≃2.00009080.\begin{aligned} {\Delta\over E_*}&={2\over1-e^{-2/|g_{C0}|}},\\ {E_*\over\Lambda_F}&=e^{-5}\simeq0.00673795,\\ {\Delta\over\Lambda_F}&={1\over\sinh5}\simeq0.01347651, \qquad {\Delta\over E_*}\simeq2.00009080. \end{aligned}

Using a spin-summed 2ν2\nu with unchanged V0V_0 incorrectly doubles ∣gC0∣|g_{C0}|. It predicts E∗/ΛF=e−2.5≃0.0820850E_*/\Lambda_F=e^{-2.5}\simeq0.0820850, an exponentially different scale. A notation change must preserve the product of the correctly counted pair density of states and the matched interaction. Neither normal-state pole alone establishes a paired ground state or phase coherence.

  • Bardeen, John, Leon N. Cooper, and J. Robert Schrieffer. “Theory of Superconductivity.” Physical Review 108, no. 5 (1957): 1175–1204. DOI. Open PDF.
  • Shankar, Ramamurti. “Renormalization Group Approach to Interacting Fermions.” Reviews of Modern Physics 66 (1994): 129–192. DOI. Cited edition: June 1993 preprint, arXiv:cond-mat/9307009v2. Open PDF.
  • Coleman, Sidney, and Erick Weinberg. “Radiative Corrections as the Origin of Spontaneous Symmetry Breaking.” Physical Review D 7, no. 6 (1973): 1888–1910.
  • Gell-Mann, Murray, and Francis E. Low. “Quantum Electrodynamics at Small Distances.” Physical Review 95, no. 5 (1954): 1300–1312.
  • Gross, David J., and Frank Wilczek. “Ultraviolet Behavior of Non-Abelian Gauge Theories.” Physical Review Letters 30, no. 26 (1973): 1343–1346.
  • Politzer, H. David. “Reliable Perturbative Results for Strong Interactions?” Physical Review Letters 30, no. 26 (1973): 1346–1349.
  • Polyakov, Alexander M. Gauge Fields and Strings. Chur: Harwood Academic Publishers, 1987, Chapter 2.
  • Schwartz, Matthew D. Quantum Field Theory and the Standard Model. Cambridge: Cambridge University Press, 2014, Chapters 23 and 26.
  • Srednicki, Mark. Quantum Field Theory. Cambridge: Cambridge University Press, 2007, Sections 28, 73, and 82.
  • Weinberg, Steven. The Quantum Theory of Fields. Volume II: Modern Applications. Cambridge: Cambridge University Press, 1996, Chapter 18.
  • Wilson, Kenneth G., and John Kogut. “The Renormalization Group and the ϵ\epsilon Expansion.” Physics Reports 12, no. 2 (1974): 75–199.
  • Zinn-Justin, Jean. Quantum Field Theory and Critical Phenomena. 4th ed. Oxford: Clarendon Press, 2002.

Original QFT.org content:CC BY 4.0, unless an item supplies different terms. Third-party material retains its own terms.