Heat-kernel spectral bound
ProvedMarkovMixing.heat_kernel_spectral_boundLet be an irreducible Markov chain on a finite state space with stationary distribution , reversible with respect to it (detailed balance: ). The heat kernel at real time is — the continuous-time chain driven by a rate-one Poisson clock. Among the eigenvalues of (real with for some nonzero ), is the largest eigenvalue different from , and is the spectral gap (Mission VII).
The theorem (Theorem 20.6 of Levin–Peres–Wilmer) asserts the pointwise estimate: for all states and all ,
In continuous time the convergence rate is governed by the spectral gap itself — not the absolute gap of the discrete theory: negative eigenvalues, which cause the discrete chain to oscillate, are killed by the Poisson smoothing, since 's eigenvalue becomes for the heat kernel. The prefactor comes from expanding the transition kernel in the -orthonormal eigenbasis of Mission VII's spectral representation. This estimate is the continuous-time engine behind the product-chain theorem of this mission.
import Definitions.Def_mm_continuous import Mathlib.Analysis.SpecialFunctions.Sqrt
namespace MarkovMixing
/-- **Theorem 20.6** (LPW): for a reversible irreducible chain with spectral
gap `γ`, the heat kernel satisfies
`|H_t(x,y) − π(y)| ≤ √(π(y)/π(x)) e^{-γt}`. -/
theorem heat_kernel_spectral_bound {V : Type*} [Fintype V] [DecidableEq V]
[Nonempty V] (P : Matrix V V ℝ) (hP : IsStochastic P)
(hirr : Irreducible P) (π : V → ℝ) (hπ : IsStationary P π)
(hrev : DetailedBalance P π) (x y : V) (t : ℝ) (ht : 0 ≤ t) :
|heatKernel P t x y - π y| ≤
Real.sqrt (π y / π x) * Real.exp (-(spectralGap P) * t) := by
sorry
end MarkovMixing