Global Dispersive Analysis of ππ Scattering

• The paper presents a data-driven global dispersive analysis that rigorously extracts ππ scattering amplitudes using N/D decomposition and conformal mapping techniques.
• It employs Roy and Roy–Steiner equations along with chiral EFT constraints to ensure a model-independent determination of low-energy thresholds and resonance poles.
• The analysis delivers precise S-wave phase shifts and resonance parameters, forming a robust framework for light scalar meson spectroscopy and QCD studies.

Global dispersive analysis of $ππ$ scattering refers to the use of analyticity, unitarity, and crossing symmetry, supplemented by data-driven and chiral EFT constraints, to reconstruct the partial-wave amplitudes for $\pi\pi\to\pi\pi$ from threshold up to the onset of Regge behavior, and to perform analytic continuation for rigorous extraction of resonance pole parameters. Modern implementations employ the $N/D$ decomposition, conformal mapping techniques for the left-hand cut, Roy and Roy-Steiner equations, and multi-channel unitarity, yielding parameterizations consistent with both direct experimental data and dispersive sum rules. This framework provides precise, model-independent determinations of low-energy phases, threshold parameters, Omnès functions, and resonance properties, and forms the backbone of contemporary analyses of hadron spectroscopy, meson-meson interactions, and constraints on QCD dynamics.

1. Analytic Structure and Partial-Wave Dispersion Relations

Within the global dispersive framework, each S-wave amplitude $t_{ab}(s)$ (with $a,b$ labeling channels such as $\pi\pi$, $K\bar{K}$) is governed by a once-subtracted dispersion relation,

$$t_{ab}(s) = t_{ab}(0) + \frac{s}{\pi} \int_{-\infty}^{s_L} \frac{ds'}{s'} \frac{\text{Disc}[t_{ab}(s')]}{s' - s} + \frac{s}{\pi} \int_{s_{\text{th}}}^{\infty} \frac{ds'}{s'} \frac{\text{Disc}[t_{ab}(s')]}{s' - s}\,,$$

where $s_{\text{th}}$ is the threshold and $s_L$ marks the onset of the left-hand cut. Unitarity dictates the discontinuity on the right-hand cut,

$$\text{Disc}\,t_{ab}(s) = \sum_{c} t_{ac}(s)\,\rho_c(s)\,t_{cb}^*(s),\qquad \rho_c(s)=\frac{\sqrt{1-4m_c^2/s}}{16\pi s} \Theta(s-4m_c^2)\,.$$

The $N/D$ representation isolates analytic structures: $t_{ab}(s) = \sum_c D^{-1}_{ac}(s)\,N_{cb}(s),$ where $N_{ab}(s)$ has only left-hand cuts and $D_{ab}(s)$ only right-hand unitarity cuts. These satisfy coupled integral equations, with $U_{ab}(s)$ driving the left-cut via

$$U_{ab}(s)=t_{ab}(0)+\frac{s}{\pi} \int_{-\infty}^{s_L} \frac{ds'}{s'} \frac{\text{Disc}[t_{ab}(s')]}{s'-s}\,.$$

This structure is foundational in the recent analyses by Danilkin et al. and underpins the semi-analytic solution of the amplitude up to the $K\bar{K}$ threshold (Deineka et al., 2022).

2. Roy and Roy–Steiner Constraints

Crossing symmetry and analyticity are enforced through the Roy equations, which provide exact relations between partial waves and their imaginary parts in $s$ and $t$ channels up to $\sim$1.1 GeV: $$t^I_\ell(s) = k^I_\ell(s) + \sum_{I'=0}^2 \sum_{\ell'=0}^\infty \int_{4m_\pi^2}^\infty \frac{dt}{t^2} K^{II'}_{\ell\ell'}(s,t) \Im\,t^{I'}_{\ell'}(t)\,.$$ The subtraction polynomials $k^I_\ell(s)$ are determined by the S-wave scattering lengths $a^0_0,a^2_0$, themselves fixed to high-precision values from NNLO chiral perturbation theory, $a^0_0\approx0.220\pm0.005$, $a^2_0\approx-0.0444\pm0.0010$. Adler zeros $s^I_A$ are fixed to NLO $\chi$PT values with conservative uncertainties, and slope parameters match NNLO theory. These constraints ensure model-independent consistency of the global fit with the full crossing relations and the low-energy expansion (Deineka et al., 2022, Kolesár et al., 31 Jul 2025).

3. Conformal Mapping and Parameterization of the Left-Hand Cut

Unknown left-hand cut contributions $U_{ab}(s)$ are represented by an expansion in a conformal variable,

$$\omega_{ab}(s) = \frac{\sqrt{s_L^{(ab)}-s_E}-\sqrt{s_L^{(ab)}-s}} {\sqrt{s_L^{(ab)}-s_E}+\sqrt{s_L^{(ab)}-s}},$$

where $s_L^{(ab)}$ is the closest branch point and $s_E$ is the expansion point. $U_{ab}(s)$ is then truncated,

$$U_{ab}(s) = \sum_{n=0}^N C_{ab,n}[\omega_{ab}(s)]^n,$$

with typically $N=3$ or 4. The coefficients $C_{ab,n}$ are treated as free parameters in fits to experimental phase shifts, inelasticities, and Roy/Roy–Steiner pseudo-data, absorbing uncertainties from short-distance physics or omitted high-order LHC dynamics (Deineka et al., 2022, Biloshytskyi et al., 2023).

4. Numerical Solution and Data-Driven Fitting Strategies

In practice, the coupled integral equations are solved via iterative procedures until convergence. In the isoscalar ($I=0$) sector, $\pi\pi$ and $K\bar{K}$ amplitudes are fitted together to experimental S-wave phase shifts $\delta^0_0(s)$ (up to 1.2 GeV), inelasticity $\eta^0_0(s)$ above $K\bar{K}$ threshold, and modulus/phase of $\pi\pi\to K\bar{K}$ scattering. Roy–Steiner solutions serve as pseudo-data for both $\pi\pi$ and $\pi\pi\to K\bar{K}$. For the $I=2$ channel, a single-channel fit suffices. All fits use a $\chi^2$ procedure with full error propagation, including bootstrap resampling and systematic studies of the conformal-map expansion point. When Roy–like pseudo-data are used, $\chi^2$/d.o.f. is typically below unity, reflecting large data correlations (Deineka et al., 2022).

5. Analytic Continuation and Extraction of Resonance Poles

Analytic continuation of the amplitude into the complex $s$-plane provides rigorous resonance pole determination. Poles correspond to zeros of $\det D(s)$ on unphysical Riemann sheets: $$\det D(s_p)=0,\qquad t_{ab}(s)\sim\frac{r_{ab}}{s_p-s},\quad(s\to s_p),$$ For $I=0$ S-wave, the pole positions determined in MeV are: $$\sqrt{s_\sigma}=458(10)^{+7}_{-15} - i\,256(9)^{+5}_{-8},\qquad \sqrt{s_{f_0}}=993(2)^{+2}_{-1} - i\,21(3)^{+2}_{-4},$$ in agreement with Roy averages and alternative dispersive analyses. For $\pi K$ ($I=1/2$), $$\sqrt{s_\kappa}=702(12)^{+4}_{-5}-i\,285(16)^{+8}_{-13}$$. The extracted residues yield the $\pi\pi$ and $K\bar{K}$ couplings, which are in mutual agreement across Roy-like frameworks (Deineka et al., 2022).

6. S-wave Phase Shifts, Threshold Parameters, and Omnès Functions

Global dispersive fits yield continuous S-wave phase shifts $\delta_0^0(s),\delta_0^2(s)$ and inelasticity $\eta_0^0(s)$, reproducing both experimental and Roy pseudo-data up to 1.2 GeV. The threshold scattering lengths extracted are: $a_0^0=0.220(5),\qquad a_0^2=-0.0444(10),$ again fully consistent with NNLO $\chi$PT and Roy solutions. Channel-by-channel Omnès matrices are defined via,

$\Omega_{ab}(s)\equiv D_{ab}^{-1}(s),$

Reducing in the single-channel case (e.g., $I=2$) to

$$\Omega(s)=\exp\left[\frac{s}{\pi} \int_{4m_\pi^2}^\infty \frac{ds'}{s'} \frac{\delta_0^2(s')}{s'-s}\right],$$

which reconstructs the relevant analytic form factors and production amplitudes (Deineka et al., 2022).

7. Model Independence, Consistency, and Application Scope

The N/D dispersive analysis, enforced by Roy and $\chi$PT constraints and informed by direct experimental data, yields a globally consistent, model-independent representation of $\pi\pi$ S-wave scattering. It matches Roy/Roy-Steiner results in both the real and complex $s$-plane, achieves excellent agreement up to $~$1.2 GeV, and delivers pole parameters for light scalar resonances compatible across high-precision dispersive studies. These results form solid input for the low-energy QCD sector, enable rigorous lattice QCD extrapolations, and facilitate the construction of hadronic form factors for weak and electromagnetic processes. The flexible framework is readily generalized to multi-channel scenarios and to alternative two-body scattering reactions (Deineka et al., 2022, Danilkin et al., 2020).

---

Feature	Treatment in Global Dispersive Analysis	Reference
S-wave dispersion	Once-subtracted, coupled $N/D$, conformal LHC	(Deineka et al., 2022)
Crossing/threshold	Roy and Roy-Steiner equations, NNLO $\chi$PT	(Kolesár et al., 31 Jul 2025)
Analytic continuation	Zeros of $\det D(s)$ for pole extraction	(Deineka et al., 2022)
Model independence	Direct fit, LHC mapping, pseudo-data consistency	(Deineka et al., 2022, Pelaez et al., 2019)

The consistency and precision achieved in these analyses position global dispersive approaches as the standard for rigorous, quantitative studies of hadronic $\pi\pi$ dynamics and light scalar meson spectroscopy.