diff b1f9bc4..56bc728 o4_model_mass_bounds_v2.txt
diff --git a/o4_model_mass_bounds_v2.txt b/o4_model_mass_bounds_v2.txt
index 9ab2de3..4de2c60 100644
--- a/o4_model_mass_bounds_v2.txt
+++ b/o4_model_mass_bounds_v2.txt
@@ -43,26 +43,26 @@ newcommand{gap}{m_{lat}}
setcounter{tocdepth}{1}
allowdisplaybreaks[2]
-# Sharp mass bounds for the two-dimensional $O(4)$ model
+# Sharp mass bounds for the two-dimensional O(4) model
maketitle
Abstract:
-We prove sharp mass bounds for the two-dimensional nearest-neighbor $O(4)$
+We prove sharp mass bounds for the two-dimensional nearest-neighbor O(4)
model. For every sufficiently large inverse temperature (inverse coupling)
-$beta$, the full transfer gap -- meaning the exponential decay rate
+beta, the full transfer gap -- meaning the exponential decay rate
of correlations of every local observable, not only the spin -- is
-sandwiched between two positive, $beta$-independent multiples of the
-asymptotically free scale $sqrtbeta,e^{-pibeta}$:
+sandwiched between two positive, beta-independent multiples of the
+asymptotically free scale sqrtbeta,e^{-pibeta}:
[
c,sqrtbeta,e^{-pibeta}
;le; gap(beta) ;le;
C,sqrtbeta,e^{-pibeta},
- qquad betagebeta_0.
+ betagebeta_0.
]
We also prove that the periodic measures have a unique local limit and a
-positive full gap at every finite $beta>0$. Whether the lattice theories
+positive full gap at every finite beta>0. Whether the lattice theories
converge to a continuum quantum field theory, and what its exact spectrum
there looks like, is not decided here; we state separately the
extra convergence hypotheses that would turn our lattice bounds into a
@@ -82,56 +82,56 @@ newpage
## The model in one paragraph
-Place a spin $sigma_xin S^3subsetR^4$ at every site $x$ of a large
-periodic square lattice $Λ=(Z/nZ)times(Z/wZ)$. Give the spins
+Place a spin sigma_xin S^3subsetR^4 at every site x of a large
+periodic square lattice Λ=(Z/nZ)times(Z/wZ). Give the spins
the Boltzmann weight
ddmu_{beta;n,w}(sigma) = Z_beta(n,w)^{-1}
exp!(beta sum_{< xy>}sigma_xcdotsigma_y)
prod_{xinΛ}ddomega(sigma_x),
-where $< xy>$ runs over nearest-neighbor pairs, each counted
-once, $ddomega$ is the normalised uniform measure on $S^3$, and
-$Z_beta(n,w)$ is the normalisation. This is the two-dimensional $O(4)$
+where < xy> runs over nearest-neighbor pairs, each counted
+once, ddomega is the normalised uniform measure on S^3, and
+Z_beta(n,w) is the normalisation. This is the two-dimensional O(4)
nonlinear sigma model on the lattice with nearest-neighbor couplings and
-bare inverse coupling $beta$.
+bare inverse coupling beta.
-The energy $betasigma_xcdotsigma_y = tfracbeta2 - tfracbeta2
-|sigma_x-sigma_y|^2$ favours neighbouring spins pointing in the same
-direction, and $beta$ is how strongly. There is no spin-ordering
+The energy betasigma_xcdotsigma_y = tfracbeta2 - tfracbeta2
+|sigma_x-sigma_y|^2 favours neighbouring spins pointing in the same
+direction, and beta is how strongly. There is no spin-ordering
transition in two dimensions (Mermin--Wagner), but there is a length
scale beyond which correlations die exponentially. That length scale is
what we want to bracket.
## What ``mass gap'' means
-On a cylinder of spatial width $w$ the transfer operator $K_w$ advances
+On a cylinder of spatial width w the transfer operator K_w advances
the configuration one time step. Its kernel is
K_w(s,s')=exp!( (beta)/(2)V_w(s)
+betasum_{i=1}^w s_icdot s'_i
+(beta)/(2)V_w(s') ),
- qquad V_w(s)=sum_i s_icdot s_{i+1}.
+ V_w(s)=sum_i s_icdot s_{i+1}.
-$K_w$ is a positive, self-adjoint, compact integral operator. Its largest
-eigenvalue $lambda_1$ is simple and its positive ground eigenvector is
-the vacuum. Closing time periodically gives $Z_beta(n,w)=Tr K_w^n$.
+K_w is a positive, self-adjoint, compact integral operator. Its largest
+eigenvalue lambda_1 is simple and its positive ground eigenvector is
+the vacuum. Closing time periodically gives Z_beta(n,w)=Tr K_w^n.
-Fix $beta$, take $wto∞$, and take a subsequential local limit. On
+Fix beta, take wto∞, and take a subsequential local limit. On
the resulting reflected Hilbert space there is a positive one-step
-translation contraction $T_beta$ with $T_betaΩ = Ω$. The
+translation contraction T_beta with T_betaΩ = Ω. The
full transfer gap is
gap(beta) = -lognorm{T_betabig|_{Ω^perp}}.
-Equivalently, autocorrelations $< F, T_beta^d F>$ of every
-centred local observable $F$ decay at least as fast as $e^{-gap d}$,
-where $d$ is the time separation. The gap is the mass; $gap=0$ means
+Equivalently, autocorrelations < F, T_beta^d F> of every
+centred local observable F decay at least as fast as e^{-gap d},
+where d is the time separation. The gap is the mass; gap=0 means
critical or ordered.
The subtlety: controlling the spin two-point function is not enough.
There is a second, rotation-invariant sector of observables -- the
-bond energies $sigma_xcdotsigma_y$ -- and $T_beta$ has spectral
+bond energies sigma_xcdotsigma_y -- and T_beta has spectral
support in that sector as well. The full gap includes both. Our
theorem bounds this bigger, smaller gap. That is why the analysis has to
carry along every local observable, not just the spin.
@@ -139,8 +139,8 @@ carry along every local observable, not just the spin.
## The answer
[Main result]
-There are finite constants $beta_0,c,C$, with $0<c<C$, such that for
-every $betagebeta_0$ the periodic measures (eq:model-ov) have a
+There are finite constants beta_0,c,C, with 0<c<C, such that for
+every betagebeta_0 the periodic measures (eq:model-ov) have a
unique local limit, and its full transfer gap obeys
c,sqrtbeta,e^{-pibeta};le;gap(beta);le;
@@ -148,56 +148,55 @@ unique local limit, and its full transfer gap obeys
-Three comments on the formula $sqrtbeta,e^{-pibeta}$.
-
-item The exponent is not a guess; it is forced by the renormalisation
+Three comments on the formula sqrtbeta,e^{-pibeta}.
+ The exponent is not a guess; it is forced by the renormalisation
group (RG). Asymptotic freedom says the running coupling satisfies
- $betamapstobeta - gamma - gamma/(2pibeta)+...$ with
- $gamma=(log L)/pi$, and integrating over many scale factors gives
- $log(inverse length)simpibeta$. So it is a two-loop
- calculation that produces $e^{-pibeta}$ with a $sqrtbeta$
+ betamapstobeta - gamma - gamma/(2pibeta)+... with
+ gamma=(log L)/pi, and integrating over many scale factors gives
+ log(inverse length)simpibeta. So it is a two-loop
+ calculation that produces e^{-pibeta} with a sqrtbeta
correction.
-item Up to bounded positive factors, $sqrtbeta,e^{-pibeta}$ is
+ Up to bounded positive factors, sqrtbeta,e^{-pibeta} is
the natural asymptotically free scale for this model in these
units. That is why the bounds are sharp.
-item The bounds are not asymptotic (the ratio need not tend to $1$, and
- in fact the exact $O(4)$ prefactor is a known transcendent constant
+ The bounds are not asymptotic (the ratio need not tend to 1, and
+ in fact the exact O(4) prefactor is a known transcendent constant
we do not touch here). They give the correct order of magnitude.
-[Full gap at every finite $beta$]
-For every finite $beta>0$, the periodic measures on squares have a unique
+[Full gap at every finite beta]
+For every finite beta>0, the periodic measures on squares have a unique
local limit. Its full transfer gap is positive, with
[
gap(beta);ge; c,exp!(-pibeta
- C√((1+beta)log(2+beta)))>0
]
-for some absolute constants $c,C>0$.
+for some absolute constants c,C>0.
Corollary~cor:allT-ov does not use the sharp two-loop
comparison; it follows from a preliminary mixing estimate alone. The
gap is therefore positive at every temperature, but the formula is too
-weak in the exponent to distinguish $sqrtbeta,e^{-pibeta}$ from
-$beta^{-1/2}e^{pibeta+sqrtbeta}$. Theorem~thm:main-ov repairs
-that at large $beta$.
+weak in the exponent to distinguish sqrtbeta,e^{-pibeta} from
+beta^{-1/2}e^{pibeta+sqrtbeta}. Theorem~thm:main-ov repairs
+that at large beta.
-## From $gap$ to a physical mass
+## From gap to a physical mass
-Choose a physical reference length $ell_{ref}>0$ and declare
-that $L^N$ lattice spacing (where $L$ is a fixed block scale we will
-choose later) has physical length $ell_{ref}$. The lattice
-spacing is $a_N=ell_{ref}L^{-N}$. Then
-$gap(beta_N)/a_N$ is the mass of the cutoff theory in physical units.
+Choose a physical reference length ell_{ref}>0 and declare
+that L^N lattice spacing (where L is a fixed block scale we will
+choose later) has physical length ell_{ref}. The lattice
+spacing is a_N=ell_{ref}L^{-N}. Then
+gap(beta_N)/a_N is the mass of the cutoff theory in physical units.
Theorem~thm:main-ov upgrades to
frac{log 2}{M_0ell_{ref}}
;le; (gap(beta_N))/(a_N)
;le; frac{log 8}{ell_{ref}}
- qquad(N large),
+ (N large),
-with $M_0$ a fixed positive integer and $beta_N$ chosen so that exactly
-$N$ RG steps have been done at terminal coupling $H$. Bounded positive
+with M_0 a fixed positive integer and beta_N chosen so that exactly
+N RG steps have been done at terminal coupling H. Bounded positive
multiplicative factors sandwich the continuum-normalised mass. Whether
the cutoff measures converge to an actual continuum theory is a separate
question; we spell out in Section~sec:physical the additional
@@ -208,15 +207,15 @@ statement about an existing continuum spectrum.
Absence of ordinary ordering is classical (Mermin--Wagner, and the
McBryan--Spencer quantitative bounds). Exponential decay of the spin
-correlation was proved separately for every finite $beta$ and every
-internal dimension $ge 3$. Perturbative renormalisation predicts the
-scale $sqrtbeta,e^{-pibeta}$; Bethe ansatz $S$-matrix theory fixes
-the exact continuum mass for $O(3)$ and $O(4)$, and for $O(n)$ in
+correlation was proved separately for every finite beta and every
+internal dimension ge 3. Perturbative renormalisation predicts the
+scale sqrtbeta,e^{-pibeta}; Bethe ansatz S-matrix theory fixes
+the exact continuum mass for O(3) and O(4), and for O(n) in
general. None of that gave a full transfer gap (including invariant
observables) for the specified lattice model. Our theorem closes that
-gap for $O(4)$ at large $beta$, and yields positivity at every finite
-$beta$ as a corollary. The proof of the exact continuum mass
-(including the ``$Gamma$-function'' prefactor) is a separate problem,
+gap for O(4) at large beta, and yields positivity at every finite
+beta as a corollary. The proof of the exact continuum mass
+(including the ``Gamma-function'' prefactor) is a separate problem,
and the proof of lattice-to-continuum convergence is also separate.
## The proof in six steps
@@ -224,28 +223,28 @@ and the proof of lattice-to-continuum convergence is also separate.
The proof has a clear overall shape.
{bf Step 1} (Section~sec:finite-volume). Reduce the gap to a
finite-volume partition-function defect
-$Delta_beta(n)=4log Z_beta(n,n)-log Z_beta(2n,2n)$.
-{bf Step 2} (Section~sec:preliminary). Show $Delta_beta(n)$ is
-already small at some scale $nsim e^{pibeta+C√(betalogbeta)}$,
+Delta_beta(n)=4log Z_beta(n,n)-log Z_beta(2n,2n).
+{bf Step 2} (Section~sec:preliminary). Show Delta_beta(n) is
+already small at some scale nsim e^{pibeta+C√(betalogbeta)},
using Ward identities and a stiffness iteration.
{bf Step 3} (Sections~sec:free-gap-ov-run--sec:calibration).
Construct an exact real-space renormalisation group (RG) step that
-integrates out $Ltimes L$ blocks and rewrites the coarse lattice in the
-same controlled form. Show it contracts past a coupling of size $Hsim$
+integrates out Ltimes L blocks and rewrites the coarse lattice in the
+same controlled form. Show it contracts past a coupling of size Hsim
constant.
{bf Step 4} (Section~sec:trajectories). Show that, at every
-depth $N$, one can pick a bare coupling $beta_N$ so the $N$-step
-trajectory ends exactly at $H$/c. Matching at the coarse end absorbs
+depth N, one can pick a bare coupling beta_N so the N-step
+trajectory ends exactly at H/c. Matching at the coarse end absorbs
the difference between two different depths.
{bf Step 5} (Section~sec:comparison). Compare the positive
coarse densities (partition functions/correlation functions) rather
than the signed expansions; positivity of the original model does the
hard work.
{bf Step 6} (Sections~sec:assembly--sec:upper). Choose one
-reference box once at a coarse depth $K$, show Step 2 already makes
-$Delta$ small there, and Step 4 keeps it small at every finer cutoff.
+reference box once at a coarse depth K, show Step 2 already makes
+Delta small there, and Step 4 keeps it small at every finer cutoff.
Then apply Step 1 to get the lower mass bound. Prove the upper bound by
-noticing that the RG produces a definite correlation at distance $L^N$.
+noticing that the RG produces a definite correlation at distance L^N.
The rest of the paper expands each step.
@@ -254,52 +253,51 @@ The rest of the paper expands each step.
%---------------------------------------------------------------------
-Let $Λ_{n,w}=(Z/nZ)times(Z/wZ)$ be a periodic rectangle and
-$sigma_xin S^3subsetR^4$ with normalised uniform measure $ddomega$
-($int_{S^3}ddomega=1$). Periods are even and at least $4$. The
+Let Λ_{n,w}=(Z/nZ)times(Z/wZ) be a periodic rectangle and
+sigma_xin S^3subsetR^4 with normalised uniform measure ddomega
+(int_{S^3}ddomega=1). Periods are even and at least 4. The
probability measure is (eq:model-ov). Every unoriented
nearest-neighbour bond is counted once.
## Transfer operator and reflection positivity
-Fix a spatial circumference $w$. A time-slice configuration is
-$s=(s_1,...,s_w)in(S^3)^w$; the one-step interaction is
-$V_w(s)=sum_{i=1}^w s_icdot s_{i+1}$ with periodic indices. The transfer
-operator $K_w$ has kernel (eq:transfer-kernel-ov). It is:
-
-item compact and self-adjoint: $beta$ finite, $S^3$ compact,
+Fix a spatial circumference w. A time-slice configuration is
+s=(s_1,...,s_w)in(S^3)^w; the one-step interaction is
+V_w(s)=sum_{i=1}^w s_icdot s_{i+1} with periodic indices. The transfer
+operator K_w has kernel (eq:transfer-kernel-ov). It is:
+ compact and self-adjoint: beta finite, S^3 compact,
kernel continuous;
-item positive semidefinite: expand
- $e^{beta scdot s'}=sum_m beta^m (scdot s')^m/m!$ and each
- $(scdot s')^m$ is a positive tensor-product kernel, then sandwich
- with $e^{beta V_w/2}$;
-item trace class: the diagonal integral of the positive kernel
+ positive semidefinite: expand
+ e^{beta scdot s'}=sum_m beta^m (scdot s')^m/m! and each
+ (scdot s')^m is a positive tensor-product kernel, then sandwich
+ with e^{beta V_w/2};
+ trace class: the diagonal integral of the positive kernel
series converges;
-item strictly positive: strictly positive kernel; hence by
- Perron--Frobenius, $lambda_1(w)$ is simple and the top
+ strictly positive: strictly positive kernel; hence by
+ Perron--Frobenius, lambda_1(w) is simple and the top
eigenvector (``vacuum'') is positive.
-Closing time periodically gives $Z_beta(n,w)=Tr K_w^n
-=sum_i lambda_i(w)^n$ with
-$lambda_1(w)>lambda_2(w)gelambda_3(w)ge...ge0$.
+Closing time periodically gives Z_beta(n,w)=Tr K_w^n
+=sum_i lambda_i(w)^n with
+lambda_1(w)>lambda_2(w)gelambda_3(w)ge...ge0.
The nice structural property is reflection positivity, following
from the site-reflection form
-$< F,G>_OS=E[(Θ barbar F)G]$F)G]$ where $Θ$
+< F,G>_OS=E[(Θ barbar F)G]F)G] where Θ
reflects in a time plane. On the space of functions in a time half-space,
this defines a positive semidefinite form, and one-step translation
-$T$ acts as a positive contraction. We record the facts we use, then
+T acts as a positive contraction. We record the facts we use, then
prove them in Section~sec:finite-volume.
[Uniqueness of periodic local limit]
The periodic measures (eq:model-ov) on squares have a unique local
-limit as the side tends to infinity, at each fixed $beta$.
+limit as the side tends to infinity, at each fixed beta.
[Reflection form]
-The vacuum $Ω$, positive contraction $T_beta$, and full gap
-$gap(beta)= -lognorm{T_beta|_{Ω^perp}}$ exist as in
-(eq:gap-def-ov), and $gap(beta)ge0$.
+The vacuum Ω, positive contraction T_beta, and full gap
+gap(beta)= -lognorm{T_beta|_{Ω^perp}} exist as in
+(eq:gap-def-ov), and gap(beta)ge0.
The proofs are in Section~sec:finite-volume (§sec:pre-part-12-2)
@@ -307,35 +305,35 @@ and use the preliminary mixing estimate of Section~sec:preliminary.
## Two scales: block depth and bare coupling
-Fix once and for all a small dyadic block size $L$ (later chosen large)
-and a large terminal inverse coupling $H$ (later chosen large). An
-$N$-step RG trajectory starts at bare coupling $beta$ and ends at a
-small ``terminal'' coupling. Call $N$ the depth. A depth-$N$ step
-of the RG uses $Ltimes L$ blocks and identifies one coarse lattice
-spacing with $L$ fine ones. After $N$ steps one coarse cell is $L^N$
+Fix once and for all a small dyadic block size L (later chosen large)
+and a large terminal inverse coupling H (later chosen large). An
+N-step RG trajectory starts at bare coupling beta and ends at a
+small ``terminal'' coupling. Call N the depth. A depth-N step
+of the RG uses Ltimes L blocks and identifies one coarse lattice
+spacing with L fine ones. After N steps one coarse cell is L^N
microscopic cells.
-Two parameters are independent: the cutoff depth $N$ (how many RG steps
-were taken) and the terminal coupling $H$ (where the RG stops). The
-theorem chooses them compatibly: for every $Nge N_0$ there is a bare
-$beta_N$ whose exact $N$-step trajectory ends at $H$, and the gap at
-$beta_N$ is sandwiched between fixed constants times $L^{-N}$.
+Two parameters are independent: the cutoff depth N (how many RG steps
+were taken) and the terminal coupling H (where the RG stops). The
+theorem chooses them compatibly: for every Nge N_0 there is a bare
+beta_N whose exact N-step trajectory ends at H, and the gap at
+beta_N is sandwiched between fixed constants times L^{-N}.
-Now note: $L^{-N}$ and $sqrtbeta,e^{-pibeta}$ are the same
+Now note: L^{-N} and sqrtbeta,e^{-pibeta} are the same
scale up to constants, because the RG integrates the running coupling
-from $beta$ down to $H$ in $N$ steps of size $gamma=(log L)/pi$.
+from beta down to H in N steps of size gamma=(log L)/pi.
Section~sec:calibration makes that identification precise.
The rest of the paper is the proof.
## Six-step map, with cross-references
-noindent
-Step 1: $Delta_beta$ criterion dotfill Ssec:finite-volume
+
+Step 1: Delta_beta criterion dotfill Ssec:finite-volume
Step 2: preliminary mixing dotfill Ssec:preliminary
Step 3: linear block dotfill Ssec:linear
-Step 3$'$: exact nonlinear RG dotfill Ssec:rg
-Step 3$''$: calibration dotfill Ssec:calibration
+Step 3': exact nonlinear RG dotfill Ssec:rg
+Step 3'': calibration dotfill Ssec:calibration
Step 4: admission and matching dotfill Ssec:trajectories
Step 5: positive-density comparison dotfill Ssec:comparison
Step 6: assembly and upper bound dotfill Ssec:assembly, Ssec:upper
@@ -343,13 +341,13 @@ Physical units dotfill Ssec:physical
## Notation
-$gap(beta)$ = full transfer gap. $rho(beta)=sqrtbeta e^{-pibeta}$
-= asymptotically free inverse length. $L$ = blocking factor (dyadic,
-large). $gamma=(log L)/pi$. $Λ_{n,w}=(Z/nZ)times(Z/wZ)$;
-$v=|Λ_{n,w}|$. $Z_beta(n,w)$ = partition function. $E$ = expectation.
-$norm{cdot}$ = operator or supremum norm, context clear. $C,c>0$ =
+gap(beta) = full transfer gap. rho(beta)=sqrtbeta e^{-pibeta}
+= asymptotically free inverse length. L = blocking factor (dyadic,
+large). gamma=(log L)/pi. Λ_{n,w}=(Z/nZ)times(Z/wZ);
+v=|Λ_{n,w}|. Z_beta(n,w) = partition function. E = expectation.
+norm{cdot} = operator or supremum norm, context clear. C,c>0 =
generic constants that may change between formulas but are
-$beta$-independent unless stated.
+beta-independent unless stated.
%---------------------------------------------------------------------
# A finite-volume criterion for the full gap
@@ -359,7 +357,7 @@ $beta$-independent unless stated.
The first step is a purely algebraic criterion. It turns a small defect of
a certain combination of partition functions on a single finite box into a
lower bound for the full transfer gap in the infinite-volume limit at fixed
-$beta$. The criterion sees every spectral sector, because it uses
+beta. The criterion sees every spectral sector, because it uses
the complete trace of the transfer operator; and the specific combination
of partition functions we use is invariant under the addition of any free
energy proportional to area, which is what will make it a perfect ``test''
@@ -367,31 +365,31 @@ for the RG comparison later. No renormalisation is used in this section.
## Eigenvalues, purity, and a single number per side
-Fix a spatial circumference $w$ and let $K_w$ act on $(S^3)^w$
+Fix a spatial circumference w and let K_w act on (S^3)^w
as in (eq:transfer-kernel-ov). Its eigenvalues satisfy
-$lambda_1(w)>lambda_2(w)gelambda_3(w)ge...ge0$, and
-$Z_beta(n,w)=Tr K_w^n=sum_ilambda_i(w)^n$.
+lambda_1(w)>lambda_2(w)gelambda_3(w)ge...ge0, and
+Z_beta(n,w)=Tr K_w^n=sum_ilambda_i(w)^n.
-Suppress $beta$ and define
+Suppress beta and define
p(n,w) &= 2log Z(n,w) - log Z(2n,w),
q(n,w) &= 2log Z(n,w) - log Z(n,2w).
The same numbers, one with time doubled and one with space doubled.
-[Both $p$ and $q$ are non-negative]
-For every $n,w$,
+[Both p and q are non-negative]
+For every n,w,
[
- p(n,w)ge0,qquad q(n,w)ge 0.
+ p(n,w)ge0, q(n,w)ge 0.
]
Proof.
-Consider non-negative eigenvalues $mu_i$ of a positive self-adjoint
-operator and their $n$-th powers. The Cauchy-Schwarz bound
-$sum_imu_i^{2n}le(sum_imu_i^n)^2$ gives
-$p(n,w)ge0$. Applying the same argument to the space direction gives
-$q(n,w)ge0$.
+Consider non-negative eigenvalues mu_i of a positive self-adjoint
+operator and their n-th powers. The Cauchy-Schwarz bound
+sum_imu_i^{2n}le(sum_imu_i^n)^2 gives
+p(n,w)ge0. Applying the same argument to the space direction gives
+q(n,w)ge0.
Introduce the normalised trace weights
@@ -401,7 +399,7 @@ Introduce the normalised trace weights
e^{-p(n,w)} = sum_i t_i^2.
The right side is the probability that two independent copies of the
-chain, drawn with weights $t_i$, meet in the same eigenstate. It is $le1$
+chain, drawn with weights t_i, meet in the same eigenstate. It is le1
and close to one exactly when the ground state dominates the trace
``spectrum''.
@@ -410,20 +408,20 @@ Set
[
s_n(w)=sum_{ige2}(lambda_i(w)/lambda_1(w))^n.
]
-Then $s_n(w)le e^{p(n,w)}-1$. If $p(n,w)le1/4$, then
+Then s_n(w)le e^{p(n,w)}-1. If p(n,w)le1/4, then
((lambda_2(w))/(lambda_1(w)))^nle 2p(n,w),
- qquad p(2n,w)le 4p(n,w)^2.
+ p(2n,w)le 4p(n,w)^2.
-The analogous statements hold for $q$ when the other period is doubled.
+The analogous statements hold for q when the other period is doubled.
Proof.
-$t_1=(1+s_n)^{-1}$ and $sum_i t_i^2le t_1$ give $s_nle e^{p}-1$.
-For $ple1/4$, $e^p-1le pe^ple2p$, which bounds the first excited term.
-The largest eigenvalue cancels from the definition of $p$, so
-$p(2n,w)=2log(1+s_{2n})-log(1+s_{4n})le2s_{2n}le2s_n^2
-le2p^2e^{2p}le4p^2$.
+t_1=(1+s_n)^{-1} and sum_i t_i^2le t_1 give s_nle e^{p}-1.
+For ple1/4, e^p-1le pe^ple2p, which bounds the first excited term.
+The largest eigenvalue cancels from the definition of p, so
+p(2n,w)=2log(1+s_{2n})-log(1+s_{4n})le2s_{2n}le2s_n^2
+le2p^2e^{2p}le4p^2.
## Why space directions need a separate argument
@@ -433,7 +431,7 @@ low-energy states from appearing as the circumference grows. The
following exact identity lets us control that limit.
[Rectangle identity]
-For every $n,w$,
+For every n,w,
p(n,2w) = 2p(n,w) - 2q(n,w) + q(2n,w).
@@ -441,7 +439,7 @@ For every $n,w$,
Proof.
Expand the definitions (eq:p-def)--(eq:q-def). Both sides
-equal $4log Z(n,w) - 2log Z(2n,w) - 2log Z(n,2w) + log Z(2n,2w)$.
+equal 4log Z(n,w) - 2log Z(2n,w) - 2log Z(n,2w) + log Z(2n,2w).
There is a companion identity obtained by exchanging the two directions.
@@ -452,38 +450,38 @@ Now define the square-box test
Delta_beta(n) = 4log Z_beta(n,n) - log Z_beta(2n,2n).
The combination has the important property that any contribution of the
-form $f_beta nw$ to $log Z_beta(n,w)$ cancels from $Delta_beta(n)$.
+form f_beta nw to log Z_beta(n,w) cancels from Delta_beta(n).
Bulk free energy generated by the RG therefore does not enter the test.
-[Two ways to write $Delta_beta$]
-For every even $nge4$,
+[Two ways to write Delta_beta]
+For every even nge4,
[
Delta_beta(n) = 2p(n,n) + q(2n,n) ge 0.
]
-Consequently $Delta_beta(n)$ is a sum of two non-negative purities
+Consequently Delta_beta(n) is a sum of two non-negative purities
combined, and vanishes at infinite inverse size in the trivial regime.
Proof.
-By direct expansion, $4log Z(n,n) - log Z(2n,2n)$ equals
-$2(2log Z(n,n)-log Z(2n,n)) + (2log Z(2n,n)-log Z(2n,2n)$, which is
-$2p(n,n)+q(2n,n)$. Non-negativity is the previous lemma.
+By direct expansion, 4log Z(n,n) - log Z(2n,2n) equals
+2(2log Z(n,n)-log Z(2n,n)) + (2log Z(2n,n)-log Z(2n,2n), which is
+2p(n,n)+q(2n,n). Non-negativity is the previous lemma.
## The criterion
[Finite-volume criterion]
-Let $nge4$ be even. If
+Let nge4 be even. If
Delta_beta(n)leeta_*:=2^{-10},
-then the transfer operators on circumferences $2^kn$, $kge0$, satisfy
+then the transfer operators on circumferences 2^kn, kge0, satisfy
(lambda_2(2^kn))/(lambda_1(2^kn))
le e^{-(log 2)/n}.
-Every local limit of the square tori with sides $2^kn$ has full gap at
-least $(log2)/n$. In particular, if the periodic infinite-plane state
+Every local limit of the square tori with sides 2^kn has full gap at
+least (log2)/n. In particular, if the periodic infinite-plane state
is unique, then
gap(beta)ge(log2)/n.
@@ -491,51 +489,51 @@ is unique, then
Proof.
-Write $Delta_r=Delta_beta(r)$. Besides $q(2r,r)leDelta_r$, the
-definitions give $p(r,2r)=Delta_r-2q(r,r)leDelta_r$.
-Lemma~lem:squaring therefore gives $p(2r,2r)le4Delta_r^2$. The
+Write Delta_r=Delta_beta(r). Besides q(2r,r)leDelta_r, the
+definitions give p(r,2r)=Delta_r-2q(r,r)leDelta_r.
+Lemma~lem:squaring therefore gives p(2r,2r)le4Delta_r^2. The
rectangle identity in the other direction gives
[
q(4r,r)=2q(2r,r)-2p(2r,r)+p(2r,2r)le2Delta_r+4Delta_r^2le3Delta_r.
]
-Squaring in that direction gives $q(4r,2r)le36Delta_r^2$. Consequently
+Squaring in that direction gives q(4r,2r)le36Delta_r^2. Consequently
[
Delta_beta(2r)=2p(2r,2r)+q(4r,2r)le44Delta_r^2le64Delta_r^2.
]
-Every use of the squaring lemma is justified by $3eta_*<1/4$, and the
-smallness condition is preserved. Set $n_k=2^kn$. By induction,
+Every use of the squaring lemma is justified by 3eta_*<1/4, and the
+smallness condition is preserved. Set n_k=2^kn. By induction,
[
Delta_beta(n_k)lefrac1{64}(64eta_*)^{2^k}lefrac1{64}2^{-2^k}.
]
-Since $( lambda_2(n_k)/lambda_1(n_k))^{n_k} le 2p(n_k,n_k)leDelta_beta(n_k)$
-by the squaring lemma, taking $n_k$-th roots proves
+Since ( lambda_2(n_k)/lambda_1(n_k))^{n_k} le 2p(n_k,n_k)leDelta_beta(n_k)
+by the squaring lemma, taking n_k-th roots proves
(eq:width-gap).
-Now pass to the limit. Let $Ω_w$ be the normalised positive ground
-eigenvector of $K_w$. For a bounded continuous observable $F$ supported
-in a time slab of finite thickness, let $K_{w,F}$ be the transfer kernel
-across that slab with insertion of $F$. Pointwise
-$|K_{w,F}|lenorm F_∞ K_w^r$ where $r$ is the number of transfer
-steps across the slab; set $A_F=K_{w,F}/lambda_1(w)^r$ for $rge0$
-(and $r=0$ gives multiplication by $F$). Its operator norm is at most
-$norm F_∞$. The infinite-time cylinder expectation is
-$omega_w(F)=inner{Ω_w}{A_FΩ_w}$. Applying (eq:width-gap)
+Now pass to the limit. Let Ω_w be the normalised positive ground
+eigenvector of K_w. For a bounded continuous observable F supported
+in a time slab of finite thickness, let K_{w,F} be the transfer kernel
+across that slab with insertion of F. Pointwise
+|K_{w,F}|lenorm F_∞ K_w^r where r is the number of transfer
+steps across the slab; set A_F=K_{w,F}/lambda_1(w)^r for rge0
+(and r=0 gives multiplication by F). Its operator norm is at most
+norm F_∞. The infinite-time cylinder expectation is
+omega_w(F)=inner{Ω_w}{A_FΩ_w}. Applying (eq:width-gap)
between two such insertions yields
|Cov_{omega_w}(F,G)|
le norm F_∞norm G_∞ e^{-(log2)d/n},
-where $d$ is the number of transfer edges separating their supports. This
+where d is the number of transfer edges separating their supports. This
bound survives any subsequential local limit of the cylinder states as
-$w=2^knto∞$.
+w=2^knto∞.
Now identify cylinder limits with square-torus limits. Let
-$widehatwidehat K_w=K_w/K_w=K_w/lambda_1(w)$. For $w=n_k$ and $rle w/2$, the
+widehatwidehat K_w=K_w/K_w=K_w/lambda_1(w). For w=n_k and rle w/2, the
square-torus expectation equals
-$Tr(A_Fwidehatwidehat K_w^K_w^{w-r})/Trwidehatwidehat K_w^w$.K_w^w$. Splitting off the ground
+Tr(A_Fwidehatwidehat K_w^K_w^{w-r})/Trwidehatwidehat K_w^w.K_w^w. Splitting off the ground
projection in each trace shows the difference from the cylinder
-expectation is at most $2norm F_∞ s_{w/2}(w)$. But
-$p(w/2,w)leDelta_beta(w/2)$ and $s_{w/2}(w)le e^{Delta_beta(w/2)}-1to0$.
+expectation is at most 2norm F_∞ s_{w/2}(w). But
+p(w/2,w)leDelta_beta(w/2) and s_{w/2}(w)le e^{Delta_beta(w/2)}-1to0.
Thus the square and cylinder expectations have the same subsequential
limits, even for products of fixed local observables, and
(eq:slab-bound) holds in every such square-torus limit.
@@ -545,7 +543,7 @@ by a local observable in the reflection-positive Hilbert space, the
one-step transfer autocorrelation has a positive spectral measure
(reflection positivity). Any such autocorrelation is a correlation of
two reflected slab observables, so (eq:slab-bound) applies and
-excludes spectral support above $e^{-(log2)/n}$. Such vectors are dense
+excludes spectral support above e^{-(log2)/n}. Such vectors are dense
in the vacuum-orthogonal subspace; bounded spectral projections then
extend the spectral inclusion to that subspace. This proves
(eq:criterion-gap).
@@ -557,8 +555,8 @@ on the size of a slab is needed. That is useful when the upper mass
bound comes up in Section~sec:upper.
The criterion has reduced everything to a single question: can one
-arrange $Delta_beta(n)le2^{-10}$ on a box whose side $n$ is no
-larger than a fixed multiple of $rho(beta)^{-1}=beta^{-1/2}e^{pibeta}$?
+arrange Delta_beta(n)le2^{-10} on a box whose side n is no
+larger than a fixed multiple of rho(beta)^{-1}=beta^{-1/2}e^{pibeta}?
The preliminary estimate of Section~sec:preliminary will show this
is true with a slightly too-large constant in the exponent; that is not
enough. The point of the RG that follows is to remove the excess.
@@ -571,43 +569,43 @@ enough. The point of the RG that follows is to remove the excess.
%---------------------------------------------------------------------
In this section we prove the first source of exponential decay, with a
-length scale $Xi(beta)$ that is bigger than the asymptotically free
+length scale Xi(beta) that is bigger than the asymptotically free
scale by a subexponential factor. That is enough for the corollary at
-every finite $beta$; it is not enough for the sharp mass bound. But the
-$Delta_beta$ bound from this section will be used at one finite
+every finite beta; it is not enough for the sharp mass bound. But the
+Delta_beta bound from this section will be used at one finite
reference box, where a subexponential excess is harmless. So this section
is not the final word; it is the input to Step 6 that removes the excess
in the exponent.
-Throughout, spins take values in the unit sphere $S^3subsetR^4$,
+Throughout, spins take values in the unit sphere S^3subsetR^4,
identified with the unit quaternions, with normalised Haar measure.
Distances between supports use maximum norm (or its periodic version).
## What we will prove
[Preliminary mixing]
-There are absolute constants $C,c,p>0$ such that, with
+There are absolute constants C,c,p>0 such that, with
Xi(beta)=exp!Bigl[pibeta+C√((1+beta)log(2+beta))Bigr],
-the following hold for every $betage0$.
+the following hold for every betage0.
For every unpinned rectangular torus with each side at least
-$CXi(beta)$, every choice $0lebeta_elebeta$ of nearest-neighbour
-couplings, and any local observables $F,G$ supported on finite sets
-$A,A'$ at distance $d$,
+CXi(beta), every choice 0lebeta_elebeta of nearest-neighbour
+couplings, and any local observables F,G supported on finite sets
+A,A' at distance d,
|Cov(F,G)|le C(1+beta+|A|+|A'|+d)^{20}L_F L_G
e^{-cd/Xi(beta)}.
-Here $L_F=norm F_∞+max_{xin A}Lip_x(F)$, the sum of
+Here L_F=norm F_∞+max_{xin A}Lip_x(F), the sum of
sup and individual Lipschitz constants. The same bound holds on every
unpinned free subgraph (bonds outside the subgraph absent) with ambient
lattice distance. In the homogeneous model the local limit of periodic
squares is unique.
-Finally, if $Z_beta(n,m)$ is the homogeneous torus partition function,
-then for every even $nge CXi(beta)$,
+Finally, if Z_beta(n,m) is the homogeneous torus partition function,
+then for every even nge CXi(beta),
0leDelta_beta(n):=4log Z_beta(n,n)-log Z_beta(2n,2n)
le C(1+beta)^pn^p e^{-cn/Xi(beta)}.
@@ -616,40 +614,40 @@ then for every even $nge CXi(beta)$,
The bound (eq:pre:theorem-delta) together with Step 1
(Proposition~prop:criterion) already gives
-$gap(beta)ge(log2)/n$ for $napproxXi(beta)$, and hence
+gap(beta)ge(log2)/n for napproxXi(beta), and hence
Corollary~cor:allT-ov. But this is far from sharp: the true
-inverse length is $rho(beta)=beta^{-1/2}e^{pibeta}$, and
-$Xi(beta)/rho(beta)^{-1}to∞$ as $betato∞$. The RG
+inverse length is rho(beta)=beta^{-1/2}e^{pibeta}, and
+Xi(beta)/rho(beta)^{-1}to∞ as betato∞. The RG
comparison below will remove this excess.
-subsection*{Structure of the proof of Theorem~thm:preliminary}
+Structure of the proof of Theorem~thm:preliminary
The proof has four macroscopic stages that we will now run in order.
-noindent
+
Stage A (Ward identities + stiffness iteration).
-Show that, at the cost of a multiplicative increase by a factor $L>1$ in
+Show that, at the cost of a multiplicative increase by a factor L>1 in
side length, the ``stiffness'' of the model -- a measure of how much the
free energy gets pushed back when you twist the boundary spins -- is
-reduced by $log L/pi$ minus an error. Iterating decreases stiffness at
-least like $log(side length)/pi$. Start with stiffness $beta$ at
-side length one. After reaching a side of length $exp(pibeta+C√(betalogbeta))$
+reduced by log L/pi minus an error. Iterating decreases stiffness at
+least like log(side length)/pi. Start with stiffness beta at
+side length one. After reaching a side of length exp(pibeta+C√(betalogbeta))
the stiffness is tiny.
-noindent
+
Stage B (small stiffness ⟹ cheap rotations).
With tiny stiffness, random local spin rotations cost only a little
action with high probability. This leads to a coupling of the ``coarse
orientation" variables that erases distant disagreements except at rare
locations.
-noindent
+
Stage C (disagreement percolation).
Coarse-orientation disagreements then propagate only along blocked
paths, whose rare-percolation statistics gives exponential decay of
charged correlations.
-noindent
+
Stage D (all observables, not just charged ones).
Use quaternion algebra: writing a unit quaternion as a pair of circle
valued fields and then decomposing further into Ising variables lets us
@@ -658,36 +656,36 @@ observables in ferromagnetic subsystems. Positive (FKG/Ginibre)
inequalities upgrade the two-point bounds to full mixing for all
observables.
-Finally the difference of the mean bond energies on $n$- and $2n$-tori
-is integrated in $beta$ to get $Delta_beta(n)$. This is
+Finally the difference of the mean bond energies on n- and 2n-tori
+is integrated in beta to get Delta_beta(n). This is
(eq:pre:theorem-delta).
The rest of the section proves each stage.
## Stage A, part 1: the local class
-Work on a pinned square $Q(S)={0,...,S}^2$. For each positively
-oriented edge $e=xto y$, write
+Work on a pinned square Q(S)={0,...,S}^2. For each positively
+oriented edge e=xto y, write
- t_e + s_e = q_y,U_e,q_x^{-1},qquad t_einR,quad s_einImmathbb H,bbH,
+ t_e + s_e = q_y,U_e,q_x^{-1}, t_einR, s_einImmathbb H,bbH,
-where $U_e$ is a fixed unit quaternion encoding a possibly twisted edge.
-The scalar $t_e=q_ycdot q_x$ is the bond energy up to twist, and the
-imaginary quaternion $s_e$ is the ``current'' along the edge: its three
+where U_e is a fixed unit quaternion encoding a possibly twisted edge.
+The scalar t_e=q_ycdot q_x is the bond energy up to twist, and the
+imaginary quaternion s_e is the ``current'' along the edge: its three
components measure the three perpendicular rotation axes. The pinning is
encoded by taking a fixed boundary spin.
-Give each edge a positive coefficient $b_elebeta$ (interior edges), with
-$b_e/2$ for edges tangent to the boundary (the half allocation makes
+Give each edge a positive coefficient b_elebeta (interior edges), with
+b_e/2 for edges tangent to the boundary (the half allocation makes
actions add correctly across cells). The local density is
-$Z^{-1}exp(sum_e b_e t_e)$.
+Z^{-1}exp(sum_e b_e t_e).
-For a $3times2$ matrix-valued macroscopic test field $F_e$ on edges,
+For a 3times2 matrix-valued macroscopic test field F_e on edges,
define the (macroscopic) current
X_Q(F)=S^{-1}sum_e b_e,s_ecdot F_{mu(e)}(z_e),
-where $z_e$ is the edge midpoint and $mu(e)$ the edge direction. Its
+where z_e is the edge midpoint and mu(e) the edge direction. Its
second variation under infinitesimal spin twists gives the stiffness
session
@@ -695,19 +693,19 @@ session
-Cov(X_Q(F),X_Q(G)).
This is ``direct quadratic cost minus current covariance''. The
-covariance is real quadratic; on constant $3times2$ matrices, the sum
-of $|b_e|$ along any line in either spatial direction is at most
-$bS^2$, so
+covariance is real quadratic; on constant 3times2 matrices, the sum
+of |b_e| along any line in either spatial direction is at most
+bS^2, so
H_Qle bI.
-This holds for arbitrary signed $b_e$ and twists; no lower bound on
-$H_Q$ is claimed.
+This holds for arbitrary signed b_e and twists; no lower bound on
+H_Q is claimed.
-For a profile $eta$ of infinitesimal imaginary quaternions vanishing on
+For a profile eta of infinitesimal imaginary quaternions vanishing on
the boundary, Haar integration by parts produces a Ward identity. The
chief technical point is the noncommutativity: rotating one bond by a
-unit imaginary $a$ changes not only the (1,a)-plane but also the
+unit imaginary a changes not only the (1,a)-plane but also the
perpendicular imaginary plane, with a cross-product correction. The
concrete formulas are:
@@ -715,7 +713,7 @@ concrete formulas are:
D_f^2mathcal AcalA &= sum_e b_e t_e(f_y-f_x)^2
-for a one-quaternion-axis profile $a f_x$, and the exact nonabelian
+for a one-quaternion-axis profile a f_x, and the exact nonabelian
derivative of the current
D_eta s_e = t_e(eta_y-eta_x) + (eta_y+eta_x)times s_e.
@@ -725,7 +723,7 @@ cross term is the obstruction to making the argument work for
non-abelian targets.
The Ward identity expressing invariance of the measure under rotations
-by $eta$ (with $eta=0$ on the pinned boundary) has the form
+by eta (with eta=0 on the pinned boundary) has the form
E[barbar X_Q(d_SX_Q(d_Seta),X_Q(F)]
= S^{-2}Esum_e b_e t_e,(d_Seta)_ebarcdot F_e
@@ -738,22 +736,22 @@ the rescaled macroscopic variables.
## Stage A, part 2: uniform moments of the current
-We want the moments of $X_Q(F)$ uniformly in the side length and couplings,
+We want the moments of X_Q(F) uniformly in the side length and couplings,
so that local errors propagate without accumulating uncontrollably.
[Uniform current moments]
-For every fixed $1le p<∞$,
+For every fixed 1le p<∞,
|X_Q(F)|_p le C_p(1+b)log(2+S)|F|_{C^1}.
The constant is independent of pins, signs, twists, and side length.
-Suppose throughout $log(2+S)le16B$.
+Suppose throughout log(2+S)le16B.
Proof.
-Consider first a compact score $f$ (a one-axis profile supported in a
-patch, zero on the boundary). Taylor's formula for $D_f^2mathcal A$calA$ is
-bounded in absolute value by $bsum_e(f_y-f_x)^2$, so Haar invariance
+Consider first a compact score f (a one-axis profile supported in a
+patch, zero on the boundary). Taylor's formula for D_f^2mathcal AcalA is
+bounded in absolute value by bsum_e(f_y-f_x)^2, so Haar invariance
gives
Eexp(tD_fmathcal A)calA)
@@ -761,57 +759,57 @@ gives
This is the exponential-tail estimate for a compact score. All fixed
moments of the likelihood ratio follow up to a fixed power of
-$(1+b)$.
+(1+b).
Next extract a transverse current on a buffered interior patch. Take a
-scalar test $h$ on a patch of side $d$, pick a rotation axis $a$
-perpendicular to the current colour to extract; let $θ$ be an
-affine compact potential with slope $1/d$. Rotate the score by opposite
-profiles $f=tfrac12arcsin(h/N)$ with $N=C√(1+b),H$ and take the
+scalar test h on a patch of side d, pick a rotation axis a
+perpendicular to the current colour to extract; let θ be an
+affine compact potential with slope 1/d. Rotate the score by opposite
+profiles f=tfrac12arcsin(h/N) with N=C√(1+b),H and take the
difference. This gives the transverse current multiplied by
-$2sin(f_x+f_y)$, and multiplication by $N$ extracts $h$ up to $CH/d$.
-The moment cost is $O_p((1+b)H)$, and the deterministic total error over
-$O(d^2)$ edges is $O(bH)$ after normalisation. Hence
+2sin(f_x+f_y), and multiplication by N extracts h up to CH/d.
+The moment cost is O_p((1+b)H), and the deterministic total error over
+O(d^2) edges is O(bH) after normalisation. Hence
Bigl|d^{-1}sum_{e∥mu}b_e(s_e)_k,h_eBigr|_p
le C_p(1+b)H.
Finally, cover the pinned square by buffered patches at every dyadic
-scale. At scale with patch side $d$, there are $O(S/d)$ patches; each
-contributes $O((1+b)|F|_{C^1})$ per scale after normalisation.
-There are $O(log(2+S))$ scales, giving (eq:pre:2-2). Boundary
+scale. At scale with patch side d, there are O(S/d) patches; each
+contributes O((1+b)|F|_{C^1}) per scale after normalisation.
+There are O(log(2+S)) scales, giving (eq:pre:2-2). Boundary
layers contribute deterministically.
## Stage A, part 3: conditioning into cells
-Tile the pinned square by $L^2$ cells of side $R=S/L$ ($L$ is an
+Tile the pinned square by L^2 cells of side R=S/L (L is an
auxiliary integer, independent of the block factor of the same name). Let
-$omega$ be the union of the boundary rings of these cells, and let $j$
+omega be the union of the boundary rings of these cells, and let j
be the current field:
-Condition on $omega$. Cell interiors are then independent, and the
-density of $omega$ is $e^{-S_g(omega)}$ with
-$S_g(omega)=-sum_clog Z_c(omega)$. Set $z_c$ to be the lower-left
-corner of cell $c$ and use the notation
+Condition on omega. Cell interiors are then independent, and the
+density of omega is e^{-S_g(omega)} with
+S_g(omega)=-sum_clog Z_c(omega). Set z_c to be the lower-left
+corner of cell c and use the notation
[
- N_c(F)=sup_{c^+}(|F|+L^{-1}|nabla F|),quad
+ N_c(F)=sup_{c^+}(|F|+L^{-1}|nabla F|),
D_c(F)=L^{-1}N_c(nabla F),
]
-where $c^+$ is a fixed enlargement of $c$. Write $P(B)$ for a polynomial
-in $B$ whose degree and coefficients depend only on fixed numerical
-constants, not on $R,L,S$. Then ((eq:pre:2-2) and
+where c^+ is a fixed enlargement of c. Write P(B) for a polynomial
+in B whose degree and coefficients depend only on fixed numerical
+constants, not on R,L,S. Then ((eq:pre:2-2) and
(eq:pre:2-1))
- |E_c X_c(F)|&le P(B)N_c(F),qquad
+ |E_c X_c(F)|&le P(B)N_c(F),
|mathcal H_c(F,G)|calH_c(F,G)|le P(B)N_c(F)N_c(G),
|mathcal H_c(F,G)-calH_c(F,G)-barbar F_cF_c H_c G_c|
&le P(B)[D_c(F)N_c(G)+N_c(F)D_c(G)].
The second equation says: ``up to slowly varying corrections, the cell
stiffness is a constant-field stiffness'', which is a big simplification.
-Under a scalar ring rotation $u$, with $O_c(u)$ the frame rotation that
+Under a scalar ring rotation u, with O_c(u) the frame rotation that
this induces on the cell,
|O_c(u)^TH_c^uO_c(u)-H_c|le P(B)D_c(u),
@@ -826,49 +824,49 @@ Finally, the Hessian identity for the ring action:
## Stage A, part 4: augmenting with Gaussian noise
The cell stiffness matrices need not be positive; adding fake Gaussian
-noise makes the desired bound on $H_Q$ into a lower bound on the
+noise makes the desired bound on H_Q into a lower bound on the
variance of an augmented current.
-Suppose throughout the pinned cell class that $H_cle hI$, and
-$B^{-D}le hle B$. Choose $v$ and $c_0$ with
+Suppose throughout the pinned cell class that H_cle hI, and
+B^{-D}le hle B. Choose v and c_0 with
- h<vle 2B,quad c_0ge B^{-D'},quad
+ h<vle 2B, c_0ge B^{-D'},
C_c:=vI-H_c-c_0P_age B^{-D'}I,
-where $P_a$ is a projection onto the two spatial components in colour $a$.
-Let $n_c$ be independent standard six-component Gaussians and $n_{0,c}$
-independent two-component Gaussians with covariance $c_0 I$, all
+where P_a is a projection onto the two spatial components in colour a.
+Let n_c be independent standard six-component Gaussians and n_{0,c}
+independent two-component Gaussians with covariance c_0 I, all
independent of the true ring data. Define the augmented current
J(F)=j(F)+L^{-1}sum_c
(√(C_c),n_c+a,n_{0,c})cdot F_c.
-For every constant real matrix $A$,
+For every constant real matrix A,
H_Q(A,A)=v|A|^2-Cov(J(A),J(A)).
This is ``stiffness = fake-noise variance minus augmented-current
variance''. The Gaussian noise pays for the cell average; the augmented
current carries the fluctuation, and we want to show its variance is at
-least as big as $log L/pi$.
-
-Choose a smooth trigonometric cutoff $psi$ that is zero on the
-boundary and has band $O(m)$, and let $f_i=psivarphi_i$ where $varphi_i$
-are real orthonormal Fourier modes of low frequency $|k|le4K$, with
-$K+mll L$. Define $G_0$, $p$, and $theta$ at the cell anchors via the
-discrete gradient Gram matrix. Substitute $p$ into the density by an
-exact change of variables; the new density is proportional to $e^{-V(p)}$
+least as big as log L/pi.
+
+Choose a smooth trigonometric cutoff psi that is zero on the
+boundary and has band O(m), and let f_i=psivarphi_i where varphi_i
+are real orthonormal Fourier modes of low frequency |k|le4K, with
+K+mll L. Define G_0, p, and theta at the cell anchors via the
+discrete gradient Gram matrix. Substitute p into the density by an
+exact change of variables; the new density is proportional to e^{-V(p)}
with
V(p) = tfrac12 c_0 p^TG_0 p + S_g(e^{atheta}omega^0)
+tfrac12sum_c|n_c^0+√(C_c^0),a,nablatheta(z_c)/L|^2.
The change of variables is exact: the ring rotation preserves product
-Haar, the $n_0$ substitution is Gaussian, and the $n_c$ substitution is
+Haar, the n_0 substitution is Gaussian, and the n_c substitution is
an orthogonal map.
-Now imagine a family of directions $u_i$ in the potential span and set
-$∂_i = ∂/∂ t_i$ at fixed base data. These commute,
+Now imagine a family of directions u_i in the potential span and set
+∂_i = ∂/∂ t_i at fixed base data. These commute,
so we can integrate by parts normally. Then
E(sum_i∂_i^* Y_i)^2
@@ -882,59 +880,59 @@ these two sums.
## Stage A, part 5: three band estimates
-Fix any $N_*>100(1+D+D')$ and demand an accuracy $B^{-10N_*}$; make the
-underlying $log Lle16B$ dominate every polynomial in $L$ that will
-appear. Choose $M,K$ with $M=B^J$ and $K=L/B^J$ and
-$Lge B^{10J}$ for one fixed $J$. Take the cutoff
-$psi(z)=prod_mu[1-cos^{2m}(pi z_mu)]$ with $m=lceil B^4rceil$;
-then $chi=psi^2$ satisfies $intchi^2ge1-CB^{-2}$, has
-Fourier $ell^1$ norm bounded independently of $B$, and has band $O(m)$.
+Fix any N_*>100(1+D+D') and demand an accuracy B^{-10N_*}; make the
+underlying log Lle16B dominate every polynomial in L that will
+appear. Choose M,K with M=B^J and K=L/B^J and
+Lge B^{10J} for one fixed J. Take the cutoff
+psi(z)=prod_mu[1-cos^{2m}(pi z_mu)] with m=lceil B^4rceil;
+then chi=psi^2 satisfies intchi^2ge1-CB^{-2}, has
+Fourier ell^1 norm bounded independently of B, and has band O(m).
Sampling, gradients, and good event.
-For a band-$O(K)$ field $F$ with $Kll L$:
+For a band-O(K) field F with Kll L:
- |N(F)|_QED&le C|F|_2,qquad
+ |N(F)|_QED&le C|F|_2,
|D(F)|_QEDle C(K/L)|F|_2,
|d_Su-nabla u|_2&le C(K/S)^2|nabla u|_2.
-Here $N$ and $D$ are the discrete cell norms. The gradient subspace
-projection $c_0^{-1}$ also implies the Gaussian bound
+Here N and D are the discrete cell norms. The gradient subspace
+projection c_0^{-1} also implies the Gaussian bound
rho := L^{-1}|nablatheta|_∞+L^{-2}|nabla^2theta|_∞
le B^{-n}
-on an event of probability $e^{-B^3}$, for every fixed $n$, once $J$ is
+on an event of probability e^{-B^3}, for every fixed n, once J is
large.
Hessian estimate.
-Replace $d_Su$ by $nabla u(z_c)$ using (eq:pre:3-2), (eq:pre:5-1),
+Replace d_Su by nabla u(z_c) using (eq:pre:3-2), (eq:pre:5-1),
(eq:pre:5-2); conjugate each cell stiffness back using (eq:pre:3-3).
-The total error is $P(B)[K/L+(K/S)^2+rho]|nabla u|_2^2$, and the
-remaining matrix sum is exactly $H_c^0+C_c^0+c_0P_a=vI$. So
+The total error is P(B)[K/L+(K/S)^2+rho]|nabla u|_2^2, and the
+remaining matrix sum is exactly H_c^0+C_c^0+c_0P_a=vI. So
∂_u^2Vle(v+B^{-N})|nabla u|_2^2
-on the good event. Everywhere else the bound holds with $P(B)$
-replacing $v+B^{-N}$; the exceptional event has probability $e^{-B^3}$,
-so it contributes an inverse power of $B$ to any finite moment.
+on the good event. Everywhere else the bound holds with P(B)
+replacing v+B^{-N}; the exceptional event has probability e^{-B^3},
+so it contributes an inverse power of B to any finite moment.
Nonabelian derivative estimate.
-For at most $CK^2$ real directions $u_i$ satisfying
-$|nablasum_i d_iu_i|_2le C|d|$ and $N_c(nabla u_i)le C$
-((eq:pre:5-6)), and for a real band-$O(K)$ $F$ perpendicular to
-the colour $a$, then
+For at most CK^2 real directions u_i satisfying
+|nablasum_i d_iu_i|_2le C|d| and N_c(nabla u_i)le C
+((eq:pre:5-6)), and for a real band-O(K) F perpendicular to
+the colour a, then
Esum_i|∂_i J(F)-2J(u_iF')|^2
le B^{-N}|F|_2^2,
- qquad F'=-atimes F.
+ F'=-atimes F.
-The proof isolates the two parts of $2J(uF')$: the anchor-conjugation
-part (bounded by $P(B)rho$) and the non-anchor square-root derivative
-part (bounded by $P(B)K^2/L^2$ after summing). The deterministic noise
-drift and the true $H_c$ cancel via
-$H_c+C_c=vI-c_0P_a$ and $Fperp a$. There are exactly four sources of
+The proof isolates the two parts of 2J(uF'): the anchor-conjugation
+part (bounded by P(B)rho) and the non-anchor square-root derivative
+part (bounded by P(B)K^2/L^2 after summing). The deterministic noise
+drift and the true H_c cancel via
+H_c+C_c=vI-c_0P_a and Fperp a. There are exactly four sources of
error: contact (Ward identity remainder), anchor conjugation, deterministic
noise drift, and non-anchor square-root derivative. Equation
(eq:pre:5-7) is the sum of their bounds.
@@ -944,10 +942,10 @@ For the uncentred defect
[
mathcal D(F,G)=vL^calD(F,G)=vL^{-2}sum_cbarbar F_cF_ccdot G_c-EJ(F)barJ(G),
]
-and for a compact $eta$ of band $O(K)$ with a fixed Fourier symbol,
+and for a compact eta of band O(K) with a fixed Fourier symbol,
apply the raw Ward identity (eq:pre:1-5) before replacing the
gradient by its leading symbol. The contact term cancels the raw defect;
-the bracket term gives a $C^1$-norm bound $P(B)/|k|$. Hence
+the bracket term gives a C^1-norm bound P(B)/|k|. Hence
|mathcal D(d_ScalD(d_Seta,G)|le P(B)/|k|+B^{-N}.
@@ -957,41 +955,41 @@ this is the critical sequencing point.
## Stage A, part 6: the two-colour closure
Close a two-colour estimate in two steps. For perpendicular colours
-$r,r'$, and for frequencies $2Mle|k_p|le4M$, put
-$b_p=|k_p|^{-2}$ and $G_p^r=r,d_S(psitildevarphi_p/|k_p|)$. The
+r,r', and for frequencies 2Mle|k_p|le4M, put
+b_p=|k_p|^{-2} and G_p^r=r,d_S(psitildevarphi_p/|k_p|). The
Ward bound (eq:pre:5-14) applied to the two annuli in the two
colours produces a closed inequality
[
sum_p b_pE D_p^2
- le CvU^r+CT_wU^{r'}+O(B^{-N_*}),qquad
+ le CvU^r+CT_wU^{r'}+O(B^{-N_*}),
U^c=sum_pb_psum_iE J(G_p^cu_i)^2,
]
with constants independent of the number of directions. Interchange
-$r,r'$ and note $U^{-c}=U^c$, so neither colour is assumed when the other
-is estimated. Let $U=max(U^r,U^{r'})$. The signal lower bound
-$U^2le Cv(v+T_w)U+CvB^{-N_*}$ gives, after solving the quadratic and
-using $vge B^{-D}$,
+r,r' and note U^{-c}=U^c, so neither colour is assumed when the other
+is estimated. Let U=max(U^r,U^{r'}). The signal lower bound
+U^2le Cv(v+T_w)U+CvB^{-N_*} gives, after solving the quadratic and
+using vge B^{-D},
U^r+U^{r'}le Cv(v+T_w).
## Stage A, part 7: extracting all spatial derivatives
-At a fixed total complex frequency $n=k_p+k_i$, $|n|le4K_1$, we
-take $w(n-k_p)ge c(M+|n|)^{-2}$ uniformly in the annulus. Write
-$psi_mu^+(z)=tfrac12[psi(z+e_mu/(2S))+psi(z-e_mu/(2S))]$. The
-exact gradient factor in spatial direction $mu$ is
+At a fixed total complex frequency n=k_p+k_i, |n|le4K_1, we
+take w(n-k_p)ge c(M+|n|)^{-2} uniformly in the annulus. Write
+psi_mu^+(z)=tfrac12[psi(z+e_mu/(2S))+psi(z-e_mu/(2S))]. The
+exact gradient factor in spatial direction mu is
[
widehatwidehat k_k_{pmu}sinc(k_{pmu}/(2S))psi_mu^+
+frac{Scos(k_{pmu}/(2S))}{i|k_p|}
[psi(z+e_mu/(2S))-psi(z-e_mu/(2S))].
]
-Under $k_pmapsto-k_p$ the first term changes sign and the second does
+Under k_pmapsto-k_p the first term changes sign and the second does
not. Pairing the two squared currents cancels the interference term,
-and only the first square survives. Coordinate reflections of $k_p$
+and only the first square survives. Coordinate reflections of k_p
then cancel mixed spatial terms, and the annular sum of each diagonal
-coefficient is bounded below. Replacing $psipsi_mu^+$ by $psi^2$
-has Fourier error $P(B)S^{-2}$, which the crude moment bound makes
+coefficient is bounded below. Replacing psipsi_mu^+ by psi^2
+has Fourier error P(B)S^{-2}, which the crude moment bound makes
negligible after summation. This yields the weighted auxiliary estimate
W^r:=sum_{|n|le4K_1}(M+|n|)^{-2}sum_{mu=1}^2
@@ -1005,13 +1003,13 @@ Notice the boxed estimate is not assumed in any of the previous steps.
The identity (eq:pre:4-5) has a crossed derivative sum, and
(eq:pre:5-7) turns it into a bilinear pairing of the derivative
array with itself, up to a summed-square error. The complex-basis
-convention: for a real matrix $D$ transformed by a unitary matrix $U$ on
-each index, $D'=UDU^T$, and
-$sum_{ij}D'_{ij}overlineoverline{D'_D'_{ji}}=Tr(D^2)$ because $U^Toverlineoverline U=I$.U=I$.
+convention: for a real matrix D transformed by a unitary matrix U on
+each index, D'=UDU^T, and
+sum_{ij}D'_{ij}overlineoverline{D'_D'_{ji}}=Tr(D^2) because U^Toverlineoverline U=I.U=I.
So the crossed pairing is a sesquilinear pairing of the two arrays at the
-same total frequency $n=k+l$, not a bilinear pairing in the complex
-basis without conjugation. For $n$ in the main annulus, let
-$H_0(n)=sum_{k+l=n}|k|^{-2}|l|^{-2}$. Elementary annular summation gives
+same total frequency n=k+l, not a bilinear pairing in the complex
+basis without conjugation. For n in the main annulus, let
+H_0(n)=sum_{k+l=n}|k|^{-2}|l|^{-2}. Elementary annular summation gives
H_0(n)le C(M+|n|)^{-2}.
@@ -1021,60 +1019,60 @@ term satisfies
|sum_{ij}E(∂_jY_i)(∂_iY_j)|
le C[vT^2+v(v+T_w)].
-The $vT^2$ piece is the direct cost of the crossed term; the $v(v+T_w)$
+The vT^2 piece is the direct cost of the crossed term; the v(v+T_w)
piece is the paired signal.
## Stage A, part 9: closing the stiffness iteration
-Fix any sufficiently large $D$, then $D'>D$, and fix $J$ to make all
+Fix any sufficiently large D, then D'>D, and fix J to make all
preceding errors as small as required. Start with the uniform upper
-bound $H_1le bI$. If $ble B^{-D}$, no step is needed. Otherwise run
+bound H_1le bI. If ble B^{-D}, no step is needed. Otherwise run
the iteration over two regimes, a depletion regime
-$hge Alog B$ and a shallow regime $hle Clog B$.
+hge Alog B and a shallow regime hle Clog B.
+
-smallskip
-noindentDepletion regime. For $hge Alog B$, set
+noindentDepletion regime. For hge Alog B, set
[
- v=h+1,quad c_0=tfrac12,quad
+ v=h+1, c_0=tfrac12,
L=bigllceil e^{√(hlog B)}bigrrceil,
]
-with $A$ large compared with $J^2$ and the constants in
-(eq:pre:7-9). Then $Lge B^{10J}$, $log Lle h$, and
+with A large compared with J^2 and the constants in
+(eq:pre:7-9). Then Lge B^{10J}, log Lle h, and
(eq:pre:4-1), (eq:pre:4-3), (eq:pre:7-9) decrease the
stiffness bound by
[
delta hge(log L)/pi-Clog B
ge c√(hlog B).
]
-There are at most $C√(B/log B)$ such steps. Summing
-$log Llepidelta h+Clog B$ over these steps gives a total logarithmic
+There are at most C√(B/log B) such steps. Summing
+log Llepidelta h+Clog B over these steps gives a total logarithmic
side length at most
pi b+C_D√(Blog B).
-This is what produces the $pibeta$ leading term of $Xi(beta)$ in
+This is what produces the pibeta leading term of Xi(beta) in
(eq:pre:Xi).
-smallskip
-noindentShallow regime. Stop at $hle Clog B$ and take a
-small fixed $delta<c_{rm sh}/4$ with
+
+noindentShallow regime. Stop at hle Clog B and take a
+small fixed delta<c_{rm sh}/4 with
[
- v=h+deltamin(h,1),quad c_0=tfrac12deltamin(h,1),quad
+ v=h+deltamin(h,1), c_0=tfrac12deltamin(h,1),
L=lceil B^{11J}rceil.
]
The same bookkeeping bounds the shallow remainder by an additional
-exponential in $C√(Blog B)$, so the total logarithmic side length
-is exactly the $Xi(beta)$ in (eq:pre:Xi) up to constants.
+exponential in C√(Blog B), so the total logarithmic side length
+is exactly the Xi(beta) in (eq:pre:Xi) up to constants.
## Stage B: coupling of orientations via cheap rotations
-With stiffness reduced to $B^{-D}$ on a square of any side up to
-$Xi(beta)$, the probability that a fixed local spin rotation costs more
-than $t$ in action is at most $e^{-ct^2/H_j}$, so most rotations are
+With stiffness reduced to B^{-D} on a square of any side up to
+Xi(beta), the probability that a fixed local spin rotation costs more
+than t in action is at most e^{-ct^2/H_j}, so most rotations are
cheap. Coarse orientation variables are then coupled via a sequence of
local rotations that erase distant boundary disagreements one cell at a
time; each disagreement step fails with probability at most
-$e^{-ε}$ for a fixed $ε>0$, uniformly in side length and
+e^{-ε} for a fixed ε>0, uniformly in side length and
couplings.
## Stage C: disagreement percolation
@@ -1082,25 +1080,23 @@ couplings.
Seed the coupling to enforce agreement on a fixed-width boundary layer.
Iteratively reveal the configuration at the boundary of the current
agreement region, coupling the two marginal laws by their common
-$e^{-ε}$ Haar component and always allowing failure at bad
+e^{-ε} Haar component and always allowing failure at bad
vertices. Then:
-
-
-item Every successful vertex is assigned the same value in both laws.
-item Any remaining disagreement is separated from the pinned boundary
+ Every successful vertex is assigned the same value in both laws.
+ Any remaining disagreement is separated from the pinned boundary
by a chain of failures.
-item A disagreement which propagates distance $d$ from the pinned
- boundary must pass through a chain of at least $cd$ failures.
+ A disagreement which propagates distance d from the pinned
+ boundary must pass through a chain of at least cd failures.
A union bound over such chains gives the following uniform bound: for
-every bounded original observable $O$ on a finite set $A$ whose average
+every bounded original observable O on a finite set A whose average
under a common left rotation is identically zero,
|E O|le C|O|_∞|A|,e^{-cd_{rm pin}/s},
-with $d_{rm pin}$ the distance from $A$ to the pinned boundary and $s$
-a microscopic constant. Extending to pairs $O(q_x)O'(q_y)$ in
+with d_{rm pin} the distance from A to the pinned boundary and s
+a microscopic constant. Extending to pairs O(q_x)O'(q_y) in
sufficiently large unpinned tori gives
|E O(q_x)O'(q_y)|
@@ -1115,139 +1111,139 @@ signed couplings survive because no association is used.
So far the decay is proved only for observables with zero common-left
average. To extend to all observables, use three structural reductions.
-smallskip
+
noindentQuaternions as four circles. A unit quaternion is a
-pair of unit complex numbers: unit quaternions are $S^3$, identified with
-$SU(2)$, which is a double cover of $SO(3)$. In the
+pair of unit complex numbers: unit quaternions are S^3, identified with
+SU(2), which is a double cover of SO(3). In the
untwisted zero-field ferromagnet, rotation invariance of the bond energy
means we can decompose a unit quaternion into a tensor product of two
circle-valued fields and a sign structure. Working this out concretely
gives a decomposition into
-$XY$ spin variables plus a discrete ($Z_2timesZ_2$) variable;
+XY spin variables plus a discrete (Z_2timesZ_2) variable;
each subsystem is genuinely ferromagnetic once the magnitudes are returned
to the weight.
-smallskip
+
noindentConditional ferromagnets. Condition on the sign
variables. The two resulting XY subsystems are independent ferromagnetic
-$XY$ models. Condition further on the $XY$-spin magnitudes
+XY models. Condition further on the XY-spin magnitudes
and phases; the remaining sign degrees of freedom form independent
ferromagnetic Ising models. At each conditioning step, positive
correlation (association / FKG) inequalities apply.
-smallskip
+
noindentAssociation + Ginibre. For a finite collection of
-associated scalar random variables $X_i$ with finite second
-moments, and functions $F,G$ with individual Lipschitz constants
-$ell_{F,i},ell_{G,j}$, one has
+associated scalar random variables X_i with finite second
+moments, and functions F,G with individual Lipschitz constants
+ell_{F,i},ell_{G,j}, one has
|Cov(F,G)|lesum_{ij}ell_{F,i}ell_{G,j}Cov(X_i,X_j).
-(The proof: set $U=sum_iell_{F,i}X_i$, $V=sum_jell_{G,j}X_j$; the
-functions $Upm F,Vpm G$ are increasing, so their covariance is
+(The proof: set U=sum_iell_{F,i}X_i, V=sum_jell_{G,j}X_j; the
+functions Upm F,Vpm G are increasing, so their covariance is
non-negative. Summing the two matching-sign bounds gives the lower bound;
summing the opposite-sign bounds gives the upper bound.)
-The $O(4)$ plane-rotator form of Ginibre's inequality
+The O(4) plane-rotator form of Ginibre's inequality
([cite]) applies here. Two-replica substitution
-$theta=varphi+psi,theta'=varphi-psi$ is an onto homomorphism of
+theta=varphi+psi,theta'=varphi-psi is an onto homomorphism of
the product tori that pushes product Haar to product Haar; expanded in
-$varphi,psi$, each term is a difference of sines with non-negative
+varphi,psi, each term is a difference of sines with non-negative
coefficient. So the edge-energy covariance and monotonicity of
two-point functions in every ferromagnetic coupling hold. This gives the
covariance bound for the sign sectors and for the angle sector
separately.
-smallskip
+
noindentA quantitative arbitrary-observable bound inside XY.
Split an XY spin as
-$(tau^1cosgamma,tau^2singamma)$, $gammain[0,pi/2]$. Given the
+(tau^1cosgamma,tau^2singamma), gammain[0,pi/2]. Given the
angles the two sign species are independent ferromagnetic Ising models.
-The $gamma$ marginal is associated by a mixed-derivative argument.
-Therefore, with $c_{ij}=inner{s_i}{s_j}$ the XY two-point function,
+The gamma marginal is associated by a mixed-derivative argument.
+Therefore, with c_{ij}=inner{s_i}{s_j} the XY two-point function,
- E[tau_i^htau_j^h]le Cc_{ij}^{1/3},qquad
+ E[tau_i^htau_j^h]le Cc_{ij}^{1/3},
|Cov(gamma_i,gamma_j)|le Cc_{ij}^{1/3}.
-smallskip
-noindentUnconditional observables. Bound original $O$ in
-terms of its conditional expectations $f(gamma)$ and $g(gamma)$, then
+
+noindentUnconditional observables. Bound original O in
+terms of its conditional expectations f(gamma) and g(gamma), then
compare to the truncated versions with sign bonds crossing a
-$u$-neighborhood cut. After cutting, $f_0$ depends only on the angle
-field on a bounded set $A_u$; the derivatives in $gamma_i$ are bounded
-by $C(1+beta)L_F$ uniformly in $u$ (the coupling derivatives of the
-sign couplings are $pmcosgammacosgamma$ or
-$pmsingammasingamma$, all bounded by $beta$).
-Applying (eq:pre:11-1) and (eq:pre:11-2) to $f_0,g_0$, then
-averaging over $u$ via $|E_alpha c^t_{ij}|le Ce^{-cd(i,j)/s}$ and
-$Jensen$ gives
+u-neighborhood cut. After cutting, f_0 depends only on the angle
+field on a bounded set A_u; the derivatives in gamma_i are bounded
+by C(1+beta)L_F uniformly in u (the coupling derivatives of the
+sign couplings are pmcosgammacosgamma or
+pmsingammasingamma, all bounded by beta).
+Applying (eq:pre:11-1) and (eq:pre:11-2) to f_0,g_0, then
+averaging over u via |E_alpha c^t_{ij}|le Ce^{-cd(i,j)/s} and
+Jensen gives
|Cov(F,G)|
le C(1+beta+|A|+|A'|+d)^{20}L_FL_Ge^{-cd/s},
-uniformly on unpinned tori with every $0lebeta_elebeta$, including
+uniformly on unpinned tori with every 0lebeta_elebeta, including
intermediate bond weakenings.
-noindent
+
The resulting global bound (eq:pre:11-5) in the previous section is
-the $Delta_beta$ bound: it holds for every pair $F,G$ on such tori.
+the Delta_beta bound: it holds for every pair F,G on such tori.
-## Stage D, part 2: $gap$ from a spin gap within XY
+## Stage D, part 2: gap from a spin gap within XY
The bound (eq:pre:11-4) plus reflection positivity
(eq:pre:11-5) gives the true finite-volume full gap within the
reflection-positive Hilbert space. We sketch the passage, then turn to
-$Delta_beta$.
+Delta_beta.
Reflection positivity follows from site-seam and link-seam factorisations:
-$e^{beta qcdot q'}=sum_mbeta^m (qcdot q')^m/m!$ is a sum of positive
+e^{beta qcdot q'}=sum_mbeta^m (qcdot q')^m/m! is a sum of positive
tensor products. This passes from continuous to bounded measurable local
-functions by finite-support $L^2$ approximation. Then
-$inner{F}{F}=E[(Θ F)F]ge0$ on bounded local functions, the
+functions by finite-support L^2 approximation. Then
+inner{F}{F}=E[(Θ F)F]ge0 on bounded local functions, the
nullspace is quotiented out, and the completion is a Hilbert space. The
-one-step translation $T$ acts as a self-adjoint contraction.
+one-step translation T acts as a self-adjoint contraction.
+
-smallskip
-Actually, bounding $p(n,w)$ requires the full-gap spectrum in
-both coordinate directions, because $Delta_beta(n)=2p(n,n)+q(2n,n)$
-and $q$ involves the full-gap spectrum in the space direction too. This is
+Actually, bounding p(n,w) requires the full-gap spectrum in
+both coordinate directions, because Delta_beta(n)=2p(n,n)+q(2n,n)
+and q involves the full-gap spectrum in the space direction too. This is
why the argument in this section carries every local observable: the
-$SO(4)$-invariant bond energies are observables whose
+SO(4)-invariant bond energies are observables whose
autocorrelations must also be controlled. Once (eq:pre:11-4) holds
for all local observables and (eq:pre:11-5) gives
reflection-positive transfer, the full-gap argument of Step 1 applies as
stated.
-## Stage D, part 3: $Delta_beta$ bound and the proof of the theorem
+## Stage D, part 3: Delta_beta bound and the proof of the theorem
-Let $e_t(n)=E_{t;n,n}(q_0cdot q_{e_1})$ and let $b$ be the maximum
+Let e_t(n)=E_{t;n,n}(q_0cdot q_{e_1}) and let b be the maximum
coupling. Use (eq:pre:11-4) with the actual maximum coupling, not
-a varying threshold. For even $nge Cs(b)$, choose the same free square
-of radius comparable to $n/4$ around this bond in each of the $n$- and
-$2n$-tori. Cut its crossing bonds by interpolation. There are $O(n)$
-such edges; their distance from the fixed bond is at least $cn$,
-their coefficients are at most $b$, and the edge observables have
+a varying threshold. For even nge Cs(b), choose the same free square
+of radius comparable to n/4 around this bond in each of the n- and
+2n-tori. Cut its crossing bonds by interpolation. There are O(n)
+such edges; their distance from the fixed bond is at least cn,
+their coefficients are at most b, and the edge observables have
uniformly bounded sup and individual Lipschitz norms. Therefore
|e_t(n)-e_t(2n)|le C(1+beta+n)^Pe^{-cn/s(b)},
- quad 0le tlebeta.
+ 0le tlebeta.
-Differentiate the normalised partition functions. There are $2n^2$ bonds
-on the $n$-square torus and $8n^2$ on the $2n$-square torus, so
+Differentiate the normalised partition functions. There are 2n^2 bonds
+on the n-square torus and 8n^2 on the 2n-square torus, so
- (d)/(dt)Delta_t(n)=8n^2[e_t(n)-e_t(2n)],quad Delta_0(n)=0.
+ (d)/(dt)Delta_t(n)=8n^2[e_t(n)-e_t(2n)], Delta_0(n)=0.
-Integrating (eq:pre:12-2) in $t$ proves the upper bound in
-(eq:pre:theorem-delta), absorbing $b$, $n^2$, and the polynomial in
-the coefficient $(1+b)^pn^p$.
+Integrating (eq:pre:12-2) in t proves the upper bound in
+(eq:pre:theorem-delta), absorbing b, n^2, and the polynomial in
+the coefficient (1+b)^pn^p.
-Non-negativity of $Delta_t(n)$ on the right side of
+Non-negativity of Delta_t(n) on the right side of
(eq:pre:theorem-delta) follows from purity; we already proved it in
-Lemma~lem:Delta-forms. Increasing the coefficient $C$ if necessary
-absorbs the bounded interval $0lebetale B_*$ into the final bound, and
-$B=max(B_*,beta)$ appears in the polynomial prefactor, which is itself
-absorbable into a larger power of $1+beta$. This completes the proof of
+Lemma~lem:Delta-forms. Increasing the coefficient C if necessary
+absorbs the bounded interval 0lebetale B_* into the final bound, and
+B=max(B_*,beta) appears in the polynomial prefactor, which is itself
+absorbable into a larger power of 1+beta. This completes the proof of
Theorem~thm:preliminary.
%---------------------------------------------------------------------
@@ -1265,65 +1261,65 @@ contains the Fourier bookkeeping; here we record the structural facts.
## Setup: precision and observation
-At a scale $h$ we are given a precision $K_h$ on the fine lattice (a
-positive-definite edge operator: $K_h=d^*B_hd$ with $B_h$ a positive
-edge lift, normalised $B_h(0)=I$). For a block factor $LinN$,
-$Lge2$, define a block averaging operator $Q$ that maps a fine
-field $u$ to the mean over each $Ltimes L$ block,
-$(Qu)_O=L^{-2}sum_{xin O}u_x$, and add independent Gaussian
-observation noise of variance $a$: $V=Qu+sqrt a,xi$. The coarse
+At a scale h we are given a precision K_h on the fine lattice (a
+positive-definite edge operator: K_h=d^*B_hd with B_h a positive
+edge lift, normalised B_h(0)=I). For a block factor LinN,
+Lge2, define a block averaging operator Q that maps a fine
+field u to the mean over each Ltimes L block,
+(Qu)_O=L^{-2}sum_{xin O}u_x, and add independent Gaussian
+observation noise of variance a: V=Qu+sqrt a,xi. The coarse
precision is
K_{h+1} = (Q^*K_h^{-1}Q + a, I)^{-1}
-(fine-scale subtraction).
-In practice one writes $K_{h+1}=QK_hQ^*-...$ with the Schur identity
-below making the subtraction explicit. Constants $C_L,c_L>0$ may
-depend on $L$ but not on the scale.
+In practice one writes K_{h+1}=QK_hQ^*-... with the Schur identity
+below making the subtraction explicit. Constants C_L,c_L>0 may
+depend on L but not on the scale.
## Conditional mean and covariance
[Free conditional mean]
-Let $k_l$ and $p$ be the Fourier momenta of the fine and coarse lattices,
-$d_mu(k_l)=1-e^{-ik_{l,mu}}$ the fine difference symbol,
-$d'(p)=1-e^{-ip_mu}$ the coarse one, $b(k_l)$ the block Fourier symbol,
-$K(k_l)=K_h$ evaluated. The conditional mean operator has kernel
+Let k_l and p be the Fourier momenta of the fine and coarse lattices,
+d_mu(k_l)=1-e^{-ik_{l,mu}} the fine difference symbol,
+d'(p)=1-e^{-ip_mu} the coarse one, b(k_l) the block Fourier symbol,
+K(k_l)=K_h evaluated. The conditional mean operator has kernel
[
P_l=L^{-2}K(k_l)^{-1}b(k_l)^*K'(p),
]
-which obeys $KP=Q^*K'$ and $a(I-QP)=K'$, preserves constants and
-centred affine fields, and, for every fixed integer $rge0$,
+which obeys KP=Q^*K' and a(I-QP)=K', preserves constants and
+centred affine fields, and, for every fixed integer rge0,
sup_xsum_y e^{c|x/L-y|}|nabla_f^rP(x,y)|le C_rL^{-r}.
-Here $nabla_f$ is a nearest-neighbor difference on the fine lattice.
+Here nabla_f is a nearest-neighbor difference on the fine lattice.
The conditional covariance
C_h=(K_h+aQ^*Q)^{-1}
has exponentially decaying kernels in coarse units, with constants
-allowed to depend on the fixed $L$. All these assertions retain a
-factor $L^{-h}$ in history differences.
+allowed to depend on the fixed L. All these assertions retain a
+factor L^{-h} in history differences.
Proof.
The normal equations for the conditional mean give
(eq:free:P) and both stated identities. On a nonprincipal alias
the smooth block bump gives arbitrary fixed inverse-power decay in
-$l$; on the principal alias the apparent pole is removable:
+l; on the principal alias the apparent pole is removable:
[
P_0=b(k_0)^{-1}Bigl[1-(a^{-1}+L^{-2}sum_{lne0}
|b(k_l)|^2K(k_l)^{-1})K'Bigr].
]
-At $p=0$ the other aliases vanish to second order, and the zero and
+At p=0 the other aliases vanish to second order, and the zero and
first jets of this expression give preservation of constants and
affine fields. Each fine difference contributes
-$|d_mu(k_l)|le C(1+|l|)/L$; with bump decay larger than the required
+|d_mu(k_l)|le C(1+|l|)/L; with bump decay larger than the required
number of differences plus the summation dimension, Fourier inversion
on a smaller coarse strip gives (eq:free:derivative). Invert
-$K+aQ^*Q$ in aliases: the nonprincipal diagonal block is invertible
+K+aQ^*Q in aliases: the nonprincipal diagonal block is invertible
with the preceding bounds; its rank-one update has positive
-denominator; the principal Schur complement multiplied by $L^2$ is
+denominator; the principal Schur complement multiplied by L^2 is
bounded below by the principal bump. Fourier inversion proves the
covariance localization.
@@ -1331,77 +1327,77 @@ covariance localization.
## The Schur energy decomposition
The conditional mean gives an exact algebraic identity for the free
-energy. For real ambient fields $u,V$ on the two lattices,
+energy. For real ambient fields u,V on the two lattices,
tfrac12inner{u}{K_hu}+tfrac a2|V-Qu|^2
=tfrac12inner{V}{K_{h+1}V}
+tfrac12inner{u-P_hV}{C_h^{-1}(u-P_hV)}.
-Indeed $C_h^{-1}P_h=aQ^*$ and $a(I-QP_h)=K_{h+1}$ make the cross and
-constant terms agree after expansion. Although $K_h$ vanishes on
-constants, $C_h^{-1}=K_h+aQ^*Q$ is strictly positive, because $Q$
+Indeed C_h^{-1}P_h=aQ^* and a(I-QP_h)=K_{h+1} make the cross and
+constant terms agree after expansion. Although K_h vanishes on
+constants, C_h^{-1}=K_h+aQ^*Q is strictly positive, because Q
preserves constants. This identity shows that the correct fluctuation
-field is $u-P_hV$, and it has a strictly positive precision quadratic
+field is u-P_hV, and it has a strictly positive precision quadratic
near rough regions. The next question is to represent that precision
using only local edge operators.
## The edge intertwiner
-The edge representation of $K_h=d^*B_hd$ is not unique: adding a curl
-correction to $B_h$ leaves $d^*B_hd$ unchanged. We use this freedom to
+The edge representation of K_h=d^*B_hd is not unique: adding a curl
+correction to B_h leaves d^*B_hd unchanged. We use this freedom to
find a positive completion whose fine-output rows obey a strong
-$L^{-1}$ fine-difference bound, so that nothing accumulates across the
-$L^2$ fine anchors per coarse one.
+L^{-1} fine-difference bound, so that nothing accumulates across the
+L^2 fine anchors per coarse one.
[Linear block stack]
-There is a coarse-edge to fine-edge map $T$ satisfying
+There is a coarse-edge to fine-edge map T satisfying
- T^*d=d'Q,qquad TB'd'=BdP,
+ T^*d=d'Q, TB'd'=BdP,
-whose aliases are bounded by $C_JL^{-1}< l>^{-J}$, and whose
-rows, columns, and fine-difference bounds depend only on $L$. One can
-find a positive scalar $B_2$ with an inequality
+whose aliases are bounded by C_JL^{-1}< l>^{-J}, and whose
+rows, columns, and fine-difference bounds depend only on L. One can
+find a positive scalar B_2 with an inequality
[
widetilde B'U_*widetilde B'=mathsf c'^*B_2sfc'^*B_2mathsf c'sfc'
]
-such that the Schur complement $U_*$ is nonnegative; consequently there
+such that the Schur complement U_* is nonnegative; consequently there
exists an ambient-field identity
U_*+mathcal M^*calM^*mathcal M+calM+mathcal N^*calN^*mathcal N=I,calN=I,
-with $mathcal M,calM,mathcal N$calN$ defined in terms of $T,B',mathsf A_s$sfA_s$
-and $S'=mathsf A'(sfA'(widetilde B')^{-1}(mathsf A')^*$sfA')^*$ an orthogonal
+with mathcal M,calM,mathcal NcalN defined in terms of T,B',mathsf A_ssfA_s
+and S'=mathsf A'(sfA'(widetilde B')^{-1}(mathsf A')^*sfA')^* an orthogonal
projection.
Proof.
-Start with $T_{0,l}=L^{-2}diag(omega_{mu,l}^*)b(k_l)^*$
-with $omega_{mu,l}=sum_{j=0}^{L-1}e^{ijk_{l,mu}}$. The identity
-$omega_{mu,l}d_mu(k_l)=d'_mu(p)$ proves the first equation. Put
-$R_l=B_ld_lP_l-T_{0,l}B'd'$; its divergence vanishes by $KP=Q^*K'$, and
-its aliases have bounds $C_sL^{-1}< l>^{-s}$ for every
-required fixed $s$, with $R_0$ vanishing through first order at $p=0$.
+Start with T_{0,l}=L^{-2}diag(omega_{mu,l}^*)b(k_l)^*
+with omega_{mu,l}=sum_{j=0}^{L-1}e^{ijk_{l,mu}}. The identity
+omega_{mu,l}d_mu(k_l)=d'_mu(p) proves the first equation. Put
+R_l=B_ld_lP_l-T_{0,l}B'd'; its divergence vanishes by KP=Q^*K', and
+its aliases have bounds C_sL^{-1}< l>^{-s} for every
+required fixed s, with R_0 vanishing through first order at p=0.
Apply the analytic alias-division lemma
(Lemma~app:alias-division) to produce
-$R_l=mathsf c_l^*sfc_l^*gamma_l$ and $H_ld'=gamma_l$ with
-$|gamma_l|+|H_l|le C_s< l>^{-s}$. Set
-$T_l=T_{0,l}+mathsf c_l^*H_l(B')^sfc_l^*H_l(B')^{-1}$, proving
+R_l=mathsf c_l^*sfc_l^*gamma_l and H_ld'=gamma_l with
+|gamma_l|+|H_l|le C_s< l>^{-s}. Set
+T_l=T_{0,l}+mathsf c_l^*H_l(B')^sfc_l^*H_l(B')^{-1}, proving
(eq:free:Tidentities) with the same alias bounds. Fourier
inversion gives rows, columns, and fine differences.
-For positivity, add a scalar correction to $widetilde B^{-1}$:
+For positivity, add a scalar correction to widetilde B^{-1}:
widetilde B^{-1}=B^{-1}
+lambda B^{-1}mathsf c^*sfc^*mathsf csfc B^{-1}
-at $B$ and at $B'$. This leaves $widetilde Bd=Bd$, hence $K$, unchanged.
-Choose $lambdage2(C+1)$ and then $L$ large enough that
-$widetilde B'U_*widetilde B'$ is positive and, after removing the
+at B and at B'. This leaves widetilde Bd=Bd, hence K, unchanged.
+Choose lambdage2(C+1) and then L large enough that
+widetilde B'U_*widetilde B' is positive and, after removing the
principal positivity, becomes
-$mathsf c'^*B_2sfc'^*B_2mathsf c'$sfc'$ with $B_2$ a positive analytic scalar.
-Positivity extends through $p=0$. Choose $c_0$ below the common lower
-bound of $widetilde B$ and set
-$ell=(widetilde B-c_0I)^{1/2}$. Taking square roots: all remain
+mathsf c'^*B_2sfc'^*B_2mathsf c'sfc' with B_2 a positive analytic scalar.
+Positivity extends through p=0. Choose c_0 below the common lower
+bound of widetilde B and set
+ell=(widetilde B-c_0I)^{1/2}. Taking square roots: all remain
uniformly analytic on a smaller strip. Define
S'&=mathsf A'(sfA'(widetilde B')^{-1}(mathsf A')^*,sfA')^*,
@@ -1411,13 +1407,13 @@ uniformly analytic on a smaller strip. Define
B_2^{1/2}mathsf c'(sfc'(widetilde B')^{-1}(mathsf A')^*z).sfA')^*z).
Substitution of (eq:free:Tidentities) proves
-(eq:free:ambient) and the cokernel property; $mathsf c'd'=0$sfc'd'=0$
-proves that $mathcal N$'scalN$'s second factor has the requested form. Using
-the definition of $U_*$ gives $mathcal M^*calM^*mathcal M+calM+mathcal N^*calN^*mathcal N=I$.calN=I$.
+(eq:free:ambient) and the cokernel property; mathsf c'd'=0sfc'd'=0
+proves that mathcal N'scalN's second factor has the requested form. Using
+the definition of U_* gives mathcal M^*calM^*mathcal M+calM+mathcal N^*calN^*mathcal N=I.calN=I.
The extra factors are uniformly summable convolution kernels, so they
-preserve the proved bounds for $T$. At depth zero take
-$widetilde B=I$ instead of the fine curl correction; this deletes a
-non-negative subtraction from $U_*$ and thus only improves positivity.
+preserve the proved bounds for T. At depth zero take
+widetilde B=I instead of the fine curl correction; this deletes a
+non-negative subtraction from U_* and thus only improves positivity.
All finite-period operators are obtained by folding these same kernels.
@@ -1430,31 +1426,31 @@ integration.
## Exponential positivity and uniform analyticity
The analytic-linear-algebra facts of this section do not degrade at
-each iteration because the constants $c_L,C_L$ are only allowed to
-depend on the fixed block size $L$ chosen once and for all.
+each iteration because the constants c_L,C_L are only allowed to
+depend on the fixed block size L chosen once and for all.
[Uniform shells and strips]
For the sampled smooth block weights used in Section~sec:rg,
-there are $delta,c,C>0$ and $L_0$ with the following properties for
-$Lge L_0$, $hge0$. The precisions $K_h$ and elliptic edge lifts
-$B_h$ are analytic on the complex $delta$-neighborhood of the momentum
+there are delta,c,C>0 and L_0 with the following properties for
+Lge L_0, hge0. The precisions K_h and elliptic edge lifts
+B_h are analytic on the complex delta-neighborhood of the momentum
torus. On the real torus,
[
- cDle K_hle CD,quad cIle B_hle CI,quad K_h=d^*B_hd,quad
+ cDle K_hle CD, cIle B_hle CI, K_h=d^*B_hd,
B_h(0)=I,
]
-where $D(p)=sum_mu4sin^2(p_mu/2)=|d|^2$. Their fixed derivatives
-are bounded, $K_h-D$ has a zero of order four at zero, and their
+where D(p)=sum_mu4sin^2(p_mu/2)=|d|^2. Their fixed derivatives
+are bounded, K_h-D has a zero of order four at zero, and their
history differences, including inverses of edge lifts, are bounded by
-$CL^{-h}$ when compared with any depth $barbar hhge h$. The same
+CL^{-h} when compared with any depth barbar hhge h. The same
assertions hold on a common smaller strip.
The proof is in Appendix~app:analytic-block; the essential point
-is the one-sided sequence of $L$-shells, so that at scale $h$ the bump
-product $W_h=prod_{j=1}^h b(xi/L^j)$ is bounded below on the real
-torus uniformly in $h$, and the nonprincipal shell contribution
-decays like $L^{-h}$.
+is the one-sided sequence of L-shells, so that at scale h the bump
+product W_h=prod_{j=1}^h b(xi/L^j) is bounded below on the real
+torus uniformly in h, and the nonprincipal shell contribution
+decays like L^{-h}.
The linear block is now fully understood. Section~sec:rg builds
the nonlinear renormalisation step on this foundation.
@@ -1477,7 +1473,7 @@ linearise around the Gaussian. Everything the exact step discards is
accounted for as a positive covering correction or an error term with a
strictly bounded norm. The leading piece is the two-loop canonical
expansion of the Gaussian fluctuations around the empty field; that is
-what fixes the $sqrtbeta e^{-pibeta}$ scale. The rest is the
+what fixes the sqrtbeta e^{-pibeta} scale. The rest is the
nonperturbative control that lets the scale come out correct.
Let me present the step in four parts.
@@ -1493,57 +1489,57 @@ Sections~sec:rg-proof and Appendices~app:analytic-block--app:rg-geometry.
## (a) Setup: interpolation, kinetic energy, positivity
-We work on a torus of side $m=L^j$ at scale $j$. The coarser torus has
-side $m'=L^{j-1}=m/L$. Both are periodic. Fine sites are $x$; coarse
-blocks are $O$ of side $L$ with $L$ odd and $ge3$. Set
-$L^2=|block|$.
+We work on a torus of side m=L^j at scale j. The coarser torus has
+side m'=L^{j-1}=m/L. Both are periodic. Fine sites are x; coarse
+blocks are O of side L with L odd and ge3. Set
+L^2=|block|.
noindentObservation kernel.
-Declare that the coarse spin $V_bin S^3$ in block $b$ is distributed
-given the fine spins ${q_x : xin b}$ as a normalised von Mises--
+Declare that the coarse spin V_bin S^3 in block b is distributed
+given the fine spins {q_x : xin b} as a normalised von Mises--
Fisher-lite kernel that will be the "observation": for a coarse spin
-$V$ and a block configuration ${U_1,...,U_{L^2}}subset S^3$ with
-normalised mean $barbar u$,u$, put
+V and a block configuration {U_1,...,U_{L^2}}subset S^3 with
+normalised mean barbar u,u, put
Q_b(Vmid u)= z_{ab}^{-1}exp!(ab,Vcdot barbar u),u),
-with $a,L^{-2}le ale 1$, $b>0$ a known coupling, $z_{ab}$ its
-normalisation. This is an $O(4)$-symmetric positive kernel; the block
-configuration enters only through the mean spin $barbar u$.u$. Conditional on
-all block interiors, the $V_b$ are independent.
+with a,L^{-2}le ale 1, b>0 a known coupling, z_{ab} its
+normalisation. This is an O(4)-symmetric positive kernel; the block
+configuration enters only through the mean spin barbar u.u. Conditional on
+all block interiors, the V_b are independent.
noindentKinetic energy.
-Given a fine configuration $q=(q_x)_{x}$, define
+Given a fine configuration q=(q_x)_{x}, define
mathcal E_calE_{h,t}(q) = sum_e S_t(r_e)
+ tfrac12sum_e m_e(q),|(ell_h dq)_e|^2,
-with $r_e=|d_eq|$, $d_eq = q_y-q_x$ for $e=xto y$, and where
+with r_e=|d_eq|, d_eq = q_y-q_x for e=xto y, and where
S_t(r)&=c_0r^2/2,&rle t,d_0tr,&r>t,
m_e(q)&=mathbf1bf1{r_fle te^{dist(e,f)/R_*}
for all f},
-with $0<d_0<c_0/16$, $R_*$ large compared with the exponential
-localisation length of $ell_h$. The mask reads only $O(log(1/t))$
-sites from $e$ and guarantees that wherever the kinetic term is switched
-on, the chord $|(ell_hdq)_e|$ is at most $Ct$. The kernels $ell_h$ are
+with 0<d_0<c_0/16, R_* large compared with the exponential
+localisation length of ell_h. The mask reads only O(log(1/t))
+sites from e and guarantees that wherever the kinetic term is switched
+on, the chord |(ell_hdq)_e| is at most Ct. The kernels ell_h are
the exponentially-localised rows of Lemma~lem:free-stack; each in
-turn depends only on fixed $L$.
+turn depends only on fixed L.
-The key positivity fact for the whole construction is: $mathcal E_calE_{h,t}$
-is a sum of non-negative terms (because $m_ein{0,1}$). This is what
+The key positivity fact for the whole construction is: mathcal E_calE_{h,t}
+is a sum of non-negative terms (because m_ein{0,1}). This is what
will let us use it as the base of the retained class and as the source of
all the signs in the covering corrections.
noindentObservation and the microscopic model.
-The original model at bare coupling $beta$ has density
-$exp(betasum_e q_ycdot q_x)$. We will treat this as the scale-$N$
+The original model at bare coupling beta has density
+exp(betasum_e q_ycdot q_x). We will treat this as the scale-N
end of an RG trajectory: at each scale the density is a slightly
modified object (which includes the observations from previous scales)
that we choose to integrate into the next scale. The observations
-already present introduce new degrees of freedom $V$ that are integrated
+already present introduce new degrees of freedom V that are integrated
out at the next step. This is an exact interchange of order of
integration because the observation kernel (eq:obs-kernel) is
positive.
@@ -1551,61 +1547,59 @@ positive.
## (b) The class of retained densities
We want the class to be:
-
-item capable of expressing the density produced by the previous step,
-item self-closed under one more RG step,
-item formally close to the two-loop canonical expansion (so that the
+ capable of expressing the density produced by the previous step,
+ self-closed under one more RG step,
+ formally close to the two-loop canonical expansion (so that the
running coupling is computed correctly).
Fix once and for all:
-
-item a block factor $L$ (odd, $ge3$) chosen later;
-item a scale-dependent coupling $b_j>0$ for the kinetic part;
-item a cutoff $0<t_j<pi$ (the kinetic smoothing scale);
-item a mask radius parameter $R_*$ (fixed);
-item microscopic constants $A=64$, $k=8$, $P_0$ (fixed after the
+ a block factor L (odd, ge3) chosen later;
+ a scale-dependent coupling b_j>0 for the kinetic part;
+ a cutoff 0<t_j<pi (the kinetic smoothing scale);
+ a mask radius parameter R_* (fixed);
+ microscopic constants A=64, k=8, P_0 (fixed after the
running-coupling bookkeeping);
-item a scaling schedule $H_0 = beta$, $H_{j+1}=H_j$, chosen to
+ a scaling schedule H_0 = beta, H_{j+1}=H_j, chosen to
contract.
[Class of density]
-Fix the parameters above and give the scale-$j$ lattice
-$Λ_j=(Z/L^jZ)^2$ on the torus of side $L^j$. A density
-$rho$ on ambient configurations $Q=(q_x)_{xinΛ_j}$ is in the
-scale-$j$ class $mathscr C_j=scrC_j=mathscr C_j(b_j,t_j,scrC_j(b_j,t_j,delta_j,h)$
+Fix the parameters above and give the scale-j lattice
+Λ_j=(Z/L^jZ)^2 on the torus of side L^j. A density
+rho on ambient configurations Q=(q_x)_{xinΛ_j} is in the
+scale-j class mathscr C_j=scrC_j=mathscr C_j(b_j,t_j,scrC_j(b_j,t_j,delta_j,h)
if:
[label=(roman*)]
-item It has the form
+ It has the form
rho(q)=e^{vm^2}expBigl{-b_jmathcal E_calE_{h,t_j}(q)
-sum_{slots}mathbf1_Xbf1_X P_X(q)
-sum_Xmathbf1_Xbf1_X f_X(q)Bigr}Xi(q),
- with $v$ a scalar (per unit volume); the second sum has strictly
- bounded support on each site; $Xi$ is a
+ with v a scalar (per unit volume); the second sum has strictly
+ bounded support on each site; Xi is a
mandatory covering gas (defined below);
-item Each canonical coefficient polynomial $P_X(q)$ is a
- finite polynomial in the relative chart $y_o(r)=q_x q_o^{-1}$
+ Each canonical coefficient polynomial P_X(q) is a
+ finite polynomial in the relative chart y_o(r)=q_x q_o^{-1}
at the anchor, invariant under simultaneous spin and coordinate
symmetries, and its Hessian vanishes on unit-affine fields at the
identity;
-item Each regular error $f_X$ is a $C^k$ function (with
- $k=8$) on the strong local domain ${r_e<4t_j:ein X}$, its
+ Each regular error f_X is a C^k function (with
+ k=8) on the strong local domain {r_e<4t_j:ein X}, its
first value and affine Hessian at the identity vanish on
bulk/nonwinding labels, and its norm satisfies
|f|_{k,j,A}:=sup_osum_{X at o}e^{As_X}[f_X]_{k,j}
- le delta_j:=H_j^{-2.05},qquad A=64,
+ le delta_j:=H_j^{-2.05}, A=64,
with the derivatives taken in the chart
- $q_x(theta)=exp(sum_itheta_i v_{i,x})q_x$ and normalised
- directions $|v_x|le t_j(1+dist(x,o))^8$;
-item Each covering label $ell$ has a bounded measurable
- weight $k_ell(q)$, a finite complete support $P_ell$, and a
- non-empty inventory $D_ell$ of bonds contained in $P_ell$,
- with $k_ellequiv0$ unless $D_ellsubset D(q)={e:r_e>t_j}$.
+ q_x(theta)=exp(sum_itheta_i v_{i,x})q_x and normalised
+ directions |v_x|le t_j(1+dist(x,o))^8;
+ Each covering label ell has a bounded measurable
+ weight k_ell(q), a finite complete support P_ell, and a
+ non-empty inventory D_ell of bonds contained in P_ell,
+ with k_ellequiv0 unless D_ellsubset D(q)={e:r_e>t_j}.
The covering gas
Xi(q)=sum_{substack{Gamma: disjoint support
@@ -1613,12 +1607,12 @@ item Each covering label $ell$ has a bounded measurable
prod_{ellinGamma}k_ell(q)
has the anchored norm bound
- $|k|_{j,A}:=sup_xsum_{ell:xin P_ell}e^{As_ell}|k_ell|_∞
- le w_j:=e^{-p_j^{1/4}}$;
-item The ambient configuration enters the labels' supports only
- through $N_{50M}$-neighbourhoods defined in
+ |k|_{j,A}:=sup_xsum_{ell:xin P_ell}e^{As_ell}|k_ell|_∞
+ le w_j:=e^{-p_j^{1/4}};
+ The ambient configuration enters the labels' supports only
+ through N_{50M}-neighbourhoods defined in
Appendix~app:rg-geometry;
-item The volume scalar $v$ and the terminal coupling $b_j$ are the
+ The volume scalar v and the terminal coupling b_j are the
same on every admitted period (translation-invariance plus
asymptotic freedom; see Section~sec:calibration).
@@ -1626,63 +1620,63 @@ item The volume scalar $v$ and the terminal coupling $b_j$ are the
The class is intentionally large: it accommodates the density of the
previous RG step and it is closed under one more. The point of the
-class is not that $rho$ has a simple form --- the form has many slots
+class is not that rho has a simple form --- the form has many slots
and one arbitrary mandatory gas --- but that all its slots are
bounded in a uniform scale-dependent way. That is what allows
contraction.
## (c) The exact integration
-The exact integration of a class-$mathscr C_j$scrC_j$ density to
-scale $j-1$ proceeds in three exact steps.
+The exact integration of a class-mathscr C_jscrC_j density to
+scale j-1 proceeds in three exact steps.
noindentStep 1: block observation.
Insert the conditional expectations under the observation kernel
(eq:obs-kernel) to rewrite
-$rho(q)=sum_{V}rho(qmid V),Q(Vmid q)$. This is an exact
+rho(q)=sum_{V}rho(qmid V),Q(Vmid q). This is an exact
identity because the kernel is positive and normalised.
noindentStep 2: exact MAP substitution.
-The exact dependence of the kinetic energy on $V$ is reorganised by
+The exact dependence of the kinetic energy on V is reorganised by
writing the kinetic part as
b_jmathcal E_calE_{h,t_j}(q)
= b_jmathcal E_calE_{h,t_j}(q)-tfrac12 a,b_j|V-Qu|^2
+tfrac12 a,b_j|V-Qu|^2,
-i.e. we insert and remove the $V$-dependent Gaussian term. Combined
+i.e. we insert and remove the V-dependent Gaussian term. Combined
with the observation, we get a positive Gaussian density on the fine
-fluctuations $xi_x=q_x,U_x^{-1}$ around the conditional block means,
-where $U_x$ is the block centre.
+fluctuations xi_x=q_x,U_x^{-1} around the conditional block means,
+where U_x is the block centre.
The residual after the conditional mean is bounded by
|q_x-q_x^{cond}|le C_L,(Mt_j)^2
-for $x$ in the block at distance at least $M$ from any core seed
-($M$ from the cutoff scale in Appendix~app:analytic-block).
+for x in the block at distance at least M from any core seed
+(M from the cutoff scale in Appendix~app:analytic-block).
noindentStep 3: Taylor to quartic + canonical map.
Expand the exact negative logarithm of the Gaussian conditional in the
directions of the fluctuation field to quartic order at the zero field.
-The bits of total homogeneous degree $ge5$ and the bits of quartic
+The bits of total homogeneous degree ge5 and the bits of quartic
order that do not fit the ``canonical coefficients'' list go into
covering corrections (positive) and regular errors (bounded by the cap
-$delta_j$).
+delta_j).
-The bits of order $le4$ are matched term by term with the canonical
+The bits of order le4 are matched term by term with the canonical
two-loop expansion, which reads
beta'=beta+alpha_h+(kappa_h)/(beta)+Delta_h,
with (matching the classical normalization of the O(4) sigma model)
- alpha_h=-gamma+O_L(rho^h),quad
+ alpha_h=-gamma+O_L(rho^h),
kappa_h=-(gamma)/(2pi)+O_L(rho^h),
-where $gamma=(log L)/pi$ and $rho<1$ a fixed constant depending only
-on $L$. Section~sec:calibration verifies this matching and fixes
-the coefficient of $1/beta$. The residual error $Delta_h$ is bounded
+where gamma=(log L)/pi and rho<1 a fixed constant depending only
+on L. Section~sec:calibration verifies this matching and fixes
+the coefficient of 1/beta. The residual error Delta_h is bounded
in Theorem~thm:rg-closure.
Step 3 is what turns the exact density into one of the retained form.
@@ -1694,12 +1688,12 @@ Definition~def:rg-class.
## (d) The closure theorem
[Exact RG step]
-Fix the parameters $L,P_0$ large and $H$ large enough so that the
+Fix the parameters L,P_0 large and H large enough so that the
hypotheses of Appendix~app:finite-choices hold. There are
-$H_j$-dependent constants $C_L,c_L>0$, and there is a scale-dependent
-data $(b_j,t_j,delta_j,h_j)$ such that, for every
-$j$-level density in $mathscr C_j(b_j,t_j,scrC_j(b_j,t_j,delta_j,h_j)$ and every
-period $m=L^j$ large enough, the exact integration to scale $j-1$ gives
+H_j-dependent constants C_L,c_L>0, and there is a scale-dependent
+data (b_j,t_j,delta_j,h_j) such that, for every
+j-level density in mathscr C_j(b_j,t_j,scrC_j(b_j,t_j,delta_j,h_j) and every
+period m=L^j large enough, the exact integration to scale j-1 gives
a density in
[
mathscr C_scrC_{j-1}(b_{j-1},t_{j-1},delta_{j-1},h_{j-1})
@@ -1707,25 +1701,25 @@ a density in
with
|b_j-b_{j-1}-alpha_{j-1}-kappa_{j-1}/b_{j-1})| le C_L b_j^{-2},
- qquad
+
frac{delta_{j}}{delta_{j-1}} ge H_j^{-1/20}.
-Here $alpha_{j-1},kappa_{j-1}$ are as in (eq:rg:kicks), with
-error $O_L(rho^{j-1})$.
+Here alpha_{j-1},kappa_{j-1} are as in (eq:rg:kicks), with
+error O_L(rho^{j-1}).
The output is invariant under simultaneous spin and spatial
symmetries, reflection-positive at any seam, translation-covariant,
-and equal to $e^{v m'^2}times(microscopic)$ up to the
+and equal to e^{v m'^2}times(microscopic) up to the
corrections listed in the class definition.
-The rate $delta_j/delta_{j-1}ge H_j^{-1/20}$ is what makes the
-induction close: each step loses a fractional power of $H_j$ but
+The rate delta_j/delta_{j-1}ge H_j^{-1/20} is what makes the
+induction close: each step loses a fractional power of H_j but
accumulates only a finite number of steps in any finite total of scale
-under control. In particular, the final error cap at scale $0$
-satisfies $delta_0ge H_0^{-2.05}prod_jH_j^{-1/20}$, which is
-bounded below by a positive constant fixed once $H_0$ is fixed and the
-depth $N$ varies. The exact proof of Theorem~thm:rg-closure
+under control. In particular, the final error cap at scale 0
+satisfies delta_0ge H_0^{-2.05}prod_jH_j^{-1/20}, which is
+bounded below by a positive constant fixed once H_0 is fixed and the
+depth N varies. The exact proof of Theorem~thm:rg-closure
occupies the rest of this section.
## Proof structure
@@ -1737,7 +1731,7 @@ Use the free-energy Schur identity of Lemma~free:bounds to split
the exact density into the free kernel part plus a positive fluctuation
piece. Perform the observation exactly and record the block mean
residual (eq:rg:residual) together with its derivatives up to
-$k=8$.
+k=8.
noindentLayer 2: perturbation and the Gaussian-ratio identity.
Expand the exact negative-log density to quartic order in the
@@ -1753,13 +1747,13 @@ signed expansion is not the same as expanding the exact (positive)
density. But the exact density is positive, so its negative logarithm
is real. We retain the exact quadratic form (bounded by the free
kernel), and register each extra contribution as a regular error
-$f_X$ satisfying the cap (eq:rg:error-norm). The Taylor
+f_X satisfying the cap (eq:rg:error-norm). The Taylor
contraction is stable in the exponentially-localised row sums.
noindentLayer 4: covering corrections.
-All Taylor terms of degree $ge5$, and all the physical predicates
+All Taylor terms of degree ge5, and all the physical predicates
(bad-edge masks, cutoff support) that would otherwise be invisible in a
-signed expansion, are reorganised as a mandatory covering gas $Xi$
+signed expansion, are reorganised as a mandatory covering gas Xi
with bounded weight. The gas is a strict positive correction to the
density, so it can be summed by the covering argument of
Appendix~app:joint-majorants (in preparation, but the mechanism
@@ -1767,10 +1761,10 @@ is the same as a standard polymer-gas argument).
noindentLayer 5: regrouping and contraction.
Regroup the gas into connected components; the exponential weight of
-each component is paid by the exponential localisation of $ell_h$,
+each component is paid by the exponential localisation of ell_h,
the classical direct cost of a bad bond, and the exceptional masses
provided by the canonical coefficients. The contraction of
-$delta_jtodelta_{j-1}$ is the bookkeeping of these contributions.
+delta_jtodelta_{j-1} is the bookkeeping of these contributions.
The full analysis of each layer is in
Sections~sec:rg-proof--sec:rg-canonical; the analytic
@@ -1784,90 +1778,88 @@ the geometric bookkeeping in Appendix~app:rg-geometry.
In this section we supply the layer-by-layer proof of
Theorem~thm:rg-closure. The proof is organised as six lemmas. Three
-of them (Schur observation, Taylor contraction, $delta$-contraction)
+of them (Schur observation, Taylor contraction, delta-contraction)
constitute the bulk of the bookkeeping; the other three (gas regrouping,
symmetry transport, assembly) are lighter.
-Let us set up notation for the section. Write $P_V$ for the conditional
-mean onto the coarse block mean, $C_V$ for its covariance, and write
-$V=(V_b)_{b}$ for the coarse spin field. Constants $C_L,c_L>0$ depend
-only on the block size $L$ (integers $Lge L_0$) and on the microscopic
-parameters $a$, $R_*$; whenever we say that $L$ is fixed we mean these
+Let us set up notation for the section. Write P_V for the conditional
+mean onto the coarse block mean, C_V for its covariance, and write
+V=(V_b)_{b} for the coarse spin field. Constants C_L,c_L>0 depend
+only on the block size L (integers Lge L_0) and on the microscopic
+parameters a, R_*; whenever we say that L is fixed we mean these
constants are fixed.
## Lemma 1: The Schur observation split
[Schur observation]
-Let $rhoinmathscr C_j(b,t,scrC_j(b,t,delta,h)$ and let $Q(Vmid u)$ be the
-positive observation kernel (eq:obs-kernel). Write $q_x^*$ for
-the conditional mean of $q_x$ given $V$, and let $C_V(x,y)$ be the
-conditional covariance matrix at pair $(x,y)$. Then the exact
+Let rhoinmathscr C_j(b,t,scrC_j(b,t,delta,h) and let Q(Vmid u) be the
+positive observation kernel (eq:obs-kernel). Write q_x^* for
+the conditional mean of q_x given V, and let C_V(x,y) be the
+conditional covariance matrix at pair (x,y). Then the exact
integration over the fine spins decomposes as
int rho(q),Q(Vmid q),dd q
= e^{A_0(V)}cdot R(V),
where
-
-item $A_0(V)$ is a quadratic-plus-linear expression in $V$
+ A_0(V) is a quadratic-plus-linear expression in V
given by the free-energy Schur complement, whose quadratic
- form is $C_V^{-1}$ (bounded above by $c_L^{-1}$ and below by
- $C_L$);
-item $R(V)ge0$ is a remainder whose logarithm is a sum of
+ form is C_V^{-1} (bounded above by c_L^{-1} and below by
+ C_L);
+ R(V)ge0 is a remainder whose logarithm is a sum of
localised Taylor remainder clauses, each sitting inside a
- complete support $N_{14M}$ of an existing slot of the class.
+ complete support N_{14M} of an existing slot of the class.
Proof.
The exact identity (eq:free:schur-energy) writes
-$tfrac12inner u{K_hu}+tfrac a2|V-Qu|^2$ as a free energy in $V$
+tfrac12inner u{K_hu}+tfrac a2|V-Qu|^2 as a free energy in V
plus the positive fluctuation cost
-$tfrac12inner{u-P_hV}{C_h^{-1}(u-P_hV)}$. Writing $u_x$ for the
-fine spin inside block $b$ and $P_hV$ for its block-mean proxy, the
-difference $u_x-(P_hV)_x$ is directly controlled by the residual
+tfrac12inner{u-P_hV}{C_h^{-1}(u-P_hV)}. Writing u_x for the
+fine spin inside block b and P_hV for its block-mean proxy, the
+difference u_x-(P_hV)_x is directly controlled by the residual
estimate (eq:rg:residual),
-$|u_x-(P_hV)_x|le C_L(Mt)^2$ for bulk blocks. Consequently the
+|u_x-(P_hV)_x|le C_L(Mt)^2 for bulk blocks. Consequently the
quadratic-plus-linear piece is the Schur complement of the free energy
-on $V$, and the remainder is the exponential of the negative
+on V, and the remainder is the exponential of the negative
fluctuation cost plus the observation weight. Completing the square
-in $V$ gives the stated quadratic coefficient as exactly $C_h^{-1}$
+in V gives the stated quadratic coefficient as exactly C_h^{-1}
restricted to the block; the coarse-scale precision after observation
is bounded above and below uniformly because
-$C_h^{-1}=K_h+aQ^*Q$ and $K_h$ is bounded by the block Fourier analysis
+C_h^{-1}=K_h+aQ^*Q and K_h is bounded by the block Fourier analysis
of Lemma~free:bounds.
Every Taylor-remainder clause arising from this split is supported on
a full neighbourhood of an existing slot because (i) the observation
-kernel enters each block only through $barbar u$,u$, and (ii) the free-energy
+kernel enters each block only through barbar u,u, and (ii) the free-energy
identity only brings in the free-kernel rows, which are exponentially
localised. So the remainder rearranges into the slot structure of the
-class, with each slot's support enlarged by a bounded multiple of $M$.
+class, with each slot's support enlarged by a bounded multiple of M.
## Lemma 2: Taylor contraction to quartic order
[Taylor contraction]
-Let $F(q)$ be the negative-log of the class-$j$ density
+Let F(q) be the negative-log of the class-j density
(eq:rg:density) expanded in the fluctuation field
-$xi_x=q_xbarbar u_x^u_x^{-1}$ around the block centre $barbar u_x$u_x$ of the
+xi_x=q_xbarbar u_x^u_x^{-1} around the block centre barbar u_xu_x of the
coarse cell. Write
-$F(xi) = F_0 + F_1cdotxi + frac12xi^TF_2xi
+F(xi) = F_0 + F_1cdotxi + frac12xi^TF_2xi
+frac16 F_3(xi,xi,xi) +frac1{24}F_4(xi,xi,xi,xi)
-+ mathcal R_5(calR_5(xi)$
++ mathcal R_5(calR_5(xi)
for the Taylor expansion to quartic order at zero with remainder
-$mathcal R_5$.calR_5$. Under the small-set residual bound
-(eq:rg:residual) for some sufficiently large $M$, the coefficients
+mathcal R_5.calR_5. Under the small-set residual bound
+(eq:rg:residual) for some sufficiently large M, the coefficients
satisfy:
-
-item $F_1$ is $O(C_LMt)$, and its derivatives to order $k=8$ are
- $O(C_LM^Dt)$;
-item $F_2$ has a lower bound $c_L > 0$ (uniform in $L$) inside the
+ F_1 is O(C_LMt), and its derivatives to order k=8 are
+ O(C_LM^Dt);
+ F_2 has a lower bound c_L > 0 (uniform in L) inside the
chart;
-item $F_3$ has a bound $O(C_LM^D t)$;
-item $F_4$ is $O(C_LM^Dt^0)$ (bounded independently of $t$ and $L$);
-item $mathcal R_5$calR_5$ has the exponential-small bound
- $|mathcal R_5|calR_5|le C_Lg^5poly(M,p)$ in the
+ F_3 has a bound O(C_LM^D t);
+ F_4 is O(C_LM^Dt^0) (bounded independently of t and L);
+ mathcal R_5calR_5 has the exponential-small bound
+ |mathcal R_5|calR_5|le C_Lg^5poly(M,p) in the
free-kernel exponential norm.
@@ -1876,81 +1868,80 @@ Proof.
The coefficients are read off by the classical computation of the local
chart expansion around the aligned configuration; this is the standard
block-type Taylor expansion used in the perfect-action literature
-[cite]. Concretely: identify $S^3$
+[cite]. Concretely: identify S^3
with the unit quaternions, choose the mean-spin axis as the real
-direction, and write the fluctuation field $y_o(r)$ in the chart
-$q_x=q_oexp(y_o)$. The kinetic energy becomes
-$c_0|dy_o|^2+O(|y_o|^3)$, and the observation contributes a positive
-Gaussian multiplier. The bounds on $F_1,F_2,F_3,F_4$ follow from direct
-computation of the chart Jacobian and Hessian: $F_2ge c_L$ because the
+direction, and write the fluctuation field y_o(r) in the chart
+q_x=q_oexp(y_o). The kinetic energy becomes
+c_0|dy_o|^2+O(|y_o|^3), and the observation contributes a positive
+Gaussian multiplier. The bounds on F_1,F_2,F_3,F_4 follow from direct
+computation of the chart Jacobian and Hessian: F_2ge c_L because the
kinetic-energy Hessian plus the Schmidt-rank bound on the block
-observation provides it; $F_1,F_3$ are controlled by the residual
-(eq:rg:residual); $F_4$ is a polynomial coefficient of bounded
+observation provides it; F_1,F_3 are controlled by the residual
+(eq:rg:residual); F_4 is a polynomial coefficient of bounded
degree.
-The remainder $mathcal R_5$calR_5$ is bounded by direct use of the fourth-order
+The remainder mathcal R_5calR_5 is bounded by direct use of the fourth-order
remainder formula
-$mathcal R_5(calR_5(xi)=int_0^1((1-s)^4)/(4!)F^{(5)}(sxi)(xi^{ ⊗ 5})ds$,
-with the free-kernel exponential-sum bound controlling $F^{(5)}$.
-Since $F_2$ is roughly the free kinetic operator with Hessian $c_0$,
+mathcal R_5(calR_5(xi)=int_0^1((1-s)^4)/(4!)F^{(5)}(sxi)(xi^{ ⊗ 5})ds,
+with the free-kernel exponential-sum bound controlling F^{(5)}.
+Since F_2 is roughly the free kinetic operator with Hessian c_0,
the free-kernel exponential-sum bound for its inverse (which is our
-exponentially-localised covariance $C_h$) plus the locality convention
-of Appendix~app:rg-geometry gives the stated polynomial-in-$M$
+exponentially-localised covariance C_h) plus the locality convention
+of Appendix~app:rg-geometry gives the stated polynomial-in-M
bound.
-## Lemma 3: The $delta$-contraction
+## Lemma 3: The delta-contraction
[Error contraction]
-Let $f_X$ be a regular error slot of the class-$j$ density with cap
-$delta$ (in the norm (eq:rg:error-norm)). Under one exact RG step,
+Let f_X be a regular error slot of the class-j density with cap
+delta (in the norm (eq:rg:error-norm)). Under one exact RG step,
its image through the Taylor map plus the Schur split satisfies
|mathcal TcalT f_X|_{k,j-1,A}
le (C/L+o(1)),delta,.
-More generally, if $A_j$ denotes the error coefficient in the class at
-scale $j$, then
+More generally, if A_j denotes the error coefficient in the class at
+scale j, then
A_{j-1}le H_j^{-1/20}A_j,
-with $H_j$ the scale-$j$ coupling.
+with H_j the scale-j coupling.
Proof.
-The composition operator $mathcal T$calT$ is bounded on the norm
-$|cdot|_{k,j,A}$ by a constant $q_c<1$ provided $L$ is chosen large
+The composition operator mathcal TcalT is bounded on the norm
+|cdot|_{k,j,A} by a constant q_c<1 provided L is chosen large
enough. The bound comes from two sources:
-
-item the free-kernel row sum bound
- (eq:free:derivative) contributes a factor $C/L$ because
- there are $L^2$ old anchors per new anchor and each contributes
- around $L^{-2}$ mass;
-item the residual (eq:rg:residual) contributes a factor
- $(Mt_j)^2$, which is $o(1)$ as $H$ grows because
- $t_j=H_j^{-2}/p_H$ (that is, $t_j$ is polynomially small).
+ the free-kernel row sum bound
+ (eq:free:derivative) contributes a factor C/L because
+ there are L^2 old anchors per new anchor and each contributes
+ around L^{-2} mass;
+ the residual (eq:rg:residual) contributes a factor
+ (Mt_j)^2, which is o(1) as H grows because
+ t_j=H_j^{-2}/p_H (that is, t_j is polynomially small).
The same argument applied to the difference of two trajectories gives the
-stated contraction rate up to the residual $H_j^{-1/20}$ (this is the
-worst-case combination of the $C/L$ kernel contribution with polynomial
-$Mt_j$ residuals; see Appendix~app:finite-choices).
+stated contraction rate up to the residual H_j^{-1/20} (this is the
+worst-case combination of the C/L kernel contribution with polynomial
+Mt_j residuals; see Appendix~app:finite-choices).
## Lemma 4: Regrouping the covering gas
[Covering regrouping]
-Under one exact RG step, the mandatory covering gas $Xi_j$ produced by
+Under one exact RG step, the mandatory covering gas Xi_j produced by
the Taylor expansion plus the physical masks corresponds bijectively to
-a mandatory covering gas $Xi_{j-1}$ of the same form as
+a mandatory covering gas Xi_{j-1} of the same form as
(eq:rg:gas), with the anchored weight bound preserved up to a
-factor $H_j^{-1/20}$.
+factor H_j^{-1/20}.
Proof.
Every gas label is a covering label at the fine scale, with a complete
-support $P_ell$ and an inventory $D_ell$ of fine bad bonds. A fine
-bad bond at scale $j$ corresponds to a fine coarse-cell bad bond after
-folding once: the coarse bond $(O,O+e_mu)$ is bad precisely when at
+support P_ell and an inventory D_ell of fine bad bonds. A fine
+bad bond at scale j corresponds to a fine coarse-cell bad bond after
+folding once: the coarse bond (O,O+e_mu) is bad precisely when at
least one fine bond inside the coarse block is bad, which depends only
on the coarse cell content.
@@ -1960,15 +1951,15 @@ inventory that coarse bad bond and weight the sum of the products, which
is bounded by the anchored norm of the fine gas. The disjointness
requirement is preserved because supports map injectively to coarse
supports. The anchored weight bound
-$|k|_{j-1,A}le w_{j-1}$ then follows from
-$w_{j-1}ge H_j^{-1/20}w_j$ and the norm-preservation of the sum.
+|k|_{j-1,A}le w_{j-1} then follows from
+w_{j-1}ge H_j^{-1/20}w_j and the norm-preservation of the sum.
## Lemma 5: Symmetry and positivity transport
[Symmetry and reflection positivity]
-The output density of any class-$j$ density under one exact RG step is
+The output density of any class-j density under one exact RG step is
invariant under simultaneous spin and spatial symmetries, and is
reflection-positive at any seam.
@@ -1977,14 +1968,14 @@ Proof.
Each input slot is by the class definition invariant under simultaneous
spin/spatial symmetries, and the exact integration commutes with these
symmetries because the observation kernel
-(eq:obs-kernel) is $O(4)$-invariant and the block geometry is
+(eq:obs-kernel) is O(4)-invariant and the block geometry is
symmetric. Thus the output is invariant.
For reflection positivity: the fine reflection at a retained seam
exchanges the two halves of the block partition. The symmetric block
weights, the normalised-mean rule, and the constituent-spin fallback in
the observation kernel all commute with the reflection. For a function
-$F$ of output spins on one side, let $mathcal QF$calQF$ be its conditional
+F of output spins on one side, let mathcal QFcalQF be its conditional
expectation given the input spins on that side. Conditional independence
and reflection covariance give
[
@@ -2000,12 +1991,12 @@ covariance; that is the whole content.
## Lemma 6: Assembly
[Assembly]
-Let $rhoinmathscr C_j(b,t,scrC_j(b,t,delta,h)$. Then under one exact RG step,
-$rho$ is mapped to $rho'inmathscr C_scrC_{j-1}(b',t',delta',h')$ with
+Let rhoinmathscr C_j(b,t,scrC_j(b,t,delta,h). Then under one exact RG step,
+rho is mapped to rho'inmathscr C_scrC_{j-1}(b',t',delta',h') with
- b'-b = alpha+kappa/b+...,qquad delta' le H_j^{-1/20}delta,
+ b'-b = alpha+kappa/b+..., delta' le H_j^{-1/20}delta,
-where the dots absorb $O_L(b^{-2})$ errors.
+where the dots absorb O_L(b^{-2}) errors.
Proof.
@@ -2014,22 +2005,22 @@ Lemma~lem:rg-Taylor to expand to quartic order, then
Lemma~lem:rg-error-contraction and
Lemma~lem:rg-cover-regroup to control the higher-order remainder
and the physical masks, then Lemma~lem:rg-symmetry to transport
-symmetry and reflection positivity. The new coupling $b'$ is read from
+symmetry and reflection positivity. The new coupling b' is read from
the coefficient of the quadratic kinetic term in the output:
-$b'=b+F_2^{(2)}/mathcal E_calE_{h,t}^{(2)}$, where $F_2^{(2)}$ is the
+b'=b+F_2^{(2)}/mathcal E_calE_{h,t}^{(2)}, where F_2^{(2)} is the
quadratic part of the exact negative logarithm and
-$mathcal E_calE_{h,t}^{(2)}$ is the quadratic kinetic term. Matching
-$F_2^{(2)}$ against the classical two-loop computation gives the
+mathcal E_calE_{h,t}^{(2)} is the quadratic kinetic term. Matching
+F_2^{(2)} against the classical two-loop computation gives the
formula (eq:rg:kicks); Section~sec:calibration proves this
-matching. The new cutoff $t'$ and the new error cap $delta'$ are
+matching. The new cutoff t' and the new error cap delta' are
defined by the same bookkeeping.
Proof. [Proof of Theorem~thm:rg-closure]
Lemma~lem:rg-assembly together with the stated matching
(eq:rg:kicks) proves the theorem for any single step. Composing
-$n$ steps and taking $rho_0$ to be the terminal scale-$0$ density
-produces a class-$0$ density with the accumulated coupling and error
+n steps and taking rho_0 to be the terminal scale-0 density
+produces a class-0 density with the accumulated coupling and error
cap. The statement that the output is equal to the microscopic model
up to corrections in the class is contained in the construction
(Lemma~lem:rg-Schur) and verified in the terminal comparison
@@ -2042,69 +2033,69 @@ up to corrections in the class is contained in the construction
%---------------------------------------------------------------------
We now prove the two-loop matching (eq:rg:kicks) for the classical
-model, and then integrate the running coupling over $N$ RG steps to
+model, and then integrate the running coupling over N RG steps to
obtain the asymptotically free scale
-$L^Napproxsqrtbeta,e^{pibeta}$.
+L^Napproxsqrtbeta,e^{pibeta}.
## Preliminary: the spin chart and Laplace expansion
-Identify $S^3$ with the unit quaternions. For a spin configuration
-$q_xin S^3$ near an aligned configuration, we choose the mean-spin
+Identify S^3 with the unit quaternions. For a spin configuration
+q_xin S^3 near an aligned configuration, we choose the mean-spin
axis and write
q_x = (√(1-|u_x|^2),u_x)in S^3subsetRtimesR^3,
- qquad sum_xu_x=0,qquad z_x=sqrtbeta,u_x.
+ sum_xu_x=0, z_x=sqrtbeta,u_x.
-Under the nearest-neighbor measure $ddmu_beta$ the sum
-$sum_e q_ycdot q_x$ becomes
-$sum_e[1-tfrac12|u_y-u_x|^2+O(|u|^3)]$, so the quadratic fluctuation
-is a massive Gaussian of mass $sim$ nearest-neighbor Laplacian, and the
-cubic and higher terms are suppressed by $1/sqrtbeta$ per spin.
+Under the nearest-neighbor measure ddmu_beta the sum
+sum_e q_ycdot q_x becomes
+sum_e[1-tfrac12|u_y-u_x|^2+O(|u|^3)], so the quadratic fluctuation
+is a massive Gaussian of mass sim nearest-neighbor Laplacian, and the
+cubic and higher terms are suppressed by 1/sqrtbeta per spin.
-For a finite torus $Λ$ of side $n$, the Laplace expansion gives
+For a finite torus Λ of side n, the Laplace expansion gives
log Z_beta(n,n)=beta|Λ|-frac32logdet K_n+O(1/beta),
-where $K_n$ is the free Laplacian on $Λ$ (three tangent
-components $timesfrac12$) and $O(1/beta)$ is the quartic coefficient.
+where K_n is the free Laplacian on Λ (three tangent
+components timesfrac12) and O(1/beta) is the quartic coefficient.
Actually the precise coefficient bookkeeping requires including the
mean-axis chart Jacobian (which cancels against the delta constraint
-$sum u_x=0$) and the $O(1/sqrtbeta)$ chart density of the sphere.
+sum u_x=0) and the O(1/sqrtbeta) chart density of the sphere.
We take this as given; the calculation is standard (see
[cite]).
For our purposes, the two quantities we need are the coefficients
-$A(n)$ of $beta^{-1}$ and $Z_2(n)$ of $beta^{-2}$ in
-$log Z_beta(n,n)$ at fixed $n$, after the constant $n^2beta$ is
+A(n) of beta^{-1} and Z_2(n) of beta^{-2} in
+log Z_beta(n,n) at fixed n, after the constant n^2beta is
discarded. The doubling operators
[
- F_d(n)=4F(n)-F(2n),qquad D_0=3D/2=6,
+ F_d(n)=4F(n)-F(2n), D_0=3D/2=6,
]
-isolate the RG-relevant parts: the $4$ against the doubled-$2n$
+isolate the RG-relevant parts: the 4 against the doubled-2n
subtracts the (known) free-energy-with-shift and leaves only the
RG-relevant coefficient from the bulk.
## The calibration calculation
-We compute the sequence of RG kicks $alpha_h,kappa_h$ by the finite
+We compute the sequence of RG kicks alpha_h,kappa_h by the finite
torus Laplace expansion. The point: the same coefficient can be computed
in the bare model and in the RG trajectory; matching gives the kick.
[Calibration]
-Let $alpha_h,kappa_h$ be the canonical kinetic kicks at depth $h$
+Let alpha_h,kappa_h be the canonical kinetic kicks at depth h
starting from zero canonical interactions, so that a single step has
form
[
beta'=beta+alpha_h+(kappa_h)/(beta)+Delta_h,
]
-with $Delta_h$ the many-body error supplied by the exact RG step.
+with Delta_h the many-body error supplied by the exact RG step.
Define
[
- B_s=sum_{h<s}alpha_h,qquad
+ B_s=sum_{h<s}alpha_h,
C_s=sum_{h<s}kappa_h.
]
-Then, for the nearest-neighbor $O(4)$ normalization of
-(eq:model-ov) and $gamma=(log L)/pi$,
+Then, for the nearest-neighbor O(4) normalization of
+(eq:model-ov) and gamma=(log L)/pi,
B_s&=-gamma s+O_L(1),&
C_s&=-(gamma)/(2pi)s+O_L(1),
@@ -2112,15 +2103,15 @@ Then, for the nearest-neighbor $O(4)$ normalization of
kappa_h&=-(gamma)/(2pi)+O_L(rho^h)
-for some $rho<1$ depending only on $L$. The constants are independent
+for some rho<1 depending only on L. The constants are independent
of the number of subsequent block transformations.
Proof.
-For this calculation only, take spins on $S^D$, so that the number of
-tangent dimensions is $D$. Let $v=n^2$ (side $n$ on the torus). Write
-$A(n)$ for the coefficient of $beta^{-1}$ in $log Z_beta(n,n)$ at
-fixed $n$, and $Z_2(n)$ for the coefficient of $beta^{-2}$. A
+For this calculation only, take spins on S^D, so that the number of
+tangent dimensions is D. Let v=n^2 (side n on the torus). Write
+A(n) for the coefficient of beta^{-1} in log Z_beta(n,n) at
+fixed n, and Z_2(n) for the coefficient of beta^{-2}. A
Laplace expansion around the aligned configuration on the chart
(eq:cal:chart) gives, for the mean-axis integration,
@@ -2128,91 +2119,91 @@ Laplace expansion around the aligned configuration on the chart
(v^{-1}sum_x√(1-|u_x|^2))^D.
Indeed the delta constraint on the mean tangent component has axis
-Jacobian $|barbar q|^q|^{-D}$, so the factor cancels that Jacobian. The
+Jacobian |barbar q|^q|^{-D}, so the factor cancels that Jacobian. The
root with positive mean component is selected. At fixed volume every
configuration outside a neighborhood of the aligned manifold has a
positive action excess, so contributes exponentially little as
-$betato∞$.
+betato∞.
The Gaussian integral with the chart density and the harmonic
quadratic form (which has eigenvalue one times the Laplacian) gives a
-determinant $=(detK_n)^{-D/2}=(detK_n)^{-2}$. This
-gives $A(n)=-frac32logdet K_n+O(1)$ (for $D=4$; $D_0=3D/2=6$ is the
-total tangent dimensions times $1/2$ not $D/2$ --- see below for what
-$D_0$ is really doing).
-
-The precise value $A(n)$ requires the standard scaling computation for
-the free energy of $D$ bosonic modes with the chart density
-(eq:cal:chart-density): the coefficient of $beta^{-1}$ equals
-$-tfrac12Dlog n^2 + tfrac12$ times a universal finite part plus $o(1)$.
-We denote this by $A(n)=-tfrac32log n+C_A+o(1)$.
-
-Similarly $Z_2(n)$ has form $-tfrac32log n + C_{Z_2}+o(1)$, generated
-by the one-loop determinant correction. Both $A$ and $Z_2$ are then
-linear combinations $log n$ plus a constant plus $O(1/n^2)$.
-
-Applying the doubling $F_d(n)=4F(n)-F(2n)$, the $log n$ cancels:
-$F_d(n)=3F_{const}+O(1/n^2)$. Given the chart density (which
+determinant =(detK_n)^{-D/2}=(detK_n)^{-2}. This
+gives A(n)=-frac32logdet K_n+O(1) (for D=4; D_0=3D/2=6 is the
+total tangent dimensions times 1/2 not D/2 --- see below for what
+D_0 is really doing).
+
+The precise value A(n) requires the standard scaling computation for
+the free energy of D bosonic modes with the chart density
+(eq:cal:chart-density): the coefficient of beta^{-1} equals
+-tfrac12Dlog n^2 + tfrac12 times a universal finite part plus o(1).
+We denote this by A(n)=-tfrac32log n+C_A+o(1).
+
+Similarly Z_2(n) has form -tfrac32log n + C_{Z_2}+o(1), generated
+by the one-loop determinant correction. Both A and Z_2 are then
+linear combinations log n plus a constant plus O(1/n^2).
+
+Applying the doubling F_d(n)=4F(n)-F(2n), the log n cancels:
+F_d(n)=3F_{const}+O(1/n^2). Given the chart density (which
must be integrated against a small-mean constraint) the explicit
coefficient is
[
- alpha_h = -gamma+O_L(rho^h),quad
+ alpha_h = -gamma+O_L(rho^h),
kappa_h = -(gamma)/(2pi)+O_L(rho^h),
]
with the exponential corrections from the finite-size eigenvalues of
the Laplacian at the boundary of the small-(aligned) set.
-The key observation is: $alpha$ and $kappa$ are the RG-relevant
+The key observation is: alpha and kappa are the RG-relevant
parts of the bare Laplace expansion under doubling. Because the
-doubling operation cancels the volume-scaling terms with $log n$
+doubling operation cancels the volume-scaling terms with log n
exactly, the RG kicks depend only on the finite-size corrections, and
-those are precisely $-gamma$ and $-gamma/(2pi)$ with error
-$O_L(rho^h)$ since the finite-size correction at depth $h$ is
-exponentially small in $h$ for the modes integrated out.
+those are precisely -gamma and -gamma/(2pi) with error
+O_L(rho^h) since the finite-size correction at depth h is
+exponentially small in h for the modes integrated out.
-Integrating: $sum_{h<s}alpha_h=-gamma s+O(1)$,
-$sum_{h<s}kappa_h=-(gamma)/(2pi)s+O(1)$. This proves
+Integrating: sum_{h<s}alpha_h=-gamma s+O(1),
+sum_{h<s}kappa_h=-(gamma)/(2pi)s+O(1). This proves
(cal:sums) and~(cal:kicks-thm).
## From running coupling to scale
-Apply the calibration to the RG trajectory: at depth $r$ the coupling
+Apply the calibration to the RG trajectory: at depth r the coupling
satisfies
beta_r = beta_N - gamma(N-r)
-(gamma)/(2pi)sum_{h=r}^{N-1}(1)/(beta_h)
+ O(1)
-where $beta_h$ is the coupling at chronological depth $h$ (so
-$beta_0$ is the terminal and $beta_N$ the bare coupling).
+where beta_h is the coupling at chronological depth h (so
+beta_0 is the terminal and beta_N the bare coupling).
-For sufficiently large $beta$, the sum on the right-hand side is
+For sufficiently large beta, the sum on the right-hand side is
approximately
-$(N-r)/(beta_{r+N})$, comparable to two-loop running. Explicitly,
-setting $beta_0=H$ at the terminal, we solve for the required bare
+(N-r)/(beta_{r+N}), comparable to two-loop running. Explicitly,
+setting beta_0=H at the terminal, we solve for the required bare
coupling:
beta_N=H+gamma N+(gamma)/(2pi)sum_{h=0}^{N-1}(1)/(beta_h)+O(1).
-Choose $N$ such that $beta_N=beta$, the desired bare inverse coupling.
+Choose N such that beta_N=beta, the desired bare inverse coupling.
Then to leading order,
N=(beta-H)/(gamma)-(1)/(2pigamma)sum_{h=0}^{N-1}(1)/(beta_h)+O(1)
=(beta)/(gamma)-frac12(logbeta)/(gamma)cdot(1)/(pi)+O(1),
-using the Euler-Maclaurin asymptotics $sum_{h=0}^{N-1}1/beta_h
+using the Euler-Maclaurin asymptotics sum_{h=0}^{N-1}1/beta_h
sim(gamma)/(gamma^2)int ds/beta(s)
-simfracgamma{2pi}(something)$ for the two-loop running.
+simfracgamma{2pi}(something) for the two-loop running.
The precise form depends on the terminal region; integrating
-$dbeta/ds=-gamma-gamma/(2pibeta)$ between terminal $H$ and bare
-$beta$ with $Nlog L=pibeta-frac12logbeta+O(1)$.
+dbeta/ds=-gamma-gamma/(2pibeta) between terminal H and bare
+beta with Nlog L=pibeta-frac12logbeta+O(1).
-Equivalently: the total RG scale $L^N$ satisfies
+Equivalently: the total RG scale L^N satisfies
L^Nasympsqrtbeta,e^{pibeta},e^{O_{L,H}(1)}.
This is the asymptotically free scale formula. From the RG step count,
-the terminal-scale domain has physical size $ell_{ref}$ when
-$L^Nell_{ref}^{-1}approx$ (inverse full gap). Using
-$gap(beta)sim L^{-N}$ gives the sharp order of $gap$.
+the terminal-scale domain has physical size ell_{ref} when
+L^Nell_{ref}^{-1}approx (inverse full gap). Using
+gap(beta)sim L^{-N} gives the sharp order of gap.