{"id":"7f37eec9-c5e6-44c7-8fa4-e83a6d5d1ed1","arxiv_id":"2411.17955","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Toponium masses from the Salpeter equation with Cornell and Coulomb potentials lie near 343.6-345.3 GeV with zero spin splittings, and predicted decay widths to photons, gluons, and leptons differ by up to a factor of two between the two potentials.","lead":"Using the Bethe-Salpeter equation with two different force laws, the authors computed masses and wave functions for toponium, a hypothetical bound state of a top quark and top antiquark. The results give concrete masses near 344 GeV and decay rates that could be tested at the LHC if toponium exists.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Decay-width formulas (Eqs. 24 and 27) are dimensionally inconsistent with the stated normalization (Eq. 12); as printed they cannot yield the keV/MeV values in Table III.","rationale":"The most load-bearing flaw is not the physical criticism (width vs binding energy), which the authors explicitly acknowledge, but the internal inconsistency of the printed decay formulas. A dimensionally wrong formula cannot be a minor typo because the numerical results are quoted to three significant figures; the code must be using a different (unstated) normalization or formula. This places the central decay predictions in Table III in doubt. The mass spectrum and degeneracy claims may survive, but the paper's headline decay widths are unsupported unless corrected. Therefore the verdict should remain CONDITIONAL, pending a corrected derivation and confirmation of the numbers.","tokens_in":10147,"tokens_out":20771,"duration_ms":171583,"concrete_test":"Compute the trace in Eq. 23 and the two-body phase space to derive the exact width formula from the Salpeter amplitude, keeping the normalization of Eq. 12. Compare the resulting M-dependence with Eq. 24; also check the gg replacement. If the correct formula is not proportional to M^-3 times the squared bracket, Table III is unsupported and must be recalculated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Eq. 12 implies φ_S has dimension GeV^-2 because d⃗q/(2π)^3 contributes GeV^3 and 4m_t M/ω contributes GeV. In Eq. 24 the integral ∫ d⃗q/(2π)^3 φ_S [1/(p1−k1)^2 + 1/(p1−k2)^2] therefore has dimension GeV^3 · GeV^-2 · GeV^-2 = GeV^-1, so its square is GeV^-2. The prefactor 12πα^2e_q^4/M^3 adds GeV^-3, giving Γ a total dimension GeV^-5 instead of GeV. The same issue applies to Eq. 27, whose prefactor is also 8/(3π)α_s^2/M^3, whereas the stated replacement e_q^4α^2→(2/9)α_s^2 in Eq. 24 would produce 8π/3 α_s^2/M^3 — an additional π^2 discrepancy. These errors mean the quoted widths (7.56 keV, 15.9 keV, etc.) cannot be reproduced from the printed equations unless φ_S or the prefactors are defined differently than stated.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper solves the instantaneous Bethe-Salpeter equation for toponium using the Cornell potential and, for comparison, a pure Coulomb potential, with mt = 172.7 GeV. It reports mass spectra for S-, P-, D-, F-, and G-wave states, finds that singlet-triplet and triplet spin splittings vanish exactly, and provides nonrelativistic wave functions for S-, P-, and D-wave toponia. It then uses these wave functions to compute the two-photon and two-gluon widths of the pseudoscalar eta_t and the dilepton width of the vector Theta. The main numerical outputs are masses in the 343-345 GeV range, with the 1S state around 343.6 GeV, and ground-state widths that differ by roughly a factor of two between the Cornell and Coulomb potentials.","tokens_in":10360,"tokens_out":13577,"duration_ms":130505,"significance":"If the central results are correct, the paper provides a useful complementary study of toponium spectroscopy and decays in a regime where the top quark is extremely heavy and the Coulomb interaction dominates. The comparison between Cornell and Coulomb potentials is transparent, and the Coulomb-potential results for the 1S mass and the decay widths are broadly consistent with earlier threshold and NRQCD estimates. The explicit normalization conditions and wave function ansaetze are helpful for reproducibility. However, the paper does not provide code or numerical details of the Salpeter solver, and the printed decay formulas have a serious dimensional inconsistency that prevents the reader from reproducing the quoted widths. The significance of the decay predictions is therefore conditional on correcting the formulas and re-deriving the numbers.","major_comments":[{"comment":"Equation (24) is dimensionally inconsistent as printed. From the normalization condition Eq. (12), phi_S has dimension GeV^-2. The integral in Eq. (24) then has dimension GeV^3 * GeV^-2 * GeV^-2 = GeV^-1, so its square is GeV^-2. The prefactor 12 pi alpha^2 e_q^4 / M^3 has dimension GeV^-3, making the right-hand side dimension GeV^-5 rather than GeV. As written, Eq. (24) cannot produce the keV values in Table III. Equation (27) inherits the same problem. In addition, the replacement e_q^4 alpha^2 -> (2/9) alpha_s^2 stated in the text would convert the prefactor of Eq. (24) to 8 pi/3 alpha_s^2 / M^3, not 8/(3 pi) alpha_s^2 / M^3 as printed in Eq. (27); the printed coefficient is smaller by a factor pi^2. The authors should provide a full derivation of the decay formulas, correct the prefactors and dimensions, and recompute or verify every entry in Table III.","section":"III, Eq. (24) and Eq. (27)"},{"comment":"The calculation treats toponium as a stable bound state: the Salpeter equation is solved without any top-quark width term, although Sec. I quotes the top width as 1.42 GeV. The computed 1S binding energy is 2 mt - 343.62 GeV = 1.78 GeV, comparable to the top width, and the 2S-1S splitting is only about 0.97 GeV. A state whose constituent decays on a timescale comparable to the binding dynamics cannot be described as a stationary eigenstate without additional justification or a complex-mass/optical-potential treatment. This issue is load-bearing for the entire mass spectrum and for the decay widths, which would in practice be smeared or suppressed. The paper should either include the width in the formalism or clearly state that all results are for a hypothetical stable top, with a quantitative discussion of the expected effect of the finite width.","section":"Sec. I and Sec. IV"},{"comment":"The exact degeneracies M(n1S0)=M(n3S1), M(n3P0)=M(n3P1)=M(n3P2)=M(n1P1), and their analogues in Eqs. (5)-(9) are not dynamical predictions of the model; they follow by construction because the input Cornell and Coulomb potentials contain no spin-dependent or spin-orbit terms. In QCD, spin splittings are nonzero, albeit small, with a parametric size of order alpha_s^4 m_t. The paper should present these equalities as approximations of a central-potential model and estimate the size of the omitted spin-dependent corrections, rather than stating that the splittings vanish. This is important because the abstract and Sec. IV present the vanishing splittings as one of the main findings.","section":"Sec. II, Eqs. (5)-(9)"}],"minor_comments":[{"comment":"The text refers to the 'color-triplet ground state toponium' when discussing Theta -> l+ l-; the vector toponium is a color singlet. This appears to be a typo and should be corrected.","section":"Sec. III, after Eq. (31)"},{"comment":"The figure legends label the curves phi1_S(q) and phi2_S(q), but the text defines a single radial function phi_S(q). The relation between phi1, phi2 and phi_S should be stated explicitly.","section":"Figs. 1 and 2"},{"comment":"There are several typographical errors, including 'exce edingly' in the abstract, 'dileton' instead of 'dilepton' in Sec. III, and inconsistent formatting such as 'm t¯t'. A careful proofreading pass is needed.","section":"Throughout"},{"comment":"The masses and widths are quoted without uncertainties or sensitivity estimates. Given that the Cornell and Coulomb potentials differ by roughly a factor of two in the ground-state decay widths, a discussion of the dependence on lambda, Lambda, and the scale of alpha_s would help the reader assess the robustness of the predictions.","section":"Tables I-III"}],"recommendation":"major_revision","confidential_remarks":"The dimensional inconsistency in Eq. (24) is the most serious issue. If the authors cannot provide a corrected derivation and reproducible numerical procedure, the decay-width claims should not be accepted as reliable. The mass spectrum itself may be salvageable, but the paper as it stands does not meet the standard for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a workmanlike application of the authors' established Salpeter solver to toponium with mt = 172.7 GeV. The mass spectrum (1S near 343.6 GeV, nearly degenerate multiplets) and the factor-of-two difference between Cornell and Coulomb decay widths are worth knowing, especially since the Coulomb results agree with earlier estimates. The paper is honest about prior literature and the experimental hints.\n\nThe serious problem is that the decay-width formulas as printed cannot be right. With the normalization in Eq. (12), φ_S has dimension GeV^-2. The integral in Eq. (24) is then GeV^-1, the prefactor GeV^-3, so Γ comes out GeV^-5, not GeV. Eq. (27) has the same issue plus a π^2 error: the stated replacement e_q^4 α^2 → (2/9) α_s^2 in Eq. (24) would give 8π/3 α_s^2/M^3, not 8/(3π) α_s^2/M^3. So the quoted 7.56 keV, 15.9 keV, etc., cannot be reproduced from the equations as written. This looks like typos in prefactors or an unstated different wave-function normalization, but it has to be fixed before the numbers are usable.\n\nOther soft spots: the exact degeneracies M(n1S0)=M(n3S1), etc., are stated as equalities without error bars; they are approximate, and while plausible for a heavy quark onium, they deserve a caveat. The top width (1.4 GeV) is comparable to the binding energy (~1.8 GeV), so treating toponium as a stable bound state is a real approximation; the paper acknowledges this but does not estimate its effect. The strong potential dependence of the widths is both a feature (discriminating) and a warning about model sensitivity.\n\nWho this is for: people actively searching for toponium or working on t-tbar threshold effects. The mass and width predictions are a useful benchmark, but only after the formulas are corrected. The central mass calculation is plausible within the model, and the citation pattern is fine.\n\nRecommendation: it deserves a serious referee, but not as-is. I would send to peer review with a clear request to fix the decay-width equations, add uncertainties to the mass degeneracies, and comment on the width-to-binding ratio.\n\nBest,\n[You]","headline":"Useful toponium potential-model calculation, but the printed decay-width equations don't match the quoted numbers—fix before publish.","tokens_in":10892,"tokens_out":4795,"would_cite":false,"duration_ms":45395,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Solving the Salpeter equation for toponium yields fully degenerate spin multiplets.","keywords":["toponium","Bethe-Salpeter equation","Salpeter equation","Cornell potential","Coulomb potential","mass spectrum","decay width","top quark"],"falsifier":"A high-statistics scan of the $t\\bar t$ invariant mass near 343–345 GeV that resolves a narrow $1S$ resonance and measures the $^1S_0$–$^3S_1$ mass separation; any spin splitting above roughly 1 MeV would contradict the claimed degeneracy. Alternatively, a measured $\\Gamma(\\eta_t\\to\\gamma\\gamma)$ well outside the predicted 7.56–15.9 keV band, or an observed excited-state mass pattern different from the degenerate tables, would falsify the wave-function and potential choices.","tokens_in":9933,"feed_emoji":"⚛️","tokens_out":2869,"duration_ms":27436,"temperature":0.7,"pith_summary":"The paper claims that toponium, a top-antitop bound state, can be described by solving the instantaneous Bethe-Salpeter equation with either the Cornell potential or the pure Coulomb potential at a top-quark mass of 172.7 GeV. Because the top quark is so heavy, the paper finds that all spin-dependent splittings vanish: the pseudoscalar and vector S-wave states share one mass, and the P-, D-, F-, and G-wave triplet and singlet states are each pairwise degenerate. The paper then uses the resulting wave functions to predict the two-photon, two-gluon, and dilepton decay widths of the ground-state toponia, finding values that differ by about a factor of two depending on whether the linear confinement term is included. These predictions matter because recent LHC data near the top-pair threshold show hints that a pseudoscalar toponium state may be contributing to the observed events.","feed_headline":"Toponium masses computed: spin splittings vanish","feed_subtitle":"New Salpeter solutions give degenerate S, P, D multiplets and concrete two-photon, two-gluon, and dilepton widths.","key_machinery":"The machinery is the instantaneous Bethe-Salpeter (Salpeter) equation, solved with a specified relativistic wave-function ansatz for each $J^{PC}$ channel and a running strong coupling $\\alpha_s(\\vec q)$ with $\\Lambda = 0.10$ GeV. The Cornell potential $V_{\\text{Cornell}}(r)=\\lambda r - \\frac{4}{3}\\frac{\\alpha_s}{r}$ supplies the confining and Coulomb interaction, while the Coulomb-only case drops the linear term. The numerical solution yields mass eigenvalues and radial wave functions $\\varphi_S$, $\\varphi_P$, $\\varphi_D$; because the top quark is so heavy, the paper reduces the relativistic wave functions to nonrelativistic forms and uses them, together with standard triangle and current matrix-element formulas, to compute the decay widths.","core_discovery":"The central claim is that, for toponium with $m_t = 172.7$ GeV, the mass spectrum is spin-degenerate: $M(n^1S_0)=M(n^3S_1)$, $M(n^3P_0)=M(n^3P_1)=M(n^3P_2)=M(n^1P_1)$, and analogous degeneracies for D, F, and G waves. Solving the Salpeter equation with the Cornell potential ($\\lambda = 0.18$ GeV$^2$) gives a $1S$ mass of 343.62 GeV, while the Coulomb-only potential gives 343.59 GeV; excited states lie between roughly 344.3 and 345.3 GeV, all below the $2m_t$ threshold. The paper further claims that the decay widths of the ground states are $\\Gamma(\\eta_t\\to\\gamma\\gamma)=7.56$ keV (Cornell) or 15.9 keV (Coulomb), $\\Gamma(\\eta_t\\to gg)=1.69$ MeV or 3.54 MeV, and $\\Gamma(\\Theta\\to\\ell^+\\ell^-)=6.09$ keV or 12.9 keV. The large spread between the Cornell and Coulomb results is attributed to the linear potential substantially reshaping the wave function, and the paper argues that future experiments can discriminate between the two potentials.","pith_inferences":["The paper treats toponium as stable, but the top-quark width (1.42 GeV) is comparable to the binding energy of about 1.8 GeV; a natural extension is to solve the Salpeter equation with a complex mass or width term and see how much the masses and widths shift, or whether the bound state dissolves.","If the spin degeneracy holds, the production amplitudes for $^1S_0$ and $^3S_1$ toponia near threshold become essentially equal at leading order, which could be used as a consistency check in $t\\bar t$ threshold scans.","The wave functions provided could be plugged into production cross-section simulations to estimate whether the predicted narrow $1S$ state is visible above the $t\\bar t$ continuum at the LHC and future colliders.","A decisive test would be a dedicated search for the diphoton final state from a 343.6 GeV resonance; a measured width outside the 7.6–15.9 keV band would force a change in the potential or the bound-state formalism."],"forward_implications":["If the spectrum is correct, only the $1S$ toponium state lies near 343.6 GeV, matching the 343–344 GeV peak region reported by threshold-production studies of $t\\bar t$ pairs.","The claimed spin degeneracy means that spin-dependent relativistic corrections can be dropped in toponium calculations, simplifying future wave-function and production studies.","The predicted two-photon, two-gluon, and dilepton widths provide concrete targets: a ground-state pseudoscalar toponium near 343.6 GeV should have a two-photon partial width of roughly 7.6–16 keV depending on the potential.","Because the Cornell and Coulomb results differ by nearly a factor of two for the ground-state decays, a sufficiently precise measurement of any one of these channels would distinguish which potential better describes the toponium system."],"supporting_citations":[{"why":"Introduces the Bethe-Salpeter equation that defines the relativistic bound-state wave function.","marker":"[39]"},{"why":"Derives the instantaneous Salpeter equation from the Bethe-Salpeter equation, the central equation solved here.","marker":"[40]"},{"why":"Provides the pseudoscalar wave-function ansatz and the method for solving the Salpeter equation for mesons.","marker":"[41]"},{"why":"Supplies the vector wave-function ansatz and the solving methodology for vector mesons.","marker":"[42]"},{"why":"Defines the Cornell potential $\\lambda r - \\frac{4}{3}\\alpha_s/r$ used for the main spectrum calculation.","marker":"[43]"},{"why":"Gives the running strong-coupling form in momentum space used in both potentials.","marker":"[44]"},{"why":"Provides a Coulomb-potential estimate of the diphoton width, 16.865 keV, used as the comparison point for the paper's Coulomb result of 15.9 keV.","marker":"[46]"},{"why":"Supplies the nonrelativistic QCD prediction of 22.6 keV for the diphoton width, which the paper compares against its own running-coupling result.","marker":"[52]"},{"why":"Reports a two-electron partial width of $13\\pm1$ keV that the paper cites as excellent agreement with its Coulomb-potential dilepton width.","marker":"[25]"}],"fun_headline_variants":["Toponium spin splittings vanish: all S, P, D states degenerate","Toponium masses solved: spin degeneracy across all multiplets","Cornell vs Coulomb toponium: decay widths to discriminate potentials","Toponium S, P, D states degenerate: masses from Salpeter equation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The computation treats toponium as a stable bound state and never includes the top quark's decay width, even though that width (1.42 GeV) is comparable to the computed binding energy; if the top quark decays before a bound state forms, the whole spectrum and width predictions lose physical meaning.","fun_headline_variants_meta":{"raw":{"variants":["Toponium spin splittings vanish: all S, P, D states degenerate","Toponium masses solved: spin degeneracy across all multiplets","Cornell vs Coulomb toponium: decay widths to discriminate potentials","Toponium S, P, D states degenerate: masses from Salpeter equation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000593,"raw_usage":{"total_tokens":2780,"prompt_tokens":949,"completion_tokens":1831,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":1753}},"tokens_in":565,"tokens_out":1831,"duration_ms":11360,"temperature":1.0,"reasoning_tokens":1753,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:40:07.052888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics scan of the $t\\bar t$ invariant mass near 343–345 GeV that resolves a narrow $1S$ resonance and measures the $^1S_0$–$^3S_1$ mass separation; any spin splitting above roughly 1 MeV would contradict the claimed degeneracy. Alternatively, a measured $\\Gamma(\\eta_t\\to\\gamma\\gamma)$ well outside the predicted 7.56–15.9 keV band, or an observed excited-state mass pattern different from the degenerate tables, would falsify the wave-function and potential choices.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Bethe-Salpeter equation that defines the relativistic bound-state wave function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the instantaneous Salpeter equation from the Bethe-Salpeter equation, the central equation solved here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the pseudoscalar wave-function ansatz and the method for solving the Salpeter equation for mesons."},{"cited_title":"Wang, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the vector wave-function ansatz and the solving methodology for vector mesons."},{"cited_title":"Eichten, K","cited_arxiv_id":null,"evidence_quote":"Defines the Cornell potential $\\lambda r - \\frac{4}{3}\\alpha_s/r$ used for the main spectrum calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides a Coulomb-potential estimate of the diphoton width, 16.865 keV, used as the comparison point for the paper's Coulomb result of 15.9 keV."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports a two-electron partial width of $13\\pm1$ keV that the paper cites as excellent agreement with its Coulomb-potential dilepton width."}],"review_version":1}