Documentation

LeanPool.Superdiffusion.Algsuperdiff.MainTheorems

MainTheorems #

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

Main results #

The three headline theorems of the formalization, stated here in full: Theorems A, B and C of the paper Superdiffusion and anomalous regularization in self-similar random incompressible flows (Armstrong--Bou-Rabee--Kuusi), in their introduction-level forms. Each theorem below restates its certified counterpart verbatim and is proved by direct application, so the statements displayed in this file are byte-faithful to the certified ones.

All three reduce to the standard axioms (propext, Classical.choice, Quot.sound).

The coefficient field, and the statements with an infrared cutoff #

Theorems B and C are stated here for the coefficient field a = ν I + k of the model itself, read on the full sample of the disorder, in their introduction-level forms. The paper proves them first for the infrared truncations a_L of that field, uniformly in the cutoff L ≥ m, and those uniform statements are formalized as well and remain part of the public surface, as Algsuperdiff.Frozen.Section4.generator_renormalization and Algsuperdiff.Frozen.Section4.anomalous_regularity. They are the more general form; the two theorems below follow from them in the limit L → ∞, and outside the law and the coefficient field the two forms share every token.

The asymptotics of the effective diffusivity #

The two-sided asymptotics of the effective diffusivity -- at every scale it tracks the superdiffusive profile (ν² + c⋆ γ⁻¹ 3^{2γm})^{1/2} predicted by the renormalization-group heuristic, with relative error of order √γ |log γ| -- is the first display of Theorem B below, the comparison between the effective diffusivity sigmaBar_m and that profile. The quantitative form proved in Section 3, which adds the two-sided bounds on sigmaBar_m at the top scale of the induction, stays public as Algsuperdiff.Frozen.Section3.diffusivity_asymptotics.

Non-vacuity #

The sample space carrying each statement is inhabited, its law being a probability measure. The model clause is another matter: this development constructs no ABKModel d, so each theorem is certified as an implication, holding for every model that satisfies the paper's assumptions. The paper's Gaussian example is the construction that would exhibit such a model, and it is not formalized here. The shell carrier also requires qualitative C² regularity, stronger than the paper's local C^{1,1} assumption. These statements preserve the paper's quantitative conclusions under that carrier regularity assumption.

Theorem A (quenched superdiffusivity). For every model whose coupling γ lies below a threshold, every time t > 0 and every exponent p in the range p ≤ C⁻¹ γ⁻¹ |log γ|⁻⁶, the annealed p-th moments of the deviation of the quenched second moment of the displacement at time t from 2 d Rtilde(t)², and of the squared quenched mean displacement, are bounded by (C (√p + √|log γ|) √γ |log γ|⁴ Rtilde(t)²)^p and (C (p + |log γ|) γ |log γ|⁷ Rtilde(t)²)^p respectively, where Rtilde(t) is the intrinsic length scale of the paper (intrinsicScale). The displacement is read through the compactified process, on which the retraction onePointRetract sends the added point to the origin.

theorem Algsuperdiff.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

Theorem B (renormalization of the generator). For every model whose coupling γ lies below a threshold and every scale m there are an effective diffusivity sigmaBar_m > 0, whose relative distance to (ν² + c⋆ γ⁻¹ 3^{2γm})^{1/2} is at most C √γ |log γ|, and a nonnegative random variable E_B on the sample of the coefficient field a = ν I + k, whose p-th moments over the disorder are bounded by (C (√p + √|log γ|) √γ |log γ|³)^p for every p in [1, C⁻¹ γ⁻¹ |log γ|⁻¹], such that almost surely the following holds for every choice of data at the scale m. If u solves the Dirichlet problem for a on the cube □_m with force g and boundary datum h, and v solves the homogenized problem with the constant coefficient sigmaBar_m I and the same data, then the rescaled uniform difference 3^{-m} |u - v| is at most E_B times sigmaBar_m⁻¹ 3^{m/2} [g] + (‖∇h‖_∞ + 3^{m/2} [∇h]), and the two Dirichlet energies ⨍ ν |∇u|² and ⨍ sigmaBar_m |∇v|² differ by at most E_B times the square of √(sigmaBar_m⁻¹) 3^{m/2} [g] + √(sigmaBar_m) (‖∇h‖_∞ + 3^{m/2} [∇h]). Here [·] is the 1/2-Hölder seminorm over the cube, for which the statement carries the explicit bounds Kg, Kh and the sup-norm bound KhInf, and the uniform estimate is read as an almost-everywhere bound on the cube. A single exceptional set serves the whole family of data; the moment bound on E_B is the separate annealed clause.

theorem Algsuperdiff.anomalous_regularity (d : ℕ) (cstar : ℝ) (_hcstar : 0 < cstar) :
∃ (gamma0 : ℝ) (C : ℝ), 0 < gamma0 ∧ 0 < C ∧ ∀ (M : Section3.ABKModel d), Section3.Disorder.cstar M = cstar → M.gamma ≤ gamma0 → ∀ (alpha : ℝ), 0 < alpha → alpha ≤ 1 - C * √M.gamma → ∀ (m : ℤ), ∃ (X : Section5.Field.FullSample d M.gamma → ℕ∞), Measurable X ∧ (∀ (N : ℕ), ↑(Section5.Field.fullSampleLaw M) {omega : Section5.Field.FullSample d M.gamma | ↑N ≤ X omega} ≤ ENNReal.ofReal (C * Real.exp (-((1 - alpha) ^ 2 * (↑N - C)) / (C * M.gamma)))) ∧ ∀ᵐ (omega : Section5.Field.FullSample d M.gamma) ∂↑(Section5.Field.fullSampleLaw M), ∀ (u h : HCPolySupport.H1Function (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) (g : HCPolySupport.Vec d → HCPolySupport.Vec d) (Kg Kh : ℝ), Section4.Support.IsDirichletSolutionOn (Section5.Field.streamCoefficient M.nu omega) (HCPolySupport.originCube d m) u 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 → (∀ y ∈ HCPolySupport.openCubeSet (HCPolySupport.originCube d m), ‖h.grad y‖ ≤ Real.rpow 3 (↑m / 2) * Kh) → Section4.Support.HasGradientOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m)) h.toFun h.grad → ∀ x ∈ HCPolySupport.openCubeSet (HCPolySupport.originCube d m), ∀ n ≤ m, X omega ≤ ↑(m - n).toNat → ENNReal.ofReal √M.nu * MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => √(HCPolySupport.vecNormSq (u.grad y))) 2 (Section4.Support.normalizedVolumeMeasureOn ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) ≤ ENNReal.ofReal (C * Real.rpow 3 ((1 - alpha) * (↑m - ↑n))) * (ENNReal.ofReal √M.nu * MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => √(HCPolySupport.vecNormSq (u.grad y))) 2 (Section4.Support.normalizedVolumeMeasureOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) + ENNReal.ofReal (√(↑(Section3.Annealed.sigmaBar M m))⁻¹ * Real.rpow 3 (↑m / 2) * Kg) + ENNReal.ofReal (√↑(Section3.Annealed.sigmaBar M m) * Real.rpow 3 (↑m / 2) * Kh)) ∧ (x ∈ HCPolySupport.openCubeSet (HCPolySupport.originCube d (m - 1)) → ENNReal.ofReal √M.nu * MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => √(HCPolySupport.vecNormSq (u.grad y))) 2 (Section4.Support.normalizedVolumeMeasureOn ((fun (y : HCPolySupport.Vec d) => x + y) '' HCPolySupport.openCubeSet (HCPolySupport.originCube d n) ∩ HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) ≤ ENNReal.ofReal (C * Real.rpow 3 ((1 - alpha) * (↑m - ↑n))) * (ENNReal.ofReal √M.nu * MeasureTheory.eLpNorm (fun (y : HCPolySupport.Vec d) => √(HCPolySupport.vecNormSq (u.grad y))) 2 (Section4.Support.normalizedVolumeMeasureOn (HCPolySupport.openCubeSet (HCPolySupport.originCube d m))) + ENNReal.ofReal (√(↑(Section3.Annealed.sigmaBar M m))⁻¹ * Real.rpow 3 (↑m / 2) * Kg)))

Theorem C (anomalous regularity). For every model whose coupling γ lies below a threshold, every Hölder exponent α in (0, 1 - C √γ] and every scale m there is a random minimal scale X, with values in ℕ∞ and with the exponentially decaying tail P[X ≥ N] ≤ C exp(-(1 - α)² (N - C) / (C γ)), such that almost surely the following excess decay holds for every choice of data at the scale m. If u solves the Dirichlet problem for the coefficient field a = ν I + k on the cube □_m with force g and boundary datum h, then at every point x of the cube and every scale n ≤ m with X ≤ m - n, the normalized energy √ν ‖∇u‖_{Lnormalized²((x + □_n) ∩ □_m)} is at most C 3^{(1-α)(m-n)} times the sum of the same energy over the whole cube, of √(sigmaBar_m⁻¹) 3^{m/2} Kg and of √(sigmaBar_m) 3^{m/2} Kh, the last term being present only when x lies outside the smaller cube □_{m-1}; that indicator is read here as the conjunction of the estimate with the term, valid at every x, and the estimate without it, valid on □_{m-1}. Here Kg bounds the 1/2-Hölder seminorm of g, while Kh bounds both the 1/2-Hölder seminorm of h.grad and its scaled size: ‖h.grad y‖ ≤ 3^{m/2} Kh throughout the cube. The energy norms are taken against the normalized volume measure of the set, and sigmaBar_m is the annealed diffusivity that Theorem B produces, quoted here by name. The decay 3^{(1-α)(m-n)} is the one a C^{0,α} function has, uniformly in the molecular diffusivity ν.