Pith. sign in

REVIEW 2 major objections 5 minor 61 references

Diffusion Monte Carlo study of deuteron-like fully light hexaquarks

T0 review · 2 major / 5 minor · reviewed 2026-07-13 · grok-4.5

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

Pith's one-line read 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.

desk verdict 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. read the letter →

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

M.C. Gordillo , J. Segovia This is my paper · ORCID
classification hep-phhep-exhep-latnucl-exnucl-th
keywords hexaquarksdibaryonsdiffusionMonteCarloconstituentquarkmodeldeuteron-likestatesmultiquarkspectroscopyAL1potential
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

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.

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.

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

Extended reading notes

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.

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.

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

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 · score 0.0 of 10

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.

Assumptions & free parameters 2 free parameters · 4 assumptions · 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.
assumptions (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.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Diffusion Monte Carlo study of deuteron-like fully light hexaquarks." pith.science (2026). https://pith.science/paper/FDPTRHE4

@misc{pith2026260709288,
  author       = {Pith},
  title        = {Pith review of: Diffusion Monte Carlo study of deuteron-like fully light hexaquarks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FDPTRHE4}},
  note         = {Machine review of arXiv:2607.09288}
}
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 the authors.

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. FIG. 2. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

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)

    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)

    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]

    arXiv 2021

  4. [4]

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

    arXiv 2016

  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]

    arXiv 2017

  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]

    arXiv 2018

  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

    arXiv 2019

  8. [8]

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

    arXiv 2020

Show all 61 references
  1. [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]

    2021 arXiv

  2. [10]

    Chen, Phys

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

    2022 arXiv

  3. [11]

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

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

    2023 arXiv

  4. [12]

    Mai, U.-G

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

    2023 arXiv

  5. [13]

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

    2022 arXiv

  6. [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]

    2023 arXiv

  7. [15]

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

    2022 arXiv

  8. [16]

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

    2021 arXiv

  9. [17]

    Huang, C

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

    2023

  10. [18]

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

    2023 arXiv

  11. [19]

    Zou, Sci

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

    2021 arXiv

  12. [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]

    2021 arXiv

  13. [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]

    2024 arXiv

  14. [22]

    Johnson, I

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

    2024 doi arXiv

  15. [23]

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

    2025 arXiv

  16. [24]

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

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

    2025 arXiv

  17. [25]

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

    2025 arXiv

  18. [26]

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

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

    2025 arXiv

  19. [27]

    Francis, Prog

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

    2025 arXiv

  20. [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]

    2025 arXiv

  21. [29]

    H¨ usken, E

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

    2025 arXiv

  22. [30]

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

    1977

  23. [31]

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

    1980

  24. [32]

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

    1990

  25. [33]

    Leandri and B

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

    1995

  26. [34]

    Pepin and F

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

    1998 arXiv

  27. [35]

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

    2001 arXiv

  28. [36]

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

    2003 arXiv

  29. [37]

    Valcarce, H

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

    2005

  30. [38]

    Ikeda and T

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

    2007 arXiv

  31. [39]

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

    2015 arXiv

  32. [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]

    2019 arXiv

  33. [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]

    2023 arXiv

  34. [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]

    2022 arXiv

  35. [43]

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

    2024 arXiv

  36. [44]

    Hammond, W

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

    1994

  37. [45]

    Foulkes, L

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

    2001

  38. [46]

    Nightingale and C

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

    2014

  39. [47]

    Ceperley, Rev

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

    1995

  40. [48]

    Giorgini, L

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

    2008 arXiv

  41. [49]

    Carlson, J

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

    1983

  42. [50]

    Carlson, J

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

    1983

  43. [51]

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

    2019 arXiv

  44. [52]

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

    2020 arXiv

  45. [53]

    Semay and B

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

    1994

  46. [54]

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

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

    1996

  47. [55]

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

    2026 arXiv

  48. [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

  49. [57]

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

    2023 arXiv

  50. [58]

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

    1964

  51. [59]

    Silvestre-Brac and J

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

    1992

  52. [60]

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

    2024 arXiv

  53. [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)

    1999

Pith tools

Reviewed July 13, 2026 · model on record in the stance chip above.