arxiv: 2604.10626 · v1 · submitted 2026-04-12 · ✦ hep-ph · hep-ex· hep-lat

Recognition: unknown

Resonance X(6600)

S. S. Agaev , K. Azizi , H. Sundu

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:36 UTC · model grok-4.3

classification ✦ hep-ph hep-exhep-lat

keywords all-charm tetraquarkX(6600)QCD sum rulesdiquark-antidiquarktensor statedecay widthscharmoniumexotic hadron

0 comments

The pith

QCD sum rules identify an all-charm tensor tetraquark as a main component of the X(6600) resonance.

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

The paper models the resonance X(6600) as an all-charm tetraquark with J^{PC}=2^{++} in the diquark-antidiquark picture. It builds an interpolating current from an axial-vector diquark and antidiquark, each made of two charm quarks, and applies two-point QCD sum rules to extract the mass while using three-point sum rules for partial decay widths. The calculation produces a mass of 6609 MeV and a total width of 165 MeV that match experimental values for X(6600). This agreement leads the authors to conclude that the tetraquark state forms an essential part of the observed resonance. The work also derives a lower limit on the mass of the state's first radial excitation.

Core claim

The resonance X(6600) is explored as the all-charm tetraquark structure with spin-parities J^{PC}=2^{++}. It is considered in the diquark-antidiquark picture and modeled as a tensor state X composed of the axial-vector diquark cCγ_μ c and antidiquark c̄ γ_ν C c̄. The mass and decay width of X are evaluated in the framework of QCD sum rule methods. The two-point SR approach is applied to find its spectroscopic parameters, while three-point SRs are used to calculate partial widths of different decay channels of X. Comparison of the mass m=(6609 ± 50) MeV and width Γ[X]=(165 ± 23) MeV of the tensor diquark-antidiquark state X with experimental data allows us to interpret it as an essential of

What carries the argument

The axial-vector diquark-antidiquark interpolating current for the J^{PC}=2^{++} all-charm tetraquark, used in QCD two-point and three-point sum rules to compute the mass and decay widths.

If this is right

Leading decays proceed to J/ψ J/ψ, η_c η_c, and χ_c1(1P) η_c channels where all four charm quarks form final-state mesons.
Subleading decays to D_{(s)} meson pairs arise from c c̄ annihilation inside the tetraquark.
The first radial excitation of the state has a mass above the lower limit obtained from the sum-rule analysis.
The computed width and mass together support identifying the calculated state as an essential component of the observed X(6600).

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the assignment holds, similar diquark-antidiquark currents could be applied to predict masses and widths of other fully charmed or bottomed exotic states.
Branching-fraction measurements in the predicted charmonium-pair channels at current or future colliders would provide a direct test of the model.
The approach may help distinguish pure tetraquark interpretations from molecular or hybrid alternatives for nearby resonances.
Observation of the radial excitation above the stated lower mass limit would further support the overall picture.

Load-bearing premise

The chosen diquark-antidiquark interpolating current accurately represents the physical state without significant mixing from other configurations, and the QCD sum rule approximations hold reliably for this heavy tetraquark system.

What would settle it

A precise experimental determination of the X(6600) mass or width lying outside the calculated ranges of 6609 ± 50 MeV and 165 ± 23 MeV, or the absence of dominant decays into J/ψ J/ψ, would falsify the tetraquark assignment.

Figures

Figures reproduced from arXiv: 2604.10626 by H. Sundu, K. Azizi, S. S. Agaev.

**Figure 2.** Figure 2: FIG. 2: Dependence of the mass [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3: Left panel: QCD data and fit function [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4: Left panel: QCD data and fit functions for strong coupl [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗

read the original abstract

The resonance $X(6600)$ is explored as the all-charm tetraquark structure with spin-parities $J^{\mathrm{PC}}=2^{++}$. It is considered in the diquark-antidiquark picture and modeled as a tensor state $X$ composed of the axial-vector diquark $cC\gamma_{\mu}c$ and antidiquark $\overline{c}% \gamma_{\nu}C\overline{c}$ with $C$ being the charge conjugation matrix. The mass and decay width of $X$ are evaluated in the framework of QCD sum rule (SR) methods. The two-point SR approach is applied to find its spectroscopic parameters, while three-point SRs used to calculate partial widths of different decay channels of $X$. We study its leading decays $X \to J/\psi J/\psi$, $X \to \eta_{c}\eta_{c}$ and $\chi _{c1}(1P)\eta _{c}$ in which all four $c$-quarks constitute final-state mesons. We consider also the subleading channels $X\to D_{(s)}^{(\ast )+}D_{(s)}^{(\ast )-}$ and $% D_{(s)}^{(\ast )0}\overline{D}_{(s)}^{(\ast )0}$ generated by annihilation of $\overline{c}c$ quarks in the tetraquark. Comparison of the mass $m=(6609 \pm 50)~ \mathrm{MeV}$ and width $\Gamma[X]=(165 \pm 23)~ \mathrm{MeV}$ of the tensor diquark-antidiquark state $X$ with experimental data allows us to interpret it as an essential component of the resonance $X(6600)$. We also provide a lower limit for the mass of the first radial excitation of $X$.

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.

Desk Editor's Note private letter to a colleague

Standard QCD sum rule calculation assigns X(6600) as 2++ all-charm tetraquark with mass and width match, but heavy-quark OPE issues remain the usual caveat.

read the letter

This paper models X(6600) as a 2++ tensor state built from axial-vector diquark and antidiquark, then runs two-point sum rules for the mass and three-point ones for the widths. The central numbers come out at 6609 ± 50 MeV and 165 ± 23 MeV, which they compare directly to the observed resonance, and they add a lower limit on the first radial excitation mass. They also compute partial widths for the leading channels J/ψJ/ψ, ηcηc, χc1ηc plus the subleading open-charm modes from cc annihilation. That set of explicit decay calculations is the concrete new piece; the rest follows the same diquark-antidiquark framework they and others have used before for lighter tetraquarks. The numerics are presented with the standard variation over Borel window and continuum threshold, so the error bars are at least internally consistent. The soft spot is exactly the one the stress-test flags: four heavy quarks make gluon condensates larger and the Borel window narrower, so OPE truncation after dimension 6 or 8 and the assumption that the chosen current has dominant overlap are less secure than in lighter systems. Without the stability plots and explicit higher-dimension terms in the full text it is hard to judge how much those choices affect the central values. The claim is therefore a reasonable fit within the method, but the method itself carries larger systematic uncertainty here than the quoted errors suggest. This is for the narrow group that follows sum-rule work on exotics. A reader already comfortable with the technique will find a competent, self-contained application and a testable prediction. I would send it to referees; the calculation is detailed enough that a review can check the parameter windows and mixing assumptions directly.

Referee Report

3 major / 2 minor

Summary. The paper claims that the resonance X(6600) can be interpreted as an all-charm tetraquark state with J^PC=2++ in the diquark-antidiquark picture. Using two-point QCD sum rules, the mass is calculated as m=(6609±50) MeV, and three-point sum rules yield a total width Γ[X]=(165±23) MeV by considering decays to J/ψJ/ψ, η_cη_c, χ_c1(1P)η_c, and subleading channels involving D mesons. This agreement with experimental data supports the interpretation, and a lower limit for the mass of the first radial excitation is provided.

Significance. If the QCD sum rule results are robust, this provides a concrete theoretical interpretation for the observed X(6600) resonance as a tensor tetraquark, aiding in the classification of exotic hadrons. The explicit computation of partial decay widths for several channels offers predictions that can be tested experimentally. The inclusion of a lower bound for the radial excitation mass adds value by extending the analysis beyond the ground state. Strengths include the systematic use of two- and three-point sum rules for both spectroscopy and decays.

major comments (3)

[Two-point SR analysis] The mass prediction m=(6609±50) MeV and width Γ=(165±23) MeV rely on the stability of the Borel parameter M^2 and continuum threshold s0 in the two-point SR. In all-charm tetraquark systems, gluon condensates scale with m_c^4 and the Borel window is typically narrow; the manuscript must demonstrate that OPE truncation after dimension-6/8 terms does not introduce systematics larger than the quoted uncertainties, as this directly affects whether the central values can be compared to data.
[Three-point SR for decays] The total width is obtained from three-point SRs assuming the axial-vector diquark-antidiquark interpolating current has dominant overlap with the physical state. This assumption is load-bearing for the claim that the calculated width matches experiment and allows interpretation as a component of X(6600); the paper should provide a sensitivity study or argument ruling out large mixing with meson-meson or alternative tetraquark currents.
[Comparison with experiment] The comparison to experimental mass and width data for the interpretation of X(6600) carries moderate circular content because s0 and M^2 are chosen for stability and can indirectly incorporate experimental information. The manuscript should quantify how independent the final m and Γ predictions remain from this choice.

minor comments (2)

[Abstract] The abstract contains LaTeX formatting artifacts (e.g., stray '%' characters and incomplete subscripts in the current definition) that should be corrected for readability.
[General] Ensure consistent numbering and referencing of all equations and tables throughout the text.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and outline the changes we will make in the revised version.

read point-by-point responses

Referee: [Two-point SR analysis] The mass prediction m=(6609±50) MeV and width Γ=(165±23) MeV rely on the stability of the Borel parameter M^2 and continuum threshold s0 in the two-point SR. In all-charm tetraquark systems, gluon condensates scale with m_c^4 and the Borel window is typically narrow; the manuscript must demonstrate that OPE truncation after dimension-6/8 terms does not introduce systematics larger than the quoted uncertainties, as this directly affects whether the central values can be compared to data.

Authors: We agree that explicit demonstration of OPE convergence is important for all-charm systems. In the original analysis the Borel window was fixed by the standard criteria of pole dominance (>50%) and OPE convergence (highest retained term <5% of the sum), with dimension-6 and dimension-8 gluon condensates included. The quoted uncertainties already incorporate the variation of M^2 and s0 inside that window. To make the truncation systematics fully transparent we will add a supplementary table (or figure) that lists the fractional contribution of each OPE term evaluated at the central working point; this will show that omitted higher-dimensional operators remain well below the quoted error. The central mass and width values themselves are unchanged. revision: partial
Referee: [Three-point SR for decays] The total width is obtained from three-point SRs assuming the axial-vector diquark-antidiquark interpolating current has dominant overlap with the physical state. This assumption is load-bearing for the claim that the calculated width matches experiment and allows interpretation as a component of X(6600); the paper should provide a sensitivity study or argument ruling out large mixing with meson-meson or alternative tetraquark currents.

Authors: The current was chosen because it carries the correct J^PC=2++ quantum numbers and is the standard interpolator in the diquark-antidiquark picture for a tensor state. While a quantitative mixing analysis with molecular or other tetraquark currents would require a coupled-channel or lattice-QCD treatment outside the present scope, we note that the same current reproduces both the mass and the dominant decay channels consistently. In the revision we will insert a short paragraph discussing why alternative currents (e.g., color-singlet meson-meson operators) would produce different decay hierarchies or mass predictions that are less compatible with the observed X(6600) data. A full numerical sensitivity scan, however, is not feasible within the present framework. revision: partial
Referee: [Comparison with experiment] The comparison to experimental mass and width data for the interpretation of X(6600) carries moderate circular content because s0 and M^2 are chosen for stability and can indirectly incorporate experimental information. The manuscript should quantify how independent the final m and Γ predictions remain from this choice.

Authors: The continuum threshold s0 is fixed by the phenomenological expectation s0 ≈ (m_X + 0.5 GeV)^2 for the first radial excitation, and M^2 is chosen solely from the requirements of pole dominance and OPE convergence; neither parameter is adjusted to reproduce the experimental mass of X(6600). The resulting mass is therefore a genuine prediction. To quantify the residual dependence we will add a short subsection (or an appendix figure) that displays the predicted mass and total width as functions of M^2 and s0 throughout the allowed window. The plot will demonstrate that the central values remain stable inside the quoted uncertainties and that the agreement with experiment is not the result of fine-tuning. revision: partial

Circularity Check

0 steps flagged

QCD sum-rule mass and width for the all-charm tensor tetraquark are computed from standard OPE inputs without reduction to experimental X(6600) data

full rationale

The derivation proceeds from the two-point correlator for the chosen diquark-antidiquark current, applies the OPE up to dimension-8 or -10 terms with known quark masses and condensates, imposes quark-hadron duality with a continuum threshold s0 chosen for Borel stability, and extracts m and the decay constants. Three-point sum rules for the widths follow the same pattern with additional meson distribution amplitudes. None of these steps inserts the experimental mass or width of X(6600) as an input; the comparison occurs only after the numerical results are obtained. No self-citation chain supplies a uniqueness theorem or ansatz that forces the final numbers, and the Borel/s0 window is fixed by internal stability rather than by fitting to the resonance under study. The calculation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The central claim rests on standard QCD SR machinery plus the postulated tetraquark structure; free parameters include Borel mass and continuum threshold chosen for stability, while the diquark-antidiquark current is an invented modeling choice without independent falsifiable evidence beyond the mass-width agreement.

free parameters (2)

Borel parameter M^2
Selected within a stability window to suppress higher states and continuum contributions.
Continuum threshold s0
Determined from mass stability or duality requirements.

axioms (2)

domain assumption Quark-hadron duality
Equates OPE side to hadronic spectral density above the continuum threshold.
standard math Convergence of operator product expansion
Truncation at finite dimension assumed sufficient for numerical accuracy.

invented entities (1)

Axial-vector diquark-antidiquark tetraquark X with J^PC=2++ no independent evidence
purpose: To represent the internal structure of the resonance as a bound state of cC gamma_mu c and conjugate antidiquark.
Postulated configuration chosen to match quantum numbers; no independent experimental handle such as a predicted production cross-section or decay branching ratio outside the mass-width fit.

pith-pipeline@v0.9.0 · 5634 in / 1623 out tokens · 71871 ms · 2026-05-10T16:36:25.574289+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Predictions for the scalar partner of the LHC tetraquark $X(6600)$
hep-ph 2026-05 accept novelty 6.0

The paper predicts a scalar cc-ccbar tetraquark state (X(6400)) and identifies it as the partner to the recently observed tensor state X(6600).

Reference graph

Works this paper leans on

50 extracted references · 3 canonical work pages · cited by 1 Pith paper · 1 internal anchor

[1]

Resonance $X(6600)$

44− 6. 47 GeV [6]. In Ref. [9] fully-heavy states were an- alyzed in the diquark-antidiquark and molecular frame- works, and the authors interpreted X(6900) as a hadronic molecule χc0χc0 or/and a tetraquark composed of pseu- doscalar components. It should be noted that alternative explanations of the observed structures were suggested as well. Indeed, new...

work page internal anchor Pith review Pith/arXiv arXiv 2026
[2]

The function Π 1(M2, s0,q 2) is determined as Π 1(M2, s0,q 2) = ∫ s0 16m2 c ds ∫ s′ 0 4m2 c ds′ρ1(s,s ′,q 2) ×e− s/M 2 1 − s′/M 2

is related to J/ψ channel. The function Π 1(M2, s0,q 2) is determined as Π 1(M2, s0,q 2) = ∫ s0 16m2 c ds ∫ s′ 0 4m2 c ds′ρ1(s,s ′,q 2) ×e− s/M 2 1 − s′/M 2
[3]

Our calculations prove that regions Eq

(23) Constraints imposed on M2 and s0 are standard for the sum rule investigations and have been detailed in the previous section. Our calculations prove that regions Eq. (13) for the parameters ( M 2 1,s 0) and M 2 2 ∈ [4, 5] GeV 2, s ′ 0 ∈ [12, 13] GeV 2. (24) for (M 2 2,s ′
[4]

The sum rule for the form factor g1(q2) is applicable in the region q2< 0

meet these conditions. The sum rule for the form factor g1(q2) is applicable in the region q2< 0. But g1(q2) determines the coupling g1 at the mass shell q2 =m2 J/ψ . For that reason, we employ the function g1(Q2) where Q2 = −q2 and apply it in our studies. The SR predictions for g1(Q2) are plotted in Fig. 3, where Q2 changes in the interval Q2 = 2 − 20 G...

2010
[5]

1) MeV are the mass and decay constant of mesons D± [39, 44], respectively. After some manipulations, one gets ˆΠ Phys µν (p,p ′) = g5(q2)Λf 2 Dm4 D m2 c (p2 − m2) (p′2 − m2 D) (q2 − m2 D) × [m4 − 2m2(m2 D +q2) + (m2 D − q2)2 12m2 gµν +p′ µp′ ν + other contributions ] . (67) For ˆΠ OPE µν (p,p ′), we ﬁnd ˆΠ OPE µν (p,p ′) = − 4 3 ∫ d4xd4yeip′ye− ipx⟨cc⟩ ×...

1969
[6]

Aaij et al

R. Aaij et al. (LHCb Collaboration), Sci. Bull. 65, 1983 (2020)

1983
[7]

Aad et al

G. Aad et al. (ATLAS Collaboration), Phys. Rev. Lett. 131, 151902 (2023)

2023
[8]

Hayrapetyan et al

A. Hayrapetyan et al. (CMS Collaboration), Phys. Rev. Lett. 132, 111901 (2024)

2024
[9]

Hayrapetyanet al.(CMS), (2026), arXiv:2602.02252 [hep-ex]

A. Hayrapetyan et al. (CMS Collaboration), arXiv:2602.02252 [hep-ex]

work page arXiv
[10]

M. A. Bedolla, J. Ferretti, C. D. Roberts and E. San- topinto, Eur. Phys. J. C 80, 1004 (2020)

2020
[11]

J. R. Zhang, Phys. Rev. D 103, 014018 (2021)

2021
[12]

Z. G. Wang, Chin. Phys. C 44, 113106 (2020)

2020
[13]

Z. G. Wang, Int. J. Mod. Phys. A 36, 2150014 (2021)

2021
[14]

R. M. Albuquerque, S. Narison, A. Rabemananjara, D. Rabetiarivony, and G. Randriamanatrika, Phys. Rev. D 102, 094001 (2020)

2020
[15]

B. C. Yang, L. Tang, and C. F. Qiao Eur. Phys. J. C 81, 324 (2021)

2021
[16]

M. C. Gordillo, F. De Soto, and J. Segovia, Phys. Rev. D 102, 114007 (2020)

2020
[17]

X. K. Dong, V. Baru, F. K. Guo, C. Hanhart, and A. Nefediev, Phys. Rev. Lett. 126, 132001 (2021); 127, 119901(E) (2021)

2021
[18]

X. K. Dong, V. Baru, F. K. Guo, C. Hanhart, A. Nefediev, and B. S. Zou, Sci. Bull. 66, 2462 (2021)

2021
[19]

Z. R. Liang, X. Y. Wu, and D. L. Yao, Phys. Rev. D 104, 034034 (2021)

2021
[20]

G. J. Wang, L. Meng, M. Oka, and S. L. Zhu, Phys. Rev. D 104, 036016 (2021)

2021
[21]

C. Deng, H. Chen, and J. Ping, Phys. Rev. D 103, 014001 (2021)

2021
[22]

Z. G. Wang, Nucl. Phys. B 985, 115983 (2022)

2022
[23]

R. N. Faustov, V. O. Galkin, and E. M. Savchenko, Sym- metry 14, 2504 (2022)

2022
[24]

P. Niu, Z. Zhang, Q. Wang, and M. L. Du, Sci. Bull. 68, 800 (2023)

2023
[25]

W. C. Dong and Z. G. Wang, Phys. Rev. D 107, 074010 (2023)

2023
[26]

G. L. Yu, Z. Y. Li, Z. G. Wang, J. Lu, and M. Yan, Eur. Phys. J. C 83, 416 (2023)

2023
[27]

H. T. An, S. Q. Luo, Z. W. Liu, and X. Liu, Eur. Phys. J. C 83, 740 (2023)

2023
[28]

S. Q. Kuang, Q. Zhou, D. Guo, Q. H. Yang, and L. Y. Dai, Eur. Phys. J. C 83, 383 (2023)

2023
[29]

M. S. liu, F. X. Liu, X. H. Zhong and Q. Zhao, Phys. Rev. D 109, 076017 (2024)

2024
[30]

Malekhosseini, S

M. Malekhosseini, S. Rostami, A. R. Olamaei and K. Az- izi, Nucl. Phys. B 1018, 116977 (2025)

2025
[31]

S. S. Agaev, K. Azizi, B. Barsbay, and H. Sundu, Phys. Lett. B 844, 138089 (2023)

2023
[32]

S. S. Agaev, K. Azizi, B. Barsbay and H. Sundu, Eur. Phys. J. Plus 138, 935 (2023)

2023
[33]

S. S. Agaev, K. Azizi, B. Barsbay and H. Sundu, Nucl. Phys. A 844, 122768 (2024)

2024
[34]

S. S. Agaev, K. Azizi, B. Barsbay and H. Sundu, Eur. Phys. J. C 83, 994 (2023)

2023
[35]

Hayrapetyan et al

A. Hayrapetyan et al. (CMS Collaboration), Nature 648, 58 (2025)

2025
[36]

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979)

1979
[37]

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 448 (1979)

1979
[38]

R. M. Albuquerque, J. M. Dias, K. P. Khemchandani, A. Martinez Torres, F. S. Navarra, M. Nielsen and C. M. Zanetti, J. Phys. G 46, 093002 (2019)

2019
[39]

S. S. Agaev, K. Azizi, and H. Sundu, Turk. J. Phys. 44, 95 (2020)

2020
[40]

Z. G. Wang, Front. Phys. 21, 016300 (2026)

2026
[41]

Becchi, A

C. Becchi, A. Giachino, L. Maiani, and E. Santopinto, Phys. Lett. B 806, 135495 (2020)

2020
[42]

Becchi, A

C. Becchi, A. Giachino, L. Maiani, and E. Santopinto, Phys. Lett. B 811, 135952 (2020)

2020
[43]

S. S. Agaev, K. Azizi, B. Barsbay, and H. Sundu, Phys. Rev. D 109, 014006 (2024)

2024
[44]

Navas et al

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

2024
[45]

Lakhina, and E

O. Lakhina, and E. S. Swanson, Phys. Rev. D 74, 014012 (2006)

2006
[46]

S. S. Agaev, K. Azizi, and H. Sundu, Phys. Lett. B 856, 138886 (2024)

2024
[47]

Hatton et al

D. Hatton et al. [HPQCD], Phys. Rev. D 102, 054511 (2020)

2020
[48]

Lucha, D

W. Lucha, D. Melikhov, and S. Simula, EPJ Web Conf. 80, 00043 (2014)

2014
[49]

J. L. Rosner, S. Stone, and R. S. Van de Water, arXiv:1509.02220

work page Pith review arXiv
[50]

Lubicz, A

V. Lubicz, A. Melis, and S. Simula, PoS LAT- TICE2016, 291 (2017)

2017