{"id":"55f3a358-72bd-43c1-9ca6-8c63f3668d90","arxiv_id":"2504.19909","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"First analytic reconstruction of the one-loop amplitude for 0 to anti-q, q, top, anti-top, Higgs using massive spinor-helicity, with expressions about twice as fast as the automatic code OpenLoops.","lead":"Physicists at Fermilab, Edinburgh, and Durham present compact analytic formulas for a one-loop quantum process that contributes to Higgs production with a top quark pair at the LHC. The work extends a powerful technique for extracting simple expressions from numerical calculations to particles with mass, which could speed up future high-precision collider predictions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Incomplete Gröbner basis leaves ansatz completeness empirical: numerical rank reduction could drop a generically independent monomial, so the 'exact' reconstruction is not fully established.","rationale":"The reader's weakest assumption identifies the same point: the reconstruction relies on an ansatz built without a complete Gröbner basis, with leftover redundancies removed numerically. I agree that this is the most load-bearing assumption, because the paper's headline result is an exact analytic form. If the numerical reduction accidentally removed a monomial that is independent on a part of the phase space not covered by the samples, then the published expressions would be wrong in that region, contradicting the central claim. The concern is not that the authors have been careless: they perform extensive cross-checks, including OpenLoops agreement, infrared pole checks, renormalization counterterms, and automated tests on physical and finite-field points. These checks make an undetected error unlikely but do not eliminate it, because none of them systematically tests the completeness of the ansatz on the full variety. The paper's own admission that no complete Gröbner basis was obtained makes this the one place where an algebraic gap could hide. The proposed test would close that gap by comparing the minimal ansatz reconstruction against a deliberately over-complete ansatz on a large, independent sample. If the two agree, the concern is resolved and acceptance is justified; if they disagree, the expressions need correction. Since the test is feasible with the tools already provided in the paper (Lips, Antares, p-adic phase-space generation) and directly probes the load-bearing assumption, a conditional verdict asking for this verification is appropriate.","tokens_in":31039,"tokens_out":16894,"duration_ms":177110,"concrete_test":"Pick one independent coefficient, e.g. c_{12×34×mmm} in Eq. (4.4). Construct a second, guaranteed-complete ansatz: enumerate all monomials in the invariants of Eq. (2.18) up to the same mass dimension and little-group weights, but without imposing any Schouten or Gram-identity reductions (so the set spans the full polynomial space up to that degree). Fit this enlarged ansatz to the same p-adic numerical data used in the paper, or to a fresh set of 10^4 random phase-space points on the variety of Eq. (2.16), including p-adic limits approaching every branch of the primary decompositions in Eqs. (2.31)–(2.48). Then compare the enlarged-fit expression and the published coefficient at another 10^4 independent finite-field points. If the values agree at every point, the minimal ansatz was complete and the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that the reconstructed expressions reproduce the one-loop amplitude exactly requires that the monomial ansatz of Sec. 2.2 spans the space of numerator polynomials for every integral coefficient. The authors explicitly state that they are 'unable to obtain a complete Gröbner basis' for the invariant ring, so the ansatz is reduced 'numerically before or during the ansatz fitting'. This numerical reduction only establishes that the retained monomials are linearly independent on the sampled phase-space points; it does not prove that they span the full five-point massive variety. If the finite-field or p-adic samples happen to lie on a subvariety that is not fully covered by the primary decompositions of Sec. 2.3, a generically independent monomial could be misclassified as redundant and dropped. The resulting coefficient would then agree with the true amplitude on the sampled regions but deviate on the complementary region. The reported OpenLoops validation is described only as 'full agreement', without specifying how many points or how they are distributed, so it does not constitute a systematic completeness test. The exactness of the analytic results therefore rests on an unproven spanning assumption, a limitation the paper itself acknowledges in Sec. 2.2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytic reconstruction of the one-loop QCD amplitude for 0 to q qbar Q Qbar H, where Q denotes a massive quark, using the massive spinor-helicity (spin-spinor) formalism. The method embeds the five-point massive kinematics into an eight-point fully massless phase space while constructing the fitting ansatz directly in five-point massive variables, and combines this with iterative pole subtraction, primary decompositions of relevant ideals, and partial fraction identities. The authors provide analytic expressions for the integral coefficients, a full ultraviolet renormalization and infrared check, and a comparison with the automatic code OpenLoops. Machine-readable expressions and evaluation code are released in Fortran and Python packages, with finite-field and p-adic evaluations and GitHub Actions tests.","tokens_in":31303,"tokens_out":5129,"duration_ms":56988,"significance":"If correct, this is the first application of analytic reconstruction to amplitudes with massive external fermions in the spin-spinor formalism, extending a line of work that previously handled massive scalars and vectors. The paper is methodologically interesting for future two-loop applications, and the release of the code, the use of finite-field and p-adic phase-space points, and the independent OpenLoops comparison are notable strengths. The claimed numerical efficiency gain (about a factor of two over OpenLoops) is modest but credible. The main risk is the completeness of the monomial ansatz used in the reconstruction, which the paper itself acknowledges in Section 2.2; this is a load-bearing point for the claimed exactness of the analytic results.","major_comments":[{"comment":"The manuscript asserts that the invariant bracket set X in Eq. (2.18) is sufficient and that the fitted ansatz is minimal, but it also states that a complete Grobner basis could not be obtained and that leftover redundancies are removed numerically. Numerical rank reduction on sampled points proves linear independence only on those points, not spanning on the full five-point massive variety. If a generically independent monomial is misclassified as redundant because the samples lie on a subvariety not covered by the primary decompositions, the reconstructed coefficient would agree with the true amplitude only on the sampled region. Since the exactness of the reconstruction is the central claim, this unproven spanning assumption needs to be either removed by a completeness argument or supported by substantially stronger validation.","section":"Section 2.2, Eqs. (2.18) and (2.20)-(2.23)"},{"comment":"The validation of the reconstructed expressions is described only as 'full agreement' with OpenLoops, without specifying the number of phase-space points, their distribution, or whether they cover the relevant branches and codimension-two limits; the GitHub Actions test evaluates the coefficients at a single physical and a single finite-field point. For a rational reconstruction in which the ansatz completeness is not proven, such a sparse check is not a systematic completeness test. The authors should either report the number and nature of the validation points (including p-adic points near all denominator poles) or soften the exactness claims to 'validated at the tested phase-space points'.","section":"Section 4 and Appendix B"}],"minor_comments":[{"comment":"The display of Eq. (2.20) is difficult to parse because of the way the four ansatz combinations are laid out with 'delta' and 'times' symbols; a clearer presentation of the four cases would help the reader.","section":"Section 2.2, Eq. (2.20)"},{"comment":"The heading 'Massive Spinors inlips and ant ares' appears to be a typo for 'Massive Spinors in Lips and Antares'.","section":"Appendix B title"},{"comment":"The phrase 'we can make see that it originates' should be 'we can see that it originates'.","section":"Section 2.3, text near Eq. (2.36)"},{"comment":"Several of the printed coefficient expressions are long and difficult to verify by eye; since the authors already provide machine-readable forms, it would be helpful to state explicitly in each case which expressions are exactly reproduced in the ancillary files and which are only representative.","section":"Section 4.1-4.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within scope for JHEP and the release of reproducible code is a substantial strength. My main concern is that the abstract's exactness claim rests on an ansatz whose completeness is explicitly not proven. If the authors can supply a credible completeness argument or a systematic validation (e.g., many random points over multiple finite-field primes and p-adic points on all relevant branches), I would be happy to support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first analytic reconstruction of a one-loop amplitude with external massive fermions in the spin-spinor formalism, for q qbar t tbar H, and it mostly delivers. The authors extend their reconstruction program to massive quarks with explicit SU(2) covariance, embed the massive five-point kinematics in a massless eight-point phase space while building the ansatz in five-point space, and develop partial-fraction and primary-decomposition tools that are genuinely new. The resulting coefficient expressions are compact enough to ship, and the paper is unusually honest about the limit: the speed gain over OpenLoops is about a factor of two, not the order of magnitude seen for massless cases.\n\nThe technical core is credible. Coefficients are reconstructed from numerical evaluations of the actual amplitude via unitarity and Passarino-Veltman reduction, so nothing is circular. The shipped code is tested in finite fields, p-adics, and complex kinematics, and the full amplitude agrees with OpenLoops for all spin choices at multiple phase-space points. That is real evidence.\n\nThe soft spot is exactly where the stress-test note points. Section 2.2 states plainly that a complete Gröbner basis for the invariant ring could not be obtained, and that leftover redundancies are removed numerically before or during the fit. That means the ansatz is not proven to span the full space of numerator polynomials. If the numerical rank reduction happened to drop a generically independent monomial, the coefficients would be wrong away from the tested points. The OpenLoops comparison is described only as 'full agreement,' with no point count or coverage detail, so it is a strong sanity check but not a systematic completeness test. I don't think this is a load-bearing flaw: random finite-field and p-adic sampling makes the missing-monomial scenario unlikely, and the reconstructed coefficients are fit to true values rather than to the ansatz. But it is a real gap in the proof of exactness, and the authors should say so more sharply or close it.\n\nCitation pattern is fine. They build on their own reconstruction papers, but that is fair, and the NLO ttH literature is cited properly.\n\nWho should read this: people working on analytic amplitudes, especially anyone planning two-loop ttH or other massive-fermion processes. It deserves a serious referee. If I were refereeing I would ask for details on the OpenLoops checks and a clearer statement about ansatz completeness, but this is not a desk reject.","headline":"First spin-spinor analytic reconstruction for one-loop amplitudes with external massive fermions; credible and useful, with a real but non-fatal caveat about the incomplete Gröbner basis.","tokens_in":31740,"tokens_out":3492,"would_cite":true,"duration_ms":38009,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The one-loop amplitude for $0 \\to \\bar{q} q t \\bar{t} H$ can be written exactly in compact massive spinor-helicity form.","keywords":["analytic reconstruction","massive spinor-helicity","one-loop amplitudes","ttH production","integral coefficients","p-adic sampling","finite-field reconstruction","partial fractions"],"falsifier":"Compute the reconstructed coefficients and the full renormalized amplitude at a previously unused random phase-space point over finite fields and p-adic numbers, and compare with an independent numerical evaluation from a second automatic generator; any discrepancy away from the construction branches would show the incomplete Gröbner basis left a missed redundancy.","tokens_in":30872,"feed_emoji":"⚛️","tokens_out":7367,"duration_ms":65953,"temperature":0.7,"pith_summary":"The paper shows that one-loop scattering amplitudes can be reconstructed analytically even when external particles are massive. It does this for $0 \\to \\bar{q} q t \\bar{t} H$, a subprocess of Higgs production with a top-quark pair at the LHC, writing every integral coefficient in compact massive spinor-helicity form. The calculation embeds the massive five-point kinematics in a massless eight-point phase space for numerical probing, while building the ansatz directly in five-point kinematics to avoid over-parameterization. Partial-fraction decompositions and common numerator factors are identified through primary decompositions in limits where pairs of denominators vanish. The resulting expressions evaluate about twice as fast as automatic numerical code, and the method is aimed at two-loop amplitudes where numerical efficiency is more important.","feed_headline":"One-loop ttH amplitudes reconstructed analytically","feed_subtitle":"A massive spinor-helicity ansatz writes the full amplitude in compact form and runs twice as fast as automatic code.","key_machinery":"The central object is the massive spinor-helicity (spin-spinor) formalism combined with analytic reconstruction in redundant variables. Massive fermion momenta are written as pairs of massless spinors with $SU(2)$ little-group indices, and the five-point massive kinematics is embedded in a fully massless eight-point phase space so that points over finite fields and $p$-adic numbers can be generated; the monomial ansatz for numerators is nevertheless built directly from five-point invariants to keep it minimal. Reconstruction proceeds by iterated pole subtraction: univariate Thiele interpolation over finite fields fixes denominator factors, $p$-adic evaluations near codimension-two varieties reveal partial-fraction structure, and primary decompositions of denominator ideals determine when numerators factor. This yields compact forms for coefficients such as $c_{13\\times 24}^{m0m}$ and organizes pole structures involving the Gram determinants $\\Delta_{12|3|4|5}$ and $\\Delta_{12|34|5}$.","core_discovery":"The paper establishes that the one-loop amplitude for $0 \\to \\bar{q} q t \\bar{t} H$ can be expressed exactly as compact analytic functions in the massive spinor-helicity (spin-spinor) formalism, with the $SU(2)$ little-group covariance of the top-quark spin states made explicit. Box, triangle, and bubble integral coefficients are reconstructed, and after subtracting the divergent box integrals the remaining triangle coefficients are proportional to the tree amplitude, so the infrared structure is fully controlled. Ultra-violet renormalization and the known pole structure are imposed explicitly, and the complete amplitude agrees with an automatic numerical evaluation at all tested phase-space points and spin choices.","pith_inferences":["The factor-of-two speedup is probably not the ceiling: the same authors' massless reconstructions achieved larger gains, so further primary decompositions and numerator-factor extractions could narrow the gap for massive processes.","The incomplete Gröbner basis is the fragile point; a useful robustness check would be to evaluate the final amplitude on a dense random sample distributed across the full phase space, including points where the numerical redundancy removal was not exercised.","The equal-mass constraint on the two heavy quarks was imposed to reduce redundancies; relaxing it should be feasible and would make the method apply to processes with distinct masses such as $t\\bar{t}W$ or $t\\bar{t}Z$ production.","The branch-dependent pole-order analysis used here is a general ansatz-reduction tool: it could be automated to shrink reconstruction problems in other multi-scale amplitudes without deriving every primary decomposition by hand."],"forward_implications":["The complete one-loop coefficient set for the $\\bar{q}q t \\bar{t}H$ channel is now available in compact analytic form; evaluating the supplied code reproduces the amplitude at any phase-space point and spin choice.","The analytic expressions are about twice as fast to evaluate as automatic one-loop generators for this process.","Because the reconstruction is built at five-point level while sampling in the embedded eight-point space, the same pipeline can be applied to the other colour structure and to $gg \\to t \\bar{t} H$, the remaining subprocess of $t\\bar{t}H$ production.","After subtraction of infrared-divergent boxes, all subtracted triangle coefficients are proportional to the tree amplitude, which fixes the infrared behaviour of the one-loop amplitude in a compact form.","The method opens a route to two-loop amplitudes with massive external quarks, where ansatz size and evaluation speed are limiting factors."],"supporting_citations":[{"why":"Supplies univariate Thiele interpolation over finite fields used to match denominator factors during reconstruction.","marker":"[1]"},{"why":"Provides the iterative pole-subtraction strategy that the whole reconstruction pipeline is built on.","marker":"[6]"},{"why":"Supplies the p-adic and algebraic-geometric ansatz framework, including the embedding variety used to generate valid phase-space points.","marker":"[18]"},{"why":"Defines the spin-spinor massive formalism used to keep explicit SU(2) covariance.","marker":"[19]"},{"why":"Provides the massive-quark helicity amplitude conventions adopted for the spin-spinors.","marker":"[20]"},{"why":"Supplies the minimal-parametrization and primary-decomposition techniques used to derive the partial-fraction forms.","marker":"[23]"},{"why":"Provides the conjectured decomposition, later proven here, that controls numerators of terms with Gram-determinant poles.","marker":"[25]"},{"why":"The automatic numerical generator used as the benchmark for correctness and evaluation speed.","marker":"[47]"}],"fun_headline_variants":["Massive spinor-helicity yields compact one-loop ttH amplitudes","Analytic reconstruction compresses one-loop ttH amplitudes","One-loop ttH amplitudes get compact analytic form","Reconstructed one-loop ttH amplitudes run faster","Compact spinor-helicity form for one-loop ttH amplitudes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reconstruction assumes that the monomial ansatz, built without a complete Gröbner basis, already contains every numerator polynomial that can appear, with the leftover redundancies removed numerically rather than by proof.","fun_headline_variants_meta":{"raw":{"variants":["Massive spinor-helicity yields compact one-loop ttH amplitudes","Analytic reconstruction compresses one-loop ttH amplitudes","One-loop ttH amplitudes get compact analytic form","Reconstructed one-loop ttH amplitudes run faster","Compact spinor-helicity form for one-loop ttH amplitudes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001463,"raw_usage":{"total_tokens":5844,"prompt_tokens":863,"completion_tokens":4981,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":4900}},"tokens_in":479,"tokens_out":4981,"duration_ms":34502,"temperature":1.0,"reasoning_tokens":4900,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:39:50.660601+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the reconstructed coefficients and the full renormalized amplitude at a previously unused random phase-space point over finite fields and p-adic numbers, and compare with an independent numerical evaluation from a second automatic generator; any discrepancy away from the construction branches would show the incomplete Gröbner basis left a missed redundancy.","supporting_citations":[],"review_version":1}