{"id":"b5aab9c0-a9f4-4925-86ae-dc871df3badd","arxiv_id":"2411.10856","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The author obtains first numerical values for the two-loop gg to t tbar g amplitude finite remainder at leading colour, taming elliptic integrals with a special function basis.","lead":"This paper reports first numerical results for the two-loop quantum chromodynamics (QCD) calculation that predicts top-quark pair production with an extra jet at the Large Hadron Collider. It is a step toward more precise theoretical predictions for backgrounds and for measuring the top quark mass.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite remainder values in Table 1 rely on an unverified order-epsilon^0 evaluation; the stated pole and gauge checks do not constrain it.","rationale":"The reader's weakest_assumption focuses on the mapping of hexagon-triangle master integrals to the pentagon-box families and on the correctness of the published master integrals. That is a legitimate concern, but I think it is not the single most load-bearing one: a mistake in that mapping, or in the master integrals, would most likely affect the pole structure as well, and the paper's pole checks would then have a reasonable chance of catching it. The deeper gap is that the finite remainder itself — the order-epsilon^0 coefficient, which is the only new numerical output — has no independent verification. The Ward identity checks gauge invariance of the whole amplitude, and the IR/UV pole checks constrain only the singular parts. Neither check is sensitive to an error that changes the epsilon^0 term. The paper explicitly states that elliptic contributions appear only at order epsilon^0, so the construction of the special function basis and its numerical evaluation by DiffExp are precisely the untested components that determine Table 1. The finite-field reconstruction of rational coefficients is also performed numerically and could contain reconstruction errors that again affect only the finite part. A second, independent evaluation of the same amplitude at one benchmark point would settle this. Since the paper is a proceedings contribution with clearly labelled preliminary results and appropriate caveats, the conditional verdict remains appropriate; no change to the reader's verdict is needed. My concern sharpens the reason for the condition rather than altering it.","tokens_in":7281,"tokens_out":2953,"duration_ms":32562,"concrete_test":"Reproduce one row of Table 1 independently. Using the rationalised phase-space values (currently 'available on request'), evaluate the two-loop leading-colour gg -> t tbar g amplitude by an independent route: generate diagrams, reduce with NeatIBP, and compute the pentagon-box master integrals with a different method, e.g. AMFlow or pySecDec sector decomposition with epsilon expansion. Recombine into A^(2)_LC(+++++; n_t n_tbar) and compare real and imaginary parts with Table 1 (first row: 19.03 - 3.108 i). Agreement within numerical integration errors (say <1%) would settle the finite-remainder evaluation; any significant deviation shows the concern lands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the numerical finite remainder A^(2)_LC (+++++) in Table 1. The two stated checks — Ward identity and analytic UV/IR pole matching — are necessary but not sufficient: both are insensitive to an error in the order-epsilon^0 coefficient. The elliptic special functions enter exactly at epsilon^0 (Sec. 2), and the values in Table 1 are produced by evaluating those functions with DiffExp. A bug in the newly constructed special-function basis, in the DiffExp evaluation, or in the finite-field reconstruction of the epsilon^0 rational coefficients would change every entry in Table 1 while leaving the stated checks unchanged. The asserted mapping of hexagon-triangle topologies to the pentagon-box families (Sec. 2, Figs. 1d–1f) is not demonstrated; if the mapping is wrong only in the finite part, the pole check would not reveal it. There is no independent numerical cross-check of the finite remainder (e.g., sector decomposition or a second differential-equation solver). This is a correctness risk on the paper's only new numerical result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports on progress toward the two-loop QCD amplitudes for pp -> tbar t j, focusing on the gg -> tbar t g channel at leading colour. After colour decomposition and a spinor-helicity setup, the amplitude is reduced to master integrals of three pentagon-box families, with hexagon-triangle topologies claimed to be mapped to them. The paper describes a special-function basis in which elliptic functions start only at order epsilon^0, allowing analytic identification of UV and IR poles, and presents a table of five numerical evaluations of the finite remainder A^(2)_LC(+++++; n_t n_tbar) obtained by evaluating the basis with DiffExp at benchmark phase-space points.","tokens_in":7428,"tokens_out":5760,"duration_ms":59064,"significance":"If correct, these are the first numerical two-loop finite remainders for the leading-colour gg -> tbar t g helicity amplitude, a key ingredient for NNLO predictions for ttbar+jet production. The paper deserves credit for the analytic structure: the pole structure is checked against established UV renormalization and IR subtraction, the Ward identity is used for gauge invariance, and the construction of a basis where elliptic sectors enter only at epsilon^0 is a useful methodological step. However, the numerical finite part, which is the only new quantitative result, lacks independent verification, and the central topology mapping is not documented; the result is therefore promising but not yet validated.","major_comments":[{"comment":"The central claim of the paper is the set of finite remainders A^(2)_LC(+++++; n_t n_tbar) in Table 1. These are values at order epsilon^0, but the checks reported in Section 3 (Ward identity and analytic matching of UV/IR poles) are insensitive to an error in the epsilon^0 coefficient; a bug in the special-function basis, in the DiffExp evaluation, or in the finite-field reconstruction would change every entry in Table 1 while leaving these checks unchanged. The manuscript should provide an independent numerical cross-check of the finite remainder at the benchmark points, for example by sector decomposition or by a second evaluation method, and should state the numerical precision of the quoted digits.","section":"Section 3 / Table 1"},{"comment":"The reduction of the hexagon-triangle master integrals to the pentagon-box families of Refs. [24,25] is asserted without proof. This mapping is load-bearing for the finite part because elliptic contributions are stated to enter only at order epsilon^0, so the analytic pole check cannot detect a mapping error at finite order. The mapping should be documented in an appendix or supplied as a transformation table, or validated independently.","section":"Section 2, Figs. 1d-1f"},{"comment":"The caption states that the rationalised kinematic values are 'available on request'. Since the numerical values in Table 1 are the only new quantitative result, the exact kinematics (rationalised invariants, signs, and values of tr5) and the corresponding numerical outputs should be provided in an ancillary file for reproducibility. 'Available on request' is not sufficient for a published result.","section":"Table 1 caption"}],"minor_comments":[{"comment":"The variable x_vec is introduced as 'a set of momentum twistor variables or Mandelstam invariants'; the manuscript should specify which choice is actually used in the numerical evaluations in Table 1.","section":"Section 2, Eq. (1)"},{"comment":"The computation of the mass counterterm is mentioned but not described; a brief statement of the renormalization scheme (on-shell vs MSbar) and how the counterterm is implemented in the finite remainder would be helpful.","section":"Section 3"},{"comment":"The phrase 'first numerical evaluations' should be accompanied by an estimate of the numerical accuracy of the DiffExp evaluation; Table 1 quotes four significant figures without stating how many are stable.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"This is a short proceedings contribution, but the claims go beyond a status report: Table 1 is a numerical result. I would ask the editors to require the missing validation or at least an explicit statement that the numbers are preliminary and should not be cited quantitatively. The reliance on the authors' own Refs. [24,25] is not circular, but the missing mapping details make the result difficult to audit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a proceedings paper, so the frame matters: its only new physics content is Table 1, five numerical finite-remainder values for the two-loop all-plus gg->tbar tg amplitude at leading colour. Compared to the cited literature, which stops at master integrals, that is genuinely new. The paper is clearly written and openly describes itself as a status report; the method—finite-field reconstruction, a pentagon-function-style special-function basis, DiffExp evaluation—is standard for this group and is applied sensibly. The mapping of the hexagon-triangle families to the previously computed pentagon-box master integrals is plausible and builds on published independent work. The Ward identity and UV/IR pole checks are nontrivial, and passing them is real evidence that the pole structure is right.\n\nThe soft spot is exactly where the stress-test points. The finite part is not constrained by either check. An error in the construction of the elliptic special-function basis, in the DiffExp evaluation at order epsilon^0, or in the finite-field reconstruction of the rational coefficients would change every entry in Table 1 and leave the stated checks unchanged. The mapping of the hexagon-triangle topologies is asserted in a single sentence; if it were wrong only in the finite part, the pole check would not notice. There is no independent evaluation—no second differential-equation solver, no sector decomposition, no public code or data. \"Rationalised values available on request\" is not a reproducibility artifact. None of this makes the numbers wrong; it does mean Table 1 is currently an unverified claim about the hardest part of the computation. The paper also doesn't state the numerical precision of the DiffExp evaluations, which for a table of numbers should be included.\n\nWho benefits? People computing tt+j at two loops, and the broader multi-loop community watching for a first benchmark. As a proceedings contribution it does its job: it records a milestone and describes the method. For a full journal paper I would want the special-function basis and the topology mapping spelled out, plus at least one independent cross-check or a public code. For a proceedings, this is enough to warrant peer review, with the expectation that the numbers are preliminary. I'd take the result as promising but not yet something to build on.","headline":"First finite-remainder numbers for gg->tbar tg at two loops, but the finite part is exactly what is not yet independently constrained.","tokens_in":7963,"tokens_out":3250,"would_cite":true,"duration_ms":36103,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Bx"],"model":"deepseek-v4-flash","headline":"The paper reports first numerical evaluations of the two-loop leading-colour helicity amplitude for $gg\\to t\\bar{t}g$, with the infrared and ultraviolet poles identified analytically, supplying the hardest missing double-virtual…","keywords":["two-loop amplitudes","QCD corrections","top-antitop-jet production","helicity amplitudes","leading colour","elliptic Feynman integrals","finite remainders"],"falsifier":"Recompute the five phase-space points of Table 1 with a different integral-reduction strategy and a different numerical evaluator; any disagreement larger than the quoted precision would show that the diagram mapping or the special-function evaluation is wrong.","tokens_in":7018,"feed_emoji":"⚛️","tokens_out":12196,"duration_ms":121730,"temperature":0.7,"pith_summary":"Top-antitop-jet production is the largest associated-production process at hadron colliders, but its complete next-to-next-to-leading-order (NNLO) QCD prediction is still out of reach because the two-loop amplitudes have not been computed. This paper focuses on the hardest channel, $gg\\to t\\bar{t}g$, and presents the framework and first numerical results: a special-function basis in which the elliptic integrals enter only at order $\\varepsilon^0$, so the ultraviolet and infrared poles can be identified analytically and the finite remainder extracted. The concrete result is Table 1, which lists five benchmark numerical values of the leading-colour finite remainder for the all-plus helicity configuration. If these values are correct, they remove the gluon-channel double-virtual blockage for NNLO $pp\\to t\\bar{t}j$ predictions.","feed_headline":"Gluon channel of top-antitop-jet gains first two-loop values","feed_subtitle":"The gluon channel is the hardest missing ingredient for NNLO top-antitop-jet predictions.","key_machinery":"The load-bearing machinery is the special-function basis built on the pentagon-box master integrals, where a pentagon-box topology is a two-loop diagram formed by a pentagon and a box sharing a loop. The six leading-colour integral families (three hexagon-triangle and three pentagon-box topologies) are reduced to the pentagon-box set, and the amplitude is written as $\\sum_i \\sum_{k=-4}^{0} \\varepsilon^k r_{ki}(\\vec x) F_i(\\vec x)$. The basis is chosen so that most $F_i$ are multiple polylogarithms and the few elliptic functions contribute only from $\\varepsilon^0$ onward; this separation is what makes the pole structure analytically accessible. The finite remainders are then obtained by evaluating the special functions through one-dimensional series expansions, using differential equations.","core_discovery":"The author establishes that the two-loop leading-colour helicity amplitude for $gg\\to t\\bar{t}g$ can be brought to a numerically evaluable form: all master integrals of the hexagon-triangle topologies are mapped to the previously computed pentagon-box master integrals, and the amplitude is expanded in a basis of special functions where the elliptic parts only contribute starting at order $\\varepsilon^0$. With this structure the ultraviolet and infrared poles are identified analytically and checked against UV renormalization and IR subtraction, and the finite remainder in the 't Hooft–Veltman scheme is extracted. The central concrete result is Table 1, which lists five phase-space values of $A^{(2)}_{\\mathrm{LC}}(+++++; n_t n_{\\bar t})$ in the physical region; these are reported as the first numerical evaluations of this finite remainder.","pith_inferences":["Beyond the paper: the same elliptic-sector separation could be applied to other 2→3 processes with internal massive propagators, such as the quark-initiated channels of $pp\\to t\\bar{t}j$, where the pole structure is likely to be handled in the same way.","Beyond the paper: if the hexagon-to-pentagon mapping is exact and the special-function basis is reconstructed analytically rather than numerically, the amplitude could be evaluated much faster over full phase space, opening the way to differential distributions at NNLO.","Beyond the paper: the benchmark values could be used to quantify the size of subleading-colour corrections once a full-colour computation becomes available, telling whether the leading-colour approximation is sufficient for precision phenomenology."],"forward_implications":["The gluon channel is the most difficult production channel; with its finite remainder in hand, the double-virtual contribution at leading colour can be combined with one-loop and real-radiation contributions to move toward NNLO predictions for $pp\\to t\\bar{t}j$.","The five benchmark phase-space points in Table 1 give independent groups a concrete numerical target for checking their own calculations of this two-loop amplitude.","The demonstrated separation of elliptic functions away from the pole order shows that analytic UV and IR pole identification remains feasible even when elliptic Feynman integrals are present.","The framework produces the finite remainder in the 't Hooft–Veltman scheme, the scheme-dependent object needed for a future NNLO phenomenology analysis."],"supporting_citations":[{"why":"computes the pentagon-box master integrals for a planar topology contributing to pp→tt j, the reduction target for the hexagon-triangle integrals.","marker":"[24]"},{"why":"provides the leading-colour two-loop integrals for tt+jet production, on which the amplitude reduction and the finite remainder depend.","marker":"[25]"},{"why":"supplies the one-loop helicity amplitudes and the spinor and structure definitions used to build the two-loop amplitude.","marker":"[23]"},{"why":"gives the two-loop leading-colour helicity amplitude framework for top-pair production in gluon fusion, including the massive-spinor formalism.","marker":"[39]"},{"why":"supplies the infrared-singularity subtraction for amplitudes with massive partons, used to check the analytic pole structure.","marker":"[40,41]"},{"why":"develops the generalised series-expansion method used to evaluate the special functions at the benchmark points.","marker":"[42]"},{"why":"provides the differential-equations-based numerical evaluation tool used to obtain the Table 1 values.","marker":"[43]"}],"fun_headline_variants":["First two-loop numbers for gluon-induced top-antitop-jet","Two-loop top-antitop-jet: gluon channel evaluated","Elliptic integrals tamed for NNLO top-antitop-jet gluon channel","Gluon channel of top-antitop-jet gets first two-loop values","Two-loop gg→ttbar g amplitude: first finite remainders"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes without proof that the more complicated loop diagrams can be rewritten in terms of the previously published pentagon-box integrals, and that those published integrals are correct.","fun_headline_variants_meta":{"raw":{"variants":["First two-loop numbers for gluon-induced top-antitop-jet","Two-loop top-antitop-jet: gluon channel evaluated","Elliptic integrals tamed for NNLO top-antitop-jet gluon channel","Gluon channel of top-antitop-jet gets first two-loop values","Two-loop gg→ttbar g amplitude: first finite remainders"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000699,"raw_usage":{"total_tokens":3100,"prompt_tokens":828,"completion_tokens":2272,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":444,"completion_tokens_details":{"reasoning_tokens":2173}},"tokens_in":444,"tokens_out":2272,"duration_ms":17434,"temperature":1.0,"reasoning_tokens":2173,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:13:09.481245+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the five phase-space points of Table 1 with a different integral-reduction strategy and a different numerical evaluator; any disagreement larger than the quoted precision would show that the diagram mapping or the special-function evaluation is wrong.","supporting_citations":[{"cited_title":"Badger, M","cited_arxiv_id":null,"evidence_quote":"computes the pentagon-box master integrals for a planar topology contributing to pp→tt j, the reduction target for the hexagon-triangle integrals."},{"cited_title":"Badger, M","cited_arxiv_id":null,"evidence_quote":"supplies the one-loop helicity amplitudes and the spinor and structure definitions used to build the two-loop amplitude."},{"cited_title":"Badger, E","cited_arxiv_id":null,"evidence_quote":"gives the two-loop leading-colour helicity amplitude framework for top-pair production in gluon fusion, including the massive-spinor formalism."},{"cited_title":"Moriello,Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops, JHEP 01, 150 (2020), doi:10.1007 /JHEP01(2020)150","cited_arxiv_id":null,"evidence_quote":"develops the generalised series-expansion method used to evaluate the special functions at the benchmark points."}],"review_version":1}