{"id":"fd7a8810-0efd-4d4e-ac3f-af5ff2d066fd","arxiv_id":"2505.10381","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"First invariant-mass threshold resummation for ttbar ttbar production, giving NLO+NLL' predictions for the invariant-mass distribution and total cross section.","lead":"This paper calculates improved predictions for four top quark production at the LHC by adding soft-gluon resummation at NLL' accuracy to the next-to-leading-order result. The result is a more precise Standard Model benchmark for a rare process used to probe new physics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Total cross section changes by ~40% between the AMT and IMT resummation schemes (Table 2); this scheme spread is not in the quoted uncertainties, so the 'most precise' claim is premature.","rationale":"The paper is a careful NLL' resummation calculation. It contains a useful internal validation: the expansion of the NLL' result to NLO reproduces NLO excluding qg channels within a few percent (Appendix C), and the delta-dependence of the Coulomb-limit prescription (Eq. 3.6) is numerically negligible. These checks give real support to the one-loop ingredients and to the Coulomb regularization. My concern is not that the calculation is internally wrong, but that the central quantitative claim, 'most precise QCD predictions', is not supported by the uncertainty estimate. Table 2 shows that, at the same nominal NLO+NLL' accuracy and with the same mu0=M/2, the IMT-res central value is 12.38 fb while the AMT-res value from Ref. [58] is 17.36 fb. The difference is traced in Appendix A to the scale-ratio logarithms in the jet function: replacing Q^2 by M^2 in g2 (the 'IMT-res test') changes the K-factor from 0.90 to 1.30. So a roughly 40% shift in the central value is caused by a scheme choice. The 7-point scale bands and the convergence across mu0 choices within the IMT scheme do not cover this spread. Therefore, unless the scheme dependence is either shown to be spurious or included as a systematic uncertainty, the paper should soften the 'most precise' claim. My recommended verdict is CONDITIONAL: accept the calculation, but require the authors to quantify this scheme uncertainty or revise the claim. This differs from the reader's assessment, which identified the Coulomb-limit treatment as the weakest point; the delta-sensitivity test and the NLL'|NLO reproduction of NLO make that issue less damaging than the scheme ambiguity visible in the paper's own comparison.","tokens_in":19120,"tokens_out":19149,"duration_ms":217307,"concrete_test":"Recompute the NLO+NLL' total cross section with the g2 scale-ratio logarithm evaluated at the interpolated scale Q_lambda = Q (M/Q)^lambda, for lambda = 0, 0.25, 0.5, 0.75, 1, with the corresponding NLO-expansion term adjusted identically, using the same mu0=M/2 setup as Table 2. If the cross section varies by more than the 7-point scale band over lambda (the endpoints are already known: about 12.4 fb and about 17.9 fb), the AMT-IMT scheme spread must be included in the theory uncertainty and the 'most precise' claim changed; if the variation is confined to a narrow range, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's own comparison in Table 2 shows that the total cross section at the same nominal accuracy (NLO+NLL') is not stable under the choice of threshold variable. With mu0=M/2, the invariant-mass threshold result of this work is 12.38 fb (K=0.90), while the absolute-mass threshold result of Ref. [58], computed on the same footing, is 17.36 fb (K=1.26). Appendix A traces the difference to scale-ratio logarithms in the jet function g2, Eq. (A.1), and shows that replacing Q^2 by M^2 in g2 (the 'IMT-res test') changes the K-factor from 0.90 to 1.30. Thus a roughly 40% shift in the central value is generated by a scheme choice that is not included in any quoted scale uncertainty. The Section 5 claim that these results constitute 'the most precise QCD predictions' is therefore not supported by the uncertainty estimate: the 7-point bands and the convergence across mu0 choices within the IMT scheme do not cover the AMT-IMT discrepancy.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first invariant-mass threshold (IMT) resummation for the production of four top quarks at the LHC. It computes the invariant-mass distribution and the total cross section for $pp\\to t\\bar t t\\bar t$ at 13.6 and 13 TeV at NLL' accuracy, matched to NLO QCD and to NLO QCD+EW, with the resummation performed in Mellin space using a colour-space soft function. The main numerical results are NLO+NLL' total cross sections of 12.38, 12.00 and 12.25 fb for $\\mu_0=M/2$, $Q/2$ and $H_T/2$ at 13.6 TeV, with a substantial reduction of the 7-point scale uncertainty relative to NLO. The paper also validates the logarithmic expansion against exact NLO excluding $qg$ channels, tests the dependence on the Mellin contour and on the Coulomb cutoff $\\delta$, and compares the IMT scheme with the earlier absolute-mass threshold (AMT) scheme in Appendix A.","tokens_in":19401,"tokens_out":7537,"duration_ms":75415,"significance":"If the results are correct, this is a technically demanding and useful computation: it extends invariant-mass threshold resummation to a process with six coloured particles, provides differential predictions beyond NLO for a rare process, and includes several robustness checks. The strongest points are the explicit validation in Appendix C, where the NLL' expansion reproduces the NLO no-$qg$ distribution within a few percent, the demonstrated insensitivity to technical parameters (Mellin contour and Coulomb cutoff), and the use of independent NLO tools (MG5 aMC@NLO and OpenLoops) for matching. However, the headline claim that these are the most precise QCD predictions is weakened by the large scheme spread documented in Table 2, which is not reflected in the quoted scale uncertainties.","major_comments":[{"comment":"The summary claim that these results constitute \"the most precise/accurate QCD predictions\" for four-top production is not supported by the quoted uncertainties. At the same nominal accuracy (NLO+NLL') and same central scale $\\mu_0=M/2$, the IMT-res cross section obtained here is 12.38 fb, while the AMT-res cross section computed with the same framework is 17.36 fb; the fictitious \"IMT-res test\" of Appendix A, which replaces $Q^2$ by $M^2$ in $g_2$ (Eq. (A.1)), moves the IMT-res result to 17.91 fb. This roughly 40% spread in the central value is larger than the 7-point scale bands (about $\\pm 20\\%$-$24\\%$) and is traced in Appendix A to threshold-variable logarithms. The authors should either incorporate this scheme ambiguity into the quoted theoretical uncertainty or explicitly restrict the \"most precise\" claim to the invariant-mass threshold scheme.","section":"Section 5 and Appendix A, Table 2"},{"comment":"The Coulomb-limit replacement $\\beta_{IJ}\\coth\\beta_{IJ}\\to 1$ for $b_{IJ}<\\delta$ is introduced through an effective-particle/Casimir-scaling picture, but the manuscript does not derive this replacement from a controlled resummation of Coulomb exchanges for a subset of particles at rest in a $2\\to 4$ process. The numerical $\\delta$-sensitivity tests and the NLO-reproduction test in Appendix C are encouraging, but varying $\\delta$ probes only the region in which the singular term is removed, not the validity of the finite part of the replacement. The authors should either provide a derivation of the effective-particle reduction, including the relation between the two forms of the finite terms in Eq. (3.5), or explicitly label this treatment as an approximation and estimate its effect on the final uncertainties.","section":"Section 3, Eq. (3.6)"}],"minor_comments":[{"comment":"The term $C^{(1)}$ in the definition of the one-loop hard function is not defined in the text; please specify its construction and its relation to the soft-collinear endpoint contributions when it is not captured by the NLL jet functions.","section":"Section 2, Eq. (2.5)"},{"comment":"The notation $\\gamma_{\\rm cusp}^{(0)}\\times 1$ is confusing; please state explicitly that the factor 1 is the finite part of $\\beta\\coth\\beta$ in the limit $b_{IJ}\\to 0$ and that the $(-i\\pi)/b_{IJ}$ term is dropped because Coulomb contributions are already contained in the hard function.","section":"Section 3, Eq. (3.6)"},{"comment":"The construction of the \"AMT-res\" differential distribution from the IMT-res distribution is described only in words; a formula analogous to Eq. (A.2) would make the reweighting procedure unambiguous.","section":"Appendix A"},{"comment":"The statement that the NLO spread decreases from 36% to 3% refers to the shown invariant-mass range; please state that range explicitly in the caption and in the text.","section":"Figure 3 and Section 4.2"}],"recommendation":"major_revision","confidential_remarks":"The technical computation appears sound and the validation in Appendix C is convincing. The main deficiency is the framing of the \"most precise\" claim in light of the AMT-IMT scheme spread; this is an interpretation issue that can be fixed in revision. The reliance on self-citations is not problematic given the continuity of the framework."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line first: this is a solid, transparent resummation calculation, and the NLO+NLL' invariant-mass distributions are genuinely new. But the paper's own Table 2 undercuts the \"most precise QCD predictions\" claim: at mu0=M/2 the total cross section is 12.38 fb in the invariant-mass scheme versus 17.36 fb in the absolute-mass scheme from Ref. [58], both labelled NLO+NLL'. That forty-percent spread is not covered by any quoted scale-uncertainty band. The authors diagnose it well in Appendix A, tracing most of it to the log(Q^2/mu_R^2) versus log(M^2/mu_R^2) term in the jet function g2, but the abstract and Section 5 do not carry that caveat. I would want that qualification in place before calling the prediction the most accurate available.\n\nWhat is genuinely good: the first application of invariant-mass threshold resummation to a six-colored-particle final state; a new differential prediction for the four-top invariant-mass spectrum; and a careful treatment of the Coulomb limit for a subset of particles at rest. The validation is also solid. Appendix C shows the NLL' expansion reproduces NLO without qg channels within 1-4% bin-by-bin and within 6% integrated, and the tests of Mellin contour and Coulomb-cutoff sensitivity are done properly. The reliance on earlier self-authored work for the formalism is not a red flag; the framework is established and the numerical result is new.\n\nSoft spots, in order. First and main: the threshold-scheme ambiguity is a load-bearing caveat, not a nitpick. If two legitimate resummation prescriptions at the same formal accuracy differ by 40% in the integrated rate, then \"reduced theoretical uncertainties\" is only true within the chosen scheme. The paper is honest in Appendix A but overclaims in the conclusions. Second: the one-loop soft boundary function S^(1) is stated as needed but not derived; the reader has to trust the cited literature. That is probably acceptable for a JHEP paper, but a few equations would have made it self-contained. Third: no public code or ancillary numerical tables are provided, so independent reproduction is not possible from the manuscript alone; that is a minor operational gap, not a scientific flaw.\n\nWho this is for: four-top phenomenologists and anyone using threshold resummation in multi-particle final states. It deserves a serious referee. My recommendation: send it to review, and ask the authors to acknowledge the AMT/IMT scheme ambiguity in the abstract and introduction, and ideally to include it as an additional uncertainty. With that revision, the paper is a useful and citable contribution.","headline":"A careful, transparent first invariant-mass threshold resummation for four-top production, whose own Table 2 shows a 40% scheme ambiguity that the 'most precise' claim does not cover.","tokens_in":19903,"tokens_out":2961,"would_cite":true,"duration_ms":31981,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By resumming soft-gluon logarithms in the invariant-mass threshold, this paper claims to obtain the most precise QCD predictions for four-top-quark production, with the spread across scale choices cut from 36% to 3%.","keywords":["four-top-quark production","threshold resummation","soft-gluon resummation","next-to-leading logarithmic accuracy","invariant-mass distribution","QCD precision predictions","top quark"],"falsifier":"Recompute the one-loop soft anomalous dimension for $2\\to 4$ massive final states in a colour basis that handles the $\\beta_{IJ}\\to 0$ singularities analytically, without the effective-particle cutoff, and re-evaluate the NLO+NLL$'$ cross section and invariant-mass distribution; a difference larger than the quoted scale bands (about 15-25%) would show the Coulomb-limit treatment is numerically significant rather than benign.","tokens_in":18929,"feed_emoji":"⚛️","tokens_out":11586,"duration_ms":107014,"temperature":0.7,"pith_summary":"Soft-gluon emissions near the invariant-mass threshold of the $t\\bar{t}t\\bar{t}$ system produce large logarithmic corrections, and this paper tries to establish that they can be resummed to all orders at NLL$'$ accuracy, giving the most precise QCD predictions for four-top-quark production at the LHC. The central quantitative payoff is a large reduction of scale uncertainty: the spread between predictions obtained with three different central scale choices drops from 36% at NLO to about 3% after matching the resummed result to NLO. If the claim holds, future LHC measurements of this rare process can be compared with a theoretical baseline whose main parametric ambiguity has been substantially reduced, which matters because four-top production is a direct probe of the top-Higgs coupling and of new physics.","feed_headline":"Resummation shrinks four-top-quark scale spread from 36% to 3%","feed_subtitle":"NLL' soft-gluon corrections shrink the scale spread and give the most precise four-top QCD cross section to date.","key_machinery":"The machinery is invariant-mass threshold resummation: in the limit where the invariant mass $Q$ of the $t\\bar{t}t\\bar{t}$ system approaches the partonic centre-of-mass energy, soft-gluon emissions produce logarithms that are resummed to all orders in Mellin space. The resummed partonic cross section factorises into a hard function $H$, a soft function $S$ obtained from the renormalisation-group evolution of the soft anomalous dimension matrix $\\Gamma$, and collinear jet functions $\\Delta_i$ for the incoming partons (Eq. (2.3)). At NLL$'$ accuracy the one-loop hard function and one-loop soft boundary condition are needed; the one-loop soft anomalous dimension for $2\\to n$ massive final states is computed from the cusp anomalous dimension $\\gamma^{(0)}_{\\rm cusp}(\\beta_{IJ}) = \\gamma^{(0)}_{\\rm cusp}\\,\\beta_{IJ}\\coth\\beta_{IJ}$, with the Coulomb limit $\\beta_{IJ}\\to 0$ handled by an effective-particle replacement below a cutoff $\\delta$ (Eq. (3.6)).","core_discovery":"The central claim is that invariant-mass threshold resummation—resumming logarithms of $1-\\hat\\rho_Q$ where $\\hat\\rho_Q = Q^2/\\hat s$ is the partonic threshold variable and $Q$ the invariant mass of the four-top final state—can be applied to a $2\\to 4$ process with six coloured particles, and that at NLL$'$ accuracy it removes most of the scale-choice dependence of the fixed-order prediction. Working in Mellin space, the resummed cross section is matched additively to exact NLO QCD and electroweak results. For $\\sqrt{S}=13.6$ TeV the matched NLO+NLL$'$ total cross sections are 12.38 fb, 12.00 fb and 12.25 fb for $\\mu_0=M/2$, $\\mu_0=Q/2$ and $\\mu_0=H_T/2$, respectively, with the spread among the three central scales falling from 36% at NLO to 3% at NLO+NLL$'$; the predicted NLL$'$ corrections themselves range from $-10\\%$ to $+18\\%$ depending on the scale choice, and the inclusion of electroweak corrections raises the cross section by 5-8%.","pith_inferences":["A further step the paper leaves implicit: the same scale-convergence improvement is likely to show up in other $2\\to n$ heavy processes (e.g. $t\\bar{t}b\\bar{b}$ or associated four-top-plus-jet production), which would make this formalism a general tool for sharpening BSM search backgrounds.","The sign flip of the NLL$'$ correction between $\\mu_0=M/2$ and $\\mu_0=Q/2$ is a residual signature that scale-log ambiguities are not fully removed; comparing the predicted shape of the invariant-mass distribution against Run 3 data could choose the preferred dynamical scale empirically.","The drop in scale spread from 36% to 3% should not be read as a 3% total theoretical error: the per-scale uncertainties quoted remain 15-24%, so the real achievement is removing the choice of central scale as a dominant source of disagreement, not eliminating higher-order uncertainty."],"forward_implications":["These NLO+NLL$'$ cross sections (about 12 fb at 13.6 TeV) provide the sharpest QCD baseline for comparing with the four-top measurements at the LHC.","The NLL$'$ corrections alter the shape of the invariant-mass distribution by up to about 25% (with sign and size depending on the scale choice), so differential data will be able to test the resummation rather than only the total rate.","The expansion of the NLL$'$ result to NLO reproduces the exact NLO QCD invariant-mass distribution within 1-4% (outside the small quark-gluon channel), so the resummation calculation doubles as a fast approximate-NLO description of this process.","Because the results are essentially independent of the Mellin contour parameters and of the Coulomb cutoff $\\delta\\in\\{0.1,0.01,0.001\\}$, the numerical implementation is stable under the choices the paper varied.","The same Mellin-space invariant-mass threshold machinery can now be applied to other multi-heavy-particle processes with more than two massive coloured final-state particles."],"supporting_citations":[{"why":"This reference provides the previous absolute-mass threshold resummation for four-top production, defines the colour-space dimension, and serves as the comparison baseline in Appendix A.","marker":"[58]"},{"why":"This reference supplies the general one-loop soft anomalous dimension for $2\\to n$ processes with massive partons, which is the starting point of the Coulomb-limit treatment.","marker":"[75]"},{"why":"This reference establishes threshold resummation for heavy-particle pair production including Coulomb corrections, motivating the effective-particle treatment of the Coulomb limit.","marker":"[71]"},{"why":"This reference provides the NLL$'$-type resummation formalism for top-quark pair production, including the one-loop soft boundary condition and the $\\bar N$ Mellin-space conventions used here.","marker":"[42]"},{"why":"This reference gives the exact NLO QCD and electroweak results for $t\\bar{t}t\\bar{t}$ production to which the resummed predictions are matched.","marker":"[31]"},{"why":"This reference defines the Minimal Prescription used for the inverse Mellin transform that returns the resummed results to momentum space.","marker":"[64]"}],"fun_headline_variants":["Four-top scale spread cut from 36% to 3% via resummation","NLL' resummation tames four-top scale uncertainty to 3%","Resummation cuts four-top scale uncertainty from 36% to 3%","Four-top production: resummation reduces scale spread from 36% to 3%","Most precise four-top cross section: NLL' resummation shrinks scale spread to 3%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that when a pair of top quarks in the final state is produced almost at rest, soft gluons effectively see the pair as a single particle with the combined colour charge, so the Coulomb-singular part of the one-loop soft anomalous dimension can be replaced by a Casimir-scaled effective-particle term below a small cutoff; if that replacement is wrong, the NLL$'$ predictions would shift outside the quoted uncertainties.","fun_headline_variants_meta":{"raw":{"variants":["Four-top scale spread cut from 36% to 3% via resummation","NLL' resummation tames four-top scale uncertainty to 3%","Resummation cuts four-top scale uncertainty from 36% to 3%","Four-top production: resummation reduces scale spread from 36% to 3%","Most precise four-top cross section: NLL' resummation shrinks scale spread to 3%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000594,"raw_usage":{"total_tokens":2782,"prompt_tokens":944,"completion_tokens":1838,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":1724}},"tokens_in":560,"tokens_out":1838,"duration_ms":12293,"temperature":1.0,"reasoning_tokens":1724,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:10:07.198645+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the one-loop soft anomalous dimension for $2\\to 4$ massive final states in a colour basis that handles the $\\beta_{IJ}\\to 0$ singularities analytically, without the effective-particle cutoff, and re-evaluate the NLO+NLL$'$ cross section and invariant-mass distribution; a difference larger than the quoted scale bands (about 15-25%) would show the Coulomb-limit treatment is numerically significant rather than benign.","supporting_citations":[],"review_version":1}