⚛️ phenomenology

Clarifying the puzzling mass shift of the $\psi(4160)$ via a reanalysis of $R$ -value data with unquenched charmonium spectroscopy

By reanalyzing BESII $R$ -value data within an unquenched charmonium framework that includes six vector states, this study resolves the puzzling upward mass shift of the $\psi(4160)$ by demonstrating that the observed enhancement near 4.19 GeV arises from coherent interference among nearby resonances rather than a genuine shift in the $\psi(4160)$ mass itself.

Original authors: Tian-Cai Peng, Xiang Liu

Published 2026-09-18

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Tian-Cai Peng, Xiang Liu

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the subatomic world, particles known as charmonium act as a unique laboratory for understanding the strong force, the invisible glue that binds matter together. These particles are not elementary; they are composed of a charm quark and its antimatter partner, a charm antiquark, locked in a tight embrace. For decades, physicists mapped out the energy levels of these particles much like astronomers chart the orbits of planets, expecting a neat, predictable pattern based on how the two quarks vibrate and rotate. However, as experiments moved into higher energy ranges, the pattern began to blur. The simple picture of two quarks orbiting each other started to fail, unable to explain the strange masses and behaviors observed in the data. This discrepancy signaled that the quarks were not acting in isolation; they were interacting with a sea of other particles, a complex environment that required a more sophisticated way of looking at the universe's building blocks.

For years, a specific puzzle has lingered in this high-energy region: the mass of a particle called the $\psi(4160)$ . Early measurements placed this particle at a mass of about 4.16 billion electron volts, but later, more detailed scans of the data seemed to push its mass up to roughly 4.19 billion electron volts. This shift of about 30 million electron volts was confusing because it did not fit with theoretical expectations derived from the simpler models. It suggested that either the particle was heavier than thought, or something fundamental was missing from the way scientists were interpreting the data. The question was whether the particle itself had changed, or if the method used to weigh it was flawed.

A team of researchers at Lanzhou University has now revisited this puzzle by looking at the raw data from a previous experiment with fresh eyes. Instead of relying on the older, simplified model that assumed only three main particles existed in this energy range, they applied a newer, more complete framework known as unquenched spectroscopy. In this updated view, the energy landscape is not empty but filled with a richer variety of states. Specifically, the researchers included six distinct vector charmonium states in their analysis, rather than the traditional three. They also accounted for a nearby particle called the $\psi(3770)$ , which sits at the lower edge of the energy range being studied. By treating all seven of these particles as a single, interconnected system where their signals overlap and interfere with one another, the team re-examined the historical data collected by the BESII experiment.

The results of this reanalysis offer a clear resolution to the decades-old mystery. The researchers found that the data could be perfectly explained without moving the mass of the $\psi(4160)$ from its original, lower value of 4.16 billion electron volts. The apparent bump in the data that had previously been interpreted as evidence for a heavier particle was actually the result of a complex interplay between three neighboring particles: the $\psi(4040)$ , the $\psi(4160)$ , and a newly emphasized state called the $\psi(4220)$ . When these three signals overlap, they create a combined wave that peaks at a higher energy, mimicking the appearance of a single, heavier particle. The study demonstrates that the "heavier" mass was an illusion created by the interference of nearby states, not a property of the $\psi(4160)$ itself.

This finding also sheds light on the strange, broad structures seen at even higher energies, around 4.4 billion electron volts. In the old model, this region was attributed to a single particle, the $\psi(4415)$ . However, the new analysis shows that this shape is better described by the combined contributions of three different particles: the $\psi(4380)$ , the $\psi(4415)$ , and the $\psi(4500)$ . Just as with the lower-energy bump, the complex shape of the data arises from the way these multiple signals blend together. The researchers noted that because the data can be fitted in many mathematically equivalent ways, there are multiple possible sets of values for the properties of these particles. However, by comparing these possibilities with known experimental values for the most stable particles and recent theoretical predictions, they identified a specific set of parameters that fits the physical reality best.

The study concludes that the historical shift in the measured mass of the $\psi(4160)$ was not a sign of a new physical phenomenon but a consequence of using an incomplete model to interpret the data. By embracing the complexity of the unquenched framework, which acknowledges the dynamic interaction between quarks and the surrounding particle sea, the researchers have restored the $\psi(4160)$ to its expected mass. While the current data allows for several mathematically valid solutions, the work highlights that the true picture of these particles is far more intricate than previously thought. Future experiments with higher precision will be needed to pin down the exact details of these overlapping signals, but this reanalysis provides a crucial step toward a more accurate map of the subatomic world.

Technical Summary: Clarifying the puzzling mass shift of the $\psi(4160)$ via a reanalysis of $R$ -value data with unquenched charmonium spectroscopy

Problem Statement
The paper addresses a longstanding discrepancy in charmonium spectroscopy regarding the mass of the $\psi(4160)$ resonance. Early measurements placed the mass near 4.16 GeV, whereas later inclusive $R$ -value analyses (specifically those utilizing BESII data) extracted a value around 4.19 GeV, representing an upward shift of approximately 30 MeV. This shift has remained puzzling within the context of the conventional "quenched" potential model, which describes the 4.0–4.5 GeV region using only three vector states: $\psi(4040)$ , $\psi(4160)$ , and $\psi(4415)$ . Additionally, there is a noted tension between the di-electron width extracted from inclusive $R$ -value analyses ( $\Gamma_{e^+e^-} \approx 0.48$ keV) and the di-muon width from BESIII exclusive measurements ( $\Gamma_{\mu^+\mu^-} \approx 2.45$ keV), suggesting that the extracted electromagnetic couplings may be sensitive to the assumed resonance content.

Methodology
The authors reanalyze the historical BESII inclusive $R$ -value data within an "unquenched" charmonium spectroscopy framework. Unlike the conventional quenched picture, the unquenched framework accounts for the coupling of bare quark cores to open-flavor meson-meson continua, resulting in a richer spectrum.

Resonance Content: The analysis incorporates a characteristic six-state vector-charmonium spectrum in the 4.0–4.5 GeV range: $\psi(4040)$ , $\psi(4160)$ , $\psi(4220)$ , $\psi(4380)$ , $\psi(4415)$ , and $\psi(4500)$ . These are organized into three $S$ – $D$ pairs. The $\psi(3770)$ is also included to cover the lower-energy threshold.
Fitting Framework: The total hadronic cross-section is modeled as the sum of a smooth non-resonant background (parameterized as a polynomial) and the coherent sum of seven Breit–Wigner amplitudes. The masses and total widths of the resonances are fixed based on experimental measurements (PDG, BESIII) and unquenched theoretical predictions, while the di-electron widths ( $\Gamma_{e^+e^-}$ ) and relative phases are treated as free parameters determined by the fit.
Ambiguity Analysis: The authors investigate the discrete ambiguity inherent in fitting coherent sums of multiple resonances. For seven interfering amplitudes, the analytic structure formally allows for $2^{7-1} = 64$ mathematically equivalent solutions that yield identical line shapes but possess different di-electron widths and relative phases.

Key Results

Resolution of the Mass Shift: The reanalysis demonstrates that the BESII $R$ -value line shape can be well reproduced over the full energy range (3.7–4.8 GeV) while retaining the lower-mass assignment for $\psi(4160)$ ( $M \approx 4.159$ GeV). The enhancement previously attributed to a shifted $\psi(4160)$ mass near 4.19 GeV is instead reproduced through the coherent interference of the nearby $\psi(4040)$ , $\psi(4160)$ , and $\psi(4220)$ amplitudes.
4.4 GeV Region Structure: Similarly, the asymmetric enhancement observed around 4.4 GeV is described by the coherent contributions of $\psi(4380)$ , $\psi(4415)$ , and $\psi(4500)$ , rather than a single isolated $\psi(4415)$ resonance. This suggests that the large spread in experimentally extracted widths for $\psi(4415)$ across different processes may stem from varying interference patterns.
Discrete Ambiguity and Solution Selection: The fit yields 16 numerically distinguishable solutions (derived from the 64 formal solutions) that are indistinguishable by the $R$ $R$ -value line shape alone. These solutions exhibit distinct patterns in di-electron widths. By comparing these solutions against:
- PDG reference values for well-established states ( $\psi(3770)$ , $\psi(4040)$ ).
- Recent coupled-channel theoretical predictions for $\psi(4040)$ and $\psi(4160)$ .
- Representative theoretical expectations for higher states ( $\psi(4220)$ , $\psi(4380)$ , $\psi(4500)$ ).
  The authors identify "Solution No. 1" as the most phenomenologically consistent. This solution provides the best agreement with established experimental references while remaining broadly compatible with theoretical expectations for higher states.
Fit Quality: The seven-state coherent fit achieves a $\chi^2/\text{d.o.f.} = 0.827$ , an improvement over the value of 1.05 reported in the earlier BESII analysis which utilized only three resonances.

Significance and Claims
The paper claims that the historical upward shift of the $\psi(4160)$ mass is not a direct indication of a physical mass shift but rather a line-shape effect arising from incomplete resonance content in previous analyses. By incorporating the unquenched spectrum, the inclusive $R$ -value data are shown to be compatible with a lower-mass $\psi(4160)$ consistent with quark-model expectations.

The authors emphasize that while the $R$ -value data alone cannot uniquely determine the di-electron widths of individual states due to the discrete ambiguity, the inclusion of the unquenched spectrum provides a natural framework for the observed structures. The study highlights the necessity of complementary measurements, specifically more precise inclusive $R$ -value data and di-muon channel measurements ( $e^+e^- \to \mu^+\mu^-$ ) at BESIII, to further constrain the resonance parameters and resolve the remaining discrete ambiguities among the equivalent solutions. The work underscores the transition of hadron spectroscopy into a high-precision era where unquenched effects are essential for interpreting data in the open-charm threshold region.

More like this