Skip to content

Scheme Transformations and RG Invariants

A renormalization scheme is a choice of finite coordinates for the same renormalized physics. When two mass-independent schemes are related by an analytic, locally invertible redefinition of couplings and fields, their beta functions are different component descriptions of one RG vector field. Coupling values, field anomalous dimensions away from fixed points, higher beta-function coefficients, and the conventional normalization of Λ\Lambda can change; consistently translated observables cannot.

This page derives those transformation laws and works a one-coupling map through three beta-function coefficients. It also marks the assumptions that are often suppressed: the map must be regular on the domain, all quantities must be re-expanded to the same order, and mass-dependent prescriptions require explicit threshold variables. These qualifications are essential before calling a fixed point, exponent, or transmuted scale “scheme independent.”

Required background. Beta Functions, Running Masses, and Field Anomalous Dimensions fixes the signs of β\beta, γ\gamma, and mass running. Renormalization Conditions, Schemes, and Finite Parts explains how finite counterterms define a prescription.

Helpful background. Operator Anomalous-Dimension Matrices supplies the operator-basis conventions used below.

Let gig^i be dimensionless renormalized couplings in one scheme and let

gi=fi(g),Jij(g)figj.g'^i=f^i(g), \qquad J^i{}_{j}(g) \equiv \frac{\partial f^i}{\partial g^j}.

Assume that ff is analytic near the perturbative point, preserves the declared coupling normalization, has no explicit μ\mu dependence, and satisfies

detJ(g)0\det J(g)\ne0

throughout the comparison domain. Differentiation at fixed bare data gives

βi(g)=Jij(g)βj(g).\boxed{ \beta'^i(g') = J^i{}_{j}(g)\,\beta^j(g). }

Thus beta functions transform as a vector field. If g(t)g(t) solves dgi/dt=βi(g)d g^i/dt=\beta^i(g), then

g(t)=f(g(t))g'(t)=f\bigl(g(t)\bigr)

solves the primed flow. The numerical coordinates differ, but the two curves are related point by point.

Finite local counterterms generate maps of this kind in perturbation theory. Their composition, identity, and local inverse form the physical finite-renormalization slice of the more general Stückelberg–Petermann description of renormalization freedom. The theorem-level construction is not needed here; the present application uses only the analytic map of a finite set of renormalized parameters and fields. Collins derives the chain-rule transformation directly from finite changes of prescription in Collins 1984/2023, § 7.8, pp. 200–202.

The conditions are substantive, not cosmetic:

ConditionWhat it protectsFailure mode
f(g)=g+O(g2)f(g)=g+O(g^2) in the chosen coupling variablesthe same weak-coupling normalizationa rescaled leading coefficient is mistaken for scheme dependence
analytic expansion near the reference pointorder-by-order perturbative translationfractional powers or essential singularities destroy the loop expansion
detJ0\det J\ne0 on the domaina one-to-one coordinate patchdistinct physical points are collapsed or a false fixed point is created
consistent field and parameter translationequality of complete predictionsonly the running coupling is changed while amplitudes are left unconverted
common active degrees of freedomcomparison of the same theorya threshold change is mislabeled as a finite scheme change

For example, g=g3g'=g^3 is not an admissible coordinate near the Gaussian point: dg/dg=0dg'/dg=0 at g=0g=0, and its inverse is nonanalytic. The formula β=Jβ\beta'=J\beta still holds algebraically, but conclusions that require a regular perturbative coordinate map do not.

One-coupling coefficients through three orders

Section titled “One-coupling coefficients through three orders”

Take the asymptotically free convention

β(g)=b0g3b1g5b2g7+O(g9)\beta(g) = -b_0g^3-b_1g^5-b_2g^7+O(g^9)

and an analytic finite redefinition tangent to the identity,

g=g+a1g3+a2g5+O(g7).g' = g+a_1g^3+a_2g^5+O(g^7).

Its inverse is

g=ga1g3+(3a12a2)g5+O(g7).g = g'-a_1g'^3 +(3a_1^2-a_2)g'^5 +O(g'^7).

The chain rule first gives

β(g)=(1+3a1g2+5a2g4+O(g6))β(g).\beta'(g') = \left( 1+3a_1g^2+5a_2g^4+O(g^6) \right)\beta(g).

Re-expanding the right-hand side in gg' yields

β(g)=b0g3b1g5b2g7+O(g9),\beta'(g') = -b'_0g'^3-b'_1g'^5-b'_2g'^7 +O(g'^9),

with

b0=b0,b1=b1,b2=b22a1b1+(2a23a12)b0.\boxed{ \begin{aligned} b'_0&=b_0, \\ b'_1&=b_1, \\ b'_2&= b_2-2a_1b_1 +(2a_2-3a_1^2)b_0. \end{aligned} }

The cancellation in b0b'_0 and b1b'_1 is the universal result for a single coupling whose beta function starts at g3g^3, compared between mass-independent analytic schemes with the same leading normalization. The next coefficient is a coordinate choice. Collins verifies the first-two-coefficient statement and its assumptions in Collins 1984/2023, § 7.8, p. 202.

As a concrete transformation, set a1=1/2a_1=1/2 and a2=0a_2=0. Then

b2=b2b134b0.b'_2 = b_2-b_1-\frac34b_0.

Nothing physical has acquired the extra term b13b0/4-b_1-3b_0/4. It is the third coefficient of the same vector field in the gg' coordinate.

For b0>0b_0>0, the leading transmuted scale in the two coordinates is

Λ=μexp[12b0g2](b0g2)b1/(2b02)[1+O(g2)],\Lambda = \mu\exp\left[-\frac{1}{2b_0g^2}\right] \bigl(b_0g^2\bigr)^{-b_1/(2b_0^2)} \left[1+O(g^2)\right],

with the same two-loop convention as on the preceding page. Since

1g2=1g22a1+O(g2),\frac{1}{g'^2} = \frac{1}{g^2}-2a_1+O(g^2),

the consistently normalized scales satisfy

ΛΛ=exp(a1b0).\boxed{ \frac{\Lambda'}{\Lambda} = \exp\left(\frac{a_1}{b_0}\right). }

For the example a1=1/2a_1=1/2, Λ/Λ=e1/(2b0)\Lambda'/\Lambda=e^{1/(2b_0)}. The scale is constant along either RG trajectory but changes by a fixed conversion factor between schemes. Collins derives this relation and explains what a measurement of Λ\Lambda means in Collins 1984/2023, § 7.9, pp. 205–206.

Let the renormalized field coordinate change by a finite factor,

Φr=Cr(g)Φr.\Phi'_r=C_r(g)\Phi_r.

With the convention dΦr/dt=γrΦrd\Phi_r/dt=-\gamma_r\Phi_r, differentiation gives

γr(g)=γr(g)βi(g)ilnCr(g).\boxed{ \gamma'_r(g') = \gamma_r(g) -\beta^i(g)\partial_i\ln C_r(g). }

If β=O(g3)\beta=O(g^3) and Cr=1+crg2+O(g4)C_r=1+c_rg^2+O(g^4), the leading O(g2)O(g^2) field anomalous dimension is unchanged, while higher coefficients can move. An elementary-field anomalous dimension can also depend on gauge fixing, so it should not be promoted to a physical observable merely because its first coefficient is stable under this narrower class of transformations.

For a multiplicatively running mass mm with

ηm1mdmdt,m=D(g)m,\eta_m \equiv \frac{1}{m}\frac{dm}{dt}, \qquad m'=D(g)m,

the sign is instead

ηm=ηm+βiilnD.\eta'_m = \eta_m +\beta^i\partial_i\ln D.

The difference follows because mm itself, rather than its inverse field normalization, is being rescaled.

For a column of renormalized operators satisfying

dOdt=γO,\frac{d\mathbf O}{dt} = -\boldsymbol\gamma\,\mathbf O,

let O=R(g)O\mathbf O'=R(g)\mathbf O. Then

γ=RγR1(βiiR)R1.\boxed{ \boldsymbol\gamma' = R\boldsymbol\gamma R^{-1} -\left(\beta^i\partial_iR\right)R^{-1}. }

Wilson coefficients transform as C=RTC\mathbf C'=R^{-T}\mathbf C, so the contraction CTO\mathbf C^T\mathbf O is unchanged. This is the matrix version of the same principle: neither an operator component nor a coefficient component is separately invariant under a finite basis change.

If gg_\star is an exact fixed point and ff is regular there, then

β(g)=0β(f(g))=0.\beta(g_\star)=0 \quad\Longrightarrow\quad \beta'\bigl(f(g_\star)\bigr)=0.

The coordinate location moves from gg_\star to g=f(g)g'_\star=f(g_\star), but the zero survives. In several couplings, define the stability matrix

Bijβigjg.B^i{}_{j} \equiv \left. \frac{\partial\beta^i}{\partial g^j} \right|_{g_\star}.

At the fixed point, terms involving derivatives of JJ multiply β(g)\beta(g_\star) and vanish. Therefore

B=JBJ1.\boxed{ B'=J_\star B J_\star^{-1}. }

The matrix entries and eigenvectors are coordinate dependent; the eigenvalues are invariant under a regular map. For operators, the derivative term in the transformed anomalous-dimension matrix also vanishes at β=0\beta=0, leaving a similarity transformation. Scaling dimensions in a closed physical operator sector are consequently invariant.

Two cautions prevent overstatement:

  1. A zero of a truncated beta function is not an exact zero. Re-expansion after a scheme change can shift, create, or remove an apparent strong-coupling root by terms of the first omitted order.
  2. A singular map does not preserve the fixed-point argument. If JJ_\star is noninvertible, the similarity relation does not exist.

A reproducible calculation uses these as adversarial checks: a regular analytic map must preserve invariant data to the declared tolerance, while a singular map or a trajectory outside the perturbative box must fail explicitly.

RG-invariant combinations and effective charges

Section titled “RG-invariant combinations and effective charges”

An invariant I(μ,g)I(\mu,g) satisfies

(μμ+βigi)I=0.\left( \mu\frac{\partial}{\partial\mu} +\beta^i\frac{\partial}{\partial g^i} \right)I=0.

The transmuted Λ\Lambda parameter is one example. For a multiplicatively running mass,

m^=m(μ)exp[g(μ)ηm(g)β(g)dg]\widehat m = m(\mu) \exp\left[ -\int^{g(\mu)} \frac{\eta_m(g')}{\beta(g')}\,dg' \right]

is constant along a one-coupling trajectory. Its conventional normalization can still change under m=D(g)mm'=D(g)m, just as Λ\Lambda changes under a coupling redefinition. A physical mass expressed in the same inputs does not.

Another useful coordinate is an effective charge defined by an observable. Let a=g2a=g^2 and normalize an infrared-safe dimensionless observable so that

R(Q)=a(μ)+r1(Q/μ)a2(μ)+.R(Q)=a(\mu)+r_1(Q/\mu)a^2(\mu)+\cdots.

Defining aR(Q)R(Q)a_R(Q)\equiv R(Q) makes the observable itself a coupling coordinate wherever dR/da0dR/da\ne0. Its beta function is

βR(aR)=dRdaβa(a),\beta_R(a_R) = \frac{dR}{da}\,\beta_a(a),

which is again the vector-field chain rule. This removes an arbitrary intermediate coupling from the final statement, but aRa_R remains observable specific and its finite-order extraction retains truncation and power-correction uncertainties. The effective-charge construction was introduced for perturbative QFT observables in Grunberg 1980, pp. 70–74.

Mass-dependent schemes and finite-order comparisons

Section titled “Mass-dependent schemes and finite-order comparisons”

The preceding coefficient claims require mass independence. If a scheme map depends on

ρmμ,g=f(g,ρ),\rho\equiv\frac{m}{\mu}, \qquad g'=f(g,\rho),

then

dρdt=(ηm1)ρ\frac{d\rho}{dt} = (\eta_m-1)\rho

and the transformed beta function is

β(g,ρ)=β(g,ρ)fg+(ηm1)ρfρ.\boxed{ \beta'(g',\rho) = \beta(g,\rho)\frac{\partial f}{\partial g} +(\eta_m-1)\rho \frac{\partial f}{\partial\rho}. }

For ρ=m2/μ2\rho=m^2/\mu^2, replace ηm1\eta_m-1 by ηm22\eta_{m^2}-2. The extra term is physical bookkeeping for the subtraction point moving relative to a mass threshold. One cannot compare its coefficients with a mass-independent beta function by applying the b0=b0b'_0=b_0, b1=b1b'_1=b_1 rule globally.

At finite perturbative order, translate in this order:

  1. express the primed couplings, masses, fields, operators, and matching inputs in terms of the unprimed ones;
  2. substitute into the complete truncated prediction;
  3. re-expand consistently through the declared order;
  4. compare only after using the same physical boundary data.

Keeping selected unexpanded pieces of the finite map inserts a subset of higher-order terms. That may define a resummation prescription, but it is not an order-by-order proof of scheme independence. The difference between consistently translated truncations begins at the first omitted order; its numerical size is a diagnostic, not a universal confidence level.

Coordinate-dependent data and invariant claims

Section titled “Coordinate-dependent data and invariant claims”

The following comparison table is reused across this chapter. Each row separates what can move under a finite convention change from the claim that survives a complete translation.

ItemWhat may changeWhat survives a consistent translationRequired qualification or check
Renormalized gig^i, masses, and field normalizationsNumerical values under finite scheme or basis changesA prediction expressed in the same physical inputsTranslate every parameter and field factor through the retained order
Beta function away from a fixed pointComponents and higher-order coefficientsThe integral curves as geometric trajectories under a nonsingular coordinate mapCompare transformed vector fields, not coefficients at equal numerical coupling
Elementary-field anomalous dimensionFinite field rescaling; gauge parameter in a gauge theoryScaling of a gauge-invariant observable after all factors are combinedNever identify a gauge-dependent elementary-field exponent with an observable
Exact fixed pointCoordinate location gig_\star^iExistence of the zero under a regular mapExclude singular redefinitions and verify the fixed point lies in the method’s domain
Fixed-point stability dataMatrix representation and basisEigenvalues in a closed physical sectorInclude operator mixing and redundant directions before diagonalizing
Transmuted scaleIts conventional normalizationMatched dimensionless ratios or predictionsState the scheme and reference condition defining the scale
Zero or singularity of a truncated beta functionLocation and even apparent existence at insufficient orderOnly the demonstrated breakdown of the stated approximationVary scheme/order and stop before couplings become large
Wilson coefficient versus power correctionFactorization scheme and, for an asymptotic series, summation prescriptionTheir consistently defined sum in an observableMatch the ambiguity of the perturbative term to the operator matrix element; developed on the renormalon page
Residual μ\mu or scheme dependenceNumerical size at finite orderVanishing in the exact consistently matched predictionTreat the residual as a diagnostic, not a universal probability law

The broader classification of finite local ambiguities is developed in Microlocal Renormalization Ambiguities and Their Classification. That theorem-level result supplies a structural setting; the table above remains the practical contract for the perturbative calculations in this volume.

Take the exact benchmark from the preceding page,

dgdt=g3,g(0)=0.4,\frac{dg}{dt}=-g^3, \qquad g(0)=0.4,

and apply the regular map

h=f(g)=g+12g3.h=f(g)=g+\frac12g^3.

Since dh/dg=1+3g2/2>0dh/dg=1+3g^2/2>0, the map is invertible on the entire positive real branch. The exact transformed equation is

dhdt=(1+32g2)g3,g=f1(h),\frac{dh}{dt} = -\left(1+\frac32g^2\right)g^3, \qquad g=f^{-1}(h),

while its perturbative re-expansion is

βh(h)=h3+34h7+O(h9).\beta_h(h) = -h^3+\frac34h^7+O(h^9).

The invariant in the new coordinate is not 1/h22t1/h^2-2t. It is

Ih(t)=1[f1(h(t))]22t=6.25.I_h(t) = \frac{1}{[f^{-1}(h(t))]^2}-2t =6.25.

Exact reference values expose the difference:

ttg(t)g(t)h(t)h(t)Corrected Ih(t)I_h(t)Naive 1/h22t1/h^2-2t
2-20.6666666666670.6666666666670.8148148148150.8148148148156.256.255.506198347115.50619834711
000.4000000000000.4000000000000.4320000000000.4320000000006.256.255.358367626895.35836762689
220.3123475237770.3123475237770.3275839883520.3275839883526.256.255.318685776105.31868577610
550.2480694691780.2480694691780.2557023759220.2557023759226.256.255.294330585885.29433058588

A numerical implementation should map the trajectory point by point and preserve the corrected invariant to relative tolerance 101010^{-10}. Comparing gg and hh at equal numerical values, or reusing the untransformed formula with hh substituted for gg, is precisely the error this test is designed to catch.

Comparing coefficients at equal numerical coupling. Equal numbers g=gg=g' generally label different physical points. Convert the coordinate first, then compare the vector fields or predictions.

Calling all beta-function coefficients universal. For one mass-independent coupling with the stated normalization, b0b_0 and b1b_1 survive; b2b_2 and higher coefficients do not. Multiple couplings and mass-dependent schemes require separate analysis.

Using a singular redefinition to erase a fixed point. Fixed-point preservation assumes an invertible Jacobian. A map that fails this test is not evidence that the original zero was unphysical.

Transforming the running but not the observable. Finite parts of amplitudes, masses, fields, operator bases, and matching conditions must change with the coupling. A partial conversion manufactures scheme dependence.

Treating residual variation as a statistical interval. Scheme and scale scans probe selected higher-order directions. Without an explicit uncertainty model they are diagnostics, not probabilities.

Starting from the general gg' map above, verify the expression for b2b'_2.

Solution

First multiply by the Jacobian:

β(g)=b0g3(b1+3a1b0)g5(b2+3a1b1+5a2b0)g7+O(g9).\begin{aligned} \beta'(g') &= -b_0g^3 -(b_1+3a_1b_0)g^5 \\ &\quad -(b_2+3a_1b_1+5a_2b_0)g^7 +O(g^9). \end{aligned}

Use

g3=g33a1g5+(12a123a2)g7+O(g9),g5=g55a1g7+O(g9),g7=g7+O(g9).\begin{aligned} g^3&=g'^3-3a_1g'^5+(12a_1^2-3a_2)g'^7+O(g'^9), \\ g^5&=g'^5-5a_1g'^7+O(g'^9), \\ g^7&=g'^7+O(g'^9). \end{aligned}

Collecting powers gives b0=b0b'_0=b_0, b1=b1b'_1=b_1, and

b2=b22a1b1+(2a23a12)b0.b'_2=b_2-2a_1b_1+(2a_2-3a_1^2)b_0.

Show that the slope of a one-coupling beta function at an exact fixed point is invariant under a regular redefinition.

Solution

For g=f(g)g'=f(g),

β(g)=f(g)β(g).\beta'(g')=f'(g)\beta(g).

Differentiate with respect to gg':

dβdg=1f(g)ddg[f(g)β(g)]=f(g)f(g)β(g)+dβdg.\frac{d\beta'}{dg'} = \frac{1}{f'(g)} \frac{d}{dg}\left[f'(g)\beta(g)\right] = \frac{f''(g)}{f'(g)}\beta(g) +\frac{d\beta}{dg}.

At g=gg=g_\star, the first term vanishes because β(g)=0\beta(g_\star)=0. Hence dβ/dgg=dβ/dggd\beta'/dg'\rvert_{g'_\star}=d\beta/dg\rvert_{g_\star}, provided f(g)0f'(g_\star)\ne0.

  • Collins, John C. Renormalization: An Introduction to Renormalization, the Renormalization Group, and the Operator-Product Expansion. Cambridge Monographs on Mathematical Physics. Cambridge: Cambridge University Press, 1984; open-access digital edition, 2023. DOI. Open PDF.
  • Grunberg, Georges. “Renormalization Group Improved Perturbative QCD.” Physics Letters B 95 (1980): 70–74; erratum 110 (1982): 501. DOI.