{"id":"73750d2d-c1d5-4b51-ba2f-210d483f9eae","arxiv_id":"1909.01361","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New 1/m_t^2 expansion terms, up to 1/m_t^8 for the box form factors and 1/m_t^14 for the triangle form factor, are computed for the three-loop gg to HH amplitude.","lead":"This paper computes new three-loop quantum corrections to the production of two Higgs bosons from gluons, assuming a very heavy top quark. The analytic results are building blocks for more precise LHC predictions of Higgs pair production and for measurements of the Higgs self-coupling.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No independent check anchors the new three-loop box coefficients; pole cancellation and self-convergence cannot catch a finite-term error.","rationale":"The reader's ACCEPT verdict is reasonable: the calculation is internally consistent and unusually well documented. However, the weakest_assumption identified by the reader—the domain of validity of the large-m_t asymptotic expansion—is not the load-bearing issue for the central claim. The expansion coefficients can be exactly correct even where the series is not convergent, and the paper's stated purpose is to supply input for Padé-type approximations. The load-bearing assumption is computational correctness of the new finite terms. The paper contains no independent validation of these terms, and the internal checks are insensitive to errors in finite pieces. For that reason I propose moving the verdict to CONDITIONAL: accept only if an independent computation reproduces at least one non-trivial new coefficient. This is a verification gap, not evidence of error; if the proposed check passes, ACCEPT would be appropriate. I mark agreement_with_reader as disagree because the reader's stated weakest assumption is about the asymptotic domain, whereas the concern I find load-bearing is the absence of an independent anchor for the new finite coefficients.","tokens_in":15545,"tokens_out":14580,"duration_ms":147166,"concrete_test":"Independently recompute one non-trivial new box coefficient from resFF.m—for concreteness the 1/m_t^2 term of F_box1^(2) for the colour structure C_F^2, or the C_F n_h structure if colour algebra makes that simpler—using the same 1/m_t asymptotic expansion but a different reduction path: generate the three-loop diagrams, apply the expansion with exp, reduce with LiteRed or Reduze instead of FIRE, and evaluate the resulting master integrals analytically or numerically at one Euclidean phase-space point. Compare the finite constant and logarithmic pieces with resFF.m. Agreement to numerical precision would retire the concern; any mismatch in a finite term would invalidate the new box coefficients as presented.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the five expansion coefficients of the new three-loop box form factors F_box1^(2) and F_box2^(2), up to 1/m_t^8, as provided in resFF.m, are correct. The checks reported for these genuinely new quantities are (1) cancellation of IR poles after the Catani subtraction of Section 4.2 and (2) apparent convergence of the 1/m_t expansion in Figures 3 and 4. Check (1) only constrains the divergent parts of the amplitude; the finite terms containing zeta(3), pi^4, and Li2 constants are untouched. Check (2) compares truncations of one and the same calculation, so a systematic error in the 14,580 projection tasks or in the FIRE reduction tables described in Section 3 would shift all curves coherently and remain invisible. The one- and two-loop comparisons with exact results in Figure 3 validate the method through two loops, but the two-loop box form factors are not compared with exact results in this paper, and the three-loop box coefficients are entirely new. Thus the load-bearing assumption is that the very large, unpublished computational pipeline is error-free; this is the least secure part of the argument.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the three-loop (NNLO) virtual form factors for gg -> HH in the large top quark mass limit. Using qgraf, q2e/exp, FORM-based expansions, and FIRE reduction, the authors obtain five expansion terms (up to 1/m_t^8) for the two box-type form factors F_box1 and F_box2 and eight expansion terms (up to 1/m_t^14) for the triangle form factor F_tri. After UV renormalization and Catani IR subtraction, the finite expressions are provided in a Mathematica ancillary file (resFF.m). The paper also describes substantial computational optimizations, including projection of three-loop vacuum integrals onto an ansatz, graph symmetries, ArgToExtraSymbol, and gzip compression of intermediate results. Numerical plots show the convergence of the 1/m_t expansion and comparisons with exact results where available.","tokens_in":15698,"tokens_out":9296,"duration_ms":90652,"significance":"If correct, these are the first three-loop results beyond the leading power for the box-type gg -> HH form factors, and they extend the triangle form factor by several powers. The results are useful input for NNLO approximations of Higgs-pair production, including Padé-based constructions and soft-virtual approximations. Strengths of the paper include: no parameters are fitted to the target result; the computational pipeline is described in unusual detail; machine-readable analytic results are supplied; and the method is validated at one loop (all form factors) and at two loops (triangle form factor) against exact results. The technical optimizations described are likely to benefit future multiloop calculations in the same framework.","major_comments":[{"comment":"The three-loop box-type form factors are entirely new, and the only reported checks for them are cancellation of IR poles after the Catani subtraction and the apparent convergence of the same 1/m_t expansion. These checks do not constrain the finite transcendental parts of the new coefficients. Moreover, while the two-loop triangle form factor is compared with the exact expression (text below Fig. 3), no comparison with exact results is shown for the two-loop box form factors, so the validation chain skips the complexity level closest to the new three-loop box calculation. Please add an independent check, for example a comparison of the two-loop box form factors with the exact NLO results from Refs. [4-6] at representative kinematic points, or a numerical evaluation of a subset of the new three-loop master integrals. If no independent check is feasible, the conclusions should explicitly state that the three-loop box coefficients currently rest on internal-consistency checks alone.","section":"Section 5, Figs. 3-4 and Section 4.2"}],"minor_comments":[{"comment":"The phrase 'We present analytic results for the form factors in the soft-virtual approximation' is misleading: the results contain full dependence on s, t and u through the logarithms and dilogarithms defined in Section 5, not only the soft-virtual limit. Please reword to say 'in the large top mass limit' or 'with full kinematic dependence within the large-mt expansion'.","section":"Abstract"},{"comment":"There is a typo: 'three-loop for factors' should read 'three-loop form factors'.","section":"Section 5, first paragraph"},{"comment":"The expression '1/m12 t ' appears to be a typesetting issue; it should be 1/m_t^12.","section":"Introduction, paragraph 1"},{"comment":"Since the exact two-loop comparison is made only for the triangle form factor, the caption should state this explicitly so that readers do not infer that all two-loop curves are checked against exact results.","section":"Figure 3 caption"},{"comment":"The full results are not printed in the manuscript but reside in an external ancillary file. Please ensure that the file is included as a permanent ancillary file with the published version (e.g., arXiv ancillary file or journal supplementary material) and that Ref. [49] remains accessible; also state the file format and any restrictions on use.","section":"Appendix A and Ref. [49]"}],"recommendation":"major_revision","confidential_remarks":"This is a strong technical paper and the calculation is likely correct, but the new three-loop box coefficients are the central product of the work and they lack any external anchor. I would support acceptance once the authors either provide a concrete independent check (for example, using the exact NLO amplitudes from the literature) or clearly state this limitation in the conclusions. The ancillary-file reproducibility issue should also be resolved before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nBottom line: this is a solid technical paper that produces genuinely new three-loop expansion coefficients for gg→HH form factors in the large-m_t limit, and it deserves a serious referee. The headline claim is exactly what it appears to be: five 1/m_t^2 expansion terms for the two box-type form factors and eight for the triangle, with the triangle going beyond the known 1/m_t^8. The ancillary file makes the results usable. The projection-operator ansatz in Section 3.2 is a real optimization, and the description of the computation is unusually detailed for this genre.\n\nThe checks are reasonable as far as they go. IR poles cancel after Catani subtraction; one-loop box and two-loop triangle agree with exact results where those exist; the 1/m_t convergence plots look plausible. I agree with the reader's ACCEPT verdict. The stress-test note identifies the genuine soft spot: for the genuinely new three-loop box coefficients, the only anchors are pole cancellation and self-convergence of the same series. Pole cancellation does not constrain finite zeta(3), pi^4, and Li2 terms, and self-convergence cannot catch a systematic error that shifts all truncations coherently. That is a real limitation, though it is the standard limitation of large automated loop calculations. It doesn't undermine the claim; it just means the coefficients should be treated as 'computed, not yet independently confirmed.' If I were referee I'd ask the authors to provide one extra numerical cross-check, e.g., a single coefficient evaluated with an independent reduction or a different expansion strategy, or at least make the stored intermediate expressions available.\n\nCitation pattern looks appropriate. The authors clearly state which terms are new and which come from Refs. [31,32], and they point to Ref. [17] as already using the triangle results. No concerns there.\n\nWho is this for? Anyone working on NNLO Higgs pair production approximations, especially Padé/hybrid constructions in the spirit of Gröber-Maier-Rauh. It's a technical input paper, not a phenomenological breakthrough, but it is exactly the kind of result that downstream approximations rely on. I'd take it to reading group, and I'd cite it if I worked on this process.\n\nRecommendation: send it to peer review. It's a legitimate calculation with a clear statement of what's new, and the main open question—independent confirmation of the box coefficients—is a matter for the community, not a reason to desk reject.","headline":"Genuinely new three-loop large-top-mass form factor coefficients for gg→HH, solidly presented with the usual caveat that the new box-type terms lack independent confirmation.","tokens_in":16272,"tokens_out":1972,"would_cite":true,"duration_ms":19619,"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":"This paper computes five expansion terms for the three-loop box-type form factors and eight for the triangle form factor in large-top-mass Higgs boson pair production, giving analytic results in the soft-virtual approximation.","keywords":["Higgs boson pair production","three-loop form factors","large top mass limit","asymptotic expansion","NNLO corrections","gluon fusion","soft-virtual approximation","analytic results"],"falsifier":"One concrete test is to compute the next order in the expansion, the $1/m_t^{10}$ box term (or the $1/m_t^{16}$ triangle term), and check that it is suppressed relative to the deepest computed term by roughly the external squared momentum divided by $m_t^2$; failure of that suppression would show the series is not converging where it is used. A future exact numerical evaluation of the three-loop virtual amplitude with finite top mass at center-of-mass energies near 300 to 400 GeV would settle the same question directly.","tokens_in":15308,"feed_emoji":"⚛️","tokens_out":17188,"duration_ms":158373,"temperature":0.7,"pith_summary":"Higgs boson pair production via gluon fusion is the main way to probe the Higgs self-coupling, but at the three-loop level the full amplitude is too hard to compute exactly. This paper extends the large-top-mass expansion of the three-loop virtual amplitude, computing five expansion terms (through $1/m_t^{8}$) for the two box-type form factors and eight terms (through $1/m_t^{14}$) for the triangle form factor. The triangle form factor is the part tied to the triple-Higgs coupling, which makes the deeper triangle expansion directly relevant to the physics goal. The paper presents finite, analytic expressions after ultraviolet renormalization and infrared subtraction, and explains the computational restructuring that made the calculation feasible.","feed_headline":"Three-loop Higgs-pair amplitude gains deeper top-mass expansions","feed_subtitle":"Deeper analytic expansions for box and triangle form factors improve Higgs-pair predictions at NNLO.","key_machinery":"The argument rests on two techniques. The first is the asymptotic expansion in $1/m_t^2$ defined by taking $m_t$ much larger than every external momentum: it splits each three-loop diagram into hard subgraphs, which are expanded in the external momenta and the co-subgraph loop momenta, and co-subgraphs, reducing the integrals to products of vacuum integrals and massless one- or two-loop integrals with at most one additional scale. The second is a projection method: each diagram is written as a polynomial in the external scalar products, and derivative operators extract the coefficients one at a time, so that three-loop tensor vacuum integrals never have to be computed directly. Applying the derivatives after the $1/m_t^2$ expansion, and storing the expanded super-diagrams in compressed form, keeps the huge intermediate expressions manageable. This is what allows the fifth-order box expansion and the eighth-order triangle expansion to be obtained.","core_discovery":"The central claim is that the three-loop amplitude for Higgs boson pair production via gluon fusion, in the limit where the top quark mass is much larger than all external momenta, can now be given analytically to order $1/m_t^{8}$ for the two box-type form factors $F_{box1}$ and $F_{box2}$, and to order $1/m_t^{14}$ for the triangle form factor $F_{tri}$. The box-type results are new; the triangle result extends the previously known expansion to three more orders. The finite forms are obtained by ultraviolet renormalization and by subtracting infrared poles with a standard prescription, leaving finite form factors expressed in terms of the colour factors, the number of light and heavy quark flavours, and the kinematical variables. The analytic results are presented for the soft-virtual approximation and supplied in computer-readable form for direct use in approximation procedures.","pith_inferences":["If the three-loop series converges like the one- and two-loop series shown below threshold, the new terms could make rational-approximant predictions trustworthy well above the two-top-quark threshold, where no exact three-loop result exists; this goes beyond what the paper explicitly demonstrates.","The observation that one box-type form factor starts only at order $1/m_t^2$ for most colour structures suggests the box corrections are more suppressed than the triangle corrections, so the effective-theory uncertainty in NNLO approximations may come mainly from the triangle sector.","A useful stress test would be comparing the deepest box expansion against the approximate NNLO cross-section predictions after the new terms are included, to see whether the size of the shift decreases with each added order.","If the projection method's efficiency holds, the same approach could target the next $1/m_t^{10}$ box terms rather than waiting for exact finite-top-mass integrals."],"forward_implications":["The new box-type form-factor terms can be inserted into existing approximation schemes that combine exact next-to-leading-order results with effective-theory NNLO building blocks, improving their kinematic coverage.","The deeper triangle expansion provides additional input for rational-approximant constructions of the full top-mass dependence, a route the paper notes has already been used for the related Higgs-gluon form factor.","Because the curves for the deepest expansions lie close together, rescaling the expanded higher-order corrections by the exact leading-order ratio should give numerically stable approximations across the plotted phase space.","The projection technique for bypassing three-loop tensor vacuum integrals can be reused in other multiloop calculations with a large internal mass, reducing both memory and CPU requirements."],"supporting_citations":[{"why":"Earlier work that computed expansion terms up to $1/m_t^4$ for the partonic cross section; this paper extends that to the individual form factors.","marker":"[21]"},{"why":"Provides the one- and two-loop form factors used for comparison and for the rescaling approximation.","marker":"[12]"},{"why":"Supplies the rational-approximant method for turning expansion terms into amplitudes valid over a wider kinematic range.","marker":"[16]"},{"why":"Uses the triangle form factor results of this paper for the related Higgs-gluon form factor.","marker":"[17]"},{"why":"Describes the approximate NNLO prediction scheme that these new virtual corrections could improve.","marker":"[23]"},{"why":"Earlier NNLO results in the heavy-top limit and part of the framework this work builds on.","marker":"[8]"},{"why":"Provides the infrared subtraction formalism used to obtain finite form factors.","marker":"[47]"},{"why":"The supplementary file containing the analytic results in computer-readable form.","marker":"[49]"}],"fun_headline_variants":["Three-loop Higgs-pair form factors reach deeper expansion orders","New analytic expansions for three-loop Higgs-pair form factors","Three-loop Higgs-pair amplitude: new top-mass expansions for form factors","Deeper top-mass expansions for three-loop Higgs-pair form factors"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the top quark is heavy enough that a power-series expansion in external momenta divided by the top mass is an accurate representation of the amplitude in the kinematic region where the form factors are applied, which includes energies above the two-top-quark threshold.","fun_headline_variants_meta":{"raw":{"variants":["Three-loop Higgs-pair form factors reach deeper expansion orders","New analytic expansions for three-loop Higgs-pair form factors","Three-loop Higgs-pair amplitude: new top-mass expansions for form factors","Deeper top-mass expansions for three-loop Higgs-pair form factors"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00095,"raw_usage":{"total_tokens":3981,"prompt_tokens":803,"completion_tokens":3178,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":419,"completion_tokens_details":{"reasoning_tokens":3106}},"tokens_in":419,"tokens_out":3178,"duration_ms":25897,"temperature":1.0,"reasoning_tokens":3106,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:20:05.905997+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One concrete test is to compute the next order in the expansion, the $1/m_t^{10}$ box term (or the $1/m_t^{16}$ triangle term), and check that it is suppressed relative to the deepest computed term by roughly the external squared momentum divided by $m_t^2$; failure of that suppression would show the series is not converging where it is used. A future exact numerical evaluation of the three-loop virtual amplitude with finite top mass at center-of-mass energies near 300 to 400 GeV would settle the same question directly.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The supplementary file containing the analytic results in computer-readable form."}],"review_version":1}