Pith. sign in

REVIEW 2 major objections 5 minor 61 references

In a constituent-quark model, compact fully light hexaquarks sit well above baryon thresholds, while dibaryon-like states are near-threshold and molecular; one is deuteron-like but unbound by about 12 MeV.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-13 01:19 UTC pith:FDPTRHE4

load-bearing objection Clean DMC survey of light hexaquarks: compact states sit well above thresholds, dibaryon-like ones are molecular and near-threshold, and the deuteron-like (1,0) channel remains unbound by ~12 MeV inside AL1. the 2 major comments →

arxiv 2607.09288 v1 pith:FDPTRHE4 submitted 2026-07-10 hep-ph hep-exhep-latnucl-exnucl-th

Diffusion Monte Carlo study of deuteron-like fully light hexaquarks

classification hep-ph hep-exhep-latnucl-exnucl-th
keywords hexaquarksdibaryonsdiffusion Monte Carloconstituent quark modeldeuteron-like statesmultiquark spectroscopyAL1 potential
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper asks whether six light quarks (three up and three down) can form stable or resonant hexaquarks, either as compact six-quark bags or as loosely bound baryon pairs. Using diffusion Monte Carlo to solve the six-body Schrödinger equation with a standard constituent-quark potential, the authors find that every fully antisymmetric compact configuration lies hundreds of MeV above the relevant two-baryon thresholds, so deeply bound compact states are not favored. By contrast, several dibaryon-like arrangements, which allow hidden-color components, sit close to but above the NN, NΔ and ΔΔ thresholds and display clear molecular spatial structure. The (S,I)=(1,0) state in particular shows two nucleon-like clusters separated by several femtometers, closely resembling the deuteron, yet its mass remains slightly unbound relative to the model NN threshold. The result matters because it supplies a clean, numerically controlled prediction for the light-hexaquark sector and indicates that any real deuteron-like state, if it exists, is likely molecular rather than compact.

Core claim

Within the AL1 constituent-quark model solved by diffusion Monte Carlo, all symmetry-allowed compact uudddd hexaquarks lie well above their baryon-baryon thresholds (lowest compact mass 2716±7 MeV versus 2MΔ=2614 MeV), while several dibaryon-like configurations lie close but above the NN, NΔ and ΔΔ thresholds and exhibit molecular radial distributions; the (S,I)=(1,0) state is spatially deuteron-like yet unbound by 12±9 MeV relative to the model NN threshold.

What carries the argument

Diffusion Monte Carlo projection of the six-body ground state with a trial wave function that is either fully antisymmetric (compact) or only intra-cluster antisymmetric (dibaryon-like, allowing hidden color), combined with radial distribution functions that distinguish compact from molecular geometries.

Load-bearing premise

The non-relativistic quark-quark potential previously fitted to ordinary mesons and baryons remains quantitatively reliable for six light quarks near two-baryon thresholds.

What would settle it

A lattice-QCD or alternative-potential calculation of the same (S,I)=(1,0) uudddd system that finds a mass below the model NN threshold, or a high-precision experimental search that rules out any near-threshold molecular hexaquark signal in the deuteron channel.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The manuscript applies the diffusion Monte Carlo method to the non-relativistic six-body Hamiltonian of Eq. (1) with the AL1 constituent-quark potential for fully light (uuuddd) hexaquarks. Two classes of antisymmetric trial wave functions are used: fully antisymmetric compact configurations and partially antisymmetrized dibaryon-like arrangements that can include hidden-color components. Table I reports the resulting energies for all symmetry-allowed (S,I) channels. Compact states lie well above the model baryon–baryon thresholds (lowest compact mass 2716 ± 7 MeV versus 2M_Δ = 2614 MeV). Several dibaryon-like states lie close but above the NN, NΔ and ΔΔ thresholds and exhibit molecular radial distributions (Figs. 1–2). The (S,I)=(1,0) configuration shows two nucleon-like subclusters separated by several fm, yet remains unbound by 12 ± 9 MeV relative to the model NN threshold of 2060 MeV.

Significance. If the AL1 results hold, the work supplies a clean, internally consistent demonstration that deeply bound compact six-quark states are disfavored in the light sector while near-threshold molecular configurations can appear. The explicit construction of both compact and dibaryon-like trial functions, the DMC projection, and the direct comparison of radial distributions to free nucleon and Δ profiles constitute a technically solid contribution to multiquark spectroscopy. The deuteron-like (S,I)=(1,0) state is a particularly useful benchmark: its spatial structure is recovered while the small binding is not, furnishing a quantitative test of the model near the two-nucleon threshold. The calculation is fully specified (Hamiltonian, potential parameters fixed by earlier hadron spectroscopy, statistical DMC errors) and therefore falsifiable within the same framework.

major comments (2)
  1. The central quantitative claim that the (S,I)=(1,0) state is unbound by only 12 ± 9 MeV rests on the AL1 potential (Eq. (1) and Refs. [53,54]) remaining accurate for six light quarks near the NN threshold. Because the same Hamiltonian generates both the hexaquark energy and the model thresholds 2M_N and 2M_Δ, the hierarchy is internally consistent, yet the absolute scale of the residual binding is model-dependent. A short discussion of how the known under-binding of the deuteron in AL1 (or a comparison with an alternative potential) would affect the 12 MeV figure would strengthen the interpretation without altering the numerical results.
  2. Table I lists several dibaryon-like channels as “· · ·” (infinite separation). The text states that DMC can split colorless clusters, but it is not shown whether the trial-function form of Eq. (4) or the partial-antisymmetrization procedure artificially suppresses or enhances dissociation in the borderline cases (e.g., (2,3), (3,2)). A brief numerical check that the same algorithm recovers free two-baryon energies when the clusters are forced far apart would remove residual doubt about the near-threshold assignments.
minor comments (5)
  1. Abstract and introduction: “fully light hexaquark containing three u quarks and three d quarks” is slightly awkward; “fully light hexaquarks composed of three u and three d quarks” is clearer.
  2. Eq. (3): the operator A is written with N_p and P_α; a short remark that the sum runs only over the allowed permutations for each ansatz (full versus partial) would help readers reconstruct the two sectors.
  3. Fig. 2 caption: the four panels are labeled by (S,I) but the text refers to (3,0,1^+); the J^P assignment should be made consistent.
  4. Reference [55] is cited as “M. C. Gordillo, (2026), arXiv:2604.18174”; if this is a companion paper still under review, a brief note on the status would be useful.
  5. Typographical: “baron-baryon” (p. 2) should be “baryon-baryon”; “Schrodinger” should be “Schrödinger” throughout.

Circularity Check

0 steps flagged

No significant circularity: hexaquark energies and model thresholds are independent outputs of the same fixed AL1 Hamiltonian solved by DMC; self-citations supply only the numerical method and prior baryon masses.

full rationale

The derivation chain is: (i) non-relativistic Hamiltonian Eq. (1) with the external AL1 potential of Refs. [53,54] (parameters fixed long ago to ordinary meson/baryon spectra, not to any hexaquark observable); (ii) two classes of antisymmetric trial functions (compact vs. dibaryon-like) constructed via the projector A of Eq. (3) and the standard Jastrow radial form Eq. (4); (iii) DMC projection of the ground-state energy for each (S,I) channel, yielding the numbers in Table I; (iv) comparison of those energies to the model two-baryon thresholds 2MN = 2060 MeV and 2MΔ = 2614 MeV that are themselves eigenvalues of the identical Hamiltonian on three-quark systems. Because both sides of every binding comparison are generated from the same fixed operator, the statements “compact states lie well above thresholds” and “the (1,0) configuration is unbound by 12 ± 9 MeV” are internally consistent numerical results, not tautologies or fitted predictions. Self-citations ([52], [55], [57], [60]) merely document the DMC implementation and the prior three-quark spectrum; they introduce no uniqueness theorem, no ansatz that forces the hexaquark spectrum, and no parameter re-tuned to the six-quark data. Radial distributions (Figs. 1–2) are post-processed observables, not inputs. Hence the paper contains no circular step of any of the six enumerated kinds.

Axiom & Free-Parameter Ledger

2 free parameters · 4 axioms · 0 invented entities

The central claim rests on a fixed non-relativistic Hamiltonian whose two-body potential was fitted elsewhere, on a particular choice of trial wave function and partial antisymmetrization that defines the two families of states, and on the assumption that DMC with that trial function yields the physical ground-state energy to within the quoted statistical errors. No new particles or forces are introduced; the free parameters are inherited from the AL1 model and the cusp-condition coefficients of the trial function.

free parameters (2)
  • AL1 potential parameters (quark masses, Coulomb strength, confinement and hyperfine couplings) = literature values of Refs. [53,54]
    Taken from Semay & Silvestre-Brac (1994) and Silvestre-Brac (1996); originally fitted to ordinary meson and baryon spectra. All absolute masses and thresholds in the paper inherit these values.
  • trial-wave-function cusp coefficients a_ij
    Chosen by hand to satisfy the two-body cusp conditions of the Coulomb-like part of AL1 for each quark pair (Eq. 4). They affect the statistical efficiency and, potentially, residual bias of the DMC projection.
axioms (4)
  • domain assumption The six-quark system is adequately described by the non-relativistic Hamiltonian of Eq. (1) with only two-body AL1 interactions.
    Stated at the opening of the methods section; relativistic, three-body, and continuum-coupling effects are omitted without quantitative estimate.
  • domain assumption A totally symmetric radial trial function times an antisymmetric spin–isospin–color function (full or partial) is sufficient for DMC to project the physical ground state of each (S,I) channel.
    Eqs. (2)–(4) and the paragraph defining the antisymmetrizer A; no multi-determinant or explicitly correlated alternative is tested.
  • ad hoc to paper Partial antisymmetrization (first three quarks among themselves, last three among themselves) correctly isolates dibaryon-like configurations that may contain hidden color.
    Introduced in the paragraph after Eq. (3) as the operational definition of the dibaryon-like sector; the mapping to physical molecular states is assumed rather than derived from a cluster decomposition theorem.
  • domain assumption Statistical DMC error bars fully capture the uncertainty relevant to the ~12 MeV comparison with the NN threshold.
    Only Monte-Carlo statistical uncertainties are quoted; time-step, population-control, and fixed-node (or fixed-phase) biases are not separately quantified.

pith-pipeline@v1.1.0-grok45 · 13895 in / 3334 out tokens · 38644 ms · 2026-07-13T01:19:54.129141+00:00 · methodology

0 comments
read the original abstract

We perform a Diffusion Monte Carlo study of fully light hexaquark containing three $u$ quarks and three $d$ quarks within a constituent-quark model. Both compact and baryon--baryon-like arrangements were considered separately. All compact hexaquark configurations are found well above their corresponding baryon--baryon thresholds, suggesting that deeply bound compact six-quark states are not favored in the light-quark sector within this model. By contrast, several dibaryon-like configurations lie close, but above, to the $NN$, $N\Delta$, and $\Delta\Delta$ thresholds and show spatial structures compatible with molecular states. One configuration exhibits two well-defined nucleon-like subclusters separated by several femtometers, closely resembling the spatial structure of the deuteron, although its calculated energy remains slightly above the corresponding threshold.

Figures

Figures reproduced from arXiv: 2607.09288 by J. Segovia, M.C. Gordillo.

Figure 1
Figure 1. Figure 1: FIG. 1. Radial distribution functions [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

61 extracted references · 41 linked inside Pith

  1. [1]

    Gell-Mann, Phys

    M. Gell-Mann, Phys. Lett.8, 214 (1964)

  2. [2]

    Zweig, CERN Report No.8182/TH.401, CERN Re- port No.8419/TH.412 (1964)

    G. Zweig, CERN Report No.8182/TH.401, CERN Re- port No.8419/TH.412 (1964)

  3. [3]

    Dong, F.-K

    X.-K. Dong, F.-K. Guo, and B.-S. Zou, Phys. Rev. Lett. 126, 152001 (2021), arXiv:2011.14517 [hep-ph]

  4. [4]

    H.-X. Chen, W. Chen, X. Liu, and S.-L. Zhu, Phys. Rept.639, 1 (2016), arXiv:1601.02092 [hep-ph]

  5. [5]

    H.-X. Chen, W. Chen, X. Liu, Y.-R. Liu, and S.-L. Zhu, Rept. Prog. Phys.80, 076201 (2017), arXiv:1609.08928 [hep-ph]

  6. [6]

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

  7. [7]

    Liu, H.-X

    Y.-R. Liu, H.-X. Chen, W. Chen, X. Liu, and S.- L. Zhu, Prog. Part. Nucl. Phys.107, 237 (2019), arXiv:1903.11976 [hep-ph]. 5

  8. [8]

    G. Yang, J. Ping, and J. Segovia, Symmetry12, 1869 (2020), arXiv:2009.00238 [hep-ph]

  9. [9]

    Dong, F.-K

    X.-K. Dong, F.-K. Guo, and B.-S. Zou, Commun. Theor. Phys.73, 125201 (2021), arXiv:2108.02673 [hep-ph]

  10. [10]

    Chen, Phys

    H.-X. Chen, Phys. Rev. D105, 094003 (2022), arXiv:2103.08586 [hep-ph]

  11. [11]

    Cao, (2023), arXiv:2301.11253 [hep-ph]

    X. Cao, (2023), arXiv:2301.11253 [hep-ph]

  12. [12]

    Mai, U.-G

    M. Mai, U.-G. Meißner, and C. Urbach, Phys. Rept. 1001, 1 (2023), arXiv:2206.01477 [hep-ph]

  13. [13]

    L. Meng, B. Wang, G.-J. Wang, and S.-L. Zhu, (2022), arXiv:2204.08716 [hep-ph]

  14. [14]

    H.-X. Chen, W. Chen, X. Liu, Y.-R. Liu, and S.-L. Zhu, Rept. Prog. Phys.86, 026201 (2023), arXiv:2204.02649 [hep-ph]

  15. [15]

    F.-K. Guo, H. Peng, J.-J. Xie, and X. Zhou, (2022), arXiv:2203.07141 [hep-ph]

  16. [16]

    P. G. Ortega and D. R. Entem, Symmetry13, 279 (2021), arXiv:2012.10105 [hep-ph]

  17. [17]

    Huang, C

    H. Huang, C. Deng, X. Liu, Y. Tan, and J. Ping, Sym- metry15, 1298 (2023)

  18. [18]

    R. F. Lebed, PoSFPCP2023, 028 (2023), arXiv:2308.00781 [hep-ph]

  19. [19]

    Zou, Sci

    B.-S. Zou, Sci. Bull.66, 1258 (2021), arXiv:2103.15273 [hep-ph]

  20. [20]

    M.-L. Du, V. Baru, F.-K. Guo, C. Hanhart, U.-G. Meißner, J. A. Oller, and Q. Wang, JHEP08, 157 (2021), arXiv:2102.07159 [hep-ph]

  21. [21]

    Liu, Y.-W

    M.-Z. Liu, Y.-W. Pan, Z.-W. Liu, T.-W. Wu, J.-X. Lu, and L.-S. Geng, (2024), arXiv:2404.06399 [hep-ph]

  22. [22]

    Johnson, I

    D. Johnson, I. Polyakov, T. Skwarnicki, and M. Wang, (2024), 10.1146/annurev-nucl-102422-040628, arXiv:2403.04051 [hep-ex]

  23. [23]

    D. R. Entem, F. Fern´ andez, P. G. Ortega, and J. Segovia, (2025), arXiv:2504.07897 [hep-ph]

  24. [24]

    Hanhart, (2025), arXiv:2504.06043 [hep-ph]

    C. Hanhart, (2025), arXiv:2504.06043 [hep-ph]

  25. [25]

    X. Wang, X. Liu, and Y. Gao, (2025), arXiv:2502.15117 [hep-ex]

  26. [26]

    Wang, (2025), arXiv:2502.11351 [hep-ph]

    Z.-G. Wang, (2025), arXiv:2502.11351 [hep-ph]

  27. [27]

    Francis, Prog

    A. Francis, Prog. Part. Nucl. Phys.140, 104143 (2025), arXiv:2502.04701 [hep-lat]

  28. [28]

    J.-H. Chen, J. Chen, F.-K. Guo, Y.-G. Ma, C.-P. Shen, Q.-Y. Shou, Q. Shou, Q. Wang, J.-J. Wu, and B.-S. Zou, Nucl. Sci. Tech.36, 55 (2025), arXiv:2411.18257 [hep-ph]

  29. [29]

    H¨ usken, E

    N. H¨ usken, E. S. Norella, and I. Polyakov, Mod. Phys. Lett. A40, 2530002 (2025), arXiv:2410.06923 [hep-ph]

  30. [30]

    R. L. Jaffe, Phys. Rev. Lett.38, 195 (1977), [Erratum: Phys.Rev.Lett. 38, 617 (1977)]

  31. [31]

    P. J. Mulders, A. T. M. Aerts, and J. J. De Swart, Phys. Rev. D21, 2653 (1980)

  32. [32]

    V. B. Kopeliovich, B. Schwesinger, and B. E. Stern, Phys. Lett. B242, 145 (1990)

  33. [33]

    Leandri and B

    J. Leandri and B. Silvestre-Brac, Phys. Rev. D51, 3628 (1995)

  34. [34]

    Pepin and F

    S. Pepin and F. Stancu, Phys. Rev. D57, 4475 (1998), arXiv:hep-ph/9710528

  35. [35]

    Q. B. Li, P. N. Shen, Z. Y. Zhang, and Y. W. Yu, Nucl. Phys. A683, 487 (2001), arXiv:nucl-th/0009038

  36. [36]

    C. P. Herzog and J. McKernan, JHEP08, 054 (2003), arXiv:hep-th/0305048

  37. [37]

    Valcarce, H

    A. Valcarce, H. Garcilazo, F. Fernandez, and P. Gon- zalez, Rept. Prog. Phys.68, 965 (2005), arXiv:hep- ph/0502173

  38. [38]

    Ikeda and T

    Y. Ikeda and T. Sato, Phys. Rev. C76, 035203 (2007), arXiv:0704.1978 [nucl-th]

  39. [39]

    S. L. Olsen, Front. Phys. (Beijing)10, 121 (2015), arXiv:1411.7738 [hep-ex]

  40. [40]

    Francis, J

    A. Francis, J. R. Green, P. M. Junnarkar, C. Miao, T. D. Rae, and H. Wittig, Phys. Rev. D99, 074505 (2019), arXiv:1805.03966 [hep-lat]

  41. [41]

    Amarasinghe, R

    S. Amarasinghe, R. Baghdadi, Z. Davoudi, W. Detmold, M. Illa, A. Parreno, A. V. Pochinsky, P. E. Shana- han, and M. L. Wagman, Phys. Rev. D107, 094508 (2023), [Erratum: Phys.Rev.D 110, 119904 (2024)], arXiv:2108.10835 [hep-lat]

  42. [42]

    Wan, S.-Q

    B.-D. Wan, S.-Q. Zhang, and C.-F. Qiao, Phys. Rev. D 105, 014016 (2022), arXiv:2109.07130 [hep-ph]

  43. [43]

    B. Wang, K. Chen, L. Meng, and S.-L. Zhu, Phys. Rev. D110, 014038 (2024), arXiv:2406.06993 [hep-ph]

  44. [44]

    Hammond, W

    B. Hammond, W. Lester, and P. Reynolds,Monte Carlo Methods in ab Initio Quantum Chemistry(World Scien- tific, Singapore, 1994)

  45. [45]

    Foulkes, L

    W. Foulkes, L. Mitas, R. Needs, and G. Rajagopal, Rev. Mod. Phys.73, 33 (2001)

  46. [46]

    Nightingale and C

    M. Nightingale and C. J. Umrigar,Quantum Monte Carlo Methods in Physics and Chemistry(Springer, Vienna, 2014)

  47. [47]

    Ceperley, Rev

    D. Ceperley, Rev. Mod. Phys.67, 279 (1995)

  48. [48]

    Giorgini, L

    S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys.80, 1215 (2008), arXiv:0706.3360 [cond-mat.other]

  49. [49]

    Carlson, J

    J. Carlson, J. B. Kogut, and V. Pandharipande, Phys. Rev. D27, 233 (1983)

  50. [50]

    Carlson, J

    J. Carlson, J. Kogut, and V. Pandharipande, Phys. Rev. D28, 2807 (1983)

  51. [51]

    Y. Bai, S. Lu, and J. Osborne, Phys. Lett. B798, 134930 (2019), arXiv:1612.00012 [hep-ph]

  52. [52]

    M. C. Gordillo, F. De Soto, and J. Segovia, Phys. Rev. D102, 114007 (2020), arXiv:2009.11889 [hep-ph]

  53. [53]

    Semay and B

    C. Semay and B. Silvestre-Brac, Z. Phys. C61, 271 (1994)

  54. [54]

    Silvestre-Brac, Few Body Syst.20, 1 (1996)

    B. Silvestre-Brac, Few Body Syst.20, 1 (1996)

  55. [55]

    M. C. Gordillo, (2026), arXiv:2604.18174 [hep-ph]

  56. [56]

    Topics in quantum mechanics,

    D. Tong, “Topics in quantum mechanics,” Lecture notes, University of Cambridge, sec. 5.4.2, “Generalised Mea- surements”, pp. 176–177

  57. [57]

    M. C. Gordillo and J. M. Alcaraz-Pelegrina, Phys. Rev. D108, 054027 (2023), arXiv:2307.08408 [hep-ph]

  58. [58]

    F. J. Dyson and N.-H. Xuong, Phys. Rev. Lett.13, 815 (1964)

  59. [59]

    Silvestre-Brac and J

    B. Silvestre-Brac and J. Leandri, Phys. Rev. D45, 4221 (1992)

  60. [60]

    M. C. Gordillo, J. Segovia, and J. M. Alcaraz-Pelegrina, Phys. Rev. D110, 094024 (2024), arXiv:2409.04130 [hep- ph]

  61. [61]

    E. G. J. Kessler, M. S. Dewey, R. D. Deslattes, A. Henins, H. G. B¨ orner, M. Jentschel, C. Doll, and H. Lehmann, Physics Letters A255, 221 (1999)