Litcius/Paper detail

Mixing time guarantees for unadjusted Hamiltonian Monte Carlo

Nawaf Bou‐Rabee, Andreas Eberle

2022Bernoulli20 citationsDOIOpen Access PDF

Abstract

We provide quantitative upper bounds on the total variation mixing time of the Markov chain corresponding to the unadjusted Hamiltonian Monte Carlo (uHMC) algorithm. For two general classes of models and fixed time discretization step size h, the mixing time is shown to depend only logarithmically on the dimension. Moreover, we provide quantitative upper bounds on the total variation distance between the invariant measure of the uHMC chain and the true target measure. As a consequence, we show that an ε-accurate approximation of the target distribution μ in total variation distance can be achieved by uHMC: (i) for a broad class of models with Od3∕4ε−1∕2log(d∕ε) gradient evaluations; and (ii) for mean field models with weak interactions with Od1∕2ε−1∕2log(d∕ε) gradient evaluations. The proofs are based on the construction of successful couplings for uHMC that realize the upper bounds.

Topics & Concepts

MathematicsMarkov chain Monte CarloMixing (physics)Statistical physicsUpper and lower boundsMarkov chainMonte Carlo methodDiscretizationCombinatoricsApplied mathematicsMathematical analysisStatisticsPhysicsQuantum mechanicsMarkov Chains and Monte Carlo MethodsMedical Imaging Techniques and ApplicationsAdvanced MRI Techniques and Applications
Mixing time guarantees for unadjusted Hamiltonian Monte Carlo | Litcius