{"id":"9f181153-eaa2-4b18-b4a2-1d6cf5289119","arxiv_id":"2411.11080","paper_version":5,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New four-loop coefficient functions are presented for the HQET heavy-light current correlator with light-quark mass terms up to second order.","lead":"The paper calculates a very precise quantum chromodynamics correction to the behavior of heavy-light mesons like B mesons, at the four-loop level of perturbation theory. The results supply new perturbative inputs for sum-rule and lattice calculations of meson properties.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The new 4-loop constants c(mn)_3, r(mn)_3, and c(sum m_i^2)_1 rest entirely on the external four-loop master integrals of Ref. [15] and on LiteRed IBP reduction; the paper supplies no independent check of the finite parts, so any error there shifts every new coefficient.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing point: the four-loop result depends on master integrals from [15] and LiteRed reduction, with no independent numerical check. My stress-test confirms this is the right soft spot because the new coefficients c3 are finite parts of four-loop integrals, and the internal checks (gauge invariance and pole cancellation) constrain only lower-order or pole parts, not the finite pieces that define c3. The lower-order agreement with [7] is necessary but not sufficient. Therefore the central claim is conditional on the external computation being correct. The abstract/full-text mismatch about the large-beta0 and naive nonabelianization analysis is a real completeness issue, but it does not bear on the correctness of the four-loop coefficients; it can be settled by inspecting the ancillary file or by removing the claim. Because the reader already assigned CONDITIONAL for this reason, my assessment leaves the verdict unchanged. No ad hominem is intended; this is a standard call for independent verification in a calculation with no machine-checked proof and no numerical cross-check.","tokens_in":6073,"tokens_out":3797,"duration_ms":98104,"concrete_test":"Independently recompute the master integrals needed for one color component of c(1)_3, for example the CF^3 term in Eq. (5), using a different method such as AMFlow or sector decomposition, and compare the finite parts with the epsilon expansions quoted from Ref. [15] at the precision needed for c3. In parallel, or as a fallback, re-run the IBP reduction of the four-loop diagrams with an independent reduction tool such as FIRE or Kira. If either the master integrals or the reduced linear combinations differ at finite order, the central coefficients are wrong; if they agree, the load-bearing assumption is verified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that c(mn)_3 in Eqs. (5)-(7), r(mn)_3 in Eqs. (12)-(14), and c(sum m_i^2)_1 in Eq. (9) are correct. The finite parts of these coefficients are obtained by IBP reduction (LiteRed v2) of four-loop diagrams to master integrals whose epsilon expansions were taken from Lee and Pikelner [15]. No independent numerical or analytic verification of these finite parts is presented. The checks cited in the paper—gauge-parameter cancellation up to xi^1 at four loops and cancellation of 1/epsilon poles in the renormalized Pi(tau) and rho(omega)—do not fix the finite parts: pole terms are determined by lower-loop coefficients and anomalous dimensions, whereas c3 depends on the finite parts of the new four-loop integrals. Agreement of c1 and c2 with [7] likewise does not constrain c3, since c3 receives contributions from integrals absent at lower orders. Thus the entire new result inherits the correctness of every master integral in [15] and of every IBP identity used by LiteRed, and that inheritance is currently an untested assumption. A secondary completeness issue is that the abstract announces a large-beta0 all-orders analysis and a statement about naive nonabelianization, but the supplied full text contains no such analysis; if that material lives only in the attached file, it needs to be inspectable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the four-loop perturbative contribution to the correlator of two HQET heavy-light currents, with the light-quark mass expansion truncated at quadratic order. It reports new four-loop constants c_3 for the unit, m, and m^2 Wilson coefficients in Eqs. (5)-(7), the corresponding spectral-density constants r_3 in Eqs. (12)-(14), and the new three-loop coefficient c(sum m_i^2)_1 in Eq. (9). The calculation uses qgraf, FORM, LiteRed v2, and the four-loop master integrals of Ref. [15], and it verifies that the one- and two-loop coefficients agree with Ref. [7]. The abstract additionally advertises a leading large-beta0 all-orders analysis and a statement about naive nonabelianization, but that analysis is not present in the supplied text.","tokens_in":6319,"tokens_out":7357,"duration_ms":77439,"significance":"If the new finite parts are correct, this is a useful four-loop advance for a correlator used in QCD sum rules and lattice comparisons of heavy-light mesons; in particular, the m^2 and sum m_i^2 terms are relevant for SU(3)-flavor-breaking quantities such as B_s^* versus B^*. The fixed-order part is internally coherent: it has no fitted parameters, it reproduces the published two- and three-loop coefficients, and it passes pole-cancellation and partial gauge-invariance checks. The main value of the paper is the new finite constants, and precisely those are the least verified part of the calculation. The lower-order inputs are external and independent, so there is no circularity in the RG extraction of the logarithmic structure. The abstract's all-orders claims, if substantiated elsewhere, would raise significance, but they cannot be evaluated from the supplied text.","major_comments":[{"comment":"The new four-loop coefficients c(mn)_3 in Eqs. (5)-(7), r(mn)_3 in Eqs. (12)-(14), and c(sum m_i^2)_1 in Eq. (9) are finite parts that rest on the epsilon expansions of the four-loop master integrals of Ref. [15] and on the completeness of the LiteRed v2 IBP reduction. The manuscript provides no independent numerical or analytical check of these finite parts. The checks that are reported are not sensitive to the new finite terms: cancellation of 1/epsilon poles is fixed by lower-loop coefficients and anomalous dimensions, the gauge-parameter check is performed only to order xi^1 (see next comment), and agreement of c1 and c2 with Ref. [7] is insensitive to the new four-loop integrals. I ask for at least one independent cross-check (for example, numerical evaluation of the most complicated master integrals, a second reduction code, or a subset of diagrams evaluated by a different method) or a clear statement that the result is contingent on the external input without further verification.","section":"Eqs. (5)-(7), (12)-(14)"},{"comment":"The four-loop gauge-invariance check is reported only to first order in the covariant-gauge parameter xi. Because the xi dependence of individual diagrams is polynomial in the number of gluon propagators, a partial check through O(xi) does not exclude O(xi^2) and higher contributions; a truly gauge-invariant result would require those to cancel as well. If the higher-order xi terms are zero for structural reasons, that argument should be stated; otherwise the check should be described as partial rather than 'strong'.","section":"Section 1, calculation setup"},{"comment":"The abstract claims that the leading large-beta0 limit is considered, that Borel images contain renormalon poles, and that naive nonabelianization works surprisingly poorly, but none of these results appears in the supplied full text. The only referenced attachment is said to contain momentum-space Wilson coefficients, not a large-beta0 analysis. This mismatch makes a headline result uninspectable. The all-orders analysis should be included in the submission, or the abstract and title should be revised to describe only the fixed-order four-loop result.","section":"Abstract vs. full text"}],"minor_comments":[{"comment":"There are several typographical errors: 'Cancelation' should be 'Cancellation', 'simler diaghrams' should be 'simpler diagrams', 'up to xi^1' has a missing space, and 'recently is has been calculated' should be 'it has recently been calculated'.","section":"Throughout"},{"comment":"The title of the version under review says 'and beyond' while the body and conclusion only claim a four-loop calculation; align the title with the content unless the large-beta0 section is added.","section":"Title and abstract"},{"comment":"The notation in Eq. (3) mixes superscripts and subscripts for the coefficients c(mn)_i; a short table or a compact notation for the three operators would improve readability.","section":"Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"Please verify that the submitted file set contains the attachment referenced in the text and any file implementing the large-beta0 analysis advertised in the abstract; as supplied, the full text does not allow those claims to be checked."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is the fixed-order four-loop calculation: the new coefficients c(mn)_3, r(mn)_3, and c(sum m_i^2)_1 in Eqs. (5)-(7), (9), (12)-(14). These are not reducible to the earlier three-loop results, and the paper reproduces the two- and three-loop limits, passes gauge-parameter cancellation at four loops, and shows 1/epsilon pole cancellation. That is a solid, useful result for the HQET sum-rule and lattice-comparison community, and the author says plainly which parts are new and which come from [7]. Credit where due: the RG structure is laid out carefully, and the note about the typo in [7] is helpful.\n\nThe soft spots, in proportion: first, the four-loop finite parts inherit the correctness of the Lee-Pikelner master integrals and the LiteRed IBP reduction, and the paper gives no independent numerical or analytic check of those finite parts. That is a common and often unavoidable feature of four-loop calculations, not a fatal flaw, but it is worth stating explicitly. Second, and more important: the abstract announces a leading large-beta0 all-orders analysis, renormalon poles, and a statement about naive nonabelianization. None of that appears in the supplied full text. If it lives in the attached coefficient file, the file needs to be inspectable; if it does not exist, the abstract overclaims and should be trimmed. This is a legitimate completeness issue, not a manufactured one.\n\nI disagree with any suggestion that the calculation is circular or that the self-citations are a problem. The input anomalous dimensions and master integrals are external published results, and the new constants are fixed by the bare calculation plus standard RG equations. The stress-test concern about unverified finite parts is correct but not disqualifying; it is the usual state of the art at this loop order.\n\nBottom line: this paper deserves a serious referee. The fixed-order result is likely correct and citable, and the abstract mismatch is fixable by either supplying the large-beta0 material or removing it. I would send it to review with a request to make the attachment available and clarify the abstract's scope.","headline":"A credible four-loop HQET correlator calculation with genuinely new constants, but the abstract promises a large-beta0 analysis that the supplied text does not contain.","tokens_in":6900,"tokens_out":1073,"would_cite":true,"duration_ms":13101,"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":"Four-loop coefficients for the HQET heavy-light current correlator are computed, extending the perturbative expansion to order α_s^3 and confirming the renormalization-group structure.","keywords":["HQET","heavy-light quark currents","correlator","operator product expansion","four loops","light-quark masses","QCD sum rules","renormalons"],"falsifier":"Recompute the four-loop master integrals with an independent numerical method (for example, sector decomposition) at a fixed value of the space-time dimension $\\varepsilon$, insert them into the reduced diagrams, and compare the resulting coefficient $c^{(1)}_3$ with Eq. (5); a mismatch beyond the integration errors would disprove the new result.","tokens_in":5823,"feed_emoji":"⚛️","tokens_out":6934,"duration_ms":61464,"temperature":0.7,"pith_summary":"Four-loop (order $α_s^{3}$) Wilson coefficients in the operator product expansion of the correlator of two HQET heavy-light currents are computed, together with the corresponding spectral-density coefficients. The new results are the constants $c^{(1)}_3$, $c^{(m)}_3$, $c^{(m^2)}_3$ and the two-loop coefficient $c^{(\\sum m_i^2)}_1$, given as closed analytic expressions in Riemann zeta values and powers of π. They extend the earlier three-loop calculation and agree with lower-loop results, while passing gauge-invariance and pole-cancellation checks. These coefficients are the input needed to push QCD sum-rule and lattice comparisons for heavy-light mesons to next order, and the mass-dependent terms describe SU(3) flavour breaking between $B_s$ and $B$ mesons.","feed_headline":"Four-loop HQET heavy-light current correlator coefficients calculated","feed_subtitle":"New α_s^3 terms refine QCD sum rules for heavy-light mesons and lattice comparisons.","key_machinery":"The renormalization-group exponentiation formula (3), which expresses each coefficient function as an exponentiated series in $\\alpha_s/(4\\pi)$ with logarithmic terms $L_\\tau$ fixed by the anomalous dimensions $\\gamma_n$ and the $\\beta$ function, is the machinery that turns loop computation into finite constants. The divergent parts (poles in $1/\\varepsilon$) of the four-loop diagrams fix the logarithmic terms through the anomalous dimensions and $\\beta$, and the constants $c_3$ are then read off from the finite remainder after subtracting these logarithms. The actual integrals are handled by IBP reduction with LiteRed v2 and the $\\varepsilon$-expanded master integrals of Lee and Pikelner.","core_discovery":"The central claim is that the renormalization-group structure of the Euclidean correlator, together with four-loop master integrals, fixes the finite parts of the Wilson coefficients at order $\\alpha_s^3$. Explicitly, Eqs. (5)–(7) give the four-loop constants $c^{(1)}_3$, $c^{(m)}_3$, $c^{(m^2)}_3$, Eq. (9) gives the two-loop constant $c^{(\\sum m_i^2)}_1$, and Eqs. (12)–(14) give the spectral-density analogues $r^{(mn)}_3$. The coefficients at two and three loops reproduce the known results of the earlier three-loop OPE calculation, and the four-loop results are new. The calculation is performed in a general covariant gauge; the final expressions are independent of the gauge parameter, and the renormalized correlator and spectral density are free of $1/\\varepsilon$ poles, both serving as internal consistency checks. The abstract also states that in the leading large-$\\beta_0$ limit the highest-$n_f$ terms are summed to all orders and that the Borel image exhibits renormalon poles, with naive nonabelianization failing to approximate the results well.","pith_inferences":["If the four-loop constants are numerically large compared with lower orders, sum-rule extractions may need to account for a slowly converging series; the paper does not assess this.","The same RG-plus-master-integral framework can be extended to the dimension-3 operators (light-quark condensate and cubic mass combinations), which the paper leaves out because of operator mixing.","The poor performance of naive nonabelianization reported in the abstract suggests that resummation of these correlators at higher orders should not rely on the standard large-$n_f$ approximation without testing other schemes."],"forward_implications":["The $\\alpha_s^3$ terms provide the next order in the OPE used by QCD sum rules for heavy-light meson masses and decay constants.","The light-quark-mass terms quantify SU(3) flavour breaking, so they sharpen predictions for the $B_s$–$B$ mass splitting.","The spectral-density coefficients $r^{(mn)}_3$ allow a direct order-by-order comparison with lattice QCD results for the same correlator.","The claimed all-order large-$\\beta_0$ analysis yields renormalon poles in the Borel image, which can be used to assess the asymptotic behaviour of the perturbative series."],"supporting_citations":[{"why":"Supplies the $\\varepsilon$-expanded four-loop master integrals that all new coefficients are built from.","marker":"[15]"},{"why":"Provides the LiteRed IBP reduction tool used to reduce the four-loop diagrams to master integrals.","marker":"[13]"},{"why":"Supplies the LiteRed version 2 upgrade used in the present calculation.","marker":"[14]"},{"why":"Gives the three-loop anomalous dimension of the heavy-light current, which fixes the logarithmic structure in the RG formulas.","marker":"[16]"},{"why":"Provides the earlier three-loop OPE calculation whose two- and three-loop coefficients are reproduced and whose notation and renormalization-group structure are adopted.","marker":"[7]"}],"fun_headline_variants":["Four-loop HQET heavy-light current correlator computed","New α_s^3 Wilson coefficients in HQET heavy-light","Renormalon poles found in HQET heavy-light correlator","Naive nonabelianization fails: HQET renormalons at all orders"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The four-loop constants are only as reliable as the four-loop master integrals from Lee-Pikelner and the LiteRed IBP reduction, and the known three-loop anomalous dimension of the heavy-light current; the paper provides no independent numerical check of the new four-loop results.","fun_headline_variants_meta":{"raw":{"variants":["Four-loop HQET heavy-light current correlator computed","New α_s^3 Wilson coefficients in HQET heavy-light","Renormalon poles found in HQET heavy-light correlator","Naive nonabelianization fails: HQET renormalons at all orders"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000275,"raw_usage":{"total_tokens":1607,"prompt_tokens":874,"completion_tokens":733,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":490,"completion_tokens_details":{"reasoning_tokens":660}},"tokens_in":490,"tokens_out":733,"duration_ms":121813,"temperature":1.0,"reasoning_tokens":660,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:55:19.974547+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the four-loop master integrals with an independent numerical method (for example, sector decomposition) at a fixed value of the space-time dimension $\\varepsilon$, insert them into the reduced diagrams, and compare the resulting coefficient $c^{(1)}_3$ with Eq. (5); a mismatch beyond the integration errors would disprove the new result.","supporting_citations":[{"cited_title":"Four-loop HQET propagators from the DRA method","cited_arxiv_id":"2211.03668","evidence_quote":"Supplies the $\\varepsilon$-expanded four-loop master integrals that all new coefficients are built from."}],"review_version":1}