REVIEW 3 major objections 6 minor 32 references
$\Omega$-deuteron Interaction in Folding Model
T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read A single folding model, using a lattice QCD Omega-nucleon potential and four deuteron wave functions, produces two-body binding energies of 6.7 to 19.2 MeV and supports the conclusion that the Omega-deuteron system is deeply bound.
desk verdict Folding-model Omega-d calculation is transparent and useful, but the paper's 'upper bound' has the inequality reversed—those numbers are lower bounds on the three-body binding. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the Watanabe single-folding ansatz, in which the $\Omega$-deuteron potential is built by integrating the $\Omega$-nucleon potential against the deuteron density constructed from four analytic S-state deuteron wave functions: two Hulthén forms and modified Reid93 and Argonne v18 forms. The folded potential is then fitted to a Woods-Saxon shape, and the bound state of that fitted potential is computed from the two-body Schrödinger equation. The upper-bound interpretation comes from the argument that folding with an unperturbed deuteron wave function omits distortion and other attractive effects, so the true three-body potential is more attractive than the folded one.
What would settle it
A full three-body calculation with the same lattice QCD $\Omega$-nucleon potential but retaining the deuteron D-state, tensor force, and nucleon distortion would produce a definite binding energy; if that binding fell below the smallest folding value of about 6.7 MeV, the paper's upper-bound claim would collapse, since the true potential must be more attractive than the folding approximation. Alternatively, a measurement of $\Omega$-deuteron correlations in heavy-ion collisions that showed no bound state near the predicted energy range would falsify the deeply-bound conclusion.
Extended reading notes
Core claim
The central claim is that the effective central potential for an $\Omega$ orbiting a deuteron in the maximal-spin (0)5/2+ state is well approximated by a Woods-Saxon form once the lattice QCD $\Omega$-nucleon potential is folded over a deuteron S-state wave function. Solving the two-body Schrödinger equation with that potential gives binding energies of 6.7 MeV for the Reid93-based wave function, 6.8 MeV for the Argonne v18-based wave function, and 15.7 and 19.2 MeV for the two Hulthén forms. The paper asserts that, since the real potential must be more attractive than the folding potential, each of these values is an upper bound for the three-body OmegaNN binding energy, and it compares them with the 20.9 MeV reported by a full three-body Faddeev calculation. The conclusion is that the $\Omega$-deuteron system in this channel is deeply bound, but the two-body folding approximation is too crude for precise quantitative predictions.
Load-bearing premise
The calculation assumes the deuteron remains an unperturbed free deuteron, using only its S-state component and only the spin-independent central part of the $\Omega$-nucleon interaction, and it takes one particular parameterization of the lattice QCD potential as representative; if deuteron distortion, D-state or tensor correlations, or a different parameter set materially change the short-range interaction, the computed binding energies (6.7-19.2 MeV) could shift by several MeV.
Editorial extensions
If this is right
- The Omega-deuteron system in the (0)5/2+ state is deeply bound, with a binding energy between roughly 7 and 21 MeV depending on the interaction model.
- The two-body folding model provides cheap estimates for three-body binding but cannot replace full three-body calculations when precision is needed.
- The spread from 6.7 to 19.2 MeV across deuteron wave functions quantifies how strongly the short-range deuteron structure controls the binding energy.
- The Woods-Saxon form for the effective Omega-deuteron potential offers a simple input for future scattering and correlation studies.
- If the upper-bound logic holds, even the smallest computed value (6.7 MeV) implies a genuine bound state, not just a virtual one.
Reading between the lines
- A decisive test would be to repeat the same folding procedure with the full deuteron wave function including the D-state; the resulting upper bound should move upward relative to the S-state-only values if the paper's reasoning is correct.
- The strong wave-function sensitivity suggests that the short-range shape of the deuteron, not the long-range Omega-nucleon tail, dominates the binding; comparing the depths and radii of the folded potentials across models could reveal whether a common scale emerges.
- The same folding prescription could be applied to other multi-strange dibaryon candidates, such as Xi-deuteron or Omega-alpha, to see whether the two-body bound-state energies track the three-body results in the same pattern.
- A practical extension is to use the actual folded potential rather than its Woods-Saxon fit in a three-body calculation; if the resulting binding differs from the 20.9 MeV Faddeev value by more than the fit's accuracy, the upper-bound claim would need revision.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript constructs an effective Omega-deuteron potential by single folding the HAL QCD spin-2 Omega-N potential over four S-state deuteron wave functions. The Omega-N input is the analytic fit of Eq. (1) with the P1 parameter set of Table I. The folded central potential of Eq. (3) is fitted to a Woods-Saxon form, Eq. (10), and the two-body Schrodinger equation is solved to obtain binding energies of 15.7, 19.2, 6.7, and 6.8 MeV for the Hulthen-1, Hulthen-2, Reid93, and Argonne v18 deuteron wave functions, respectively. These are compared with the three-body Faddeev result of 20.9 MeV from Ref. [9]. The paper concludes that the two-body folding values are upper bounds on the three-body binding energy and that the Omega-d system is deeply bound, while also emphasizing that the single-folding reduction is inadequate for quantitative accuracy.
Significance. The calculation is transparent and uses external lattice QCD input; it demonstrates that a simple two-body folding model can reproduce the qualitative feature of deep binding and that the result depends strongly on the short-range behavior of the deuteron wave function. These are useful reference points for more refined few-body calculations. However, the central interpretation that the computed binding energies are upper bounds is mathematically reversed: they are variational lower bounds. This weakens the logical force of the comparison and requires correction in the abstract, the results section, and the conclusions. The numerical values themselves remain consistent with deep binding, so the qualitative conclusion survives, but the paper's central claim needs to be reframed.
major comments (3)
- [Abstract; Section IV (after Table II); Section V] The direction of the variational inequality is reversed. For the trial wave function Psi(R,r) = phi(R) psi_d(r), the expectation value is E_trial = E_d + E_2b = -B_d - B_2b, and the variational principle gives E_trial >= E_exact. Therefore B_exact = -E_exact >= -E_trial = B_d + B_2b. Since the Table II values are quoted with respect to the Omega-N-N threshold, they are lower bounds, not upper bounds, on the three-body binding energy. This correction applies to the abstract, to the sentence 'Since the real potential must be more attractive than the folding potential, the resulting energy can only be an upper bound...' in Section IV, and to the analogous statement in Section V. The reported numbers remain consistent with deep binding, but the logical status of the comparison to the 20.9 MeV Faddeev value changes.
- [Table II; Eq. (10)] The fitted diffuseness parameters for the Hulthen-1 and Hulthen-2 cases are negative (c = -0.19 fm and c = -0.31 fm). With negative c, Eq. (10) does not describe a standard Woods-Saxon potential: instead of a monotone attractive pocket, it produces a potential whose magnitude increases from r = 0 up to r approximately R and then saturates. Since the binding energies are obtained after solving the Schrodinger equation with this fitted potential, the claim that the effective potential has a 'simple Wood-Saxon form' is not supported for those two cases. The authors should justify the negative-diffuseness fits, use a different functional form, or demonstrate that the binding energies are insensitive to the fitting ansatz.
- [Section III; Table II] All numerical results are generated with only the P1 parameter set of Table I and with the deuteron restricted to its S-state component. The HAL QCD input contains four parameter sets, P1-P4, and the D-state probability is discarded after Eq. (5). A single-parameter-set, S-wave-only calculation cannot establish the quoted 6.7-19.2 MeV spread as a robust theoretical range. At minimum, the authors should repeat the folding and binding calculation for P2-P4 and provide an estimate of the D-state contribution, for example by folding the full u^2(r) + w^2(r) density. Without such a sensitivity check, the comparison with the 20.9 MeV Faddeev result is not as strong as the text implies.
minor comments (6)
- [Abstract; Section II] The abstract describes the Omega-N interaction as 'separable,' but Eq. (1) is a local central potential with Gaussian and Yukawa terms. The word 'separable' appears to be a misnomer and should be replaced by 'analytic' or 'local.'
- [Eq. (9)] The sign convention for the Coulomb term is not specified: Eq. (9) uses a '+-' symbol. The following sentence about the Coulomb interaction increasing binding for systems containing a proton is ambiguous, since the Omega-proton force is attractive but the Omega-neutron force has no Coulomb term. Please specify the signs and clarify the sentence.
- [Table II; Section IV] The table heading uses B_Omega D while the text refers to B_2b, Omega d, and the caption states the values are with respect to the Omega-N-N threshold. Please clarify once whether the quoted numbers include the deuteron binding energy, and use a single notation consistently.
- [References] Several references have incomplete bibliographic data: Refs. [7], [13], and [19] contain '000' placeholders, and Ref. [8] lists page '0'. These should be completed before publication.
- [Fig. 2; Eq. (10)] The quality of the Woods-Saxon fits is not quantified. Reporting the chi-squared per degree of freedom or a similar measure for each of the four fits would help the reader judge whether the fitted form is actually adequate.
- [Throughout] The spelling 'Wood-Saxon' should be corrected to 'Woods-Saxon,' following standard usage and Ref. [27].
Circularity Check
No significant circularity: the Ωd binding energies are computed from an external HAL QCD ΩN potential and standard deuteron wave functions, not fitted to the three-body value; the upper/lower-bound wording error is a physics mistake, not a circular derivation.
full rationale
The derivation chain is not circular. The ΩN interaction is an external input taken from the HAL QCD lattice calculation (Refs. [17,18]), with parameters listed in Table I and the P1 parameterization selected as representative. Although one of the present authors is a coauthor of those lattice papers, the lattice potential is an independent, externally obtained input: it does not include the Ωd binding energy, is not fitted in this paper, and is not equivalent to the target result. The folded potential is constructed by Eqs. (2) and (3) using standard deuteron wave functions (Hulthen, Reid93, Av18), and the Woods-Saxon parameters in Eq. (10) are fitted to the folded potential itself, not to any binding-energy datum. The binding energies in Table II are then obtained by solving the Schrödinger equation and compared with the independent three-body Faddeev result of 20.9 MeV from Ref. [9]. No parameter is fitted to the three-body binding energy, and no equation reduces by construction to an input. The claim that the folding energy is an upper bound on the three-body binding is physically misstated (the variational inequality actually makes these binding energies lower bounds), but that is a correctness error, not a circularity. Thus the paper is self-contained with respect to circularity.
Assumptions & free parameters
free parameters (4)
- Woods-Saxon potential depth V0 =
538, 675, 143, 148 MeV for Hulthen-1, Hulthen-2, Reid93, Av18
- Woods-Saxon radius R =
0.98, 0.94, 0.89, 0.90 fm
- Woods-Saxon diffuseness c =
-0.19, -0.31, 1.39, 1.34 fm
- Coulomb screening radius r0 =
50 fm
assumptions (5)
- standard math The Omega-deuteron relative motion is described by the non-relativistic Schrodinger equation with the folded potential (Section IV).
- ad hoc to paper The HAL QCD Omega-nucleon 5S2 potential (Eq. 1, Table I) is a faithful input and the P1 parameter set is representative.
- domain assumption The deuteron is described by free-space S-wave wave functions only; D-state probability is neglected.
- domain assumption Watanabe folding gives an upper bound on the true three-body binding energy.
- domain assumption The (0)5/2+ maximal-spin channel is the relevant state and spin-independent central interactions suffice.
Cite this review
Pith. "Pith review of $\Omega$-deuteron Interaction in Folding Model." pith.science (2026). https://pith.science/paper/LDBT53XM
@misc{pith2026190811484,
author = {Pith},
title = {Pith review of: $\Omega$-deuteron Interaction in Folding Model},
year = {2026},
howpublished = {\url{https://pith.science/paper/LDBT53XM}},
note = {Machine review of arXiv:1908.11484}
}
abstract
A simple single folding model for $\Omega$-deuteron with maximal spin $\left(I\right)J^{P}=\left(0\right)5/2^{+}$ is investigated. $\Omega$ is assumed to orbit an unperturbed deuteron in a $\Omega$-deuteron potential based on a separable $\Omega$-nucleon potential from lattice QCD. We show that the effective central folding potential of $\Omega d$ in the $^{5}S_{2}$ channel has a simple Wood-Saxon form, and approximate the upper bound for the binding energy of $\Omega$ particle on a deuteron. In order to investigate how changes in the wave functions affect the results, we consider four analytical forms for $S$-state deuteron wave functions, i.e., two widely used Hulth\'en forms, as well as the modified Reid93 and Argonne v18 forms. Our calculations of binding energy from simple two-body approximation are compared with the results reported for the three-body problem; it is confirmed that the $\Omega d$ system is deeply bound. Although the single folding model reduces the three-body problem to a two-body problem, this simplification is inadequate.
Figures
Reference graph
Works this paper leans on
- [9]
- [1]
-
[2]
H. Garcilazo, A. Valcarce, and J. Vijande, Phys. Rev. C 94, 024002 (2016)
work page 2016
- [3]
-
[4]
N. N. Kolesnikov and V. A. Kopylov, Sov. Phys. J. 31, 210 (1988)
work page 1988
-
[5]
S. P. Verma and D. P. Sural, Phys. Rev. C 20, 781 (1979)
work page 1979
- [6]
- [7]
Show all 32 references
-
[8]
S. M. M. Nejad and A. Armat, International Journal of Modern Physics E 0, 0 (0), https://doi.org/10.1142/S0218301319500873
-
[10]
Etminan and M
F. Etminan and M. M. Firoozabadi, Modern Physics Letters A 29, 1450177 (2014), https://doi.org/10.1142/S0217732314501776
2014 doi
-
[11]
Nemura (HAL QCD), International Journal of Modern Physics E 23, 1461006 (2014), https://doi.org/10.1142/S0218301314610060
H. Nemura (HAL QCD), International Journal of Modern Physics E 23, 1461006 (2014), https://doi.org/10.1142/S0218301314610060
2014 doi
-
[12]
Sasaki, S
K. Sasaki, S. Aoki, T. Doi, S. Gongyo, T. Hatsuda, Y. Ikeda, T. Inoue, T. Iritani, N. Ishii, D. Kawai, K. Mu- rano, and H. Nemura (HAL QCD), EPJ Web Conf. 175, 05010 (2018)
2018
-
[13]
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 (HAL QCD), arXiv:1912.08630 [hep-lat] , 000 (2019), arXiv:1912.08630 [hep-lat]
2019 arXiv
-
[14]
J. G. Congleton, J. Phys. G: Nucl. Part. Phys. 18, 339 (1992)
1992
-
[15]
Morita, A
K. Morita, A. Ohnishi, F. Etminan, and T. Hatsuda, Phys. Rev. C 94, 031901(R) (2016)
2016
-
[16]
Adam (STAR), Phys
J. Adam (STAR), Phys. Lett. B 790, 490 (2019)
2019
-
[17]
Iritani, S
T. Iritani, S. Aoki, T. Doi, F. Etminan, S. Gongyo, T. Hatsuda, Y. Ikeda, T. Inoue, N. Ishii, T. Miyamoto, and K. Sasaki (HAL QCD), Phys. Lett. B 792, 284 (2019)
2019
-
[18]
Etminan, H
F. Etminan, H. Nemura, S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, T. Inoue, N. Ishii, K. Murano, and K. Sasaki (HAL QCD), Nucl. Phys. A 928, 89 (2014)
2014
- [19]
-
[20]
Watanabe, Nucl
S. Watanabe, Nucl. Phys. 50, 484 (1958). 5
1958
-
[21]
P. W. Keaton and D. D. Armstrong, Phys. Rev. C 8, 1692 (1973)
1973
-
[22]
J. M. Blatt and V. F. Weisskopf, in Theoretical Nuclear Physics (Springer, New York, NY, 1958) Chap. 2, pp. 99–101
1958
-
[23]
V. I. Zhaba, Modern Physics Letters A 31, 1650139 (2016)
2016
-
[24]
R. J. Adler, T. K. Das, and A. F. Filho, Phys. Rev. C 16, 1231 (1977)
1977
-
[25]
Mr´ owczy´ nski, Physics Letters B277, 43 (1992)
S. Mr´ owczy´ nski, Physics Letters B277, 43 (1992)
1992
-
[26]
Gongyo, K
S. Gongyo, K. Sasaki, S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, T. Inoue, T. Iritani, N. Ishii, T. Miyamoto, and H. Nemura, Phys. Rev. Lett. 120, 212001 (2018)
2018
-
[27]
M. A. Melkanoff, J. S. Nodvik, D. S.Saxon, and D. Woods, Phys. Rev. 106, 793 (1957)
1957
-
[28]
Miyamoto, S
T. Miyamoto, S. Aoki, T. Doi, S. Gongyo, T. Hatsuda, Y. Ikeda, T. Inoue, T. Iritani, N. Ishii, D. Kawai, K. Mu- rano, H. Nemura, and K. Sasaki (HAL QCD), Nucl. Phys. A 971, 113 (2018)
2018
-
[29]
S. Aoki, B. Charron, T. Doi, T. Hatsuda, T. Inoue, and N. Ishii, Phys. Rev. D 87, 034512 (2013)
2013
-
[30]
H. D. Vries, C. D. Jager, and C. D. Vries, At. Data Nucl. Data Tables 36, 495 (1987)
1987
-
[31]
M. E.-A. Farid and M. Hassanain, Nucl. Phys. A 678, 39 (2000)
2000
-
[954]
5 MeV/c 2, respectively [17, 26])
7 MeV/c 2 and 1711. 5 MeV/c 2, respectively [17, 26]). One can see from Table II that the binding energies cal- culated by Hulth´ en-1 (15. 7 MeV) and Hulth´ en-2 (19 . 2 4 MeV) DWFs are more closer to the binding energy from three-body calculations. Since the real potential m...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.