---
title: Hadronic Vacuum Polarization (HVP) in QCD
url: https://www.emergentmind.com/topics/hadronic-vacuum-polarization-hvp
type: topic
---

# Hadronic Vacuum Polarization (HVP) in QCD

Hadronic vacuum polarization (HVP) is a central nonperturbative QCD effect that modifies the photon propagator via a virtual quark–antiquark polarization cloud. Its leading contribution arises from the insertion of the hadronic vacuum polarization tensor into the photon two-point function, and it dominates the theoretical uncertainty in precision Standard Model (SM) predictions for the anomalous magnetic moments of the muon and electron. The two-pion channel, governed by low-energy ππ dynamics and resonant ρ-meson exchange, provides the largest contribution and encodes the dominant isospin-1 light-quark dependence. Control of the chiral (quark-mass) extrapolation of the ππ contribution to $a_\mu^{\mathrm{HVP}}$ is indispensable for reaching sub-percent precision in lattice-QCD and phenomenological evaluations.

## 1. Omnès Representation for the Pion Vector Form Factor

The core hadronic input to the two-pion part of HVP is the pion vector form factor,
\[
F_\pi^V(s) = P(s)\exp\Bigg(\frac{s}{\pi}\int_{4m_\pi^2}^\infty \frac{\delta_1^1(s')\,ds'}{s'(s'-s-i\varepsilon)}\Bigg) \equiv P(s)\,\Omega_1^1(s)
\]
where  $\Omega_1^1(s)$ is the Omnès function built from the $\pi\pi$ $P$-wave phase shift $\delta_1^1(s)$, and $P(s)$ is a real low-degree polynomial (typically $P(s) = 1 + \beta s$) absorbing subtraction constants and residual inelastic effects. The only free parameter of $P(s)$ in this approach is the pion charge radius, $\langle r_\pi^2\rangle$, or equivalently $\beta$. This analytic representation enforces unitarity, analyticity, and correct threshold behavior.

## 2. Inverse-Amplitude Method for $\pi\pi$ Phase Shifts

The elastic $P$-wave phase shift $\delta_1^1(s)$ encodes both low-energy pion dynamics and the ρ resonance and is parametrized nonperturbatively via the inverse-amplitude method (IAM). For a partial wave $t(s)$, elastic unitarity implies
\[
\operatorname{Im} t(s) = \sigma_\pi(s)\,|t(s)|^2, \qquad \sigma_\pi(s) = \sqrt{1-4m_\pi^2/s}
\]
The IAM produces a nonperturbative approximation using chiral perturbation theory amplitudes $t_2(s)$, $t_4(s)$, $t_6(s)$ at orders $p^2$, $p^4$, $p^6$:
\[
t^{\mathrm{NLO}}_{\mathrm{IAM}}(s) = \frac{t_2^2(s)}{t_2(s) - t_4(s)} \,, \qquad
t^{\mathrm{NNLO}}_{\mathrm{IAM}}(s) = \frac{t_2^2(s)}{t_2(s) - t_4(s) + t_4^2(s)/t_2(s) - t_6(s)}
\]
and the phase shift is extracted as $\delta_1^1(s) = \arg t_1^1(s)$. The NLO (one-loop) IAM involves low-energy constants (LECs) $l_2^r-2l_1^r$, while the NNLO (two-loop) extension incorporates higher-order LECs $r_{a,b,c}$ and $r_F^r$.

## 3. Calculation of the Two-Pion HVP Contribution

With $F_\pi^V(s)$ constructed as above, the two-pion contribution to $a_\mu^{\mathrm{HVP}}$ is given by the standard dispersion relation:
\[
a_\mu^{\mathrm{HVP},\,\pi\pi} = \left(\frac{\alpha}{\pi}\right)^2 \int_{4m_\pi^2}^{\infty} ds\; \frac{K(s)}{12\pi}\, \beta_\pi^3(s)\, |F_\pi^V(s)|^2
\]
where $K(s)$ is a calculable kernel maximizing the weight near $s \sim m_\mu^2$ and $\beta_\pi(s) = \sqrt{1-4m_\pi^2/s}$. This representation ensures that the dominant long-distance, low-energy hadronic physics is correctly incorporated, and the dispersively determined function agrees with all relevant low-energy constraints.

## 4. Chiral Extrapolation in Lattice QCD and Phenomenology

To control the dependence on unphysical pion masses (as used in lattice ensembles above or below the physical point), the following procedure is adopted:
- The Omnès input is fixed at the physical point by the empirical phase shifts and measured $\langle r_\pi^2\rangle=0.429(4)\,{\rm fm}^2$.
- The IAM parametrization for $\delta_1^1(s)$ is expressed in terms of the chiral parameters ($F$, $m_\pi$, $l_i^r$, $r_i^r$), with the two-loop pion charge radius formula,
  \[
  \langle r_\pi^2 \rangle = \frac{1}{16\pi^2F^2}\left\{ R_4 + \frac{m_\pi^2}{16\pi^2F^2}R_6 \right\}
  \]
  $R_4$, $R_6$ are functions of the relevant LECs.
- The $m_\pi$ dependence of the IAM input is propagated into $F_\pi^V(s)$, the Omnès function, and $a_\mu^{\mathrm{HVP},\,\pi\pi}(m_\pi)$ constructed.
- Inelastic (non-elastic) effects are included as a mild polynomial correction to $P(s)$.

Numerical implementation yields:
\[
\begin{align*}
\bar a_\mu^{\mathrm{HVP}}[\pi\pi,\leq 1\,\mathrm{GeV}]|_{\rm NNLO,IAM} &= 508.0(28.8)_{\rm fit}(1.6)_\chi(8.5)_{r_\pi} \times 10^{-10} \\
\bar a_\mu^{\mathrm{HVP}}[\pi\pi,\leq 1\,\mathrm{GeV}]|_{\rm NNLO,IAM+\delta} &= 482.4(0.1)(0.7)(8.0)\times 10^{-10}
\end{align*}
with the errors from fit parameters, chiral truncation, and $\langle r_\pi^2\rangle$ uncertainties. Imposing the precise physical $\delta_1^1(s)$ drastically reduces uncertainty. Varying the cutoff to infinity modifies the result slightly, with chiral-truncation errors remaining below $5\times10^{-10}$ for all $m_\pi \leq 250$ MeV.

## 5. Systematic Uncertainties and Chiral Fit Strategies

The IAM-driven chiral extrapolation enables robust quantification of the $m_\pi$ dependence. Systematic uncertainties from inelastic corrections and higher-order LECs are assessed to contribute at most a few $10^{-10}$, which is orders of magnitude below present statistical and systematic errors from lattice QCD and R-ratio phenomenology. Practical fitting strategies for lattice calculations include:
- Fitting the chiral LECs $\{F, l_i^r, r_i^r\}$ directly to multi-ensemble lattice data for the phase shift and decay constant, with $\beta$ (or $\langle r_\pi^2\rangle$) as the only free normalization parameter for $F_\pi^V(s)$; inelastic effects are included as simple polynomials in $m_\pi^2$.
- Alternatively, working in the spacelike domain and fitting the subtracted vacuum polarization function as
  \[
  \frac{\bar\Pi(-Q^2)}{Q^2} = \frac{a + b Q^2}{1 + c Q^2 + d Q^4}
  \]
  with the $m_\pi$ dependence determined by the IAM, typically requiring only a $m_\pi^{-2}$ and possibly a $\log m_\pi^2$ term.

## 6. Physical Implications, Precision, and Lattice Impact

The dominance of the $\pi\pi$ channel in the isospin-1 correlator ensures that IAM-based chiral referencing directly constrains the leading light-quark HVP contribution with controlled and negligible uncertainties relative to the experimental or phenomenological targets. For lattice QCD determinations aiming at sub-percent uncertainties in $a_\mu^{\mathrm{HVP}}$, the IAM-guided chiral extrapolation provides an essential and reliable bridge from heavier-than-physical pion masses to the physical point. Uncertainties from inelastic channels and higher-order effects remain well below the precision goal for the dominant component.

A plausible implication is that IAM-based parametrizations will continue to serve as the standard reference for chiral and mass extrapolations in high-precision lattice studies of the HVP, especially as statistical errors decrease and the controlling of systematic uncertainties from pion-mass dependence becomes ratio-limiting. The methodology uniquely ensures that low-energy QCD constraints, resonance dynamics, and the correct analytic structure are implemented ab initio in the extrapolation needed for HVP precision studies.

Source: https://www.emergentmind.com/topics/hadronic-vacuum-polarization-hvp