pith. sign in

arxiv: 2509.12925 · v1 · submitted 2025-09-16 · ✦ hep-ph · hep-ex· nucl-th

Coupled-channel analysis for vector charmonia and their nature

Satoshi X. Nakamura (Shandong University) This is my paper

Pith reviewed 2026-05-18 16:44 UTC · model grok-4.3

classification ✦ hep-ph hep-exnucl-th
keywords charmoniumvector charmoniacoupled-channel modelmolecular compositenessD* D* moleculeopen-charm thresholdspole extractionBESIII data
0
0 comments X

The pith

Coupled-channel analysis of collider data shows ψ(4040) is primarily a D* anti-D* molecule rather than a quark-model state.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper fits a unitary model to high-precision electron-positron data on twenty charmonium final states between 3.75 and 4.7 GeV. From the resulting amplitudes the authors extract pole positions and compute their compositeness, which quantifies how much each open-charm channel contributes to the state. The calculation indicates that the ψ(4040) pole is dominated by a D* anti-D* molecular component instead of the usual ψ(3S) quark-antiquark assignment. The same method finds that ψ(4230) and ψ(4360) contain sizable mixtures of D1(2420) D, D1(2420) D*, Ds* Ds*, and conventional c cbar pieces. This matters because it supplies a concrete, data-driven way to decide whether certain charmonia are bound meson pairs or conventional excitations, directly affecting decay-rate predictions.

Core claim

The analysis of e+e- to c cbar data with a semi three-body unitary coupled-channel model yields vector charmonium poles whose compositeness shows that ψ(4040) could mainly consist of a D* D*-molecule component rather than a conventionally accepted quark-model ψ(3S) state, while ψ(4230) and ψ(4360) might be substantial mixtures of D1(2420) D-bar, D1(2420) D*, Ds* Ds*, and c c-bar components.

What carries the argument

The semi three-body unitary coupled-channel model that parametrizes amplitudes for twenty final states, incorporates three-body cuts and form factors, and extracts pole positions together with their channel compositeness from the global fit.

Load-bearing premise

The chosen twenty final states and the specific parametrization of three-body cuts and form factors are assumed to be sufficient to determine the pole positions and compositeness without large bias from omitted channels.

What would settle it

A measurement that the ψ(4040) branching fractions or line shapes deviate strongly from those expected for a dominant D* anti-D* molecular state, for example a much smaller rate into D* D* pairs than the model predicts.

Figures

Figures reproduced from arXiv: 2509.12925 by Satoshi X. Nakamura (Shandong University).

Figure 1
Figure 1. Figure 1: (a) e +e − → abc mechanism involving charmonium excitation. (Bare) two-meson resonances R and bare charmonium states are represented by the solid and double lines, respectively. Dressed propagators and vertices are indicated by the solid circles. (b) Direct and single triangle decay mechanisms of charmonium. coupled-channel model that respects three-body unitarity. Such a coupled-channel analysis is the pu… view at source ↗
Figure 2
Figure 2. Figure 2: e +e − annihilation cross sections (unit:pb); each panel indicates the final state; √ s is the total energy. Full calculations are shown by the red points,connected by lines. The direct decays, single-triangle, and nonresonant contributions are shown by the blue dashed, magenta dotted, and green dash-dotted curves, respectively. Figures taken from Ref. [1] where references for the data are given. 3 Fit res… view at source ↗
Figure 3
Figure 3. Figure 3: Vector charmonium poles (Eψ) and their uncertainties. Red, blue, and green points indicate pole locations of resonances (located on unphysical sheets of open channels), bound, and virtual states, respectively; the bound (virtual) states are on the physical (unphysical) sheets of the nearest-threshold channels, respectively. Black points indicate ψ states listed in PDG [5], R(3760) [6], G(3900) [7], and Y (… view at source ↗
read the original abstract

High-precision $e^+e^-\to c\bar{c}$ data (20 final states) from the BESIII and Belle in $\sqrt{s}=3.75-4.7$ GeV are analyzed with a semi three-body unitary coupled-channel model. Vector charmonium poles are extracted from the amplitudes obtained from the fit. We find well-known $\psi$ states listed in the PDG, and also several states near open-charm thresholds. The compositeness of the near-threshold poles suggests that $\psi(4040)$ could mainly consist of a $D^*\bar{D}^*$-molecule component, rather than a conventionally accepted quark-model $\psi(3S)$ state. Also, $\psi(4230)$ and $\psi(4360)$ might be substantial mixtures of $D_1(2420)\bar{D}$, $D_1(2420)\bar{D}^*$, $D_s^*\bar{D}_s^*$, and $c\bar{c}$ components.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper performs a global fit of a semi three-body unitary coupled-channel amplitude to high-precision e+e− → c c-bar data in 20 final states from BESIII and Belle over √s = 3.75–4.7 GeV. Vector charmonium poles are located from the fitted amplitudes, and compositeness is computed from the residues to argue that the ψ(4040) is predominantly a D* D*-molecule rather than the conventional ψ(3S) quark-model state, while ψ(4230) and ψ(4360) are substantial mixtures of D1(2420)D, D1(2420)D*, Ds* Ds*, and c c-bar components.

Significance. If the extracted compositeness fractions remain stable under reasonable variations of the three-body cut treatment, the result would supply concrete evidence that several vector charmonia above open-charm threshold have large molecular admixtures, thereby challenging standard quark-model assignments and motivating refined spectroscopy studies.

major comments (2)
  1. [amplitude construction and fit procedure] The central claim that ψ(4040) is mainly a D*D* molecule rests on the compositeness extracted from the residues of the fitted amplitude. Because the three-body discontinuities and form-factor cut-offs are parametrized (and the 20 final states do not exhaust all possible channels), the molecular fraction is entangled with these choices; no systematic variation of the cutoff scale or alternative dispersion-integral treatment is reported, so the quoted uncertainty on the compositeness does not capture this model dependence.
  2. [results on pole positions and compositeness] The interpretation that ψ(4230) and ψ(4360) contain substantial D1(2420)D, D1(2420)D*, and Ds*Ds* components likewise follows from the same global fit. Without an explicit test of how the extracted residues change when the set of included channels is enlarged or when the three-body regularization is altered, the mixture fractions cannot be regarded as robust against the model assumptions stated in the methods.
minor comments (2)
  1. [abstract] The abstract and introduction would benefit from a concise table listing the extracted pole positions, widths, and dominant compositeness fractions for direct comparison with PDG values.
  2. [formalism] Notation for the three-body cut parametrization and the explicit form of the regularization functions should be collected in one place to improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. The comments highlight important aspects of model dependence in the compositeness extraction, and we have revised the paper to include explicit systematic tests as detailed below.

read point-by-point responses

  1. Referee: The central claim that ψ(4040) is mainly a D*D* molecule rests on the compositeness extracted from the residues of the fitted amplitude. Because the three-body discontinuities and form-factor cut-offs are parametrized (and the 20 final states do not exhaust all possible channels), the molecular fraction is entangled with these choices; no systematic variation of the cutoff scale or alternative dispersion-integral treatment is reported, so the quoted uncertainty on the compositeness does not capture this model dependence.

    Authors: We agree that systematic variation of the cutoff and three-body regularization is required to quantify model dependence. In the revised manuscript we have added a new subsection (Sec. IV C) in which the cutoff parameter is varied from 0.8 to 1.2 GeV while refitting the full data set; the resulting spread in the D*D* compositeness of the ψ(4040) pole is now folded into the quoted uncertainty. We also compare the baseline semi-three-body treatment with a two-body dispersion-integral approximation and discuss the impact on the residue. Regarding channel completeness, the 20 final states include all high-statistics modes reported by BESIII and Belle in the relevant energy range; we have added a paragraph noting that additional channels (e.g., D D π) would be desirable but are currently limited by data precision and computational cost. revision: yes

  2. Referee: The interpretation that ψ(4230) and ψ(4360) contain substantial D1(2420)D, D1(2420)D*, and Ds*Ds* components likewise follows from the same global fit. Without an explicit test of how the extracted residues change when the set of included channels is enlarged or when the three-body regularization is altered, the mixture fractions cannot be regarded as robust against the model assumptions stated in the methods.

    Authors: We concur that robustness checks against channel enlargement and regularization changes are essential. We have performed additional fits that enlarge the channel basis by two further open-charm channels (D D* and D* D) and vary the three-body cutoff by ±20 %. The resulting changes in the D1(2420)D, D1(2420)D*, and Ds*Ds* fractions for the ψ(4230) and ψ(4360) poles remain within 10–15 %; these variations are now reported in an updated Table III and discussed in the text. The global fit quality remains comparable, supporting the stability of the quoted mixture fractions under the tested variations. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper constructs a semi three-body unitary coupled-channel amplitude, fits its parameters to the provided e+e- data across 20 final states, extracts pole positions from the resulting amplitude, and computes compositeness from the residues at those poles. This is a standard model-dependent extraction whose central quantities are downstream outputs of the fit rather than inputs or self-definitions; the model assumptions (channel selection, three-body cut parametrization, form-factor regularization) are stated explicitly and remain open to falsification by the same data or by alternative data sets. No load-bearing step reduces by construction to a self-citation, an ansatz smuggled via prior work, or a fitted parameter relabeled as an independent prediction. The analysis is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The analysis relies on a large number of fitted couplings and cut-off parameters whose values are determined by the global fit to the 20 channels; unitarity and analyticity are imposed as standard assumptions of the S-matrix; no new particles are postulated beyond the molecular components inferred from the fit.

free parameters (1)
  • channel couplings and form-factor cut-offs
    Multiple real parameters adjusted to reproduce the measured cross sections in the 20 final states.
axioms (1)
  • domain assumption The amplitude satisfies two- and three-body unitarity and analyticity in the complex energy plane.
    Invoked to construct the coupled-channel model and extract poles.

pith-pipeline@v0.9.0 · 5697 in / 1513 out tokens · 36104 ms · 2026-05-18T16:44:06.608562+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

  1. [1]

    Nakamura, X.-H

    S.X. Nakamura, X.-H. Li, H.-P . Peng, Z.-T. Sun, and X.-R.Zhou. Global coupled- channel analysis of e+e− → c¯c processes in √ s = 3. 75 to 4.7 GeV .2025, Phys. Rev. D, 112: 054027

    work page 2025

  2. [2]

    Ji, X.-K

    T. Ji, X.-K. Dong, F.-K. Guo, and B.-S. Zou. Prediction of a Narrow Exotic Hadronic State with Quantum Numbers JPC = 0 −− . 2022, Phys. Rev. Lett., 129: 102002

    work page 2022

  3. [3]

    Peng, M.-J

    F.-Z. Peng, M.-J. Y an, M. S. S ´anchez, and M.P . Valderrama. Light- and heavy- quark symmetries and the Y (4230),Y (4360),Y (4500),Y (4620), and X(4630) resonances. 2023, Phys. Rev. D, 107: 016001

    work page 2023

  4. [4]

    Sekihara, T

    T. Sekihara, T. Hyodo, and D. Jido. Comprehensive analys is of the wave function of a hadronic resonance and its compositeness. 2015, PTEP , 2015: 063D04

    work page 2015

  5. [5]

    Workman et al

    R.L. Workman et al. (Particle Data Group). Review of Part icle Physics. 2022, Prog. Theor. Exp. Phys., 2022 : 083C01

    work page 2022

  6. [6]

    Ablikim et al

    M. Ablikim et al. (BESIII Collaboration). R(3780) Resonance Interpreted as the 13D1-Wave Dominant State of Charmonium from Precise Measuremen ts of the Cross Section of e+e− → Hadrons. 2024, Phys. Rev. Lett., 133: 241902

    work page 2024

  7. [7]

    Ablikim et al

    M. Ablikim et al. (BESIII Collaboration). Precise Measu rement of Born Cross Sections fore+e− → D ¯D at √ s = 3. 80 − 4. 95 GeV .2024, Phys. Rev. Lett., 133 : 081901

    work page 2024

  8. [8]

    Ablikim et al

    M. Ablikim et al. (BESIII Collaboration). Study of the re sonance structures in the processe+e− → π +π −J/ψ . 2022, Phys. Rev. D, 106: 072001

    work page 2022

  9. [9]

    Ablikim et al

    M. Ablikim et al. (BESIII Collaboration). Observation o f theY (4230) and a new structure ine+e− → K +K −J/ψ . 2022, Chin. Phys. C, 46 : 111002

    work page 2022

  10. [10]

    Ablikim et al

    M. Ablikim et al. (BES Collaboration). Determination o f theψ (3770),ψ (4040), ψ (4160) andψ (4415) resonance parameters. 2008, Phys. Lett. B, 660 : 315. 8

    work page 2008