English

Irreversible Langevin MCMC on Lie Groups

Statistics Theory 2019-03-22 v1 Differential Geometry Statistics Theory

Abstract

It is well-known that irreversible MCMC algorithms converge faster to their stationary distributions than reversible ones. Using the special geometric structure of Lie groups G\mathcal G and dissipation fields compatible with the symplectic structure, we construct an irreversible HMC-like MCMC algorithm on G\mathcal G, where we first update the momentum by solving an OU process on the corresponding Lie algebra g\mathfrak g, and then approximate the Hamiltonian system on G×g\mathcal G \times \mathfrak g with a reversible symplectic integrator followed by a Metropolis-Hastings correction step. In particular, when the OU process is simulated over sufficiently long times, we recover HMC as a special case. We illustrate this algorithm numerically using the example G=SO(3)\mathcal G = SO(3).

Keywords

Cite

@article{arxiv.1903.08939,
  title  = {Irreversible Langevin MCMC on Lie Groups},
  author = {Alexis Arnaudon and Alessandro Barp and So Takao},
  journal= {arXiv preprint arXiv:1903.08939},
  year   = {2019}
}
R2 v1 2026-06-23T08:14:53.333Z