REVIEW 2 major objections 5 minor 36 references
No evidence found for strangeness-rich dibaryons in bottomonium decays, with first upper limits on their production set at the 10^-7 to 10^-6 level in branching fraction.
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 · deepseek-v4-flash
2026-08-02 12:51 UTC pith:TTE6VT2X
load-bearing objection First limits on Ξ0p, Ω−p, Ω−n dibaryons in Υ(1S)/Υ(2S) decays — a careful, honest null result whose quoted limits carry an unquantified model dependence from the assumed phase space for the accompanying hadrons. the 2 major comments →
Search for Xi⁰p, Ω^- p, and Ω^- n dibaryons in Upsilon(1S) and Upsilon(2S) decays at Belle
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors find that in 102 million Upsilon(1S) and 158 million Upsilon(2S) decays, there is no evidence for a Ξ0p, Ω−p, or Ω−n dibaryon in the mass window extending from 30 MeV below to 30 MeV above each baryon-pair threshold. All observed yields are consistent with background. The first upper limits on the branching fractions of these decays are set at 90% confidence, ranging from about 10^-7 to 10^-6, thereby constraining the production of multistrange dibaryons in gluon-rich bottomonium hadronization.
What carries the argument
The search reconstructs each dibaryon candidate from its exclusive decay chain (e.g., Ξ0p → π0Λp, Ω−p → Ξ0Λ, Ω−p → Ω−p) using kinematic fits, then scans the invariant mass relative to the two-baryon threshold in 2 MeV steps from -30 to +30 MeV. The conversion of observed yields into branching fractions relies on B = N_sig/(2N_Υ ε_rec B_sub), with reconstruction efficiencies computed under an assumed phase-space model in which the dibaryon is produced accompanied by a 4ππ0 hadronic system; the profile-likelihood q_μ test statistic converts the null yields into 90% confidence upper limits. This machinery is what turns a visually flat invariant-mass spectrum into quantitative production-rate bo
Load-bearing premise
The reported limits assume that a produced dibaryon would decay through the specific reconstructed chains and that its production is accompanied by exactly the pion multiplicity assumed in simulation; if real dibaryons prefer different decay or production patterns, the limits would not apply.
What would settle it
Taking ten times more data at the same collision energies—or else looking in hadronic collisions—and observing a narrow resonant structure in the Ξ0p, Ω−p, or Ω−n invariant mass spectrum with a production rate above the reported 90% limits would directly falsify the null claim.
If this is right
- If correct, these limits rule out production of near-threshold Ξ0p, Ω−p, and Ω−n dibaryons in bottomonium decays with branching fractions above roughly 10^-6.
- The results provide the first experimental constraints on multistrange dibaryon formation in gluon-rich hadronization, complementing searches in hadronic and nuclear reactions.
- The achieved sensitivity is comparable to the measured antideuteron production rate in the same decays, suggesting the search would have seen a similar multibaryon signal if such dibaryons were produced at comparable rates.
- Models predicting shallow bound states with binding energies up to a few MeV are now constrained at the branching-fraction level, offering a new benchmark for baryon-baryon interaction models.
Where Pith is reading between the lines
- Because the search window covers only ±30 MeV around threshold, a dibaryon bound more deeply than about 30 MeV would fall outside the scan and would not be excluded by these limits.
- The limits are tied to the assumed signal model: if the true production mechanism populates a different number of additional pions or a different momentum distribution, the reconstruction efficiency—and hence the limits—could shift by a factor that is not included in the quoted systematic uncertainties.
- The Ω−n limit applies only to a bound state decaying strongly to ΞΛ; a neutron-star-relevant Ω−n state with a different dominant decay mode would be invisible to this search, and a dedicated search with neutron detection would be needed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a search for dibaryon states with strangeness content (Ξ0p, Ω−p, Ω−n) in Υ(1S) and Υ(2S) decays using the full Belle dataset (102M Υ(1S) and 158M Υ(2S) decays). Both bound-state (below-threshold) and unbound near-threshold hypotheses are reconstructed through hyperon decay chains, with optimized selections developed on simulation and background control regions before signal inspection. No significant signal is observed in any channel, and 90% confidence-level upper limits on the branching fractions are set as a function of the mass difference from the corresponding baryon-pair threshold, with quoted values in the O(10^-7)–O(10^-6) range. The limits are obtained with a profile-likelihood construction with Gaussian-constrained nuisance parameters, and the total systematic uncertainties are 5–6%.
Significance. If the result holds, these are the first upper limits for these multistrange dibaryons in bottomonium decays, providing constraints on attractive ΞN and ΩN interactions that are relevant for neutron-star matter and for discriminating between theoretical models. The analysis is careful and methodical: selection is optimized before signal regions are examined, the systematic budget is documented in detail, and observed limits are consistent with expectations within 2σ. Presenting limits as a function of mass difference from threshold enables direct comparison with theory and with other experiments. The main caveat is that the reconstruction efficiency—and therefore the numerical limits—depends on an assumed signal-production model that is not fully varied in the systematic budget.
major comments (2)
- [Simulated-samples paragraph; Eq. (1); Table I] The signal simulation assumes a fixed final state 'H′ B1 B2 4ππ0', described as a 'representative five-body hadronic system' to restrict phase space. The reconstruction efficiency ε_rec enters the branching-fraction limit linearly in Eq. (1), yet the systematic assigned to 'Reconstruction efficiency' in Table I (0.1–0.7%) is derived only from varying the antibaryon combinations, not from the pion multiplicity or the phase-space model. If the true production has a different hadronic multiplicity (e.g., 2π or 6π) or different momentum distribution, ε_rec could change by more than the quoted 5–6% total systematic, shifting the central numerical limits. The authors should add a systematic explicitly accounting for the assumed multiplicity/kinematics, or quote limits under a conservative efficiency choice, or clearly state in the abstract and conclusions that the limits are conditional on the
- [ΩN bound-state reconstruction paragraph] For the Ω−p and Ω−n bound-state hypotheses, the reconstruction assumes a strong decay to Ξ0Λ (Ω−p) or Ξ−Λ (Ω−n). This is a model assumption, motivated by the threshold difference, but an alternative dominant decay mode (e.g., involving Ξ* resonances or other strangeness-conserving channels) would remove sensitivity entirely. The manuscript should explicitly state that the reported limits are upper limits on the product B(Υ→Dibaryon+X)×B(Dibaryon→ΞΛ) for these two channels, and ideally discuss the theoretical basis for the assumed decay mode and its possible branching fraction.
minor comments (5)
- [Abstract] The abstract claims 'first 90% confidence-level upper limits' without noting the model dependence of the signal efficiency. A short qualifier such as 'under the assumed production model' would be more precise.
- [Figure 2 and text after Eq. (1)] The 'gap at zero mass difference in the Ω−p channel' is explained as reflecting 'the reduced reconstruction efficiency of the bound-state hypothesis... relative to the all-charged unbound final state.' This explanation is not self-evident; clarify whether the point is omitted due to a fit failure, an efficiency cut, or an intentionally excluded scale.
- [All figures] Several figures carry the label 'Preliminary'. If this manuscript is intended for journal submission, these labels should be removed or clearly explained.
- [Supplemental Tables 1 and 2] The best-fit branching fractions include negative values. While this is an expected artifact of unconstrained background-only fits, a one-sentence note in the supplemental material would help readers interpret these entries as unphysical.
- [Reference [29]] The supplemental-material link is given as 'XXXX-link-provided-by-PRL'. Ensure the URL is filled in before submission.
Circularity Check
No significant circularity: the limits follow from data through an explicitly defined branching-fraction formula with independently measured inputs.
full rationale
This is a search measurement, not a derivation. The central formula, B = Nsig/(2N_Upsilon eps_rec B_sub) (Eq. 1), is the definition of a branching fraction: Nsig comes from an unbinned extended maximum-likelihood fit to the invariant-mass spectra, N_Upsilon is taken from a dedicated Belle measurement [36], eps_rec is determined from signal MC, and B_sub is taken from the PDG [26]. The upper limits are produced by the q_mu profile-likelihood construction; no parameter is fitted to a subset of the data and then relabeled as a prediction, and the signal and background yields are explicitly allowed to float while all shape parameters are fixed from simulation. The signal-model assumption 'a representative five-body hadronic system, 4 pi pi0, is assumed for the additional particles in order to restrict the available phase space' is an acknowledged external modeling choice affecting eps_rec, and the paper assigns a systematic for the variation over antibaryon combinations allowed by charge and strangeness conservation. The unquantified dependence on the pion multiplicity is a model-robustness caveat that would affect the numerical limits, but it is not circularity: eps_rec is not derived from the branching fraction being constrained. Self-citations in the paper are to published Belle calibration studies (Lambda selection, pi0 selection, Upsilon counting) and software, which are independent external inputs, not load-bearing arguments that reduce to an unverified self-citation or to a uniqueness theorem from the same authors. The central claim—'No significant signals are observed'—is a direct data observation, and the reported first upper limits are new outputs conditional on stated assumptions, not restatements of those assumptions by construction.
Axiom & Free-Parameter Ledger
free parameters (3)
- Assumed width of ΩN bound-state signal =
5 MeV/c²
- Representative binding energies in signal MC =
2 MeV (Ξ⁰p), 10 MeV (ΩN)
- Phase-space model for additional final-state particles =
4ππ⁰ (5-body)
axioms (6)
- domain assumption Monte Carlo simulation (EvtGen/Pythia + Geant3) accurately models signal and background shapes and efficiencies
- domain assumption Υ(1S) and Υ(2S) decays share similar hadronization dynamics via the three-gluon intermediate state, justifying sample combination
- domain assumption The e⁺e⁻→qq continuum sample at √s = 10.52 GeV models the background under the resonances
- domain assumption Bound ΩN states decay strongly to ΞΛ (Ξ⁰Λ for Ω⁻p, Ξ⁻Λ for Ω⁻n)
- standard math CP symmetry holds, allowing baryon and antibaryon channels to be combined
- domain assumption PDG branching fractions and the Belle-measured Υ(1S)/Υ(2S) event counts are correct
Cite this review
Pith. "Pith review of Search for $\Xi^0p$, $\Omega^- p$, and $\Omega^- n$ dibaryons in $\Upsilon(1S)$ and $\Upsilon(2S)$ decays at Belle." pith.science (2026). https://pith.science/paper/TTE6VT2X
@misc{pith2026260529778,
author = {Pith},
title = {Pith review of: Search for $\Xi^0p$, $\Omega^- p$, and $\Omega^- n$ dibaryons in $\Upsilon(1S)$ and $\Upsilon(2S)$ decays at Belle},
year = {2026},
howpublished = {\url{https://pith.science/paper/TTE6VT2X}},
note = {Machine review of arXiv:2605.29778}
}
read the original abstract
We search for $\Xi^0p$, $\Omega^-p$, and $\Omega^-n$ dibaryon states in $\Upsilon(1S)$ and $\Upsilon(2S)$ decays, probing mass regions near the corresponding baryon-pair thresholds. Multistrange baryon-baryon interactions are relevant to neutron-star matter but remain largely unconstrained. Experimental and theoretical studies supporting attractive $\Xi N$ and $\Omega N$ interactions, where $N$ denotes a nucleon, motivate searches for weakly bound states. We use samples of $102$ million $\Upsilon(1S)$ and $158$ million $\Upsilon(2S)$ decays collected with the Belle detector at the KEKB asymmetric-energy $e^+e^-$ collider. No significant signals are observed, and the first $90$% confidence-level upper limits are set on the branching fractions of $\Upsilon(1S)$ and $\Upsilon(2S)$ decays to $\Xi^0p$, $\Omega^-p$, and $\Omega^-n$ dibaryon states, at the level of $O(10^{-7})$-$O(10^{-6})$, depending on the channel and the assumed mass difference from the corresponding baryon-pair threshold.
Figures
Reference graph
Works this paper leans on
-
[1]
Demorestet al., Nature4671081 (2010)
P. Demorestet al., Nature4671081 (2010)
2010
-
[2]
Djapo, B.-J
H. Djapo, B.-J. Schaefer, and J. Wambach, Phys. Rev. C81035803 (2010)
2010
-
[3]
Lonardoni, A
D. Lonardoni, A. Lovato, S. Gandolfi, and F. Pederiva, Phys. Rev. Lett.114092301 (2015)
2015
-
[4]
Acharyaet al.(ALICE Collaboration), Phys
S. Acharyaet al.(ALICE Collaboration), Phys. Rev. Lett.123112002 (2019)
2019
-
[5]
Acharyaet al.(ALICE Collaboration), Nature588 232 (2020), [Erratum: Nature 590, E13 (2021)]
S. Acharyaet al.(ALICE Collaboration), Nature588 232 (2020), [Erratum: Nature 590, E13 (2021)]
2020
-
[6]
Morita, A
K. Morita, A. Ohnishi, F. Etminan, and T. Hatsuda, Phys. Rev. C94031901 (2016), [Erratum: Phys.Rev.C 100, 069902 (2019)]
2016
-
[7]
Iritaniet al.(HAL QCD Collaboration), Phys
T. Iritaniet al.(HAL QCD Collaboration), Phys. Lett. B792284 (2019)
2019
-
[8]
Sekihara and T
T. Sekihara and T. Hashiguchi, Phys. Rev. C108 065202 (2023)
2023
-
[9]
Sasakiet al.(HAL QCD Collaboration), Nucl
K. Sasakiet al.(HAL QCD Collaboration), Nucl. Phys. A998121737 (2020)
2020
-
[10]
Haidenbauer and U.-G
J. Haidenbauer and U.-G. Meißner, Eur. Phys. J. A55 23 (2019)
2019
-
[11]
M. M. Nagels, T. A. Rijken, and Y. Yamamoto, Phys. Rev. C102054003 (2020)
2020
-
[12]
Fujiwara, Y
Y. Fujiwara, Y. Suzuki, and C. Nakamoto, Prog. Part. Nucl. Phys.58439 (2007)
2007
-
[13]
Hiyamaet al., Phys
E. Hiyamaet al., Phys. Rev. Lett.124092501 (2020)
2020
-
[14]
J. P. Leeset al.(BaBar Collaboration), Phys. Rev. D 89111102 (2014)
2014
-
[15]
D. M. Asneret al.(CLEO Collaboration), Phys. Rev. D 75012009 (2007)
2007
-
[16]
Adachiet al.(Belle and Belle II Collaborations), Phys
I. Adachiet al.(Belle and Belle II Collaborations), Phys. Rev. Lett.135041901 (2025)
2025
-
[17]
B. H. Kimet al.(Belle Collaboration), Phys. Rev. Lett. 110222002 (2013)
2013
-
[18]
J. P. Leeset al.(BaBar Collaboration), Phys. Rev. Lett.122072002 (2019)
2019
-
[19]
Abashianet al., Nucl
A. Abashianet al., Nucl. Instr. and Meth. A479117 (2002)
2002
-
[20]
D. J. Lange, Nucl. Instrum. Meth. A462152 (2001)
2001
-
[21]
Sjostrand, S
T. Sjostrand, S. Mrenna, and P. Z. Skands, JHEP05 026 (2006)
2006
-
[22]
Brunet al.CERN-DD-EE-84-1 (1987)
R. Brunet al.CERN-DD-EE-84-1 (1987)
1987
-
[23]
Kuhret al.(Belle II Framework Software Group), Comput
T. Kuhret al.(Belle II Framework Software Group), Comput. Softw. Big Sci.31 (2019)
2019
-
[24]
M. Gelb, T. Keck, M. Primet al., Comput. Software Big Sci.29 (2018)
2018
-
[25]
Krohnet al.(Belle II analysis software Group), Nucl
J.-F. Krohnet al.(Belle II analysis software Group), Nucl. Instrum. Meth. A976164269 (2020)
2020
-
[26]
Navaset al.(Particle Data Group), Phys
S. Navaset al.(Particle Data Group), Phys. Rev. D 110030001 (2024)
2024
-
[27]
Abeet al.(Belle Collaboration), Phys
K. Abeet al.(Belle Collaboration), Phys. Rev. D65 091103 (2002)
2002
-
[28]
Punzi, eConfC030908MODT002 (2003)
G. Punzi, eConfC030908MODT002 (2003)
2003
-
[29]
See Supplemental Material at XXXX-link-provided-by-PRL for additional plots and tables
-
[30]
J. E. Gaiser, PhD thesis, SLAC-R-255 (1982), Appendix F; T. Skwarnicki, PhD thesis, DESY F31-86-02 (1986), Appendix E
1982
-
[31]
Ikedaet al.(Belle Collaboration), Nucl
H. Ikedaet al.(Belle Collaboration), Nucl. Instrum. Meth. A441401 (2000)
2000
-
[32]
Cowan, K
G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Eur. Phys. J. C711554 (2011), [Erratum: Eur.Phys.J.C 73, 2501 (2013)]
2011
-
[33]
Rodrigueset al., EPJ Web Conf.24506028 (2020)
E. Rodrigueset al., EPJ Web Conf.24506028 (2020)
2020
-
[34]
Katoet al.(Belle Collaboration), Phys
Y. Katoet al.(Belle Collaboration), Phys. Rev. D94 032002 (2016)
2016
-
[35]
Bhardwajet al.(Belle Collaboration), Phys
V. Bhardwajet al.(Belle Collaboration), Phys. Rev. D 99111101 (2019)
2019
-
[36]
Jiaet al.(Belle Collaboration), Phys
S. Jiaet al.(Belle Collaboration), Phys. Rev. D95 012001 (2017). 1 Supplemental material Figure 1 shows the reconstruction and selection efficiency for all channels as a function of the mass difference from the corresponding baryon-pair threshold. 30 20 10 0 10 20 30 M (MB1 + MB2) (MeV) 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 Efficiency Preliminary 0...
2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.