---
title: Renormalon-Free Gluon Condensate Scheme
url: https://www.emergentmind.com/topics/renormalon-free-gluon-condensate-scheme
type: topic
---

# Renormalon-Free Gluon Condensate Scheme

The renormalon-free gluon-condensate scheme is a class of prescriptions for defining the dimension-four gluon condensate in a way that removes the \(O(\Lambda^4)\) ambiguity associated with the leading \(u=2\) infrared renormalon of perturbative Wilson coefficients. In these schemes, the perturbative contribution and the condensate term are rearranged so that the operator-product expansion (OPE) becomes unambiguous term by term, and the condensate is promoted from a prescription-dependent remainder to a well-defined nonperturbative matrix element. Concrete realizations include large-\(\beta_0\) dispersive constructions, minimal-term and principal-value prescriptions for the lattice plaquette, pNRQCD factorization for the static potential, \(R\)-dependent subtraction schemes for \(\tau\) observables, and gradient-flow formulations in which the flow time \(1/\sqrt{t}\) serves as a gauge-invariant factorization scale [1807.10064, 1502.00086, 2009.01285, 1712.05435, 2202.10957, 2309.10853].

## 1. Renormalon ambiguity at dimension four

In a standard single-scale OPE,
\[
X(Q^2)=c_{1,X}(Q^2)\cdot \langle 1\rangle
+c_{FF,X}(Q^2)\cdot \frac{\langle (\alpha/\pi)\{F^2\}\rangle}{Q^4}
+O(1/Q^6),
\]
the coefficient \(c_{1,X}\) is computed in perturbation theory but is afflicted by the \(u=2\) infrared renormalon, which induces an intrinsic uncertainty of order \((\Lambda^2/Q^2)^2\sim \Lambda^4/Q^4\), precisely the same order as the condensate term. This prevents a reliable determination of the gluon condensate, because any subtraction of \(c_{1,X}\) at finite order leaves an \(O(\Lambda^4/Q^4)\) remnant [1807.10064].

The same structural problem appears in lattice regularization. For the plaquette,
\[
P(a)=P_{\rm pert}(a)+C_G(a)\,\langle G^2\rangle+\cdots,
\]
the perturbative tail is asymptotic, and the identity operator mixes with the plaquette, producing a power-divergent contribution \(Z(\beta)\sim a^{-4}\). In this setting, the separation of scales in the OPE does not correspond to a separation of perturbative and non-perturbative contributions, so the ambiguity attached to the perturbative series has to be absorbed into the definition of the condensate itself [1811.05427].

For the Adler function in the conventional \(\overline{\rm MS}\) OPE, the leading non-perturbative correction is the dimension-4 gluon condensate, while the perturbative coefficients contain the same \(u=2\) singularity that generates factorial growth and an \(O(\Lambda_{\rm QCD}^4)\) ambiguity in the standard condensate. This is the ambiguity targeted by renormalon-free gluon-condensate schemes [2202.10957].

## 2. General subtraction logic and scheme architectures

A common strategy is to introduce an auxiliary scale or prescription that isolates the renormalon-sensitive part of the perturbative coefficient and then absorb the matching ambiguity into a redefined condensate. In the large-\(\beta_0\) framework of Suzuki and Takaura, the perturbative Wilson coefficient is written in a dispersive form,
\[
c_{1,X}(Q^2)=\int_0^\infty \frac{d\tau}{2\pi\,\tau}\,w_X(\tau/Q^2)\,\alpha_{\beta_0}(\tau),
\]
with
\[
\alpha_{\beta_0}(\tau)=\frac{4\pi/\beta_0}{\ln(\tau/(e^{5/3}\Lambda^2))}.
\]
Introducing an infrared cutoff \(\mu_f\) yields \(c_{1,X}(Q^2;\mu_f)\), whose \(\mu_f\)-dependent part can be matched against the ultraviolet-cutoff dependence of the condensate. This leads to
\[
\langle(\alpha/\pi)\{F^2\}\rangle^{\rm RF}
\equiv
\langle(\alpha/\pi)\{F^2\}\rangle(\mu_f)
-
A\,{\rm Im}\!\int_{C_b}\frac{d\tau}{\pi\,\tau}\;i\,\tau^2\,\alpha_{\beta_0}(\tau),
\]
which is \(\mu_f\)-independent, free of \(O(\Lambda^4)\) renormalon ambiguity, and universal in the sense that no observable label \(X\) appears in its definition [1807.10064].

Different implementations realize the same logic in different variables.

| Framework | Observable | Characteristic subtraction |
|---|---|---|
| Large-\(\beta_0\) OPE [1807.10064] | \(X(Q^2)\) | \(\mu_f\)-dependent split of \(c_{1,X}\) and \(\langle(\alpha/\pi)\{F^2\}\rangle\) |
| Plaquette hyperasymptotics [2009.01285] | Wilson plaquette | PV Borel sum \(P_{\rm PV}\) plus terminant \(\Omega_{G^2}\) |
| Minimal-term lattice schemes [1502.00086, 1811.05427] | Plaquette | truncation at the minimal term \(S_P\) |
| pNRQCD static potential [1712.05435] | \(V_{\rm QCD}(r)\) | separation into \(V_S^{\rm RF}\), \(\delta E_{\rm us}^{\rm RF}\), and \(\langle g^2G^2\rangle_{\rm RF}\) |
| \(R\)-scheme for \(\tau\) moments [2202.10957] | Adler function and moments | subtraction \(R^4\sum_\ell N_g r_\ell^{(4,0)}\bar a_R^\ell\) |
| Gradient flow [2309.10853, 2510.12193] | Adler function, \(\tau\) width | subtraction with flowed \(1/t^2\) coefficient \( \tilde C_1(t)\) |

These formulations differ in the auxiliary quantity that carries the subtraction—\(\mu_f\), the truncation point, a principal-value prescription, the pNRQCD cutoffs \(\mu_1,\mu_2\), the infrared factorization scale \(R\), or the gradient-flow time \(t\)—but each aims to cancel the same dimension-four ambiguity between perturbation theory and the condensate term.

## 3. Plaquette-based realizations on the lattice

The plaquette has supplied the most explicit numerical determinations of a renormalon-free gluon condensate. In four-dimensional SU(3) pure gauge theory with the Wilson action, the perturbative expansion
\[
P_{\rm pert}(\alpha)=\sum_{n=0}^\infty p_n\,\alpha^{n+1}
\]
was computed to 35 loops. The large-order coefficients exhibit the expected asymptotic behavior associated with the dimension-4 operator at \(u=2\), and the natural prescription in early lattice work was truncation at the minimal term. Defining
\[
S_P(\alpha)=\sum_{j=0}^{n_0(\alpha)} p_j\,\alpha^{j+1},
\]
with \(n_0\) chosen so that \(p_{n_0}\alpha^{n_0+1}\) is minimal, one extracts
\[
\langle G^2\rangle
=
\frac{36}{\pi^2 a^4}\,C_G^{-1}(\alpha)\,
\bigl[\langle P\rangle_{\rm MC}(a)-S_P(\alpha)\bigr]
+O(a^2).
\]
Using high-precision Monte Carlo data for \(\langle P\rangle\) and perturbative coefficients up to \(n=35\), Bali et al. obtained
\[
r_0^4\,\langle G^2\rangle = 3.18(29),
\qquad
\langle G^2\rangle = 0.077(8)\,{\rm GeV}^4,
\]
with an intrinsic “renormalon-prescription” ambiguity
\[
\delta\langle G^2\rangle = 27(11)\Lambda^4 \simeq 0.087\,{\rm GeV}^4
\]
in pure gluodynamics [1502.00086].

A later NSPT study extended the same logic to \(N_f=2\) massless QCD with massless staggered fermions and twisted boundary conditions. In that computation, the coefficients \(p_n\) were measured up to \(n\simeq 35\), the ratios \(p_n/(n\,p_{n-1})\) approached the renormalon-predicted constant, and the minimal-term subtraction yielded a condensate that exhibited clean \(a^4\) scaling and remained stable under modest variations of the prescription point \(\bar n\) or of the scale \(\mu\). The quoted result was
\[
\left\langle \frac{\alpha_s}{\pi}G^2\right\rangle_{\rm RF}
=
(0.015\pm0.003_{\rm stat}\pm0.004_{\rm syst})\,{\rm GeV}^4
\]
for \(N_f=2\) massless QCD [1811.05427].

The hyperasymptotic formulation of Ayala et al. replaced optimal truncation by a principal-value Borel prescription and an explicit terminant associated with the leading renormalon. The perturbative sum is defined as
\[
P_{\rm PV}(g)\equiv {\rm PV}\int_0^\infty du\, e^{-u/[\beta_0 g^2]}\,B[P](u),
\]
and the OPE becomes
\[
\langle P\rangle_{\rm MC}
=
P_{\rm PV}(\alpha)
+
(\pi^2/36)\,C_G(\alpha)\,a^4\,\langle G^2\rangle_{\rm PV}
+
O(a^6).
\]
Using the hyperasymptotic expansion
\[
S_{\rm PV}(\alpha)=
\sum_{n=0}^{N_P}p_n\alpha^{n+1}
+\Omega_{G^2}(\alpha)
+\sum_{n=N_P+1}^{N'} [p_n-p_n^{(\rm as)}]\alpha^{n+1}
+\cdots,
\]
they obtained
\[
\langle G^2\rangle_{\rm PV}(n_f=0)=3.15(18)\,r_0^{-4}.
\]
The final result showed a very flat dependence on \(a^2\), repeating the analysis in the \(\overline{\rm MS}\) scheme gave the same central value within the quoted error, and the scheme and scale dependence in \(C_G\), \(\Lambda\), and \(Z_P\) canceled in the combination \(\langle P\rangle_{\rm MC}-P_{\rm PV}\) [2009.01285].

## 4. pNRQCD and the static QCD potential

In pNRQCD, the static QCD potential is organized as
\[
V_{\rm QCD}(r)=V_S(r)+\delta E_{\rm us}(r)+\cdots.
\]
The soft singlet potential \(V_S\) contains renormalons at positive half-integers \(u=1/2,3/2,5/2,\ldots\), and the leading \(r\)-dependent ambiguity is the \(u=3/2\) renormalon \(O(\Lambda^3 r^2)\). The next-to-leading ultrasoft term \(\delta E_{\rm us}(r)\), which scales as \(\alpha_s^3 r^2\) at small \(r\), contains the same \(u=3/2\) renormalon. With two cutoffs satisfying
\[
\Lambda \ll \mu_2 \ll \mu_1 \ll 1/r,
\]
one writes
\[
V_S(r;\mu_1)=V_S^{\rm RF}(r)+C_2(\mu_1)\,r^2+O(r^3),
\]
\[
\delta E_{\rm us}(r;\mu_1,\mu_2)=\delta E_{\rm us}^{\rm RF}(r)-C_2(\mu_1)\,r^2+O(\mu_2^4 r^3),
\]
so that the \(O(\mu_1^3r^2)\) ambiguities cancel in \(V_S+\delta E_{\rm us}\) [1712.05435].

At still longer wavelengths, the Wilson-line correlator in \(\delta E_{\rm us}\) admits a local expansion whose leading term is the local gluon condensate. In the same large-\(\beta_0\) framework, the condensate contribution carries an \(O(\mu_2^4 r^2/\Lambda^2)\) ambiguity from the \(u=2\) renormalon, but this dependence cancels against the corresponding term in \(\delta E_{\rm us}(r;\mu_1,\mu_2)\). The resulting renormalon-free condensate is
\[
\langle g^2G^2\rangle_{\rm RF}
\equiv
\langle g^2G^2\rangle_{k<\mu_2}
-
\frac{12N_c}{T_F}\cdot \frac{C_F}{8\pi}
\int_0^{\mu_2^2} d\tau\,\tau\,\alpha_{\rm eff}(\tau),
\]
and the final short-distance expansion reads
\[
V_{\rm QCD}(r)
=
V_S^{\rm RF}(r)
+
\delta E_{\rm us}^{\rm RF}(r)
+
\frac{T_F}{12N_c}\frac{r^2}{\Lambda^2}\langle g^2G^2\rangle_{\rm RF}
+
O(r^4\Lambda^4).
\]

This factorization is designed for very short distances \(\Lambda_{\rm QCD}r\lesssim0.1\). In that regime, \(\delta E_{\rm us}^{\rm RF}(r)\) is negative, of order \(\alpha_s^3 r^2\), decreases as \(r\to 0\), and at \(r=0.05\,\Lambda^{-1}\) is only a few MeV for \(\Lambda\sim 300\,{\rm MeV}\); its effect on the force is at the \(0.1\%\) level [1712.05435].

## 5. Infrared-subtracted condensates in \(\tau\) observables

The renormalon-free gluon-condensate scheme acquired particular prominence in analyses of the Adler function and hadronic \(\tau\) decay, where the discrepancy between contour-improved perturbation theory (CIPT) and fixed-order perturbation theory (FOPT) was traced to infrared sensitivity associated with the gluon-condensate renormalon. In the large-\(\beta_0\) test case, the Adler-function Borel transform has a simple pole at \(u=2\), and one isolates the pure gluon-condensate piece
\[
B_{G^2}(u)=\frac{8}{\beta_0}\,\frac{e^{5u/3}}{2-u}.
\]
Subtracting this pole defines a subtracted perturbative Adler function and, at the same time, a renormalon-free condensate
\[
\langle G^2\rangle_{\rm RF}
=
\langle G^2\rangle_{\overline{\rm MS}}
+
\delta\langle G^2\rangle(s_0).
\]
For the kinematic \(\tau\) weight \(W_\tau(x)=(1-x)^3(1+x)\), one has \(J_{2,W_\tau}=0\), so the explicit condensate term vanishes, yet the subtraction removes the bulk of the FOPT–CIPT discrepancy: a standard-scheme difference of about \(+0.010\) around orders \(n=5\)–\(7\) is reduced to \(\lesssim 2\times 10^{-3}\) in the renormalon-free scheme [2111.09614].

The full-QCD \(R\)-scheme introduced an infrared factorization scale \(R\) and defined the order-dependent standard condensate by
\[
\langle \bar G^2\rangle^{(n)}
=
\langle G^2\rangle(R^2)
-
R^4\sum_{\ell=1}^n N_g\,r_\ell^{(4,0)}\,\bar a_R^\ell,
\]
with
\[
r_\ell^{(4,0)}
=
\Bigl(\frac12\Bigr)^{\ell+4\hat b_1}
\frac{\Gamma(\ell+4\hat b_1)}{\Gamma(1+4\hat b_1)}.
\]
This subtraction removes the \(u=2\) renormalon from the perturbative Adler series and repackages it into the scale-dependent condensate. In this framework, FOPT contributes identically zero to gluon-condensate-suppressed moments, while the CIPT subtraction is nonzero and cancels the piece responsible for the earlier asymptotic separation. In the large-\(\beta_0\) toy model, in large-\(\beta_0\), and in the 5-loop plus renormalon model for full QCD, FOPT and CIPT in the renormalon-free scheme approach the same “true” Borel sum, and moments that were poorly convergent in the \(\overline{\rm MS}\) scheme become well behaved as well [2202.10957].

The phenomenological implementation depends on the renormalon norm \(N_g\). A follow-up analysis used three methods to determine \(N_g\) for \(n_f=3\)—a multi-renormalon Borel model, a conformal-mapping method, and an optimal-subtraction \(\chi^2\) method—and adopted
\[
N_g=0.57\pm0.23.
\]
Applying the scheme to state-of-the-art \(\alpha_s(m_\tau^2)\) determinations, the truncated-OPE analysis changed from \(0.3169\pm0.0065\) versus \(0.3366\pm0.0070\) in the standard scheme to \(0.3169\pm0.0068\) versus \(0.3197\pm0.0100\) in the renormalon-free scheme, with combined average
\[
\alpha_s(m_\tau^2)=0.3183\pm0.0069.
\]
In the DV-model analysis, the standard-scheme values \(0.3083\pm0.0072\) and \(0.3256\pm0.0097\) became \(0.3081\pm0.0072\) and \(0.3159\pm0.0096\), with combined average
\[
\alpha_s(m_\tau^2)=0.3120\pm0.0082.
\]
In both cases, the FOPT result was essentially unchanged and the CIPT result moved toward FOPT [2207.01116]. A related presentation using a realistic high-order Borel model quoted \(N_g\equiv N_2\approx 0.64\) and emphasized that, for the kinematic and \(x^2\)-suppressed moments, the renormalon-free scheme yields compatible FOPT and CIPT extractions near the input \(\alpha_s(m_\tau^2)=0.315\), while moments with unsuppressed gluon-condensate contributions also become more stable [2207.06081].

## 6. Gradient flow as a gauge-invariant renormalon subtraction

The gradient-flow approach replaces the auxiliary infrared cutoff by the flow time \(t\), which acts as a gauge-invariant factorization scale. One introduces flowed fields \(B_\mu(t,x)\) satisfying
\[
\partial_t B_\mu(t,x)=D_\nu G_{\nu\mu}(t,x),
\qquad
B_\mu(0,x)=A_\mu(x),
\]
and defines the flowed action density
\[
E(t)\equiv \Big\langle \frac{g_s^2}{4}\,G^a_{\mu\nu}(t,x)\,G^{a,\mu\nu}(t,x)\Big\rangle.
\]
Its small-\(t\) OPE is
\[
E(t)=\pi^2\Bigl[\frac{\tilde C_1(t)}{t^2}
+\tilde C_{G^2}(t)\,\langle (\alpha_s/\pi)G^2\rangle
+O(t)\Bigr],
\]
or, in the notation of the later paper,
\[
\frac{E(t)}{\pi^2}=\frac{Y_1(t)}{t^2}+Y_{FF}(t)\,\Big\langle \frac{\alpha_s}{\pi}G^2\Big\rangle+O(t^2).
\]
Solving this relation for the condensate and substituting it into the Adler-function OPE produces a subtracted perturbative coefficient,
\[
D_{\rm ren.free}(Q^2)=C_1(Q^2)-\frac{r}{t^2Q^4}\,\tilde C_1(t),
\qquad
r(Q^2,t)=\frac{C_{G^2}(Q^2)}{\tilde C_{G^2}(t)},
\]
so that the same infrared renormalon singularity in \(C_1\) and in \(\tilde C_1\) cancels exactly, order by order in \(\alpha_s\), without introducing or determining the non-perturbative Stokes constant that would normalize the renormalon [2309.10853].

Applied to hadronic \(\tau\) decay, this construction modifies only the contour-improved series in practice. Using \(\mu=m_\tau\), \(\alpha_s(m_\tau)=0.34\), and \(t\) chosen so \(1/\sqrt{t}\approx 1.5\,{\rm GeV}\), the quoted partial sums through \(O(\alpha_s^4)\) are
\[
\delta_{\rm FO}^{(0)}\approx 0.2200(+0.0088),\qquad
\delta_{\rm CI}^{(0)}\approx 0.1984(+0.0038),\qquad
\delta_{\rm CI,RS}^{(0)}\approx 0.2207(+0.0053),
\]
so the subtraction leaves FO essentially unchanged but lifts CI into agreement with FO already by the fourth order [2309.10853].

The later gradient-flowed OPE formulation used lattice data on the action density to estimate the flowed condensate and reported a marked improvement in convergence. At \(Q=m_\tau\) with \(\alpha_s(m_\tau)=0.314\), the unsubtracted perturbative coefficients for \(\Delta_D=\sum d_n a_Q^n\) were
\[
d_n:\ 1.00,\ 1.64,\ 6.37,\ 49.08,\ldots,
\]
while after subtraction, with \(t=20/(8Q^2)\),
\[
\tilde C_n:\ 0.92,\ 0.94,\ 0.69,\ 4.63,\ldots.
\]
The partial sums
\[
S_N=\sum_{n=1}^N \tilde C_n a_Q^n
\]
stabilized already at \(N\ge 2\) to \(S_N\approx 0.136(1)\). Repeating the analysis for \(1.2\,{\rm GeV}\le Q\le 3\,{\rm GeV}\), the paper reported significantly reduced theoretical uncertainty and extended reliable predictions down to \(Q\approx 1.2\,{\rm GeV}\); for \(\tau\) decay, the FO and CI prescriptions converged to the same result, in favour of the fixed-order treatment [2510.12193].

## 7. Shared structural features, limitations, and scope

Across the existing literature, a renormalon-free gluon-condensate scheme has three recurring elements. First, the perturbative contribution is split into a renormalon-free part and a subtraction term tied to the \(u=2\) singularity. Second, the same subtraction is inserted into the condensate definition, producing a nonperturbative quantity that is independent of the auxiliary prescription variable. Third, the OPE is rewritten so that the perturbative coefficient and the condensate matrix element are separately free of the leading \(O(\Lambda^4)\) ambiguity. This structure is explicit in the \(\mu_f\)-dependent large-\(\beta_0\) OPE [1807.10064], in the PV and hyperasymptotic plaquette formulation [2009.01285], in the \(R\)-dependent \(\tau\) scheme [2202.10957], and in the gradient-flowed Adler-function expansion [2510.12193].

The practical limitations differ by implementation. In the \(R\)-scheme for \(\tau\) observables, the subtraction depends on the renormalon norm \(N_g\), which has to be supplemented independently and was assigned a \(40\%\) uncertainty in the dedicated determination [2207.01116]. In the plaquette hyperasymptotic analysis, the error budget included a \(\pm 40\%\) error in the renormalon normalization \(Z_P\) [2009.01285]. In the original large-\(\beta_0\) proposal based on gradient-flow energy-density data, the difference
\[
\Delta(t)=\hat E_{\rm latt}(t)-c_{1,\hat E}^{\rm RF}(t)
\]
did not yet display a clean \(t^2\) behavior in the available \(t\)-window, and higher-power or higher-log terms reflecting truncation errors spoiled a stable extraction; the authors therefore emphasized the need to go beyond the large-\(\beta_0\) approximation and beyond existing three-loop information [1807.10064].

These results suggest that “renormalon-free gluon-condensate scheme” is best understood not as a single unique prescription but as a family of OPE rearrangements adapted to different observables and regulators. What is common is the demand that the perturbative subtraction and the condensate definition be matched so that the leading renormalon ambiguity is canceled explicitly. Within that shared framework, the lattice plaquette, the static potential, \(\tau\) moments, and gradient-flowed operators provide complementary realizations of the same underlying principle.

Source: https://www.emergentmind.com/topics/renormalon-free-gluon-condensate-scheme