Instability in the Molecular Dynamics Step of Hybrid Monte Carlo in Dynamical Fermion Lattice QCD Simulations

Abstract

We investigate instability and reversibility within Hybrid Monte Carlo simulations using a non-perturbatively 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 optimisation experiments with respect to tolerance prarameters are also presented.

arXiv commentRevTeX, Some results here were presented at Vielat 99, Vienna, Austria, Sept 1999 22 Pages, 10 figures, to be submitted to PRD

References in corpus (5)

• A New Way to Set the Energy Scale in Lattice Gauge Theories and its Application to the Static Force and ${\alpha }_s$ in SU(2) Yang--Mills Theory
• Instabilities and Non-Reversibility of Molecular Dynamics Trajectories
• Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms
• The non-perturbative O(a)-improved action for dynamical Wilson fermions
• Implementing the Generalised Hybrid Monte-Carlo Algorithm

Cited by in corpus (28)

• 2+1 Flavor Lattice QCD toward the Physical Point
• Accelerating Dynamical Fermion Computations using the Rational Hybrid Monte Carlo (RHMC) Algorithm with Multiple Pseudofermion Fields
• Higher representations on the lattice: numerical simulations. SU(2) with adjoint fermions
• Speeding up Lattice QCD simulations with clover-improved Wilson Fermions
• Effects of non-perturbatively improved dynamical fermions in QCD at fixed lattice spacing
• Two-flavor QCD simulation with exact chiral symmetry
• Hadron Spectroscopy with Dynamical Chirally Improved Fermions
• Numerical integrators for the Hybrid Monte Carlo method
• Cost of the Generalised Hybrid Monte Carlo Algorithm for Free Field Theory
• Improved Wilson QCD simulations with light quark masses
• Polynomial Hybrid Monte Carlo algorithm for lattice QCD with an odd number of flavors
• Exploring the HMC trajectory-length dependence of autocorrelation times in lattice QCD
• Algorithms for Lattice QCD with Dynamical Fermions
• Light hadron spectroscopy in two-flavor QCD with small sea quark masses
• R-symmetry in the Q-exact (2,2) 2d lattice Wess-Zumino model
• Roberge-Weiss transition in $N_{\text{f}}=2$ QCD with Wilson fermions and $N_{\tau }=6$
• QCD simulations with staggered fermions on GPUs
• Simulating the Schroedinger functional with two pseudo-fermions
• Wilson fermions with chirally twisted mass
• Progress in lattice algorithms
• A Numerical Study of the 2-Flavour Schwinger Model with Dynamical Overlap Hypercube Fermions
• Suppressing dislocations in normalized hypercubic smearing
• Impact of large cutoff-effects on algorithms for improved Wilson fermions
• Upper and lower Higgs boson mass bounds from a chirally invariant lattice Higgs-Yukawa model
• The RHMC algorithm for theories with unknown spectral bounds
• Hessian-free force-gradient integrators
• Cutoff-effects in the spectrum of dynamical Wilson fermions
• Numerical stability of force-gradient integrators and their Hessian-free variants in lattice QCD simulations