Answers are grounded in this paper and cite their sources.
We present a simple recursive iteration of the leapfrog discretization of Newton's equations which leads to a removal of the finite-step-size error to any desired order. This is done in a manner that preserves phase-space areas and reversibility, as required for use in the hybrid Monte Carlo method for simulating fermionic fields. The resulting asymptotic volume dependence is exp[(lnV)1/2]. We test the scheme on the (2+1)-dimensional Hubbard model.
No takes yet. Share an insight, caveat, or question.
Creutz et al. (1989) studied this question.
Synapse has enriched 5 closely related papers on similar clinical questions. Consider them for comparative context: