Documentation

LeanPool.Superdiffusion.Algsuperdiff.Frozen.Introduction.GeneratorRenormalization

GeneratorRenormalization #

Ported from scottnarmstrong/Superdiffusion at 7f20476385d4af418f110303085ad0740c9eba46. Module imports, support namespace and local option directives adapted for Lean Pool.

Generator renormalization for the coefficient field — [ABK] Theorem B #

Source: algsuperdiff.tex of arXiv:2601.22142v2, Theorem B of the introduction (t.homogenization).

At every scale m the coarse-grained generator of the coefficient field a = nu I + k is a constant multiple sigmaBar of the identity, quantitatively close to sqrt (nu^2 + cstar gamma^{-1} 3^{2 gamma m}), and the disorder solution u is compared with the homogenized solution v of the same Dirichlet problem: a single random error amplitude with Gaussian-type moments controls both the uniform difference u - v and the difference of the two Dirichlet energies.

The carrier and the coefficient field #

The disorder is read on Algsuperdiff.Section5.Field.FullSample d M.gamma under Algsuperdiff.Section5.Field.fullSampleLaw M, the carrier the quenched superdiffusivity estimates use: the subtype of samples on which the field's two-sided scale rates and gradient bounds hold. That event has full measure for the unrestricted law, so the restriction loses nothing, and it is what makes the coefficient field a total function of the sample rather than a partially defined one.

The field is Algsuperdiff.Section5.Field.streamCoefficient M.nu omega, that is nu I + k_omega with k_omega x = Σ_{n ∈ ℤ} (j_n x - j_n 0). The shell sum is normalized at the origin, which fixes the additive constant the source leaves free. The choice is invisible to the two Dirichlet problems of the statement: the shift is by a constant antisymmetric matrix, and Algsuperdiff.Section4.Provider.ExcessDecay.isDirichletSolutionOn_sub_const_of_skew says that such a shift moves no solution of the weak problem.

The form with an infrared cutoff #

Algsuperdiff.Frozen.Section4.generator_renormalization asserts the same estimate for the truncated field (coefficientCutoff nu L omega).toCoeffField, uniformly in the cutoff L ≥ m, on the unrestricted sample space. That form is the more general one and stays part of the public surface; the statement here is the one the introduction prints, and follows from it in the limit L → ∞. Outside the carrier, the law and the coefficient field the two statements share every token.

Reading of the statement #

Proved; reduces to the standard axioms.

theorem Algsuperdiff.Frozen.Introduction.generator_renormalization (d : ℕ) (cstar : ℝ) (_hcstar : 0 < cstar) :
∃ (gamma0 : ℝ) (C : ℝ), 0 < gamma0 ∧ 0 < C ∧ ∀ (M : Section3.ABKModel d), Section3.Disorder.cstar M = cstar → M.gamma ≤ gamma0 → ∀ (m : ℤ), ∃ (sigmaBarM : ℝ), 0 < sigmaBarM ∧ |sigmaBarM - √(M.nu ^ 2 + cstar * M.gamma⁻¹ * Real.rpow 3 (2 * M.gamma * ↑m))| ≤ C * √M.gamma * |Real.log M.gamma| * sigmaBarM ∧ ∃ (EB : Section5.Field.FullSample d M.gamma → ℝ), (∀ (omega : Section5.Field.FullSample d M.gamma), 0 ≤ EB omega) ∧ Measurable EB ∧ (∀ (p : ℝ), 1 ≤ p → p ≤ C⁻¹ * M.gamma⁻¹ * |Real.log M.gamma|⁻¹ → ∫⁻ (omega : Section5.Field.FullSample d M.gamma), ENNReal.ofReal (EB omega) ^ p ∂↑(Section5.Field.fullSampleLaw M) ≤ ENNReal.ofReal (C * (√p + √|Real.log M.gamma|) * √M.gamma * |Real.log M.gamma| ^ 3) ^ p) ∧ ∀ᵐ (omega : Section5.Field.FullSample d M.gamma) ∂↑(Section5.Field.fullSampleLaw M), ∀ (u v h : HCPolySupport.H1Function (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) (g : HCPolySupport.Vec d → HCPolySupport.Vec d) (Kg Kh KhInf : ℝ), Section4.Support.IsDirichletSolutionOn (Section5.Field.streamCoefficient M.nu omega) (HCPolySupport.originCube d m) u h g → Section4.Support.IsDirichletSolutionOn (fun (x : HCPolySupport.Vec d) => sigmaBarM • 1) (HCPolySupport.originCube d m) v h g → Section4.Support.HolderSeminormBoundOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) (1 / 2) Kg g → Section4.Support.HolderSeminormBoundOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) (1 / 2) Kh h.grad → (∀ x ∈ HCPolySupport.openCubeSet (HCPolySupport.originCube d m), ‖h.grad x‖ ≤ KhInf) → Section4.Support.HasGradientOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) h.toFun h.grad → (∀ᵐ (x : HCPolySupport.Vec d) ∂MeasureTheory.volume.restrict (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)), Real.rpow 3 (-↑m) * |u.toFun x - v.toFun x| ≤ EB omega * (sigmaBarM⁻¹ * Real.rpow 3 (↑m / 2) * Kg + (KhInf + Real.rpow 3 (↑m / 2) * Kh))) ∧ |(HCPolySupport.volumeAverage (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) fun (y : HCPolySupport.Vec d) => M.nu * HCPolySupport.vecNormSq (u.grad y)) - HCPolySupport.volumeAverage (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) fun (y : HCPolySupport.Vec d) => sigmaBarM * HCPolySupport.vecNormSq (v.grad y)| ≤ EB omega * (√sigmaBarM⁻¹ * Real.rpow 3 (↑m / 2) * Kg + √sigmaBarM * (KhInf + Real.rpow 3 (↑m / 2) * Kh)) ^ 2