arxiv: 2604.05044 · v1 · submitted 2026-04-06 · ✦ hep-ph · hep-th

Recognition: 2 theorem links

· Lean Theorem

The chromomagnetic moment of a heavy quark with hyperasymptotic precision

Cesar Ayala , Antonio Pineda

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:14 UTC · model grok-4.3

classification ✦ hep-ph hep-th

keywords infrared renormalonchromomagnetic momentheavy quarkhyperfine splittingB mesonD mesonoperator product expansionperturbative series

0 comments

The pith

The leading infrared renormalon normalization for the chromomagnetic moment of a heavy quark is fixed at 0.507(7) GeV² by fitting hyperfine splittings in B and D mesons.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper determines the normalization of the leading infrared renormalon in the perturbative series for the chromomagnetic moment of a heavy quark. It incorporates the associated terminant into a hyperasymptotic expansion to compute the hyperfine splitting of ground-state B and D mesons. Experimental data are then fitted to the operator product expansion prediction, leaving hat mu squared G PV as the single free parameter. This yields the value 0.507(7) GeV². A reader would care because the result supplies a concrete non-perturbative input for heavy-quark systems that improves the reliability of meson mass and splitting predictions.

Core claim

We determine the normalization of the leading infrared renormalon of the chromomagnetic moment of a heavy quark. Estimates of higher order coefficients of the perturbative series are given. We compute the hyperfine splitting of the B and D mesons for the ground state with hyperasymptotic precision by including the leading terminant, associated with the first infrared renormalon. We fit the experimental data to the operator product expansion theoretical prediction with hat mu²_{G,PV} as the free parameter. We obtain hat mu²_{G,PV}=0.507(7) GeV² for the ground state.

What carries the argument

The leading terminant associated with the first infrared renormalon in the perturbative expansion of the chromomagnetic moment, which is added to produce the hyperasymptotic approximation.

If this is right

The extracted normalization supplies estimates for all higher-order coefficients in the perturbative series of the chromomagnetic moment.
The same hyperasymptotic method can be applied to other observables governed by the same infrared renormalon structure.
The value 0.507(7) GeV² can be inserted directly into operator product expansion predictions for additional heavy-meson properties.
The approach yields a concrete benchmark against which alternative resummation schemes for the same series can be tested.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The extracted parameter offers a target for direct comparison with non-perturbative lattice determinations of the chromomagnetic operator matrix element.
The method suggests that analogous terminant insertions could reduce theoretical uncertainties in related heavy-quark effective theory observables such as decay rates or form factors.
Consistency of the fitted value across different meson systems would strengthen the assumption that the leading renormalon dominates the series.
The result supplies a fixed input that can be propagated into global fits of heavy-quark parameters used in flavor physics.

Load-bearing premise

The perturbative series for the chromomagnetic moment is dominated by its leading infrared renormalon, allowing the terminant to produce a reliable approximation whose only remaining free parameter is the fitted value of hat mu²_{G,PV}.

What would settle it

A new lattice QCD computation or experimental measurement of the ground-state hyperfine splitting in B or D mesons that lies more than two standard deviations outside the range predicted by the fitted value 0.507(7) GeV².

Figures

Figures reproduced from arXiv: 2604.05044 by Antonio Pineda, Cesar Ayala.

**Figure 2.** Figure 2: Plot of the scale dependence of the ratio of [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗

**Figure 3.** Figure 3: Determination of ZLS for nl = 0 and nl = 3. The dotted, dashed and solid black lines correspond to Eq. (31) for n = 0, 1 and 2 respectively, where x = ν/mPV. The dash-dotted red line corresponds to Z fit LS cn(x) c (as) n,fit(x) using the normalizations obtained in Eqs. (32) and (33). The horizontal line and the blue band yield our prediction for ZLS and the error estimate listed in [PITH_FULL_IMAGE:figur… view at source ↗

**Figure 4.** Figure 4: Charm quark case. Upper figure: Determinations of c (c) PV at O(αs) (dotted-black line), O(α 2 s ) (dashed-black line), O(α 3 ) (dash-dotted black line), at order (1,0) in the hyperasymptotic expansion counting (continuous-black line), and at order (1,N = 2) in the hyperasymptotic expansion counting (dash-dotted red line). Lower figure: Zoom of the previous figure adding the determination of c (c) PV at … view at source ↗

**Figure 5.** Figure 5: Bottom quark case. Upper figure: Determinations of c (b) PV at O(αs) (dotted-black line), O(α 2 s ) (dashed-black line), O(α 3 ) (dash-dotted black line), at order (1,0) in the hyperasymptotic expansion counting (continuous-black line), and at order (1,N = 1) in the hyperasymptotic expansion counting (dash-dotted red line). Lower figure: Zoom of the previous figure adding the determination of c (b) PV at… view at source ↗

**Figure 6.** Figure 6: Determination of ˆc (b) PV/cˆ (c) PV at order (1,0) in the hyperasymptotic expansion counting (continuous black line), at order (1,N) (where N = 1 for ˆc (b) PV and N = 2 for ˆc (c) PV) in the hyperasymptotic expansion counting (dash-dotted red line), and at order (1,0) in the hyperasymptotic expansion counting (dotted-red line) using the normalization of the renormalons obtained in Eq. 33. The horizonta… view at source ↗

read the original abstract

We determine the normalization of the leading infrared renormalon of the chromomagnetic moment of a heavy quark. Estimates of higher order coefficients of the perturbative series are given. We compute the hyperfine splitting of the $B$ and $D$ mesons for the ground state with hyperasymptotic precision by including the leading terminant, associated with the first infrared renormalon. We fit the experimental data to the operator product expansion theoretical prediction with $\hat \mu^2_{G,\rm PV}$ as the free parameter. We obtain $\hat \mu^2_{G,\rm PV}=0.507(7)$ GeV$^2$ for the ground state.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper fits the chromomagnetic matrix element after adding the leading renormalon terminant to the hyperfine splitting series, giving a data-driven value rather than an ab initio result.

read the letter

The main point is that the authors fix the normalization of the leading infrared renormalon for the chromomagnetic moment operator and then use the associated terminant to reach hyperasymptotic precision on the ground-state hyperfine splittings of B and D mesons. They fit the remaining parameter hat mu^2_G,PV directly to experimental data and report 0.507(7) GeV². That fitted number is the central output. They also supply estimates for higher-order perturbative coefficients in the series for the moment itself. The work sits squarely inside standard renormalon resummation techniques applied to heavy-quark effective theory. It improves control over the perturbative piece of the operator product expansion for these splittings and produces a concrete number with a small quoted uncertainty. The approach is internally consistent with how the field handles IR renormalons in similar operators. The soft spot is that the result is obtained by treating the matrix element as the single free parameter and adjusting it to match data. This makes the quoted central value a fit by construction, not an independent theoretical prediction. The claim of reliable hyperasymptotic accuracy therefore rests on the assumption that the leading renormalon dominates the remaining series once the terminant is included. That assumption is common but would be stronger with explicit checks on subleading renormalon effects or comparisons to lattice results. The paper is aimed at people who need precise non-perturbative inputs for heavy-meson spectroscopy and decays inside HQET. Readers who already work with renormalon methods or who compare effective-theory predictions to B and D data will get direct use from the number and the method. It is worth sending for peer review. The technical steps are clear, the application is focused, and the result supplies a usable input even if the value itself comes from the fit.

Referee Report

1 major / 2 minor

Summary. The paper determines the normalization of the leading infrared renormalon in the chromomagnetic moment of a heavy quark. It provides estimates of higher-order perturbative coefficients, computes the hyperfine splitting of ground-state B and D mesons via an operator product expansion that includes the leading terminant for hyperasymptotic precision, and fits the resulting theoretical expression to experimental data with a single free parameter, obtaining the value hat mu^2_{G,PV} = 0.507(7) GeV^2.

Significance. If the assumption that the perturbative series is dominated by the leading infrared renormalon holds and the terminant yields a reliable approximation, the result supplies a high-precision value for a key matrix element in heavy-quark effective theory. This can serve as input for meson spectroscopy, lattice QCD benchmarks, and other phenomenological applications, while illustrating the practical use of hyperasymptotic expansions to control renormalon ambiguities in QCD.

major comments (1)

The central result is obtained by treating hat mu^2_{G,PV} as a free parameter and fitting it to hyperfine data after including the terminant (abstract). While this is a standard procedure for extracting non-perturbative parameters, the quoted uncertainty of 0.007 GeV^2 requires explicit demonstration that the remainder after the leading terminant (including possible subleading renormalon contributions) lies below this level; the abstract mentions estimates of higher-order coefficients but does not quantify the associated truncation error.

minor comments (2)

The abstract states that estimates of higher-order coefficients are given, but the main text should tabulate or plot these estimates alongside the known low-order terms to allow readers to assess the convergence pattern before the terminant is added.
Notation for the principal-value prescription and the definition of the terminant should be introduced with a brief equation in the introductory section on the OPE to improve readability for readers unfamiliar with hyperasymptotic techniques.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the recommendation of minor revision. We address the single major comment below and will revise the manuscript accordingly to strengthen the justification of our quoted uncertainty.

read point-by-point responses

Referee: The central result is obtained by treating hat mu^2_{G,PV} as a free parameter and fitting it to hyperfine data after including the terminant (abstract). While this is a standard procedure for extracting non-perturbative parameters, the quoted uncertainty of 0.007 GeV^2 requires explicit demonstration that the remainder after the leading terminant (including possible subleading renormalon contributions) lies below this level; the abstract mentions estimates of higher-order coefficients but does not quantify the associated truncation error.

Authors: We agree that an explicit demonstration is needed to support the quoted uncertainty. The manuscript already supplies estimates of higher-order perturbative coefficients (Section 3) that indicate dominance by the leading infrared renormalon, and the hyperasymptotic calculation incorporates the associated leading terminant. However, we acknowledge that a dedicated quantitative bound on the remainder (including subleading renormalons) was not provided. In the revised version we will add a concise paragraph in the hyperfine-splitting section that converts the coefficient estimates into an explicit upper limit on the truncation error after the leading terminant, showing that this contribution lies below 0.007 GeV² and thereby justifying the precision of the fitted value. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper explicitly computes perturbative coefficients and the leading IR renormalon normalization for the chromomagnetic moment, incorporates the associated terminant into the hyperasymptotic OPE for hyperfine splittings, and then fits the single remaining parameter hat mu^2_{G,PV} directly to B/D meson experimental data, reporting the fitted value as obtained from that comparison. This is a transparent parameter extraction using the theoretical framework rather than a first-principles prediction that reduces to its own inputs by construction. No self-definitional relations, mislabeled fitted quantities presented as predictions, load-bearing self-citations, or ansatz smuggling appear in the derivation chain. The approach is self-contained against external experimental benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The work rests on standard QCD assumptions about the existence and dominance of infrared renormalons plus the validity of the operator-product expansion for hyperfine splittings; the only explicit free parameter is the fitted scale.

free parameters (1)

hat mu^2_{G,PV} = 0.507(7) GeV^2
Single free parameter adjusted to match experimental hyperfine splittings of ground-state B and D mesons.

axioms (1)

domain assumption The leading infrared renormalon dominates the large-order behavior of the perturbative series for the chromomagnetic moment.
Used to justify truncation plus terminant resummation for hyperasymptotic accuracy.

pith-pipeline@v0.9.0 · 5404 in / 1425 out tokens · 35957 ms · 2026-05-10T19:14:57.718167+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We determine the normalization of the leading infrared renormalon of the chromomagnetic moment... hyperasymptotic precision by including the leading terminant... fit ... with μ̂²_G,PV as the free parameter.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Borel transform near the closest infrared renormalon singularity... Z_LS ... terminant

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Renormalon-based resummation for B(D) Mesons
hep-ph 2026-05 unverdicted novelty 4.0

Renormalon resummation with a holomorphic timelike coupling yields estimates for b and c quark pole masses and related hadronic parameters from MS-bar inputs and known series coefficients.

Reference graph

Works this paper leans on

32 extracted references · 25 canonical work pages · cited by 1 Pith paper

[1]

Heavy Quark Symmetry

M. Neubert, Phys. Rept.245, 259-396 (1994) [arXiv:hep-ph/9306320 [hep-ph]]

work page Pith review arXiv 1994
[2]

I. I. Y. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D52, 196-235 (1995) [arXiv:hep-ph/9405410 [hep-ph]]

work page arXiv 1995
[3]

M. E. Luke and A. V. Manohar, Phys. Lett. B286, 348-354 (1992) [arXiv:hep- ph/9205228 [hep-ph]]

work page arXiv 1992
[4]

Amoros, M

G. Amoros, M. Beneke and M. Neubert, Phys. Lett. B401, 81-90 (1997) [arXiv:hep- ph/9701375 [hep-ph]]

work page arXiv 1997
[5]

HQET chromomagnetic interaction at two loops,

A. Czarnecki and A. G. Grozin, Phys. Lett. B405, 142-149 (1997) [erratum: Phys. Lett. B650, 447 (2007)] [arXiv:hep-ph/9701415 [hep-ph]]

work page arXiv 1997
[6]

A. G. Grozin, P. Marquard, J. H. Piclum and M. Steinhauser, Nucl. Phys. B789, 277-293 (2008) [arXiv:0707.1388 [hep-ph]]

work page arXiv 2008
[7]

A. I. Davydychev and A. G. Grozin, Phys. Rev. D59, 054023 (1999) [arXiv:hep- ph/9809589 [hep-ph]]

work page arXiv 1999
[8]

Peset, A

C. Peset, A. Pineda and J. Segovia, JHEP09, 167 (2018) [arXiv:1806.05197 [hep-ph]]

work page arXiv 2018
[9]

The bottom quark mass from theΥ(1S) system at NNNLO,

C. Ayala, G. Cvetiˇ c and A. Pineda, JHEP09, 045 (2014) [arXiv:1407.2128 [hep-ph]]

work page arXiv 2014
[10]

A. G. Grozin and M. Neubert, Nucl. Phys. B508, 311-326 (1997) [arXiv:hep-ph/9707318 [hep-ph]]

work page arXiv 1997
[11]

Parisi, Phys

G. Parisi, Phys. Lett. B76, 65 (1978)

1978
[12]

A. H. Mueller, Phys. Lett. B308, 355 (1993)

1993
[13]

’t Hooft, inThe Whys of Subnuclear Physics, edited by A

G. ’t Hooft, inThe Whys of Subnuclear Physics, edited by A. Zichichi (Plenum, New York, 1978). 25

1978
[14]

Superasymptotic and hyperasymptotic approximation to the operator product expansion,

C. Ayala, X. Lobregat and A. Pineda, Phys. Rev. D99, no.7, 074019 (2019) [arXiv:1902.07736 [hep-th]]

work page arXiv 2019
[15]

Takaura, JHEP10, 039 (2020) [arXiv:2002.00428 [hep-ph]]

H. Takaura, JHEP10, 039 (2020) [arXiv:2002.00428 [hep-ph]]

work page arXiv 2020
[16]

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147, 385-447 (1979)

1979
[17]

G. S. Bali, C. Bauer and A. Pineda, Phys. Rev. Lett.113, 092001 (2014) [arXiv:1403.6477 [hep-ph]]

work page arXiv 2014
[18]

Ayala, X

C. Ayala, X. Lobregat and A. Pineda, JHEP12, 093 (2020) [arXiv:2009.01285 [hep-ph]]

work page arXiv 2020
[19]

Ayala, X

C. Ayala, X. Lobregat and A. Pineda, Nucl. Part. Phys. Proc.309-311, 77-86 (2020) [arXiv:1910.04090 [hep-ph]]

work page arXiv 2020
[20]

Hyperasymptotic approximation to the top, bottom and charm pole mass,

C. Ayala, X. Lobregat and A. Pineda, Phys. Rev. D101, no.3, 034002 (2020) [arXiv:1909.01370 [hep-ph]]

work page arXiv 2020
[21]

M. V. Berry and C. J. Howls, Hyperasymptotics, Proc. Roy. Soc. London A, 430 (1990), pp. 653-668

1990
[22]

Dingle,Asymptotic Expansions: Their Derivation and Interpretation(Academic Press, London, 1973)

R.B. Dingle,Asymptotic Expansions: Their Derivation and Interpretation(Academic Press, London, 1973)

1973
[23]

Ayala, X

C. Ayala, X. Lobregat and A. Pineda, JHEP09, 016 (2020) [arXiv:2005.12301 [hep-ph]]

work page arXiv 2020
[24]

Ayala and A

C. Ayala and A. Pineda, Phys. Rev. D106, no.5, 056023 (2022) [arXiv:2208.07389 [hep-ph]]

work page arXiv 2022
[25]

G. S. Bali, C. Bauer, A. Pineda and C. Torrero, Phys. Rev. D87, 094517 (2013) [arXiv:1303.3279 [hep-lat]]

work page arXiv 2013
[26]

Ayala and A

C. Ayala and A. Pineda, JHEP10, 176 (2025) [arXiv:2505.18827 [hep-ph]]

work page arXiv 2025
[27]

Navaset al.[Particle Data Group], Phys

S. Navaset al.[Particle Data Group], Phys. Rev. D110, no.3, 030001 (2024)

2024
[28]

Version 3 of {\tt RunDec} and {\tt CRunDec}

F. Herren and M. Steinhauser, Comput. Phys. Commun.224, 333-345 (2018) [arXiv:1703.03751 [hep-ph]]

work page Pith review arXiv 2018
[29]

Extraction of heavy-quark-expansion parameters from unquenched lattice data on pseudoscalar and vector heavy-light meson masses,

P. Gambino, A. Melis and S. Simula, Phys. Rev. D96, no.1, 014511 (2017) [arXiv:1704.06105 [hep-lat]]

work page arXiv 2017
[30]

Bazavovet al.(Fermilab Lattice, MILC, TUMQCD), Up-, down-, strange-, charm-, and bottom-quark masses from four-flavor lattice QCD, Phys

A. Bazavovet al.[Fermilab Lattice, MILC and TUMQCD], Phys. Rev. D98, no.5, 054517 (2018) [arXiv:1802.04248 [hep-lat]]

work page arXiv 2018
[31]

Hayashi, Y

Y. Hayashi, Y. Sumino and H. Takaura, Phys. Lett. B829, 137068 (2022) [arXiv:2202.01434 [hep-ph]]

work page arXiv 2022
[32]

Renormalon-based resummation for spacelike and timelike QCD quantities whose perturbation expansion has general form,

C. Ayala, G. Cvetiˇ c and R. K¨ ogerler, J. Phys. G53(2026) no.2, 025005 [arXiv:2509.09042 [hep-ph]]. 26

work page arXiv 2026