Skip to content

Vacuum Polarization, Running Charge, and Screening

Vacuum polarization inserts a charged-fermion loop into the photon propagator. Current conservation makes the insertion transverse; renormalization turns its scalar coefficient into either a momentum-dependent effective charge or a running scheme parameter. Both descriptions encode screening, but they are not numerically identical until their subtraction conventions are matched.

Required background. Electron and photon renormalization supplies Z3Z_3 and the subtraction conditions. Beta functions and anomalous dimensions supplies the general renormalization-group definition.

Helpful background. Subtracted dispersion relations explains how the timelike spectral cut determines the spacelike function.

For one Dirac fermion of mass mm and charge magnitude ee, the one-loop proper photon two-point function is

iΠμν(q)=(ie)2μ2ϵdd(2π)dtr[γμi( ⁣ ⁣ ⁣/+m)2m2+i0γνi( ⁣ ⁣ ⁣/+q ⁣ ⁣ ⁣/+m)(+q)2m2+i0].i\Pi^{\mu\nu}(q) =-(-ie)^2\mu^{2\epsilon} \int\frac{\mathrm d^d\ell}{(2\pi)^d} \operatorname{tr}\left[ \gamma^\mu\frac{i(\ell\!\!\!/+m)}{\ell^2-m^2+i0} \gamma^\nu\frac{i(\ell\!\!\!/+q\!\!\!/+m)}{(\ell+q)^2-m^2+i0} \right].

The initial minus sign is the closed-fermion-loop sign. A symmetry-preserving regulator and consistent momentum routing give

qμΠμν(q)=0,q_\mu\Pi^{\mu\nu}(q)=0,

so Lorentz invariance fixes a single scalar function. Choose the positive-spectral current-correlator convention

Πμν(q)=(qμqνq2ημν)ΠM(q2).\Pi^{\mu\nu}(q) =(q^\mu q^\nu-q^2\eta^{\mu\nu})\Pi_M(q^2).

Any residual ημν\eta^{\mu\nu} term not accompanied by the required q2q^2 is a photon-mass artifact and violates the Ward identity. The loop calculation and its Coulomb-potential interpretation are developed in Schwartz 2014, § 16.2, pp. 304–314 and Weinberg 1995, § 11.2, pp. 473–482.

Let Q2=q2>0Q^2=-q^2>0 be spacelike. To make the sign convention explicit, define the positive Euclidean screening function

ΠE(Q2)[ΠM(Q2)ΠM(0)].\Pi_E(Q^2) \equiv-\bigl[\Pi_M(-Q^2)-\Pi_M(0)\bigr].

Zero-momentum subtraction then gives the finite one-loop result

ΠE(Q2)=2απ01dxx(1x)ln[1+Q2m2x(1x)]\boxed{ \Pi_E(Q^2)=\frac{2\alpha}{\pi} \int_0^1\mathrm dx\,x(1-x) \ln\left[1+\frac{Q^2}{m^2}x(1-x)\right]}

with α=e2/(4π)\alpha=e^2/(4\pi). Dyson summation of transverse insertions motivates the zero-momentum-subtracted Euclidean effective charge

αeff(Q2)=α1ΠE(Q2).\alpha_{\rm eff}(Q^2) =\frac{\alpha}{1-\Pi_E(Q^2)}.

This definition is tied to the chosen subtraction and spacelike kinematics. It should not be silently identified with an MS\overline{\mathrm{MS}} coupling.

At low momentum, expand the logarithm:

ΠE(Q2)=2απQ2m201dxx2(1x)2+O(Q4/m4)=α15πQ2m2+O(Q4/m4).\begin{aligned} \Pi_E(Q^2) &=\frac{2\alpha}{\pi}\frac{Q^2}{m^2} \int_0^1\mathrm dx\,x^2(1-x)^2 +O(Q^4/m^4)\\ &=\frac{\alpha}{15\pi}\frac{Q^2}{m^2} +O(Q^4/m^4). \end{aligned}

The Q2/m2Q^2/m^2 suppression is decoupling: a massive charged species cannot produce logarithmic running far below its pair threshold. At the independent checkpoint Q2=m2Q^2=m^2,

ΠE(m2)=α9π[7+5ln(945)].\Pi_E(m^2) =\frac{\alpha}{9\pi} \left[7+\sqrt5\ln(9-4\sqrt5)\right].

At Q2m2Q^2\gg m^2,

ΠE(Q2)=α3π[lnQ2m253]+O ⁣(m2Q2lnQ2m2).\Pi_E(Q^2) =\frac{\alpha}{3\pi} \left[\ln\frac{Q^2}{m^2}-\frac53\right] +O\!\left(\frac{m^2}{Q^2}\ln\frac{Q^2}{m^2}\right).

The positive logarithm makes αeff\alpha_{\rm eff} increase as shorter distances are resolved.

The Euclidean integral has no imaginary part. Continuing the underlying Minkowski scalar ΠM\Pi_M to timelike s=q2>0s=q^2>0 requires the +i0+i0 prescription and a declared logarithm branch. The first physical cut begins at

s=4m2,s=4m^2,

where the photon can create a real fermion–antifermion pair. With the scalar convention above, the spectral density is

ρ(s)=1πImΠM(s+i0)=α3π(1+2m2s)14m2sθ(s4m2).\rho(s)=\frac{1}{\pi}\operatorname{Im}\Pi_M(s+i0) =\frac{\alpha}{3\pi} \left(1+\frac{2m^2}{s}\right) \sqrt{1-\frac{4m^2}{s}}\, \theta(s-4m^2).

A once-subtracted dispersion relation reconstructs the spacelike function:

ΠE(Q2)=Q24m2dss(s+Q2)ρ(s).\Pi_E(Q^2) =Q^2\int_{4m^2}^{\infty} \frac{\mathrm ds}{s(s+Q^2)}\rho(s).

This formula supplies an independent sign check: ρ(s)0\rho(s)\ge0 implies ΠR(Q2)0\Pi_R(Q^2)\ge0 for Q2>0Q^2>0. Using the Euclidean logarithm directly above timelike threshold would miss both the imaginary part and the correct sheet.

In MS\overline{\mathrm{MS}}, the one-loop photon counterterm for NfN_f active unit-charge Dirac fermions is

Z3=1Nfα3πϵˉ+O(α2).Z_3=1-\frac{N_f\alpha}{3\pi\bar\epsilon}+O(\alpha^2).

Together with e0=μϵZ31/2ee_0=\mu^\epsilon Z_3^{-1/2}e, holding the bare coupling fixed gives

dαdlnμ=2Nf3πα2+O(α3).\boxed{ \frac{\mathrm d\alpha}{\mathrm d\ln\mu} =\frac{2N_f}{3\pi}\alpha^2+O(\alpha^3)}.

For charges QfQ_f in units of ee, replace NfN_f by fQf2\sum_f Q_f^2. At fixed NfN_f, integration gives the one-loop invariant

α1(μ)+2Nf3πlnμμref=constant,\alpha^{-1}(\mu) +\frac{2N_f}{3\pi}\ln\frac{\mu}{\mu_{\rm ref}} =\text{constant},

or

α1(μ)=α1(μ0)2Nf3πlnμμ0.\alpha^{-1}(\mu) =\alpha^{-1}(\mu_0) -\frac{2N_f}{3\pi}\ln\frac{\mu}{\mu_0}.

The factor is written for lnμ\ln\mu. Differentiating with respect to lnμ2\ln\mu^2 halves it; mixing the two conventions is a common error. The general counterterm derivation appears in Schwartz 2014, §§ 23.1–23.2, pp. 417–426.

A mass-independent scheme does not automatically decouple a heavy field. Below a mass MM, match the full theory onto an EFT without that active species and run with the lower NfN_f. At leading order the coupling is continuous at μ=M\mu=M; higher orders add a finite decoupling coefficient.

For the exact benchmark

Nf=1(μ<M),Nf=2(μ>M),α(M)=110,N_f=1\quad(\mu<M), \qquad N_f=2\quad(\mu>M), \qquad \alpha(M)=\frac{1}{10},

one-loop running with leading matching gives

α1(M/e)=10+23π,α1(eM)=1043π.\alpha^{-1}(M/e)=10+\frac{2}{3\pi}, \qquad \alpha^{-1}(eM)=10-\frac{4}{3\pi}.

Running through MM with one unchanged beta coefficient fails the threshold check.

To compare αMS(μ)\alpha_{\overline{\rm MS}}(\mu) with αeff(Q2)\alpha_{\rm eff}(Q^2), choose a matching point, compute the finite conversion there, and only then compare slopes. Their raw values belong to different definitions.

Virtual charged pairs polarize the vacuum so that a distant probe sees a more screened charge than a short-distance probe. Equivalently, the beta function is positive and the renormalized coupling grows toward the ultraviolet. This interpretation is gauge invariant when formulated through the current–current amplitude or a physical effective charge, not through a gauge-dependent photon field normalization alone.

Formally extending the one-loop solution produces a scale where its denominator vanishes. That “Landau pole” lies outside the regime in which the one-loop approximation can justify itself. It signals that perturbative QED does not supply a controlled ultraviolet completion; it is not a prediction of a physical pole at a trustworthy energy.

An MS\overline{\mathrm{MS}} coupling is not a measured potential by definition. Match it to a physical effective charge before comparing numbers.

Euclidean and timelike polarization are not the same real function. Above 4m24m^2, the +i0+i0 prescription and branch cut are essential.

Mass-independent running does not implement decoupling automatically. Change the active theory and match at each threshold.

The one-loop Landau pole is not controlled evidence of a state. It occurs where the approximation used to derive it has already failed.

  • Matthew D. Schwartz, Quantum Field Theory and the Standard Model, Cambridge University Press (2014), §§ 16.2 and 23.1–23.2, doi:10.1017/9781139540940.
  • Steven Weinberg, The Quantum Theory of Fields, Volume I: Foundations, Cambridge University Press (1995), § 11.2, doi:10.1017/CBO9781139644167.