{"id":"b5a020ce-9a39-48a0-ab3a-e51219afedce","arxiv_id":"2411.15402","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A characteristic-operator reformulation of the PMC scale-setting method gives a scale-invariant N4LO prediction Gamma(H->bbar)=2.3819 MeV with Q*=55.2916 GeV.","lead":"Physicists propose a new mathematical operator to make QCD predictions independent of the arbitrary renormalization scale, and use it to compute the Higgs decay width into bottom quarks at four-loop accuracy. The method yields a scale-invariant width of 2.3819 MeV, but the improvement over the standard result is small and the choice of which terms to absorb is debated.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim depends on the contested choice to treat nf terms in gamma_SS as conformal coefficients; the paper provides no numerical check of the alternative beta-expansion treatment.","rationale":"The paper is a serious and internally consistent application of the PMCs framework to H→bbar. The CO formalism is a compact way to organise the scale-displacement and degeneracy relations, and the numerical results reproduce the conventional N4LO width (Eq. 4.4) when standard inputs are used, indicating no gross algebraic error. The claimed outputs that the abstract highlights, Q*=55.2916 GeV and Gamma=2.3819 MeV, are direct consequences of the procedure for deciding which nf terms are non-conformal and should be absorbed into the running coupling and mass. That decision is not a theorem but a convention: the paper's own Appendix D acknowledges a dispute with Ref. [90] and argues its position by analogy with the QED Adler function, an analogy the opposing literature does not accept. Because the final series and scale are sensitive to this convention, the central claim is conditional on it. The massless-quark approximation, by contrast, only affects the numerical value at the permille level and is secondary. I found no evidence that the CO derivation itself contains an algebraic error; the proposed test therefore probes the one convention on which the result hinges. If the alternative β-expansion treatment yields a numerically similar Q* and width, the concern is lifted; if not, the quoted central values would need revision. The reader's CONDITIONAL verdict is therefore appropriate, and no change is needed.","tokens_in":24719,"tokens_out":16324,"duration_ms":137277,"concrete_test":"Perform an independent PMC single-scale calculation for the same N4LO H→bbar series in which the nf terms of gamma_SS are first converted into {beta_i}-terms according to the β-expansion recipe of Ref. [90], then re-solve Eq. (2.24) for Q* and recompute Gamma using Eq. (3.10). Compare the new Q* and Gamma with Eqs. (4.1) and (4.13): if Q* shifts by more than ~1 GeV or Gamma shifts by more than the ±0.021 MeV alpha_s error, the paper's central numerical claim is not robust to the disputed treatment. As a secondary check, include the finite bottom-mass O(alpha_s^3) corrections of Ref. [59] to quantify the neglected massless-quark approximation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the CO/PMCs procedure yields a scheme-and-scale invariant series with Q*=55.2916 GeV and Gamma(H→bbar)|_PMC = 2.3819 MeV—rests on the convention, stated in Sec. 3 and defended in App. D, that the nf terms in the anomalous dimension gamma_SS (App. C) are 'not RGE-involved' and must be kept as conformal coefficients r_{i,0}, rather than re-expressed as {beta_i}-terms. This choice is load-bearing because Q* is fixed by Eq. (2.24) using only the j>=1 non-conformal sector, while all r_{i,0} are evaluated at Q* in Eq. (3.10). If the alternative β-expansion treatment advocated in Ref. [90] is correct, the degeneracy relations (2.17)-(2.21) and the resulting Q* and width would all shift. The paper's App. D argues by analogy with the QED photon anomalous dimension in D_ns, but this analogy is contested (Refs. [89]-[92]), and no numerical estimate of the shift under the alternative treatment is provided. Thus the uniqueness and correctness of the quoted Q* and Gamma remain an unsettled convention, not a derived fact. The massless-quark approximation is a secondary permille-level numerical issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a 'characteristic operator' Dhat_{n_gamma,n_beta} (Eq. (2.6)) that packages the beta-function and the quark-mass anomalous dimension into a single differential operator, allowing a compact derivation of scale-displacement relations, QCD degeneracy relations, and PMC single-scale formulas that run alpha_s and the MS mass simultaneously. The formalism is applied to Gamma(H -> b bbar) at N4LO in the massless approximation: the PMC scale is fixed to Q* = 55.2916 GeV (Eq. (4.1)) and the scale-invariant width is Gamma(H -> b bbar)|_PMC = 2.3819^{+0.0230}_{-0.0231} MeV (Eq. (4.13)), with uncertainties from Delta alpha_s(M_Z), Delta M_H, Delta m_b, and a Bayesian estimate of the N5LO contribution. All operator coefficients and input series are listed in Appendices A-C.","tokens_in":1,"tokens_out":10811,"duration_ms":168564,"significance":"The CO formalism is a useful technical generalization of PMCs: it unifies the treatment of running coupling and running mass, simplifies the degeneracy relations, and is presented with enough coefficient-level detail that the derivation can be checked step by step. The numerical outcome is transparent and the claimed scale independence is explicitly demonstrated (Fig. 2, Eq. (4.13)); the result is consistent with the conventional N4LO prediction within the quoted uncertainties. The framework would be a valuable benchmark if the two load-bearing caveats identified below are resolved; as it stands, the central value and its precision statement depend on a contested convention for how n_f terms in the anomalous dimension gamma_SS are classified. The paper deserves credit for listing all input coefficients and for giving a clear Bayesian procedure for the uncalculated N5LO term.","major_comments":[{"comment":"The central numerical result depends on the decision, stated in Sec. 3 and defended in App. D, that the n_f terms in the anomalous dimension gamma_SS are 'not RGE-involved' and are therefore retained as conformal coefficients r_{i,0} (Eq. (3.9)) instead of being converted into {beta_i}-terms before solving Eq. (2.24) for Q*. This choice is load-bearing: if one instead follows the beta-expansion treatment advocated in Ref. [90] and related literature, the degeneracy relations (2.17)-(2.21), the resulting Q*, and therefore Eq. (4.13) would all shift. App. D argues by analogy with the QED photon anomalous dimension in the Adler function, but that analogy is precisely the point disputed in Refs. [89]-[92]. Please provide a quantitative sensitivity check: repeat the PMC scale-setting with the alternative treatment of the gamma_SS n_f terms and report the shift in Q* and Gamma(H -> b bbar), or give a direct RGI-based proof specific to gamma_SS showing why these terms cannot enter Eq. (2.24). Without this, the uniqueness of Q* = 55.2916 GeV as 'the' PMC scale is not established.","section":"Sec. 3, App. D, Eqs. (2.24), (3.9), (3.10)"},{"comment":"The analysis is restricted to massless quarks in the QCD corrections (Sec. 1), yet Eq. (4.13) is quoted as the total decay width Gamma(H -> b bbar) with a permille-level central value and a detailed error budget. The known O(alpha_s^3) finite-bottom-mass corrections of Ref. [59] are not included and their numerical size is not estimated. These corrections could be comparable to the Bayesian UHO error Delta Gamma = +/- 0.0001 MeV quoted in Sec. 4.2, and they are certainly relevant for a four-significant-digit central value. Please either include the mass effects, provide a numerical estimate of their contribution to the error budget, or state explicitly throughout the abstract and Sec. 4 that the quoted width applies only to the massless-quark approximation.","section":"Sec. 1, Sec. 4, Ref. [59]"},{"comment":"The paper states in the abstract and Sec. 5 that the CO framework 'achieves a scheme-and-scale invariant pQCD series.' The numerical evidence in Sec. 4 concerns scale independence only: Eq. (3.10) is constructed to be independent of the initial renormalization scale, and Fig. 2 displays this explicitly. Scheme independence is not demonstrated for the H -> b bbar series; the coefficients r_{i,0} are MS-scheme quantities and the CO is built from the MS-scheme beta and gamma_m. Please state explicitly which previous general PMC result (e.g., Ref. [22]) guarantees scheme invariance in the present case, or provide a numerical scheme-variation check. As written, the scheme-invariance claim goes beyond what this paper demonstrates.","section":"Abstract, Sec. 5"}],"minor_comments":[{"comment":"Eq. (2.12) uses binomial coefficients C_i^k, while C_k^j is defined just below with reversed indices; the notation should be made consistent.","section":"Eq. (2.12)"},{"comment":"Eq. (4.8): the two cases of the Bayesian coefficient are garbled in the typeset formula; the conditions should read DoB <= p/(p+1) and DoB >= p/(p+1).","section":"Eq. (4.8)"},{"comment":"The appendix heading reads 'Appdendix'; this should be corrected.","section":"Appendix heading"},{"comment":"Sec. 4.2 contains 'a truncated perturbation series cannot not automatically satisfy'; the double negative should be removed.","section":"Sec. 4.2"},{"comment":"The abstract says 'errors are squared averages'; the quoted uncertainties are combined in quadrature, so the phrase should be 'quadrature sum' or 'root-sum-square', not 'squared averages'.","section":"Abstract"},{"comment":"Figure 3 has a dense legend; please spell out in the caption which symbols correspond to which scale choice and order, and ensure that the text's description of 'red solid squares' matches the figure.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the paper is technically competent and the derivation is reproducible from the appendices. The main risk is the contested treatment of gamma_SS; I believe this is fixable in revision by adding a quantitative sensitivity analysis rather than by rewriting the paper. The paper is within JHEP's scope as a methods-oriented phenomenological study, and the self-citation pattern, while noticeable, is conventional for this research program."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things before you read it. The characteristic operator is a genuinely useful formal device: it gives a compact derivation of the PMC single-scale formulas when alpha_s and an MS-scheme mass run simultaneously, and it extends cleanly to other RGE-running quantities. That part is solid. The numerical H->bbar example, though, is not a big departure--2.3819 MeV sits 0.1% from the conventional 2.3842 MeV, within uncertainties--and the 'scheme-and-scale invariant' claim rests on a convention that the authors themselves flag but do not defend quantitatively.\n\nCredit where due. The derivation in Sec. 2 is coherent; the degeneracy relations reduce to earlier PMCs forms in the right limit; Appendix C supplies all input coefficients and the quoted quadrature errors match the stated inputs. I checked the key relations and found no internal numerical inconsistency. The paper is also honest about the massless-quark approximation and about the contested treatment of nf terms in gamma_SS. That is more than many papers in this area do.\n\nSoft spots. The load-bearing choice is the treatment of the nf terms in gamma_SS. The authors keep them as conformal coefficients, arguing by analogy with the QED photon anomalous dimension in the Adler function. That analogy is exactly what Kataev and collaborators dispute, and it matters: Q* is fixed by Eq. (2.24) using only the j>=1 non-conformal sector, so if those gamma_SS nf terms should be beta-expanded, Q* and the extracted width shift. App. D gives a conceptual argument but no numerical estimate of the alternative. This is the one thing I would want resolved before trusting the central claim. The second issue is minor: the O(alpha_s^3) finite-mb effect from Ref. [59] is neglected with a massless-quark argument and no permille-level estimate. Third, the practical gain is modest because the N4LO conventional series already has small scale dependence and the dominant uncertainty is external alpha_s. That is not a flaw, but it means the paper's value is mainly methodological.\n\nBottom line: this is a serious, checkable paper for people who work on renormalization-scale setting. It deserves a referee, not a desk rejection. Ask the authors to quantify the alternative gamma_SS treatment and the massless-quark error; then the conditional verdict can be resolved.","headline":"A solid, internally consistent method paper that extends PMC single-scale setting to running masses via a characteristic operator; the H->bbar number is not numerically new, and the scheme-invariance claim depends on a contested convention about gamma_SS nf terms that the paper does not quantify.","tokens_in":25525,"tokens_out":3367,"would_cite":true,"duration_ms":31663,"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":"A characteristic operator makes the four-loop Higgs-to-bottom decay width independent of the renormalization scale.","keywords":["characteristic operator","principle of maximum conformality","renormalization-scale ambiguity","Higgs decay to bottom quarks","running quark mass","N4LO QCD corrections","Bayesian uncertainty estimate"],"falsifier":"Recompute the N$^4$LO width with the alternative $\\beta$-expansion prescription for the $\\gamma_{SS}$ $n_f$ terms (the approach of Ref. [90]); if $Q_*$ and $\\Gamma(H\\to b\\bar b)$ move outside the quoted $\\pm0.023$ MeV error band, the non-absorption rule is the decisive assumption. Independently, adding the O($\\alpha_s^3$) finite bottom-mass effect of Ref. [59] would test the massless-quark approximation at the permille level.","tokens_in":24472,"feed_emoji":"⚛️","tokens_out":9520,"duration_ms":82430,"temperature":0.7,"pith_summary":"The paper tries to prove that a single mathematical operator, the characteristic operator, lets the principle of maximum conformality fix the strong coupling and the running bottom-quark mass at the same time, producing a perturbative prediction that is independent of the arbitrary renormalization scale. It applies this to the dominant decay of the Standard-Model Higgs boson, $H \\to b\\bar b$, at fourth-order QCD accuracy. The result is a single effective scale $Q_* = 55.2916$ GeV and a decay width $\\Gamma(H\\to b\\bar b) = 2.3819^{+0.0230}_{-0.0231}$ MeV with the scale ambiguity removed. If true, this would make the theory prediction genuinely scheme-and-scale invariant and give a cleaner foundation for estimating missing higher-order terms.","feed_headline":"One operator fixes Higgs-to-bottom width at 2.3819 MeV","feed_subtitle":"Renormalization-scale ambiguity vanishes at four loops, with the effective QCD scale fixed to 55.2916 GeV","key_machinery":"The characteristic operator is defined as $\\hat D_{n_\\gamma,n_\\beta} = n_\\gamma\\gamma_m + n_\\beta\\,\\beta/\\alpha_s + \\beta\\,\\partial/\\partial\\alpha_s$, where $\\beta(\\alpha_s)$ and $\\gamma_m(\\alpha_s)$ are the renormalization-group functions for the strong coupling and the $\\overline{\\rm MS}$ quark mass. It packages the scale-running of both quantities in a single object, so that repeated application generates the scale-displacement relation and the new QCD degeneracy relations used to split the series into conformal and non-conformal parts. In the $H\\to b\\bar b$ application, the coefficients of the two-point correlator and the anomalous dimension $\\gamma_{SS}$ are inserted into this machinery; demanding that all non-conformal terms vanish fixes $Q_*$ and leaves a conformal series $\\sum r_{i,0}\\alpha_s^{i-1}(Q_*)$ multiplying $\\overline m_b^2(Q_*)$.","core_discovery":"Working in the $\\overline{\\rm MS}$ scheme and treating the QCD corrections as massless, the paper demonstrates that all non-conformal $\\{\\beta_i\\}$-terms in the N$^4$LO series for $\\Gamma(H\\to b\\bar b)$ can be absorbed into a single effective coupling $\\alpha_s(Q_*)$ and a single effective running mass $\\overline m_b(Q_*)$ simultaneously. With the characteristic-operator form of the PMC single-scale equations, the effective scale is fixed to N$^3$LL accuracy, $Q_* = 55.2916$ GeV, and the resulting conformal series is independent of the initial renormalization scale. The central prediction is $\\Gamma(H\\to b\\bar b)|_{\\rm PMC} = 2.3819^{+0.0230}_{-0.0231}$ MeV, with errors obtained by adding in quadrature the effects of $\\Delta\\alpha_s(M_Z)$, $\\Delta M_H$, $\\Delta\\overline m_b(\\overline m_b)$, and a Bayesian estimate of the uncalculated N$^5$LO contribution ($\\pm0.0001$ MeV at 95.5% degree of belief).","pith_inferences":["The operator construction is not limited to $\\alpha_s$ and the quark mass: any renormalization-group-evolving input, such as a parton distribution or fragmentation function with an anomalous dimension of the same form, could in principle be handled the same way, though the paper only demonstrates the quark-mass case.","Applying the same operator machinery to $H\\to gg$, or to the ratio $\\Gamma(H\\to b\\bar b)/\\Gamma(H\\to gg)$, would test whether the simultaneous-running scheme remains scale invariant when the mass factor does not enter the leading order in the same way.","If a future high-precision measurement of $\\Gamma(H\\to b\\bar b)$ approaches the quoted uncertainty, the scale-invariant series could be inverted to extract $\\alpha_s(M_Z)$ more precisely, a direction the paper only hints at in its closing remark."],"forward_implications":["The N$^4$LO total width no longer depends on the initial renormalization scale, with the effective scale fixed at $Q_* = 55.2916$ GeV to N$^3$LL accuracy.","The PMC series converges substantially faster than the conventional series: the higher-order fractional contributions drop from $\\{20.3\\%, 3.7\\%, 0.19\\%, 0.14\\%\\}$ at $\\mu_r = M_H$ to $\\{6.8\\%, 0.67\\%, 0.04\\%, 0.01\\%\\}$.","The total uncertainty is dominated by $\\Delta\\alpha_s(M_Z)$, which contributes about $\\pm0.021$ MeV of the quoted $\\pm0.023$ MeV, so the precision benefit depends mainly on how well the strong coupling is known.","The Bayesian estimate of the uncalculated N$^5$LO term is $\\pm0.0001$ MeV for the PMC series, compared with $+0.0004/-0.0024$ MeV for the conventional series over $\\mu_r \\in [M_H/2, 2M_H]$.","Because the PMC and conventional central values agree within conventional scale errors, the practical consequence is a scale-independent prediction with a defensible error budget rather than a numerically very different width."],"supporting_citations":[{"why":"Supplies the O($\\alpha_s^4$) massless scalar-correlator coefficients that define the fixed-order series for $\\Gamma(H\\to b\\bar b)$.","marker":"[51]"},{"why":"Provides the explicit analytical O($\\alpha_s^4$) correlator results used as the N$^4$LO input.","marker":"[52]"},{"why":"Defines the single-scale PMC approach that the characteristic-operator framework generalizes.","marker":"[20]"},{"why":"Gives the QCD degeneracy relations that the paper rederives in operator form to separate conformal from non-conformal terms.","marker":"[17]"},{"why":"Introduces the renormalization-group treatment combining the $\\beta$-function and the quark-mass anomalous dimension, which is the basis for running $\\alpha_s$ and $m_b$ simultaneously.","marker":"[33]"},{"why":"Supplies the earlier PMC multi-scale analysis of $H\\to b\\bar b$ that this paper improves by also running the $b$-quark mass.","marker":"[70]"},{"why":"Presents the competing $\\beta$-expansion treatment of $n_f$ terms in anomalous dimensions that the paper argues should not be absorbed.","marker":"[90]"},{"why":"Computes the O($\\alpha_s^3$) finite bottom-mass effect that the massless-quark approximation neglects.","marker":"[59]"},{"why":"Provides the world-average input values for $\\alpha_s(M_Z)$, $M_H$, and $m_b$ used in the numerical results.","marker":"[42]"}],"fun_headline_variants":["Scale-free Higgs width: 2.3819 MeV at N4LO","One operator erases renormalization-scale ambiguity in H→bb","Characteristic operator makes Higgs width scale-invariant","Four-loop Higgs width now independent of renormalization scale","New operator fixes QCD scale to 55.2916 GeV for H→bb width"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands or falls on whether the flavour-number terms in the anomalous dimension $\\gamma_{SS}$ should stay out of the effective coupling; the paper keeps them as fixed conformal coefficients, whereas an alternative $\\beta$-expansion treatment would absorb them and shift both $Q_*$ and the width.","fun_headline_variants_meta":{"raw":{"variants":["Scale-free Higgs width: 2.3819 MeV at N4LO","One operator erases renormalization-scale ambiguity in H→bb","Characteristic operator makes Higgs width scale-invariant","Four-loop Higgs width now independent of renormalization scale","New operator fixes QCD scale to 55.2916 GeV for H→bb width"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000744,"raw_usage":{"total_tokens":3840,"prompt_tokens":1237,"completion_tokens":2603,"prompt_tokens_details":{"cached_tokens":1152},"prompt_cache_hit_tokens":1152,"prompt_cache_miss_tokens":85,"completion_tokens_details":{"reasoning_tokens":2525}},"tokens_in":85,"tokens_out":2603,"duration_ms":203571,"temperature":1.0,"reasoning_tokens":2525,"cache_read_input_tokens":1152,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:21:18.706114+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the N$^4$LO width with the alternative $\\beta$-expansion prescription for the $\\gamma_{SS}$ $n_f$ terms (the approach of Ref. [90]); if $Q_*$ and $\\Gamma(H\\to b\\bar b)$ move outside the quoted $\\pm0.023$ MeV error band, the non-absorption rule is the decisive assumption. Independently, adding the O($\\alpha_s^3$) finite bottom-mass effect of Ref. [59] would test the massless-quark approximation at the permille level.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the O($\\alpha_s^4$) massless scalar-correlator coefficients that define the fixed-order series for $\\Gamma(H\\to b\\bar b)$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the explicit analytical O($\\alpha_s^4$) correlator results used as the N$^4$LO input."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the single-scale PMC approach that the characteristic-operator framework generalizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the renormalization-group treatment combining the $\\beta$-function and the quark-mass anomalous dimension, which is the basis for running $\\alpha_s$ and $m_b$ simultaneously."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the earlier PMC multi-scale analysis of $H\\to b\\bar b$ that this paper improves by also running the $b$-quark mass."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the competing $\\beta$-expansion treatment of $n_f$ terms in anomalous dimensions that the paper argues should not be absorbed."},{"cited_title":"Analytic decay width of the Higgs boson to massive bottom quarks at order $\\alpha_s^3$","cited_arxiv_id":"2411.07493","evidence_quote":"Computes the O($\\alpha_s^3$) finite bottom-mass effect that the massless-quark approximation neglects."}],"review_version":1}