High-Order Langevin Diffusion Yields an Accelerated MCMC Algorithm
Abstract
We propose a Markov chain Monte Carlo (MCMC) algorithm based on third-order Langevin dynamics for sampling from distributions with log-concave and smooth densities. The higher-order dynamics allow for more flexible discretization schemes, and we develop a specific method that combines splitting with more accurate integration. For a broad class of -dimensional distributions arising from generalized linear models, we prove that the resulting third-order algorithm produces samples from a distribution that is at most in Wasserstein distance from the target distribution in steps. This result requires only Lipschitz conditions on the gradient. For general strongly convex potentials with -th order smoothness, we prove that the mixing time scales as .
Cite
@article{arxiv.1908.10859,
title = {High-Order Langevin Diffusion Yields an Accelerated MCMC Algorithm},
author = {Wenlong Mou and Yi-An Ma and Martin J. Wainwright and Peter L. Bartlett and Michael I. Jordan},
journal= {arXiv preprint arXiv:1908.10859},
year = {2020}
}
Comments
Changes from v1: improved algorithm with $O (d^{1/4} / \varepsilon^{1/2})$ mixing time