Observation of the Electromagnetic Dalitz Transition h_c rightarrow e^+e^-η_c
Pith reviewed 2026-05-23 23:54 UTC · model grok-4.3
The pith
The electromagnetic Dalitz decay h_c to e+e- eta_c is observed for the first time at 5.4 sigma significance.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using (27.12 ± 0.14) × 10^8 ψ(3686) decays and data samples of e⁺e⁻ collisions with √s from 4.130 to 4.780 GeV, the first observation of the electromagnetic Dalitz transition h_c → e⁺e⁻ η_c is reported with a statistical significance of 5.4σ. The ratio of branching fractions B(h_c → e⁺e⁻ η_c)/B(h_c → γ η_c) is measured separately in two production channels and averaged to (0.59 ± 0.10(stat.) ± 0.04(syst.))%.
What carries the argument
The electromagnetic Dalitz transition h_c → e⁺e⁻ η_c, which proceeds through a virtual photon and is compared directly to the radiative decay h_c → γ η_c to extract the branching-fraction ratio.
If this is right
- The measured ratio of 0.59 percent supplies a direct experimental benchmark for calculations of electromagnetic Dalitz transitions in charmonium.
- The observation confirms that the h_c state participates in this three-body electromagnetic decay mode in addition to its known radiative decay.
- The two production channels (via psi(3686) and via direct e+e- production) yield consistent ratio values, supporting the robustness of the extraction.
- The result opens the possibility of differential studies of the e+e- eta_c invariant-mass spectrum with larger data sets.
Where Pith is reading between the lines
- The ratio value can be compared with quark-model or effective-field-theory predictions to constrain the h_c wave function at the origin.
- Similar Dalitz transitions may exist for other vector or pseudoscalar charmonium states and could be searched for with the same technique.
- If the ratio remains stable with increased statistics, it may serve as a normalization channel for other rare h_c decays.
- The observation constrains the size of higher-order electromagnetic corrections in charmonium decays.
Load-bearing premise
Background contributions are correctly modeled and subtracted when the signal yield is extracted from the invariant-mass distributions.
What would settle it
An independent reanalysis that employs a different background model and finds the fitted signal yield consistent with zero at less than 5 sigma significance would falsify the observation.
Figures
read the original abstract
Using $(27.12\pm 0.14)\times10^8$ $\psi(3686)$ decays and data samples of $e^+e^-$ collisions with $\sqrt{s}$ from 4.130 to 4.780~GeV collected with the BESIII detector, we report the first observation of the electromagnetic Dalitz transition $h_c\to e^+e^-\eta_c$ with a statistical significance of $5.4\sigma$. We measure the ratio of the branching fractions $\frac{\mathcal{B}(h_c\rightarrow e^+e^-\eta_c)}{\mathcal{B}(h_c\rightarrow \gamma \eta_c)}$ separately for the $h_c$ samples produced via $\psi(3686)\to\pi^0h_c$ and $e^+e^-\to\pi^+\pi^-h_c$. The average ratio is determined to be $(0.59\pm0.10(\text{stat.})\pm0.04(\text{syst.}))\%$, where the uncertainty includes both statistical and systematic components.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports the first observation of the electromagnetic Dalitz transition h_c → e⁺e⁻ η_c using (27.12 ± 0.14) × 10^8 ψ(3686) decays and e⁺e⁻ collision data at √s = 4.130–4.780 GeV collected with BESIII. It extracts a 5.4σ statistical significance from fits to the e⁺e⁻η_c invariant mass and measures the ratio B(h_c → e⁺e⁻ η_c)/B(h_c → γ η_c) separately in the two production modes ψ(3686) → π⁰ h_c and e⁺e⁻ → π⁺π⁻ h_c, reporting an average value of (0.59 ± 0.10(stat.) ± 0.04(syst.))%.
Significance. If the background modeling and efficiency corrections hold, this constitutes the first observation of the Dalitz transition, providing a direct experimental test of electromagnetic matrix elements in the charmonium system. The use of two independent production modes and the separation of statistical and systematic uncertainties on the ratio are strengths that allow cleaner comparison with QED-based predictions for the Dalitz decay width.
major comments (2)
- [Results section describing the invariant-mass fits and background parametrization] The 5.4σ significance and the extracted signal yields that enter the ratio are determined entirely from fits to the e⁺e⁻η_c invariant mass distributions. The manuscript must demonstrate that the combinatorial and peaking background shapes (from continuum, ψ(3686) feed-down, and other sources) are independently validated for each production mode; any under-subtraction would directly inflate the significance and the reported ratio.
- [Section on branching-fraction ratio extraction and averaging procedure] The averaging of the ratio over the two production modes requires a clear statement of how the efficiency corrections, acceptance, and systematic uncertainties are combined; the current description leaves open whether the two samples are treated as fully independent or whether correlated systematics are properly accounted for in the final average.
minor comments (2)
- [Abstract and significance paragraph] The abstract states the significance is statistical only; the text should explicitly confirm whether any systematic component was considered in the significance calculation or whether the 5.4σ is purely from the likelihood ratio test on the signal yield.
- [Throughout the results and discussion sections] Notation for the two production modes should be standardized throughout (e.g., consistent use of subscripts or superscripts) to avoid ambiguity when comparing the separate ratios.
Simulated Author's Rebuttal
We thank the referee for the positive assessment and for identifying areas where additional clarity on background validation and uncertainty treatment would strengthen the manuscript. We address both major comments below and have revised the text and supplementary material accordingly.
read point-by-point responses
-
Referee: [Results section describing the invariant-mass fits and background parametrization] The 5.4σ significance and the extracted signal yields that enter the ratio are determined entirely from fits to the e⁺e⁻η_c invariant mass distributions. The manuscript must demonstrate that the combinatorial and peaking background shapes (from continuum, ψ(3686) feed-down, and other sources) are independently validated for each production mode; any under-subtraction would directly inflate the significance and the reported ratio.
Authors: We agree that explicit validation of the background parametrization for each production mode is required. In the revised manuscript we have added a dedicated subsection (now Section IV.C) that presents: (i) sideband studies in the e⁺e⁻η_c mass spectrum for both the ψ(3686)→π⁰h_c and e⁺e⁻→π⁺π⁻h_c samples, (ii) comparison of the combinatorial background shape extracted from data sidebands versus the shape obtained from inclusive Monte Carlo samples normalized to the same luminosity, and (iii) a cross-check using an alternative background parametrization (second-order polynomial versus exponential) that yields consistent signal yields within statistical uncertainties. These additions confirm that the background model does not under-subtract in either channel. revision: yes
-
Referee: [Section on branching-fraction ratio extraction and averaging procedure] The averaging of the ratio over the two production modes requires a clear statement of how the efficiency corrections, acceptance, and systematic uncertainties are combined; the current description leaves open whether the two samples are treated as fully independent or whether correlated systematics are properly accounted for in the final average.
Authors: We have expanded the description of the averaging procedure (now in Section V.B and the associated appendix). The two data samples are treated as statistically independent. Efficiency corrections and acceptance factors are evaluated separately for each mode using dedicated Monte Carlo samples. Systematic uncertainties are classified as fully correlated (e.g., tracking efficiency, photon reconstruction, and the η_c branching fractions) or uncorrelated (e.g., fit-range variations and background shape choices). The weighted average is computed with the full covariance matrix that incorporates the correlated components; the resulting total uncertainty on the averaged ratio is reported with the correlated and uncorrelated contributions shown separately. The revised text now includes the explicit formula and the breakdown of the covariance terms. revision: yes
Circularity Check
No circularity: direct experimental ratio from data yields
full rationale
The paper is a pure experimental measurement reporting the first observation of hc→e+e−ηc at 5.4σ and the ratio B(hc→e+e−ηc)/B(hc→γηc)=(0.59±0.10(stat.)±0.04(syst.))% extracted from event yields in e+e−ηc invariant-mass distributions. The derivation chain consists of data collection, background subtraction via standard parametrizations, and yield extraction; none of these steps reduce by the paper's own equations or self-citations to the reported ratio or significance. No self-definitional, fitted-input-called-prediction, or load-bearing self-citation patterns are present. The result is self-contained against external benchmarks (observed counts in two independent production modes).
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Standard model of particle physics and quantum electrodynamics govern electromagnetic transitions between charmonium states.
- domain assumption BESIII detector efficiency and resolution are accurately modeled by Monte Carlo simulation for event reconstruction.
Reference graph
Works this paper leans on
- [1]
-
[2]
R. L. Workman et al. (Particle Data Group), Prog. Theor. Exp. Phys. 2022, 083C01 (2022)
work page 2022
- [3]
- [4]
-
[5]
M. Ablikim et al. (BESIII Collaboration), Phys. Rev. Lett. 111, 242001 (2013)
work page 2013
- [6]
-
[7]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 39, 093001 (2015)
work page 2015
-
[8]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 45, 103001 (2021)
work page 2021
-
[9]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 46, 113002 (2022)
work page 2022
-
[10]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 46, 113003 (2022)
work page 2022
-
[11]
M. Ablikim et al. (BESIII Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 614, 345 (2010)
work page 2010
-
[12]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 44, 040001 (2020)
work page 2020
-
[13]
S. Agostinelli et al. (GEANT4 Collaboration), Nucl. In- strum. Methods Phys. Res., Sect. A 506, 250 (2003)
work page 2003
-
[14]
D. J. Lange, Nucl. Instrum. Methods Phys. Res., Sect. A 462, 152 (2001)
work page 2001
-
[15]
R. G. Ping, Chin. Phys. C 32, 599 (2008)
work page 2008
-
[16]
R. L. Yang, R. G. Ping, and H. Chen, Chin. Phys. Lett. 31, 061301 (2014)
work page 2014
-
[17]
Z. R. Xu and K. L. He, Chin. Phys. C 36, 742 (2012)
work page 2012
-
[18]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 106, 072007 (2022)
work page 2022
-
[19]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 96, 111101(R) (2017)
work page 2017
-
[20]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 83, 112005 (2011)
work page 2011
-
[21]
M. Ablikim et al . (BESIII Collaboration), Phys. Rev. D 89, 092008 (2014)
work page 2014
-
[22]
M. Ablikim et al . (BESIII Collaboration), Eur. Phys. J. C 76, 369 (2016)
work page 2016
-
[23]
M. Ablikim et al. (BESIII Collaboration), Chin. Phys. C 40, 113001 (2016)
work page 2016
- [24]
-
[25]
J. L. Fu, H. B. Li, X. S. Qin and M. Z. Yang, Mod. Phys. Lett. A 27, 1250223 (2012)
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.