Subelliptic Random Walks on Riemannian Manifolds and Their Convergence to Equilibrium
Davide Tramontana
TL;DR
This work constructs a global $(h,\rho)$-subelliptic diffusion on a closed Riemannian manifold driven by a subelliptic operator $A$, and proves convergence to a stationary distribution using a detailed spectral analysis of the associated Markov operator $S_{h,\rho}$. The authors deploy Fefferman–Phong localization to reduce $A$ locally to a universal model $\tilde{A}$, then patch these local models into a global diffusion via a finite subunit atlas. They establish a lower bound for the Markov kernel, perform a high/low frequency decomposition to obtain Weyl-type spectral estimates, and derive a Weyl-type infinitesimal generator $\mathsf{A}_{0,\rho}$, leading to a discrete spectrum near eigenvalue 1 and a spectral gap $g_{h,\rho}$ of order $h^2$. Consequently, the diffusion converges to the equilibrium distribution at an exponential rate $e^{-k g_{h,\rho}}$ in total variation, with the rate improving as $h\to0$. The results connect subelliptic diffusion on manifolds to Gibbs–Metropolis-type dynamics, providing quantitative, geometry-aware convergence guarantees governed by the subellipticity order $\varepsilon$.
Abstract
The aim of this work is to study the convergence to equilibrium of an $(h,ρ)$-subelliptic random walk on a closed, connected Riemannian manifold $(M,g)$ associated with a subelliptic second-order differential operator $A$ on $M$. In such a random walk, $h$ roughly represents the step size and $ρ$ the speed at which it is carried out. To construct the random walk and prove the convergence result, we employ a technique due to Fefferman and Phong, which reduces the problem to the study of a constant-coefficient operator $\tilde{A}$ that is locally equivalent to our second-order subelliptic operator $A$, in the sense that the diffusion generated by $\tilde{A}$ induces a local diffusion for $A$. By using the compactness of $M$ this local diffusion can be lifted to a global diffusion, and the convergence result is then obtained via the spectral theory of the associated Markov operator.
