Abstract
We investigate instability and reversibility within hybrid Monte Carlo simulations using a nonperturbatively improved Wilson action. We demonstrate the onset of instability as tolerance parameters and molecular dynamics step sizes are varied. We compare these findings with theoretical expectations and present limits on simulation parameters within which a stable and reversible algorithm is obtained for physically relevant simulations. Results of optimization experiments with respect to tolerance parameters are also presented.
- Received 26 May 2000
DOI:https://doi.org/10.1103/PhysRevD.62.114501
©2000 American Physical Society