Force-Gradient Nested Multirate Methods for Hamiltonian System
Abstract
Force-gradient decomposition methods are used to improve the energy preservation of symplectic schemes applied to Hamiltonian systems. If the potential is composed of different parts with strongly varying dynamics, this multirate potential can be exploited by coupling force-gradient decomposition methods with splitting techniques for multi-time scale problems to further increase the accuracy of the scheme and reduce the computational costs. In this paper, we derive novel force-gradient nested methods and test them numerically. Such methods can be used to increase the acceptance rate for the molecular dynamics step of the Hybrid Monte Carlo algorithm (HMC) and hence improve its computational efficiency.
Cite
@article{arxiv.1312.3113,
title = {Force-Gradient Nested Multirate Methods for Hamiltonian System},
author = {Dmitry Shcherbakov and Matthias Ehrhardt and Michael Günther and Michael Peardon},
journal= {arXiv preprint arXiv:1312.3113},
year = {2013}
}
Comments
16 pages 3 figures, submitted to Applied Numerical Mathematics