Pith. sign in

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 →

arxiv 2607.09539 v2 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 correctionsfinite-width effectstoponiumLHC physicsspin correlationsNLO event generators
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.

This paper tries to establish that the enhancement in top-antitop production near threshold seen at the LHC is a perturbative QCD effect, not a sign of a genuine toponium bound state. The authors construct three next-to-leading-order event generators that resum threshold-enhanced Coulomb corrections to all orders and account for the finite width of the top quark. In the threshold bin defined by m_ttbar < 350 GeV and p* < 50 GeV, they obtain a cross section near ten picobarns, about four picobarns above the cross section from standard NLO generators. They also find that corrections beyond the first three perturbative orders contribute less than one picobarn, so the effect is under perturbative control. A reader should care because this reframes the reported 'toponium excess' as an ordinary QCD cross-section enhancement that was missing from the baseline simulations.

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.

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

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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.

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

Referee Report

3 major / 4 minor

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)
  1. [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
  2. [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
  3. [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)
  1. [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.
  2. [Conclusions] Typo: 'When running with effects are included' should be 'When running-width effects are included'.
  3. [Table captions] The captions of Tables 1–8 contain a typo: 'T able 1' should be 'Table 1', etc.
  4. [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

1 steps flagged

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
  1. 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

4 free parameters · 7 axioms · 0 invented entities

The central calculation rests on the standard non-relativistic/Coulomb Green's-function treatment of the ttbar threshold (with finite-width insertions), plus several technical and heuristic choices (regularization parameters, colour decomposition, scale prescription, and an approximate finite-width recipe). No new entities are postulated.

free parameters (4)
  • E_cut (delta-regulator width) = 1 GeV (varied by factor 2)
    Regularizes the bound-state delta function in eq. (2.7); the paper states output is insensitive to this choice.
  • h (minimum virtuality offset) = 2 GeV
    Technical lower-bound parameter for sampling top/W virtualities in thr2 (eqs. 3.2–3.3).
  • threshold-scale factor = 1 (varied 0.5 and 2)
    Multiplicative factor on the Coulomb-coupling scale defined around eq. (4.1); scale variations in table 4 quantify the resulting uncertainty.
  • double-top classification window = +/-15 GeV
    Defines the dt/st/nt event categories in §4; affects the tabulated cross-section decomposition but not the inclusive results.
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.
    Used throughout §2, from eq. (2.1) onward.
  • 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.
    Assumed in Appendix A (fig. 9 and surrounding text) and in the NRR-based formula (2.1).
  • 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).
    Needed to define the (alpha_s/v)^4 terms in §2.2.
  • 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.
    Section 3.4: 'assuming a 2/7 colour singlet and a 5/7 colour octet fraction'.
  • 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.
    Section 3.5 and Conclusions; authors state they cannot be certain the double-top phase-space suppression survives the exact treatment.
  • 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.
    Section 3.3: modification of sigcollsoft2.f to restore soft cancellations.
  • 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.
    Paragraph after eq. (4.1): 'we proceed in a heuristic way'.

pith-pipeline@v1.3.0-alltime-deepseek · 30071 in / 16942 out tokens · 181436 ms · 2026-08-02T07:31:59.972689+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

34 extracted references · 23 linked inside Pith

  1. [3]

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

  2. [4]

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

  3. [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]

  4. [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]

  5. [7]

    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. [8]

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

  7. [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

  8. [10]

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

  9. [11]

    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. [12]

    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. [13]

    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. [14]

    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. [15]

    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]. – 37 –

  14. [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]

  15. [17]

    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]

  16. [18]

    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]

  17. [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]

  18. [20]

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

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

  19. [21]

    SakharovJETP18(1948) 631

    A. SakharovJETP18(1948) 631

  20. [22]

    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. [23]

    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]

  22. [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]

  23. [26]

    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]

  24. [27]

    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]

  25. [28]

    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]

  26. [29]

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

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

  27. [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

  28. [31]

    Cacciari, G

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

  29. [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

  30. [33]

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

  31. [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]

  32. [35]

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

  33. [36]

    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]

  34. [37]

    Sj¨ ostrand, V

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