Skip to content

QED and Yang–Mills Beta Functions

The previous pages developed three views of vacuum polarization: a momentum-space transverse tensor, a running charge, and a proper-time determinant in a background magnetic field. We now combine them into the one-loop beta functions of Abelian and non-Abelian gauge theory.

The central physical question is deceptively simple: why does QED screen charge while Yang–Mills theory antiscreens it? Charged scalars and fermions polarize the vacuum in a way that weakens long-distance electric fields. Gauge bosons in a non-Abelian theory also carry charge, but their spin-one magnetic moment produces a stronger paramagnetic response with the opposite sign. The final result is asymptotic freedom: at short distances, the Yang–Mills coupling becomes small.

The background-field method computes this coefficient while preserving gauge covariance of the background. One separates a smooth background from the quantum fluctuation and extracts the local ultraviolet correction to its F2F^2 term. This page derives that coefficient in background Feynman gauge, explains the spin decomposition, and states the general one-loop result with matter. The calculation concerns perturbative ultraviolet running; it does not establish a stable chromomagnetic vacuum or an infrared mass gap.

Required background. Effective Actions in Background Fields supplies the proper-time determinants and the orbital-versus-spin decomposition used below. Running Charge, Screening, and Antiscreening fixes the inverse-coupling interpretation of screening. All one-loop results below are translated into b0b_0, defined by β(g)=−b0g3/(16π2)+O(g5)\beta(g)=-b_0g^3/(16\pi^2)+O(g^5). One Dirac or Weyl fermion in representation RR contributes respectively −4T(R)/3-4T(R)/3 or −2T(R)/3-2T(R)/3; one complex or real scalar contributes −T(R)/3-T(R)/3 or −T(R)/6-T(R)/6; and pure Yang–Mills contributes +11CA/3+11C_A/3. Thus Abelian matter gives b0<0b_0<0 and β(e)>0\beta(e)>0.

Gauge-field and beta-function normalization. The background-field calculation is Euclidean. We use the rescaled gauge-field normalization

Γ[A]⊃14∫d4x 1g2FμνaFμνa\boxed{ \Gamma[A]\supset {1\over4}\int d^4x\,{1\over g^2}F^a_{\mu\nu}F^a_{\mu\nu} }

or, equivalently, Γ[A]⊃(2g2)−1∫tr⁡FμνFμν\Gamma[A]\supset (2g^2)^{-1}\int \operatorname{tr}F_{\mu\nu}F_{\mu\nu} when tr⁡(TaTb)=δab/2\operatorname{tr}(T^aT^b)=\delta^{ab}/2. The coupling sits in front of the gauge kinetic term; covariant derivatives in the fluctuation operator use the rescaled background field.

For Abelian matter of unit charge,

Dμ=∂μ−iAμ,[Dμ,Dν]=−iFμν.D_\mu=\partial_\mu-iA_\mu, \qquad [D_\mu,D_\nu]=-iF_{\mu\nu}.

For adjoint non-Abelian fluctuations,

DˉμX=∂μX+[Aˉμ,X],[Dˉμ,Dˉν]X=[Fˉμν,X].\bar D_\mu X=\partial_\mu X+[\bar A_\mu,X], \qquad [\bar D_\mu,\bar D_\nu]X=[\bar F_{\mu\nu},X].

Here, only in the adjoint fluctuation formulas, Aˉμ\bar A_\mu denotes the anti-Hermitian connection Aˉμ=−iAˉμ,HaTa\bar A_\mu=-i\bar A^a_{\mu,H}T^a, and Fˉμν=−iFˉμνaTa\bar F_{\mu\nu}=-i\bar F^a_{\mu\nu}T^a. Thus DˉμX=∂μX−i[Aˉμ,H,X]\bar D_\mu X=\partial_\mu X-i[\bar A_{\mu,H},X] is precisely the Hermitian-generator convention above. The components Fˉμνa\bar F^a_{\mu\nu} are real; group invariants below are defined with the Hermitian matrices TaT^a. Relative to an unrescaled Hermitian connection, AˉH=gAphys\bar A_H=gA_{\mathrm{phys}} and Fˉa=gFphysa\bar F^a=gF^a_{\mathrm{phys}}, so the fluctuation curvature carries no additional factor of gg in our formulas.

The beta function is

β(g)=μdgdμ.\beta(g)=\mu {dg\over d\mu}.

We write the Yang–Mills result as

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

Thus b0>0b_0>0 means asymptotic freedom.

In the rescaled normalization, the two-point function of the background gauge field defines the running inverse coupling. For an Abelian field in Euclidean momentum space,

Γ(2)[A]=12∫qAμ(−q)(q2δμν−qμqν)1e2(q)Aν(q)+⋯ .\Gamma^{(2)}[A] ={1\over2}\int_q A_\mu(-q) \left(q^2\delta_{\mu\nu}-q_\mu q_\nu\right) {1\over e^2(q)}A_\nu(q)+\cdots.

The transverse tensor is forced by gauge invariance. The entire one-loop question is the coefficient of the logarithm in 1/e2(q)1/e^2(q).

For NfN_f Dirac fermions and NsN_s charged complex scalars of unit Abelian charge,

1e2(q)=1e2(μ)+116π2(43Nf+13Ns)log⁡μ2q2.{1\over e^2(q)} ={1\over e^2(\mu)} +{1\over16\pi^2} \left({4\over3}N_f+{1\over3}N_s\right) \log{\mu^2\over q^2}.

Only the one-loop logarithm is displayed; finite matching terms depend on the definition of the effective charge, and higher-loop corrections begin at the next order in e2e^2. Differentiating at fixed bare coupling gives

βQED(e)=e316π2(43Nf+13Ns)+O(e5).\boxed{ \beta_{\mathrm{QED}}(e) ={e^3\over16\pi^2} \left({4\over3}N_f+{1\over3}N_s\right)+O(e^5). }

The sign is positive. As the renormalization scale μ\mu is raised, the charge increases. Equivalently, at longer distances the charge is screened.

For Yang–Mills theory it is more natural to write

1g2(q)=1g2(μ)+b016π2log⁡q2μ2,{1\over g^2(q)} ={1\over g^2(\mu)} +{b_0\over16\pi^2}\log{q^2\over\mu^2},

so that

β(g)=−b016π2g3+O(g5).\boxed{ \beta(g)=-{b_0\over16\pi^2}g^3+O(g^5). }

The sign has reversed relative to QED if b0>0b_0>0: the inverse coupling grows in the ultraviolet, so the coupling itself shrinks.

Opposite one-loop flows of Abelian and non-Abelian inverse couplings

QED screening and Yang–Mills antiscreening in the inverse-coupling language. At one loop, 1/e21/e^2 decreases toward the ultraviolet, while 1/g21/g^2 in pure Yang–Mills increases toward the ultraviolet.

The coefficient being computed is that of the transverse background-field kinetic term. Background gauge covariance constrains its form; the full off-shell effective action can still depend on the gauge parameter used for the quantum fluctuation. It is not itself a physical observable. The distinction and the relation to coupling renormalization are developed in Background-Field Yang–Mills Effective Action.

Here is the sign dictionary used in the rest of the course:

theoryUV behavior of 1/g2(μ)UV behavior of g(μ)QED matterdecreasesincreasespure Yang–Millsincreasesdecreases\begin{array}{c|c|c} \text{theory} & \text{UV behavior of }1/g^2(\mu) & \text{UV behavior of }g(\mu) \\ \hline \text{QED matter} & \text{decreases} & \text{increases} \\ \text{pure Yang–Mills} & \text{increases} & \text{decreases} \end{array}

Most sign mistakes in this calculation come from silently switching between gg, g2g^2, and 1/g21/g^2.

This statement is also why the background-field method is more than a shortcut. In ordinary gauges one separately computes wavefunction, vertex, and coupling counterterms. In background-field gauge, background gauge invariance ties the renormalization of the background two-point function directly to the coupling renormalization.

The matter coefficients can be read directly from the proper-time result of the previous page. A complex scalar in a constant magnetic field gives

Δ(1e2)s=116π213log⁡Λ2m2.\Delta\left({1\over e^2}\right)_s ={1\over16\pi^2}{1\over3}\log {\Lambda^2\over m^2}.

A Dirac fermion gives

Δ(1e2)f=116π243log⁡Λ2m2.\Delta\left({1\over e^2}\right)_f ={1\over16\pi^2}{4\over3}\log {\Lambda^2\over m^2}.

The scalar contribution is orbital. The fermion contribution contains both orbital motion and the Pauli spin coupling. In the magnetic-field heat kernel, this distinction appears in the expansion

4xsinh⁡xcosh⁡x=4(1−x26+⋯ )(1+x22+⋯ )=4+43x2+⋯ .4{x\over\sinh x}\cosh x =4\left(1-{x^2\over6}+\cdots\right) \left(1+{x^2\over2}+\cdots\right) =4+{4\over3}x^2+\cdots.

The −2/3-2/3 orbital piece and the +2+2 spin piece add to 4/34/3:

−23+2=43.-{2\over3}+2={4\over3}.

Both scalar and fermion matter screen. In a non-Abelian theory, matter fields still screen. The only new ingredient is the group-theory weight.

Let RR be a representation of the gauge group, with Hermitian generators normalized by

tr⁡R(TaTb)=TRδab.\operatorname{tr}_R(T^aT^b)=T_R\delta^{ab}.

Then one Dirac fermion in RR contributes

Δb0Dirac=−43TR,\Delta b_0^{\mathrm{Dirac}}=-{4\over3}T_R,

and one complex scalar in RR contributes

Δb0complex scalar=−13TR.\Delta b_0^{\mathrm{complex\ scalar}}=-{1\over3}T_R.

The minus signs appear because b0b_0 is defined by β(g)=−b0g3/(16π2)\beta(g)=-b_0g^3/(16\pi^2). Matter has the QED sign, so it reduces b0b_0 and pushes the theory away from asymptotic freedom. This is often the cleanest way to remember the formula: gauge bosons add to b0b_0, ordinary matter subtracts from b0b_0.

For a real scalar the coefficient is half as large,

Δb0real scalar=−16TR,\Delta b_0^{\mathrm{real\ scalar}}=-{1\over6}T_R,

and for a Weyl fermion it is half a Dirac fermion,

Δb0Weyl=−23TR.\Delta b_0^{\mathrm{Weyl}}=-{2\over3}T_R.

These factors are often the easiest way to check supersymmetric cancellations.

Background-field expansion of Yang–Mills theory

Section titled “Background-field expansion of Yang–Mills theory”

Now consider pure Yang–Mills theory. Write the gauge field as a background plus a quantum fluctuation,

Aμ=Aˉμ+aμ.A_\mu=\bar A_\mu+a_\mu.

The field strength expands as

Fμν[Aˉ+a]=Fˉμν+Dˉμaν−Dˉνaμ+[aμ,aν].F_{\mu\nu}[\bar A+a] =\bar F_{\mu\nu} +\bar D_\mu a_\nu-\bar D_\nu a_\mu +[a_\mu,a_\nu].

A background gauge transformation acts by

Aˉμ↦h−1Aˉμh+h−1∂μh,aμ↦h−1aμh.\bar A_\mu\mapsto h^{-1}\bar A_\mu h+h^{-1}\partial_\mu h, \qquad a_\mu\mapsto h^{-1}a_\mu h.

Thus the background is a connection, while the fluctuation transforms homogeneously. The gauge condition

Dˉμaμ=0\bar D_\mu a_\mu=0

preserves background gauge invariance. This is the reason the method is so efficient: the effective action generated for Aˉμ\bar A_\mu must be a gauge-invariant functional of the background.

Expanding the Yang–Mills action to quadratic order in aμa_\mu and adding background Feynman gauge fixing gives the vector fluctuation operator

Lμνvec=−Dˉ2δμν−2ad⁡(Fˉμν),\mathcal L_{\mu\nu}^{\mathrm{vec}} =-\bar D^2\delta_{\mu\nu}-2\operatorname{ad}(\bar F_{\mu\nu}),

where

ad⁡(Fˉμν)X=[Fˉμν,X].\operatorname{ad}(\bar F_{\mu\nu})X=[\bar F_{\mu\nu},X].

The Faddeev–Popov ghosts are Grassmann scalar fields in the adjoint representation with operator

Lgh=−Dˉ2.\mathcal L_{\mathrm{gh}}=-\bar D^2.

Therefore the one-loop effective action is schematically

ΓYM(1)[Aˉ]=12Tr⁡1,adjlog⁡(−Dˉ2δμν−2ad⁡Fˉμν)−Tr⁡0,adjlog⁡(−Dˉ2).\boxed{ \Gamma^{(1)}_{\mathrm{YM}}[\bar A] ={1\over2}\operatorname{Tr}_{1,\mathrm{adj}}\log \left(-\bar D^2\delta_{\mu\nu}-2\operatorname{ad}\bar F_{\mu\nu}\right) -\operatorname{Tr}_{0,\mathrm{adj}}\log(-\bar D^2). }

The factor 1/21/2 is the Gaussian determinant of a real bosonic vector field. The ghost and antighost form one independent Grassmann pair per adjoint component, giving a determinant with weight −1-1, not −1/2-1/2. The vector trace includes Lorentz-vector indices and adjoint color indices; the ghost trace includes only adjoint color indices. These are the unrestricted vector and ghost determinants in background Feynman gauge, as in Vassilevich 2003, §3.4, pp.27–28, Eqs.(3.46), (3.49), (3.51)–(3.53), preprint PDF. A strictly transverse vector determinant has a different ghost power and must not be substituted into this formula.

Background-field expansion with vector and ghost determinants

The background-field split A=Aˉ+aA=\bar A+a preserves manifest gauge covariance in the background. Gauge fluctuations contribute a spin-one determinant, while ghosts contribute scalar adjoint determinants that remove unphysical gauge modes.

The CP-even logarithmic divergence quadratic in a slowly varying background must be proportional to

∫d4x FˉμνaFˉμνa,\int d^4x\,\bar F^a_{\mu\nu}\bar F^a_{\mu\nu},

up to total derivatives. The other familiar dimension-four gauge invariant, tr⁡FˉμνFˉ~μν\operatorname{tr}\bar F_{\mu\nu}\widetilde{\bar F}_{\mu\nu}, is CP odd and topological, so it is not generated by this CP-even fluctuation determinant. Extracting the coefficient of Fˉ2\bar F^2 gives the beta function.

Computing the ultraviolet heat coefficient

Section titled “Computing the ultraviolet heat coefficient”

Take a smooth background in flat Euclidean space and a compact simple gauge group. Use a periodic box, or sufficient falloff to discard integrated total derivatives. We need only the short-proper-time expansion. A finite upper cutoff sIRs_{\mathrm{IR}} separates it from possible zero or negative modes of the vector operator; no convergence of the full integral at s→∞s\to\infty is assumed. The distinction between ultraviolet coefficients and infrared convergence is explicit in Vassilevich 2003, §1, pp.7–8, Eqs.(1.16)–(1.21), preprint PDF.

Write an operator of Laplace type as

P=−(Dˉ2+E),Ωμν=[Dˉμ,Dˉν].P=-(\bar D^2+E), \qquad \Omega_{\mu\nu}=[\bar D_\mu,\bar D_\nu].

For its diagonal heat trace in four flat dimensions,

Tr⁡e−sP∼1(4πs)2∫d4x tr⁡[1+sE+s2(E22+ΩμνΩμν12+Dˉ2E6)+⋯ ].\begin{aligned} \operatorname{Tr}e^{-sP} \sim {1\over(4\pi s)^2}\int d^4x\, \operatorname{tr}\Bigl[\mathbf 1+sE +s^2\Bigl({E^2\over2} +{\Omega_{\mu\nu}\Omega_{\mu\nu}\over12} +{\bar D^2E\over6}\Bigr)+\cdots\Bigr]. \end{aligned}

Here Tr⁡\operatorname{Tr} includes the spacetime integration, whereas tr⁡\operatorname{tr} acts on the finite internal indices. This local expansion, with exactly the sign convention for PP above, is Vassilevich 2003, §2.1, p.11, Eq.(2.2), and §4.1, pp.39–40, Eqs.(4.25), (4.28), preprint PDF. It is the mathematical input we use; a general proof for curved spaces or boundaries is outside this lesson. Heat Kernels, Zeta Functions, and Spectral Determinants explains the broader operator assumptions.

The two coefficients relevant here also have simple checks. A constant commuting EE multiplies the free heat kernel by esEe^{sE}, fixing the E2/2E^2/2 term. For a unit-charge scalar with F12=BF_{12}=B, the previous lesson gave

Bssinh⁡Bs=1−B2s26+O(s4).{Bs\over\sinh Bs}=1-{B^2s^2\over6}+O(s^4).

Since Ω12=−iB\Omega_{12}=-iB and both ordered index pairs occur, ΩμνΩμν=−2B2\Omega_{\mu\nu}\Omega_{\mu\nu}=-2B^2. Thus the curvature coefficient is 1/121/12. Gauge covariance and dimensional counting leave precisely these two quadratic invariants, plus a total derivative, in this flat CP-even coefficient.

Define

Q=CAFˉμνaFˉμνa,tr⁡adj(TadjaTadjb)=CAδab.\mathcal Q=C_A\bar F^a_{\mu\nu}\bar F^a_{\mu\nu}, \qquad \operatorname{tr}_{\mathrm{adj}}(T^a_{\mathrm{adj}}T^b_{\mathrm{adj}}) =C_A\delta^{ab}.

Because Ωμν=−iFˉμνaTadja\Omega_{\mu\nu}=-i\bar F^a_{\mu\nu}T^a_{\mathrm{adj}}, its color trace has a minus sign:

tr⁡adjΩμνΩμν=−Q.\operatorname{tr}_{\mathrm{adj}} \Omega_{\mu\nu}\Omega_{\mu\nu}=-\mathcal Q.

For the vector operator, Eμν=2ΩμνE_{\mu\nu}=2\Omega_{\mu\nu}. The connection acts identically on each of the four vector components, but EE mixes them. Therefore

tr⁡1,adjΩ2=−4Q,tr⁡1,adjE2=4∑μ,νtr⁡adjΩμνΩνμ=4Q.\begin{aligned} \operatorname{tr}_{1,\mathrm{adj}}\Omega^2&=-4\mathcal Q,\\ \operatorname{tr}_{1,\mathrm{adj}}E^2 &=4\sum_{\mu,\nu}\operatorname{tr}_{\mathrm{adj}} \Omega_{\mu\nu}\Omega_{\nu\mu} =4\mathcal Q. \end{aligned}

The second sign uses Ωνμ=−Ωμν\Omega_{\nu\mu}=-\Omega_{\mu\nu}; treating E2E^2 as four unrelated scalar squares would miss it. For the ghost, E=0E=0 and only the color trace remains. After the total derivative is integrated out, denote the coefficient without (4π)−2(4\pi)^{-2} by h4h_4:

h4vec=−4Q12+4Q2=53Q,h4gh=−Q12,H4≡12h4vec−h4gh=1112Q.\begin{aligned} h_4^{\mathrm{vec}} &=-{4\mathcal Q\over12}+{4\mathcal Q\over2} ={5\over3}\mathcal Q,\\ h_4^{\mathrm{gh}}&=-{\mathcal Q\over12},\\ H_4\equiv {1\over2}h_4^{\mathrm{vec}}-h_4^{\mathrm{gh}} &={11\over12}\mathcal Q. \end{aligned}

This is the decisive vector-and-ghost trace calculation. Vassilevich 2003, §4.2.1, pp.40–41, Eqs.(4.30)–(4.34), preprint PDF obtains the same result, but defines his integrated total heat coefficient as a4vec−2a4gha_4^{\mathrm{vec}}-2a_4^{\mathrm{gh}}. That convention is twice (16π2)−1∫d4x H4(16\pi^2)^{-1}\int d^4x\,H_4, because our coefficient includes the real-vector Gaussian weight.

The proper-time identity contributes a further minus sign:

ΓUV(1)=−∫0sIRdss(12Tr⁡e−sLvec−Tr⁡e−sLgh).\Gamma^{(1)}_{\mathrm{UV}} =-\int_0^{s_{\mathrm{IR}}}{ds\over s} \left({1\over2}\operatorname{Tr}e^{-s\mathcal L_{\mathrm{vec}}} -\operatorname{Tr}e^{-s\mathcal L_{\mathrm{gh}}}\right).

This expression denotes its ultraviolet part, with a regulator at the lower endpoint and background-independent terms subtracted. In d=4−2ϵd=4-2\epsilon, the F2F^2 term contains ∫0sIRds s−1+ϵ\int_0^{s_{\mathrm{IR}}}ds\,s^{-1+\epsilon}, whose pole is 1/ϵ1/\epsilon. Four vector components suffice for this pole; their O(ϵ)O(\epsilon) correction contributes only to finite terms. Hence

Γdiv(1)=−μ−2ϵ16π2ϵ11CA12∫ddx FˉμνaFˉμνa.\Gamma^{(1)}_{\mathrm{div}} =-{\mu^{-2\epsilon}\over16\pi^2\epsilon} {11C_A\over12} \int d^dx\,\bar F^a_{\mu\nu}\bar F^a_{\mu\nu}.

Canceling it with the bare kinetic term requires

1g02=μ−2ϵ[1g2+b016π2ϵ+⋯ ],b0=113CA.{1\over g_0^2} =\mu^{-2\epsilon} \left[{1\over g^2}+{b_0\over16\pi^2\epsilon}+\cdots\right], \qquad b_0={11\over3}C_A.

The factor four comes from the normalization F2/(4g02)F^2/(4g_0^2). Differentiate at fixed g0g_0, retaining the dimensional part of βd=μ dg/dμ\beta_d=\mu\,dg/d\mu:

0=−2ϵg2−2b016π2−2βdg3+O(g2),βd=−ϵg−b0g316π2+O(g5).\begin{aligned} 0&=-{2\epsilon\over g^2}-{2b_0\over16\pi^2} -{2\beta_d\over g^3}+O(g^2),\\ \beta_d&=-\epsilon g-{b_0g^3\over16\pi^2}+O(g^5). \end{aligned}

Taking ϵ→0\epsilon\to0 yields the advertised negative Yang–Mills beta function. As a diagrammatic check, the separate background-gauge ghost and vector contributions give CA/3C_A/3 and 10CA/310C_A/3 to b0b_0 in Schwartz 2014, §34.3.2, pp.756–757, Eqs.(34.96), (34.101). His dimensional parameter is ϵS=2ϵ\epsilon_{\mathrm S}=2\epsilon because he uses d=4−ϵSd=4-\epsilon_{\mathrm S}; the pole and fixed-bare calculation above specify our normalization directly.

Spin-one paramagnetism and the coefficient eleven-thirds

Section titled “Spin-one paramagnetism and the coefficient eleven-thirds”

The operator

−Dˉ2δμν-\bar D^2\delta_{\mu\nu}

is the orbital part. If this were the whole story, a gauge boson would resemble several charged scalar fields in the adjoint representation and would screen the charge. The crucial term is

−2ad⁡(Fˉμν).-2\operatorname{ad}(\bar F_{\mu\nu}).

This is the spin-one magnetic-moment coupling. It is the vector analog of the Pauli term in the squared Dirac operator. Its effect has the opposite sign and is larger.

A clean comparison uses one quantity throughout. Write the one-loop contribution of a field species as

β(g)=c g316π2+O(g5).\beta(g)=c\,{g^3\over16\pi^2}+O(g^5).

For one complex scalar, c/TR=+1/3c/T_R=+1/3. For one Dirac fermion, the orbital and Pauli contributions to c/TRc/T_R give

−23+2=43.-{2\over3}+2={4\over3}.

For the pure Yang–Mills gauge-and-ghost sector, the orbital terms in H4H_4 are 12(−Q/3)−(−Q/12)=−Q/12\tfrac12(-\mathcal Q/3)-(-\mathcal Q/12)=-\mathcal Q/12, while the vector spin term is 12(2Q)=Q\tfrac12(2\mathcal Q)=\mathcal Q. Since c=−b0c=-b_0 and b0b_0 is four times the coefficient of FˉaFˉa\bar F^a\bar F^a in H4H_4, this gives

cCA=13⏟orbital+(−4)⏟spin-one moment=−113.{c\over C_A}=\underbrace{{1\over3}}_{\text{orbital}} +\underbrace{(-4)}_{\text{spin-one moment}} =-{11\over3}.

Multiplying by 2g2g converts this to the coefficient in the flow of g2g^2:

dg2dlog⁡μ=2gβ(g),{d g^2\over d\log\mu}=2g\beta(g),

so the same pure-gauge decomposition becomes

23−8=−223.{2\over3}-8=-{22\over3}.

Thus

βpure YM(g)=−g316π2113CA+O(g5).\boxed{ \beta_{\mathrm{pure\ YM}}(g) =-{g^3\over16\pi^2}{11\over3}C_A+O(g^5). }

Ghosts are essential in the actual calculation. They do not merely “subtract two polarizations” in a naive way; rather, they enforce gauge invariance and cancel the unphysical pieces of the vector determinant. A quick check on any derivation is that the final logarithmic divergence must be proportional to the background-gauge-invariant operator ∫Fˉ2\int \bar F^2, with no leftover gauge-parameter dependence. But after the dust settles, the physical mnemonic is reliable:

spin-one paramagnetism dominates orbital screening.\text{spin-one paramagnetism dominates orbital screening.}

Orbital and spin contributions to one-loop beta-function signs

Schematic orbital and spin contributions to β(g)=c g3/(16π2)\beta(g)=c\,g^3/(16\pi^2) in the background Feynman-gauge calculation. Each row shows the coefficient after its group factor has been removed: c/TRc/T_R for matter and c/CAc/C_A for Yang–Mills. The labeled numbers, including the row totals, use that normalization; bar lengths are approximate. Scalar and spinor matter screen, while the gauge-and-ghost sum gives c=−11CA/3c=-11C_A/3.

The group factor is

facdfbcd=CAδab,f^{acd}f^{bcd}=C_A\delta^{ab},

so for SU(N)SU(N),

CA=N.C_A=N.

This is why the pure SU(N)SU(N) beta function is

β(g)=−11N3g316π2+O(g5).\beta(g)=-{11N\over3}{g^3\over16\pi^2}+O(g^5).

Combining gauge, fermion, and scalar loops gives

β(g)=−g316π2[113CA−43∑Dirac fT(Rf)−13∑complex sT(Rs)]+O(g5).\boxed{ \beta(g)=-{g^3\over16\pi^2} \left[ {11\over3}C_A -{4\over3}\sum_{\mathrm{Dirac\ f}}T(R_f) -{1\over3}\sum_{\mathrm{complex\ s}}T(R_s) \right]+O(g^5). }

Equivalently, for Weyl fermions and real scalars,

b0=113CA−23∑Weyl fT(Rf)−16∑real sT(Rs).\boxed{ b_0={11\over3}C_A -{2\over3}\sum_{\mathrm{Weyl\ f}}T(R_f) -{1\over6}\sum_{\mathrm{real\ s}}T(R_s). }

For a chiral matter list, this coefficient presupposes a consistent anomaly-free gauge theory; the numerical beta coefficient alone does not check gauge-anomaly cancellation.

For SU(N)SU(N) with NfN_f Dirac fermions in the fundamental representation,

TF=12,CA=N,T_F={1\over2}, \qquad C_A=N,

so

b0=113N−23Nf.\boxed{ b_0={11\over3}N-{2\over3}N_f. }

The condition for one-loop asymptotic freedom is

b0>0.b_0>0.

For QCD with gauge group SU(3)SU(3), this condition is

11−23Nf>0,Nf<332.11-{2\over3}N_f>0, \qquad N_f<{33\over2}.

Thus a small enough number of quark flavors preserves asymptotic freedom. Matter screens, but the gauge field antiscreens more strongly.

Because NfN_f is an integer, an SU(3)SU(3) theory with fundamental Dirac fermions is asymptotically free at one loop for Nf≤16N_f\leq16.

These formulas assume that the matter fields are light compared with the renormalization scale. If a field has mass MM, it contributes to the running above MM and decouples from low-energy logarithms below MM after matching. In a mass-independent subtraction scheme this decoupling is not automatic in the beta function; it is implemented by an effective field theory with threshold matching.

Gauge, fermion, and scalar contributions to the one-loop non-Abelian coefficient

The one-loop coefficient b0b_0 receives a positive contribution from gauge bosons and negative contributions from matter. In the convention β(g)=−b0g3/(16π2)\beta(g)=-b_0g^3/(16\pi^2), positive b0b_0 means asymptotic freedom.

The important conceptual point is that the sign of a beta function is not determined merely by whether particles are bosons or fermions. It is determined by the full fluctuation operator: statistics, number of degrees of freedom, representation under the gauge group, and spin coupling to the background field all matter.

For reference, in the convention above:

field in representation Rcontribution to b0Dirac fermion−43T(R)Weyl fermion−23T(R)complex scalar−13T(R)real scalar−16T(R)\begin{array}{c|c} \text{field in representation }R & \text{contribution to }b_0 \\ \hline \text{Dirac fermion} & -{4\over3}T(R) \\ \text{Weyl fermion} & -{2\over3}T(R) \\ \text{complex scalar} & -{1\over3}T(R) \\ \text{real scalar} & -{1\over6}T(R) \end{array}

This small table is a useful guardrail when moving between particle-physics, supersymmetry, and statistical-field-theory normalizations.

Supersymmetric cancellations and the special role of ten dimensions

Section titled “Supersymmetric cancellations and the special role of ten dimensions”

The coefficients above give quick checks of supersymmetric gauge theories. In four-dimensional N=4\mathcal N=4 super-Yang–Mills theory all fields are in the adjoint representation. The field content can be written as

one gauge boson+4 Weyl fermions+6 real scalars.\text{one gauge boson} \quad+ 4\ \text{Weyl fermions} \quad+ 6\ \text{real scalars}.

Using T(adj)=CAT(\mathrm{adj})=C_A, the one-loop coefficient is

b0=113CA−23(4CA)−16(6CA).b_0 ={11\over3}C_A -{2\over3}(4C_A) -{1\over6}(6C_A).

Therefore

b0=(113−83−1)CA=0.\boxed{ b_0=\left({11\over3}-{8\over3}-1\right)C_A=0. }

This cancellation is not an accident of arithmetic. Four-dimensional N=4\mathcal N=4 super-Yang–Mills can be obtained by dimensional reduction of ten-dimensional minimal supersymmetric Yang–Mills theory. Here we have checked its one-loop gauge beta function; an all-orders conclusion requires additional supersymmetric arguments.

One-loop cancellation in N=4 super-Yang-Mills theory

In N=4\mathcal N=4 super-Yang–Mills theory, the gauge-boson contribution to b0b_0 is exactly canceled by four adjoint Weyl fermions and six adjoint real scalars: 11/3−8/3−1=011/3-8/3-1=0.

A related mnemonic is that the one-loop coefficient for dimensionally reduced supersymmetric Yang–Mills theories knows about the critical dimension D=10D=10. The four-dimensional N=4\mathcal N=4 theory inherits precisely the field content of ten-dimensional minimal SYM, and the one-loop gauge beta function vanishes. This remark should not be confused with a proof of all-order finiteness; it is only the one-loop statement visible from the beta-function coefficients on this page.

For a theory with

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

we have

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

Integrating between μ0\mu_0 and μ\mu gives

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

Thus g(μ)g(\mu) becomes small when μ\mu becomes large. Perturbation theory improves in the ultraviolet.

The same equation can be written in terms of the RG-invariant scale

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

Then

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

at one loop. This formula is trustworthy only when μ≫Λ\mu\gg\Lambda. As μ\mu approaches Λ\Lambda, perturbation theory breaks down. The next page develops this phenomenon as dimensional transmutation: the classical theory has a dimensionless coupling, but the quantum theory produces a scale.

The beta function of a gauge theory is encoded in the logarithmic renormalization of the background F2F^2 term. In the rescaled normalization, matter loops correct 1/g21/g^2 directly.

For QED matter,

β(e)=e316π2(43Nf+13Ns)+O(e5),\beta(e)={e^3\over16\pi^2} \left({4\over3}N_f+{1\over3}N_s\right)+O(e^5),

so the Abelian charge is screened in the infrared and grows toward the ultraviolet.

For pure Yang–Mills theory,

β(g)=−g316π2113CA+O(g5),\beta(g)=-{g^3\over16\pi^2}{11\over3}C_A+O(g^5),

so the coupling decreases at short distances. The physical origin of the sign is spin-one paramagnetism: the gauge boson magnetic-moment term dominates the orbital screening part. Ghosts are required for the gauge-invariant coefficient.

With matter,

b0=113CA−43∑Dirac fT(Rf)−13∑complex sT(Rs).b_0={11\over3}C_A -{4\over3}\sum_{\mathrm{Dirac\ f}}T(R_f) -{1\over3}\sum_{\mathrm{complex\ s}}T(R_s).

Matter reduces b0b_0. A non-Abelian theory is asymptotically free when the gauge-boson term wins.

Inferring the sign from statistics alone. The sign comes from the full quadratic operator in the background field. Spin couplings are decisive.

Dropping the ghost determinant. The background-field vector determinant by itself contains unphysical gauge modes. The ghost determinant is part of the gauge-fixed definition of the theory and is necessary for the coefficient 11/311/3.

Mixing the flows of gg and g2g^2. If

β(g)=−113CAg316π2,\beta(g)=-{11\over3}{C_Ag^3\over16\pi^2},

then

dg2dlog⁡μ=−223CAg416π2.{d g^2\over d\log\mu}=-{22\over3}{C_Ag^4\over16\pi^2}.

The extra factor of two is a common source of mismatched coefficients.

Taking “antiscreening” too literally. A color-charge cloud is not by itself a gauge-invariant observable. The coefficient derived here governs perturbative short-distance running and can enter predictions for specified physical observables. Background gauge covariance alone does not identify the full off-shell effective action with such an observable, or establish the stability of a constant chromomagnetic field.

Exercise 1: Derive the QED beta function from the inverse coupling

Section titled “Exercise 1: Derive the QED beta function from the inverse coupling”

Starting from

1e2(q)=1e2(μ)+b16π2log⁡μ2q2,{1\over e^2(q)}={1\over e^2(\mu)}+{b\over16\pi^2}\log{\mu^2\over q^2},

derive

β(e)=b16π2e3.\beta(e)={b\over16\pi^2}e^3.
Solution

Hold the physical momentum qq fixed and differentiate the stated relation with respect to log⁡μ\log\mu. The effective coupling on the left is independent of the arbitrary subtraction point, so

0=−2e3β(e)+2b16π2.0=-{2\over e^3}\beta(e)+{2b\over16\pi^2}.

Solving gives

β(e)=b16π2e3.\boxed{ \beta(e)={b\over16\pi^2}e^3. }

Exercise 2: Translate the Yang–Mills flow from the coupling to its square

Section titled “Exercise 2: Translate the Yang–Mills flow from the coupling to its square”

For pure Yang–Mills theory,

β(g)=−113CAg316π2.\beta(g)=-{11\over3}{C_Ag^3\over16\pi^2}.

Find dg2/dlog⁡μd g^2/d\log\mu and explain why the coefficient is −22CA/3-22C_A/3 rather than −11CA/3-11C_A/3.

Solution

Use

dg2dlog⁡μ=2gdgdlog⁡μ=2gβ(g).{d g^2\over d\log\mu}=2g{dg\over d\log\mu}=2g\beta(g).

Substituting the beta function,

dg2dlog⁡μ=2g[−113CAg316π2]=−223CAg416π2.{d g^2\over d\log\mu} =2g\left[-{11\over3}{C_Ag^3\over16\pi^2}\right] =-{22\over3}{C_Ag^4\over16\pi^2}.

The coefficient doubles because g2g^2 changes twice as fast logarithmically as gg:

dlog⁡g2dlog⁡μ=2dlog⁡gdlog⁡μ.{d\log g^2\over d\log\mu}=2{d\log g\over d\log\mu}.

Exercise 3: Find the asymptotic-freedom window for fundamental matter

Section titled “Exercise 3: Find the asymptotic-freedom window for fundamental matter”

For SU(N)SU(N) gauge theory with NfN_f Dirac fermions in the fundamental representation and no scalars, show that

b0=113N−23Nf.b_0={11\over3}N-{2\over3}N_f.

For what values of NfN_f is the theory asymptotically free?

Solution

For SU(N)SU(N),

CA=N,TF=12C_A=N, \qquad T_F={1\over2}

in the fundamental representation. The general formula with NfN_f Dirac fermions is

b0=113CA−43TFNf.b_0={11\over3}C_A-{4\over3}T_FN_f.

Substituting the group factors gives

b0=113N−4312Nf=113N−23Nf.b_0={11\over3}N-{4\over3}{1\over2}N_f ={11\over3}N-{2\over3}N_f.

Asymptotic freedom requires b0>0b_0>0, so

113N−23Nf>0.{11\over3}N-{2\over3}N_f>0.

Multiplying by 33 gives

11N−2Nf>0,11N-2N_f>0,

hence

Nf<112N.\boxed{N_f<{11\over2}N.}

For SU(3)SU(3) this becomes Nf<33/2N_f<33/2.

Exercise 4: Check the maximally supersymmetric one-loop cancellation

Section titled “Exercise 4: Check the maximally supersymmetric one-loop cancellation”

Check the one-loop cancellation in N=4\mathcal N=4 super-Yang–Mills theory. Use one gauge boson, four Weyl fermions, and six real scalars, all in the adjoint representation.

Solution

For adjoint fields,

T(adj)=CA.T(\mathrm{adj})=C_A.

The gauge boson contributes

113CA.{11\over3}C_A.

Each Weyl fermion contributes

−23CA,-{2\over3}C_A,

so four Weyl fermions contribute

4(−23CA)=−83CA.4\left(-{2\over3}C_A\right)=-{8\over3}C_A.

Each real scalar contributes

−16CA,-{1\over6}C_A,

so six real scalars contribute

6(−16CA)=−CA.6\left(-{1\over6}C_A\right)=-C_A.

The total coefficient is

b0=113CA−83CA−CA=(113−83−33)CA=0.b_0={11\over3}C_A-{8\over3}C_A-C_A =\left({11\over3}-{8\over3}-{3\over3}\right)C_A=0.

Thus the one-loop beta function vanishes.

Exercise 5: Construct the one-loop RG-invariant scale

Section titled “Exercise 5: Construct the one-loop RG-invariant 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 the flow of the inverse coupling:

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

Now take the logarithm of Λ\Lambda:

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

Differentiate:

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

Using the inverse-coupling flow,

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

Therefore Λ\Lambda is RG invariant at one loop.

Exercise 6: Test the ghost weight in the heat coefficient

Section titled “Exercise 6: Test the ghost weight in the heat coefficient”

Use h4vec=5Q/3h_4^{\mathrm{vec}}=5\mathcal Q/3 and h4gh=−Q/12h_4^{\mathrm{gh}}=-\mathcal Q/12 to find the inferred b0b_0 if the ghost determinant is (a) omitted or (b) mistakenly assigned weight −1/2-1/2. Explain why neither reproduces Yang–Mills theory.

Solution

Write H4=QhH_4=\mathcal Q h. The kinetic normalization gives b0=4CAhb_0=4C_Ah. Without the ghost,

h=1253=56,b0=103CA.h={1\over2}{5\over3}={5\over6}, \qquad b_0={10\over3}C_A.

With half the required ghost weight,

h=56+124=78,b0=72CA.h={5\over6}+{1\over24}={7\over8}, \qquad b_0={7\over2}C_A.

The actual Grassmann pair has weight −1-1, giving h=5/6+1/12=11/12h=5/6+1/12=11/12 and b0=11CA/3b_0=11C_A/3. The ghost operator acts on adjoint scalar components, but its determinant weight follows from Grassmann integration; it cannot be chosen by counting transverse polarizations in the unrestricted vector operator.

  • Schwartz, Matthew D. Quantum Field Theory and the Standard Model. Cambridge University Press, 2014. DOI. The diagrammatic comparison uses §34.3.2, Eqs.(34.96) and (34.101), with the dimensional-parameter translation stated above.
  • Vassilevich, Dmitri V. “Heat Kernel Expansion: User’s Manual.” Physics Reports 388 (2003): 279–360. DOI. Open PDF, arXiv:hep-th/0306138v3. All page locators above refer to this preprint’s printed pagination.
  • Abbott, L. F. “Introduction to the Background Field Method.” Acta Physica Polonica B 13, no. 1 (1982): 33–50.
  • 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.
  • Srednicki, Mark. Quantum Field Theory. Cambridge: Cambridge University Press, 2007, Sections 66, 73, and 78.
  • Weinberg, Steven. The Quantum Theory of Fields. Volume II: Modern Applications. Cambridge: Cambridge University Press, 1996, Chapters 17–18.
  • 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.