Di-Pion Correlations in Heavy Quarkonium

• The paper demonstrates that di-pion correlations exhibit a bump–dip structure near threshold, revealing interference and final-state interaction effects in quarkonium decays.
• It employs chiral unitary and coupled-channel methods to quantitatively model ππ interactions and reproduce spectral anomalies with high statistical agreement.
• The study highlights how multi-channel interference and non-perturbative QCD dynamics provide a unified framework for interpreting charmonium and bottomonium decay patterns.

Di-pion correlation in heavy quarkonium decays refers to the strong dynamical interplay between the two-pion subsystem emerging from transitions of charmonium ($\psi$) and bottomonium ($\Upsilon$) mesons. The correlations arise from final-state interaction (FSI) effects, channel couplings, and interference among multiple decay topologies, leading to pronounced anomalies and substructure in di-pion invariant mass spectra. These phenomena provide high-sensitivity probes of non-perturbative QCD dynamics, chiral symmetry breaking, and resonance formation in the meson-meson sector.

1. Experimental Features: Observation of Di-pion Substructure

Recent data from the BESIII Collaboration, representing the largest sample of $\psi(3686)$ decays ($2.7124\pm0.0014\times10^9$ events), revealed a distinct “bump–dip” feature near the di-pion mass threshold in $\psi(3686)\to J/\psi\,\pi^+\pi^-$ (Wang et al., 13 Nov 2025). The anomaly is characterized by:

• A sharp peak at $M_{\pi\pi}\simeq0.285$ GeV.
• A pronounced dip near $M_{\pi\pi}\simeq0.305$ GeV (threshold region $0.28-0.32$ GeV).
• Statistical significance: $\chi^2/\mathrm{d.o.f.}=1.02$ in the local fit region, strongly excluding purely phase-space or QCD multipole expansion explanations.

Similar near-threshold enhancements and intermediate-mass dips are consistently observed in bottomonium transitions, e.g., $\Upsilon(nS)\to\Upsilon(mS)\,\pi^+\pi^-$, and reflected in multi-channel Dalitz plot analyses (Baru et al., 2020).

2. Chiral Unitary and Coupled-Channel Formalism

Di-pion correlations are fundamentally rooted in the resummation of strong $\Upsilon$0 final-state interactions, encoded using a chiral unitary approach (Wang et al., 13 Nov 2025, Surovtsev et al., 2016):

• The effective chiral Lagrangian for pseudoscalar-meson interactions is

$\Upsilon$1

with $\Upsilon$2 and $\Upsilon$3 GeV.

• Projected onto $\Upsilon$4-wave, the coupled channel interaction kernel $\Upsilon$5 spans $\Upsilon$6, $\Upsilon$7, and $\Upsilon$8.
• Unitarisation is performed via the Bethe–Salpeter equation,

$\Upsilon$9

with $\psi(3686)$0 regularized using either a cutoff ($\psi(3686)$1 GeV) or dimensional subtraction for threshold matching.

• The approach successfully captures key resonances ($\psi(3686)$2 and $\psi(3686)$3).

In alternative formulations, analyticity and unitarity are imposed in the Surovtsev–Bydžovský et al. coupled-channel S-matrix formalism, combining $\psi(3686)$4, $\psi(3686)$5, and $\psi(3686)$6 sectors on an eight-sheeted Riemann surface and expressing the multi-channel $\psi(3686)$7 amplitudes (Surovtsev et al., 2016).

In bottomonium transitions, a dispersive Omnès formalism is utilized (Baru et al., 2020), with detailed inclusion of heavy-quark spin-symmetry (HQSS) constraints and point-like transitions between open- and hidden-bottom channels. Short-range couplings ($\psi(3686)$8) and low-energy constants ($\psi(3686)$9) are fitted to channel-specific data, while the coupled-channel $2.7124\pm0.0014\times10^9$0-$2.7124\pm0.0014\times10^9$1 $2.7124\pm0.0014\times10^9$2-matrix is fully parameterized from Roy-equation results.

3. Correlation Function and Observables

The di-pion correlation function $2.7124\pm0.0014\times10^9$3 is a central observable, encoding both FSI dynamics and spatial emission characteristics:

$2.7124\pm0.0014\times10^9$4

• $2.7124\pm0.0014\times10^9$5: relative momentum ($2.7124\pm0.0014\times10^9$6)
• $2.7124\pm0.0014\times10^9$7: two-body phase-space factor
• $2.7124\pm0.0014\times10^9$8: spatial source function (Gaussian model),

$2.7124\pm0.0014\times10^9$9

with $\psi(3686)\to J/\psi\,\pi^+\pi^-$0, $\psi(3686)\to J/\psi\,\pi^+\pi^-$1, $\psi(3686)\to J/\psi\,\pi^+\pi^-$2 source radius.

In practical fits, $\psi(3686)\to J/\psi\,\pi^+\pi^-$3 peaks at $\psi(3686)\to J/\psi\,\pi^+\pi^-$4 near $\psi(3686)\to J/\psi\,\pi^+\pi^-$5 GeV for $\psi(3686)\to J/\psi\,\pi^+\pi^-$6 fm, and damps toward unity at higher $\psi(3686)\to J/\psi\,\pi^+\pi^-$7. Variations in $\psi(3686)\to J/\psi\,\pi^+\pi^-$8 modulate the correlation, consistent with femtoscopic and source-size analyses (Wang et al., 13 Nov 2025). The coupled-channel formalism enables generalizations to all intermediate states.

Alternative definitions employ the normalized di-pion mass spectrum as a proxy for the correlation function:

$\psi(3686)\to J/\psi\,\pi^+\pi^-$9

reflecting FSI dynamics via $M_{\pi\pi}\simeq0.285$0.

4. Interference Effects and Channel Dynamics

The observed di-pion spectral features are not resonant phenomena but arise from interference among OZI-suppressed decay topologies and coupled-channel FSI (Wang et al., 13 Nov 2025, Surovtsev et al., 2016):

• Decay amplitudes are linear combinations of $M_{\pi\pi}\simeq0.285$1 ($M_{\pi\pi}\simeq0.285$2: $M_{\pi\pi}\simeq0.285$3, $M_{\pi\pi}\simeq0.285$4, $M_{\pi\pi}\simeq0.285$5), e.g.,

$M_{\pi\pi}\simeq0.285$6

where coefficients are polynomials in $M_{\pi\pi}\simeq0.285$7 and encode tree-level and loop-induced contributions.

• The expansion,

$M_{\pi\pi}\simeq0.285$8

demonstrates constructive interference near threshold and destructive interference at intermediate masses, producing bell-shaped peaks and dips in the invariant mass spectrum.

• Phenomenological fits indicate $M_{\pi\pi}\simeq0.285$9, $M_{\pi\pi}\simeq0.305$0, phase $M_{\pi\pi}\simeq0.305$1 radians for $M_{\pi\pi}\simeq0.305$2 decays, with robust statistical agreement to the BESIII data.

These mechanisms underlie spectrum shaping in both charmonium and bottomonium transitions, with channel-dependent features (e.g., dips at $M_{\pi\pi}\simeq0.305$3–$M_{\pi\pi}\simeq0.305$4 GeV and near $M_{\pi\pi}\simeq0.305$5 GeV in $M_{\pi\pi}\simeq0.305$6).

5. Numerical Fits and Comparative Analysis

Global fits to heavy quarkonium decay data across multiple collaborations (ARGUS, CLEO, CUSB, Crystal Ball, Belle, BaBar, Mark II–III, DM2, BES II) yield $M_{\pi\pi}\simeq0.305$7 (Surovtsev et al., 2016). The coupled-channel models, with a small number of parameters, reproduce not only the sharp threshold enhancements but also broad dips and spectral distortions:

• The chiral-unitary model produces the full inclusive $M_{\pi\pi}\simeq0.305$8 spectrum up to $M_{\pi\pi}\simeq0.305$9 GeV, consistent with BESIII anomaly extension (Wang et al., 13 Nov 2025).
• Dalitz-plot fits in $0.28-0.32$0 exploit short-range B-meson interactions, HQSS-mandated relative couplings, and dispersive treatments of FSI (Baru et al., 2020).
• For $0.28-0.32$1 final states, the fitted low-energy constants are channel-dependent but highly correlated, reflecting underlying dynamical constraints.
• Switching off resonant channels confirms that spectrum features result from multi-channel interference and not from isolated resonant structures.

6. Physical Interpretation, Outlook, and Broader Implications

The observed di-pion correlations and related anomalies affirm the central role of strong FSI, chiral dynamics, and unitarity in shaping heavy quarkonium decay spectra:

• The bump–dip structure near threshold is a nonresonant interference effect, not attributable to exotic states, but to the interplay of OZI-suppressed topologies feeding $0.28-0.32$2 FSI (Wang et al., 13 Nov 2025).
• The chiral-unitary and coupled-channel approaches provide a unified description across both $0.28-0.32$3 and $0.28-0.32$4 families, consistent with universality of di-pion correlations (Surovtsev et al., 2016).
• Future femtoscopic measurements, analyzing $0.28-0.32$5 dependence on emission source size $0.28-0.32$6, collision system, or decay channel, can further constrain non-perturbative sector dynamics. Upcoming data from BESIII, Belle II, and PANDA are poised to test predictions related to threshold cusps, source radii, and channel mixing effects.
• The strong agreement between dispersive, HQSS-constrained models and experiment in $0.28-0.32$7 decays provides substantive support for the molecular interpretation of $0.28-0.32$8 and $0.28-0.32$9 as B-meson bound/resonant states (Baru et al., 2020).
• The correlation formalism opens a route for quantitative studies of femtoscopy in heavy-quarkonium environments, deepening the understanding of correlated mesonic matter and non-perturbative QCD.

In summary, di-pion correlations in heavy quarkonium decays stand as a sensitive, theoretically rigorous probe of QCD at low energies, with practical implications for particle spectroscopy, source imaging, and the dynamical structure of meson-meson interactions across multiple quarkonium systems.