{"id":"d3f43ed0-8e39-4490-964a-628adeb2b2f3","arxiv_id":"1909.01370","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"PV-regulated pole mass expansions produce a 28 MeV purely theoretical error for the top quark mass conversion at 163 GeV, with lattice cross-checks for the bottom and charm sectors.","lead":"This paper develops hyperasymptotic expansions for the heavy quark pole mass defined with a principal value prescription, and applies them to charm, bottom, and top quarks. It gives a new estimate of the theoretical uncertainty in converting the MS top mass to this regulated pole mass, about 28 MeV at mbar_t = 163 GeV, before adding the separate strong-coupling uncertainty.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 28 MeV theoretical error is conditional on the u=1 IR renormalon being absent; no quantitative bound is given, and an O(Λ^2/μc) contribution at μc=5 GeV would shift the central value.","rationale":"The paper has genuine supporting evidence: the large-β0 toy model demonstrates the hyperasymptotic subtraction works, and the quenched lattice and B-meson analyses give mutually consistent results with plausible error budgets. The central numerical claim, however, depends on what is subtracted at μc=5 GeV. Eq. (59) explicitly accounts for the leading u=1/2 renormalon and leaves a remainder whose size is set by the next renormalon, and the paper's own Sec. V.C notes that the existence of the u=1 IR renormalon is a matter of debate. No quantitative limit on Z_2 is derived; the f_n signs and the zero result of [38] are suggestive but not a bound. Since an O(Λ^2/μc) contribution at μc=5 GeV is numerically of the same order as the quoted 28 MeV, this is the most load-bearing uncertainty in the central claim. The reader's verdict of CONDITIONAL already captures this, so the stress-test does not change the verdict; it sharpens the specific condition that would need to be verified. The abstract's omission of the α_s error is a presentation issue, secondary to the renormalon assumption but worth fixing in a revision.","tokens_in":34305,"tokens_out":8929,"duration_ms":91753,"concrete_test":"Fit the B-meson spin-averaged mass relation, Eq. (49) with Eq. (50), allowing an explicit mΩ_2 term with free normalization Z_2 for the u=1 IR renormalon, using the same perturbative inputs as the paper; obtain the 68% CL interval for Z_2. Then recompute Eq. (59) at μc=5 GeV with this interval and propagate to m_t,PV in Eq. (66). If the implied shift or added uncertainty in m_t,PV exceeds about 10 MeV, the 28 MeV theoretical error is understated and the central value is not robust to the u=1 renormalon assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The 28 MeV theoretical error in Eq. (66) is built from Eq. (61)'s 22 MeV higher-order estimate plus the scale and Z_m errors in Eq. (62). But the whole subtraction at μc=5 GeV removes only the leading u=1/2 infrared renormalon. Eq. (59) omits the mΩ_2 terminant and assigns the remainder O(μc e^{-2π/(β0 α)(1+ln2)}); in Sec. V.C the paper argues only that the f_n sign pattern and the lattice analysis of [38] are 'consistent' with the u=1 renormalon being zero, without deriving any bound on its normalization Z_2. If Z_2 is of typical O(1) size, the omitted term is of order Λ^2/μc ~ 10-30 MeV at μc=5 GeV, i.e. comparable to or larger than the quoted 28 MeV, and it shifts the central value rather than being covered by the scale-variation error (which probes the μc-dependence of the already-subtracted series) or by the Z_m error (which normalizes only the leading renormalon through Eq. (22)). The abstract's 28 MeV also excludes the +119/-123 MeV α_s error of Eq. (66); that omission is defensible only if α_s is treated as an external parameter, but the paper does not say so. The load-bearing issue is therefore that the headline uncertainty is conditional on an unquantified assumption about a debated renormalon.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs hyperasymptotic expansions for the heavy-quark pole mass regulated with the principal-value (PV) prescription, extending earlier work by the same group to include ultraviolet renormalons. After establishing the general formalism, the authors test it in the large-β0 approximation, where exact results are available, and demonstrate that successive hyperasymptotic corrections (mP, mΩm, subleading sums, mΩ−2) produce the expected convergence in both lattice and MS schemes. They then use the method to extract Λ̄_PV from quenched lattice data and from the B-meson mass, and compare with an independent lattice result. Finally, they construct a renormalon-free running function F(m,nf) to evolve the top mass from 163 GeV down to μc=5 GeV, decoupling bottom and charm quarks, and evaluate m_PV(μc)+μcΩm there. The central result is Eq. (66): m_{t,PV}(163 GeV)=173033 (+25/−28)(th) (+119/−123)(α) MeV, i.e., a 28 MeV theoretical uncertainty for fixed α_s, which the abstract quotes.","tokens_in":34794,"tokens_out":8909,"duration_ms":86799,"significance":"The paper is a serious, technically careful contribution. The large-β0 toy model is a genuine check of the hyperasymptotic expansion, and the consistency between lattice-scheme and MS-scheme determinations of Λ̄_PV is encouraging. The decoupling/renormalon-free-running strategy for the top mass is well motivated, and the bottom/charm finite-mass effects are treated in detail. If the error assessment is correct, the result would substantially reduce the often-quoted O(Λ_QCD) pole-mass ambiguity for the top. The main caveat is that the headline uncertainty is conditional on an unquantified assumption about the u=1 renormalon and on a heuristic truncation-error estimate; because the paper's central claim is precisely about the size of an uncertainty, these assumptions are load-bearing rather than cosmetic.","major_comments":[{"comment":"The abstract's δm_{t,PV}=28 MeV and the central value in Eq. (66) are obtained from Eq. (59), which omits the mΩ_2 terminant associated with a possible u=1 infrared renormalon and estimates the remainder as O(μc e^{−2π/(β0α)(1+ln2)}). At μc=5 GeV this remainder is numerically comparable to the quoted 28 MeV, so the quoted uncertainty is not a bound on this contribution but an assumption that its coefficient is small. Section V.C offers only plausibility arguments—the sign pattern of f_n in Table V and the lattice analysis of Ref. [38]—and no quantitative bound on the normalization Z_2. If Z_2 were of order one, the omitted mΩ_2 would shift the central value by an amount comparable to the quoted error and would not be covered by the μc-variation error in Eq. (62), which tests the already-subtracted series. Please provide a quantitative estimate or bound on mΩ_2, or explicitly present the 28 MeV as conditional on Z_2 being zero or negligible.","section":"§V.B, §V.C (Eqs. (59), (66))"},{"comment":"A major component of the theoretical error is the 22 MeV assigned in Eq. (61) to the truncation of F(m,nf), estimated as half the last computed term. This is a standard rule for sign-alternating asymptotic series, but the available series for nf=3 has only four coefficients, and the sign pattern is itself used in Section V.C as evidence about renormalons. Applying the 'half the last term' rule to such a short series is an assumption; it would be more convincing to test the stability of the error estimate under different truncation orders (e.g., excluding or including the last coefficient) and to compare with the large-β0 prediction, where the exact remainder is known.","section":"§V.B, Eq. (61)"}],"minor_comments":[{"comment":"The 28 MeV is the theory error for fixed α_s; the α_s contribution is +119/−123 MeV. Please state this explicitly in the abstract, since as written the abstract may be read as a total uncertainty.","section":"Abstract and Eq. (66)"},{"comment":"Please clarify how the quoted +25/−28 MeV theoretical error is obtained from the components +22/−22 (h.o.), +22/−22 (μ), +7/−15 (Z_m), and +9/−9 (μ_c) listed in Eq. (63); a direct quadrature of these components gives a somewhat larger interval.","section":"Eqs. (63) and (66)"},{"comment":"The unit '163MeV' in Eqs. (63) and (66) should read '163 GeV'.","section":"Eqs. (63) and (66)"},{"comment":"The word 'Elucubrative' after Eq. (42) is nonstandard; 'Speculatively' or 'As a speculative ansatz' is preferable. Similarly, 'aminorates' in Sec. IV.E should be 'ameliorates'.","section":"§IV.A and §IV.E"},{"comment":"The labels (a)–(d) in the lower panel are not legible at the printed size; consider enlarging the font or using a legend.","section":"Figure 11"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a substantial extension of the authors' earlier work and is within the scope of the journal. My main concern is that the headline 28 MeV number is likely to be cited without the caveats that appear in the body; a revision that quantifies or explicitly conditions the u=1 renormalon contribution would materially strengthen the paper. I do not see a circularity problem in the top-mass conversion: Z_m is taken from independent determinations and the top mass itself is not fitted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a read if you work on heavy quark masses and renormalons. It's the follow-up to their 2019 hyperasymptotic paper, now applied to the pole mass with ultraviolet renormalons included. The new deliverables are concrete: a lattice-based bar-Lambda_PV in quenched QCD (1.42 r0^-1), a bottom PV mass of 4.836 GeV, and a top conversion m_t,PV(163 GeV)=173033 MeV with a +25/-28 MeV theory error. The running/decoupling treatment for the top is new and looks careful; the bottom/charm finite-mass effects come out at the few MeV level.\n\nThe strongest part is the large-beta0 toy model. It shows the hyperasymptotic expansion converging down to surprisingly low scales, and it demonstrates that the u=-1 UV renormalon is suppressed in a high-scale scheme like the lattice. That gives real confidence the machinery works, not just formal manipulation.\n\nThe soft spots are real but not fatal. The 28 MeV headline is built from the truncation estimate for the F integral, taken as half the last computed term, plus scale and Z_m errors at mu_c=5 GeV. The half-last-term rule is a heuristic from sign-alternating series; fine, but not a theorem. More importantly, the subtraction at mu_c=5 GeV removes only the u=1/2 renormalon. The missing u=1 term is assigned to O(mu_c e^{-2 pi/(beta0 alpha)(1+ln2)}) but its normalization Z_2 is never bounded. Section V.C argues from the f_n sign pattern and the Fermilab/MILC lattice analysis that the u=1 renormalon is likely zero, but that is not a proof. If Z_2 were O(1), the omitted term would be order Lambda^2/mu_c ~ 10-30 MeV, the same size as the quoted error, and it would shift the central value. Scale variation does not obviously cover that, because it probes the mu_c dependence of the subtracted series, not the size of a term that was never included.\n\nThe abstract's 28 MeV also excludes the alpha_s uncertainty, which Eq. (66) gives as +119/-123 MeV. That is defensible if alpha_s is treated as an external input, but the abstract does not say that. A sentence would fix it.\n\nOn citation practice: they lean on their own earlier papers and on [13] for Z_m, but that is the standard source for the renormalon normalization, and their lattice result is consistent with earlier superasymptotic determinations. No red flags there.\n\nBottom line: the paper deserves peer review. A referee should push for an explicit statement that the 28 MeV assumes no u=1 renormalon and should ask for a bound or a dedicated error term if any evidence for Z_2 exists. The toy model and the lattice analysis are solid enough that the paper should be published after those clarifications.","headline":"A careful extension of the authors' hyperasymptotic method to the pole mass, with a plausible 28 MeV theory error for the top PV mass once you accept that the u=1 renormalon is absent; the abstract's headline omits the larger alpha_s error and the conditionality.","tokens_in":35220,"tokens_out":6357,"would_cite":true,"duration_ms":55434,"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 paper claims that the heavy-quark pole mass, when regulated by the principal-value prescription, is a well-defined mass, and that the MS-bar-to-PV conversion for the top quark carries a theoretical uncertainty of 28 MeV.","keywords":["heavy quark pole mass","hyperasymptotic expansion","principal value prescription","renormalons","top quark mass","MS-bar mass","Lambda-bar","B meson"],"falsifier":"Decisive tests: compute the renormalon normalization $Z_m$ independently, for instance from the static potential with reduced uncertainty or from a sum rule free of the leading renormalon, since $m\\Omega_m$ is proportional to $Z_m$ and a shift of $Z_m$ by its current uncertainty changes the $173033$ MeV central value by more than 20 MeV. Also determine whether the $u=1$ infrared renormalon normalization is nonzero: if it is, the missing terminant contributes an additional $O(\\Lambda_{\\rm QCD})$ shift in the conversion, and the 28 MeV theoretical error estimate understates the uncertainty.","tokens_in":34155,"feed_emoji":"⚛️","tokens_out":8817,"duration_ms":79326,"temperature":0.7,"pith_summary":"The paper claims that the longstanding 'ambiguity' of the heavy-quark pole mass is not an ambiguity at all: once the divergent perturbative series is resummed with the principal-value (PV) prescription, quantities such as the PV pole mass are well defined with arbitrary achievable accuracy, and the practical question is only how well the conversion between the MS-bar mass and a specific pole-mass definition is known. For the top quark with MS-bar mass 163 GeV, the paper finds this conversion to have a theoretical uncertainty of 28 MeV, with a central value $m_{t,\\rm PV} = 173033$ MeV. The argument runs through hyperasymptotic expansions, in which the truncated perturbative sum is supplemented by 'terminants' that encode factorial divergences (renormalons) in the Borel plane. The method is tested in the large-$\\beta_0$ approximation, applied to lattice and $B/D$-meson determinations of $\\bar\\Lambda$, and then used for the top quark.","feed_headline":"Top quark's pole-mass 'ambiguity' shrinks to 28 MeV","feed_subtitle":"A Borel-resummed principal-value mass replaces the pole-mass ambiguity with a computable 28 MeV theoretical error.","key_machinery":"The central object is the hyperasymptotic expansion of the pole mass, a Borel-resummed principal-value version of the divergent perturbative series. The workhorse is the terminant $\\Omega_d$ (and its mass analogue $m\\Omega_m$), the completion that restores the part of the series associated with a singularity at $u=d/2$ in the Borel plane; for the leading infrared renormalon at $u=1/2$ it has the form $\\Omega_m \\sim \\sqrt{\\alpha(\\mu)}\\, K_X^{(P)}\\, (\\mu/m)\\, e^{-2\\pi/(\\beta_0\\alpha(\\mu))}\\, (\\beta_0\\alpha(\\mu)/4\\pi)^{-b}$, with $K_X^{(P)}$ proportional to the normalization $Z_m$. The argument is carried by matching this terminant to the known asymptotic behaviour of the perturbative coefficients $r_n$, and by a renormalon-free running function $F(m,n_f)$ that lowers the top quark mass to scales around 5 GeV where the hyperasymptotic expansion can be evaluated, plus explicit decoupling of bottom and charm finite-mass effects.","core_discovery":"The central claim is that the PV-regulated pole mass $m_{\\rm PV}(m)$ can be computed from the MS-bar mass $m$ by a hyperasymptotic expansion whose leading non-perturbative term is the terminant $m\\Omega_m$, proportional to $\\sqrt{\\alpha}\\,(\\Lambda_{\\rm QCD}/m)$ times the renormalon normalization $Z_m$, and that this expansion is accurate enough that for the top quark the relation $\\bar m_t = 163$ GeV to $m_{t,\\rm PV}$ has a theoretical error of only 28 MeV (plus a much larger $+119/-123$ MeV from the strong coupling). A corollary is that the often-quoted pole-mass ambiguity of order $\\Lambda_{\\rm QCD}$ is replaced by a computable, definition-dependent difference: different legitimate pole-mass definitions differ by $O(\\Lambda_{\\rm QCD})$ with an arbitrary coefficient, so statements about 'the' pole mass ambiguity are ill-posed without a specific definition. The paper also finds evidence for the ultraviolet renormalon at $u=-1$ in the sign-alternating coefficients of the renormalon-free running function $F(m,n_f)$, while finding no evidence for an infrared renormalon at $u=1$.","pith_inferences":["If the $u=1$ infrared renormalon is truly absent, the PV prescription becomes a particularly clean mass definition for electroweak-precision fits, because the dominant theoretical error would then come from $\\alpha_s$ and higher-order coefficients, both improvable.","Because the large-$\\beta_0$ tests show the OPE picture holding down to scales around 667 MeV, the same machinery might give charm quark mass determinations with controlled errors despite the lower scale.","The observed sign-alternating coefficients in $F(m,n_f)$ imply that the $u=-1$ ultraviolet renormalon, usually thought to be subleading, may set the practical truncation point in the MS scheme; a lattice-scheme computation with a high effective scale could push the expansion further.","Relating the PV mass to experimentally measured top cross sections rather than assuming the measured mass is the pole mass would remove an $O(\\Lambda_{\\rm QCD})$ theoretical ambiguity from top-quark mass extractions; the paper makes this concrete by giving the PV-to-MS-bar shift with a 28 MeV error."],"forward_implications":["The top quark pole mass can be defined and used in a scheme-and-scale controlled way: with current perturbative input, $m_{t,\\rm PV}(163\\,{\\rm GeV}) = 173033^{+25}_{-28}({\\rm th})^{+119}_{-123}(\\alpha)$ MeV, so the purely theoretical uncertainty is well below the typical $\\Lambda_{\\rm QCD}$ ambiguity quoted in the literature.","The same hyperasymptotic machinery gives $m_{b,\\rm PV} = 4836^{+8}_{-17}(Z_m)^{+12}_{-11}(\\alpha)^{+8}_{-9}$ MeV for the bottom quark and $\\bar\\Lambda_{\\rm PV} = 477\\pm 46$ MeV from the $B$ meson mass, allowing cross-checks between lattice and $B$-physics determinations.","The renormalon-free running function $F(m,n_f)$ shows a sign-alternating pattern consistent with an ultraviolet renormalon at $u=-1$ and no sign of an infrared renormalon at $u=1$; if the latter is truly absent, the heavy-quarkonium mass analysis closes consistently.","The error in the conversion is dominated by the strong coupling, not by renormalon ambiguity, and it is systematically improvable order by order in perturbation theory."],"supporting_citations":[{"why":"Establishes the hyperasymptotic expansion framework for observables with an OPE and the PV summation that this paper extends to include ultraviolet renormalons.","marker":"[1]"},{"why":"Supplies the renormalon normalization $Z_m$ and its uncertainty that set the size of the leading terminant $m\\Omega_m$.","marker":"[13]"},{"why":"Provides the Borel transform of the pole mass in the large-$\\beta_0$ approximation, the starting point for identifying the $u=1/2$ and $u=-1$ renormalons.","marker":"[19]"},{"why":"Gives the large-$\\beta_0$ pole-mass Borel representation and the divergent-series structure used to build the asymptotic coefficients $r_n^{(as)}$.","marker":"[20]"},{"why":"Supplies the high-order lattice perturbative coefficients $c_n$ and the renormalon-based determination of $Z_m$ in the lattice scheme used for $\\bar\\Lambda_{\\rm PV}$ fits.","marker":"[22]"},{"why":"Provides the lattice-scheme normalization $Z_{\\rm latt}^m = 17.9(1.0)$ and beta-function estimates that set the $O(1/n)$ corrections in the asymptotic coefficients.","marker":"[23]"},{"why":"The alternative top pole-mass ambiguity estimate that this paper compares against and finds less sensitive to $Z_m$.","marker":"[45]"},{"why":"Supplies the four-loop coefficient of the pole-to-MS-bar relation used in the running and in the large-$\\beta_0$ comparison.","marker":"[9]"}],"fun_headline_variants":["Hyperasymptotics trim top pole-mass error to 28 MeV","Top pole mass uncertainty: 28 MeV via hyperasymptotics","PV-resummed pole mass: top error just 28 MeV","Hyperasymptotic expansion sets top mass error at 28 MeV","Pole-mass ambiguity resolved to 28 MeV for top quark"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the size of the leading 'runaway growth' of the perturbative series—the renormalon normalization—is known to the accuracy of the cited determination, and that there is no competing divergence at twice that location; if either assumption fails, the central value and the 28 MeV error both move.","fun_headline_variants_meta":{"raw":{"variants":["Hyperasymptotics trim top pole-mass error to 28 MeV","Top pole mass uncertainty: 28 MeV via hyperasymptotics","PV-resummed pole mass: top error just 28 MeV","Hyperasymptotic expansion sets top mass error at 28 MeV","Pole-mass ambiguity resolved to 28 MeV for top quark"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001076,"raw_usage":{"total_tokens":4493,"prompt_tokens":927,"completion_tokens":3566,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":3476}},"tokens_in":543,"tokens_out":3566,"duration_ms":25298,"temperature":1.0,"reasoning_tokens":3476,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:19:49.480608+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Decisive tests: compute the renormalon normalization $Z_m$ independently, for instance from the static potential with reduced uncertainty or from a sum rule free of the leading renormalon, since $m\\Omega_m$ is proportional to $Z_m$ and a shift of $Z_m$ by its current uncertainty changes the $173033$ MeV central value by more than 20 MeV. Also determine whether the $u=1$ infrared renormalon normalization is nonzero: if it is, the missing terminant contributes an additional $O(\\Lambda_{\\rm QCD})$ shift in the conversion, and the 28 MeV theoretical error estimate understates the uncertainty.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the hyperasymptotic expansion framework for observables with an OPE and the PV summation that this paper extends to include ultraviolet renormalons."},{"cited_title":"QCD perturbation theory at large orders with large renormalization scales in the large $\\beta_0$ limit","cited_arxiv_id":"hep-ph/0307070","evidence_quote":"Provides the Borel transform of the pole mass in the large-$\\beta_0$ approximation, the starting point for identifying the $u=1/2$ and $u=-1$ renormalons."}],"review_version":1}