A simple method to optimize HMC performance

Preprint number: CP3-Origins-2016-38 DNRF90
Authors: Andrea Bussone (CP3-Origins & DIAS), Michele Della Morte (CP3-Origins & DIAS), Vincent Drach (CERN), Martin Hansen (CP3-Origins & DIAS), Ari Hietanen (CP3-Origins & DIAS), Claudio Pica (CP3-Origins & DIAS), and Jarno Rantaharju (CP3-Origins & DIAS)
External link: arXiv.org

We present a practical strategy to optimize a set of Hybrid Monte Carlo parameters in simulations of QCD and QCD-like theories. We specialize to the case of mass-preconditioning, with multiple time-step Omelyan integrators. Starting from properties of the shadow Hamiltonian we show how the optimal setup for the integrator can be chosen once the forces and their variances are measured, assuming that those only depend on the mass-preconditioning parameter.