pith. sign in

arxiv: 2509.10658 · v2 · pith:IHUUBJSZnew · submitted 2025-09-12 · ✦ hep-lat · hep-ph· hep-th· nucl-th

Renormalization Group Approach to Confinement

Pith reviewed 2026-05-21 22:54 UTC · model grok-4.3

classification ✦ hep-lat hep-phhep-thnucl-th
keywords renormalization groupgluon condensaterunning couplingconfinementinfrared fixed pointgradient flowQCDscale invariance
0
0 comments X

The pith

A scale-invariant gluon condensate forces the strong coupling to run as Lambda squared over mu squared, driving an infrared fixed point at infinite coupling.

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

The paper applies renormalization group transformations via gradient flow to evolve the strong coupling from short-distance perturbative behavior to long-distance confinement. It argues that the gluon condensate stays constant across scales due to self-similar properties, which directly produces the functional form alpha_S(mu) approximately equal to Lambda_S squared over mu squared. This running leads to a fixed point where the inverse coupling reaches zero, matching the expectation of infrared slavery. A reader would care because the result isolates the condensate as the single essential ingredient for confinement, a feature QCD shares with many other models, and backs the derivation with numerical checks.

Core claim

Building on the scale invariance of the gluon condensate, the running coupling satisfies alpha_S(mu) approximately Lambda_S squared over mu squared. This evolves toward the infrared fixed point where one over alpha_S equals zero, the signature of infrared slavery. The derivation uses renormalization group transformations applied through the gradient flow, and the only essential input is the presence of the condensate itself, a universal feature of QCD and related theories. The analytical form is supported by numerical simulations.

What carries the argument

The scale-invariant gluon condensate, which fixes the infrared running of the strong coupling through renormalization group flow with gradient flow.

If this is right

  • Confinement follows from the running coupling alone once the condensate is present.
  • The infrared fixed point at infinite coupling emerges automatically from the condensate's scale invariance.
  • The same mechanism applies to any model sharing the gluon condensate as a universal feature.
  • Numerical simulations on the lattice confirm the predicted running without additional assumptions.

Where Pith is reading between the lines

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

  • The method could be applied to other gauge theories that possess a similar scale-invariant condensate to predict their confining behavior.
  • It offers a route to extract the confinement scale directly from the condensate value without introducing extra parameters.
  • If the scale invariance holds, the result might generalize to theories with varying numbers of colors or fermion flavors, providing a testable prediction for how confinement strength changes.

Load-bearing premise

The gluon condensate keeps the same value no matter which scale is used to measure it.

What would settle it

A direct lattice measurement showing that the gluon condensate changes value when the renormalization scale is varied, or that the running coupling at low energies deviates from the one-over-mu-squared form.

Figures

Figures reproduced from arXiv: 2509.10658 by Gerrit Schierholz.

Figure 1
Figure 1. Figure 1: The energy density L 2E(t)/8π 2 in units of p 2χt /π as a function of t/L 2 , where χt = hQ 2 i/V has been taken from [18]. The curves are a fit of the form c/t to the left-hand data (t/L 2 . 0.2), and a linear fit to the right-hand data. fall on a single line. This tells us two things. First, hE(t)i ∝ 1/t for t/L 2 . 0.2, in accord with αGF ≃ 8Λ 2 GFt. At t/L 2 ≈ 0.21, corresponding to √ 8t ≈ 1.25 L, the … view at source ↗
Figure 2
Figure 2. Figure 2: The gluon condensate on the 324 lattice in units of the flow parameter w0 as a function of the dimensionless quantity t/L 2 in the gradient flow scheme. is not new [21], but the reference to a physical scale and the onset of classical behavior is. In [18] the energy density has been broken down according to the topological charge |Q|. We now turn to the gluon condensate, on which our derivations are essent… view at source ↗
Figure 3
Figure 3. Figure 3: The running coupling αGF(µ) on the 324 lattice as a function of t/L 2 . On this lattice µ ≈ 60 MeV at the border value. which is subject to considerable theoretical uncertainty however. See the Appendix. To conclude, in [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The ratio of the gluon condensate GQ=0 for topological charge zero to the full gluon condensate G as a function of t/L 2 . The solid curve is a Gaussian fit. Substituting √ 3G/16 for hE(t)it finally gives the result for G in the MS scheme √ G ≈ 0.21 × 41.5 √ χt = 8.7 √ χt ∴ G = 75.7 χt . (20) Given that χt = (180(4) MeV)4 [25], and our result (18) for G, the expected ratio is G/χt = (535/180)4 = 78, which … view at source ↗
Figure 5
Figure 5. Figure 5: The gluon condensate in units of w0 as a function of t/w 2 0 for three values of β on the 323 × 64 (β = 5.65) and 483 × 96 lattice (β = 5.8 and 5.95). exposed in [39]. In previous work we generated flowed configurations with 2 + 1 quark flavors, in an attempt to compute the gradient flow scales √ t0 and w0 [40], to which I refer here. The calculations were done at the SU(3) flavor symmetric point [41], wit… view at source ↗
Figure 6
Figure 6. Figure 6: The running coupling αGF(µ) as a function of t/w 2 0 for our three values of β. 11 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
read the original abstract

While we have several complementary models of confinement, some of which are phenomenologically appealing, we do not have the ability to calculate analytically even simple aspects of confinement, let alone have a framework to eventually prove confinement. The problem we are facing is to evolve the theory from the perturbative regime to the long distance confining regime. This is generally achieved by renormalization group transformations. With the gradient flow we now have a technique to address the problem from first principles. The primary focus is on the running coupling $\alpha_S(\mu)$, from which confinement can be concluded alone. A central point is that the gluon condensate is scale invariant, which reflects its self-similar behavior across different scales. Building on that, we derive $\alpha_S(\mu) \simeq \Lambda_S^2/\mu^2$, which evolves to the infrared fixed point $1/\alpha_S = 0$ in accordance with infrared slavery. The only important factor appears to be the presence of the gluon condensate, which is a universal feature that QCD shares with many other models. The analytical results are supported by numerical simulations.

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.

Referee Report

2 major / 1 minor

Summary. The manuscript proposes a renormalization group framework for confinement in QCD based on the gradient flow. It identifies the scale-invariant gluon condensate as the key ingredient and derives the running coupling α_S(μ) ≃ Λ_S²/μ², which is stated to flow to an infrared fixed point 1/α_S = 0 consistent with infrared slavery. The analytic result is claimed to be supported by numerical simulations, with the gluon condensate presented as a universal feature shared across models.

Significance. If the central derivation were shown to be consistent with both the perturbative ultraviolet regime and the confining infrared regime, the approach would supply a compact analytic link between them that depends primarily on the gluon condensate. The manuscript does not, however, demonstrate such consistency or provide the required matching conditions.

major comments (2)
  1. Abstract and central paragraph on the gluon condensate: the claimed form α_S(μ) ≃ Λ_S²/μ² produces the beta function β(α) = dα/d ln μ = −2α. This linear dependence is incompatible with the perturbative QCD beta function, which begins at quadratic order, β(α) ∼ −(11N_c/3 − 2N_f/3)α²/(2π) + O(α³), at large μ. The scale-invariance assumption on the condensate directly enforces the power-law running, so the ultraviolet matching required by the paper’s stated goal of evolving from the perturbative to the confining regime is not established.
  2. Abstract: the statement that numerical simulations support the analytic results is given without error analysis, explicit fitting procedures, or data-exclusion criteria, leaving the quantitative evidence for the functional form unverified.
minor comments (1)
  1. The manuscript would benefit from an explicit step-by-step derivation of the running coupling from the gradient-flow equation and the scale-invariance assumption.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the insightful comments. We provide point-by-point responses to the major comments below.

read point-by-point responses
  1. Referee: Abstract and central paragraph on the gluon condensate: the claimed form alpha_S(mu) simeq Lambda_S^2/mu^2 produces the beta function beta(alpha) = d alpha / d ln mu = -2 alpha. This linear dependence is incompatible with the perturbative QCD beta function, which begins at quadratic order, beta(alpha) sim -(11 N_c/3 - 2 N_f/3) alpha^2/(2 pi) + O(alpha^3), at large mu. The scale-invariance assumption on the condensate directly enforces the power-law running, so the ultraviolet matching required by the paper's stated goal of evolving from the perturbative to the confining regime is not established.

    Authors: The referee correctly notes that the derived form alpha_S(mu) simeq Lambda_S^2/mu^2 implies a beta function beta(alpha) = -2 alpha. This linear dependence follows directly from the assumption of scale invariance of the gluon condensate within the gradient flow approach. The manuscript centers on the infrared regime and the resulting fixed point consistent with confinement. We acknowledge that the current derivation does not include explicit matching conditions to the perturbative ultraviolet beta function. In the revised version we will add a clarifying paragraph on the regime of validity of the scale-invariance assumption and note that a full ultraviolet-infrared connection remains a topic for future work. revision: yes

  2. Referee: Abstract: the statement that numerical simulations support the analytic results is given without error analysis, explicit fitting procedures, or data-exclusion criteria, leaving the quantitative evidence for the functional form unverified.

    Authors: We agree that the abstract statement on numerical support requires more detail. The manuscript contains comparisons to lattice data, but we will revise both the abstract and the numerical section to include error analysis, a description of the fitting procedures used to test the functional form, and explicit criteria for data selection and exclusion. revision: yes

Circularity Check

1 steps flagged

Power-law form α_S(μ) ≃ Λ_S²/μ² follows directly from scale-invariance assumption on gluon condensate

specific steps
  1. self definitional [Abstract, central paragraph on the gluon condensate]
    "A central point is that the gluon condensate is scale invariant, which reflects its self-similar behavior across different scales. Building on that, we derive α_S(μ) ≃ Λ_S²/μ², which evolves to the infrared fixed point 1/α_S = 0 in accordance with infrared slavery."

    The scale invariance is taken as given input; the specific power-law form α_S(μ) ≃ Λ_S²/μ² (with Λ_S tied to the condensate) and the resulting IR fixed point are then obtained directly from that assumption, reducing the claimed derivation to a restatement of the input rather than an independent prediction from the RG flow.

full rationale

The paper states that the gluon condensate is scale invariant as a central point reflecting self-similar behavior, then builds the explicit functional form of the running coupling on that assumption. This makes the derived expression and its IR fixed point equivalent to the input invariance by construction, rather than an independent RG outcome. The derivation chain is otherwise presented as first-principles via gradient flow, with no evident self-citation load-bearing or parameter fitting in the abstract. The result is partially circular because the key functional dependence is fixed by the declared invariance.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the gluon condensate being universal and scale-invariant; Λ_S appears as a scale parameter whose numerical value is not independently derived in the abstract.

free parameters (1)
  • Λ_S
    Scale parameter appearing in the derived coupling; its value is set by the gluon condensate or by matching to data.
axioms (1)
  • domain assumption The gluon condensate is scale invariant across different renormalization scales.
    Stated as the central point that allows the functional form α_S(μ) ≃ Λ_S²/μ² to be written down.

pith-pipeline@v0.9.0 · 5715 in / 1305 out tokens · 40380 ms · 2026-05-21T22:54:03.466164+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

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 echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    A central point is that the gluon condensate is scale invariant, which reflects its self-similar behavior across different scales. Building on that, we derive α_S(μ) ≃ Λ_S²/μ², which evolves to the infrared fixed point 1/α_S = 0

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.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages · 17 internal anchors

  1. [1]

    Narayanan and H

    R. Narayanan and H. Neuberger, JHEP 0603 (2006) 064 [hep-th /0601210]

  2. [2]

    Properties and uses of the Wilson flow in lattice QCD

    M. L¨ uscher, JHEP 1008 (2010) 071 [Erratum: JHEP 1403 (2014) 092] [arXiv:1006.4518 [hep-lat]]

  3. [3]

    R. V . Harlander and T. Neumann, JHEP 06 (2016) 161 [arXiv:1606.03756 [hep-ph]]

  4. [4]

    M. A. Shifman, A. I. V ainshtein and V . I. Zakharov, Nucl. P hys. B 147 (1979) 385; Nucl. Phys. B 147 (1979) 448

  5. [5]

    Perturbative analysis of the gradient flow in non-abelian gauge theories

    M. L¨ uscher and P . Weisz, JHEP02 (2011) 051 [arXiv:1101.0963 [hep-th]]. 14

  6. [6]

    J. Artz, R. V . Harlander, F. Lange, T. Neumann and M. Praus a, JHEP 06 (2019) 121 [erra- tum: JHEP 10 (2019) 032] [arXiv:1905.00882 [hep-lat]]

  7. [7]

    Banks, R

    T. Banks, R. Horsley, H. R. Rubinstein and U. Wol ff, Nucl. Phys. B 190 (1981) 692

  8. [8]

    Energy-momentum tensor from the Yang-Mills gradient flow

    H. Suzuki, PTEP 2013 (2013) 083B03 [erratum: PTEP 2015 (2015) 079201] [arXiv:1304.0533 [hep-lat]]

  9. [9]

    Kluberg-Stern and J

    H. Kluberg-Stern and J. B. Zuber, Phys. Rev. D 12 (1975) 467

  10. [10]

    Tarrach, Nucl

    R. Tarrach, Nucl. Phys. B 196 (1982) 45

  11. [11]

    V . A. Novikov, M. A. Shifman, A. I. V ainshtein and V . I. Za kharov, Nucl. Phys. B 165 (1980) 67

  12. [12]

    F. R. Graziani, Z. Phys. C 33 (1987) 397

  13. [13]

    Suzuki and H

    H. Suzuki and H. Takaura, PTEP 2019 (2019) 103B04 [arXiv:1807.10064 [hep-ph]]

  14. [14]

    Schierholz, arXiv:2410.17677 [hep-lat], to be publ ished in Eur

    G. Schierholz, arXiv:2410.17677 [hep-lat], to be publ ished in Eur. Phys. J. C

  15. [15]

    C. T. Hill, arXiv:hep-th /0510177 [hep-th]

  16. [16]

    G. K. Savvidy, Phys. Lett. B 71 (1977) 133

  17. [17]

    N. K. Nielsen and P . Olesen, Nucl. Phys. B 144 (1978) 376

  18. [18]

    Nakamura and G

    Y . Nakamura and G. Schierholz, Nucl. Phys. B 986 (2023) 116063 [arXiv:2106.11369 [hep- ph]]

  19. [19]

    Schierholz, J

    G. Schierholz, J. Phys. G 52 (2025) 04LT01 [arXiv:2403.13508 [hep-ph]]

  20. [20]

    G¨ ockeler, A

    M. G¨ ockeler, A. S. Kronfeld, M. L. Laursen, G. Schierho lz and U. J. Wiese, Phys. Lett. B 233 (1989) 192

  21. [21]

    E. M. Ilgenfritz, M. L. Laursen, G. Schierholz, M. M¨ ull er-Preussker and H. Schiller, Nucl. Phys. B 268 (1986) 693

  22. [22]

    Wilson loops to 20th order numerical stochastic perturbation theory

    R. Horsley, G. Hotzel, E. M. Ilgenfritz, R. Millo, H. Per lt, P . E. L. Rakow, Y . Nakamura, G. Schierholz and A. Schiller, Phys. Rev. D 86 (2012) 054502 [arXiv:1205.1659 [hep-lat]]

  23. [23]

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

  24. [24]

    FLAG Review 2024

    Y . Aoki et al. [Flavour Lattice Averaging Group (FLAG)], arXiv:2411.04268 [hep-lat]

  25. [25]

    M. C` e, C. Consonni, G. P . Engel and L. Giusti, Phys. Rev. D 92 (2015) 074502 [arXiv:1506.06052 [hep-lat]]

  26. [26]

    V . A. Novikov, M. A. Shifman, A. I. V ainshtein and V . I. Za kharov, Nucl. Phys. B 191 (1981) 301. 15

  27. [27]

    I. E. Halperin and A. Zhitnitsky, Mod. Phys. Lett. A 13 (1998) 1955 [arXiv:hep-ph/9707286 [hep-ph]]

  28. [28]

    The infrared behaviour of quantum chro modynamics,

    D. P . Crewther, “The infrared behaviour of quantum chro modynamics,” doi:10.7907/JYNH-9286

  29. [29]

    On the relation between QCD potentials in momentum and position space

    M. Jezabek, M. Peter and Y . Sumino, Phys. Lett. B 428 (1998) 352 [arXiv:hep-ph/9803337 [hep-ph]]

  30. [30]

    J. L. Richardson, Phys. Lett. B 82 (1979) 272

  31. [31]

    The static potential in QCD to two loops

    Y . Schr¨ oder, Phys. Lett. B447 (1999) 321 [arXiv:hep-ph/9812205 [hep-ph]]

  32. [32]

    Testing perturbation theory on the nf=0 static quark potential

    S. Necco and R. Sommer, Phys. Lett. B 523 (2001) 135 [arXiv:hep-ph/0109093 [hep-ph]]

  33. [33]

    Quark and Color Confinement Through Dynamical Higgs Mechanism,

    G. Mack, “Quark and Color Confinement Through Dynamical Higgs Mechanism,” DESY preprint DESY 77 /58 (1977), doi:10.3204 /PUBDB-2017-12281

  34. [34]

    de Forcrand and M

    P . de Forcrand and M. Pepe, Nucl. Phys. B 598 (2001) 557 [arXiv:hep-lat /0008016 [hep- lat]]

  35. [35]

    Abelian and center gauge fixing in continuum Yang-Mills-Theory for general gauge groups

    H. Reinhardt and T. Tok, Phys. Lett. B 500 (2001) 173 [arXiv:hep-th/0011068 [hep-th]]

  36. [36]

    ’t Hooft, Nucl

    G. ’t Hooft, Nucl. Phys. B 190 (1981) 455

  37. [37]

    A. S. Kronfeld, G. Schierholz and U. J. Wiese, Nucl. Phys . B 293 (1987) 461

  38. [38]

    A. S. Kronfeld, M. L. Laursen, G. Schierholz and U. J. Wie se, Phys. Lett. B 198 (1987) 516

  39. [39]

    K. U. Can, R. Horsley, Y . Nakamura, H. Perlt, P . E. L. Rako w, G. Schierholz, H. St¨ uben, R. D. Y oung and J. Zanotti PoS LATTICE2021 (2022) 055 [arXiv:2112.06355 [hep-lat]]

  40. [40]

    V . G. Bornyakov, R. Horsley, R. Hudspith, Y . Nakamura, H. Perlt, D. Pleiter, P . E. L. Rakow, G. Schierholz, A. Schiller and H. St¨ uben and J. M. Zanotti, arXiv:1508.05916 [hep-lat]

  41. [41]

    Flavour blindness and patterns of flavour symmetry breaking in lattice simulations of up, down and strange quarks

    W. Bietenholz, V . Bornyakov, M. G¨ ockeler, R. Horsley, W. G. Lockhart, Y . Nakamura, H. Perlt, D. Pleiter, P . E. L. Rakow and G. Schierholz, A. Schi ller, T. Streuer, H. St¨ uben , F. Winter and J. M. Zanotti, Phys. Rev. D 84 (2011) 054509 [arXiv:1102.5300 [hep-lat]]

  42. [42]

    G. F. de T´ eramond, A. Paul, S. J. Brodsky, A. Deur, H. G. D osch, T. Liu and R. S. Sufian, Phys. Rev. Lett. 133 (2024) 181901 [arXiv:2403.16126 [hep-ph]]

  43. [43]

    Monopole action and condensation in SU(2) QCD

    H. Shiba and T. Suzuki, Phys. Lett. B 351 (1995) 519 [arXiv:hep-lat/9408004 [hep-lat]]

  44. [44]

    Engelhardt and H

    M. Engelhardt and H. Reinhardt, Nucl. Phys. B 585 (2000) 591 [arXiv:hep-lat /9912003 [hep-lat]]

  45. [45]

    Sonoda and H

    H. Sonoda and H. Suzuki, arXiv:2506.23567 [hep-th]. 16

  46. [46]

    The Vacuum Polarization: Power Corrections beyond OPE ?

    M. G¨ ockeler, R. Horsley, W. K¨ urzinger, V . Linke, D. Pleiter, P . E. L. Rakow and G. Schier- holz, Nucl. Phys. B Proc. Suppl. 94 (2001) 571 [arXiv:hep-lat/0012010 [hep-lat]]

  47. [47]

    Narison, Z

    S. Narison, Z. Phys. C 26 (1984) 209

  48. [48]

    V . A. Novikov, M. A. Shifman, A. I. V ainshtein and V . I. Za kharov, Nucl. Phys. B 249 (1985) 445. 17