This paper provides a convergence analysis for generalized Hamiltonian Monte Carlo samplers, a family of Markov Chain Monte Carlo methods based on leapfrog integration of Hamiltonian dynamics and kinetic Langevin diffusion, that encompasses the unadjusted Hamiltonian Monte Carlo method. Assuming that the target distribution $\pi$ satisfies a log-Sobolev inequality and mild conditions on the corresponding potential function, we establish quantitative bounds on the relative entropy of the iterates defined by the algorithm, with respect to $\pi$. Our approach is based on a perturbative and discrete version of the modified entropy method developed to establish hypocoercivity for the continuous-time kinetic Langevin process. As a corollary of our main result, we are able to derive complexity bounds for the class of algorithms at hand. In particular, we show that the total number of iterations to achieve a target accuracy $\varepsilon >0$ is of order $d/\varepsilon^{1/4}$, where $d$ is the dimension of the problem. This result can be further improved in the case of weakly interacting mean field potentials, for which we find a total number of iterations of order $(d/\varepsilon)^{1/4}$.
翻译:本文针对广义哈密顿蒙特卡洛采样器提供收敛性分析,该采样器是一类基于哈密顿动力学跳蛙积分和动力学朗之万扩散的马尔可夫链蒙特卡洛方法,包含了非调整哈密顿蒙特卡洛方法。假设目标分布 $\pi$ 满足对数索博列夫不等式且对应势函数满足温和条件,我们建立了算法迭代点相对于 $\pi$ 的相对熵的定量界。我们的方法基于扰动离散版本的修正熵方法,该方法最初为建立连续时间动力学朗之万过程的亚强制性而发展。作为主要结果的推论,我们能够推导出该类算法的复杂度界。特别地,我们证明达到目标精度 $\varepsilon >0$ 所需的总迭代次数阶为 $d/\varepsilon^{1/4}$,其中 $d$ 是问题的维度。在弱相互作用平均场势的情况下,该结果可进一步优化,我们得到总迭代次数阶为 $(d/\varepsilon)^{1/4}$。