{"id":"efaa6573-ed95-4e14-a86a-eaea3081fe5e","arxiv_id":"1908.09932","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The thesis computes three-loop and two-loop QCD corrections to Higgs processes with full quark mass dependence, including a new analytic series-expansion treatment of elliptic two-loop master integrals for Higgs-plus-jet production.","lead":"A particle physics dissertation computes very high-order quantum chromodynamics corrections to Higgs boson processes that run through heavy quark loops. It delivers new results for three Higgs processes, including a new technique for the hardest two-loop Higgs-plus-jet integrals.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The matched-series evaluation of the elliptic master integrals is the load-bearing new claim, but the supplied text omits the numerical checks in Chapter 7, leaving the convergence/coverage assertion unverified.","rationale":"The reader's verdict is CONDITIONAL with low confidence, and the weakest assumption is identified as the convergence and matching coverage of the one-dimensional series-expansion method for the elliptic master integrals. My stress-test arrives at the same load-bearing point. The two non-elliptic results in the thesis are independently supported: the H→Zγ two-loop result is an analytic re-derivation confirming published numerics, and the three-loop Hb bbar form factor matches the predicted infrared pole structure. These parts do not raise a correctness concern. The genuinely novel and consequential claim is the numerical evaluation of elliptic multi-scale master integrals via matched series expansions. The supplied text contains the method description (Chapter 6) and the outline of its application (Chapter 7 table of contents and forward references), but the actual numerical checks—Section 7.9 including degree of series expansions, timings, relative deviation, and truncation error—are truncated away. Without those checks, or without accompanying code and data, the claim that the expansions cover the whole physical phase space to the required precision cannot be independently verified. This is a missing-evidence concern, not a demonstrated inconsistency; hence the appropriate verdict remains CONDITIONAL rather than ACCEPT or REJECT. A concrete, decisive test is to reproduce the matching for a representative elliptic sector and compare against independent high-precision numerics across the physical phase space, including overlap and boundary regions. This test directly targets the weakest assumption and would settle whether the concern lands.","tokens_in":61268,"tokens_out":4605,"duration_ms":48349,"concrete_test":"Implement the series-expansion matching for one representative elliptic MI (e.g., sector A6,215 from Appendix C.2) along the one-dimensional parameter λ, using the code/data from the thesis if available. Then: (1) compute the radius of convergence of each expansion from the nearest singularity in λ for a grid of (x,z,h) values spanning the physical region; (2) evaluate the matched series at a dense set of points, especially in overlap regions and near the endpoints, and compare against an independent high-precision evaluation (e.g., pySecDec sector decomposition or adaptive numerical integration of the Feynman-parameter representation). If any test point shows relative deviation larger than the truncation error claimed in Sections 7.9.3–7.9.4, the matching coverage fails and the central claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new claim is that the planar two-loop master integrals for Higgs-plus-jet production, including the elliptic sectors A6,215 and A7,247, can be evaluated numerically 'in a fast and reliable way' by series expansions derived from differential equations and matched across singular points (English Summary; Chapter 7, Sections 7.5.5, 7.6.5, 7.9.3–7.9.4). The supplied manuscript truncates before Section 7.9, so the numerical evidence for this claim is absent. The method description itself shows why this is load-bearing: the matching procedure (Sections 6.5–6.6) partitions the λ-phase space at singular points and connects local series solutions. A local series around a singular point has radius of convergence limited by the nearest other singularity in the λ-plane. For the elliptic sectors the singularity positions depend on the auxiliary kinematic invariants (x,z,h in Eq. (4.15)); as those invariants approach thresholds, two singular points can coalesce, shrinking the convergence disk to zero and leaving gaps in the claimed coverage. The thesis asserts (Sections 7.9.3–7.9.4) that overlapping intervals cover the whole physical phase space, but neither the convergence radii nor the numerical checks are shown in the provided text. This is not an accusation of error; it is a claim that the central novel result is currently unsupported by the available evidence. The H→Zγ analytic re-derivation and the three-loop Hb bbar form factor have independent checks (earlier numerics, IR pole structure), but the elliptic MI evaluation has no such external anchor in the supplied material.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This thesis-style manuscript presents three higher-order QCD calculations involving Higgs bosons and heavy-quark loops. Chapter 3 derives the three-loop QCD corrections to the Hb bbar form factor in the massless limit and verifies the infrared pole structure against known factorization formulae. Chapter 5 computes the two-loop QCD corrections to H->Z gamma with full quark mass dependence, analytically confirming earlier purely numerical results obtained with on-shell renormalization. Chapters 6 and 7 develop and apply a method to evaluate planar two-loop master integrals for Higgs-plus-jet production with full quark mass dependence, using series expansions derived from differential equations because the master integrals involve elliptic structures. The abstract claims that these elliptic integrals can be evaluated numerically in a fast and reliable way by matching multiple series expansions over the physical phase space.","tokens_in":61526,"tokens_out":5397,"duration_ms":54264,"significance":"If the matched-series method works as claimed, it represents a significant technical advance: exact top-quark-mass dependence for the two-loop Higgs-plus-jet amplitude is a recognized high-priority target, and a practical method for elliptic multi-scale master integrals would have broader applicability. The H->Z gamma analytic computation and the three-loop form factor provide strong internal consistency checks: the former reproduces known numerics from Ref. [252], and the latter matches the predicted infrared pole structure (Section 3.4, Eqs. (3.17)-(3.19)). The method introduces no fitted parameters; the only free choices are the renormalization scale and the series truncation order. However, the central new claim—the fast and reliable numerical evaluation of the elliptic master integrals—is not verifiable from the text supplied for review, because the numerical validation sections are absent.","major_comments":[{"comment":"The load-bearing claim that the planar two-loop master integrals for Higgs-plus-jet production, including the elliptic sectors A6,215 and A7,247, can be evaluated 'in a fast and reliable way' by matched series expansions is not supported by the evidence in the available manuscript. The table of contents lists the numerical checks (degree of the series expansions, timings, relative deviation, truncation error), but the supplied text truncates before these sections, so no convergence radii, overlap tests, or accuracy numbers are shown. This is not a cosmetic issue: the matching procedure in Sections 6.5-6.6 connects local series solutions only where their convergence disks overlap, and the claimed coverage of the whole physical phase space requires numerical demonstration. Please provide the missing validation or state explicitly that it is deferred to a separate publication.","section":"Sections 7.9.1-7.9.4 (with 6.5-6.6 and 7.5.5, 7.6.5)"},{"comment":"The partitioning of the phase space through singular points is asserted to give overlapping intervals covering the full physical region, but the argument is incomplete. A local series around a singular point lambda_0 has radius of convergence limited by the nearest other singularity in the lambda-plane; as the auxiliary kinematic invariants x, z, h defined in Eq. (4.15) approach thresholds, two singular points can coalesce and shrink the convergence disk to zero. The manuscript does not provide a bound on the distances between adjacent expansion centers relative to these radii, nor an explicit scan of the (x, z, h) parameter region near thresholds. I request either a proof of gap-free coverage or a numerical demonstration (for example, a grid scan of relative deviations) for the elliptic sectors.","section":"Sections 6.6.1 and 7.5.2"},{"comment":"The text supplied for review ends during Section 5.3 and does not include the announced Sections 5.4 (numerical results), 6, 7.6-7.10, or the appendices. Consequently, the claimed computation of the two-loop amplitude in terms of master integrals (Section 7.3.2), the treatment of the elliptic sectors A6,215 and A7,247 (Sections 7.6-7.7), and the numerical checks (Section 7.9) cannot be checked. If this manuscript is intended for journal publication, the complete derivations and validation must be included; a thesis may be a self-contained document, but the submitted excerpt is not.","section":"Sections 5.4-5.5 and 7.6-7.10"}],"minor_comments":[{"comment":"The sentence discussing the Lorentz structure of the external momenta contains the duplicated article in 'the the Lorentz structure'; this should be corrected.","section":"Section 2.1.1"},{"comment":"The word 'tapole' appears in the text and should read 'tadpole'.","section":"Section 4.4.3"},{"comment":"The three-loop form factor expressions are long; an electronic ancillary file with the Laurent coefficients in machine-readable form would substantially aid verification and reuse.","section":"Eqs. (3.14)-(3.16)"},{"comment":"The bibliography is not included in the supplied text, so citations such as [252] cannot be resolved; please ensure the complete reference list is part of any revised submission.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This manuscript appears to be a doctoral thesis posted to arXiv rather than a focused research article, and much of the content overlaps with the author's prior publications (indicated as Refs. [28-30]) as well as with Ref. [43]. If the journal has policies on prior publication or on the originality of submitted work, the editor should clarify the overlap with the authors. Independent of that, any revised version must contain the numerical validation of the matched-series method for the elliptic master integrals, because the current text does not support the central claim of fast and reliable numerical evaluation over the full physical phase space."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The part of this thesis you can fully check is good, and the part you cannot fully check is the part that matters most. The three-loop Hbb form factor and the analytic two-loop H→Zγ result both have independent support: the form factor reproduces the predicted infrared pole structure, and the H→Zγ derivation confirms earlier numerics. Those two sections are careful, readable, and trustworthy.\n\nThe genuinely new piece is the matched-series evaluation of the planar elliptic master integrals for Higgs-plus-jet production. The method itself is plausible and well explained: reduce the multi-scale problem to one variable, compute series expansions around singular points of the differential equations, and match them across the physical phase space. The matching idea is sensible, and the author is honest that the elliptic sectors (A6,215 and A7,247) resist the usual decoupling. What is missing is the evidence that the method works. Sections 7.9.3–7.9.4, which are supposed to show relative deviations and truncation errors, are not in the supplied text. The English Summary claims the results are evaluated “in a fast and reliable way,” but the numerical checks, the convergence radii, and the coverage arguments are absent from what I received. The stress-test concern about singular points coalescing in the λ-plane is legitimate, and without those checks the convergence claim is unsupported.\n\nI do not think this is an error; it is a hole in the available evidence. The thesis likely contains those checks in the full version, and the work is serious enough that a referee should see them. Two smaller points: no code or data is attached, which makes independent verification harder, and the pedagogical chapters (complex analysis, MPLs) could be trimmed, though they do not hurt. The citation pattern looks normal for a thesis of this type; I see no circularity.\n\nBottom line: this deserves a serious referee, but the referee should have either the complete thesis or the companion papers. If the numerical checks hold up, the elliptic method is a real contribution to multi-scale two-loop calculations. If they do not, the H→Zγ and form-factor parts still stand on their own. I would cite the H→Zγ result and would bring the thesis to a reading group as a maybe, mainly to discuss the matching strategy.","headline":"A thesis with two solid, cross-checked results and one interesting but evidentially incomplete claim about evaluating elliptic master integrals via matched series expansions.","tokens_in":62074,"tokens_out":1891,"would_cite":true,"duration_ms":23645,"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":"This thesis claims that higher-order Higgs-boson amplitudes with full quark-mass dependence can be computed, including elliptic two-loop master integrals for Higgs-plus-jet production via matched series expansions derived from…","keywords":["Higgs boson","quantum chromodynamics","heavy-quark loops","master integrals","differential equations","elliptic integrals","series expansions","Higgs-plus-jet production"],"falsifier":"Compute one of the planar elliptic master integrals at a dense grid of physical phase-space points using an independent numerical method, such as sector decomposition of the Feynman-parameter representation, and compare with the matched-series evaluation; any disagreement beyond the stated truncation error would falsify the claim of full physical-region coverage.","tokens_in":61006,"feed_emoji":"⚛️","tokens_out":7132,"duration_ms":71119,"temperature":0.7,"pith_summary":"The thesis aims to push perturbative QCD predictions for Higgs-boson processes mediated by heavy-quark loops to higher orders while keeping the full quark mass. It derives the three-loop QCD correction to the form factor for the Higgs coupling to bottom quarks in the massless-quark limit, and it analytically computes the two-loop QCD corrections to the H→Zγ decay with the full internal quark mass, confirming earlier numerical results. Its principal new method reduces multi-scale master integrals to one variable, derives series expansions around the singular points of their differential equations, and matches those expansions so that the planar two-loop master integrals for Higgs-plus-jet production—including elliptic multi-scale cases—can be evaluated numerically across the whole physical phase space. If the method is right, NLO predictions for Higgs-plus-jet production with full top-mass dependence become feasible in the high-transverse-momentum region where the usual infinite-top-mass effective theory is unreliable.","feed_headline":"Elliptic two-loop Higgs integrals computed via matched series expansions","feed_subtitle":"Full quark-mass dependence reaches NLO in a region where the usual infinite-top-mass approximation breaks down.","key_machinery":"The load-bearing mechanism is the series-expansion-from-differential-equations method presented in Chapter 6. A multi-scale problem is reduced to a single variable by a one-dimensional parametrization of the phase space; the differential equations for the master integrals are solved as series around their singular points, using homogeneous and particular solutions of second-order ordinary differential equations for the elliptic sectors; and a matching procedure joins these expansions through overlapping convergence regions, with boundary conditions fixed by regularity at pseudo-thresholds. The non-elliptic sectors are handled through canonical differential equations in d-log form, while the planar elliptic sector A6,215 is carried by the second-order equations whose homogeneous solutions are elliptic functions.","core_discovery":"The central discovery claimed is that full quark-mass dependence does not block higher-order Higgs amplitudes: the three-loop Hb bbar form factor is obtained in massless QCD with its infrared poles matching the factorization prediction; the two-loop H→Zγ amplitude is obtained analytically in terms of multiple polylogarithms with full quark-mass dependence; and the planar master integrals for two-loop Higgs-plus-jet production, including the elliptic sectors, are evaluated by a one-dimensional parametrization of the phase space, series expansions around singular points of the differential equations, and a matching procedure whose overlapping radii of convergence cover the physical region. The thesis claims the resulting numerical evaluation is fast and reliable, and that the two-loop Higgs-plus-jet scattering amplitude can be expressed in terms of these planar master integrals together with the non-planar ones.","pith_inferences":["(Editor's inference) If the matching procedure carries over to the non-planar sectors, the full NLO Higgs-plus-jet amplitude with exact top-mass dependence becomes numerically tractable, and the main remaining bottleneck would be computational, not conceptual.","(Editor's inference) The one-dimensional parametrization suggests a direct independent test: compare the matched-series values of the elliptic master integrals at selected physical points against a completely different numerical technique; the thesis's internal consistency checks are necessary but not sufficient.","(Editor's inference) The single-logarithmic small-mass behaviour of H→Zγ, in contrast to the double logarithms in H→γγ, points to a structural feature of the Z-boson coupling that could be studied separately, since the thesis reports it but does not explain it."],"forward_implications":["The analytic two-loop H→Zγ result confirms the earlier numerical computation and allows a direct study of renormalization-scheme and scale dependence of the decay width.","The three-loop Hb bbar form factor supplies a core ingredient for third-order QCD corrections to Higgs production in bottom-quark fusion and to the H→bbar decay rate.","The matched-series method evaluates the planar two-loop master integrals for Higgs-plus-jet production with full quark mass dependence in the physical region, including elliptic multi-scale integrals, in a fast and reliable way.","Because the quark mass is retained, the two-loop Higgs-plus-jet amplitude can be trusted at high Higgs transverse momentum, where the infinite-top-mass effective theory is not appropriate."],"supporting_citations":[{"why":"Supplies the previous numerical two-loop H→Zγ result that the analytic calculation confirms.","marker":"[252]"},{"why":"Provides the prior two-loop double-Higgs production calculation with full top-mass dependence in numerical form, a benchmark for the series-expansion approach.","marker":"[41,42]"},{"why":"A recent numerical two-loop Higgs-plus-jet calculation with fixed top-mass ratio, the state of the art this work extends.","marker":"[43]"},{"why":"Introduces canonical differential equations in d-log form, the framework used for the non-elliptic master integrals.","marker":"[201,202]"},{"why":"Provides the method of degenerate integration-by-parts relations in fixed dimensions used to decouple differential equations and find suitable integral bases.","marker":"[199]"},{"why":"Formulates the infrared factorization that predicts the pole structure of the three-loop form factor, used as the central check of Chapter 3.","marker":"[84,145–148]"},{"why":"Gives the multi-loop workflow and integral basis that the three-loop Hb bbar form factor calculation follows.","marker":"[116–118]"}],"fun_headline_variants":["Full quark mass dependence at three loops for Higgs form factors","Two-loop Higgs amplitudes with exact quark masses","Elliptic integrals tamed: Higgs amplitudes with full mass dependence","Series expansions crack elliptic two-loop Higgs integrals","Massive quarks in Higgs loops: higher orders computed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the series expansions around the singular points of the differential equations can be matched so that their combined radii of convergence cover every point of the physical phase space; if some region is left uncovered, the numerical evaluation of the elliptic master integrals and thus of the two-loop Higgs-plus-jet amplitude would be unreliable.","fun_headline_variants_meta":{"raw":{"variants":["Full quark mass dependence at three loops for Higgs form factors","Two-loop Higgs amplitudes with exact quark masses","Elliptic integrals tamed: Higgs amplitudes with full mass dependence","Series expansions crack elliptic two-loop Higgs integrals","Massive quarks in Higgs loops: higher orders computed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000707,"raw_usage":{"total_tokens":3133,"prompt_tokens":843,"completion_tokens":2290,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":459,"completion_tokens_details":{"reasoning_tokens":2214}},"tokens_in":459,"tokens_out":2290,"duration_ms":20072,"temperature":1.0,"reasoning_tokens":2214,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:57:45.194254+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute one of the planar elliptic master integrals at a dense grid of physical phase-space points using an independent numerical method, such as sector decomposition of the Feynman-parameter representation, and compare with the matched-series evaluation; any disagreement beyond the stated truncation error would falsify the claim of full physical-region coverage.","supporting_citations":[],"review_version":1}