Skip to content

Static Potentials, Flux Tubes, and Effective Strings

A rectangular Wilson loop is a transfer-matrix correlator for a static source–antisource pair. Taking its Euclidean-time extent to infinity at fixed separation extracts the lowest static energy; only after that projection may one test linear growth, string breaking, or the universal long-string correction. A visible tube at finite time or distance can instead be a metastable excitation.

Required background. Line operators, screening, and symmetry diagnostics fixes the genuine probe and predicts whether dynamical matter can break its string.

Helpful background. Background-field Yang–Mills supplies the gauge-covariant expansion used in effective descriptions.

Shared comparison. The diagnostic map distinguishes a static potential from the other confinement claims, and the claim–evidence table records the relevant failure limits.

Rectangular loops project onto static energies

Section titled “Rectangular loops project onto static energies”

Work at zero temperature in dd-dimensional Euclidean spacetime. Let WR(R,TE)W_R(R,T_E) be a renormalized rectangular loop of spatial width RR and Euclidean-time height TET_E in a genuine representation RR. Cutting the rectangle across a time slice prepares a gauge-invariant state containing two static sources. The transfer matrix gives

WR(R,TE)=ncn(R)2eEn(R)TE,\langle W_R(R,T_E)\rangle =\sum_{n}|c_n(R)|^2e^{-E_n(R)T_E},

where cnc_n is the overlap of that Wilson operator with the nnth state in the static-source sector. Therefore

VR(R)=E0(R)2δmR=limTE1TElogWR(R,TE)ren.V_R(R)=E_0(R)-2\delta m_R =-\lim_{T_E\to\infty}\frac{1}{T_E} \log \langle W_R(R,T_E)\rangle_{\mathrm{ren}}.

The subtraction 2δmR2\delta m_R removes the scheme-dependent static self-energies, equivalently the temporal part of the loop’s perimeter divergence. Additive constants in VRV_R are not physical; forces, level splittings, and the coefficient of asymptotic growth are.

If a large loop obeys

logW(R,TE)=σRTE+μ(R+TE)+,-\log\langle W(R,T_E)\rangle =\sigma RT_E+\mu(R+T_E)+\cdots,

then taking TET_E\to\infty yields V(R)=σR+μ+V(R)=\sigma R+\mu+\cdots. This is the precise bridge from an area law to a linear static energy, already implicit in Wilson’s strong-coupling loop analysis Wilson 1974, §§IV–V, pp. 2450–2456. The converse requires the same loop family and compatible limits; it cannot be inferred from a finite interval that merely looks straight.

Suppose the chosen probe cannot be screened, the ground state contains a stable long tube, and

Rσ1.R\sqrt{\sigma}\gg1.

At wavelengths much longer than the tube thickness, the only mandatory gapless world-sheet fields are its d2d-2 transverse displacements Xi(t,z)X^i(t,z). To quadratic order in static gauge,

Sws=σRTE+σ2i=1d2dtdzαXiαXi+,S_{\mathrm{ws}} =\sigma RT_E +\frac{\sigma}{2}\sum_{i=1}^{d-2} \int\mathrm d t\,\mathrm d z\, \partial_\alpha X^i\partial_\alpha X^i+\cdots,

with Dirichlet endpoints at the static sources. Each transverse scalar has modes ωn=πn/R\omega_n=\pi n/R. Its regularized zero-point energy is

12n=1ωn=π2Rζ(1)=π24R.\frac12\sum_{n=1}^{\infty}\omega_n =\frac{\pi}{2R}\,\zeta(-1) =-\frac{\pi}{24R}.

Summing over transverse directions gives the universal open-string asymptotics

V(R)=σR+μπ(d2)24R+O(R3).V(R)=\sigma R+\mu -\frac{\pi(d-2)}{24R} +O(R^{-3}).

The negative 1/R1/R coefficient is the Lüscher term. It depends only on d2d-2 massless transverse modes and the open-string boundary conditions, not on the microscopic flux profile Lüscher, Symanzik, and Weisz 1980, pp. 365–396. Extra massless world-sheet fields or different boundaries change the coefficient; short strings need not follow the expansion.

The same effective theory predicts excitation gaps at leading order,

En(R)E0(R)=πnR+O(R3),E_n(R)-E_0(R)=\frac{\pi n}{R}+O(R^{-3}),

organized by transverse oscillator quantum numbers. Agreement of both the ground-state correction and the excitation pattern tests a string description more sharply than observing approximate linearity alone.

String breaking changes the ground-state branch

Section titled “String breaking changes the ground-state branch”

With dynamical matter capable of screening the probe, compare a flux-tube state and a two-hadron state,

ES(R)2δm+σR,EM(R)2MQq.E_{S}(R)\simeq2\delta m+\sigma R, \qquad E_{M}(R)\simeq2M_{Qq}.

In the two-state approximation,

H(R)=(ES(R)x(R)x(R)EM(R)),H(R)= \begin{pmatrix} E_S(R)&x(R)\\ x(R)&E_M(R) \end{pmatrix},

so the eigenvalues repel. The lower eigenvalue saturates rather than remaining linear. Yet a thin Wilson loop often has cScM|c_S|\gg|c_M|; at practical TET_E it can track the upper, string-like eigenstate past the avoided crossing. Mixed correlators containing both string and two-hadron operators reveal the true ground state, as demonstrated in unquenched lattice QCD by Bali et al. 2005, §§II–IV.

Thus the order of limits is part of the claim:

project TE at fixed Rbefore testing R.\boxed{\text{project }T_E\to\infty\text{ at fixed }R \quad\text{before testing }R\to\infty.}

When the string breaks, the long-string expansion remains meaningful for a metastable level only over distances where that resonance is narrow and the world-sheet scale separation survives. It is no longer an asymptotic order parameter.

Flux profiles and discriminating measurements

Section titled “Flux profiles and discriminating measurements”

A gauge-invariant flux profile can be defined by correlating the Wilson loop with a local field-strength composite, subtracting its vacuum expectation value, and taking the same large-TET_E projection. Useful observables include the transverse energy-density profile, its width, and excitation energies. They answer different questions:

  • V(R)V'(R) measures the force;
  • the profile asks whether energy is spatially collimated;
  • the spectrum tests the world-sheet degrees of freedom;
  • mixed-operator overlaps locate string breaking.

No single one of these measurements identifies a microscopic mechanism. A dual-superconductor model, a center-vortex description, and a controlled monopole plasma can share a linear potential while predicting different field profiles or defect responses.

1. Excited-state contamination. Suppose W(T)=a0eE0T+a1eE1TW(T)=a_0e^{-E_0T}+a_1e^{-E_1T} with E1>E0E_1>E_0. Show that the effective energy Eeff(T)=TlogW(T)E_{\mathrm{eff}}(T)=-\partial_T\log W(T) approaches E0E_0 from above.

Solution

Writing r=(a1/a0)e(E1E0)Tr=(a_1/a_0)e^{-(E_1-E_0)T} gives Eeff=(E0+rE1)/(1+r)=E0+(E1E0)r/(1+r)E_{\mathrm{eff}}=(E_0+rE_1)/(1+r)=E_0+(E_1-E_0)r/(1+r). For positive spectral weights, the correction is positive and vanishes exponentially. Poor ground-state overlap means rr can remain large for a long Euclidean time.

2. Lüscher coefficient. Evaluate the universal term in d=3d=3 and d=4d=4 spacetime dimensions.

Solution

There are respectively one and two transverse fields. Hence the terms are π/(24R)-\pi/(24R) in d=3d=3 and π/(12R)-\pi/(12R) in d=4d=4.

  • Bali, Gunnar S., Hartmut Neff, Thomas Düsel, Thomas Lippert, and Klaus Schilling. “Observation of String Breaking in QCD.” Physical Review D 71 (2005): 114513. DOI. Open PDF.
  • Lüscher, Martin, Kurt Symanzik, and Peter Weisz. “Anomalies of the Free Loop Wave Equation in the WKB Approximation.” Nuclear Physics B 173 (1980): 365–396. DOI.
  • Wilson, Kenneth G. “Confinement of Quarks.” Physical Review D 10 (1974): 2445–2459. DOI.