computo ergo sum

2026-10-09 · chapter Yang–Mills and the mass gapPart 11 Further considerationsthe length at which conservation breaks, coordinates, two-dimensional test beds · Kiichi, Shiori and Rin

The whole book Contents and reading orders
  1. Part 0 Entrance — how to read this book
  2. 0-01 How to read this book
  3. 0-02 One picture
  4. 0-03 Reading the problem statement word by word
  5. 0-04 Four dimensions in mathematics and four dimensions in physics
  6. Part 1 Mathematical preparation
  7. 1-01 Vector spaces and linear maps
  8. 1-02 What a group is
  9. 1-03 Lie groups and Lie algebras
  10. 1-04 SU(2) and SU(3)
  11. 1-05 Manifolds and tangent spaces
  12. 1-06 Differential forms and the exterior derivative
  13. 1-07 Vector bundles and connections
  14. 1-08 Probability and measure
  15. 1-09 Probability in infinite dimensions
  16. 1-10 Hilbert spaces and self-adjoint operators
  17. 1-11 Fourier analysis and distributions
  18. 1-12 Detour: what happens only in four dimensions
  19. Part 2 Physical preparation
  20. 2-01 Lagrangians and the action
  21. 2-02 The idea of a field
  22. 2-03 Electromagnetism is a gauge theory
  23. 2-04 Special relativity and spacetime
  24. 2-05 The bones of quantum mechanics
  25. 2-06 The harmonic oscillator and creation and annihilation
  26. 2-07 The idea of the path integral
  27. 2-08 Statistical mechanics and phase transitions
  28. 2-09 Detour: the gap and the correlation length in the Ising model
  29. 2-10 The relaxation time and the dynamical exponent z
  30. Part 3 Yang–Mills theory (classical)
  31. 3-01 The gauge principle
  32. 3-02 Non-abelian gauge fields
  33. 3-03 The action and the equations
  34. 3-04 Gauge theory as geometry
  35. 3-05 Instantons and topology
  36. 3-06 Detour: Donaldson and exotic R⁴
  37. 3-07 Yang–Mills inside the Standard Model
  38. Part 4 Quantisation
  39. 4-01 Canonical quantisation and the Hamiltonian
  40. 4-02 The path integral and Euclideanisation
  41. 4-03 Perturbation theory and Feynman diagrams
  42. 4-04 Divergences and renormalisation
  43. 4-05 The algebra of divergences
  44. 4-06 Asymptotic freedom
  45. 4-07 Dimensional transmutation and Λ
  46. 4-08 Gauge fixing and Faddeev–Popov
  47. 4-09 The axioms of quantum field theory
  48. 4-10 The Osterwalder–Schrader axioms and reconstruction
  49. 4-11 The definition of the mass gap
  50. 4-12 Detour: the perturbation series diverges
  51. Part 5 Lattice gauge theory
  52. 5-01 Wilson's lattice
  53. 5-02 The strong-coupling expansion
  54. 5-03 Reflection positivity and the transfer matrix
  55. 5-04 The infinite-volume limit and the cluster expansion
  56. 5-05 Weak coupling and the continuum limit
  57. 5-06 How U(1) differs from the non-abelian case
  58. 5-07 The Monte Carlo method
  59. 5-08 Measuring glueballs and the string tension
  60. 5-09 What to hold fixed when taking the limit
  61. 5-10 Lattice gauge theory with finite groups
  62. 5-11 Detour: Bakry–Émery and one-link integration
  63. Part 6 Constructive field theory
  64. 6-01 What constructive field theory is
  65. 6-02 Solvability in two dimensions and the Yang–Mills measure
  66. 6-03 Constructing scalar fields
  67. 6-04 The cluster expansion
  68. 6-05 The sequence of renormalisation group steps
  69. 6-06 Ultraviolet stability of three-dimensional Yang–Mills
  70. 6-07 Four dimensions
  71. 6-08 Triviality of φ⁴ in four dimensions
  72. 6-09 Stochastic quantisation and regularity structures
  73. 6-10 Where four dimensions stops
  74. 6-11 Obstacles other than divergence
  75. 6-12 Detour: "if 1/0 is given a symbol"
  76. Part 7 From the side of physics
  77. 7-01 How physics sees it
  78. 7-02 Mechanisms of confinement
  79. 7-03 The string picture
  80. 7-04 Large N
  81. 7-05 String theory and extra dimensions
  82. 7-06 How extra dimensions become invisible
  83. 7-07 Where gauge fields come from
  84. 7-08 Holography
  85. 7-09 The mass gap becomes geometry
  86. 7-10 Back to four dimensions
  87. 7-11 Supersymmetry and Seiberg–Witten
  88. 7-12 The sentence that plays the role of an equivalence principle
  89. 7-13 Where the way physics digs parts from mathematics
  90. Part 8 A dictionary of two languages — the physical and the mathematical view
  91. 8-01 How to read the dictionary
  92. 8-02 Dictionary A
  93. 8-03 Dictionary B
  94. 8-04 Dictionary C
  95. 8-05 The physical counterparts of our assumptions
  96. Part 9 Where things stand, and open tasks
  97. 9-01 How far the world has come
  98. 9-02 Where the two ways of digging part
  99. 9-03 Candidates for new concepts
  100. 9-04 The list of open tasks
  101. 9-05 Common misconceptions
  102. Part 10 Towards a proof of the mass gap — building the first step of the staircase
  103. 10-00 Entrance — the staircase, and the random walk, the ruler and the clock
  104. 10-01 The line of construction and Assumption H
  105. 10-02 Roads tried on the four-dimensional lattice, and walls
  106. 10-03 The four steps of the staircase — why we start from 0+1 dimensions
  107. 10-04 The path measure of the double well, and kinks
  108. 10-05 The statement of the first step — four conditions and Proposition 4.5
  109. 10-06 (H1) Local Poincaré within a phase
  110. 10-07 The map of (H2) — seven routes and the seam of the ring
  111. 10-08 What computer-assisted means — the conditions on the chain inside the window
  112. 10-09 The sector framework and the three strands of route 5
  113. 10-10 The list of detour plans — types of breakage
  114. 10-11 (H3)(H4) and the constant of the core
  115. 10-12 The integrated table and the dependency graph
  116. 10-13 Lean statements and the limits of machine checking
  117. 10-14 Reconciliation — the numbers of the record and the known numbers
  118. 10-15 What remains and its cost — the distance to a theorem
  119. 10-16 What is not claimed
  120. 10-17 Insights and small questions
  121. Part 11 Further considerations — beyond Assumption H, and other views
  122. 11-00 Further considerations — entrance
  123. 11-01 One form that obstructs H — a nearly conserved global quantity
  124. 11-02 Choosing the ruler and the coordinates again — the length at which conservation breaks
  125. 11-03 Two ways to repair H — the sector H and the local H
  126. 11-04 Which H the one-sided bridge needs
  127. 11-05 Making time continuous — the Hamiltonian lattice and the vacuum measure
  128. 11-06 Paying with rarity instead of convexity — renormalisation on groups and log-Sobolev
  129. 11-07 Because it is four dimensions — the dimensionless coupling and topology
  130. 11-08 Questions, and what is not claimed
  131. 11-A1 Table of notation for Part 11
  132. Appendices
  133. A-01 Table of notation
  134. A-02 Glossary
  135. A-03 Guide to the literature
  136. A-04 Lean and machine checking
  137. A-05 How this book was made
  138. A-06 Index of assumptions

Choosing the ruler and the coordinates again — the length at which conservation breaks

What this chapter shows — Assumption H was a sentence comparing the relaxation time CPC_P with the square of one length ξ\xi. We split the rulers that could play the role of that ξ\xi into three: the static correlation length, the length at which conservation breaks ℓ∗\ell_*, and the time for the global quantity to move by one. On a hierarchical lattice we see that, even for the same equilibrium measure, whether H is true changes with the coordinates used to move it (the choice of Dirichlet form) (CP=1/βC_P=1/\beta in tree coordinates, the square of the side in link coordinates). We state that the length at which conservation breaks, measured at five points on the two-dimensional U(1)U(1) strip, is not proportional to the static length; that it comes from the same wrapping as one step of the topological charge; and that splitting the relaxation rate into a "raw rate" and a "static fluctuation" reveals the factor the paper prediction had dropped. Finally we set out a candidate correction that compares with max⁡(ξ2, the time for the global quantity to move by one)\max(\xi^2,\ \text{the time for the global quantity to move by one}), and that on the two-dimensional CPN−1CP^{N-1} test bed the freezing of the topological charge splits with NN while the ratios of local quantities did not (conditionally).

Prerequisite chapters — 11-01 (the pattern, the two accounts of (4)–(7), row (c) of Table 1), 10-01 §05 ((13) for Assumption H and the note on topological sectors), 10-01 §07 (the bookkeeping length ξRG\xi_{\rm RG} and the logarithmic factor), the detour in 2-10 (conserving dynamics), "Differences other than constants" in 2-10 (the choice of variables). The notation is collected in 11-A1.

Said first — The level is [Y3]. The numbers in this chapter are numbers for the Gauss approximation on a hierarchical lattice and for two-dimensional models (the U(1)U(1) strip, CPN−1CP^{N-1}). Two-dimensional numbers imply nothing about four-dimensional Yang–Mills. "Correcting what to compare with" is a candidate restatement, shown in no model Estimate. This chapter shows neither the mass gap nor whether Assumption H is true.

The order of this chapter
  1. Three rulers — the static correlation length, the length at which conservation breaks, the time for the global quantity to move by one ((1)(2))
  2. Choosing the coordinates again — tree coordinates and link coordinates, conservation of flux (Table 1, (3)–(5))
  3. The two-dimensional U(1)U(1) strip — the length at which conservation breaks at five points (Table 2, (6))
  4. Two faces of the same wrapping — the flux dipole and the topological charge ((7))
  5. Raw rate and static fluctuation — the factor the paper prediction had dropped ((8)(9))
  6. Correcting what to compare with — max⁡(ξ2,τglob)\max(\xi^2,\tau_{\rm glob}) ((10))
  7. The two-dimensional CPN−1CP^{N-1} test bed — freezing splits with NN, local ratios do not (Table 3, Table 4)
  8. What cannot be said
  9. Table of the five grades and the dependency graph — Table 5

01

Three rulers — what length is the ξ\xi of H?

Assumption H is the power statement lim⁡ln⁡CP/ln⁡ξ=2\lim\ln C_P/\ln\xi=2 ((13) of 10-01), comparing the relaxation time CPC_P with the square of one length ξ\xi. Part 10 took this ξ\xi to be the correlation length 1/mlat1/m_{\rm lat} and used the bookkeeping length ξRG\xi_{\rm RG} in its place (10-01 §07). But in the pattern of 11-01, what fixes the slowness is not the correlation length but "the scale of the events that break the global quantity QQ". So we split the candidates for what to compare with into three.

(a) The static correlation length ξ\xi. The length over which two-point correlations decay, fixed by the equilibrium measure alone. In four dimensions ξ=1/mlat\xi=1/m_{\rm lat}. In the two-dimensional U(1)U(1) strip, one row of the strip becomes a one-dimensional chain of rotors in axial gauge, and the length over which Wilson loops decay with area (the inverse string tension),

ξXY=−1ln⁡(I1(β)/I0(β)) ≈ 2β(β large)(1)\xi_{XY}=-\frac{1}{\ln\bigl(I_1(\beta)/I_0(\beta)\bigr)}\ \approx\ 2\beta\qquad(\beta\ \text{large}) \tag{1}

plays this role (IνI_\nu is a modified Bessel function. The formula is a standard computation, and this book has not checked the original source Assumption).

(b) The length at which conservation breaks ℓ∗\ell_*. The global quantity (charge) QRQ_R in a subregion RR changes by two routes: the route of transport, flowing out through the boundary (taking time of order R2R^2), and the route of conservation breaking inside the region (dislocations, wrapping). The length at which the two balance is ℓ∗\ell_*. The crossing length R∗R_* of (7) in 11-01 was its static version. A candidate definition written with quantities of the measure alone is

S:=sup⁡RVar⁡(QR)E(QR),ℓ∗:=inf⁡{R: Var⁡(QR)E(QR) ≥ S2}(S=∞ implies ℓ∗=∞)(2)S:=\sup_R\frac{\operatorname{Var}(Q_R)}{\mathcal E(Q_R)},\qquad \ell_*:=\inf\Bigl\{R:\ \frac{\operatorname{Var}(Q_R)}{\mathcal E(Q_R)}\ \ge\ \frac S2\Bigr\}\quad(S=\infty\ \text{implies}\ \ell_*=\infty) \tag{2}

Estimate. Putting in the estimate Var⁡/E≈χtR2/(hR2+dλW)\operatorname{Var}/\mathcal E\approx\chi_tR^2/(hR^2+d\lambda_W) of (7) in 11-01 gives S=χt/hS=\chi_t/h, and the ratio reaches half of it when hR2=dλWhR^2=d\lambda_W, that is, ℓ∗=R∗\ell_*=R_* Assumption (on top of the estimate (7)). If conservation never breaks (h=0h=0), then S=∞S=\infty and ℓ∗=∞\ell_*=\infty, which agrees with the flux of the layers in the Gauss approximation of §02 below (where conservation holds exactly). Strictly, a definition of the form "the supremum of RR with Var⁡/E≤KR2\operatorname{Var}/\mathcal E\le KR^2", comparing the ratio with R2R^2, does not fix a length: under the estimate (7) the ratio divided by R2R^2 decreases in RR, so the set of RR satisfying the condition is a half-line open upwards or empty, and the supremum is infinite or undefined. (2) avoids this by cutting at half the saturation value.

(c) The time for the global quantity to move by one τglob\tau_{\rm glob}. A time, not a length. Its representative is the integrated autocorrelation time τint(Q)\tau_{\rm int}(Q) of the topological charge QQ in the periodic box, and for reversible dynamics τint(Q)≤CP\tau_{\rm int}(Q)\le C_P ((1) of 11-01).

The three can take different values within the same model. In the 0+1-dimensional double well, the slow quantity (the number of kinks) and the correlation length are born from the same kinks, and the three differ only by logarithms (11-01 Table 3). In the two-dimensional strip of §03, (a) grows linearly in β\beta while (b)(c) grow exponentially, and the three move apart. Fixing the model does not automatically fix which ruler H is a sentence about.

What this section can say: (2) is a candidate definition, and the four-dimensional value of ℓ∗\ell_* is not in this book. The ξXY\xi_{XY} of (1) is used only as the static length of the strip; it does not stand in for the four-dimensional correlation length.

02

Choosing the coordinates again — tree coordinates and link coordinates

The denominator of the Poincaré constant CP=sup⁡Var⁡(f)/E(f)C_P=\sup\operatorname{Var}(f)/\mathcal E(f) is a Dirichlet form, and it is not fixed until one decides the variables whose gradients it measures. We saw in Part 2 that changing variables by a constant factor changes CPC_P by a constant factor ("Differences other than constants" in 2-10: a factor β\beta depending on whether one measures by link angles or by the field). In this section we give an example where the choice of coordinates changes the type of growth of CPC_P.

Stage. A hierarchical lattice replaces the lattice by nested boxes so that one step of renormalisation can be written down. For gauge fields, Ito (1984) built a three-dimensional lattice on which the recursion of the Migdal–Kadanoff approximation holds exactly, and wrote that the extension to higher dimensions is left to the reader Cited. This book defines for itself a four-dimensional version that coarsens by factors of two (each layer a two-dimensional periodic lattice of side L=2NL=2^N, with the bonds of level rr identified across 4r4^r layers), and treats only the Gauss approximation on it (U(1)U(1), or linearisation to one component of the Lie algebra of SU(2)SU(2)) Assumption (the four-dimensional construction is this book's extension).

Tree coordinates. Taking the values of the plaquettes (the holonomies at the nodes) as variables, in the Gauss approximation the configuration space is a linear subspace SS, the equilibrium measure is the density e−β∣x∣2/2e^{-\beta|x|^2/2} on SS, and the Hessian of the Langevin potential is βI\beta I on SS.

tree coordinates (Gauss approximation):λ1=β,CP=1β(every number of levels N)(3)\text{tree coordinates (Gauss approximation)}:\qquad \lambda_1=\beta,\qquad C_P=\frac1\beta\qquad(\text{every number of levels}\ N) \tag{3}

Proved here Proof. It is an Ornstein–Uhlenbeck process with unit metric and drift matrix βI\beta I, so the gap of the generator is β\beta. ∎ Lean HierLink.tree_no_decay (bundle LatticeGaugeHierarchicalLink; assumption: 0<β0\lt\beta; only the arithmetic that the ratio is 1). Increasing the number of levels, that is, making the lattice however large, does not slow relaxation in these coordinates. Even if the correlation length grows under the non-abelian flow, the coarse levels are expected to stay close to Haar with CP=O(1)C_P=O(1) Estimate, and in tree coordinates H is false or empty.

Link coordinates. What H in Part 10 names is Langevin in the link variables (10-01 §04). The link generator in the Gauss approximation is β DΠD∗\beta\,D\Pi D^* (DD is the exterior derivative from links to plaquettes, Π\Pi the projection of the identification), and its smallest positive eigenvalue is

λmin⁡(N)=4sin⁡2π2N+1,4Nλmin⁡(N) ⟶ π2(N→∞)(4)\lambda_{\min}(N)=4\sin^2\frac{\pi}{2^{N+1}},\qquad 4^N\lambda_{\min}(N)\ \longrightarrow\ \pi^2\qquad(N\to\infty) \tag{4}

Proved here (paper: splitting functions of the layers into sectors by the hierarchical projection, the block operator for sector kk becomes the Neumann Laplacian of the two-dimensional dual lattice of side 2k+12^{k+1}, and the minimum is the value for sector k=N−1k=N-1 with side 2N2^N) Reading (agreeing with the direct full spectrum for N≤2N\le2, and with a separate implementation by sectors for N=3,4N=3,4). The gap of the generator is βλmin⁡\beta\lambda_{\min}, which shrinks inversely with the square of the side. Lean HierLink.lam_le (no assumptions: 4sin⁡2(π/2N+1)≤π2/4N4\sin^2(\pi/2^{N+1})\le\pi^2/4^N; an upper bound in the direction of slowness; that λmin⁡\lambda_{\min} is equal to this expression is outside Lean), HierLink.measured_scaled (with the interval of the measured λmin⁡(4)\lambda_{\min}(4) as assumption, 44λmin⁡∈(9.83,9.84)4^4\lambda_{\min}\in(9.83,9.84)).

Table 1 Coordinates and gap (Gauss approximation, one component of U(1)U(1), this book's four-dimensional extension; the gap is divided by β\beta)
NNSide L=2NL=2^NTree coordinatesLink coordinates (metric of the original lattice)4Nλmin⁡4^N\lambda_{\min}Link coordinates (metric per class)Gauge quotient on the four-dimensional original latticeControl with layers decoupled 4sin⁡2(π/L)4\sin^2(\pi/L)
1212.0000008.00002.0000002.00004
2410.58578649.37260.87689440.58578642
3810.15224099.74340.2575235—0.5858
41610.03842949.83790.0667663—0.1522

Labels of Table 1: the tree-coordinate column is Proved here by (3). Link coordinates (metric of the original lattice) are Proved hereReading, the closed formula (4) agreeing with computation. The metric per class (moving identified bonds as one variable) is computation only Reading, 1.0–1.74 times the metric of the original lattice. The gauge quotient on the four-dimensional original lattice agrees with the metric of the original lattice to 10−910^{-9} for N≤2N\le2 Reading. The lowest eigenvector lies, for all NN, in the same sector (the contrast between the four top-level groups of layers) Reading.

The slowness comes from conservation of flux. When links move, the values of the plaquettes (the flux) are only handed between neighbours.

ddtF=d(ddtA),∑p∈blockFp=0(when the boundary bonds are fixed; Stokes)(5)\frac{d}{dt}F=d\Bigl(\frac{d}{dt}A\Bigr),\qquad \sum_{p\in\text{block}}F_p=0\quad(\text{when the boundary bonds are fixed; Stokes}) \tag{5}

Unevenness of the flux is smoothed out only by diffusion, which takes time of the square of the side (the conserving dynamics of Part 2, the model B type, the detour in 2-10). The decisive point is the last column of Table 1. If the layers are decoupled and every bond made an independent variable per layer, then statically the plaquettes are almost independent (only the constraint on the total) and the correlation length is 0, yet the smallest eigenvalue is 4sin⁡2(π/L)4\sin^2(\pi/L) and the same L−2L^{-2} type remains Proved here. Hence this slowness does not come from the correlation length. The static correlation of plaquettes in the Gauss approximation decays masslessly, of order −β−1R−4-\beta^{-1}R^{-4}; β\beta does not flow (the Gauss-approximation face of the coupling being dimensionless in four dimensions, 11-07), and the correlation length is only cut off by the volume Proved here.

Consequence. Assumption H is not a sentence about the measure alone but a sentence about the pair (measure, Dirichlet form) Estimate. Tree coordinates and link coordinates have the same equilibrium measure, yet in the former CPC_P does not depend on the number of levels, and in the latter it grows as the square of the side. In Part 2, "CPC_P changes by a constant factor with the coordinates" did not change the power of the growth, but here the power itself (N0N^0 versus 4N4^N) changes. H in Part 10 names Langevin in the link variables, so it already meets this point. Reinterpreting heat-bath or overrelaxation numbers as CPC_P of link Langevin is an assumption (2-10 §"Do not identify heat bath with Langevin"). And even in link coordinates, conservation holds exactly in the Gauss approximation, so ℓ∗\ell_* of (2) is infinite, and comparing "CP∝L2C_P\propto L^2" with the correlation length has no meaning. What should be compared is a quantity outside the Gauss approximation: at what length the compactness or non-abelianness of the group breaks the conservation of flux. The next section is the first measurement of it.

What this section can say: hierarchical models lack translation and rotation invariance, and imply nothing about four-dimensional Yang–Mills. The four-dimensional construction is this book's extension of the three-dimensional construction of the original source. The numbers are numbers for the Gauss approximation, with N≤4N\le4.

03

The two-dimensional U(1)U(1) strip — the length at which conservation breaks at five points

Model. We cut one layer of the previous section into a narrow strip. Plaquettes (x,y)(x,y), x=0,…,M−1x=0,\dots,M-1, y=0,1y=0,1 (width W=2W=2); all links on the outer boundary are fixed at angle 0, and the inner links are moved by the Langevin dynamics dθℓ=−∂ℓS dt+2 dBℓd\theta_\ell=-\partial_\ell S\,dt+\sqrt2\,dB_\ell for the action S=−β∑pcos⁡φpS=-\beta\sum_p\cos\varphi_p (unit metric). Gauge transformations at inner vertices do not change the plaquette angles φp\varphi_p, so on gauge-invariant functions the link generator reduces to a diffusion on φ\varphi, whose metric is the Neumann Laplacian of the dual graph of the plaquettes Proved here. The equilibrium measure is the product of independent eβcos⁡φe^{\beta\cos\varphi} per plaquette conditioned on ∑pφp≡0 (mod 2π)\sum_p\varphi_p\equiv0\ (\mathrm{mod}\ 2\pi), and there is almost no static correlation between plaquettes.

In the Gauss approximation conservation holds exactly. With sin⁡φ→φ\sin\varphi\to\varphi, the smallest gap is β⋅4sin⁡2(π/2M)≈βπ2/M2\beta\cdot4\sin^2(\pi/2M)\approx\beta\pi^2/M^2, inversely proportional to the square of the strip length Proved here (Lean ConsLen.neumann_gap_lower, ConsLen.neumann_gap_upper, bundle LatticeGaugeConservationLength; assumption: 1≤M1\le M; 4/M2≤4sin⁡2(π/2M)≤π2/M24/M^2\le4\sin^2(\pi/2M)\le\pi^2/M^2. That the gap equals this expression is outside Lean).

If the group is compact, conservation breaks. An event in which the angle of a single plaquette crosses ±π\pm\pi (a wrapping) changes by one the sum of the angles reduced to (−π,π](-\pi,\pi],

Q:=12π∑p[φp] ∈ Z(6)Q:=\frac1{2\pi}\sum_p[\varphi_p]\ \in\ \mathbb Z \tag{6}

(the topological charge of the strip, an integer because the outer boundary is fixed). A wrapping is a rare event crossing a barrier of height 2β2\beta. We measure two quantities: λX:=1/τint(X)\lambda_X:=1/\tau_{\rm int}(X) for the flux dipole X=∑pwp[φp]X=\sum_pw_p[\varphi_p] (wp=cos⁡(π(x+12)/M)w_p=\cos(\pi(x+\tfrac12)/M), the slowest mode of the Gauss approximation), and λQ:=1/τint(Q)\lambda_Q:=1/\tau_{\rm int}(Q) for the topological charge. For reversible dynamics τint≤CP\tau_{\rm int}\le C_P ((1) of 11-01), so both are read as upper bounds on the gap. The length at which conservation breaks is defined by fitting λX=A/M2+B\lambda_X=A/M^2+B and setting ℓc:=A/B\ell_c:=\sqrt{A/B} (the length at which the diffusion rate A/M2A/M^2 equals the saturation value BB).

Table 2 Five points for the two-dimensional U(1)U(1) strip (width 2, outer boundary frozen, gauge quotient of link Langevin, strip length M≤256M\le256; λQ\lambda_Q is a weighted mean over MM with Var⁡Q≥0.5\operatorname{Var}Q\ge0.5; E:=E [φp]2E:=\mathbb E\,[\varphi_p]^2; errors by Madras–Sokal and the covariance of the fit)
β\betaξXY\xi_{XY}ℓc\ell_cPaper estimateλQ\lambda_Qτ(Q)\tau(Q)π/EλQ\pi/\sqrt{E\lambda_Q}ρ=ℓc/πEλQ\rho=\ell_c\big/\tfrac{\pi}{\sqrt{E\lambda_Q}}ℓc/ξXY\ell_c/\xi_{XY}
22.783.29(20)3.530.769(21)1.304.100.80(5)1.2
34.759.88(39)8.280.224(13)4.4610.00.985(50)2.1
46.8224.7(9)19.60.0579(29)17.323.91.03(5)3.6
58.8753.6(3.1)47.70.0127(10)7958.50.92(6)6.0
610.90139(14)1180.00276(32)3621391.00(11)12.8

Every label in Table 2 is Reading. "Paper estimate" is the value, by numerical integration, of a capacity estimate with a test function (the length at which the two terms of an upper bound become equal, using a test function whose return from wrapping is optimised by capacity; for large β\beta, eβπ/(6β)e^\beta\sqrt{\pi/(6\beta)}) Assumption (an estimate on top of the independence approximation), and is of the same order as the measured ℓc\ell_c.

Slopes. With a weighted straight line over β=2,…,6\beta=2,\dots,6:

Reading Lean ConsExp.lc_slope_2_6 (bundle LatticeGaugeConservationExponent; assumptions: ℓc(2)∈[3.09,3.49]\ell_c(2)\in[3.09,3.49], ℓc(6)∈[125,153]\ell_c(6)\in[125,153], the measured 1σ1\sigma intervals; conclusion: a∈[3/4,1]a\in[3/4,1]), ConsExp.lamQ_slope_2_6 (with the 1σ1\sigma intervals of λQ(2)\lambda_Q(2) and λQ(6)\lambda_Q(6) as assumptions, b∈[5/4,3/2]b\in[5/4,3/2]). The bundles guarantee only interval arithmetic.

Not proportional to the static length. The static length ξXY≈2β\xi_{XY}\approx2\beta is linear in β\beta and grows only about fourfold over the five points. The length at which conservation breaks, ℓc∝e0.92β\ell_c\propto e^{0.92\beta}, grows 42-fold, and the ratio ℓc/ξXY\ell_c/\xi_{XY} opens from 1.21.2 to 12.812.8 Reading. Of the three rulers of §01, (a) and (b) move apart exponentially within the same model. When the strip is longer than ℓc\ell_c, relaxation stops depending on the strip length and is fixed by the time for QQ to move by one.

Conditions: two-dimensional U(1)U(1), width 2, outer boundary frozen, gauge quotient of link Langevin with unit metric, β=2,…,6\beta=2,\dots,6, strip length M≤256M\le256, τint\tau_{\rm int} as a reading of an upper bound on the gap, fit of the form A/M2+BA/M^2+B (for β=2\beta=2, χ2/dof=3.3\chi^2/\text{dof}=3.3, and the form does not fit). Nothing is said about the exponent as β→∞\beta\to\infty.

04

Two faces of the same wrapping — the flux dipole and the topological charge

The ρ\rho column of Table 2 shows that the two quantities are one process. The relaxation of the flux dipole XX falls as M−2M^{-2} at the rate μ1(M)/E\mu_1(M)/E (μ1≈π2/M2\mu_1\approx\pi^2/M^2), with the Gauss stiffness replaced by 1/E1/E, while the relaxation rate λQ\lambda_Q of the topological charge QQ does not depend on the strip length. The length at which the time for diffusion to cross a length ℓ\ell equals the time for QQ to move by one is

π2E ℓ2=λQ⟺ℓ=πE λQ(7)\frac{\pi^2}{E\,\ell^2}=\lambda_Q\quad\Longleftrightarrow\quad \ell=\frac{\pi}{\sqrt{E\,\lambda_Q}} \tag{7}

Proved here (algebra) Lean ConsExp.lc_eq_iff (assumptions: 0<E0\lt E, 0<λ0\lt\lambda, 0<ℓ0\lt\ell), and its ratio ρ\rho to the measured ℓc\ell_c is 0.800.80 to 1.031.03 over the five points Reading. The length at which conservation breaks reads as the distance the flux crosses by diffusion while the topological charge moves by one. The flux dipole departing from M−2M^{-2} and the topological charge relaxing are two faces of the same wrapping.

A single point at β=8\beta=8 supplements this reading. Up to strip length 64, XX keeps M−2M^{-2} (ratio to the Gauss stiffness 0.860.86 to 1.031.03), but at strip length 32, τint(Q)=5.3(1.7)×103\tau_{\rm int}(Q)=5.3(1.7)\times10^3 (run length 10610^6, Var⁡Q=0.145\operatorname{Var}Q=0.145), and a slow mode at about 1/3801/380 of the value for XX appears Reading. Because Var⁡Q\operatorname{Var}Q is small (a sector where Q≠0Q\ne0 is rare), this point is not included in the slopes.

What this section can say: ℓc\ell_c is the length of the bend in the relaxation of XX, not the length of a bend in the gap itself (including QQ, the upper bound on the gap at β=4\beta=4 already saturates at M≈8M\approx8 to 1616). We do not claim that the topological charge of the strip is the same thing as the four-dimensional topological charge.

05

Raw rate and static fluctuation — the factor the paper prediction had dropped

The paper prediction before measuring was "λQ∝βe−2β\lambda_Q\propto\beta e^{-2\beta}" (the Kramers rate for the barrier 2β2\beta), so the slope bb should have been near 2. The measured b=1.355(18)b=1.355(18) is gentler, and moreover split before and after β\beta (bottom of Table 2). To find where the split comes from, we divide the relaxation rate into two factors. QQ can be regarded as a chain of ±1\pm1 jumps with a restoring force, so from the total jump rate WMνWM\nu (ν\nu is the rate of wrapping per plaquette) and the variance Var⁡Q≈WM E/(4π2)\operatorname{Var}Q\approx WM\,E/(4\pi^2) comes

λQ=κ⋅2π2 νE,E=E [φp]2(8)\lambda_Q=\kappa\cdot\frac{2\pi^2\,\nu}{E},\qquad E=\mathbb E\,[\varphi_p]^2 \tag{8}

this form Assumption (the chain picture; κ\kappa is a correction for the raw count. It is a different quantity from the comparison constant κ\kappa of Proposition B (11-A1)). The ν\nu in the numerator is the raw rate (a quantity of the dynamics), and the EE in the denominator is the static fluctuation (a quantity of the measure). Taking logarithms gives

Δln⁡λQ=Δln⁡ν−Δln⁡E(9)\Delta\ln\lambda_Q=\Delta\ln\nu-\Delta\ln E \tag{9}

Proved here (algebra) Lean ConsExp.lamQ_log_split (assumptions: 0<κ0\lt\kappa, each quantity positive).

Counting the raw number of wrappings νraw\nu_{\rm raw} (the number of times ±π\pm\pi is crossed), we found the following Reading.

From (9), b=(slope of the raw rate)−(ln⁡(1/E) slope)b=(\text{slope of the raw rate})-(\ln(1/E)\ \text{slope}), and since the latter decreases with β\beta, bb increases with β\beta. This is what the split really is. The paper prediction had dropped the static factor 1/E1/E. The corrected form is λQ∝βe−2β/E\lambda_Q\propto\beta e^{-2\beta}/E (for large β\beta, ∝β2e−2β\propto\beta^2e^{-2\beta}, b→2−2/βb\to2-2/\beta), in the same direction as the measured bb Estimate. In the combination π/EλQ\pi/\sqrt{E\lambda_Q}, EE cancels and only ν−1/2\nu^{-1/2} remains, so the slope is half the slope of the raw rate (about 0.880.88) and constant over all intervals. That is the third slope at the bottom of Table 2.

Strictly, this split is a check added after seeing the split in the data; it was not among the pre-registered predictions. Also, nothing can be said from these five points about the exponent as β→∞\beta\to\infty (whether b→2b\to2).

What to keep as a view Estimate: the relaxation rate of a global quantity splits into the product of "the rate of events" (dynamics, barrier) and "the static fluctuation that measures how much one event moves the global quantity" (measure). The one-loop accounting that the original sources write for the four-dimensional topological charge (11-07 §02) also counts only the exponent of the density of events. Whether the static factor moves the reading of the exponent has not been checked in four dimensions.

06

Correcting what to compare with — max⁡(ξ2,τglob)\max(\xi^2,\tau_{\rm glob})

In the strip, measuring τ(Q)\tau(Q) with the static length ξXY\xi_{XY} as the ruler, τ(Q)/ξXY2\tau(Q)/\xi_{XY}^2 is about 0.40.4 at β=4\beta=4 and about 2424 at β=8\beta=8 (strip length 32) Reading. Since CP≥τint(Q)C_P\ge\tau_{\rm int}(Q), CP/ξXY2C_P/\xi_{XY}^2 also grows with β\beta. "CP≲ξ2C_P\lesssim\xi^2", with only the static correlation length as the ruler, tends to fail with β\beta inside this toy (row (c) of the pattern in 11-01). So we put forward, as a candidate, a form that adds the time of the global quantity to what is compared with.

CP ≍ max⁡(ξ2, τglob),τglob:=(the time for the global quantity to move by one)(10)C_P\ \asymp\ \max\bigl(\xi^2,\ \tau_{\rm glob}\bigr),\qquad \tau_{\rm glob}:=\bigl(\text{the time for the global quantity to move by one}\bigr) \tag{10}

Estimate. This is one direction for correcting H. The other direction is to restrict the side of test functions to families local in time and space (11-03 §02), where, as long as the box is smaller than ℓ∗\ell_*, the time of the global quantity is expected not to enter (§07).

(10) has three weaknesses. (i) If τglob\tau_{\rm glob} is defined by τint(Q)\tau_{\rm int}(Q) of the same dynamics, (10) is close to a restatement of τint(Q)≤CP\tau_{\rm int}(Q)\le C_P and has no content. To give it content, τglob\tau_{\rm glob} must be fixed independently of the dynamics. Candidates are the split of §05 (barrier rate and static fluctuation, τglob≈E/(2π2κν)\tau_{\rm glob}\approx E/(2\pi^2\kappa\nu)), or making the diffusion time ℓ∗2/D\ell_*^2/D from the static ℓ∗\ell_* of (2) Estimate. (ii) The strip numbers see only the lower side of (10), CP≥τglobC_P\ge\tau_{\rm glob}, and not a single number bounds CPC_P from above. At β=4\beta=4, τ(Q)≈17\tau(Q)\approx17 against ξXY2≈46\xi_{XY}^2\approx46, and whether CPC_P is on the side of ξXY2\xi_{XY}^2 has not been measured. (iii) In four dimensions the candidate for τglob\tau_{\rm glob} is the topological charge in the periodic box, whose measurements are the heat-bath and HMC numbers of the original sources, and reinterpreting them for Langevin is an assumption (11-01 §06).

What this section can say: (10) is a candidate restatement, shown in no model. Nor do we propose replacing H of Part 10 with (10). The order of the account is: because H was set up sharply, we could name the places where a single thing to compare with is not enough.

07

The two-dimensional CPN−1CP^{N-1} test bed — freezing splits with NN, local ratios do not

Why this model. The two-dimensional CPN−1CP^{N-1} model shares two points with four-dimensional Yang–Mills — a dimensionless coupling, and instantons with no scale (11-07) — and lets one vary NN, which plays the role of the size of the group. There are primary sources reporting that in the periodic box the autocorrelation of the topological mode grows like exp⁡(1.72 ξ)\exp(1.72\,\xi) for ξ=2\xi=2 to 66 (N=21N=21, Hasenbusch 2017), and worsens with NN (Del Debbio–Manca–Vicari 2004, Berni–Bonanno–D'Elia 2019) Cited. It is also reported that with open boundaries the growth returns to z≈1z\approx1 Cited.

Setting. The action is Hasenbusch's form S=−βN∑x,μ(zˉx+μzxλx,μ+c.c.−2)S=-\beta N\sum_{x,\mu}(\bar z_{x+\mu}z_x\lambda_{x,\mu}+{\rm c.c.}-2) (zx∈CNz_x\in\mathbb C^N, ∣zx∣=1|z_x|=1, λ∈U(1)\lambda\in U(1)), and the topological charge is the sum of plaquette angles Q=12π∑θplaqQ=\frac1{2\pi}\sum\theta_{\rm plaq}. The dynamics is heat bath (sites and links) mixed with novn_{\rm ov} overrelaxation steps; it is not Langevin. The autocorrelations measured here are in units of heat-bath sweeps, and reinterpreting them as the Langevin CPC_P is an assumption (2-10). The local quantities are three for a box of side RR: the energy EBE_B, the topological charge of the box qBq_B, and the "magnetisation" of the box mBm_B.

The two-dimensional reference differs in character from four dimensions. The reference of the ratio form KrelK_{\rm rel} of 11-03 (the box value for the massless free field) diverges in two dimensions not only with the box size RR but logarithmically in the lattice side LL (the constant part of the continuum kernel −12πln⁡∣z∣-\tfrac1{2\pi}\ln|z| sits on the box average) Proved here. Removing the box average makes it finite, about 0.05480.0548 by extrapolation in 1/R21/R^2 Reading. Moreover, even for mean-zero quantities the free-field Langevin autocorrelation grows with ln⁡L\ln L Reading (in a box with R=16R=16, from 36 to 119 as the lattice side goes from 64 to 1024). So the ratios of this test bed are ratios of autocorrelations with the lattice side fixed at L=64L=64, Krel(c;qty):=τ(qtyB(cξ))/τfree(DB(cξ))K_{\rm rel}(c;\text{qty}):=\tau(\text{qty}_B(c\xi))/\tau^{\rm free}(D_B(c\xi)) (DBD_B is the first Neumann mode of the free-field box, with the same algorithm and the same LL) and Kglob:=τ(Q2)/τfree(DB(4ξ))K_{\rm glob}:=\tau(Q^2)/\tau^{\rm free}(D_B(4\xi)), quantities different from the static ratios of 11-03.

Positive control. With two runs (run length 70000 sweeps) at N=21N=21, β=0.70\beta=0.70, L=44L=44, nov=6n_{\rm ov}=6, it passed with τint(Q2)=160(26)\tau_{\rm int}(Q^2)=160(26) (0.80 times 199(18)199(18) of Hasenbusch's Table II, inside the band [139,259][139,259] fixed in advance), an energy difference of −9×10−6-9\times10^{-6}, and a ξ\xi difference of −0.1%-0.1\% Reading.

Table 3 Core pairs and pure heat-bath pairs (L=64L=64, periodic, ξ≈4\xi\approx4; KK is three values at c=1,2,4c=1,2,4; one run each, some runs shortened from the registered length)
NNβ\betanovn_{\rm ov}ξ2nd\xi_{2\rm nd}τ(Q2)\tau(Q^2)Changes of QQ / measurementsKglobK_{\rm glob}Krel(EB)K_{\rm rel}(E_B)Krel(mB)K_{\rm rel}(m_B)Krel(qB)K_{\rm rel}(q_B)
210.7063.738(14)154(46)527 / 30600Undefined (the reference τfree\tau^{\rm free} is 0.33–0.54)
21.14064.25(8)0.5093123 / 99000Same as above
210.7003.76(6)3874(2810)120 / 41400550.32, 0.11, 0.0380.68, 0.59, 0.511.71, 1.74, 2.27
100.8004.62(8)70.9(10.3)24347 / 810001.010.35, 0.13, 0.0490.83, 0.78, 0.571.51, 1.11, 0.83
21.14004.89(24)0.73133493 / 1440000.0100.26, 0.084, 0.0300.66, 0.41, 0.250.21, 0.062, 0.022

Every label in Table 3 is Reading. τ(Q2)\tau(Q^2) for N=21N=21, nov=0n_{\rm ov}=0 has a run length of only about 11τ11\tau and is treated as a lower bound.

What can be read.

  1. The freezing of the topological charge splits with NN. At ξ≈4\xi\approx4, τ(Q2)\tau(Q^2) for N=21N=21 against N=2N=2 is 154154 against 0.500.50 with overrelaxation (about 300 times), and 38743874 against 0.730.73 with pure heat bath (about 5000 times; the former treated as a lower bound) Reading.
  2. The ratios of local quantities are small regardless of NN. For the pure heat-bath pairs, KrelK_{\rm rel} for the energy and the magnetisation is almost the same at N=2,10,21N=2,10,21 (0.03–0.35 for EBE_B, 0.25–0.83 for mBm_B), and the ratios of one-step Rayleigh quotients (the ratio of 1/(1−ρ(1))1/(1-\rho(1)) to its free-field value) are all 1.1 or below Reading. The verdict at the freezing point was on the side of "freezing of topology does not break the local ratios up to box 4ξ4\xi". Conditions: L=64L=64, a single point ξ≈4\xi\approx4, one run each, a post-hoc cell (below). There is only one freezing point, (N=21, ξ≈4)(N=21,\ \xi\approx4), and Kglob=55K_{\rm glob}=55 has a 1σ1\sigma width of 15–95 from the error of one run, so whether it exceeds the decision threshold 30 is within that error. It is a weak reading at one point.
  3. ℓ∗\ell_* grows with NN. Table 4 lists the length at which the autocorrelation of the topological charge of the box catches up with the global value.
Table 4 Autocorrelation of the topological charge of the box τ(qB(R))\tau(q_B(R)) and the global τ(Q)\tau(Q) (pure heat bath, L=64L=64)
NNR=2R=2448816163232τ(Q)\tau(Q)τ(qB(32))/τ(Q)\tau(q_B(32))/\tau(Q)Exponent for R→2RR\to2R
213.089.8238.916073840710.181.67, 1.98, 2.04, 2.21
103.018.7224.758.41132220.511.54, 1.50, 1.24, 0.95
20.981.211.391.511.531.550.990.30, 0.20, 0.12, 0.02

Reading For N=2N=2 the autocorrelation of the box charge saturates at the global value at R≈8R\approx8 to 1616; for N=10N=10 it begins to bend at R=32R=32; for N=21N=21 it still grows as R2R^2 at R=32R=32. The box charge takes time R2R^2 through transport across the boundary (the route (b) of §01), and the length at which it catches up with the global value is the reading of ℓ∗\ell_*. In this model and this range, freezing does not break the ratios for boxes smaller than ℓ∗\ell_*, and rides only on quantities on scales beyond ℓ∗\ell_* Estimate. That in Table 3 only Krel(qB)K_{\rm rel}(q_B) of the box charge splits with NN (0.022→0.83→2.270.022\to0.83\to2.27 at c=4c=4) we read as a reflection of this difference in catching-up length. In the language of 11-03, the local H asks about the transport coefficient dd and the global H about the coefficient of breaking of conservation hh (11-03 §05), so this is one two-dimensional face of the prospect that local ratios are unaffected by freezing. It also agrees in direction with the conclusion of 11-04 that what the one-sided bridge needs is the local side (11-04), but the ratios here are ratios of autocorrelations, not the ratio of the bridge's threshold (a static ratio).

What should be disclosed.

What this section can say: numbers for the two-dimensional CPN−1CP^{N-1}, this action and these algorithms (heat bath and overrelaxation), L=64L=64, a single point ξ≈4\xi\approx4; nothing is said about four dimensions, SU(N)SU(N), Assumption H or the mass gap.

08

What cannot be said

(a) Two dimensions and four. The numbers of §03–§07 are numbers for the two-dimensional U(1)U(1) strip and the two-dimensional CPN−1CP^{N-1}. Two-dimensional U(1)U(1) is a trivial theory in the continuum limit, and neither set of numbers implies anything about exponents in four dimensions, non-abelian, β→∞\beta\to\infty.

(b) Readings of upper bounds. λX\lambda_X and λQ\lambda_Q are readings as upper bounds on the gap (lower bounds on CPC_P), not the gap itself. ℓc\ell_c is the length of the bend in the relaxation of the flux dipole and depends on the form of the fit A/M2+BA/M^2+B.

(c) Candidates. The ℓ∗\ell_* of (2) and the thing to compare with of (10) are candidates, shown in no model. The strip numbers see only the lower side of (10).

(d) The hierarchical lattice. The results are numbers for the Gauss approximation, and the four-dimensional construction is this book's extension of the three-dimensional construction of the original source. Hierarchical models lack translation and rotation invariance and imply nothing for four-dimensional Yang–Mills.

(e) Conditions of the test bed. The CPN−1CP^{N-1} ratios are ratios of autocorrelations of heat-bath dynamics, conditional quantities with L=64L=64 fixed, and different from the static ratios of 11-03. For each NN there is one point in ξ\xi, one freezing point, and the verdict on local ratios is a weak reading at one point on a cell added after registration; the registered verdict could not be defined. zQ(N)z_Q(N) has not been read.

(f) Afterthought. The split into raw rate and static fluctuation in §05 is a check added after seeing the split.

(g) The analogy with SU(2)SU(2). Using the split of freezing with NN as a two-dimensional analogy of the borderline reading for SU(2)SU(2) is an estimate, and N=2N=2 carries an ultraviolet label (11-07).

09

Table of the five grades and the dependency graph

Table 5 Rows of this chapter ("—" marks rows without machine checking)
RowStatement (one line)GradeGroundsScopeMachine checking
[01]The Gauss approximation (abelian linearisation), and the four-dimensional extension of the hierarchical latticeAssumption§02N≤4N\le4—
[02]λ1=β\lambda_1=\beta in tree coordinates ((3))Proved here§02Gauss approximation, all NNHierLink.tree_no_decay (assumption: 0<β0\lt\beta; arithmetic of the ratio)
[03]λmin⁡(N)=4sin⁡2(π/2N+1)\lambda_{\min}(N)=4\sin^2(\pi/2^{N+1}) in link coordinates ((4))Proved hereSector decomposition of §02Gauss approximationHierLink.lam_le (no assumptions; only the direction of the upper bound)
[04]Cross-check for N≤4N\le4, metric per class, gauge quotient on the original latticeReadingTable 1N≤4N\le4 (gauge quotient N≤2N\le2)HierLink.measured_scaled (measured interval as assumption)
[05]L−2L^{-2} remains even with the layers decoupled (conservation of flux)Proved here(5), Table 1Gauss approximation—
[06]H is a sentence about the pair (measure, Dirichlet form)Estimate§02——
[07]Gauss gap of the strip β⋅4sin⁡2(π/2M)\beta\cdot4\sin^2(\pi/2M)Proved here§03Width 2, outer boundary frozenConsLen.neumann_gap_lower/upper (assumption: 1≤M1\le M)
[08]τint≤CP\tau_{\rm int}\le C_PProved here(1) of 11-01Reversible dynamicsLatticeGauge.FourDTrial.tau_int_le_CP (assumptions: weights ≥0\ge0, eigenvalues ≥λ1>0\ge\lambda_1\gt0)
[09]Five points of the strip: ℓc\ell_c, λQ\lambda_Q, slopes a,ba,bReadingTable 2Conditions of §03ConsExp.lc_slope_2_6, lamQ_slope_2_6 (1σ1\sigma intervals as assumptions)
[10]Algebra of (7)Proved here§04—ConsExp.lc_eq_iff (assumptions: 0<E,λ,ℓ0\lt E,\lambda,\ell)
[11]ρ=0.80\rho=0.80 to 1.031.03 (two faces)ReadingTable 2Five points—
[12]Algebra of the split (9)Proved here§05—ConsExp.lamQ_log_split (assumptions: 0<κ0\lt\kappa, positivity)
[13]The raw rate within 1–3% of the Kramers slope, κ\kappa constant (added afterwards)Reading§05β=3..6\beta=3..6—
[14]The paper prediction had dropped 1/E1/EEstimate§05Strip—
[15]τ(Q)/ξXY2\tau(Q)/\xi_{XY}^2 goes 0.4→240.4\to24Reading§06β=4,8\beta=4,8—
[16]Compare with max⁡(ξ2,τglob)\max(\xi^2,\tau_{\rm glob}) ((10))Estimate§06——
[17]Passing the positive control for CPN−1CP^{N-1}Reading§07N=21N=21, L=44L=44—
[18]τ(Q2)\tau(Q^2) splits with NNReadingTable 3L=64L=64, ξ≈4\xi\approx4—
[19]Local ratios are small regardless of NNReadingTable 3Post-hoc cell, one point, weak—
[20]The length at which the box charge catches up with the global one grows with NNReadingTable 4L=64L=64—
[21]N=2N=2 may fail to freeze because the ultraviolet fixes QQEstimate§07——
[22]Freezing rides only on quantities on scales beyond ℓ∗\ell_*Estimate§07This model, this range—

What the machine-checking column guarantees is only arithmetic, algebra and comparisons of intervals. That eigenvalues equal closed formulas, the spectral representation of τint≤CP\tau_{\rm int}\le C_P, the fits and the Monte Carlo lie outside Lean (10-13).

Dependency graph (premise numbers are smaller than conclusion numbers):

[01] ASM GA    Gauss approx., 4D extension of hier. lattice     N ≤ 4
[02] PRV TREE  λ₁ = β in tree coordinates ⟸ [01]                Lean tree_no_decay
[03] PRV LINK  λ_min(N) in link coordinates ⟸ [01]              Lean lam_le (upper bound)
[04] RDG CHK   cross-check for N ≤ 4 ⟸ [03]                     Lean measured_scaled
[05] PRV FLUX  L⁻² even with layers decoupled ⟸ [01]
[06] EST PAIR  H is a sentence about a pair (measure, form) ⟸ [02][03][05]
[07] PRV NG    Gauss gap of the strip                           Lean neumann_gap_lower/upper
[08] PRV TI    τ_int ≤ C_P ((1) of 11-01)                       Lean tau_int_le_CP
[09] RDG BAND  five points of the strip ℓ_c, λ_Q, a, b ⟸ [07][08]
               Lean lc_slope_2_6, lamQ_slope_2_6
[10] PRV LC    algebra of (7)                                   Lean lc_eq_iff
[11] RDG RHO   ρ = 0.80–1.03 ⟸ [09][10]
[12] PRV SPL   the split (9)                                    Lean lamQ_log_split
[13] RDG KR    raw rate and κ (added afterwards) ⟸ [09]
[14] EST EF    the paper prediction dropped 1/E ⟸ [12][13]
[15] RDG TQ    τ(Q)/ξ_XY² goes 0.4 → 24 ⟸ [08][09]
[16] EST MAX   compare with max(ξ², τ_glob) ⟸ [06][11][15]
[17] RDG PC    positive control for CP^{N−1}
[18] RDG SPLIT τ(Q²) splits with N ⟸ [17]
[19] RDG LOC   local ratios small regardless of N ⟸ [17]        post-hoc cell, one point
[20] RDG QB    catching-up length grows with N ⟸ [17]
[21] EST UV    is N = 2 on the side where UV fixes Q? ⟸ [18]
[22] EST ELL   freezing rides only on quantities beyond ℓ_* ⟸ [19][20]

How to check: see by eye that on every "⟸" the numbers on the right are smaller than the number on the left. There are no circles. The conclusion nodes of this chapter, [06][16][22], are all "Estimate", and there is neither a four-dimensional node nor a node for the mass gap. The two-dimensional readings [09]–[22] are premises of none of the four-dimensional rows.


What this chapter can say

SaidNot said
Proved here On the hierarchical lattice in the Gauss approximation, the gap in tree coordinates is β\beta and the gap in link coordinates is β⋅4sin⁡2(π/2N+1)\beta\cdot4\sin^2(\pi/2^{N+1}), and the slowness comes from conservation of fluxThe gap of non-linear, non-abelian hierarchical lattices. Implications for four-dimensional Yang–Mills
Reading The length at which conservation breaks, ℓc\ell_c, at five points of the two-dimensional U(1)U(1) strip with slope 0.916(21)0.916(21); that it is not proportional to the static length; ℓc≈π/EλQ\ell_c\approx\pi/\sqrt{E\lambda_Q}The exponent as β→∞\beta\to\infty. The bend of the gap itself. The four-dimensional ℓ∗\ell_*
Reading The raw rate follows the Kramers slope, and the split comes from the static factor 1/E1/E (added afterwards)The static factor in the rate of the four-dimensional topological charge
Estimate The three rulers, the candidate definition (2) of ℓ∗\ell_*, comparing with max⁡(ξ2,τglob)\max(\xi^2,\tau_{\rm glob})A sentence shown in any model. The right way to correct H
Reading In two-dimensional CPN−1CP^{N-1} (L=64L=64, ξ≈4\xi\approx4), freezing splits with NN, local ratios are small regardless of NN (post-hoc cell, one point, weak), and the catching-up length grows with NNzQ(N)z_Q(N). The registered verdict (could not be defined). Implications for four dimensions and SU(N)SU(N)

Sources and reproduction

ItemKindSource / tools
Hierarchical gauge latticesCitedT. Ito, Commun. Math. Phys. 95 (1984) 247 (§1 checked in the text; writes that the extension to higher dimensions is left to the reader). The four-dimensional version is this book's extension
Comparison values for the hierarchical Gaussian free fieldCitedR. Bauerschmidt, D. C. Brydges, G. Slade, arXiv:1907.05474 §4.1.2 (text checked)
Algebra of (3)(4), (5), (7)(9)Proved here, LeanProofs in the text. Bundles LatticeGaugeHierarchicalLink (HierLink.*), LatticeGaugeConservationLength (ConsLen.*), LatticeGaugeConservationExponent (ConsExp.*). The bundles guarantee only arithmetic, algebra and comparisons of intervals
Numbers of Table 1ReadingAll eigenvalues of the dense matrix of the class model (N≤2N\le2), a separate implementation per sector (N≤4N\le4), the gauge quotient on the constraint surface of the four-dimensional original lattice (N≤2N\le2)
Numbers of Table 2ReadingSimulation of link Langevin on the strip (step β dt≈0.04\beta\,dt\approx0.04, strip length ≤256\le256). Positive control: τ\tau of the Gauss dynamics agrees with the closed formula within 1σ1\sigma. Negative control: moving the outer links too gives O(β)O(\beta) independent of MM. Errors by Madras–Sokal, with block averages as a separate route
The formula for ξXY\xi_{XY}, the Kramers rateAssumptionUsed as standard computations (original sources not checked)
The action, topological charge and positive-control values for CPN−1CP^{N-1}CitedM. Hasenbusch, Phys. Rev. D 96 (2017) 054504, arXiv:1706.04443 (equations (1), (10)–(12) and Table II checked in the text). B. Berg, M. Lüscher, Nucl. Phys. B 190 (1981) 412 (the original source for the triangle definition, bibliographic data only)
Dependence on NN of the freezing of the topological modeCitedL. Del Debbio, G. M. Manca, E. Vicari, Phys. Lett. B 594 (2004) 315 (text), M. Berni, C. Bonanno, M. D'Elia, Phys. Rev. D 100 (2019) 114509 (text), M. Campostrini, P. Rossi, E. Vicari, Phys. Rev. D 46 (1992) 2647 (abstract)
Numbers of Tables 3 and 4ReadingSimulation of heat bath (O(3)O(3) heat bath on sites, von Mises heat bath on links) and overrelaxation, autocorrelation with Wolff's window. The free-field reference uses the same algorithm and the same LL
The ultraviolet problem of the lattice topological susceptibility for O(3)O(3)AssumptionFrom memory (original source not obtained)

Terms used in this chapter — the three rulers, the length at which conservation breaks ℓ∗\ell_* (ℓc\ell_c), the Dirichlet form (coordinates), tree coordinates and link coordinates, conservation of flux, wrapping, the topological charge of the strip, the flux dipole, raw rate and static fluctuation, the time for the global quantity to move by one, two-dimensional test beds. Their meanings are in the A-02 Glossary and 11-A1.

Next chapter to read: 11-03 The sector H and the local H — two forms that restate H to avoid the pattern. What is specific to four dimensions (why the two lengths split) is treated in 11-07.

← 11-01 One form that obstructs H — a nearly conserved global quantityContents11-03 Two ways to repair H — the sector H and the local H →

Revised 2026-10-09: first version.