Pith. sign in

REVIEW 3 major objections 5 minor 48 references

Study on charmonium(-like) mesons within a diabatic approach

T0 review · 3 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A diabatic coupled-channel model with spin-dependent charm quark interactions yields charmonium masses and component probabilities that support molecular assignments for chi_c1(3872), psi(4040), and chi_c2(3930).

desk verdict Useful incremental extension of the diabatic approach, but the chi_c2(3930) molecular claim is contradicted by the authors' own flavor-scaling test. read the letter →

arxiv 2504.12149 v1 pith:UUHWJFNT submitted 2025-04-16 hep-ph

Zi-Zhao Zhang , Rong Li , Bo-Chao Liu This is my paper · ORCID
classification hep-ph
keywords charmoniumstatesspectrumapproachcomponentsdiabaticlikeobtain
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Charmonium is the bound state of a charm quark and a charm antiquark, like a tiny atom. Some measured particles in this mass region do not look like pure quark-antiquark states, so physicists suspect they contain a pair of mesons as well. This paper uses a diabatic formalism, which sets up a matrix of potentials describing the quark-antiquark part and several meson-meson parts together, and solves the coupled equations.

The new ingredient is the inclusion of spin-dependent forces between the charm quark and antiquark. These forces split levels that were degenerate in earlier diabatic calculations. The authors take most parameters from an earlier quark model and fit one extra parameter to the known mass of chi_c1(3872). From the resulting wave functions, they compute masses and the probability of each component. They find that chi_c1(3872) is about 94 percent D D-bar-star, that psi(4040) has a sizable D_s D_s-star component, and that chi_c2(3930) has about a quarter meson-meson content.

The authors note a serious caveat: if the coupling to charmed-strange channels is weakened by a flavor factor, the 3934 MeV 2++ state disappears. So the chi_c2(3930) interpretation is fragile, while the lower states are more robust.

Extended reading notes

Core claim

The central claim is that including spin-dependent c-cbar interactions in the diabatic coupled-channel model reproduces the charmonium spectrum below 4.1 GeV and yields component probabilities identifying chi_c1(3872), psi(4040), and chi_c2(3930) as having significant molecular components, while chi_c0(3860) and psi(3770) are conventional chi_c0(2P) and psi(1D) states. The paper states: 'Our results support the arguments that the chi_c1(3872), psi(4040) and chi_c2(3930) have significant molecular components.'

Load-bearing premise

The calculation assumes the same mixing strength Delta for charmed and charmed-strange meson-meson channels. The authors test a flavor-scaled coupling with a factor mq/ms and find that the 2++ state at 3934 MeV, which they identify with chi_c2(3930), no longer exists. This assumption enters in Section IV.B and the sensitivity check is reported in the final discussion, so the chi_c2(3930) assignment depends on it.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript extends the diabatic coupled-channel approach of Bruschini and González to include spin-dependent c\bar c interactions. The authors construct a potential matrix with a Gaussian mixing potential, fix the mixing strength Δ by fitting the mass of χc1(3872), and solve the coupled-channel Schrödinger equations for J^PC = 0−+, 0++, 1+−, 1++, 1−−, and 2++ states below 4.1 GeV. The calculated masses are compared with PDG data and with previous quark-model and diabatic results, and the component probabilities are used to argue that χc1(3872), ψ(4040), and χc2(3930) have significant molecular components, while χc0(3860) and ψ(3770) correspond to χc0(2P) and ψ(1D). The paper closes with a limitations section that acknowledges the bound-state-only treatment and reports a sensitivity test in which the charmed-strange mixing strength is scaled by m_q/m_s, causing the calculated 2++ state at 3934 MeV to disappear.

Significance. If the central claims were robust, the work would be a useful contribution to the phenomenology of charmonium and charmonium-like states: it provides a unified coupled-channel description below 4.1 GeV, includes spin-dependent effects absent from earlier diabatic studies, and offers explicit, falsifiable predictions for the composition of controversial states such as χc1(3872) and χc0(3860). The spectrum in Table IV is in reasonably good agreement with experiment, and the authors are transparent about the main approximations and about the sensitivity of the 2++ state to the flavor dependence of the mixing strength. However, the abstract and conclusions overstate the certainty of the molecular assignments given the fitted nature of Δ and the results of the authors' own sensitivity test.

major comments (3)
  1. [Section IV.B, final paragraph; Tables IV and V] The central claim that χc2(3930) has a significant molecular component is not supported by the authors' own sensitivity analysis. When the mixing strength for charmed-strange channels is scaled by m_q/m_s, as suggested by the 3P0 model, the calculated 2++ bound state at 3934 MeV, which is identified with χc2(3930) in Table IV, no longer exists. This is explicitly stated in the final discussion. Since the abstract and conclusion list χc2(3930) among the states with significant molecular components, the claim must either be withdrawn or qualified, or the baseline calculation must be redone with a flavor-dependent Δ. As written, the identification is an artifact of the equal-Δ assumption.
  2. [Section IV.A, parameter fitting; Table III] The parameter Δ is obtained by fitting the mass of χc1(3872). Consequently, the 94% DD* probability for the 3871.7 MeV 1++ state in Table V is not an independent prediction of the model: a single mixing parameter tuned to place a bound state at the DD* threshold essentially forces the state to be nearly molecular. The claim that χc1(3872) has a significant molecular component is therefore circular in this model. The authors should present this result as a consistency check and examine how the composition varies with Δ, rather than presenting it as independent evidence for the molecular interpretation.
  3. [Section IV.B, last paragraph; Table V] The stated bound-state-only approximation is not reconciled with the probabilities listed in Table V for channels whose thresholds lie below the computed mass. The text says that mixing with meson pairs at energies above their thresholds is simply ignored, yet Table V reports, for example, a 20% DD* component for the 4060.1 MeV 1−− state even though the DD* threshold is 3872 MeV. The boundary conditions used for open channels need to be specified; as it stands, the molecular probability for ψ(4040) is not well defined. This affects the abstract claim that ψ(4040) has a significant molecular component.
minor comments (5)
  1. [Title] The title contains a typo: 'approa ch' should read 'approach'.
  2. [Section IV.B] The reference for the m_q/m_s factor in the 3P0 model is missing in the manuscript (the text contains '[?]'); this should be supplied.
  3. [Table V] The column alignment of Table V is difficult to read in the preprint; explicit column headers and explicit zero entries would make the component probabilities clearer.
  4. [Abstract and Section V] The abstract states that χc2(3930) has a significant molecular component, while the conclusion says it 'may have' one; the strength of this claim should be made consistent with the sensitivity analysis.
  5. [Section III.A] The regularization of the 1/r^3 terms in the spin-orbit and tensor potentials by the cutoff r_cut is mentioned but the explicit regularized forms are not given; writing them out would improve reproducibility.
Assumptions & free parameters 7 free parameters · 6 assumptions · 0 invented entities

The central calculation depends on seven fixed or fitted numbers in the potential matrix and several structural approximations. Only Delta is fitted in this paper; the others come from a prior quark model and from earlier diabatic studies. The Gaussian mixing ansatz and the neglect of meson-meson interactions are the most model-dependent choices.

free parameters (7)
  • Delta = 0.116 GeV
    Fitted in this work to the mass of chi_c1(3872); sets the overall cc to meson-meson mixing strength.
  • rho = 0.3 fm
    Taken from Ref [24]; width of the Gaussian mixing potential, argued from lattice QCD results.
  • mc = 1.4830 GeV
    Charm quark mass from the quark model fit of Ref [36] to 12 charmonium states.
  • alpha_s = 0.5461
    Strong coupling in the Cornell potential from Ref [36].
  • b = 0.1425 GeV^2
    String tension in the confining potential from Ref [36].
  • sigma = 1.1384 GeV
    Gaussian smearing width for the spin-spin contact term from Ref [36].
  • rcut = 0.202 fm
    Cutoff for the 1/r^3 singularity in the spin-orbit and tensor potentials from Ref [36].
assumptions (6)
  • domain assumption Born-Oppenheimer static limit and diabatic expansion at a fixed r0 with light-field states evaluated at r0
    Section II, Eq. (4) expands the wavefunction in states at fixed r0; this assumes the heavy-quark separation can be treated as a c-number parameter.
  • domain assumption Interactions between different meson-meson channels are neglected (V_ij = 0 for i != j)
    Section III.B, Eq. (28), justified by analogy to lattice QCD studies.
  • domain assumption Off-diagonal tensor couplings between different L,S states are neglected
    Section III.A, stated as 'very small and can be neglected'.
  • ad hoc to paper The mixing potential V_mix has a Gaussian form in (Vcc - T)
    Section III.C, Eq. (31); a phenomenological ansatz rather than a derived result.
  • ad hoc to paper The same coupling Delta is used for charmed and charmed-strange meson channels in the baseline calculation
    Section IV.B, final discussion; the authors test an alternative but use the uniform Delta for the main results.
  • domain assumption Only bound states below meson-meson thresholds are treated; mixing above thresholds is neglected
    Section IV.B, final discussion; limits the validity near and above open-charm thresholds.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Study on charmonium(-like) mesons within a diabatic approach." pith.science (2026). https://pith.science/paper/UUHWJFNT

@misc{pith2026250412149,
  author       = {Pith},
  title        = {Pith review of: Study on charmonium(-like) mesons within a diabatic approach},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UUHWJFNT}},
  note         = {Machine review of arXiv:2504.12149}
}
abstract

In this work, we study the charmonium(-like) spectrum below 4.1 GeV using the diabatic approach, which offers a unified description of conventional and unconventional heavy meson states. Compared to previous studies, we consider a more realistic $c\bar c$ potential with including the spin-dependent interactions, which allows us to obtain more states and get more insights on the charmonium spectrum. Based on our calculation, we obtain the masses of the charmonium spectrum which align with the experimental data well. We also present the probabilities of finding various components, i.e. $c\bar c$ or meson-meson pair, in those states. Our results support the arguments that the $\chi_{c1}(3872)$, $\psi(4040)$ and $\chi_{c2}(3930)$ have significant molecular components. In addition, our calculations show that the $\chi_{c0}(3860)$ and $\psi(3770)$ can be looked as the candidates for the charmonium states $\chi_{c0}(2P)$ and $\psi(1D)$, respectively.

Figures

Figures reproduced from arXiv: 2504.12149 by the authors.

Figure 1
Figure 1. FIG. 1: Radial wave function of the calculated 0 [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Radial wave function of the calculated 0 [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4: Radial wave function of the calculated 0 [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figures from the paper (2 more)
Figure 3
Figure 3. Figure 3: FIG. 3: Radial wave function of the calculated 0 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png]
Figure 5
Figure 5. Figure 5: FIG. 5: Radial wave function of the calculated 0 [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

48 extracted references · 36 canonical work pages

  1. [1]

    For the 0 ++ state χc0(3860)[30, 39], our calculated mass is 3862 MeV, which is exactly same as the cen- tral value suggested by PDG

    The rms radius can be obtained as 6.55 fm, which is significantly large compared to its c c component with a rms radius of 1.15 fm. For the 0 ++ state χc0(3860)[30, 39], our calculated mass is 3862 MeV, which is exactly same as the cen- tral value suggested by PDG. Its radial wave functions of various components are drawn in Fig

    work page

  2. [2]

    The probabilities of various /s48 /s49 /s50 /s51 /s45/s48/s46/s52 /s45/s48/s46/s50 /s48/s46/s48 /s48/s46/s50 /s48/s46/s52 /s48/s46/s54 /s32 /s32 /s32 FIG

    The result- ing rms radius is 1.05 fm. The probabilities of various /s48 /s49 /s50 /s51 /s45/s48/s46/s52 /s45/s48/s46/s50 /s48/s46/s48 /s48/s46/s50 /s48/s46/s52 /s48/s46/s54 /s32 /s32 /s32 FIG. 1: Radial wave function of the calculated 0 +(1++) state with a mass of 3871.7 MeV. /s48 /s49 /s50 /s51 /s45/s49 /s48 /s49 /s50 /s51 /s32 /s32 /s32 /s32 FIG. 2: Ra...

    work page

  3. [3]

    G. S. Abrams, D. Briggs, W. Chinowsky, C. E. Friedberg, G. Goldhaber, R. J. Hollebeek, J. A. Kadyk, A. Litke, B. Lulu, F. Pierre, et al., Phys. Rev. Lett., v. 33, no. 24, pp. 1453-1455 (1974)

    work page 1974

  4. [4]

    (6) The introduction of the diabatic potential matrix is a crucial step in the diabatic approach

    into Eq.( 2) and multiplying ⟨r| and ⟨ζj(r0)| on the left side of the equation give ∑ i ( − ℏ2 2µδji∇2 +Vji(r, r0) −Eδji ) ˜ψi(r, r0) = 0, (5) with the diabatic potential matrix being defined as Vji(r, r0) ≡ ⟨ζj (r0)|H light static(r)|ζi(r0)⟩. (6) The introduction of the diabatic potential matrix is a crucial step in the diabatic approach. Since the r0 is ...

    work page

  5. [5]

    To understand the properties of these states, vari- ous potential models were proposed to describe the inter- actions between quarks [ 8, 9]

    and others[ 6, 7] were successively discovered in experi- ments. To understand the properties of these states, vari- ous potential models were proposed to describe the inter- actions between quarks [ 8, 9]. Among them, the Cornell potential [ 10] is widely adopted to predict and interpret the charmonium spectrum. While the Cornell potential is rather succ...

    work page

  6. [6]

    (32) As mentioned above, the mixing of c¯c and M ¯M pri- marily occurs near the crossing radius rc

    as V0i(r) = Vi0(r) = V mix i (r), i ⁄= 0. (32) As mentioned above, the mixing of c¯c and M ¯M pri- marily occurs near the crossing radius rc. In this region, the long-range part of Vcc(r)(Eq.27) will play a major role in the mixing potential V mix i (r). In our model, the spin-dependent terms are mainly short-range, and the mixing of different c¯c states d...

    work page 2007

  7. [7]

    J. J. Aubert, U. Becker, P. J. Biggs, J. Burger, M. Chen, G. Everhart, P. Goldhagen, J. Leong, T. McCorriston, T. G. Rhoades, et al., Phys. Rev. Lett. 33, 1404 (1974)

    work page 1974

  8. [8]

    J. E. Augustin et al., Phys. Rev. Lett. 33, 1406 (1974)

    work page 1974

Show all 48 references
  1. [9]

    C. J. Biddick, T. H. Burnett, G. E. Masek, E. S. Miller, J. G. Smith, J. P. Stronski, M. K. Sullivan, W. Vernon, D. H. Badtke, B. A. Barnett, et al., Phys. Rev. Lett. 38, 1324 (1977)

    1977

  2. [10]

    J. S. Whitaker, W. Tanenbaum, G. S. Abrams, M. S. Alam, A. M. Boyarski, M. Breidenbach, W. Chinowsky, R. G. DeVoe, G. J. Feldman, C. E. Friedberg, et al., Phys. Rev. Lett. 37, 1596 (1976)

    1976

  3. [11]

    Tanenbaum, J

    W. Tanenbaum, J. S. Whitaker, G. S. Abrams, A. M. Boyarski, M. Breidenbach, F. Bulos, W. Chinowsky, G. J. Feldman, C. E. Friedberg, D. Fryberger, et al., Phys. Rev. Lett. 35, 1323 (1975)

    1975

  4. [12]

    Partridge, C

    R. Partridge, C. Peck, F. Porter, W. Kollmann, M. Richardson, K. Strauch, K. Wacker, D. Aschman, T. Burnett, M. Cavalli-Sforza, et al., Phys. Rev. Lett. 45, 1150 (1980)

    1980

  5. [13]

    Eichten, K

    E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T. M. Yan, Phys. Rev. D 17, 3090 (1978)

    1978

  6. [14]

    Voloshin, Progress in Particle and Nuclear Physics 61, 455 (2008)

    M. Voloshin, Progress in Particle and Nuclear Physics 61, 455 (2008)

    2008

  7. [15]

    Eichten, K

    E. Eichten, K. Gottfried, T. Kinoshita, J. B. Kogut, K. D . Lane, and T.-M. Yan, Phys. Rev. Lett. 34, 369 (1975), [Erratum: Phys.Rev.Lett. 36, 1276 (1976)]

    1975

  8. [16]

    Brambilla, S

    N. Brambilla, S. Eidelman, C. Hanhart, A. Nefediev, C.- P. Shen, C. E. Thomas, A. Vairo, and C.-Z. Yuan, Phys. Rept. 873, 1 (2020), 1907.07583

    2020 arXiv

  9. [17]

    R. L. Workman and Others (Particle Data Group), PTEP 2022, 083C01 (2022)

    2022

  10. [18]

    F.-K. Guo, C. Hanhart, U.-G. Meißner, Q. Wang, Q. Zhao, and B.-S. Zou, Rev. Mod. Phys. 90, 015004 (2018)

    2018

  11. [19]

    Ji, X.-K

    T. Ji, X.-K. Dong, M. Albaladejo, M.-L. Du, F.-K. Guo, J. Nieves, and B.-S. Zou, Sci. Bull. 68, 688 (2023), 2212.00631

    2023 arXiv

  12. [20]

    Hanhart and E

    C. Hanhart and E. Klempt, Int. J. Mod. Phys. A 35, 2050019 (2020), 1906.11971

    2020 arXiv

  13. [21]

    H.-X. Chen, W. Chen, X. Liu, Y.-R. Liu, and S.-L. Zhu, Rept. Prog. Phys. 86, 026201 (2023), 2204.02649

    2023 arXiv

  14. [22]

    B.-Q. Li, C. Meng, and K.-T. Chao, Phys. Rev. D 80, 014012 (2009)

    2009

  15. [23]

    Y. S. Kalashnikova, Phys. Rev. D 72, 034010 (2005)

    2005

  16. [24]

    M. R. Pennington and D. J. Wilson, Phys. Rev. D 76, 077502 (2007), URL https://link.aps.org/doi/10.1103/PhysRevD.76.077502

    2007 doi

  17. [25]

    Ono and N

    S. Ono and N. T¨ ornqvist, Zeitschrift f¨ ur Physik C Parti- cles and Fields 23, 59 (1984)

    1984

  18. [26]

    Born and W

    M. Born and W. Heisenberg, Original Scientific Papers Wissenschaftliche Originalarbeiten pp. 216–246 (1985)

    1985

  19. [27]

    G. S. Bali, Physics Reports 343, 1 (2001)

    2001

  20. [28]

    J. T. Castella, Physical Review D 106, 094020 (2022)

    2022

  21. [29]

    Bruschini and P

    R. Bruschini and P. Gonz´ alez, Phys. Rev. D 102, 074002 (2020)

    2020

  22. [30]

    R. F. Lebed and S. R. Martinez, Phys. Rev. D 106, 074007 (2022), 2207.01101

    2022 arXiv

  23. [31]

    R. F. Lebed and S. R. Martinez (2024), 2404.15186

    2024 arXiv

  24. [32]

    Baer, Beyond Born-Oppenheimer: Electronic Nonadi- abatic Coupling Terms and Conical Intersections, (John Wiley & Sons, New York) (2006)

    M. Baer, Beyond Born-Oppenheimer: Electronic Nonadi- abatic Coupling Terms and Conical Intersections, (John Wiley & Sons, New York) (2006)

    2006

  25. [33]

    Bruschini and P

    R. Bruschini and P. Gonz´ alez, Physical Review D 103 (2021)

    2021

  26. [34]

    A. A. Petrov, J. Phys. Conf. Ser. 9, 83 (2005)

    2005

  27. [35]

    G. L. Yu, Z. G. Wang, and Z. Y. Li, Chin. Phys. C 42, 4 (2018), 1704.06763

    2018 arXiv

  28. [36]

    Bulava, B

    J. Bulava, B. H¨ orz, F. Knechtli, V. Koch, G. Moir, C. Morningstar, and M. Peardon, Physics Letters B 793, 493–498 (2019)

    2019

  29. [37]

    Eichten, S

    E. Eichten, S. Godfrey, H. Mahlke, and J. L. Rosner, Reviews of Modern Physics 80, 1161 (2008)

    2008

  30. [38]

    Lucha and F

    W. Lucha and F. F. Schoberl, in International Summer School for Students on Development in Nuclear Theory and Particle Physics (1995), hep-ph/9601263

    1995 arXiv

  31. [39]

    Barnes, S

    T. Barnes, S. Godfrey, and E. Swanson, Physical Review D 72, 054026 (2005)

    2005

  32. [40]

    Heikkil¨ a, N

    K. Heikkil¨ a, N. A. T¨ ornqvist, and S. Ono, Phys. Rev. D 29, 110 (1984)

    1984

  33. [41]

    W.-J. Deng, H. Liu, L.-C. Gui, and X.-H. Zhong, Phys. Rev. D 95, 034026 (2017)

    2017

  34. [42]

    B.R.johnson, The Journal of Chemical Physics 69, 4678 (1978)

    1978

  35. [43]

    Navas et al

    S. Navas et al. (Particle Data Group), Phys. Rev. D 110, 030001 (2024)

    2024

  36. [44]

    Chilikin et al

    K. Chilikin et al. (Belle), Phys. Rev. D 95, 112003 (2017), 1704.01872

    2017 arXiv

  37. [45]

    P. G. Ortega, J. Segovia, D. R. Entem, and F. Fern´ andez, Physics Letters B 778, 1 (2018)

    2018

  38. [46]

    Z.-Y. Zhou, Z. Xiao, and H.-Q. Zhou, Phys. Rev. Lett. 115, 022001 (2015)

    2015

  39. [47]

    Guo and U.-G

    F.-K. Guo and U.-G. Meissner, Phys. Rev. D 86, 091501 (2012), 1208.1134

    2012 arXiv

  40. [48]

    Aaij et al

    R. Aaij et al. (LHCb), Phys. Rev. Lett. 125, 242001 (2020), 2009.00025

    2020

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.