Two-lepton tales: Dalitz decays of heavy quarkonia

We study the Dalitz decays of heavy quarkonia, which result from the internal virtual photon conversion into an $\ell^+ \ell^-$ lepton pair. Heavy-quark symmetries allow us to establish systematic relations between transitions of different quarkonium states, and to precisely determine the branching fractions for several charmonium and bottomonium decay modes. For charmonium, existing data on $χ_{cJ}(1P)\to J/ψ\ell^+ \ell^-$ and $ψ(2S)\to χ_{cJ}(1P) \ell^+ \ell^-$ enable us to determine the parameters of the transition form factors and to predict the rates of yet-unobserved modes. The Dalitz transitions of $χ_{c1}(3872)$ are important, as they can help assessing the structure of this meson. For bottomonium, recent LHCb measurements allow us to predict the branching fractions of $χ_{bJ}(nP)\to Υ(1S)\ell^+ \ell^-$ and $h_b(nP)\to η_b(1S) \ell^+ \ell^-$ ($n=1,,2)$. We also investigate the sensitivity of heavy quarkonia Dalitz modes to the contribution of a new light vector mediator, such as the putative $X(17)$.

💡 Research Summary

This paper presents a comprehensive study of Dalitz decays of heavy quarkonia—processes in which an internal virtual photon converts into a charged lepton pair (ℓ⁺ℓ⁻). The authors start by deriving a model‑independent relation between the differential Dalitz decay rate and the corresponding radiative decay width, Eq. (1):

 dB(M′→M ℓ⁺ℓ⁻)/dq² = B_rad(M′→M γ) × λ½(m²_M′, m²_M, q²) λ½(m²_M′, m²_M, 0) × F_QED(q²) × |f(q²)|²,

where λ is the usual triangular kinematic function, F_QED(q²) is the pure QED factor, and f(q²) is the transition form factor (TFF) normalized to unity at q² = 0. The TFF encodes non‑perturbative QCD dynamics and is the only unknown once the radiative branching fraction B_rad is known.

To constrain f(q²) the authors exploit heavy‑quark spin symmetry (HQSS) that emerges in the limit of infinite heavy‑quark mass. Within HQET they construct effective Lagrangians for electric‑dipole (E1) transitions between P‑wave (ℓ = 1) and S‑wave (ℓ = 0) multiplets, and for magnetic‑dipole (M1) transitions within the S‑wave doublet. The key result is that a single coupling constant δ_{nP mS} governs all E1 transitions from a given P‑wave radial level n to an S‑wave level m, irrespective of the total spin J of the initial state. This dramatically reduces the number of independent parameters.

The authors then adopt a simple vector‑meson‑dominance inspired pole ansatz for the TFF:

 f(q²) = 1 / (1 − q²/a²),

with a single pole mass a that is assumed to be the same for all 1P→1S transitions in the heavy‑quark limit. Using the measured Dalitz branching fractions for χ_{c1}(1P)→J/ψ μ⁺μ⁻ and χ_{c2}(1P)→J/ψ μ⁺μ⁻, they extract a_{c1}=0.83^{+0.17}{‑0.11} GeV and a{c2}=0.71^{+0.13}_{‑0.07} GeV, confirming the HQSS expectation that the two values are compatible. By performing a combined χ² fit they obtain a common pole mass a ≈ 0.77 ± 0.10 GeV, which they then use to predict the Dalitz branching fractions for all charmonium 1P→1S transitions, both in the electron and muon channels. Representative results include

 B(χ_{c0}→J/ψ e⁺e⁻) ≈ 4 × 10⁻⁶,
 B(χ_{c1}→J/ψ μ⁺μ⁻) ≈ 2.3 × 10⁻⁴,
 B(χ_{c2}→J/ψ μ⁺μ⁻) ≈ 2.0 × 10⁻⁴,

and analogous predictions for the yet‑unobserved χ_{c0,1,2}→ψ(2S)ℓ⁺ℓ⁻ modes.

Turning to bottomonium, the authors employ recent LHCb measurements of the radiative decays χ_{bJ}(nP)→Υ(1S) γ and h_b(nP)→η_b(1S) γ (n = 1, 2) to fix the same coupling δ_{nP 1S}. Assuming the same pole mass a, they predict Dalitz branching fractions such as

 B(χ_{b1}(1P)→Υ(1S) μ⁺μ⁻) ≈ 3.5 × 10⁻³,
 B(χ_{b2}(2P)→Υ(1S) μ⁺μ⁻) ≈ 6.6 × 10⁻³,

and corresponding electron‑channel rates. These predictions are within reach of forthcoming LHCb and Belle II analyses.

A particularly interesting application concerns the exotic state X(3872), identified here as χ_{c1}(3872). The Dalitz decays χ_{c1}(3872)→J/ψ ℓ⁺ℓ⁻ and χ_{c1}(3872)→ψ(2S) ℓ⁺ℓ⁻ are sensitive to the internal structure (compact charmonium vs. D⁰\bar{D}^{*0} molecule). By measuring the dilepton invariant‑mass spectrum and comparing it with the TFF shape inferred from the radiative widths B(χ_{c1}(3872)→J/ψ γ) and B(χ_{c1}(3872)→ψ(2S) γ), one can discriminate between competing models.

Finally, the paper investigates the potential impact of a new light vector boson, the hypothesized X(17) particle, on Dalitz spectra. If X(17) mixes with the virtual photon, it would introduce an additional resonant contribution at q² ≈ (17 MeV)², most visible in the electron channel because of the larger phase space. The authors estimate that current experimental uncertainties are too large to resolve such a tiny effect, but a future sensitivity at the 10⁻⁶ level in branching fractions would allow a meaningful search for X(17) in heavy‑quarkonium Dalitz decays.

In summary, the work provides a unified, symmetry‑based framework linking radiative and Dalitz decays of heavy quarkonia, extracts the transition form‑factor parameters from existing data, delivers a suite of concrete predictions for unmeasured charmonium and bottomonium modes, and outlines how these processes can probe both the nature of exotic states like X(3872) and possible new light gauge bosons. The methodology and predictions set the stage for a rich experimental program at LHCb, BESIII, and Belle II.