Hamiltonian Monte Carlo

Hamiltonian Monte Carlo introduces an auxiliary Gaussian momentum, approximately follows Hamiltonian dynamics by a reversible volume-preserving integrator, and applies a Metropolis correction.

Surrogate Hamiltonian Monte Carlo uses an auxiliary smooth probability density function $\nu$ to generate trajectories but applies a Metropolis–Hastings acceptance probability using the intended target $\mu$. For Gaussian momenta and an involutive Metropolis proposal generated by flipped leapfrog integration, the acceptance probability is $\min (1,\mu (x^{\prime })e^{-\Vert p^{\prime }{\Vert }^2/2}/[\mu (x)e^{-\Vert p{\Vert }^2/2}])$. Only the surrogate gradient is required during the trajectory; target density evaluations are still required at the endpoints.