⚛️ phenomenology

Mass spectrum and decay widths of charmonium-like mesons: A diabatic approach with complex scaling

This paper extends a diabatic framework for charmonium-like mesons by incorporating the complex scaling method to treat bound and resonant states on equal footing, thereby enabling parameter-free extraction of decay widths and revealing significant molecular characteristics in states such as $\chi_{c1}(3872)$ , $\psi(4040)$ , and $\psi(4230)$ while proposing $X(3940)$ as a $1^{++}$ candidate.

Original authors: Zi-Zhao Zhang, Bo-Chao Liu

Published 2026-09-25

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Zi-Zhao Zhang, Bo-Chao Liu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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

Deep within the heart of every atom lies a world of particles that behave less like solid balls and more like a complex, shifting fog. Among these, a family of particles called charmonium has long served as a testing ground for our understanding of the strong force, the invisible glue that binds the universe's smallest building blocks together. For decades, physicists have tried to map out the "spectrum" of these particles, much like a musician tuning a piano to find every possible note. The traditional view held that these particles were simple pairs of a charm quark and its antimatter partner, locked together in a neat, predictable orbit. However, in recent years, experiments have revealed a crowded and confusing landscape of new, strange particles that do not fit this simple picture. These new states often appear right at the edge of where they can break apart into other particles, suggesting that their identity is not just a single pair, but a messy, dynamic mix of different configurations. Understanding exactly what these particles are and how they decay is crucial, because it reveals the deeper rules of how matter holds itself together.

A team of researchers at Xi'an Jiaotong University has taken a fresh look at this problem by building a more complete model of how these particles interact. Instead of treating the particles as isolated pairs, their new approach explicitly includes the constant possibility that the particle can briefly transform into a pair of lighter particles before snapping back together. This continuous exchange, known as coupling to a continuum, is difficult to calculate because it involves states that are not stable and do not last long enough to be measured directly. To solve this, the researchers employed a mathematical technique called complex scaling. Imagine trying to listen to a faint, short-lived sound in a noisy room; this method effectively turns down the noise of the background while amplifying the specific signal of the fleeting particle, allowing the scientists to calculate its properties with high precision without needing to invent new, unknown variables.

Using this unified framework, the team calculated the masses and decay widths—the rate at which these particles fall apart—for a wide range of charmonium-like states below 4.3 GeV. Their results offer a clearer picture of the internal structure of these mysterious particles. They found that the famous X(3872) particle, long a puzzle because of its mass being almost identical to the threshold for breaking apart, behaves exactly as a loosely bound molecule made of two different mesons rather than a tight, conventional pair. Similarly, they identified the ψ(4040) and ψ(4230) states as having significant molecular characteristics, meaning they are heavily influenced by the surrounding cloud of lighter particles. The study also proposed a new identity for the X(3940), suggesting it is a specific type of particle with quantum numbers that allow it to decay in ways consistent with experimental observations, acting as a hybrid mixture of a traditional charm pair and a molecular state.

The researchers also revisited the nature of the χc0(3860), a candidate for the first excited scalar charmonium state. While their calculated mass aligns perfectly with the experimental value for this particle, their model predicts a much narrower decay width than what is currently observed in experiments. This discrepancy suggests that the true nature of this particle might be more complex than the current model can capture, possibly involving hidden mechanisms not yet included in the calculation. In the sector of vector particles, the model successfully reproduced the number of excited states seen in experiments, though the calculated masses and widths showed minor deviations. The study confirms that for many of these exotic states, the simple picture of a single quark pair is insufficient; instead, the particles are deeply dressed in a cloud of other particles, blurring the line between a fundamental state and a composite molecule.

One limitation of this work is that it focused only on open-charm channels, where the particles can break apart into lighter, open-flavor mesons. Several other exotic candidates observed in experiments, such as the X(3915), decay into hidden-charm modes that were not included in this specific calculation. The absence of these states in the new spectrum suggests that their formation relies on different dynamics that require further investigation. Nevertheless, this study provides a robust, parameter-free foundation for understanding the charmonium spectrum, demonstrating that by treating bound states and fleeting resonances on equal footing, physicists can begin to unravel the complex, hybrid nature of the particles that populate the subatomic world.

Technical Summary: Mass spectrum and decay widths of charmonium-like mesons: A diabatic approach with complex scaling

Problem Statement
The experimental observation of numerous charmonium-like states (XYZ states) near open-charm thresholds challenges the conventional $q\bar{q}$ description. States such as $X(3872)$ exhibit masses and decay patterns inconsistent with simple potential models. While the Cornell-type potential successfully describes low-lying states, it fails to provide a unified description for highly excited or near-threshold states where coupled-channel effects and meson-meson continuum contributions are significant. Previous diabatic frameworks addressed the spectrum but often neglected continuum contributions to decay widths or treated bound and resonant states separately, requiring perturbative corrections or additional free parameters to estimate widths.

Methodology
This work extends a previous diabatic framework by integrating the Complex Scaling Method (CSM) to treat bound states and resonances on an equal footing.

Diabatic Formalism: The system is described using a multi-channel Schrödinger equation with a matrix-valued potential. The Hamiltonian includes a heavy quark-antiquark ( $c\bar{c}$ ) core and open-flavor meson-meson channels. The potential matrix incorporates:
- Diagonal Terms: A spin-dependent Cornell potential for the $c\bar{c}$ core (including confinement, Coulomb, spin-spin, spin-orbit, and tensor terms) and constant thresholds for meson-meson channels.
- Off-Diagonal Terms: Mixing potentials ( $V_{mix}$ ) representing string breaking, parametrized as Gaussian functions of the energy difference between the $c\bar{c}$ potential and the meson-meson threshold.
- Spin-Dependence: Explicit inclusion of spin-dependent interactions derived from one-gluon exchange, regularized via Gaussian smearing to handle short-distance singularities.
Complex Scaling Method (CSM): To extract decay widths without perturbative approximations, the coordinates and momenta are rotated into the complex plane ( $r \to re^{i\theta}$ , $p \to pe^{-i\theta}$ ). This transforms the Hamiltonian into a non-Hermitian form where resonant states appear as discrete complex eigenvalues ( $E = M - i\Gamma/2$ ) separated from the rotated continuum.
- Numerical Solution: The coupled-channel equations are solved using the Gaussian Expansion Method (GEM).
- Observables: The real part of the eigenvalue yields the mass ( $M$ ), and the imaginary part yields the decay width ( $\Gamma$ ). The study also calculates complex root-mean-square radii ( $r_{RMS}$ ) and channel component weights ( $P_i$ ) using the $c$ -product (biorthogonal product) to analyze internal structures.
Parameters: The model utilizes parameters fixed by fitting established charmonium states and lattice QCD results. The mixing strength ( $\Delta$ ) is constrained by the mass of the $\chi_{c1}(3872)$ , while the mixing range ( $\rho$ ) is fixed based on previous diabatic studies. No additional free parameters are introduced for the continuum treatment.

Key Contributions and Results
The study calculates the mass spectrum and decay widths for charmonium-like states below 4.3 GeV across various $J^{PC}$ sectors:

$1^{++}$ Sector:
- $\chi_{c1}(3872)$ : Identified as a pole near the $D\bar{D}^*$ threshold with a vanishingly small width. The analysis reveals a dominant $D\bar{D}^*$ component weight (88.4%) and a large complex RMS radius, supporting its interpretation as a loosely bound hadronic molecule.
- $X(3940)$ : A pole is found at $M \approx 3941$ MeV with $\Gamma \approx 40$ MeV. The authors propose this as a candidate for the experimentally observed $X(3940)$ with $J^{PC}=1^{++}$ . The state shows significant mixing between $c\bar{c}$ and $D\bar{D}^*$ configurations, and the $1^{++}$ assignment is consistent with the observed suppression of the $D\bar{D}$ decay mode.
$0^{++}$ Sector:
- $\chi_{c0}(3860)$ : A pole is identified at $M \approx 3865$ MeV with a narrow width ( $\Gamma \approx 4$ MeV). While the mass aligns well with the experimental $\chi_{c0}(3860)$ , the calculated width is significantly narrower than the experimental central value ( $\sim 200$ MeV). The state is predominantly $c\bar{c}$ (77%). The paper notes this tension suggests the experimental nature of $\chi_{c0}(3860)$ may be more complex.
- $X(3915)$ : No pole corresponding to the $X(3915)$ is found, consistent with the view that it cannot be described as a conventional $c\bar{c}$ state coupled primarily to open-charm channels.
$2^{++}$ Sector:
- $\chi_{c2}(3930)$ : The calculated mass ( $\approx 3900$ MeV) and width ( $\approx 36$ MeV) agree well with the experimental $\chi_{c2}(3930)$ . Component analysis indicates substantial mixing, with the $c\bar{c}$ weight at $\sim 48\%$ , highlighting the importance of meson-meson clouds.
- Prediction: A second pole is predicted at $M \approx 4009$ MeV ( $\Gamma \approx 17$ MeV), tentatively assigned to the $\chi_{c2}(1F)$ state.
$1^{--}$ (Vector) Sector:
- The model reproduces the multiplicity of observed vector states ( $\psi(2S)$ , $\psi(3770)$ , $\psi(4040)$ , $\psi(4160)$ , $\psi(4230)$ ).
- $\psi(4040)$ : Identified as having a significant $D^*\bar{D}^*$ molecular component ( $c\bar{c}$ weight $\sim 34\%$ ).
- $\psi(4230)$ : Identified with a dominant $D_s\bar{D}_s^*$ molecular component, despite the non-Hermitian nature of the weights.
Other Sectors: Bound states for $\eta_c(1S, 2S)$ and $h_c(1P)$ are reproduced with high accuracy. Predictions are made for resonances in the $0^{-+}$ ( $\sim 4085$ MeV) and $1^{+-}$ ( $\sim 3937$ MeV) sectors.

Significance and Claims
The paper claims that the integration of the Complex Scaling Method into the diabatic framework provides a unified, parameter-free description of both bound and resonant charmonium-like states. Key achievements include:

Direct Extraction of Widths: Decay widths are obtained directly from complex eigenvalues without perturbative corrections or ad-hoc assumptions.
Structural Insight: The calculation of complex component weights and RMS radii elucidates the internal structure of exotic candidates, confirming the molecular nature of $\chi_{c1}(3872)$ and the mixed nature of $\chi_{c2}(3930)$ .
Candidate Identification: The work provides strong theoretical support for the $1^{++}$ assignment of the $X(3940)$ and offers a mass-consistent identification for $\chi_{c0}(3860)$ as the $\chi_{c0}(2P)$ , while highlighting discrepancies in the width that warrant further investigation.
Limitations and Future Directions: The authors explicitly note that the current model is restricted to open-charm channels. The absence of states like $X(3915)$ , $\chi_{c1}(4140)$ , and $\chi_{c1}(4274)$ in the spectrum suggests that mechanisms involving hidden-charm channels (e.g., $J/\psi\phi$ , $\omega J/\psi$ ) are necessary to fully describe these exotic candidates. The paper concludes that incorporating these hidden-charm channels is a necessary extension for a complete physical picture.

More like this