Pith. sign in

REVIEW 2 major objections 7 minor 130 references

Reviewed by Pith at T0; open to challenge.

T0 means a machine referee read the full paper against a public rubric. The mark states how deep the mechanical check went, never who wrote it. the ladder, T0–T4 →

T0 review · glm-5.2

Step scaling pins large-N Yang-Mills Λ at 0.639(36)

2026-07-09 18:14 UTC pith:XDH67IHL

load-bearing objection First step-scaling determination of the large-N Yang-Mills Λ-parameter; solid execution but the 2-loop perturbative matching is the load-bearing soft spot. the 2 major comments →

arxiv 2607.07176 v1 pith:XDH67IHL submitted 2026-07-08 hep-lat

The large-N Yang--Mills Λ-parameter from step scaling

classification hep-lat
keywords lambdamathrmlarge-overlinescriptscriptstylesqrtyang--millsdependence
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The Λ-parameter is the fundamental scale that dimensional transmutation generates in a Yang-Mills theory: it sets the confinement scale and is a renormalization-group invariant, though it depends on the chosen renormalization scheme. This paper determines how Λ depends on the number of colors N in pure SU(N) Yang-Mills theory, using the step-scaling method — a non-perturbative procedure that recursively runs a finite-volume coupling from hadronic scales up to perturbative scales by halving the box size at each step. The calculation is performed in the Twisted Gradient Flow scheme, which combines a gradient-flow-defined coupling with twisted boundary conditions on an asymmetric volume (long sides ℓ, short sides ℓ/N), so that finite-volume effects are controlled by the large extent ℓ rather than the physical volume ℓ⁴/N². This makes the scheme efficient for large N. Results at N=3, 5, and 8 give √(8t₀)Λ_MS = 0.577(23), 0.632(32), and 0.611(43), extrapolating to a large-N limit of 0.639(36) with N-dependence √(8t₀)Λ_MS(N) = 0.639(36)[1 − 0.85(62)/N² + O(1/N⁴)]. The large-N value is only about 11% above the physical N=3 result, and the 1/N² coefficient is consistent with zero within errors, meaning the confinement scale varies mildly with N. This is the first large-N determination of the Yang-Mills Λ-parameter that does not rely on asymptotic scaling — the alternative approach where one reads Λ off from the bare coupling's dependence on the lattice spacing assuming perturbative running, which suffers from a window problem requiring simultaneously large volumes and tiny lattice spacings. The results agree with previous asymptotic-scaling determinations at the ~5% level for all N studied.

Core claim

The Λ-parameter of pure SU(N) Yang-Mills theory — the dynamically generated scale of dimensional transmutation — varies only mildly with the number of colors. Using step scaling in a twisted-gradient-flow finite-volume scheme, the authors find that the large-N limit √(8t₀)Λ_MS(N=∞) = 0.639(36) exceeds the physical N=3 value by about 11%, with the N-dependence well described by a 1/N² expansion whose coefficient −0.85(62) is consistent with zero. This is the first large-N Λ determination achieved without asymptotic scaling, and the results are compatible with prior asymptotic-scaling estimates within current uncertainties.

What carries the argument

The step-scaling method in the Twisted Gradient Flow (TGF) scheme: a finite-volume renormalization scheme where the coupling is defined through the gradient flow (a smoothing procedure for gauge fields) evaluated at a flow time tied to the box size, and twisted boundary conditions on an asymmetric ℓ² × (ℓ/N)² volume suppress finite-volume effects by making the effective size ℓ rather than the physical volume. The coupling is projected to the topological sector Q=0 to avoid topological freezing at fine lattice spacings. The step-scaling function — which gives the coupling at half the box size given the coupling at full size — is extracted from a global fit of lattice data at multiple lattices

Load-bearing premise

The dominant uncertainty comes from extrapolating to the perturbative limit using only the two-loop (universal) β-function in the TGF scheme. At the weakest couplings reached (λ_pt ≈ 1.83), the extrapolation assumes that unknown higher-order perturbative corrections are small enough to be captured by a linear fit in λ_pt. If the three-loop coefficient is anomalously large, the extrapolated Λ values could shift beyond the quoted errors.

What would settle it

A computation of the three-loop β-function coefficient in the TGF scheme that is anomalously large would shift the λ_pt → 0 extrapolation and could move the Λ values outside the quoted error bars.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the three-loop β-function coefficient in the TGF scheme were computed (e.g., via numerical stochastic perturbation theory), the dominant uncertainty from the λ_pt → 0 extrapolation would drop from linear to quadratic in the coupling, potentially reaching per-cent-level precision for Λ.
  • The mild N-dependence (coefficient consistent with zero) suggests confinement-scale physics is remarkably insensitive to the number of colors, supporting the use of large-N approximations for real-world QCD phenomenology.
  • The framework could be extended to a fully volume-reduced twisted Eguchi-Kawai formulation where the step-scaling sequence is generated by varying N itself, directly accessing N=∞ without finite-N extrapolation.
  • Agreement with asymptotic-scaling results at the ~5% level for N>3 helps validate both methodologies, but per-cent-level precision would be needed to resolve whether the tensions discussed in the N=3 literature reflect real systematic differences or are statistical fluctuations.
  • The pure-gauge Λ values feed into decoupling-based strategies for extracting the strong coupling in full QCD, so improved large-N determinations would tighten a key input for precision α_s determinations.

Where Pith is reading between the lines

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

  • With only three data points (N=3, 5, 8) and a 1/N² coefficient consistent with zero, the functional form of the N-dependence is not tightly constrained — a 1/N correction or other non-analytic terms cannot be excluded at current precision.
  • The ~11% gap between N=3 and N=∞, if confirmed at higher precision, would quantify the systematic error introduced by large-N approximations in hadron phenomenology applications.
  • The agreement between step-scaling and asymptotic-scaling results at N>3, where error bars are larger, may mask real differences that are visible only at N=3 where both methods have been pushed to higher precision — the consistency at N>3 could be a consequence of larger uncertainties rather than genuine agreement.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 7 minor

Summary. This paper presents the first determination of the Yang--Mills Λ-parameter at large N using step-scaling rather than asymptotic-scaling methods. The authors implement step-scaling in the Twisted Gradient Flow (TGF) scheme with twisted boundary conditions for SU(N) at N=3, 5, and 8, extract Λ_MS in units of √(8t₀) for each N, and perform a 1/N² extrapolation to N=∞. The methodology is carefully implemented: continuum extrapolations are performed at multiple lattice spacings, lattice artifacts are parameterized and controlled, the step-scaling function is determined from global fits, and scale setting uses independent infinite-volume simulations from Ref. [109]. The resulting large-N value √(8t₀)Λ_MS(N=∞) = 0.639(36) is compatible with existing asymptotic-scaling determinations within errors.

Significance. This is a methodologically important first: all prior large-N determinations of Λ relied on asymptotic scaling, which suffers from a window problem requiring simultaneously large volumes and fine lattice spacings. The step-scaling approach sidesteps this, and the use of twisted boundary conditions provides effective volume reduction advantageous for large-N studies. The N=3 result is cross-checked against prior work (Ref. [27]) using the same data but with non-perturbative SF matching, and the SU(4) step-scaling result of Ref. [24] is shown to be consistent with the fitted N-dependence. The results are compared comprehensively with the literature in Sec. V. The work is relevant both for large-N gauge dynamics and for decoupling-based strategies for α_s extraction in full QCD.

major comments (2)
  1. Sec. IV D and Fig. 5: The dominant systematic uncertainty is the λ_pt → 0 extrapolation, performed linearly using only the 2-loop universal β-function in the TGF scheme. For N=3, the paper itself provides a quantitative anchor: the SF-matched result of Ref. [27] gives √(8t₀)Λ_MS = 0.607(17), compared to the TGF-only value of 0.577(23) — a ~5% upward shift when higher-order perturbative information is included. This shift is comparable to the quoted errors and directly affects the large-N extrapolation. The critical question for the central large-N claim is whether this ~5% truncation error is N-dependent. The 3-loop coefficient of the TGF β-function is scheme-dependent; while its leading large-N behavior is universal, subleading N-dependence could cause differential shifts at N=5 and N=8, biasing the fitted 1/N² slope A = −0.85(62) (currently only 1.4σ from zero) and the N=∞ intercept. I
  2. The authors acknowledge this issue in the outlook (Sec. V) but do not quantify the potential impact on the large-N extrapolation. Given that the large-N fit uses only 3 data points for 2 parameters with unconstrained O(1/N⁴) terms, it has limited resilience against an N-dependent systematic of this magnitude. The authors should provide a more quantitative discussion: for instance, using the N=3 SF comparison to estimate the size of the 3-loop truncation effect, and arguing (or at least discussing) whether the N-dependence of this effect is expected to be small compared to the fitted 1/N² coefficient. Without this, the reader cannot assess the robustness of the large-N intercept.
minor comments (7)
  1. Sec. IV C: The choice of fit orders d₁ = d₂ = 4 for N=3,8 and d₁=4, d₂=5 for N=5 is stated to be stable under variation of d₁, d₂ between 4 and 6, but no quantitative detail on this stability check is provided (e.g., how much the continuum step-scaling function changes). A brief statement of the range of variation observed would strengthen confidence.
  2. Sec. IV F, Eq. (53): The large-N fit assumes O(1/N⁴) corrections are negligible, but with only 3 points and 2 parameters, these are completely unconstrained. It would be useful to show the result of alternative fit ansätze (e.g., a constant fit with no 1/N² term, or a fit including a 1/N⁴ term with a prior) to demonstrate that the intercept is not an artifact of the chosen parameterization.
  3. Table II: The λ_pt → 0 extrapolated values (k=∞ row) appear to be obtained from a linear fit using only the last few points, but the number of points used and the fit quality (χ²/dof) are not stated. This information should be provided, perhaps in a caption or footnote.
  4. Sec. III, regarding the topological Q=0 projection (Eq. 28): The consistency between projected and non-projected coupling definitions was verified for N=3 in Ref. [91]. For N=5 and N=8, no such cross-check is mentioned. A brief comment on whether this has been checked, or why it is expected to hold, would be appropriate.
  5. Fig. 2: The two Padé fit ranges are described as motivated by a change in behavior around b ≈ 0.40–0.45, but the physical origin of this change is not explained. A brief comment would help the reader.
  6. Sec. IV E: The statement that √(8t₀) ≃ 0.475 fm is adopted for all N is a convention choice, but it would be worth noting explicitly that this is a definition rather than a physical measurement for N≠3, to avoid confusion.
  7. Typo in Sec. IV C: 'In the end, we quote results...' — the 'I' at the start appears to be a formatting artifact.

Circularity Check

0 steps flagged

No significant circularity: the Λ-parameter is extracted from non-perturbative step-scaling data matched to perturbation theory, with independent scale-setting; one minor self-citation for the N=3 dataset reuse is not load-bearing for the large-N claim.

full rationale

The paper's central result — the large-N extrapolation of √(8t₀)Λ_MS — is derived from a chain that is not circular. The step-scaling function is determined non-perturbatively from lattice simulations at N=3,5,8 (Sec. IV C, Eq. 46), the coupling is run to the perturbative regime and matched to the 2-loop β-function (Sec. IV D, Eq. 8), and the Λ-parameter in each scheme is extracted from a λ_pt→0 extrapolation (Fig. 5). The scale-setting conversion factor μ_had√(8t₀) comes from independent infinite-volume simulations in Ref. [109] (a companion paper by overlapping authors, but with separate simulation data and a distinct methodology using PTBC). The N=3 result reuses simulation data from Ref. [27] (Bribián 21), but the paper explicitly states the analysis differs (different μ_had choice, TGF-only matching vs. SF matching) and cross-checks against the SF-matched result (0.577 vs. 0.607), showing the reuse is not a definitional identity. The large-N fit (Eq. 53) is a standard 1/N² extrapolation over three independent data points, not a fit that recovers its inputs by construction. The one-loop TGF→MS conversion factor (Eq. 21) is a perturbative calculation, not a fitted input. The self-citations (Refs. [27], [109]) provide supporting infrastructure but the central derivation chain — lattice data → step-scaling function → perturbative matching → Λ extraction → large-N extrapolation — has independent content at each stage. The dominant systematic (2-loop truncation in the λ_pt→0 extrapolation) is a correctness concern, not a circularity: the paper does not define its output in terms of that extrapolation, it uses it as a controlled approximation with stated assumptions. No step reduces to its inputs by construction.

Axiom & Free-Parameter Ledger

5 free parameters · 5 axioms · 0 invented entities

No new physical entities are postulated. The free parameters are standard choices in lattice step-scaling calculations (scheme parameter c, starting coupling u₀, fit orders). The large-N fit parameters A and Λ(∞) are the results being determined. The axioms are all domain-standard assumptions in lattice gauge theory, not ad hoc inventions.

free parameters (5)
  • u₀ = λ_TGF(μ_had) = 44.81
    Starting coupling for the step-scaling sequence, chosen to balance scale-setting feasibility against the strong-coupling regime (Sec. IV B). Not fitted to data but selected by the authors.
  • c = 0.3
    Scheme parameter defining the renormalization scale μ=1/(cℓ) in the TGF scheme (Sec. II B). Chosen freely; each value defines a different scheme.
  • d₁, d₂ (fit orders) = d₁=d₂=4 for N=3,8; d₁=4,d₂=5 for N=5
    Orders of the polynomial expansion of the inverse step-scaling function and its lattice artifacts (Sec. IV C). Chosen after checking stability by varying between 4 and 6.
  • A (large-N slope) = -0.85(62)
    Coefficient of 1/N² in the large-N extrapolation (Eq. 53), fitted to the three data points N=3,5,8.
  • √(8t₀)Λ_MS(∞) = 0.639(36)
    Large-N limit of the Λ-parameter, fitted from the same three data points.
axioms (5)
  • domain assumption Large-N volume independence via twisted boundary conditions holds at finite N with controlled finite-volume effects
    The TGF scheme relies on twisted volume reduction being effective at finite N, with finite-volume effects controlled by the effective symmetric volume ℓ⁴ rather than ℓ⁴/N² (Sec. II B). This is supported by perturbation theory but is an assumption at finite N.
  • domain assumption Lattice artifacts are O(a²) = O(1/L²) at leading order
    Assumed in the continuum extrapolations of both the step-scaling function (Eq. 45) and the scale-setting conversion factor (Sec. IV E, Fig. 6). Wilson plaquette action is O(a)-improved by the clover discretization.
  • domain assumption The 1/N² expansion of √(8t₀)Λ_MS(N) is the correct functional form for the large-N extrapolation
    Standard large-N counting predicts O(1/N²) corrections, but with only three data points (N=3,5,8) the O(1/N⁴) term cannot be constrained (Sec. IV F, Eq. 53).
  • domain assumption The Q=0 topological projection yields a valid renormalization scheme for step-scaling
    The coupling is defined with δ_{Q,0} projection (Eq. 28). This was justified in Ref. [104] and cross-checked in Ref. [91], but it is a scheme choice that could in principle introduce systematics.
  • domain assumption Two-loop perturbative matching is sufficient with a linear λ_pt→0 extrapolation
    Only the universal two-loop β-function is known in the TGF scheme. The authors assume O(λ_pt²) corrections are negligible at their weakest couplings (λ_pt≈1.83) and extrapolate linearly (Sec. IV D, Fig. 5). This is the dominant systematic uncertainty.

pith-pipeline@v1.1.0-glm · 34998 in / 3116 out tokens · 186175 ms · 2026-07-09T18:14:18.692710+00:00 · methodology

0 comments
read the original abstract

We use the step-scaling method and results obtained at $N = 3, 5$ and $8$ to determine the $N$-dependence of the dynamically generated scale $\Lambda$ of $\mathrm{SU}(N)$ Yang--Mills theories. We implement the step-scaling method in a suitable finite-volume renormalization scheme based on twisted boundary conditions, introduced to effectively achieve large-$N$ volume independence, and on a coupling defined through the gradient flow. In the $\overline{\mathrm{MS}}$ scheme, we obtain the following values in terms of the gradient flow scale $t_0$: $\sqrt{8t_0}\Lambda_{\scriptscriptstyle{\overline{\mathrm{MS}}}} = 0.577(23)$, $0.632(32)$, and $0.611(43)$ for $N=3,5$ and $8$, respectively. They extrapolate to a large-$N$ value of: $\sqrt{8t_0}\Lambda_{\scriptscriptstyle{\overline{\mathrm{MS}}}} (N=\infty) = 0.639(36)$, and the $N$-dependence is given by $\sqrt{8t_0}\Lambda_{\scriptscriptstyle{\overline{\mathrm{MS}}}}(N)=0.639(36)[1-0.85(62)/N^2+\mathcal{O}(1/N^4)]$. This work represents the first calculation of the Yang--Mills $\Lambda$-parameter in the large-$N$ limit that does not rely on asymptotic scaling strategies.

Figures

Figures reproduced from arXiv: 2607.07176 by Andrea Giorgieri, Claudio Bonanno, Jorge Luis Dasilva Gol\'an, Margarita Garc\'ia P\'erez.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: the central panel of [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

130 extracted references · 130 canonical work pages · 87 internal anchors

  1. [1]

    e−1/[2b0λs(µ)] e−Is(λs(µ)) ,(3) Is(λ) = Z λ 0 dx 1 2βs(x) + 1 2b0x2 − b1 2b2 0x .(4) The values of Λ in two schemes s,s ′ can be matched ex- actly through a one-loop calculation: log Λs′ Λs = css′ 2b0 ,(5) with λs′ =λ s (1 +c ss′λs +. . .).(6) By introducing two scalesµ 1 andµ 2, one obtains the exact relation: Λs µ1 = Λs µ2 exp ( − Z λs(µ1) λs(µ2) dx 2βs...

  2. [2]

    (4) evaluated with then-loopβ-function

    e − 1 2b0 λs (µpt ) e−I (n) s [λs(µpt)] +O[λ n−1 s (µpt)], (8) whereI (n) s (λ) is the perturbative truncation ofI s(λ) in Eq. (4) evaluated with then-loopβ-function. Instead, the non-perturbative running fromµ had toµ pt can be controlled on the lattice. In this case, the integral in the exponent of Eq. (7) reduces to, − Z λs(µhad) λs(µpt) dx 2βs(x) =− Z...

  3. [3]

    This is discussed in Sec

    Determineu 1 =λ(2µ had), withµ had defined by the chosenu 0 =λ(µ had), by fitting theb-dependence of λon the lattice. This is discussed in Sec. IV B

  4. [4]

    This is discussed in Sec

    From a global fit of the values ofλon the lattice, determine the inverse step-scaling functionu(σ) in a range of couplingsσthat includesu 1 and extends to the perturbative regime. This is discussed in Sec. IV C

  5. [5]

    This is discussed in Sec

    Determine Λ MS/µhad by usingu(σ) to run the cou- pling to sufficiently high energies in order to match with perturbation theory. This is discussed in Sec. IV D

  6. [6]

    Non- perturbative aspects of fundamental interactions, in the Standard Model and beyond

    Convert Λ MS/µhad to Λ MS √8t0 using the results of scale setting to determineµ had √8t0. This is dis- cussed in Sec. IV E. B. Initial step fromµ had to2µ had In this section, we describe the first step of the step- scaling sequence, that is, the choice ofu 0 =λ TGF(µhad) and the determination ofu 1 =λ TGF(2µhad). As antic- ipated in Sec. IV A, this requi...

  7. [7]

    ’t Hooft, Nucl

    G. ’t Hooft, Nucl. Phys. B72, 461 (1974)

  8. [8]

    Witten, Nucl

    E. Witten, Nucl. Phys. B149, 285 (1979)

  9. [9]

    ’t Hooft, Phys

    G. ’t Hooft, Phys. Rev. Lett.37, 8 (1976)

  10. [10]

    Witten, Nucl

    E. Witten, Nucl. Phys. B156, 269 (1979)

  11. [11]

    Veneziano, Nucl

    G. Veneziano, Nucl. Phys. B159, 213 (1979)

  12. [12]

    Kawarabayashi and N

    K. Kawarabayashi and N. Ohta, Nucl. Phys. B175, 477 (1980)

  13. [13]

    Witten, Annals Phys.128, 363 (1980)

    E. Witten, Annals Phys.128, 363 (1980)

  14. [14]

    Di Vecchia and G

    P. Di Vecchia and G. Veneziano, Nucl. Phys. B171, 253 (1980)

  15. [15]

    SU(N) gauge theories at large N

    B. Lucini and M. Panero, Phys. Rept.526, 93 (2013), arXiv:1210.4997 [hep-th]

  16. [16]

    Introductory lectures to large-N QCD phenomenology and lattice results

    B. Lucini and M. Panero, Prog. Part. Nucl. Phys.75, 1 (2014), arXiv:1309.3638 [hep-th]

  17. [17]

    The Large $N_c$ limit of QCD on the lattice

    P. Hern´ andez and F. Romero-L´ opez, Eur. Phys. J. A57, 52 (2021), arXiv:2012.03331 [hep-lat]. 15

  18. [18]

    Prospects for large N gauge theories on the lattice

    M. Garc´ ıa P´ erez, PoSLA TTICE2019, 276 (2020), arXiv:2001.10859 [hep-lat]

  19. [19]

    Lucini and M

    B. Lucini and M. Teper, JHEP06(2001), 050, arXiv:hep- lat/0103027 [hep-lat]

  20. [20]

    On the running of the bare coupling in SU(N) lattice gauge theories

    C. Allton, M. Teper, and A. Trivini, JHEP07(2008), 021, arXiv:0803.1092 [hep-lat]

  21. [21]

    Rectangular Wilson Loops at Large N

    R. Lohmayer and H. Neuberger, JHEP08(2012), 102, arXiv:1206.4015 [hep-lat]

  22. [22]

    The string tension from smeared Wilson loops at large N

    A. Gonzalez-Arroyo and M. Okawa, Phys. Lett. B718, 1524 (2013), arXiv:1206.0049 [hep-th]

  23. [23]

    SU(N) gauge theories in 3+1 dimensions: glueball spectrum, string tensions and topology

    A. Athenodorou and M. Teper, JHEP12(2021), 082, arXiv:2106.00364 [hep-lat]

  24. [24]

    Testing (asymptotic) scaling in Yang-Mills theories in the large-$N_c$ limit

    P. Butti and A. Gonz´ alez-Arroyo, PoSLA TTICE2023, 381 (2024), arXiv:2311.18696 [hep-lat]

  25. [25]

    L¨ uscher, P

    M. L¨ uscher, P. Weisz, and U. Wolff, Nucl. Phys. B359, 221 (1991)

  26. [26]

    Infinite N phase transitions in continuum Wilson loop operators

    R. Narayanan and H. Neuberger, JHEP03(2006), 064, arXiv:hep-th/0601210

  27. [27]

    Continuous smearing of Wilson Loops

    R. Lohmayer and H. Neuberger, PoSLA TTICE2011, 249 (2011), arXiv:1110.3522 [hep-lat]

  28. [28]

    Trivializing maps, the Wilson flow and the HMC algorithm

    M. L¨ uscher, Commun. Math. Phys.293, 899 (2010), arXiv:0907.5491 [hep-lat]

  29. [29]

    ’t Hooft, Nucl

    G. ’t Hooft, Nucl. Phys. B153, 141 (1979)

  30. [30]

    The running of the coupling in SU(N) pure gauge theories

    B. Lucini and G. Moraitis, Phys. Lett. B668, 226 (2008), arXiv:0805.2913 [hep-lat]

  31. [31]

    The gradient flow running coupling with twisted boundary conditions

    A. Ramos, JHEP11(2014), 101, arXiv:1409.1445 [hep-lat]

  32. [32]

    E. I. Bribi´ an and M. Garc´ ıa P´ erez, JHEP03(2019), 200, arXiv:1903.08029 [hep-lat]

  33. [33]

    E. I. Bribi´ an, J. L. Dasilva Gol´ an, M. Garc´ ıa P´ erez, and A. Ramos, Eur. Phys. Jour. C81, 951 (2021), arXiv:2107.03747 [hep-lat]

  34. [34]

    Eguchi and H

    T. Eguchi and H. Kawai, Phys. Rev. Lett.48, 1063 (1982)

  35. [35]

    Gonz´ alez-Arroyo and M

    A. Gonz´ alez-Arroyo and M. Okawa, Phys. Lett. B120, 174 (1983)

  36. [36]

    Gonz´ alez-Arroyo and M

    A. Gonz´ alez-Arroyo and M. Okawa, Phys. Rev. D27, 2397 (1983)

  37. [37]

    Large $N$ reduction with the Twisted Eguchi-Kawai model

    A. Gonz´ alez-Arroyo and M. Okawa, JHEP07(2010), 043, arXiv:1005.1981 [hep-th]

  38. [38]

    Testing volume independence of SU(N) pure gauge theories at large N

    A. Gonz´ alez-Arroyo and M. Okawa, JHEP12(2014), 106, arXiv:1410.6405 [hep-lat]

  39. [39]

    Non-perturbative renormalization by decoupling

    M. Dalla Brida, R. H¨ ollwieser, F. Knechtli, T. Korzec, A. Ramos, and R. Sommer (ALPHA), Phys. Lett. B807, 135571 (2020), arXiv:1912.06001 [hep-lat]

  40. [40]

    Determination of $\alpha_s(m_Z)$ by the non-perturbative decoupling method

    M. Dalla Brida, R. H¨ ollwieser, F. Knechtli, T. Korzec, A. Nada, A. Ramos, S. Sint, and R. Sommer (ALPHA), Eur. Phys. J. C82, 1092 (2022), arXiv:2209.14204 [hep-lat]

  41. [41]

    Dalla Brida, R

    M. Dalla Brida, R. H¨ ollwieser, F. Knechtli, T. Korzec, A. Ramos, S. Sint, and R. Sommer, Nature652, 328 (2026)

  42. [42]

    Past, present, and future of precision determinations of the QCD coupling from lattice QCD

    M. Dalla Brida, Eur. Phys. J. A57, 66 (2021), arXiv:2012.01232 [hep-lat]

  43. [43]

    Lattice determinations of the strong coupling

    L. Del Debbio and A. Ramos, Phys. Rept.920, 1 (2021), arXiv:2101.04762 [hep-lat]

  44. [44]

    The strong coupling constant: State of the art and the decade ahead

    D. d’Enterriaet al., J. Phys. G51, 090501 (2024), arXiv:2203.08271 [hep-ph]

  45. [45]

    Alpha_s from the non-perturbatively renormalised lattice three-gluon vertex

    B. Alles, D. Henty, H. Panagopoulos, C. Parrinello, C. Pit- tori, and D. G. Richards, Nucl. Phys. B502, 325 (1997), arXiv:hep-lat/9605033

  46. [46]

    Capitani, M

    S. Capitani, M. L¨ uscher, R. Sommer, and H. Wittig, Nuclear Physics B544, 669 (1999)

  47. [47]

    Three-loop beta function and non-perturbative $\alpha_s$ in asymmetric momentum scheme

    P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, and C. Roiesnel, JHEP12(1998), 004, arXiv:hep-ph/9810437

  48. [48]

    Lattice calculation of $\alpha_s$ in momentum scheme

    P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, and C. Roiesnel, JHEP10(1998), 017, arXiv:hep-ph/9810322

  49. [49]

    Asymptotic scaling of the gluon propagtor on the lattice

    D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodriguez-Quintero, and C. Roiesnel, Phys. Rev. D61, 114508 (2000), arXiv:hep-ph/9910204

  50. [50]

    Asymptotic behaviour of the gluon propagator from lattice QCD

    D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodriguez-Quintero, and C. Roiesnel, Phys. Rev. D60, 094509 (1999), arXiv:hep-ph/9903364

  51. [51]

    Consistent OPE Description of Gluon Two- and Three-point Green Function?

    P. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene, and J. Rodriguez-Quintero, Phys. Lett. B493, 315 (2000), arXiv:hep-ph/0008043

  52. [52]

    Boucaudet al., JHEP04(2000), 006, arXiv:hep- ph/0003020

    P. Boucaudet al., JHEP04(2000), 006, arXiv:hep- ph/0003020

  53. [53]

    Remarks on the determination of the Landau gauge OPE for the Asymmetric three gluon vertex

    F. De Soto and J. Rodriguez-Quintero, Phys. Rev. D64, 114003 (2001), arXiv:hep-ph/0105063

  54. [54]

    Testing Landau gauge OPE on the Lattice with a $<A^2>$ Condensate

    P. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene, and J. Rodriguez-Quintero, Phys. Rev. D63, 114003 (2001), arXiv:hep-ph/0101302

  55. [55]

    G¨ ockeler, R

    M. G¨ ockeler, R. Horsley, A. C. Irving, D. Pleiter, P. E. L. Rakow, G. Schierholz, and H. St¨ uben (QCDSF-UKQCD Col- laboration), Phys. Rev. D73, 014513 (2006)

  56. [56]

    Asymptotic behavior of the ghost propagator in SU3 lattice gauge theory

    P. Boucaud, J. P. Leroy, A. Le Yaouanc, A. Y. Lokhov, J. Micheli, O. Pene, J. Rodriguez-Quintero, and C. Roiesnel, Phys. Rev. D72, 114503 (2005), arXiv:hep-lat/0506031

  57. [57]

    Ghost-gluon running coupling, power corrections and the determination of $\Lambda_{\bar {\rm MS}}$

    P. Boucaud, F. De Soto, J. P. Leroy, A. Le Yaouanc, J. Micheli, O. Pene, and J. Rodriguez-Quintero, Phys. Rev. D79, 014508 (2009), arXiv:0811.2059 [hep-ph]

  58. [58]

    QCD Lambda parameter from Landau-gauge gluon and ghost correlations

    A. Sternbeck, E. M. Ilgenfritz, K. Maltman, M. Muller- Preussker, L. von Smekal, and A. G. Williams, PoS LA T2009, 210 (2009), arXiv:1003.1585 [hep-lat]

  59. [59]

    SU(3) Landau gauge gluon and ghost propagators using the logarithmic lattice gluon field definition

    E.-M. Ilgenfritz, C. Menz, M. Muller-Preussker, A. Schiller, and A. Sternbeck, Phys. Rev. D83, 054506 (2011), arXiv:1010.5120 [hep-lat]

  60. [60]

    Precision determination of r_0Lambda_MS from the QCD static energy

    N. Brambilla, X. Garcia i Tormo, J. Soto, and A. Vairo, Phys. Rev. Lett.105, 212001 (2010), [Erratum: Phys. Rev. Lett. 108 (2012) 269903], arXiv:1006.2066 [hep-ph]

  61. [61]

    Determination of LambdaMS from the gluon and ghost propagators in Landau gauge

    A. Sternbeck, K. Maltman, M. Muller-Preussker, and L. von Smekal, PoSLA TTICE2012, 243 (2012), arXiv:1212.2039 [hep-lat]

  62. [62]

    Twisted space-time reduced model of large N QCD with two adjoint Wilson fermions

    A. Gonz´ alez-Arroyo and M. Okawa, Phys. Rev. D88, 014514 (2013), arXiv:1305.6253 [hep-lat]

  63. [63]

    Equation of State for SU(3) Gauge Theory via the Energy-Momentum Tensor under Gradient Flow

    M. Kitazawa, T. Iritani, M. Asakawa, T. Hatsuda, and H. Suzuki, Phys. Rev. D94, 114512 (2016), arXiv:1610.07810 [hep-lat]

  64. [64]

    Non-perturbative determination of the $\Lambda$-parameter in the pure SU(3) gauge theory from the twisted gradient flow coupling

    K.-I. Ishikawa, I. Kanamori, Y. Murakami, A. Naka- mura, M. Okawa, and R. Ueno, JHEP12(2017), 067, arXiv:1702.06289 [hep-lat]

  65. [65]

    SU(3) Yang Mills theory at small distances and fine lattices

    N. Husung, M. Koren, P. Krah, and R. Sommer, EPJ Web Conf.175, 14024 (2018), arXiv:1711.01860 [hep-lat]

  66. [66]

    The gradient flow coupling at high-energy and the scale of SU(3) Yang-Mills theory

    M. Dalla Brida and A. Ramos, Eur. Phys. J. C79, 720 (2019), arXiv:1905.05147 [hep-lat]

  67. [67]

    An analysis of systematic effects in finite size scaling studies using the gradient flow

    A. Nada and A. Ramos, Eur. Phys. J. C81, 1 (2021), arXiv:2007.12862 [hep-lat]

  68. [68]

    Husung, A

    N. Husung, A. Nada, and R. Sommer, PoSLA TTICE2019, 263 (2020)

  69. [69]

    $\Lambda$ parameter of the SU(3) Yang-Mills theory from the continuous $\beta$ function

    A. Hasenfratz, C. T. Peterson, J. van Sickle, and O. Witzel, Phys. Rev. D108, 014502 (2023), arXiv:2303.00704 [hep-lat]

  70. [70]

    C. H. Wong, S. Borsanyi, Z. Fodor, K. Holland, and J. Kuti, PoSLA TTICE2022, 043 (2023)

  71. [71]

    Brambilla, V

    N. Brambilla, V. Leino, J. Mayer-Steudte, and A. Vairo, Phys. Rev. D109, 114517 (2024)

  72. [72]

    FLAG Review 2021

    Y. Aokiet al.(Flavour Lattice Averaging Group (FLAG)), Eur. Phys. J. C82, 869 (2022), arXiv:2111.09849 [hep-lat]

  73. [73]

    FLAG Review 2024

    Y. Aokiet al.(Flavour Lattice Averaging Group (FLAG)), Phys. Rev. D113, 014508 (2026), arXiv:2411.04268 [hep-lat]

  74. [74]

    Volume independence for Yang-Mills fields on the twisted torus

    M. Garc´ ıa P´ erez, A. Gonz´ alez-Arroyo, and M. Okawa, Int. J. Mod. Phys. A29, 1445001 (2014), arXiv:1406.5655 [hep-th]

  75. [75]

    The $SU(\infty)$ twisted gradient flow running coupling

    M. Garc´ ıa P´ erez, A. Gonz´ alez-Arroyo, L. Keegan, and M. Okawa, JHEP01(2015), 038, arXiv:1412.0941 [hep-lat]

  76. [76]

    Mass anomalous dimension of Adjoint QCD at large N from twisted volume reduction

    M. Garc´ ıa P´ erez, A. Gonz´ alez-Arroyo, L. Keegan, and M. Okawa, JHEP08(2015), 034, arXiv:1506.06536 [hep-lat]

  77. [77]

    Large N meson masses from a matrix model

    A. Gonz´ alez-Arroyo and M. Okawa, Phys. Lett. B755, 132 (2016), arXiv:1510.05428 [hep-lat]

  78. [78]

    Meson spectrum in the large $N$ limit

    M. Garc´ ıa P´ erez, A. Gonz´ alez-Arroyo, and M. Okawa, JHEP 04(2021), 230, arXiv:2011.13061 [hep-lat]

  79. [79]

    Scale setting for large-$N$ SUSY Yang-Mills on the lattice

    P. Butti, M. Garc´ ıa P´ erez, A. Gonz´ alez-Arroyo, K.- I. Ishikawa, and M. Okawa, JHEP07(2022), 074, arXiv:2205.03166 [hep-lat]

  80. [80]

    The large-$N$ limit of the chiral condensate from twisted reduced models

    C. Bonanno, P. Butti, M. Garc´ ıa Per´ ez, A. Gonz´ alez- 16 Arroyo, K.-I. Ishikawa, and M. Okawa, JHEP12(2023), 034, arXiv:2309.15540 [hep-lat]

Showing first 80 references.