Pith. sign in

REVIEW 1 major objections 4 minor 33 references

Near-threshold top-pair excesses at the LHC are mostly ordinary higher-order QCD, not a free-standing toponium signal.

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 · grok-4.5

2026-07-13 02:17 UTC pith:Y3KKOGSM

load-bearing objection Public NLO+PS generators plus clean ATLAS-bin numbers that settle the size of threshold enhancement and the loosely defined toponium piece. the 1 major comments →

arxiv 2607.09539 v1 pith:Y3KKOGSM submitted 2026-07-10 hep-ph

Top-Antitop Production and Decay at Threshold at the LHC in QCD Perturbation Theory

classification hep-ph
keywords top-antitop productionthreshold resummationCoulomb correctionstoponiumNLO+PS generatorsfinite top widthLHC phenomenology
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

When a top and antitop are produced with mass just above twice the top mass, Coulomb-like soft-gluon exchanges generate large corrections that can be summed to all orders. The paper builds three public NLO event generators that include those corrections and asks how large they really are once standard fixed-order QCD is already accounted for, how much of them can be labelled “toponium,” and whether the top’s finite lifetime matters. In the narrow mass window used by ATLAS the extra enhancement over ordinary NLO is only a few picobarns, the loosely defined bound-state piece is smaller still, and truncating the series at third order already captures nearly everything. Finite-width effects stay modest once the experimental resolution is coarser than the top width. The generators therefore give experiments a clean way to re-interpret the recent excess claims without double-counting or over-attributing the signal to a quasi-bound state.

Core claim

In the ATLAS threshold bin the full resummed Coulomb enhancement over a baseline NLO calculation is of order four picobarns; the contribution that can be ascribed to would-be toponium bound states is only a few picobarns; and the first three orders of the expansion already contain essentially the whole effect. Finite-width corrections remain small for any mass resolution much larger than the top width.

What carries the argument

Three NLO+PS generators (thr1, thr2, bb4l) that multiply the Born matrix element by the all-order Sommerfeld/Coulomb factor (including the discrete bound-state poles) while subtracting the pieces already present at NLO, with an approximate but inclusively consistent treatment of off-shell top virtualities.

Load-bearing premise

The Monte-Carlo treatment of finite-width effects uses a simplified recipe that matches the known inclusive formula after integration over top virtualities, but does not implement the fully differential off-shell formula derived in the appendix.

What would settle it

A high-resolution measurement of the top-pair invariant-mass spectrum below the nominal 2mt threshold (or a future e+e- threshold scan) that shows a shape or size incompatible with the thr2/bb4l predictions once the approximate off-shell recipe is replaced by the exact resolvent formula of Appendix A.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 4 minor

Summary. The paper constructs three public NLO+PS Monte Carlo generators (thr1, thr2, bb4l) that incorporate all-order Coulomb-enhanced threshold corrections (Sommerfeld factor plus bound-state poles) for ttbar production and decay near the 2mt threshold at the LHC. Using these tools the authors quantify, for the ATLAS bin (mtt < 350 GeV, p* < 50 GeV) and for coarser mtt cuts, the size of the enhancement beyond ordinary NLO/NNLO, the loosely defined bound-state contribution, and the residual impact of finite top width. They conclude that O(αs^{3}) already captures the bulk of the effect (~4 pb over baseline NLO), that the bound-state piece is a few pb, and that finite-width shifts remain at the 1–2 pb level once experimental resolution is coarser than Γt.

Significance. The work supplies concrete, publicly available generators that allow experimental analyses to include threshold resummation without double-counting NLO pieces and without treating toponium as an extra free-floating resonance. The order-by-order tables, the three-generator cross-check, the semi-analytic phase-space test (Fig. 7), and the comparison with an independent NLO-NRQCD calculation constitute a solid, falsifiable clarification of the ATLAS/CMS excess. The public codes and the explicit separation of double-top / single-top / no-top samples are particularly valuable for future LHC studies.

major comments (1)
  1. [Appendix A, Eq. (A.19)] Appendix A derives the fully differential resolvent formula (A.19) at fixed top and antitop virtualities, yet the thr2 and bb4l generators implement only an approximate recipe that recovers the inclusive Fadin–Khoze spectral density after integration over virtualities. While the authors argue (and the numerical checks support) that this is adequate for LHC-resolution observables, a short quantitative estimate of the residual error on the double-top sample below the nominal threshold would strengthen the claim that the approximation does not affect the central LHC conclusions.
minor comments (4)
  1. [Abstract / §1] In the abstract and introduction the phrase “enhanced t tbar production near threshold in the pseudoscalar channel” is repeated; a single, precise statement of the experimental claim would improve readability.
  2. [Table 3, §5] Table 3 shows a large single-top contamination in the bb4l NLO sample; a brief remark on the jet-radius dependence of the dt/st classification would help users of the generator.
  3. [§2.2] The regularisation of the δ-function contribution (Eq. 2.7) and the polynomial P(v) used for the (αs/v)^{4} term are technical but essential; a short sentence confirming numerical stability under variation of Ecut and η would be useful.
  4. A few typographical slips remain (“as as a”, “to toponium”, missing spaces around equations). A final proof-reading pass is recommended.

Circularity Check

0 steps flagged

No significant circularity: threshold resummation follows standard Fadin–Khoze NRQCD; generators implement known formulae without fitting free parameters to the ATLAS/CMS excess.

full rationale

The paper’s central claims (size of Coulomb-enhanced corrections beyond NLO/NNLO, size of the loosely defined bound-state piece, and smallness of finite-width effects for LHC-resolution bins) are obtained by implementing the textbook non-relativistic spectral density of Fadin–Khoze (eqs. 2.1–2.6, A.27–A.33) inside three independent NLO+PS generators. No free parameters are fitted to the experimental excess; the only free choices are conventional scale factors that are varied and tabulated. Self-citations to the authors’ earlier NRR paper supply intermediate formulae that are re-derived or extended here (higher-order terms, finite-width recipe of Appendix A). The approximate finite-width mapping (A.24) is acknowledged as incomplete relative to the fully differential resolvent (A.19), but that incompleteness is a technical limitation, not a circular reduction of a prediction to an input. External comparison with an independent NLO-NRQCD calculation (ref. [18]) further anchors the results. Consequently the derivation chain is self-contained against external benchmarks and exhibits no circular step that forces the quoted cross-section numbers by construction.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

The calculation inherits standard QCD, the Coulomb potential, and the classic Fadin–Khoze spectral density. Free choices are limited to technical regularisation parameters, the scale prescription for αs in the Coulomb factor, colour-channel fractions, and the approximate finite-width recipe. No new dynamical entities are postulated.

free parameters (4)
  • E_cut (delta-function regularisation) = 1 GeV (central)
    Width of the polynomial approximation to δ(E) used for bound-state terms in thr1; varied by a factor of two around 1 GeV to check stability.
  • Coulomb αs scale factor = central = 1
    Heuristic scale √(m √(E²+Γt²)) multiplied by a user factor; varied by ½ and 2 in tables. Not fitted to data.
  • Colour singlet/octet fractions in bb4l = 2/7, 5/7
    Fixed 2/7 and 5/7 ansatz for Born colour decomposition when full colour information is not re-projected.
  • Top virtuality window for dt/st/nt classification = 15 GeV
    ±15 GeV cut used to separate double-top, single-top and no-top samples; analysis choice, not a fit.
axioms (4)
  • domain assumption Non-relativistic Coulomb resummation (Sommerfeld factor + bound-state poles) correctly captures the leading (αs/v)^n threshold series for colour-singlet and colour-octet channels.
    Taken from Fadin–Khoze and used throughout §§2–3 and Appendix A.
  • ad hoc to paper After integration over top and antitop virtualities the approximate MC recipe reproduces the inclusive finite-width spectral density of refs. [3,4].
    Stated in Appendix A; full differential formula (A.19) is derived but not implemented.
  • domain assumption Landau–Yang suppression of the spin-triplet gg→tt amplitude near threshold remains valid for the off-shell kinematics used here.
    Used to explain the dominance of the spin-singlet component (Introduction and §2.1).
  • domain assumption Standard POWHEG NLO matching and Les-Houches shower interface remain valid after the Born reweighting by the threshold factor.
    Underlying framework of all three generators.

pith-pipeline@v1.1.0-grok45 · 35115 in / 2982 out tokens · 36531 ms · 2026-07-13T02:17:00.308671+00:00 · methodology

0 comments
read the original abstract

In this work we consider the production of a top-antitop pair at the LHC when the mass of the pair is relatively near to the nominal threshold, that is to say to twice the top pole mass. In this regime, enhanced perturbative corrections arise that can be computed to all orders in perturbation theory. We present three generators of the NLO+PS kind (Next-to-Leading-Order that can be interfaced to parton showers) that include these threshold enhanced effects. Using these generators we address the following questions: what is the size of enhanced non-relativistic effects that are not already present in the well known NLO and NNLO perturbative results; what is the size of the contribution from these effects that can be loosely attributed to toponium production; and to what extent the finite width of the top quark affects threshold enhanced corrections. Our generators are relevant for the recent observation of enhanced $t{\bar t}$ production near threshold in the pseudoscalar channel by the ATLAS and CMS collaborations.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 22 linked inside Pith

  1. [1]

    V. S. Fadin and V. A. Khoze,Threshold Behavior of Heavy Top Production in e+ e- Collisions,JETP Lett.46(1987) 525–529

  2. [2]

    V. S. Fadin, V. A. Khoze, and T. Sjostrand,On the Threshold Behavior of Heavy Top Production,Z. Phys. C48(1990) 613–622

  3. [3]

    B. Fuks, K. Hagiwara, K. Ma, and Y.-J. Zheng,Simulating toponium formation signals at the LHC,Eur. Phys. J. C85(2025), no. 2 157, [arXiv:2411.18962]

  4. [4]

    B. Fuks, K. Hagiwara, K. Ma, and Y.-J. Zheng,Signatures of toponium formation in LHC run 2 data,Phys. Rev. D104(2021), no. 3 034023, [arXiv:2102.11281]

  5. [5]

    Sumino and H

    Y. Sumino and H. Yokoya,Bound-state effects on kinematical distributions of top quarks at hadron colliders,JHEP09(2010) 034, [arXiv:1007.0075]. [Erratum: JHEP 06, 037 (2016)]

  6. [6]

    L. D. Landau,On the angular momentum of a system of two photons,Dokl. Akad. Nauk SSSR60(1948), no. 2 207–209

  7. [7]

    Yang,Selection Rules for the Dematerialization of a Particle Into Two Photons,Phys

    C.-N. Yang,Selection Rules for the Dematerialization of a Particle Into Two Photons,Phys. Rev.77(1950) 242–245

  8. [8]

    M. A. Braun,Positronium singularities in quantum electrodynamics and perturbation theory, Zh. Eksp. Teor. Fiz.54(1968) 1220–1227

  9. [9]

    Melnikov, A

    K. Melnikov, A. Vainshtein, and M. Voloshin,Remarks on the effect of bound states and threshold in g-2,Phys. Rev. D90(2014), no. 1 017301, [arXiv:1402.5690]

  10. [10]

    Beneke and P

    M. Beneke and P. Ruiz-Femenia,Threshold singularities, dispersion relations and fixed-order perturbative calculations,JHEP08(2016) 145, [arXiv:1606.02434]

  11. [11]

    Nason, E

    P. Nason, E. Re, and L. Rottoli,Spin correlations int tproduction and decay at the LHC in QCD perturbation theory,JHEP10(2025) 149, [arXiv:2505.00096]

  12. [12]

    Maltoni, C

    F. Maltoni, C. Severi, S. Tentori, and E. Vryonidou,Quantum detection of new physics in top-quark pair production at the LHC,JHEP03(2024) 099, [arXiv:2401.08751]

  13. [13]

    Frixione, P

    S. Frixione, P. Nason, and G. Ridolfi,A Positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction,JHEP09(2007) 126, [arXiv:0707.3088]. – 36 –

  14. [14]

    Frixione, E

    S. Frixione, E. Laenen, P. Motylinski, and B. R. Webber,Angular correlations of lepton pairs from vector boson and top quark decays in Monte Carlo simulations,JHEP04(2007) 081, [hep-ph/0702198]

  15. [15]

    Jeˇ zo, J

    T. Jeˇ zo, J. M. Lindert, P. Nason, C. Oleari, and S. Pozzorini,An NLO+PS generator fort ¯t andW tproduction and decay including non-resonant and interference effects,Eur. Phys. J. C76(2016), no. 12 691, [arXiv:1607.04538]

  16. [16]

    M. V. Garzelli, G. Limatola, S. O. Moch, M. Steinhauser, and O. Zenaiev,Updated predictions for toponium production at the LHC,Phys. Lett. B866(2025) 139532, [arXiv:2412.16685]

  17. [17]

    Sommerfeld,Atombau und spektrallinien,Bd.2(Vieweg, Braunschweig, 1939)

    A. Sommerfeld,Atombau und spektrallinien,Bd.2(Vieweg, Braunschweig, 1939)

  18. [18]

    SakharovJETP18(1948) 631

    A. SakharovJETP18(1948) 631

  19. [19]

    Alwall et al.,A standard format for les houches event files,Comput

    J. Alwall et al.,A standard format for les houches event files,Comput. Phys. Commun.176 (2007) 300–304, [hep-ph/0609017]

  20. [20]

    Mazzitelli, P

    J. Mazzitelli, P. F. Monni, P. Nason, E. Re, M. Wiesemann, and G. Zanderighi,Top-pair production at the LHC with MINNLO P S,JHEP04(2022) 079, [arXiv:2112.12135]

  21. [21]

    Mazzitelli, P

    J. Mazzitelli, P. F. Monni, P. Nason, E. Re, M. Wiesemann, and G. Zanderighi, Next-to-Next-to-Leading Order Event Generation for Top-Quark Pair Production,Phys. Rev. Lett.127(2021), no. 6 062001, [arXiv:2012.14267]. [24]NNPDFCollaboration, R. D. Ball et al.,Parton distributions for the LHC Run II,JHEP 04(2015) 040, [arXiv:1410.8849]

  22. [22]

    Buckley, J

    A. Buckley, J. Ferrando, S. Lloyd, K. Nordstr¨ om, B. Page, M. R¨ ufenacht, M. Sch¨ onherr, and G. Watt,LHAPDF6: parton density access in the LHC precision era,Eur. Phys. J. C75 (2015) 132, [arXiv:1412.7420]

  23. [23]

    Hou et al.,New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC,Phys

    T.-J. Hou et al.,New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC,Phys. Rev. D103(2021), no. 1 014013, [arXiv:1912.10053]

  24. [24]

    Bailey, T

    S. Bailey, T. Cridge, L. A. Harland-Lang, A. D. Martin, and R. S. Thorne,Parton distributions from LHC, HERA, Tevatron and fixed target data: MSHT20 PDFs,Eur. Phys. J. C81(2021), no. 4 341, [arXiv:2012.04684]

  25. [25]

    Fischler,Quark - anti-Quark Potential in QCD,Nucl

    W. Fischler,Quark - anti-Quark Potential in QCD,Nucl. Phys. B129(1977) 157–174

  26. [26]

    Billoire,How Heavy Must Be Quarks in Order to Build Coulombic q anti-q Bound States, Phys

    A. Billoire,How Heavy Must Be Quarks in Order to Build Coulombic q anti-q Bound States, Phys. Lett. B92(1980) 343–347

  27. [27]

    Cacciari, G

    M. Cacciari, G. P. Salam, and G. Soyez,The anti-k t jet clustering algorithm,JHEP04 (2008) 063, [arXiv:0802.1189]

  28. [28]

    M. V. Garzelli, G. Limatola, S.-O. Moch, M. Steinhauser, and O. Zenaiev,Threshold Top-Quark Pair-Production: Cross Sections and Key Uncertainties,arXiv:2604.09485

  29. [29]

    M. J. Strassler and M. E. Peskin,The Heavy top quark threshold: QCD and the Higgs,Phys. Rev. D43(1991) 1500–1514

  30. [30]

    A. H. Hoang and T. Teubner,Top quark pair production at threshold: Complete next-to-next-to-leading order relativistic corrections,Phys. Rev. D58(1998) 114023, [hep-ph/9801397]

  31. [31]

    A. H. Hoang and M. Stahlhofen,The Top-Antitop Threshold at the ILC: NNLL QCD Uncertainties,JHEP05(2014) 121, [arXiv:1309.6323]. – 37 –

  32. [32]

    Beneke, Y

    M. Beneke, Y. Kiyo, P. Marquard, A. Penin, J. Piclum, and M. Steinhauser, Next-to-Next-to-Next-to-Leading Order QCD Prediction for the Top AntitopS-Wave Pair Production Cross Section Near Threshold ine +e− Annihilation,Phys. Rev. Lett.115(2015), no. 19 192001, [arXiv:1506.06864]

  33. [33]

    Sj¨ ostrand, V

    T. Sj¨ ostrand, V. A. Khoze, and C. T. Preuss,Top Pair Threshold Revisited, arXiv:2605.19546. – 38 –