REVIEW 4 major objections 7 minor 57 references
Mass and Decay Properties of Toponium Using SUSY QM Factorization Method
T0 review · 4 major / 7 minor · reviewed 2026-07-30 · grok-4.5
Pith's one-line read SUSY QM applied to the Cornell potential yields toponium ground-state masses near 344.1 GeV, a 2.57 MeV di-gluon width, and radii of 0.015–0.025 fm.
desk verdict Clean SUSY-QM Cornell application that reproduces its own numbers, but the MeV-precision 1S masses describe a broad threshold bump, not an isolated bound state. 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 SUSY-QM factorization of the radial Schrödinger equation for the Cornell potential, which expands the superpotential in the dimensionless string-tension parameter λ and yields the compact analytic formula |Ψ(0)|² = (1/π a₀³)(1 + 3/2 λ − 31/16 λ²) that controls every mass, radius and decay width.
What would settle it
A precision threshold scan at a future e⁺e⁻ or muon collider that measures the 1S peak position and the η_t → gg (or γγ) partial width; disagreement with 344.11 GeV and 2.57 MeV at the few-MeV level would falsify the O(λ²) SUSY-QM prediction.
Extended reading notes
Core claim
Within the SUSY-QM factorization of the Cornell potential, the ground-state toponium system has spin-averaged mass 344.134 GeV, binding energy −0.9056 GeV, hyperfine splitting 27.2 MeV (giving η_t = 344.114 GeV and Θ_t = 344.141 GeV), dominant di-gluonic width 2.57 MeV, and mean/RMS/most-probable radii all lying between 0.015 and 0.025 fm—all quantities fixed by the O(λ²) expansion of |Ψ(0)|².
Load-bearing premise
That a stable non-relativistic Schrödinger bound state still makes sense when each top quark already decays with a width larger than the computed binding energy, so the system is at best a broad quasi-bound resonance whose wave function is still given by the same expansion.
Editorial extensions
If this is right
- Pseudoscalar and vector toponium masses sit only ~27 MeV apart near 344.1 GeV, fixing the location of the near-threshold excess already seen by ATLAS and CMS.
- The di-gluon channel dominates annihilation (2.57 MeV), while electroweak modes such as Z⁰H (~595 keV) and W⁺W⁻ supply the cleanest experimental tags.
- The system is extremely compact (radii ≲ 0.025 fm), placing it firmly inside the perturbative QCD regime and making the 1S mass a sharp definition of the top-quark mass.
- Because the wave function at the origin is analytic, every partial width scales directly with |Ψ(0)|² and can be updated instantly when α_s or the string tension is refined.
Reading between the lines
- If the single-top width really exceeds the binding energy, the extracted |Ψ(0)|² should be re-interpreted as the residue of a complex pole rather than a normalizable bound-state density; that reinterpretation would systematically shift all quoted widths.
- The same O(λ²) SUSY expansion can be applied without change to the first radial excitations, giving a parameter-free prediction for the 2S–1S splitting that threshold scans could test next.
- Because the calculation is fully analytic in λ, it supplies a transparent uncertainty budget for the top Yukawa extraction once the Z⁰H channel is observed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies the SUSY-QM factorization/Supersymmetric Expansion Algorithm to the Cornell potential for the toponium ground state. Treating the linear term perturbatively through O(λ²), it obtains |Ψ(0)|², a binding energy E_B=−0.9056 GeV, pseudoscalar and vector masses of 344.114 and 344.141 GeV, a 27.2 MeV hyperfine splitting, several characteristic radii, and a broad set of annihilation widths. The reported numbers are compared with earlier potential-model calculations, and the independently decaying top-quark width is also estimated.
Significance. The study is timely given the ATLAS and CMS near-threshold tt̄ measurements and provides explicit, readily checkable potential-model predictions for the spectrum, wave function at the origin, radii, and annihilation channels. A particular strength is the analytic O(λ²) treatment: for the parameters implied by the text, λ is small and the reported numerical chain is internally reproducible. The side-by-side decay tables are also useful. If the quasi-bound-state interpretation and uncertainty treatment are made rigorous, the paper can serve as a compact benchmark for conventional NR potential-model expectations, although the method itself is an application of previously published formalism rather than a new theoretical development.
major comments (4)
- [§II, Eqs. (8)–(12) and §II.A] The numerical potential is not fully specified. Equation (10) defines λ=2σ/(μ²α³), but the manuscript never states σ, λ, or explicitly defines the reduced mass μ used in the calculation. Inverting the quoted E_B with μ=m_t/2 gives approximately λ=0.0157 and σ=0.184 GeV², but readers should not have to reverse-engineer a central input. Please state these values, show the O(λ²) energy formula used for E_B, justify the chosen σ, and give the associated truncation/parameter uncertainty.
- [§II.C and §III] There is an unresolved tension between the bound-state language and the paper's own width estimate. The manuscript obtains Γ_tt̄=2Γ_t=2.9566 GeV, which is more than three times |E_B|=0.9056 GeV and also exceeds the roughly 0.76 GeV 1S–2S spacing implied by the same Coulomb-plus-linear approximation. Thus the object is a broad quasi-bound threshold enhancement, not an isolated narrow 1S level. The tables can still be meaningful as conventional partial annihilation widths evaluated with a frozen real-potential wave function, but the authors must state that approximation, qualify the phrase “tightly bound state,” and discuss the effect of width smearing on the quoted masses and observables.
- [Abstract, Table I, and §III] The precision attached to the central results is not supported by an uncertainty analysis. The masses are quoted to 1 MeV even though the input m_t=172.52±0.33 GeV alone gives an uncertainty of about ±0.66 GeV in 2m_t; α_s(μ), σ, omitted higher-order potential terms, and relativistic corrections add further model dependence. Similar uncertainties propagate into |Ψ(0)|² and every width in Tables II–III. An uncertainty budget is needed, or the results should be rounded and presented explicitly as central values within the stated model.
- [Abstract and §III] The statement that η_t→gg is the dominant decay mode is accurate only within the restricted set of annihilation channels. The manuscript itself finds the independent weak decays of the two constituents to give Γ_tt̄≈2.96 GeV, roughly three orders of magnitude above Γ(η_t→gg)=2.57 MeV. The abstract, discussion, and table captions should consistently say “dominant annihilation channel” and, where branching fractions or experimental relevance are discussed, normalize to the total quasi-bound-state width rather than only to annihilation widths.
minor comments (7)
- [§II.A, Eq. (14)] The hyperfine notation is ambiguous. As written, Eq. (14) multiplies the coefficient by ⟨S₁·S₂⟩, which gives spin-dependent level shifts; the quoted 27.2 MeV is instead the coefficient and hence the vector–pseudoscalar mass difference. Please distinguish δE_S=C|Ψ(0)|²⟨S₁·S₂⟩ from ΔM(Θ_t−η_t)=C|Ψ(0)|².
- [Abstract and §II.B] The abstract states that all three radii lie in 0.015–0.025 fm, but the quoted RMS radius is 0.02657 fm. Either adjust the range or correct the reported value.
- [§II, Eq. (11)] The one-loop coupling formula uses log(e+μ²/Λ²), rather than the more standard log(μ²/Λ²). Although the Euler-constant offset is numerically negligible at μ=m_t, its origin and the scale uncertainty should be explained. It would also help to state explicitly that μ=m_t/2 in the two-body Schrödinger equation.
- [§II.B] The “hyperfine splittings” assigned to r_mp and r_rms (3.72 and 0.88 MeV) are introduced without a defining formula. Since the Fermi contact interaction is controlled by |Ψ(0)|², the meaning and derivation of these radius-based estimates should be given or the statements removed.
- [Tables II–III] Please clarify whether entries such as e⁺e⁻, μ⁺μ⁻, τ⁺τ⁻ and the three neutrino channels in Table III are per flavor or summed over flavors. Adding branching fractions relative to Γ_tt̄ would also make the phenomenological hierarchy clearer.
- [§I] The experimental statements should be more cautious: the collaborations observe a near-threshold cross-section enhancement, while its specific quasi-bound-state interpretation is model dependent. “Effectively proving that the top and anti-top quarks pair up briefly” overstates the inference.
- [General] There are several typographical and wording issues, including “succesfully,” “dimesionless,” “Ricatti,” “Riccati equation” formatting, “event estimated,” and inconsistent spacing around units and particle symbols. The denominator in Eq. (14) should also be typeset unambiguously as 9m_Q².
Circularity Check
No significant circularity: masses, widths, and radii are computed from external PDG inputs, the standard Cornell form, an independent SUSY expansion (Ref. [34]), and literature decay formulas.
full rationale
The derivation chain is feed-forward and non-circular. Inputs are external and independently specified: PDG top mass and couplings (Eq. 11, §II.C constants), the Cornell potential (Eq. 8), the SUSY QM λ-expansion of |Ψ(0)|² and the radial density taken from Napsuciale et al. [34] (no author overlap with the present paper), and standard annihilation-width formulas cited from the broader literature [23–25, 37–39]. From these the paper evaluates E_B, spin-averaged and hyperfine-split masses (Eqs. 13–14), radii via moments of u²(x) (Eqs. 15–16), and partial widths (Eqs. 20–33). None of the reported outputs is algebraically identical to an input by construction, nor is any parameter fitted to toponium data and then re-presented as a prediction of the same data. Mild model dependence (choice of Λ, fixed-scale α_s, implicit σ/λ) is ordinary potential-model practice, not circularity under the stated criteria. Self-citations ([26] Parmar thesis supervised by one co-author) are background and not load-bearing for the uniqueness or form of the results. Score 0 is therefore appropriate.
Assumptions & free parameters
free parameters (3)
- QCD scale Λ in α_s formula =
0.10 GeV
- string tension σ (or normalized λ) =
not reported (implicit via E_B)
- renormalization scale choice μ=m_t for α_s =
μ=m_t=172.52 GeV
assumptions (6)
- domain assumption Quarkonium dynamics are described by the non-relativistic Schrödinger equation with the Cornell potential V(r)=−(4/3)α_s/r + σr.
- domain assumption The linear term may be treated as a perturbation and the superpotential expanded to O(λ²), yielding the quoted |Ψ(0)|² and u²(x) series from the algorithm of Ref. [34].
- domain assumption Spin-averaged mass is 2m_t + E_B and hyperfine splitting is the contact formula ΔE_hf=(32πα_s/9m_Q²)|Ψ(0)|²⟨S1·S2⟩.
- domain assumption Annihilation decay widths factorize into known short-distance coefficients times |ψ(0)|² (and listed NLO α_s factors) as written in Eqs. (20)–(33).
- domain assumption Top and anti-top in the bound state decay independently, so Γ_tt̄≈2Γ_t with the LO three-body weak formula.
- standard math SUSY QM factorization, partner Hamiltonians, and Riccati equation for the superpotential are the correct solution method for the radial problem.
Cite this review
Pith. "Pith review of Mass and Decay Properties of Toponium Using SUSY QM Factorization Method." pith.science (2026). https://pith.science/paper/BC2Y4YN3
@misc{pith2026260723643,
author = {Pith},
title = {Pith review of: Mass and Decay Properties of Toponium Using SUSY QM Factorization Method},
year = {2026},
howpublished = {\url{https://pith.science/paper/BC2Y4YN3}},
note = {Machine review of arXiv:2607.23643}
}
abstract
The Supersymmetric quantum mechanics (SUSY QM) factorization method is applied to study the Cornell potential in the context of toponium $(t\bar{t})$, the heaviest quarkonium system. The method is successfully applied to determine the mass, binding energy, and spin-dependent hyperfine splitting, giving pseudoscalar $(\eta_t)$ and vector $(\Theta_t)$ state masses of 344.114 GeV and 344.141 GeV, respectively. The decay properties of both the pseudoscalar and vector states are also investigated through their respective decay widths, highlighting the crucial role of the bound-state wavefunction. Among the decay channels, the di-gluonic $(\eta_t\rightarrow gg)$ mode is found to be dominant, with an estimated decay width of 2.57 MeV. In addition, the spatial characteristics of toponium are examined by evaluating the mean radius, root-mean-square (RMS) radius, and most probable radius, all of which are found to lie within the range of $0.015\text{--}0.025\,\mathrm{fm}$.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
These transitions arise from the annihilation of the quark–antiquark pair into gluons
Gluonic Decay Width Gluonic decay processes are governed by Quantum Chro- modynamics (QCD) and typically provide the dominant con- tribution in regimes where electromagnetic decay channels are suppressed. These transitions arise from the annihilation of the quark–antiquark pair into gluons. Accordingly, the di-gluonic decay of the pseudoscalar state, (η t...
-
[2]
Photonic Decay Width Even though electromagnetic decay channels, such as pho- tonic decays, are typically suppressed, they remain important 5 because of their clean experimental signatures. Even when such decays proceed through higher-order or indirect mech- anisms, the absence of strong interaction backgrounds makes them particularly valuable for precisi...
-
[3]
Fermionic Decay Width The fermionic decay modes which involves the produc- tion of fermion–antifermion pairs, are particularly sensitive to mass-dependent couplings and thus offer insight into the interplay between bound-state dynamics and electroweak in- teractions [41]. The SUSY QM approach, through its shape- invariant construction and analytic control...
-
[4]
Bosonic decay At higher energy scales, bosonic decay channels, especially those involving massive gauge bosons become increasingly relevant [44, 45]. These bosonic decays of the pseudoscalar quarkonia are computed as [23, 25, 39] Γ(ηt →Z 0γ) = 8παemαZ 27M 2ηt (3−8 sin2 θW )2λ1/2(1, RZ ,0)|ψ(0)| 2, (26) 6 Γ(Θt →Z 0γ) = 8παemαZ(1−R 2 Z) 9M 2 Θt RZ |ψ(0)|2 (...
-
[5]
Navaset al.(Particle Data Group Collaboration), Phys
S. Navaset al.(Particle Data Group Collaboration), Phys. Rev. D110, 030001 (2024)
2024
-
[6]
V . M. Abazovet al., Physical Review D85, 10.1103/phys- revd.85.091104 (2012)
doi:10.1103/phys- 2012
-
[7]
Varma and O
M. Varma and O. K. Baker, Nuclear Physics A1042, 122795 (2024)
2024
-
[8]
R. Francener, V . P. Goncalves, and D. E. Martins, Investigating the exclusive toponium production at the lhc and fcc (2025), arXiv:2502.03295 [hep-ph]
arXiv 2025
Show all 57 references
-
[9]
Fu, Y .-J
J.-H. Fu, Y .-J. Li, H.-M. Yang, Y .-B. Li, Y .-J. Zhang, and C.-P. Shen, Physical Review D111, 10.1103/fqc9-k315 (2025)
2025 doi
-
[10]
I. I. Y . Bigi and H. Krasemann, Z. Phys. C7, 127 (1981)
1981
-
[11]
Aadet al., Journal of High Energy Physics2023, 10.1007/jhep07(2023)141 (2023)
G. Aadet al., Journal of High Energy Physics2023, 10.1007/jhep07(2023)141 (2023)
2023 doi
-
[12]
Tumasyanet al., Journal of High Energy Physics2025, 10.1007/jhep02(2025)064 (2025)
A. Tumasyanet al., Journal of High Energy Physics2025, 10.1007/jhep02(2025)064 (2025)
2025 doi
-
[13]
T. C. Collaboration, Reports on Progress in Physics88, 087801 (2025)
2025
-
[14]
Collaboration, Observation of a cross-section enhancement near thet ¯tproduction threshold in √s= 13tevppcollisions with the atlas detector (2026), arXiv:2601.11780 [hep-ex]
A. Collaboration, Observation of a cross-section enhancement near thet ¯tproduction threshold in √s= 13tevppcollisions with the atlas detector (2026), arXiv:2601.11780 [hep-ex]
2026 arXiv
-
[15]
T. C. S. Group, Cepc conceptual design report: V olume 2 - physics and detector (2018), arXiv:1811.10545 [hep-ex]
2018 arXiv
-
[16]
T. C. S. Group, Cepc technical design report – reference detec- tor (2025), arXiv:2510.05260 [hep-ex]
2025
-
[17]
Agapovet al., Future circular lepton collider fcc-ee: Overview and status (2022), arXiv:2203.08310 [physics.acc- ph]
I. Agapovet al., Future circular lepton collider fcc-ee: Overview and status (2022), arXiv:2203.08310 [physics.acc- ph]
2022 arXiv
-
[18]
Accetturaet al., Interim report for the international muon collider collaboration (imcc) (2025), arXiv:2407.12450 [physics.acc-ph]
C. Accetturaet al., Interim report for the international muon collider collaboration (imcc) (2025), arXiv:2407.12450 [physics.acc-ph]
2025
-
[19]
T. Han, D. Liu, and S. Wang, Physical Review D111, 9 10.1103/physrevd.111.035015 (2025)
2025 doi
-
[20]
Feigenbaum, Phys
J. Feigenbaum, Phys. Rev. D43, 264 (1991)
1991
-
[21]
Sehgal and P
L. Sehgal and P. Zerwas, Nuclear Physics B183, 417 (1981)
1981
-
[22]
M. B. V oloshin, Nucl. Phys. B154, 365 (1979)
1979
-
[23]
Dom `enech, M
G. Dom `enech, M. Goodsell, and C. Wetterich, Journal of High Energy Physics2021, 10.1007/jhep01(2021)180 (2021)
2021 doi
-
[24]
G. F. Giudice, R. Rattazzi, and J. D. Wells, Nuclear Physics B 544, 3–38 (1999)
1999
-
[25]
Wang, T.-F
G.-L. Wang, T.-F. Feng, and Y .-Q. Wang, Phys. Rev. D111, 096016 (2025)
2025
-
[26]
Akbar, I
N. Akbar, I. Asghar, and Z. Ahmad, Advances in High Energy Physics2026, 9454944 (2026), https://onlinelibrary.wiley.com/doi/pdf/10.1155/ahep/9454944
2026
-
[27]
Yang, Chinese Physics C50, 033101 (2026)
2026
-
[28]
Jiang, B.-Q
S.-J. Jiang, B.-Q. Li, G.-Z. Xu, and K.-Y . Liu, Study on toponium: Spectrum and associated processes (2024), arXiv:2412.18527 [hep-ph]
2024 arXiv
-
[29]
Bai, T.-K
Y . Bai, T.-K. Chen, and Y . Yang, Journal of High Energy Physics2026, 10.1007/jhep05(2026)154 (2026)
2026 doi
-
[30]
R. H. Parmar,Study of Potential Problems Using Supersymmet- ric Quantum Mechanics, Ph.D. thesis, Sardar Patel University, Vallabh Vidyanagar, Gujarat, India (2002), supervised by Prof. P. C. Vinodkumar
2002
-
[31]
Schwabl,Quantum Mechanics, 4th ed
F. Schwabl,Quantum Mechanics, 4th ed. (Springer, Berlin, Hei- delberg, 2007)
2007
-
[32]
Cooper, A
F. Cooper, A. Khare, and U. Sukhatme, Physics Reports251, 267–385 (1995)
1995
-
[33]
Tong, Supersymmetric quantum mechanics, Lecture notes available athttps://www.damtp.cam.ac.uk/user/ tong/susy/susyqm.pdf(2007), accessed: 2025-10-09
D. Tong, Supersymmetric quantum mechanics, Lecture notes available athttps://www.damtp.cam.ac.uk/user/ tong/susy/susyqm.pdf(2007), accessed: 2025-10-09
2007
-
[34]
Khare, Proceedings of the Indian Academy of Sciences – Section A49, 41 (1997), accessed: 2025-10-09
A. Khare, Proceedings of the Indian Academy of Sciences – Section A49, 41 (1997), accessed: 2025-10-09
1997
-
[35]
Kulkarni, Resonance: Journal of Science Education8, 28 (2003), accessed: 2026-03-28
A. Kulkarni, Resonance: Journal of Science Education8, 28 (2003), accessed: 2026-03-28
2003
-
[36]
Napsuciale and S
M. Napsuciale and S. Rodr´ıguez, Physics Letters B816, 136218 (2021)
2021
-
[37]
A Alkathiri, L
A. A Alkathiri, L. S. El-Sherif, A. M. Bakry, A. Ismail, and M. Allosh, Physica Scripta99, 105305 (2024)
2024
-
[38]
Napsuciale, S
M. Napsuciale, S. Rodr ´ıguez, and A. Villanueva-Guti ´errez, Physical Review D111, 10.1103/4xk3-x1qy (2025)
2025 doi
-
[39]
K. Saha, S. Upadhaya, and S. Ghosh, inProceedings of the DAE Symposium on Nuclear Physics, V ol. 56 (Department of Atomic Energy, India, 2011) p. E71, accessed: 2026-03-28
2011
-
[40]
Z. R. Najjar and K. Azizi, The European Physical Journal C86, 412 (2026)
2026
-
[41]
Kwong, P
W. Kwong, P. B. Mackenzie, R. Rosenfeld, and J. L. Rosner, Phys. Rev. D37, 3210 (1988)
1988
-
[42]
Fabiano, A
N. Fabiano, A. Grau, and G. Pancheri, Il Nuovo Cimento A (1965-1970)107, 2789 (1994)
1965
-
[43]
Barger, E
V . Barger, E. W. N. Glover, K. Hikasa, W.-Y . Keung, M. G. Olsson, C. J. Suchyta III, and X. R. Tata, Phys. Rev. D35, 3366 (1987)
1987
-
[44]
K ¨all´en,Elementary particle physics(Addison-Wesley, Reading, MA, 1964)
G. K ¨all´en,Elementary particle physics(Addison-Wesley, Reading, MA, 1964)
1964
-
[45]
Bagnaschi, L
E. Bagnaschi, L. Biermann, and M. Spira, SciPost Physics Community Reports 10.21468/scipostphyscommrep.19 (2026)
2026 doi
-
[46]
J. C. Baez and J. Huerta, The algebra of grand unified theories (2010), arXiv:0904.1556 [hep-th]
2010 arXiv
-
[47]
Hern ´andez-Tom´e, C
G. Hern ´andez-Tom´e, C. S. Kim, and G. L. Castro, The Euro- pean Physical Journal C85, 10.1140/epjc/s10052-025-14379-4 (2025)
2025 doi
-
[48]
Kilian, M
W. Kilian, M. Sekulla, T. Ohl, and J. Reuter, Physical Review D91, 10.1103/physrevd.91.096007 (2015)
2015 doi
-
[49]
Bergstrom and G
L. Bergstrom and G. Hulth, Nucl. Phys. B259, 137 (1985), [Erratum: Nucl.Phys.B 276, 744–744 (1986)]
1985
-
[50]
Khachatryanet al., Journal of High Energy Physics2014, 10.1007/jhep09(2014)087 (2014)
V . Khachatryanet al., Journal of High Energy Physics2014, 10.1007/jhep09(2014)087 (2014)
2014 doi
-
[51]
Fu, Y .-J
J.-H. Fu, Y .-J. Zhang, G.-Z. Xu, and K.-Y . Liu, Toponium: Implementation of a toponium model in feynrules (2025), arXiv:2504.12634 [hep-ph]
2025 arXiv
-
[52]
J. Gao, C. S. Li, and H. X. Zhu, Physical Review Letters110, 10.1103/physrevlett.110.042001 (2013)
2013 doi
-
[53]
Zhang, X
S.-W. Zhang, X. Luo, H.-M. Yang, and H.-X. Chen, Universe 11, 334 (2025)
2025
-
[54]
S. C. Benzahra, A. A. Evans, and R. A. Porter, Binding ener- gies and dissociations ofj/ψ,υ, and toponium in quark-gluon plasma (2005), arXiv:hep-ph/0503139 [hep-ph]
2005 arXiv
-
[55]
B. Fuks, K. Hagiwara, K. Ma, and Y .-J. Zheng, Physical Review D104, 10.1103/physrevd.104.034023 (2021)
2021 doi
-
[56]
Bejar, P
I. Bejar, P. Fessia, L. Rossi, O. Bruning, L. Tavian, G. Apolli- nari, and M. Lamont,High-Luminosity Large Hadron Collider (HL-LHC) Technical Design Report V . 0.1(2017)
2017
-
[57]
Kats and M
Y . Kats and M. D. Schwartz, Journal of High Energy Physics 2010, 10.1007/jhep04(2010)016 (2010)
2010 doi
Reviewed July 30, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.