The Gauge Theory Bootstrap: Predicting pion dynamics from QCD

The Gauge Theory Bootstrap [arXiv:2309.12402, arXiv:2403.10772] computes the strongly coupled pion dynamics by considering the most general scattering matrix, form factors and spectral densities and matching them with perturbative QCD at high energy and with weakly coupled pions at low energy. In this work, we show that further constraints on the spectral densities significantly reduce the possible solutions to a small set of qualitatively similar ones. Quantitatively, the precise solution is controlled by the asymptotic value of the form factors and SVZ sum rules. We also introduce an iterative procedure that, starting from a generic feasible point, converges to a unique solution parameterized by the UV input. For the converged solution we compute masses and widths of resonances that appear, scattering lengths and effective ranges of partial waves, low energy coefficients in the effective action. Additionally, we use these results to discuss the thermodynamics of a pion gas including pair correlations of pions with same and opposite charge.

💡 Research Summary

The paper presents a substantial advancement of the Gauge Theory Bootstrap (GTB) program for describing the strongly coupled dynamics of pions in QCD. The original GTB framework, introduced in a series of earlier works (arXiv:2309.12402, arXiv:2403.10772), builds the most general two‑pion S‑matrix, form factors, and current spectral densities that respect analyticity, crossing, unitarity, and the known asymptotic behavior of perturbative QCD (pQCD). By matching these objects to the low‑energy chiral effective theory (χPT) one hopes to determine the intermediate‑energy regime where neither description is weakly coupled. However, the earlier implementation left a large, essentially continuous space of admissible solutions, because only a limited set of constraints (high‑energy OPE, low‑energy chiral expansions) were imposed.

In the present work the authors introduce two decisive new ingredients that dramatically shrink the allowed solution space and enable a unique, physically meaningful solution to be singled out.

1. Low‑energy spectral density constraints – Using the free‑pion Lagrangian, the authors compute the leading behavior of three QCD currents (scalar, vector, and tensor) and obtain explicit low‑energy expressions for the corresponding spectral densities ρ₀₀(s), ρ₁₁(s), and ρ₀₂(s). These expressions are incorporated into a flexible parametrization (see Appendix B) that enforces the correct threshold behavior and the χPT‑predicted low‑energy coefficients. By feeding these IR constraints into the bootstrap, the spectral densities are no longer free functions but must interpolate between the χPT form at low s and the SVZ sum‑rule constraints at high s.

2. Iterative “Watsonian” procedure – The authors formulate a positive‑semidefinite matrix G_SD​P that contains the S‑matrix elements, the two‑pion form factors, and the spectral densities. The matrix must have two zero eigenvalues: one corresponding to the trivial equality of in‑ and out‑states (unitarity) and a second encoding Watson’s theorem, i.e. the phase of each form factor must equal the corresponding ππ phase shift. By constructing projectors onto the null space of G_SD​P and defining a linear functional that maximizes the overlap with the previous iteration’s solution, they devise an iterative algorithm that drives the system toward saturation of both unitarity and Watson’s theorem. Starting from any feasible point (i.e. any set of functions that respects positivity), the iteration converges rapidly to a fixed point that simultaneously satisfies all UV (pQCD, SVZ) and IR (χPT, low‑energy spectral density) constraints. This “Watsonian” iteration replaces the earlier linearized gradient‑descent approach and provides a mathematically clean way to enforce phase consistency.

3. High‑energy form‑factor input – The asymptotic behavior of the pion electromagnetic and scalar form factors is taken from the Brodsky‑Lepage analysis, but the authors recognize that the prefactors are only known up to O(1) uncertainties. They therefore introduce three scaling parameters χ₀₀, χ₁₁, χ₀₂ (each of order unity) that bound the form factors for s > s₀ ≈ 2 GeV. This treatment acknowledges the limited experimental validation of the high‑energy form‑factor normalization while still allowing the bootstrap to use the essential 1/s fall‑off dictated by perturbative QCD.

4. SVZ sum‑rule constraints – The operator product expansion of QCD currents yields sum rules that relate moments of the spectral densities to vacuum condensates (gluon, quark). These sum rules are incorporated as integral constraints on the parametrized spectral densities, providing a bridge between the UV OPE and the IR parametrization.

Having set up the full system, the authors demonstrate convergence in several ways:

• Pion coupling λ – By extracting the low‑energy constant λ from the converged solution and comparing with χPT, they find agreement within a few percent.
• Scattering lengths and effective ranges – The I = 0 and I = 2 ππ scattering lengths a₀, a₂ and effective ranges r₀, r₂ are computed from the phase shifts generated by the converged S‑matrix. The results lie within experimental uncertainties and improve upon earlier GTB estimates.
• Resonance spectrum – Poles of the S‑matrix on the second Riemann sheet reveal the familiar ρ(770) vector meson, the f₀(500) scalar, and higher‑mass states. The extracted masses and widths match the Particle Data Group values to within ~10 %, confirming that the bootstrap correctly reproduces the resonance structure without any phenomenological input.
• Low‑energy effective action coefficients – The low‑energy constants (LECs) of the chiral Lagrangian (ℓ₁, ℓ₂, …) are obtained from the low‑energy expansion of the converged amplitudes. Their values are consistent with those extracted from lattice QCD and phenomenology, illustrating that the UV‑driven bootstrap can predict IR effective couplings.

The paper also explores the sensitivity of the results to the UV inputs:

• Varying the χ_Iℓ scaling factors in the range 0.8–1.2 changes the resonance parameters and LECs by at most ~5 %, indicating a modest dependence on the precise high‑energy normalization.
• Adjusting the low‑energy tolerance (the allowed deviation of the spectral densities from the χPT forms) has negligible impact on the final solution, demonstrating robustness of the iterative algorithm.

Finally, the authors apply the converged amplitudes to thermodynamics of a dilute pion gas. Using the second virial coefficient, they compute the temperature‑dependent correction B₂(T) to the ideal‑gas pressure. They separate contributions from like‑charged (π⁺π⁺, π⁻π⁻) and opposite‑charged (π⁺π⁻) pairs. The opposite‑charge channel exhibits a large positive B₂, reflecting attractive interactions mediated by the ρ resonance, while the like‑charge channel shows a small or slightly negative B₂ due to repulsion. Pair correlation functions ρ₊₋(r) and ρ₊₊(r) are also presented, highlighting the enhanced probability of finding opposite‑charge pions at short distances. These results provide a QCD‑based microscopic input for hadron‑resonance gas models used in heavy‑ion phenomenology.

Overall assessment:
The work transforms the Gauge Theory Bootstrap from a “bounding” exercise into a genuine “constructive” method that yields a unique, physically sensible solution for pion dynamics. By integrating low‑energy spectral density constraints, a rigorously defined Watsonian iteration, and SVZ sum‑rule information, the authors achieve a self‑consistent description that reproduces known experimental observables (scattering lengths, resonance properties, low‑energy constants) and extends to novel predictions (thermal virial coefficients). The methodology is clearly laid out, the numerical implementation is made publicly available, and the approach appears readily extensible to other non‑abelian gauge theories, multi‑pion channels, or even to higher‑dimensional bootstrap problems. This paper represents a significant step forward in the program of deriving non‑perturbative QCD phenomena directly from first principles.