Principle of Maximum Conformality (PMC) in QCD

• PMC is a renormalization-scale setting method that decomposes perturbative series into conformal and non-conformal parts by absorbing β-function terms into the running coupling.
• The technique reorganizes fixed-order predictions to eliminate scale ambiguity and improve convergence by matching perturbative coefficients with renormalization-group invariance.
• PMC enhances scheme independence and practical stability in QCD calculations, with applications spanning heavy-quarkonium production, event shapes, and other observables.

Searching arXiv for recent and foundational PMC papers to ground the article. The Principle of Maximum Conformality (PMC) is a renormalization-scale setting method for perturbative quantum chromodynamics (pQCD) in which all non-conformal terms associated with the QCD $\beta$-function are identified order by order and absorbed into the running coupling, leaving a residual perturbative series identical in structure to that of a conformal theory with $\beta=0$ [(Brodsky et al., 2011); (Yan et al., 2023)]. In the PMC formulation, the renormalization scale is not treated as a guessed kinematic parameter but as a process- and, in many applications, kinematics-dependent quantity fixed by the renormalization-group structure of the perturbative coefficients themselves [(Brodsky et al., 2011); (Wang et al., 2019)]. The method is presented in the literature as the all-orders extension of the Brodsky–Lepage–Mackenzie procedure, as consistent with standard renormalization-group invariance, and as reducing to the Gell-Mann–Low prescription in the Abelian limit [(Brodsky et al., 2011); (Brodsky et al., 2012); (Yan et al., 2023)].

1. Conceptual definition and renormalization-group basis

PMC is built on the premise that the running coupling exists to resum all perturbative contributions generated by renormalization, namely the terms controlled by the QCD $\beta$-function coefficients $\{\beta_i\}$ [(Brodsky et al., 2011); (Yan et al., 2023)]. In this framework, a fixed-order pQCD prediction is first decomposed into conformal pieces and non-conformal pieces. The non-conformal pieces are those proportional to $\beta_0,\beta_1,\dots$, often visible at low orders through their $n_f$-dependence, provided that only the ultraviolet-renormalization-associated $n_f$-terms are identified as $\beta$-terms (Yan et al., 2023, Wang et al., 2018). Once those terms are absorbed into the arguments of the running couplings, the remaining coefficients are the conformal coefficients, meaning the coefficients of the corresponding $\beta=0$ theory [(Brodsky et al., 2011); (Giustino et al., 2020)].

This definition is tied directly to standard renormalization-group invariance. The PMC literature repeatedly states that physical observables should be independent of the initial renormalization scale and the renormalization scheme, and that conventional fixed-order truncation violates this requirement in practice because the scale dependence of the coupling and the coefficients is not matched order by order [(Brodsky et al., 2012); (Huang et al., 2024)]. PMC is proposed as the remedy: by resumming all known renormalization-group-controlled $\beta$-terms into the coupling, it produces a conformal series whose finite-order prediction is independent of the initial choice of renormalization scale up to unknown higher-order terms [(Brodsky et al., 2012); (Yan et al., 2023)].

A standard low-order illustration writes an observable as a perturbative series such as

$\beta=0$0

and then uses the scale-displacement relation

$\beta=0$1

to choose

$\beta=0$2

so that the $\beta=0$3-dependent term is absorbed into the coupling and the NLO coefficient becomes conformal (Yan et al., 2023). At higher orders, the same logic is applied recursively to the full $\beta=0$4-pattern (Zeng et al., 2018, Salinas-Arizmendi et al., 2022).

The method is also framed as a non-Abelian generalization of QED scale setting. In the Abelian limit, PMC is claimed to reproduce the Gell-Mann–Low prescription, where the photon virtuality sets the physical scale of the running electromagnetic coupling [(Brodsky et al., 2011); (Brodsky et al., 2012); (Yan et al., 2023)]. This Abelian-limit requirement is treated in the PMC literature as a decisive theoretical consistency condition.

2. Perturbative reorganization and scale determination

A recurrent formal pattern in PMC applications is the rewriting of a fixed-order series into conformal and non-conformal components. For example, in the NLO analysis of $\beta=0$5, the total cross section is written as

$\beta=0$6

with

$\beta=0$7

and the NLO coefficient decomposed as

$\beta=0$8

This is then reorganized into

$\beta=0$9

where $\beta$0 and $\beta$1 are conformal and $\beta$2 is the non-conformal term to be absorbed into the coupling (Huang et al., 2024). The corresponding PMC scale at leading-logarithmic accuracy is

$\beta$3

and the PMC-improved series becomes

$\beta$4

(Huang et al., 2024).

The same structure appears in event shapes. For the thrust distribution, the NLO coefficient is decomposed as

$\beta$5

and the PMC scale is determined by

$\beta$6

leading to the conformal form

$\beta$7

(Wang et al., 2019).

For higher-order applications, the $\beta$8-pattern becomes more elaborate. In the $\beta$9 analysis of $\{\beta_i\}$0, the width is reorganized as

$\{\beta_i\}$1

and then transformed to the multi-scale conformal series

$\{\beta_i\}$2

(Zeng et al., 2018). In this multi-scale realization, each perturbative order receives its own effective scale $\{\beta_i\}$3, reflecting the possibility that the relevant virtuality differs from order to order (Zeng et al., 2018).

A major methodological variant is the single-scale PMC. In this implementation, one effective scale $\{\beta_i\}$4 is chosen so that all known non-conformal terms vanish simultaneously. This version is emphasized in several applications because it suppresses residual scale dependence associated with the highest unknown-order scale and simplifies practical implementation (Yan et al., 2023, Yu et al., 2019, Huang et al., 2019).

3. Variants, extensions, and relation to other frameworks

PMC is presented in the literature as the all-orders extension of BLM, with the “PMC–BLM correspondence principle” asserting their equivalence at the level of underlying scale-setting logic [(Brodsky et al., 2012); (Brodsky et al., 2011)]. What PMC adds is a systematic all-orders language in terms of explicit $\{\beta_i\}$5-pattern identification and conformal/non-conformal separation [(Brodsky et al., 2012); (Yan et al., 2023)].

The formalism has been developed in both multi-scale and single-scale forms. The multi-scale form assigns distinct scales to different perturbative orders or subprocess structures, as in $\{\beta_i\}$6 (Zeng et al., 2018), top-pair production (Brodsky et al., 2012), and polarized double-charmonium production (Wang et al., 2013). The single-scale form is used when one seeks one effective process scale, for example in heavy-quarkonium inclusive decays (Yu et al., 2019), $\{\beta_i\}$7 (Huang et al., 2019), and the NLO analysis of $\{\beta_i\}$8 (Huang et al., 2024).

A more ambitious extension is “PMC$\{\beta_i\}$9,” which is based on the claimed existence of intrinsic Conformality ($\beta_0,\beta_1,\dots$0). In this construction, perturbative observables are decomposed into disjoint “conformal subsets,” each governed by one conformal coefficient and one intrinsic scale (Giustino et al., 2020). The observable is first written in the form

$\beta_0,\beta_1,\dots$1

$\beta_0,\beta_1,\dots$2

$\beta_0,\beta_1,\dots$3

so that each subset is scale invariant by itself (Giustino et al., 2020). This formulation claims that the all-orders scale for a subset is fixed once the lowest-order $\beta_0,\beta_1,\dots$4-logarithm is known (Giustino et al., 2020). A plausible implication is that PMC$\beta_0,\beta_1,\dots$5 is intended not only as a scale-setting prescription but also as a structural reorganization of perturbation theory.

The literature also places PMC in the context of commensurate scale relations (CSRs), where relations between effective charges of different observables are written with process-dependent relative scales chosen so that all scheme-dependent $\beta_0,\beta_1,\dots$6-terms are absorbed [(Brodsky et al., 2011); (Yan et al., 2023)]. The generalized Crewther relation is cited as a canonical example (Yan et al., 2023).

4. Phenomenological applications across QCD processes

PMC has been applied to a wide range of observables, and the applications consistently use the same core logic: identify $\beta_0,\beta_1,\dots$7-controlled terms, absorb them into the coupling, and compare the resulting conformal series with conventional fixed-order predictions.

In top-pair hadroproduction at NNLO, the partonic cross sections are decomposed into channel-dependent structures with distinct non-conformal terms, including separate Coulombic and non-Coulombic contributions (Brodsky et al., 2012). After PMC scale setting, the total $\beta_0,\beta_1,\dots$8 cross section is reported to be nearly independent of the initial choice of renormalization scale even when $\beta_0,\beta_1,\dots$9 is varied from $n_f$0 to $n_f$1 (Brodsky et al., 2012). The same framework was used for the Tevatron forward-backward asymmetry, where the dip behavior of the NLO PMC scale in the dominant $n_f$2 channel enhances the effective coupling and raises the asymmetry prediction (Brodsky et al., 2012).

In event-shape physics, the thrust distribution in $n_f$3 annihilation provides a clean illustration of a kinematics-dependent PMC scale (Wang et al., 2019). The thrust variable,

$n_f$4

is used to show that the effective scale is not constant but rises monotonically with $n_f$5, reflecting the changing virtuality of the QCD subprocess (Wang et al., 2019). The same analysis emphasizes that the PMC scale also determines the correct number of active flavors via the scale dependence of $n_f$6 (Wang et al., 2019).

In heavy-quarkonium production and decay, PMC has been extensively used within NRQCD. For exclusive double charmonium production $n_f$7, the perturbative series is separated into a pure-QCD sector and a QED-interference sector,

$n_f$8

$n_f$9

with PMC scales

$n_f$0

(Sun et al., 2018). In that analysis, both sectors yield the same numerical value,

$n_f$1

interpreted as the typical momentum flow of the process (Sun et al., 2018).

For $n_f$2, the NLO polarized cross sections are written as

$n_f$3

with

$n_f$4

and channel-dependent PMC scales

$n_f$5

(Wang et al., 2013). This application is frequently cited to illustrate that different helicity channels need not share a common effective scale (Wang et al., 2013).

For heavy-quarkonium decays, the ratio

$n_f$6

was analyzed through NNLO using the PMC single-scale method (Yu et al., 2019). The PMC-improved conformal series is written as

$n_f$7

with the PMC scale expanded as

$n_f$8

(Yu et al., 2019). A related N$n_f$9LO application to $\beta$0 yielded a single effective scale $\beta$1, substantially below the conventional guess $\beta$2 GeV (Huang et al., 2019).

In transition form factors, the NNLO NRQCD prediction for $\beta$3 is reorganized as

$\beta$4

with PMC scale

$\beta$5

(Wang et al., 2018). That paper is notable for stressing that not every $\beta$6-term is a $\beta$7-term: the ultraviolet-finite light-by-light term $\beta$8 is treated as conformal and kept out of the scale displacement (Wang et al., 2018).

PMC has also been used in Higgs phenomenology. For $\beta$9, the multi-scale analysis through $\beta=0$0 led to approximate scales

$\beta=0$1

at $\beta=0$2 (Zeng et al., 2018). The same application emphasizes that the new $\beta=0$3 input greatly suppresses the residual scale dependence that remained in the earlier $\beta=0$4 analysis (Zeng et al., 2018).

Thermal QCD has provided another nontrivial testing ground. In the free energy density at high temperature, the perturbative expansion is separated into a hard contribution $\beta=0$5 and a soft contribution $\beta=0$6,

$\beta=0$7

and PMC is applied separately to the hard-scale coefficients and to the EFT parameters $\beta=0$8 and $\beta=0$9 (Bu et al., 2018). This application makes especially explicit that renormalization-scale setting and factorization-scale setting are distinct issues, and that PMC addresses only the former (Bu et al., 2018).

Most recently, the NNLO study of the $\beta$0-wave quarkonium decays $\beta$1 combines PMC for $\beta$2 with explicit LDME evolution for the factorization scale $\beta$3 (Zhou et al., 23 Jun 2026). After the $\beta$4-dependence is canceled through LDME running, the perturbative width

$\beta$5

is rewritten as

$\beta$6

and the PMC scale is set by

$\beta$7

(Zhou et al., 23 Jun 2026). This illustrates the continued use of PMC in current quarkonium precision studies.

5. Claimed advantages and recurring empirical patterns

The most frequently stated advantage of PMC is the elimination of conventional renormalization-scale ambiguity. In many applications, the PMC prediction is reported to be flat with respect to the initial choice of $\beta$8, whereas the conventional prediction varies strongly under standard scale scans [(Brodsky et al., 2012); (Huang et al., 2024); (Huang et al., 2019)]. This is treated as practical evidence that PMC restores the renormalization-group invariance expected of a physical observable.

A second recurring claim is improved perturbative convergence. The mechanism given is that conventional coefficients contain unresummed running-coupling contributions, including factorially divergent renormalon-type terms of the form $\beta$9, whereas PMC absorbs those $\beta=0$00-dependent pieces into the running coupling and leaves a conformal residual series [(Brodsky et al., 2011); (Brodsky et al., 2012); (Yan et al., 2023)]. In applications such as $\beta=0$01, the relative importance of the NLO correction changes from $\beta=0$02 under conventional scale setting at $\beta=0$03 to $\beta=0$04 after PMC (Huang et al., 2024). In $\beta=0$05, the NNLO term becomes much smaller relative to NLO after PMC than under conventional scale setting (Zhou et al., 23 Jun 2026).

A third claimed advantage is renormalization-scheme independence. The argument is that scheme dependence in a renormalizable gauge theory is tied to the $\beta=0$06-function, so once all scheme-dependent non-conformal $\beta=0$07-terms are removed from the coefficients and absorbed into the coupling, the remaining conformal coefficients are scheme independent [(Brodsky et al., 2011); (Yan et al., 2023)]. The literature often connects this point to commensurate scale relations and to the Abelian-limit reduction to Gell-Mann–Low scale setting [(Brodsky et al., 2011); (Yan et al., 2023)].

A fourth repeated theme is that PMC can determine physically meaningful effective scales rather than guessed ones. In thrust, the scale rises monotonically with $\beta=0$08, reflecting event topology (Wang et al., 2019). In $\beta=0$09, the derived scale $\beta=0$10 GeV is interpreted as the typical momentum flow and retrospectively explains why conventional phenomenology tended to prefer $\beta=0$11-$\beta=0$12 GeV (Sun et al., 2018). In the $\beta=0$13 leptonic decay, the extracted $\beta=0$14 GeV implies a softer effective momentum flow than the conventional hard-scale guess (Huang et al., 2019).

The method is also presented as useful for estimating unknown higher-order terms. Several papers combine PMC with Padé approximation approaches, arguing that the conformalized, renormalon-free series provides a more reliable basis for extrapolating to uncalculated orders (Huang et al., 2024, Yu et al., 2019, Zhou et al., 23 Jun 2026). This suggests a broader use of PMC not only in central-value prediction but also in perturbative uncertainty quantification.

6. Limitations, subtleties, and points of dispute

Despite the strong claims made in the PMC literature, the same papers also identify several limitations and technical caveats.

The first is that residual dependence does not vanish completely at finite order. Even after all known $\beta=0$15-terms are absorbed, the highest perturbative term lacks the next-order information needed to determine its own PMC scale exactly (Zeng et al., 2018, Yan et al., 2023). Thus a residual uncertainty remains from unknown higher-order $\beta=0$16-terms and from the perturbative truncation of the PMC scale itself (Zeng et al., 2018, Huang et al., 2019).

The second is that not every $\beta=0$17-dependent term should be treated as non-conformal. Several papers explicitly warn that ultraviolet-finite $\beta=0$18-dependent terms, such as light-by-light contributions, are conformal and must not be absorbed into the running coupling (Wang et al., 2018, Yan et al., 2023). This point is technically important because a naive $\beta=0$19-counting procedure can misidentify the $\beta=0$20-pattern.

The third is that PMC does not address factorization-scale ambiguities. This distinction is emphasized in both the thermal free-energy analysis and the $\beta=0$21 study: factorization-scale dependence is associated with the separation between perturbative and nonperturbative physics and remains even in a conformal theory (Bu et al., 2018, Zhou et al., 23 Jun 2026). In such cases, separate EFT or LDME evolution equations are needed to control $\beta=0$22-dependence.

A fourth limitation is that a conformalized series can still have large conformal coefficients. This is stated explicitly in the $\beta=0$23 paper, where the QCD-interference sector becomes very well behaved after PMC, but the pure-QCD sector still has a sizable NLO correction because of a large conformal coefficient $\beta=0$24 (Sun et al., 2018). The authors stress that PMC fixes the renormalization-scale problem but cannot remove genuinely large conformal higher-order effects (Sun et al., 2018).

A fifth issue is practical. In some multi-scale implementations, the extracted scales can become very small in portions of phase space, pushing the coupling toward the nonperturbative domain. This is discussed explicitly in the analysis of $\beta=0$25 from the Bjorken sum rule, where the multi-scale PMC implemented in the $\beta=0$26 auxiliary scheme was found to be effectively applicable only for $\beta=0$27 GeV, motivating suggestions to use a different intermediate scheme or a single effective PMC scale (Deur et al., 2017). Similarly, in event-shape applications of PMC$\beta=0$28, singular or unstable extracted scales in certain kinematic regions had to be regularized (Giustino et al., 2020).

Finally, PMC remains a subject of conceptual dispute. The paper “The Principle of Maximum Conformality Correctly Resolves the Renormalization-Scheme-Dependence Problem” is explicitly framed as a rebuttal to criticisms that PMC does not solve renormalization-scheme dependence and is merely a disguised variant of the Principle of Minimal Sensitivity (Yan et al., 2023). That paper insists that PMC differs fundamentally from PMS because it does not seek a stationary point of the original truncated series but instead transforms the series into a conformal one by absorbing all known $\beta=0$29-dependent terms into the coupling (Yan et al., 2023). The existence of such rebuttal literature indicates that the broader conceptual status of PMC, especially its claims of exact scheme independence at finite order, remains debated.

In aggregate, the PMC literature presents the method as a renormalization-group-based reorganization of pQCD in which scale setting is derived from the perturbative series rather than guessed, the residual coefficients are conformal, and the resulting predictions are more stable and often phenomenologically improved [(Brodsky et al., 2011); (Brodsky et al., 2012); (Wang et al., 2019); (Huang et al., 2024)]. A plausible implication is that PMC is best understood not as a universal replacement for all perturbative uncertainty analysis but as a specific, strongly RG-motivated prescription for isolating running-coupling effects and assigning them to the coupling rather than to the coefficients.