Hessian-free force-gradient integrators
arXiv:2403.10370 · doi:10.1016/j.cpc.2024.109478
Abstract
We propose a new framework of Hessian-free force-gradient integrators that do not require the analytical expression of the force-gradient term based on the Hessian of the potential. Due to that the new class of decomposition algorithms for separable Hamiltonian systems with quadratic kinetic energy may be particularly useful when applied to Hamiltonian systems where an evaluation of the Hessian is significantly more expensive than an evaluation of its gradient, e.g. in molecular dynamics simulations of classical systems. Numerical experiments of an N-body problem, as well as applications to the molecular dynamics step in the Hybrid Monte Carlo (HMC) algorithm for lattice simulations of the Schwinger model and Quantum Chromodynamics (QCD) verify these expectations.
28 pages, 6 figures
References in corpus (9)
- Speeding up the Hybrid-Monte-Carlo algorithm for dynamical fermions
- Lattice QCD with open boundary conditions and twisted-mass reweighting
- HMC algorithm with multiple time scale integration and mass preconditioning
- Testing and tuning symplectic integrators for Hybrid Monte Carlo algorithm in lattice QCD
- On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
- Shadow Hamiltonians, Poisson Brackets, and Gauge Theories
- Scaling test of fermion actions in the Schwinger model
- Optimising creation operators for charmonium spectroscopy on the lattice
- Domain Wall Fermion QCD with the Exact One Flavor Algorithm