REVIEW 3 major objections 4 minor 34 references
Accounting for all-order Coulomb threshold corrections and the finite top-quark width, perturbative QCD predicts a top-pair production cross section of about ten picobarns in the LHC threshold bin, with an excess of about four picobarns ove
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 07:31 UTC pith:Y3KKOGSM
load-bearing objection Solid threshold-resummation Monte Carlo work whose central ATLAS-bin number is generator-dependent and should not be quoted as a robust QCD prediction. the 3 major comments →
Top-Antitop Production and Decay at Threshold at the LHC in QCD Perturbation Theory
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 central claim is that the observed threshold excess is fully contained in the resummed perturbative series. Building on a spectral-density picture in which the top-antitop pair interacts through the Coulomb potential and decays through a finite width, the paper implements a resummation of all terms of order (alpha_s/v)^n. The all-order factor multiplies the Born cross section after subtracting the alpha_s term already present at NLO to avoid double counting. The resulting prediction for the threshold bin is a cross section near ten picobarns, with an excess over standard NLO generators near four picobarns. The paper derives an exact formula for the resummed cross section at fixed top and
What carries the argument
The central machinery is the non-relativistic spectral density of the top-antitop system, obtained from the resolvent of the Coulomb Hamiltonian evaluated at complex energy E + iGamma, where the top width Gamma smears the bound-state poles into Breit-Wigner shapes and the continuum into a Sommerfeld-enhanced distribution. This is packaged as a multiplicative threshold factor applied to the Born cross section, with the NLO alpha_s part subtracted to prevent double counting. The three generators differ in how off-shell top decays are treated: one uses the narrow-width approximation, one generates off-shell phase space with a projection onto an on-shell matrix element, and one uses a full reson
Load-bearing premise
Everything hinges on the approximate treatment of the top's finite width at fixed virtuality: the paper uses a recipe that is exact only after integrating over the top and antitop virtualities, and it states that it has not implemented the exact fixed-virtuality formula, so the predicted suppression of double-top events, worth about 1 pb, could change under a more exact treatment.
What would settle it
Implement the paper's exact fixed-virtuality resummation formula (eq. A.19) in a Monte Carlo and compute the double-top restricted cross section in the bin m_ttbar < 350 GeV, p* < 50 GeV; if the suppression below the nominal 2m_t threshold disappears, or the predicted enhancement shifts by more than about 1 pb, the claimed compatibility with the observed excess would be put in doubt.
If this is right
- In the threshold bin, the full resummed prediction is about 10 pb, about 4 pb above standard NLO, so the observed pseudoscalar excess is compatible with perturbative QCD without a genuine toponium bound state.
- Threshold corrections beyond order alpha_s^3 contribute less than about 1 pb in the colour-singlet channel, so the enhancement is dominated by the first three perturbative orders.
- Finite-width effects are small for coarse bins, but become visible at the 1-2 pb level through running-width corrections and through the kinematic suppression of double-top events below the nominal 2m_t threshold.
- The 'toponium' contribution is part of the perturbative cross section and is calculable; with realistic top-virtuality cuts it is reduced by about 1 pb relative to the narrow-width treatment.
- Using a fixed scale associated with the invariant-mass cut rather than an event-by-event scale reduces the predicted threshold effect by several picobarns for larger mass cuts, so experimental analyses should be aware of the choice of scale.
Where Pith is reading between the lines
- If the exact fixed-virtuality formula the paper derives but does not implement were implemented in a Monte Carlo, the predicted double-top suppression below threshold could change; the paper's own numbers suggest the shift could be of order 1 pb, enough to matter in a precision comparison.
- A cleaner experimental test would be to fit the unfolded m_ttbar lineshape in the threshold region to a single resummed prediction with smearing, rather than adding a separate pseudoscalar signal on top of a baseline: the paper's logic implies such two-sample fits risk double counting.
- Tightening the p* cut or using higher-luminosity measurements with better invariant-mass resolution would make finite-width effects larger, and the approximate recipe would likely need to give way to the exact formula, as the paper notes for future lepton-collider studies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs three NLO+PS generators (thr1, thr2, bb4l) for top-antitop production near the ttbar threshold, implementing all-order Coulomb/threshold-enhanced corrections in QCD and, in the thr2 and bb4l versions, including finite-width effects through an approximate recipe derived in Appendix A. The authors use these generators to compute the cross section in the ATLAS bin (m_ttbar < 350 GeV, p* < 50 GeV), separate double-top, single-top and no-top contributions, estimate bound-state (toponium-like) contributions, and compare with the NRQCD calculation of ref. [19]. Their central claim is that the ATLAS/CMS near-threshold pseudoscalar excess is compatible with perturbative QCD without invoking a genuine toponium bound state: they expect the ATLAS-bin cross section to be near 10 pb, with an excess of about 4 pb over standard NLO results. The paper also explicitly states that the exact fixed-virtuality resummation formula, eq. (A.19), has been derived but not implemented, with the approximate recipe (A.24) being used instead.
Significance. If the central claim holds, this is an important result: it would reinterpret a widely discussed 'toponium' excess as an ordinary perturbative-QCD threshold effect and provide public NLO+PS tools that can be used by experimental analyses. The paper has several genuine strengths: the derivation in Appendix A is careful and reduces to the Fadin-Khoze inclusive result in the appropriate limit; the numerical outputs are validated against the independent NRQCD calculation of ref. [19] and against a semi-analytic phase-space computation (fig. 7); the codes are made public; and the authors are unusually explicit about the limitations of their finite-width treatment. The manuscript is not fitted to the data, and the comparison points are external. The main weakness is that the headline numerical prediction is not stable across the paper's own preferred generators without a gauge-invariant event definition, and the approximate finite-width recipe leaves an unquantified systematic in the quoted 10/4 pb numbers.
major comments (3)
- [Section 5, Table 3] The headline prediction 'near ten picobarn, excess near four picobarn' is not robust across the paper's own NLO generators. In Table 3, the thr2 NLO result with running-width corrections gives dt+st+nt ≈ 9.80 + 0.715 ≈ 10.5 pb, while bb4l NLO gives dt+st+nt ≈ 13.2 + 11.9 ≈ 25.1 pb. This is not a threshold-resummation effect: the baseline values already differ by a factor of about 2.7 (≈7.1 pb vs ≈19.2 pb). The paper attributes the difference to 'single top events leaking into the double top sample', but Section 3.4 explicitly states that the single/double resonant separation is not gauge invariant and is 'to a large extent arbitrary'. Since the ATLAS bin does not impose the 15 GeV virtuality window used to define 'double top', the physical cross section in that bin includes the single-top-like/off-shell events. Discarding the st column changes the bb4l prediction by about 12 pb. The quot
- [Section 3.5 and Appendix A] The finite-width treatment is the load-bearing approximation of the paper. The exact fixed-virtuality formula (A.19) is derived but not implemented; instead the paper uses the recipe (A.24), which reproduces the inclusive Fadin-Khoze result only after integration over the top and antitop virtualities. The manuscript itself concedes (Section 4, near figs. 5–7, and Conclusions) that 'we cannot be absolutely certain that this effect survives when a more exact treatment will become available.' This caveat is not merely cosmetic: the predicted double-top suppression (which reduces the bound-state contribution by about 1 pb), the running-width effects, and the event-composition tables all depend on the virtuality distribution, which is exactly what (A.19) would modify. The authors should either provide a quantitative estimate of the difference between (A.24) and (A.19) in a simplified but cont
- [Section 7, Table 7] The choice between the default event-by-event threshold scale and the fixed-scale prescription changes the central ATLAS-bin prediction by amounts comparable to the quoted excess. For example, in Table 7 the thr2-wc1 NLO full cross section with running scale is 10.65 pb, while the fixed-scale result is 10.09 pb; the bound-state-only contribution changes from 2.62 pb to 2.24 pb. This is not an error, but it shows that the 'near ten pb' and 'near four pb' figures are scale-prescription dependent at the level of about 0.5–1 pb. The paper should make this explicit in the conclusions and, ideally, quote a range that folds in the running/fixed-scale difference together with Table 4 scale variations.
minor comments (4)
- [Section 7] In the text before Table 7, 'the last column of fig. 7' should read 'the last column of table 7'; the discussion is about the table, not the figure.
- [Conclusions] Typo: 'When running with effects are included' should be 'When running-width effects are included'.
- [Table captions] The captions of Tables 1–8 contain a typo: 'T able 1' should be 'Table 1', etc.
- [Appendix A, eq. (A.22)] The notation for the Green's function is inconsistent: eq. (A.19) uses R(x=0,k_f), while eq. (A.22) writes R(0x,k_f). Please define the position-space argument clearly and use a consistent notation throughout the appendix. There are also typographical issues with vector arrows in several equations.
Circularity Check
No significant circularity: the central ATLAS-bin cross-section and excess are computed from the generators and compared with external data/NRQCD; the only by-construction element is the finite-width recipe (A.24), which is explicitly calibrated to reproduce the inclusive Fadin–Khoze result.
specific steps
-
other
[Conclusions (Section 9) and Appendix A, eq. (A.24)]
"In essence, the recipe is such that the top and antitop can be produced with different virtualities and that after integration over their virtuality one reproduces the fully inclusive result of ref. [3]."
The implemented finite-width recipe (A.24) is constructed, by its normalization, to integrate to the Fadin–Khoze spectral density of refs. [3,4]; hence the inclusive m_ttbar distribution obtained from it is the input FK result by construction, not an independent prediction. This is explicitly acknowledged. The paper's central claims are more exclusive (ATLAS-bin m_ttbar<350, p*<50, dt/st/nt decomposition) and are not calibrated to that inclusive target; they are cross-checked against external NRQCD results [19,32], so the by-construction normalization does not force the 10 pb / 4 pb conclusions.
full rationale
The derivation chain is largely self-contained. The threshold-enhanced corrections are taken from the standard Coulomb resummation (NRR/Fadin–Khoze), extended here to higher (α_s/v)^n orders, and implemented in three NLO+PS generators. The central numerical claims—~10 pb in the ATLAS bin and ~4 pb excess over standard NLO—are computed from these generators with fixed scales and compared with, not fitted to, the ATLAS/CMS excess and with the external NRQCD calculations of refs. [19,32]. No parameter is tuned to the experimental excess. Self-citations to NRR are present, but they are not load-bearing: the paper summarizes and re-derives the needed NRR formulas, adds new higher-order and finite-width content, and the main conclusions are also supported by independent external comparisons. The important caveats—the exact fixed-virtuality formula (A.19) is not implemented; the dt/st separation is not gauge invariant and is 'to a large extent arbitrary'; the bb4l vs thr2 NLO difference in table 3 is unresolved and attributed to single-top leakage—are limitations/robustness concerns, not circularity. The only element that is 'by construction' is the finite-width recipe (A.24), which is explicitly defined to reproduce the fully inclusive Fadin–Khoze integral after integrating over virtualities; that limits the inclusive m_ttbar shape as an independent test, but does not by itself force the more exclusive ATLAS-bin predictions. Hence score 1.
Axiom & Free-Parameter Ledger
free parameters (4)
- E_cut (delta-regulator width) =
1 GeV (varied by factor 2)
- h (minimum virtuality offset) =
2 GeV
- threshold-scale factor =
1 (varied 0.5 and 2)
- double-top classification window =
+/-15 GeV
axioms (7)
- domain assumption Near threshold (m_ttbar < 350 GeV, p* < 50 GeV, v ~ 0.24), the dominant QCD corrections are Coulomb-like, and the cross section is described by a non-relativistic spectral density times a smooth hard factor.
- domain assumption The production vertex is effectively pointlike; the threshold-enhanced series is obtained by dressing the Born ttbar configuration with an arbitrary number of Coulomb gluon exchanges.
- standard math Scaleless integrals of even powers of v times the plus distributions vanish (eq. 2.11), and the plus distribution can be regularized by the polynomial continuation (2.13)–(2.16).
- domain assumption In the threshold limit of gluon fusion, the colour singlet and octet fractions are 2/7 and 5/7, and these fractions also hold for the off-shell Born kinematics in bb4l.
- domain assumption The recipe of eq. (A.24), which reproduces the fully inclusive Fadin-Khoze result after integrating over top virtualities but does not implement the exact fixed-virtuality formula (A.19), is adequate for LHC observables with invariant-mass resolution much larger than the top width.
- domain assumption In the real-emission soft limit, the eikonal factors can be evaluated with the on-shell-projected top momenta (P_t, P_tbar) rather than the off-shell ones.
- domain assumption The alpha_s entering the Coulomb potential should be evaluated at the scale sqrt(m sqrt(E^2 + Gamma_t^2)), converted to the MS scheme via eq. (4.1), with an event-by-event scale for the generator default.
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.
Reference graph
Works this paper leans on
-
[3]
V. S. Fadin and V. A. Khoze,Threshold Behavior of Heavy Top Production in e+ e- Collisions,JETP Lett.46(1987) 525–529
1987
-
[4]
V. S. Fadin, V. A. Khoze, and T. Sjostrand,On the Threshold Behavior of Heavy Top Production,Z. Phys. C48(1990) 613–622
1990
-
[5]
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]
Pith/arXiv arXiv 2025
-
[6]
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]
Pith/arXiv arXiv 2021
-
[7]
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)]
Pith/arXiv arXiv 2010
-
[8]
L. D. Landau,On the angular momentum of a system of two photons,Dokl. Akad. Nauk SSSR60(1948), no. 2 207–209
1948
-
[9]
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
1950
-
[10]
M. A. Braun,Positronium singularities in quantum electrodynamics and perturbation theory, Zh. Eksp. Teor. Fiz.54(1968) 1220–1227
1968
-
[11]
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]
Pith/arXiv arXiv 2014
-
[12]
M. Beneke and P. Ruiz-Femenia,Threshold singularities, dispersion relations and fixed-order perturbative calculations,JHEP08(2016) 145, [arXiv:1606.02434]
Pith/arXiv arXiv 2016
- [13]
-
[14]
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]
Pith/arXiv arXiv 2024
-
[15]
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]. – 37 –
Pith/arXiv arXiv 2007
-
[16]
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]
Pith/arXiv arXiv 2007
-
[17]
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]
Pith/arXiv arXiv 2007
-
[18]
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]
Pith/arXiv arXiv 2016
-
[19]
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]
Pith/arXiv arXiv 2025
-
[20]
Sommerfeld,Atombau und spektrallinien,Bd.2(Vieweg, Braunschweig, 1939)
A. Sommerfeld,Atombau und spektrallinien,Bd.2(Vieweg, Braunschweig, 1939)
1939
-
[21]
SakharovJETP18(1948) 631
A. SakharovJETP18(1948) 631
1948
-
[22]
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]
Pith/arXiv arXiv 2022
-
[23]
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]
Pith/arXiv arXiv 2021
-
[24]
F. Bach, B. C. Nejad, A. Hoang, W. Kilian, J. Reuter, M. Stahlhofen, T. Teubner, and C. Weiss,Fully-differential Top-Pair Production at a Lepton Collider: From Threshold to Continuum,JHEP03(2018) 184, [arXiv:1712.02220]. [25]NNPDFCollaboration, R. D. Ball et al.,Parton distributions for the LHC Run II,JHEP 04(2015) 040, [arXiv:1410.8849]
Pith/arXiv arXiv 2018
-
[26]
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]
Pith/arXiv arXiv 2015
-
[27]
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]
Pith/arXiv arXiv 2021
-
[28]
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]
Pith/arXiv arXiv 2021
-
[29]
Fischler,Quark - anti-Quark Potential in QCD,Nucl
W. Fischler,Quark - anti-Quark Potential in QCD,Nucl. Phys. B129(1977) 157–174
1977
-
[30]
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
1980
-
[31]
M. Cacciari, G. P. Salam, and G. Soyez,The anti-k t jet clustering algorithm,JHEP04 (2008) 063, [arXiv:0802.1189]
Pith/arXiv arXiv 2008
-
[32]
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
-
[33]
M. J. Strassler and M. E. Peskin,The Heavy top quark threshold: QCD and the Higgs,Phys. Rev. D43(1991) 1500–1514. – 38 –
1991
-
[34]
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]
Pith/arXiv arXiv 1998
-
[35]
A. H. Hoang and M. Stahlhofen,The Top-Antitop Threshold at the ILC: NNLL QCD Uncertainties,JHEP05(2014) 121, [arXiv:1309.6323]
Pith/arXiv arXiv 2014
-
[36]
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]
Pith/arXiv arXiv 2015
-
[37]
T. Sj¨ ostrand, V. A. Khoze, and C. T. Preuss,Top Pair Threshold Revisited, arXiv:2605.19546. – 39 –
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.