Numerical integrators for the Hybrid Monte Carlo method

Abstract

We construct numerical integrators for Hamiltonian problems that may advantageously replace the standard Verlet time-stepper within Hybrid Monte Carlo and related simulations. Past attempts have often aimed at boosting the order of accuracy of the integrator and/or reducing the size of its error constants; order and error constant are relevant concepts in the limit of vanishing step-length. We propose an alternative methodology based on the performance of the integrator when sampling from Gaussian distributions with not necessarily small step-lengths. We construct new splitting formulae that require two, three or four force evaluations per time-step. Limited, proof-of-concept numerical experiments suggest that the new integrators may provide an improvement on the efficiency of the standard Verlet method, especially in problems with high dimensionality.

30 pages, 5 figures

Cited by in corpus (19)

• Randomized Hamiltonian Monte Carlo
• Geometric integrators and the Hamiltonian Monte Carlo method
• Compressible Generalized Hybrid Monte Carlo
• Modified Hamiltonian Monte Carlo for Bayesian inference
• Extra Chance Generalized Hybrid Monte Carlo
• Adaptive multi-stage integrators for optimal energy conservation in molecular simulations
• Palindromic 3-stage splitting integrators, a roadmap
• Convergence of unadjusted Hamiltonian Monte Carlo for mean-field models
• Multi-stage splitting integrators for sampling with modified Hamiltonian Monte Carlo methods
• Couplings for Andersen Dynamics
• Quantum Monte Carlo Simulation of the 3D Ising Transition on the Fuzzy Sphere
• Higher Order Hamiltonian Monte Carlo Sampling for Cosmological Large-Scale Structure Analysis
• Geometric Methods for Sampling, Optimisation, Inference and Adaptive Agents
• Composition Methods for Dynamical Systems Separable into Three Parts
• Adaptive multi-stage integration schemes for Hamiltonian Monte Carlo
• Unadjusted Hamiltonian MCMC with Stratified Monte Carlo Time Integration
• Precision annealing Monte Carlo methods for statistical data assimilation and machine learning
• Benchmarking field-level cosmological inference from galaxy redshift surveys
• Differentiable Fuzzy Cosmic-Web for Field Level Inference