{"id":"44154fc5-0632-433f-a4ca-515f44a9b86a","arxiv_id":"2412.06881","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"The three-loop inclusive bHQET jet function for boosted heavy quarks is computed analytically, completing the fixed-order ingredients for N3LL' top-mass observables.","lead":"This paper computes the three-loop quantum chromodynamics correction to the jet function that describes narrow radiation inside a boosted heavy quark. The result is the final missing ingredient for more precise top quark mass predictions at future lepton colliders.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The one unshown, load-bearing input is the non-planar master integral MI3L_t and the unpublished B_bare^3/MI combination; a systematic error there would change Eq. (4.7).","rationale":"The paper's central claim is a three-loop coefficient; all other results (mass scheme, R-evolution, four-loop estimate) are derived or secondary. The computation is a conventional Feynman-integral calculation: qgraf to Looping to FIRE to master integrals. The final result passes impressive checks: gauge-parameter independence, non-Abelian exponentiation in d dimensions, reproduction of Gamma_c^2 and gamma_B^2, and agreement with the n_l^2 renormalon prediction. These make a gross error unlikely, and I do not believe the claim is circular or wrong for an identifiable reason. However, the path from the 20 MIs to Eq. (4.7) is not publicly reproducible: B_bare^3 is given only on request, and the one non-planar MI is stated as a single 3F2 line without derivation. The numerical FIESTA/pySecDec checks in Appendix F are good evidence, but they evaluate the listed MIs; they do not independently establish the analytic all-orders expression for MI3L_t, nor the algebraic reduction that assembles it into b30. Because b30 is the load-bearing output, the residual risk is concentrated there. This matches the reader's weakest_assumption. My recommendation is UNCHANGED (CONDITIONAL): no new defect is established, but the conditionality due to unshipped artifacts and the unshown MI-level reduction should remain until an independent check of MI3L_t or release of the B_bare^3 combination.","tokens_in":47820,"tokens_out":11968,"duration_ms":125133,"concrete_test":"Independently re-derive MI3L_t at d=4-2*epsilon using a second Mellin-Barnes representation (e.g., starting from the Feynman-parameter form after integrating two loops, using a different contour package) and compare the resulting epsilon-expansion through O(epsilon^3) with the expansion of Eq. (F.8u). As a second leg, request the B_bare^3/MI combination of footnote 15 from the authors and substitute the Appendix F expansions symbolically to confirm Eq. (4.7); if either check changes b30 by more than 0.1%, the central number needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new number, Eq. (4.7), is obtained by inserting 20 analytic three-loop master integrals into a linear combination B_bare^3 that is not displayed (footnote 15). Among the MIs, the single most vulnerable entry is the non-planar MI3L_t: Eq. (F.8u) is presented as a closed 3F2 result, but Appendix G only sketches the Mellin-Barnes procedure and does not show the contour manipulations, residue sums, or the analytic continuation that produces the all-orders expression. The planar MIs received additional cross-checks (IBP/dimensional recurrence, quasi-finite/HyperInt), and the final result is checked against gauge invariance, non-Abelian exponentiation, and known anomalous dimensions; those checks are real but they are consistency checks on the whole combination, not an independent verification of MI3L_t. The quasi-finite integrals F.8n-F.8s are likewise expanded only to O(epsilon) or O(epsilon^2), and the paper states without showing the B_bare^3 decomposition that this is 'one order higher than necessary.' A wrong prefactor or a missed epsilon term in MI3L_t would propagate directly into the C_F C_A^2 part of b30, which is currently the largest and least cross-checked piece of the numerical result. This is a verification gap rather than a demonstrated error.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first three-loop (O(alpha_s^3)) computation of the inclusive bHQET jet function for boosted heavy quarks. The calculation is performed by direct Feynman-diagram evaluation, with IBP reduction to 20 master integrals, followed by analytic evaluation using Mellin-Barnes, integration-by-parts, and quasi-finite/HyperInt methods. The main quantitative result is Eq. (4.7), the three-loop non-logarithmic coefficient b30 of the exponent of the position-space jet function, numerically b30 = 50.054 n_l^2 - 1899.8 n_l + 12834 for N_c=3. The authors verify the result against a number of independent constraints: gauge-parameter independence, non-Abelian exponentiation in d dimensions, reproduction of the known cusp and non-cusp anomalous dimensions, and agreement with the large-beta_0 prediction for the n_l^2 term. They also use the result to derive the three-loop pole-to-jet-mass scheme conversion, compare two short-distance jet-mass schemes with the MSR mass, and estimate the four-loop non-logarithmic coefficient using renormalon dominance. The paper is careful and detailed, but one load-bearing piece of information is not shown: the decomposition of the bare three-loop matrix element B_bare^3 into the master integrals, and the derivation of the non-planar master integral MI3L_t is only sketched.","tokens_in":48053,"tokens_out":4574,"duration_ms":51564,"significance":"If correct, this is an important result. The three-loop bHQET jet function is the last missing perturbative ingredient for N3LL'-accurate self-normalized thrust-type distributions in the peak region of boosted top production, and it also provides the first direct determination of the three-loop non-cusp anomalous dimension of this jet function through the renormalization-group consistency of the calculation. The paper benefits from unusually strong internal and external consistency checks: the gauge-parameter independence is checked, non-Abelian exponentiation is verified before epsilon-expansion, the anomalous dimensions are reproduced, the n_l^2 term agrees with the large-beta_0 prediction, and all master integrals are checked numerically with FIESTA/pySecDec. These checks make an outright error in Eq. (4.7) unlikely, but they do not remove the need for the transparency concerns raised below.","major_comments":[{"comment":"The central new number, Eq. (4.7), is obtained by inserting the 20 three-loop master integrals of Appendix F into the linear combination B_bare^3, but the expression for B_bare^3 in terms of MIs is not given; footnote 15 states that it can be obtained from the authors upon request. Because an error in this combination would propagate directly into b30, this is a load-bearing transparency gap. I request that the full d-dimensional expression for B_bare^3 (or at least the epsilon-expanded coefficients needed for the renormalized result) be provided as an ancillary file or in an appendix, in a computer-readable form.","section":"Sec. 4.3 and footnote 15"},{"comment":"The non-planar master integral MI3L_t is the least cross-checked entry in the C_F C_A^2 part of b30, which is the largest and least independently verified piece of the numerical result. Eq. (F.8u) is presented as a closed all-orders 3F2 result, but Appendix G.2 only sketches the Mellin-Barnes procedure and does not show the contour choices, residue sums, or analytic continuation for this specific integral. Since the numerical checks in Appendix F are strong but not a proof, I ask the authors to provide a more complete derivation of MI3L_t, or an independent verification such as an evaluation through the quasi-finite/HyperInt method, so that this load-bearing ingredient can be audited.","section":"Appendix G.2 and Eq. (F.8u)"}],"minor_comments":[{"comment":"The quasi-finite integrals F.8n through F.8s are expanded only to O(epsilon) or O(epsilon^2), with the statement that this is one order higher than necessary. It would be helpful to state explicitly, once B_bare^3 is provided, which epsilon order is required from each master integral, so that readers can verify that the quoted truncations are sufficient.","section":"Appendix G.4"},{"comment":"The four-loop estimate relies on the assumption that the u=1/2 renormalon dominates b40. The internal consistency tests in Table 2 and Fig. 3 are reassuring, but the text could state more explicitly that this is an assumption rather than a derivation, and that the resulting uncertainties are only as reliable as that assumption.","section":"Sec. 5"},{"comment":"The caption uses \"Revolver\" while the reference and the main text use \"REvolver\"; please unify the spelling.","section":"Fig. 1 caption and reference [58]"},{"comment":"The abstract's claim that the result is \"the last missing piece\" for N3LL'-accurate self-normalized thrust distributions is carefully qualified in Sec. 1, which notes that unnormalized cross sections also require the N3LL' hard-matching coefficient. Consider adding a similar qualifier in the abstract itself.","section":"Abstract and Sec. 1"}],"recommendation":"major_revision","confidential_remarks":"My recommendation is driven by the disclosure gap, not by any demonstrated error. The cross-checks reported in the paper are strong, and the numerical checks of the master integrals substantially mitigate the risk in MI3L_t. I would expect that adding an ancillary file with B_bare^3 and a fuller derivation of MI3L_t can resolve the major comments without new physics. If the authors prefer to keep the expression 'available upon request' only, I would not consider the manuscript suitable for publication in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers the first three-loop result for the boosted heavy-quark jet function, and it is a genuinely serious calculation. The new terms in Eq. (4.7) are the CF CA^2, CF CA TF nl, and CF^2 TF nl pieces of b30; the nl^2 term was already known from renormalon calculus, but the full color structure is new. What the paper does well is cross-check itself. Gauge-parameter independence, non-Abelian exponentiation in d dimensions, reproduction of the known cusp and non-cusp anomalous dimensions, and agreement with the nl^2 prediction are all non-trivial checks, and they pass. The master integrals are given analytically in Appendix F and numerically checked with FIESTA/pySecDec. That is real evidence, not decoration.\n\nThe soft spot is the one the stress-test identifies: the non-planar MI3L_t in Eq. (F.8u) is presented as a closed 3F2 result, but Appendix G only sketches the Mellin-Barnes procedure, and the unpublished B_bare^3 combination in footnote 15 is exactly what inserts that integral into Eq. (4.7). A wrong prefactor or a missed epsilon term there would change the largest color piece of b30. This is a verification gap, not a demonstrated error. The consistency checks pass on the whole combination, but they do not independently re-verify MI3L_t. The quasi-finite integrals expanded to one order less than the claimed need are a minor related concern, not a fatal one.\n\nI do not think the reader was too harsh. A conditional verdict is fair: the physics is credible, but the reproducibility is below what a three-loop fixed-order result should ship with. The four-loop estimate in Sec. 5 is explicitly model-dependent and is honestly labeled as such; that is not a flaw.\n\nWho this is for: anybody working on top-quark mass extraction, boosted-top event shapes, or bHQET/SCET factorization. The paper earns a serious referee, and if the authors provide the full bare B^3 in d dimensions and expand the MI3L_t derivation, it is publishable as is. I would accept it for peer review, conditional on those documentation improvements.","headline":"A well-executed three-loop computation of the bHQET jet function; the main soft spot is documentation of the non-planar master integral, not the physics.","tokens_in":48647,"tokens_out":1110,"would_cite":true,"duration_ms":15150,"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":"The three-loop inclusive jet function for boosted heavy quarks in bHQET is computed, completing the last missing fixed-order ingredient for N$^3$LL$'$ resummed thrust in boosted top-pair events.","keywords":["boosted heavy-quark effective theory","jet function","three-loop QCD","top quark mass","thrust","non-Abelian exponentiation","renormalon","N3LL' resummation"],"falsifier":"Recompute the 20 master integrals of Appendix F to the required order in epsilon with an independent numerical or analytic method; any discrepancy with the quoted expansions would disprove Eq. (4.7).","tokens_in":47571,"feed_emoji":"⚛️","tokens_out":8066,"duration_ms":75597,"temperature":0.7,"pith_summary":"The paper computes the inclusive jet function for boosted heavy quarks to three loops in the strong coupling, within boosted Heavy-Quark Effective Theory. This function controls collimated radiation inside an energetic heavy-quark jet in the regime where the jet invariant mass satisfies $M^2 - m^2 \\ll m^2$, the situation relevant for boosted top-pair production. The central result is the three-loop non-logarithmic coefficient of the logarithm of the position-space jet function, Eq. (4.7), which numerically reads $50.054\\, n_\\ell^2 - 1899.8\\, n_\\ell + 12834$ for $N_c=3$. This completes the last missing fixed-order ingredient for N$^3$LL$'$ resummed self-normalized thrust distributions in the peak region of boosted top-pair events, and with it a more reliable calibration of the top quark mass parameter in parton-shower Monte Carlo generators.","feed_headline":"Three-loop top-jet function is computed","feed_subtitle":"Last missing piece for N3LL' resummed peak-region thrust in boosted top-pair events.","key_machinery":"The central object is the logarithm of the position-space jet function, $\\tilde b(x,\\mu)=\\log[m\\tilde B(x,\\mu)]$, whose non-Abelian exponentiation means only fully connected color factors appear, so the $C_F^3$ and $C_F^2 C_A$ terms vanish at three loops. The argument is carried by the analytic evaluation of the 20 three-loop master integrals, obtained from roughly 1100 Feynman diagrams via automated integral-family assignment, integration-by-parts reduction, and a combination of Feynman-parameter, Mellin-Barnes, integration-by-parts, and quasi-finite-integral techniques. These master integrals enter the renormalized result through the exponent $\\tilde b_{30}$ of Eq. (4.7).","core_discovery":"The paper's main claim is that the three-loop coefficient $\\tilde b_{30}$ of the logarithm of the position-space bHQET jet function is given by the analytic expression in Eq. (4.7), with numerical value $50.054\\, n_\\ell^2 - 1899.8\\, n_\\ell + 12834$ for $N_c=3$. The result passes several internal consistency checks: it satisfies non-Abelian exponentiation in $d$ dimensions, it reproduces the known cusp and non-cusp anomalous dimensions of the jet function through $\\mathcal{O}(\\alpha_s^3)$, and its $n_\\ell^2$ piece agrees with the large-$\\beta_0$ prediction from renormalon calculus. The exact three-loop value lies inside the uncertainty band of the earlier renormalon-based estimate. As by-products, the paper derives the relation between the pole mass and two renormalon-free short-distance jet-mass schemes at $\\mathcal{O}(\\alpha_s^3)$, and estimates the non-logarithmic four-loop coefficient of the jet function from renormalon dominance.","pith_inferences":["If the master integrals hold, the same quasi-finite and Mellin-Barnes machinery can likely be extended to extract the four-loop non-cusp anomalous dimension of the bHQET jet function, which is currently the bottleneck for N$^4$LL.","The agreement between the exact three-loop result and the renormalon-dominance estimate strengthens the case that renormalon calculus can be used to estimate unknown higher-order coefficients of other factorization functions, not just the jet function.","The scheme comparison suggests that for Monte Carlo top-mass calibrations with limited perturbative order, the non-derivative jet-mass may be a more efficient renormalon-free mass definition than the derivative jet-mass.","The same jet function is a building block for hemisphere mass, heavy jet mass and C-parameter distributions; this calculation puts those observables within reach of N$^3$LL$'$ as well once their soft functions and matching coefficients are available at that order."],"forward_implications":["The computation reproduces the three-loop non-cusp anomalous dimension of the bHQET jet function, providing its first direct derivation rather than an inference from RG consistency.","With the hard and soft functions already available, the missing jet function means N$^3$LL$'$ self-normalized 2-jettiness (thrust) distributions in the peak region of boosted top-pair events can now be constructed.","The pole-mass to jet-mass scheme relation at $\\mathcal{O}(\\alpha_s^3)$ provides two practical short-distance mass schemes; the non-derivative jet-mass is well behaved already at one loop while the derivative jet-mass requires two loops before it approaches the MSR benchmark.","The four-loop non-logarithmic coefficient estimated from renormalon dominance gives a numerical target for future explicit four-loop calculations.","Once the four-loop cusp and non-cusp anomalous dimensions are known, the new jet function enables N$^4$LL resummation for these observables."],"supporting_citations":[{"why":"supplies the previous two-loop jet function and the jet-mass scheme definition that the three-loop calculation extends.","marker":"[55]"},{"why":"provides the large-beta0 prediction for the three-loop non-log coefficient and the renormalon-dominance estimate, against which Eq. (4.7) is checked.","marker":"[65]"},{"why":"derives the bHQET peak-region factorization for boosted top jets, the physical context that the jet function enters.","marker":"[18]"},{"why":"gives the three-loop soft-function calculation with the same Wilson-line structures, whose integral-reduction strategy is adapted here.","marker":"[67]"},{"why":"introduces the quasi-finite integral and dimensional-recurrence method used to evaluate the three-loop master integrals.","marker":"[79]"},{"why":"provides the hard matching factor for boosted tops whose RG consistency determines the three-loop non-cusp anomalous dimension used as a check.","marker":"[50]"},{"why":"contributes the known lower-order cusp anomalous dimension coefficients used to verify the renormalized result.","marker":"[100]"},{"why":"supplies the three-loop cusp anomalous dimension that the calculation reproduces as a consistency check.","marker":"[101]"},{"why":"is used together with [50,67] to derive the known three-loop non-cusp anomalous dimension reproduced here.","marker":"[102]"}],"fun_headline_variants":["Three-loop top-jet function: N3LL' thrust within reach","Last missing piece for N3LL' top-mass calibration: 3-loop jet function","Three-loop bHQET jet function: complete for boosted tops","Precise top mass from N3LL' thrust: 3-loop jet function now complete","Three-loop jet function: final ingredient for N3LL' top mass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the correctness of the analytic epsilon expansions of the 20 three-loop master integrals; a mistake in any one of them would change Eq. (4.7).","fun_headline_variants_meta":{"raw":{"variants":["Three-loop top-jet function: N3LL' thrust within reach","Last missing piece for N3LL' top-mass calibration: 3-loop jet function","Three-loop bHQET jet function: complete for boosted tops","Precise top mass from N3LL' thrust: 3-loop jet function now complete","Three-loop jet function: final ingredient for N3LL' top mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00043,"raw_usage":{"total_tokens":2262,"prompt_tokens":1078,"completion_tokens":1184,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":1082}},"tokens_in":694,"tokens_out":1184,"duration_ms":9271,"temperature":1.0,"reasoning_tokens":1082,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:19:31.029570+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the 20 master integrals of Appendix F to the required order in epsilon with an independent numerical or analytic method; any discrepancy with the quoted expansions would disprove Eq. (4.7).","supporting_citations":[],"review_version":1}