About HMC
HMC augments position q with momentum p, then simulates Hamiltonian dynamics:
H(q,p) = U(q) + K(p)
The leapfrog integrator approximately conserves H, making proposed samples nearly always accepted. Unlike random-walk MCMC, trajectories explore efficiently with low correlation.