arxiv: 2507.11874 · v2 · submitted 2025-07-16 · ✦ hep-ph

Recognition: 3 theorem links

· Lean Theorem

Fully strange tetraquark states via QCD sum rules

Bing-Dong Wan , Ji-Chong Yang

Authors on Pith no claims yet

Pith reviewed 2026-05-06 17:21 UTC · model claude-opus-4-7

classification ✦ hep-ph PACS 12.38.Lg12.39.Mk14.40.Rt

keywords fully strange tetraquarkQCD sum rulesexotic hadronsX(2300)molecular tetraquarkssss̄operator product expansionBESIII

0 comments

The pith

Fully-strange ssss̄ tetraquarks with masses 2.07–3.12 GeV emerge from a QCD sum-rule analysis, and the 1^+- state matches the BESIII X(2300).

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

The paper asks whether bound or resonant states made entirely of two strange quarks and two strange antiquarks should exist, and at what masses. Using QCD sum rules with molecular-type interpolating currents — pairs of strange quark-antiquark bilinears coupled through different Dirac structures — the authors compute masses for every spin-parity channel that admits such a current, finding a spectrum stretched between roughly 2.07 and 3.12 GeV. They argue that the 1^+- entry near 2.29 GeV lines up with the X(2300) recently seen by BESIII in ψ(3686) → φηη', supporting a tetraquark reading of that resonance. They also predict an exotic 0^-- state near 2.46 GeV, which cannot arise from a conventional quark-antiquark meson and would therefore be a clean signature of multiquark dynamics if found. For each current they identify the two-meson channels it most naturally decays into, giving experiments a concrete shopping list of final states to scan.

Core claim

The authors apply QCD sum rules to a complete set of meson–meson interpolating currents built from two strange quarks and two strange antiquarks, organised by the quantum numbers J^PC = 0^++, 0^-+, 0^--, 1^--, 1^+-, 1^++. For each current they extract a Borel window in which both the operator-product expansion converges and the ground-state pole dominates, and they read off masses ranging roughly from 2.07 to 3.12 GeV. The 1^+- channel built from a γ5 ⊗ γμ molecular current comes out at 2.29 ± 0.14 GeV, which they read as a tetraquark interpretation of the X(2300) resonance recently reported in ψ(3686) → φηη'. They also tabulate decay channels (φφ, ηη, φη, φη', K K-bar*, φηπ, ηηγ, …) tied to

What carries the argument

A set of fourteen meson–meson interpolating currents of the form [s̄ Γ1 s][s̄ Γ2 s], one per allowed Dirac structure for each J^PC, fed into the standard QCD sum-rule machinery: two-point correlator, operator-product expansion through dimension-8 condensates (quark condensate, gluon condensate, mixed quark-gluon condensate, and their products), Borel transform, and a continuum threshold √s0 chosen so that the extracted mass is stationary in the Borel parameter while the pole contribution exceeds the continuum. Each current's Dirac content (γ5⊗γ5, γμ⊗γμ, σμν⊗σμν, etc.) also dictates which two-meson decay channel it couples to.

If this is right

The X(2300) observed by BESIII in ψ(3686) → φηη' is consistent with a 1^+- fully-strange molecular tetraquark rather than a conventional strangeonium.
An exotic 0^-- fully-strange state near 2.46 GeV should exist; because 0^-- is forbidden for q q-bar mesons, finding it (e.g. in φηπ or ηηγ) would be a clean signal of multiquark dynamics.
Different Dirac structures of the molecular current preferentially decay into different two-meson channels (e.g. ηη, f0 f0, φφ, f1 f1 for the 0^++ family), giving a channel-by-channel template for sorting candidates.
A series of additional fully-strange resonances should populate φη, φη', K K-bar, and K* K-bar invariant-mass spectra between roughly 2.5 and 3.1 GeV, accessible at BESIII, Belle II, and LHCb.
The partial overlap with constituent-quark-model spectra for ssss̄ supports treating molecular and compact diquark pictures as complementary rather than mutually exclusive.

Where Pith is reading between the lines

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

Because each interpolating current's Dirac structure already encodes a preferred two-meson decay channel, the table of currents is effectively a falsifiable channel-mass map: a single experimental partial-wave fit could test several entries simultaneously rather than one resonance at a time.
The exotic 0^-- prediction near 2.46 GeV is the cleanest discriminator in the set, since any signal there cannot be faked by ordinary q q-bar mixing; this is where an experiment would get the most information per event.
The wide spread of 0^++ masses (2.23 to 3.00 GeV) across different currents suggests these are not five different states but five different projections onto the same multi-resonance scalar spectrum, and a coupled-channel treatment would likely collapse them into fewer physical poles.
The agreement with constituent-quark-model masses near 2.2–3.3 GeV is non-trivial because the two methods make different structural assumptions; this convergence makes the existence of *some* ssss̄ structure in this window more plausible than the identification of any specific state.

Load-bearing premise

That for each chosen current the lowest physical state it couples to is a genuine narrow resonance — not a smeared two-meson continuum or a mixture with ordinary strangeonium — and that the continuum threshold picked to stabilise the Borel plateau really separates that resonance from everything above it.

What would settle it

A partial-wave analysis at BESIII, Belle II, or LHCb that pins down the J^PC of X(2300) and measures its dominant decay channels: if X(2300) is 1^+- and decays primarily to φη / φη' / K K-bar*, the identification stands; if its quantum numbers are different, or if no fully-strange resonances appear in φφ, φη, or φηπ within the predicted 2.07–3.12 GeV window (especially the exotic 0^-- prediction near 2.46 GeV), the molecular tetraquark spectrum claimed here is wrong.

read the original abstract

In this paper, we have systematically explored the mass spectrum of fully strange tetraquark candidates within the framework of QCD sum rules, focusing on states with quantum numbers $J^{PC}=0^{++}$, $0^{-+}$, $0^{--}$, $1^{--}$, $1^{+-}$, and $1^{++}$. The analysis reveals the existence of fully strange tetraquark states with masses ranging from approximately $2.07$ to $3.12$ GeV. These predictions are confronted with existing experimental observations of potential fully strange tetraquark resonances, notably the $X(2300)$ recently reported by the BESIII Collaboration, which may be interpreted as a fully strange tetraquark state. Furthermore, the possible decay modes of these fully strange tetraquark states are analyzed, providing guidance for their identification in current and future high energy experiments such as BESIII, Belle II, and LHCb.

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

Competent but routine QCDSR survey of ss̄ss̄ molecular currents; the X(2300) match is suggestive at best, and the factorizable-diagram worry is real but not fatal.

read the letter

Quick read: this is a workmanlike QCD sum rules paper that catalogues 14 molecular [s̄Γs][s̄Γs] interpolating currents across six J^PC channels, runs the standard Shifman–Vainshtein–Zakharov machinery up to dimension-8 condensates, and reports masses in the 2.2–3.1 GeV range (Table I). Nothing here is methodologically new — it's the same pipeline Wan and collaborators have applied repeatedly — but the per-channel coverage of the fully strange sector with molecular currents is a useful reference table, and §IV's mapping from current structure to likely decay modes (ηη, φφ, φη, K*K̄, etc.) is genuinely helpful for experimentalists doing partial-wave analyses.

What works: the standard criteria (OPE convergence, >50% pole contribution, Borel-window stability) are checked channel-by-channel in the appendix figures, the input parameters are mainstream, and the discussion section now (after revision) honestly acknowledges that different currents probe different internal structures rather than pretending they all converge on one state. The added comparison with Liu et al. (Ref. [22]) and the GEM analysis of Ma et al. ([28]) is a real improvement.

Soft spots, in order of seriousness:

1. The stress-test concern about factorizable two-meson contributions is the most substantive issue. Every current is C⊗C color-singlet, and the disconnected piece — free propagation of two s̄Γs mesons — is not subtracted, separated, or even displayed. For j^A_{1+−} the natural disconnected channel sits near φ+η, and 2.29 GeV is exactly where you'd expect the factorizable continuum to peak. This is the Kondo–Morimatsu–Nishikawa / Lucha–Melikhov critique, and the paper does not address it. The X(2300) "agreement" should be read with that caveat.

2. Spectral densities are not written out, no code is shipped, and there is no direct numerical comparison with the earlier QCDSR ss̄ss̄ analyses by Chen and collaborators ([23,24,26,27]) or Wang ([29]) — which is what one actually wants to know to judge whether this adds anything beyond a relabeling.

3. The decay discussion is qualitative. "φη likely dominant due to phase space" is not a width estimate.

4. With ~150 MeV uncertainties spread over 14 channels, hitting 2300 MeV once is weak evidence for the X(2300) interpretation.

Recommendation: send it out to referees. It's a legitimate addition to the literature in a sector that does need systematic coverage, and the factorizable-diagram issue is something a referee can press the authors to address explicitly. Not reading-group material — too incremental — and I wouldn't cite it unless I were specifically working on ss̄ss̄ phenomenology.

Referee Report

6 major / 8 minor

Summary. The authors apply QCD sum rules to a set of molecular-type interpolating currents [s̄Γ₁s][s̄Γ₂s] to estimate masses of fully-strange tetraquark candidates with J^PC = 0^{++}, 0^{-+}, 0^{--}, 1^{--}, 1^{+-}, 1^{++}. Using the OPE up to dimension 8 (perturbative + ⟨s̄s⟩, ⟨G²⟩, ⟨s̄Gs⟩, ⟨s̄s⟩², mixed), Borel transformation, and the standard pole-dominance / OPE-convergence criteria to fix s₀ and the Borel window, they extract 14 mass values (Table I) spanning roughly 2.23–3.07 GeV. They compare with the X(2300) reported by BESIII and propose, in particular, that the J^PC = 1^{+-} A-type current (mass 2.29 ± 0.14 GeV) is a viable interpretation. A largely qualitative discussion of dominant decay channels and a comparison with potential-model results (Refs. [22,28]) are provided.

Significance. If the predictions hold, the paper provides a useful, organized survey of the molecular-type ss̄ss̄ sector across six J^PC channels, including the exotic 0^{--}, complementing the previous 1^{-+} sum-rule work of the authors and recent quark-model results. The 1^{+-} match to X(2300) is testable: the proposed dominant decay modes (φη, φη', KK̄*) overlap with BESIII reconstruction channels and would be falsifiable at BESIII, Belle II, and LHCb. The systematic enumeration of independent Lorentz/Dirac-structure currents per channel, with Borel-window diagnostics shown explicitly in Figs. 2–15, is a strength and supports reproducibility. The main correctness-risk concerns are technical and well-known in the molecular-QCDSR literature; they are addressable rather than fatal, but they do bear on the X(2300) identification.

major comments (6)

[Eqs. (1)–(14), (15)–(20)] Factorizable two-meson contamination of the molecular currents is not addressed. Every current is a color-singlet × color-singlet product [s̄Γ₁s][s̄Γ₂s]; in the two-point function (Eqs. 15–16) the leading perturbative loop and the lower-dimensional condensate contributions decompose into a disconnected piece that is the product of two ordinary qq̄ correlators (free two-meson propagation) plus a connected/non-factorizable piece. As emphasized by Kondo–Morimatsu–Nishikawa and by Lucha–Melikhov–Sazdjian, the factorizable piece does not couple to a genuine tetraquark pole and tends to saturate the dispersion integral with two-meson continuum strength near the relevant threshold. The manuscript does not show the spectral density ρ^OPE(s) explicitly, does not separate connected from disconnected diagrams, and does not demonstrate that the non-factorizable contribution is comparable to the fact
[§III, Table I; §V (X(2300) interpretation)] The X(2300) identification is load-bearing for the abstract and conclusions, but the evidence as presented is weak. (a) j^A_{1^{+-}} couples to η_s ⊗ φ-type two-meson states; the φη' threshold (~1.98 GeV) and the φη(1295) threshold lie within roughly 1σ of the quoted 2.29 ± 0.14 GeV, exactly where factorizable contributions dominate. (b) The 14 currents span 2.23–3.07 GeV; given this density, picking the channel/current closest to a target observed mass is not a strong test. (c) Ref. [72], cited by the authors, prefers an s̄sqq̄ assignment for X(2300); the manuscript notes this in passing but does not quantitatively confront it. A more discriminating analysis (e.g., predicted partial widths into φη vs φη', or a finite-energy sum rule with explicit two-meson subtraction) is needed to make the 1^{+-} assignment more than coincidence.
[§IV.A; Table I, 0^{++} rows] The interpretation in §IV.A — that the five 0^{++} currents (Eqs. 1–5) correspond to five physically distinct tetraquark states with masses 2.23, 2.24, 2.50, 2.66, 3.00 GeV — is not justified within the QCDSR framework as applied here. For a fixed J^PC, different currents are not orthogonal: each in general couples to all states with that J^PC, with different overlaps. The standard procedure is either (i) to identify the lowest state with the ground state (and treat the others as projection artefacts/excited states with a duality cut) or (ii) to diagonalize the matrix of currents. The current text reads instead as if each current selects 'its own' resonance, with the J^E_{0^{++}} mass driven up to 3.00 GeV explained by 'coupling to two 2^{--} mesons' — but σ_{μν}⊗σ_{μν} does not select J=2 in a meaningful sense. Please clarify the physical content: are these claimed to be distinct narrow
[Figs. 2–15; Eqs. (26)–(27)] Pole dominance is invoked as R^PC > 0.5, but inspection of the figures suggests this is marginal at the upper end of several quoted Borel windows. For example, Fig. 2(a): at M_B² ≈ 2.2 GeV² with √s₀ = 2.9 GeV the red R^PC curve lies near or below 0.5; similar behaviour appears in Figs. 5, 7, 8, 12, 13. Please provide tabulated values of R^PC (and R^OPE) at the lower and upper boundaries of every quoted Borel window so the reader can verify that both criteria are jointly satisfied across the entire window for each current. If the plateau region with R^PC > 0.5 is narrower than reported, the mass uncertainties in Table I are underestimated.
[Eq. (20); Eq. (21)] OPE convergence is asserted but not demonstrated. With the propagator (21) and four-quark operators, contributions of dimension d ≥ 8 (e.g., ⟨s̄s⟩⟨s̄Gs⟩, ⟨G²⟩⟨s̄s⟩, and the d=10 four-quark/gluon mixed operators) are known to be sizeable for ss̄ss̄ systems, and ⟨s̄s⟩² already enters at d=6 with a strange-quark enhancement. Please show the partial OPE contributions order-by-order (or as a ratio plot) so that R^OPE = |L_{d=8}/L_total| < some threshold is genuinely a measure of convergence and not just of the highest term retained. As a minimum, a discussion of why d=10 contributions are expected to be subleading would strengthen the analysis.
[§IV.B–§IV.C; Table II] The decay-mode discussion is qualitative and in places contains errors that should be corrected: (i) §IV.C invokes 'CKM suppressed' couplings for fully-strange tetraquark strong decays, where CKM has no role; (ii) the listed channels mix S-wave/P-wave assignments without explicit angular-momentum bookkeeping (e.g., 0^{++} → φφ is S-wave only above 2 m_φ ≈ 2.04 GeV and even then has nontrivial polarization structure). Since these decay tables are the experimental hook of the paper, even a leading-order quark-pair-creation (3P0) or simple coupling-constant estimate of the ratios ϕη : ϕη' : KK̄* would substantially raise the predictive content.

minor comments (8)

[Abstract] The stated range '2.07 to 3.12 GeV' is the envelope obtained by adding/subtracting the quoted uncertainties to the extreme entries of Table I; it would be clearer to state 'central values 2.23–3.07 GeV, with uncertainties extending the envelope to 2.07–3.12 GeV'.
[§I, paragraph after Eq. introduction of currents] The statement that 'the 0^{+-} state does not admit a molecular configuration' deserves a one-line justification (C-parity of γ_μ × γ_ν combinations) rather than being asserted.
[§II, Eqs. (1)–(14)] Several typos in subscripts/superscripts of currents (e.g., 'j^{k,μ}_{JJC}' in Eq. 16, 'jD_{0^{++}}' uses γ_μγ_5 ⊗ γ^μγ_5 — please clarify the index contraction explicitly). Please proofread.
[§IV.A] The phrase 'j^E_{0^{++}} couples to two 2^{--} mesons' is misleading. σ_{μν} interpolates 1^{+-}/1^{++} mesons depending on parity, not 2^{--}. Please rephrase.
[References] Refs. [41]–[69] are heavily author-self-citations of QCDSR applications. While not improper, the introduction would benefit from a few additional independent references on factorization issues in molecular QCDSR (e.g., Kondo–Morimatsu–Nishikawa; Lucha–Melikhov–Sazdjian) given that this is the dominant methodological concern for this paper.
[§V, paragraph on Ref. [72]] The note that Ref. [72] favors an s̄sqq̄ structure for X(2300) is added almost as an aside. Given that this directly contradicts the paper's headline interpretation, a more substantive discussion (one short paragraph) of how the proposed decay-channel analysis would distinguish the two assignments would be valuable.
[Figs. 2–15] The figures are numerous and nearly identical in form. Consider consolidating into a multi-panel figure or moving some to a supplemental file; alternatively, ensure axis labels and curve labels are uniform (e.g., R^PC_{A,0++} vs R^A_PC notation appears in different figures).
[Table I, λ values] The decay constants λ are quoted but never used downstream. Either employ them (e.g., to estimate production cross sections at BESIII) or remove to streamline.

Simulated Author's Rebuttal

6 responses · 1 unresolved

We thank the referee for a careful and technically substantive report. The major comments fall into three groups: (i) a foundational issue concerning factorizable two-meson contamination of color-singlet molecular currents, which directly bears on the X(2300) interpretation; (ii) presentation/diagnostic deficiencies (Borel-window endpoints, OPE convergence shown order-by-order, multi-current interpretation in §IV.A); and (iii) errors in the decay-mode section (a stray 'CKM suppressed' reference, missing partial-wave bookkeeping). We accept essentially all of the referee's technical points. In revision we will (a) explicitly separate connected from disconnected contributions to every correlator and condition the ground-state extraction on connected-piece dominance, (b) soften the X(2300) identification to a consistency statement and add a φη/φη'/KK̄* branching-ratio estimate as a falsifiable test, (c) reinterpret multiple currents per J^PC as overlapping projections rather than distinct states and adopt the lowest extracted value as the ground-state estimate, (d) tabulate R^OPE and R^PC at every Borel-window endpoint and tighten windows where R^PC<0.5, (e) display the OPE term-by-term and add d=10 as a systematic, and (f) correct the decay-mode discussion. One limitation remains: a full current-matrix diagonalization and finite-energy sum rule with explicit two-meson subtraction are beyond this revision and will be flagged as such.

read point-by-point responses

Referee: Factorizable two-meson contamination of molecular currents: every current is color-singlet × color-singlet, so the two-point function decomposes into a disconnected (free two-meson) piece plus a connected/non-factorizable piece (Kondo–Morimatsu–Nishikawa; Lucha–Melikhov–Sazdjian). The manuscript does not show ρ^OPE(s) explicitly, does not separate connected from disconnected diagrams, and does not demonstrate that the non-factorizable contribution is comparable to the factorizable one.

Authors: We agree this is the central correctness issue for molecular-type currents and we did not address it adequately. In the revision we will (i) explicitly decompose each correlator into disconnected (factorizable) and connected (non-factorizable) Wick contractions at the diagrammatic level, displaying the relevant Feynman graphs; (ii) tabulate the relative size of L_0^conn vs. L_0^disc inside each Borel window; and (iii) restrict the quoted ground-state masses to the cases where the connected piece is not parametrically subdominant. For currents where the disconnected part dominates near the relevant two-meson threshold, we will explicitly flag the extracted value as compatible with two-meson continuum strength rather than a genuine tetraquark pole, following the prescription of Lucha–Melikhov–Sazdjian. We will also add the explicit ρ^OPE(s) expressions (currently relegated to internal notes) as an appendix so readers can reproduce the connected/disconnected split. We thank the referee for forcing this clarification. revision: yes
Referee: The X(2300) identification is load-bearing but weakly supported: the j^A_{1^{+-}} two-meson thresholds (φη', φη(1295)) sit within 1σ of 2.29 ± 0.14 GeV where factorizable contributions dominate; with 14 currents spanning 2.23–3.07 GeV, picking the closest one is not a strong test; Ref. [72] prefers s̄sqq̄. A discriminating analysis (partial widths φη vs φη', or finite-energy sum rule with explicit two-meson subtraction) is needed.

Authors: The referee's criticism is fair. In the revised version we will: (a) soften the abstract and conclusions to state only that the 1^{+-} A-current mass is *consistent* with X(2300), not that it identifies the state; (b) add the φη' and φη(1295) thresholds explicitly to the discussion and acknowledge that the central value lies in the factorization-dominated region, conditioning the interpretation on the connected-piece analysis described in our reply to comment 1; (c) include a quantitative comparison with Ref. [72] noting that the s̄sqq̄ assignment cannot be excluded by our analysis, since QCDSR with [s̄s][s̄s] currents is by construction blind to that flavor structure; and (d) compute, at the level of QCDSR three-point coupling estimates, the ratios Γ(φη)/Γ(φη') for the 1^{+-} A current as a falsifiable discriminator. A full finite-energy sum rule with explicit two-meson subtraction is a substantial separate project; we will mention it as the next step rather than claim it here. revision: yes
Referee: Interpretation in §IV.A that the five 0^{++} currents correspond to five physically distinct states is not justified: different currents are not orthogonal, each couples to all J^PC states. Standard procedure is to identify the lowest as ground state or diagonalize the current matrix. The σ_{μν}⊗σ_{μν} → 'two 2^{--} mesons' explanation for j^E is not meaningful.

Authors: The referee is correct on both points. We will rewrite §IV.A. The five 0^{++} currents are not orthogonal projectors onto five distinct resonances; each in principle couples to the full tower of 0^{++} states with different overlaps λ_k. The mass extracted from a given current is best interpreted as a weighted average dominated by whichever physical state has the largest overlap with that interpolating field. In the revision we will (i) explicitly state this caveat; (ii) identify the *lowest* extracted mass per J^PC as our ground-state estimate, treating the higher-current values as upper bounds / signals of excited or mixed components rather than independent states; (iii) remove the incorrect statement that σ_{μν}⊗σ_{μν} 'couples to two 2^{--} mesons' — the σ_{μν} bilinear contains both 1^{+-} and 1^{--} components and the J=2 wording was loose. A proper diagonalization of the 5×5 current matrix is in principle possible and we will outline it, though carrying it out fully is beyond the present scope. revision: yes
Referee: Pole dominance R^PC > 0.5 appears marginal at the upper end of several Borel windows (Figs. 2, 5, 7, 8, 12, 13). Please tabulate R^PC and R^OPE at the lower and upper boundaries of every quoted window. If the plateau is narrower than reported, mass uncertainties in Table I are underestimated.

Authors: We will add a new table giving R^OPE and R^PC at both endpoints of every Borel window for all 14 currents, plus the value at the central M_B². On re-examination several windows do indeed have R^PC close to 0.5 at the upper edge; for those cases we will tighten the upper bound of M_B² to enforce R^PC ≥ 0.5 strictly across the entire window and re-extract the masses. Where this narrows the window appreciably, we will inflate the quoted mass uncertainty accordingly. We expect the central values to shift by less than the existing error bars, but the error bars themselves will increase for a subset of currents (most likely the 0^{++} E and 0^{-+} B channels, where the OPE is most sensitive). We thank the referee for catching this. revision: yes
Referee: OPE convergence is asserted, not demonstrated. ⟨s̄s⟩² already enters at d=6 with strange-quark enhancement; d=8 mixed condensates and d=10 contributions can be sizeable for ss̄ss̄. R^OPE as defined only measures the highest term retained. Please show partial contributions order-by-order, or at minimum justify why d=10 is subleading.

Authors: We accept this. In the revised manuscript we will provide (as a figure and a supplementary table) the partial contributions of each OPE term — perturbative, ⟨s̄s⟩, ⟨G²⟩, ⟨s̄Gs⟩, ⟨s̄s⟩², ⟨G²⟩⟨s̄s⟩, ⟨s̄s⟩⟨s̄Gs⟩ — as a function of M_B² for each current, normalized to the total. This makes hierarchical convergence visible rather than asserted. Regarding d=10: a fully consistent inclusion of d=10 four-quark/gluon mixed condensates is not standard in the molecular-tetraquark QCDSR literature and would require additional vacuum-saturation assumptions. We will discuss why the leading d=10 operators are expected to be suppressed by additional factors of m_s² or by extra loop counting, and we will treat the omitted d=10 contribution as a systematic uncertainty added in quadrature to the mass error. We agree the previous R^OPE definition was a weak diagnostic and will replace it. revision: yes
Referee: Decay-mode discussion contains errors: (i) §IV.C invokes 'CKM suppressed' couplings for fully-strange strong decays where CKM plays no role; (ii) S-wave/P-wave assignments are mixed without angular-momentum bookkeeping (e.g., 0^{++}→φφ requires 2m_φ threshold and has polarization structure). A leading-order ³P_0 or coupling-constant estimate of φη:φη':KK̄* ratios would help.

Authors: Both errors are correct and will be fixed. (i) The phrase 'CKM suppressed' is a straightforward mistake — strong decays of an ss̄ss̄ state involve no flavor-changing currents — and will be removed. (ii) We will rewrite Table II and the accompanying text to make the partial-wave assignment for each channel explicit (S/P/D), include phase-space factors |p|^(2L+1), and flag channels that are kinematically closed for the lower-mass currents (φφ for any state below ~2.04 GeV; f_1 f_1 for everything in our spectrum). For the 1^{+-}, 1^{--}, and 1^{++} states we will provide leading-order ³P_0-model estimates of the ratios Γ(φη):Γ(φη'):Γ(KK̄*), which directly support the X(2300) discussion in comment 2. We will not promote these to absolute partial widths, since the ³P_0 strength parameter for a tetraquark source is poorly constrained, but the ratios are robust enough to be experimentally useful. revision: yes

standing simulated objections not resolved

A complete diagonalization of the full current-mixing matrix for each J^PC channel, and a full finite-energy sum rule with explicit two-meson continuum subtraction, are acknowledged as the methodologically correct next step but will not be carried out in the present revision; we will state this explicitly as a limitation rather than claim to have done it.

Circularity Check

1 steps flagged

Standard QCD sum rules analysis with the usual self-consistent s₀ selection; no central claim reduces to its own inputs by definition.

specific steps

fitted input called prediction [Sec. III, Eqs. (26)–(27) and Table I]
"the goal is to identify an optimal value of s₀ that yields a stable Borel window for the extracted mass of the fully strange tetraquark state. Within this window, the mass prediction should exhibit minimal dependence on the Borel parameter M²_B... In practical calculations, √s₀ is varied by 0.1 GeV around a central value"

The continuum threshold s₀ is selected by requiring a stable plateau for the very mass being predicted; across Table I √s₀ tracks M+0.5–0.7 GeV. This is the conventional QCDSR self-consistency criterion, so the 'prediction' is partially shaped by the answer. However, this is the universal QCDSR convention, not a step unique to this paper, and the masses still confront external thresholds and experimental data, so the circularity is mild.

full rationale

This paper applies a well-established QCDSR methodology to fully-strange tetraquark interpolating currents. The derivation chain is: (i) construct molecular currents in Eqs. (1)–(14); (ii) compute the OPE-side correlation function with quark propagator (21) and condensates; (iii) match to a phenomenological pole+continuum form (18); (iv) Borel transform and extract mass via Eq. (22). The inputs (quark masses, condensates ⟨s̄s⟩, ⟨G²⟩, ⟨s̄Gσ·Gs⟩, m₀²) are taken from external references [34–39], and the resulting masses are not algebraically equivalent to any single input. The principal methodological loop is the standard QCDSR self-consistency: the continuum threshold √s₀ is chosen by requiring a 'stable Borel window' for the very mass being extracted (Sec. III, Eqs. 26–27), and in practice √s₀ ≈ M + 0.5–0.7 GeV across Table I. This is a known soft circularity in all QCDSR work — s₀ is informed by the answer — but it is the universally accepted convention rather than a step unique to this paper, and the resulting masses can in principle be confronted with external data (e.g., X(2300)). The reader's skeptic attack about factorizable [s̄Γs][s̄Γs] disconnected diagrams dominating the OPE is a serious correctness/methodology concern (the paper indeed never separates connected from disconnected pieces, never displays the spectral density, and the agreement with two-meson thresholds near 2 m_φ, m_φ+m_η' is suggestive), but this is a *validity* problem, not a circular derivation: the extracted number is not by construction equal to a fitted input. Self-citations (Ref. [31] for 1⁻⁺ states, and the long list [42–69] of prior QCDSR work by the same author group) are conventional methodological references, not load-bearing uniqueness theorems being smuggled in. The interpretive claim that the 1⁺⁻ state at 2.29 GeV 'is' the X(2300) is post-hoc identification, not circular derivation. Overall: standard methodology, ordinary self-consistency, score 2.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The central claim rests on standard QCDSR machinery: SVZ duality, OPE truncation, vacuum condensate values from external sum-rule literature, and Borel-window selection. There are no invented entities. The "free parameters" are the per-current continuum threshold √s₀ (chosen for plateau stability) and the Borel window M²_B; vacuum condensates and m₀² are inputs from prior QCDSR phenomenology rather than this paper's fits. The molecular-current ansatz itself is a domain assumption that excludes diquark-antidiquark and color-adjoint structures, as the authors note.

pith-pipeline@v0.9.0 · 9498 in / 6560 out tokens · 97963 ms · 2026-05-06T17:21:02.679671+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

Unification.YangMillsMassGap spectral_gap / mass_gap_from_forcing_chain unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the masses of the fully strange tetraquark states in molecular configurations with J^PC = 0++, 0−+, 0−−, 1−−, 1+−, and 1++ are investigated within the framework of QCD sum rules
Cost.FunctionalEquation washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

ρ^OPE(s) = ρ^pert(s) + ρ⟨s̄s⟩(s) + ρ⟨G²⟩(s) + ρ⟨s̄Gs⟩(s) + ρ⟨s̄s⟩²(s) + ...
Foundation.ConstantDerivations all_constants_from_phi unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

input parameters ... m_s = (95±5) MeV, ⟨s̄s⟩ = (0.8±0.1)⟨q̄q⟩, ⟨g_s² G²⟩ = (0.88±0.25) GeV^4, m₀² = (0.8±0.2) GeV²; M = 2.23 ± 0.15 GeV ...

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.

Forward citations

Cited by 1 Pith paper

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

QCD sum rule analysis of local meson-meson currents for the $K(1690)$ state
hep-ph 2026-04 unverdicted novelty 4.0

QCD sum rules with local meson-meson currents for the K(1690) consistently predict masses around 2 GeV or above, disfavoring a molecular interpretation in favor of a compact multiquark state.