math.PRAug 4, 2026

A Direct Route to Markov Chain Convergence via Asymptotic Equivalence with the Target

Authors: Patrick Forré

Organizations: AI4Science Lab Korteweg-de Vries Institute for Mathematics University of Amsterdam

Abstract

For a Markov kernel TT with an invariant probability measure ππ, we give a self-contained proof of the Markov chain convergence theorem via a criterion called asymptotic equivalence with the target. It assumes two parts about the Lebesgue decompositions of TxnT^{n}_{x} and ππ for every starting point xx: 1.) asymptotic absolute continuity: the singular mass sing(Txnπ)(T^{n}_{x}\midπ) tends to 00; 2.) asymptotic domination of the target: the singular mass sing(πTxn)(π\mid T^{n}_{x}) tends to 00, as nn \to \infty. This criterion, on countably generated measurable spaces, is both sufficient and necessary for the Markov chain convergence. A density version of this criterion is verified on general measurable spaces in three cases: (i) TT has a positive transition density wrt ππ; (ii) TT consists of an absolutely continuous part with positive transition density together with an atom at the starting point, which covers the Metropolis--Hastings algorithm; (iii) the transition density is positive only after a finite number of steps that may depend on the starting point xx. To demonstrate our general criterion, we investigate the Gibbs sampler with random scan and the parallel tempering algorithm. Furthermore, we show that in all mentioned settings Birkhoff's ergodic theorem applies, so as to obtain the strong law of large numbers. Throughout this paper, neither irreducibility, nor aperiodicity, nor recurrence, nor couplings, nor splitting constructions, nor small sets are used. In most results, the state space is a general measurable space, which carries no structure beyond a σσ-algebra. Countable generation is only assumed where the density-free form of the criterion is stated. None of the theorems proved here is new; what is offered is a short route to a single, widely applicable Markov chain convergence criterion, which is both sufficient and necessary.

Explore similar work

Jun 26, 2025stat.ML

Gaussian Invariant Markov Chain Monte Carlo

We develop sampling methods, which consist of Gaussian invariant versions of random walk Metropolis (RWM), Metropolis adjusted Langevin algorithm (MALA) and second order Hessian or Manifold MALA. Unlike standard RWM and MALA, we show that Gaussian invariant sampling can lead to ergodic estimators with improved statistical efficiency. This is due to a remarkable property of Gaussian invariance that allows us to obtain exact analytical solutions to the Poisson equation for Gaussian targets. These solutions can be used to construct efficient and easy to use control variates for variance reduction of estimators under any intractable target. We demonstrate the new samplers and estimators in several examples, including high dimensional targets in latent Gaussian models where we compare against several advanced methods and obtain state-of-the-art results. We also provide theoretical results regarding geometric ergodicity, and an optimal scaling analysis that shows the dependence of the optimal acceptance rate on the Gaussianity of the target.
Michalis K. Titsias, Angelos Alexopoulos, Siran Liu +1
May 28, 2026stat.CO

True Self-Avoiding Walk for Accelerating Markov-Chain Monte Carlo Integration

We study true self-avoiding walk (TSAW) as a mechanism for improving empirical integral estimation via Markov chain Monte Carlo (MCMC). We consider finite-state adaptive sampling dynamics associated with an irreducible Markov kernel PP on a finite set, with stationary distribution ππ, in which the transition probabilities are penalized according to empirical overuse. Our main result is that the empirical occupation counts Lt(i)L_t(i) and transition counts Nt(i,j)N_t(i,j) of the resulting TSAW-based walk satisfy Lt(i)tπi=O(logt)andNt(i,j)tπiPij=O(logt)almost surelyL_t(i)-tπ_i = O(\sqrt{\log t}) \quad\text{and}\quad N_t(i,j)-tπ_iP_{ij}=O(\sqrt{\log t}) \qquad\text{almost surely} for every state ii and every edge (i,j)(i,j) with Pij>0P_{ij}>0. Consequently, for every bounded function f:VRf:V\to\mathbb R, the error of our integral estimator converges as 1ts=0t1f(Xs)iVπif(i)=O(logtt)almost surely.\left|\frac1t\sum_{s=0}^{t-1} f(X_s)-\sum_{i\in V}π_i f(i)\right| = O\left(\frac{\sqrt{\log t}}{t}\right) \qquad\text{almost surely}. These results show that, in contrast with the usual t1/2t^{-1/2} error scaling for empirical averages under standard random-walk-based methods, TSAW-based estimator yields empirical integral errors of order O(logt/t)O(\sqrt{\log t}/t) almost surely, thereby achieving a substantially sharper dependence on the sample size tt.
Qinghua, Ding, Venkat Anantharam
May 1, 2026stat.ML

Decentralized Proximal Stochastic Gradient Langevin Dynamics

We propose Decentralized Proximal Stochastic Gradient Langevin Dynamics (DE-PSGLD), a decentralized Markov chain Monte Carlo (MCMC) algorithm for sampling from a log-concave probability distribution constrained to a convex domain. Constraints are enforced through a shared proximal regularization based on the Moreau-Yosida envelope, enabling unconstrained updates while preserving consistency with the target constrained posterior. We establish non-asymptotic convergence guarantees in the 2-Wasserstein distance for both individual agent iterates and their network averages. Our analysis shows that DE-PSGLD converges to a regularized Gibbs distribution and quantifies the bias introduced by the proximal approximation. We evaluate DE-PSGLD for different sampling problems on synthetic and real datasets. As the first decentralized approach for constrained domains, our algorithm exhibits fast posterior concentration and high predictive accuracy.
Mohammad Rafiqul Islam, Lingjiong Zhu