Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms
arXiv:hep-lat/9708017 · doi:10.1016/S0550-3213(98)00174-6
Abstract
We study the phenomenon of lack of reversibility in molecular dynamics algorithms for the case of Wilson's lattice QCD. We demonstrate that the classical equations of motion that are employed in these algorithms are chaotic in nature. The leading Liapunov exponent is determined in a range of coupling parameters. We give a quantitative estimate of the consequences of the breakdown of reversibility due to round-off errors.
Latex2e file, 4 figures, 19 pages
References in corpus (4)
- Chronological Inversion Method for the Dirac Matrix in Hybrid Monte Carlo
- Kramers equation algorithm for simulations of QCD with two flavors of Wilson fermions and gauge group SU(2)
- Instabilities and Non-Reversibility of Molecular Dynamics Trajectories
- Study of Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms
Cited by in corpus (12)
- HMC algorithm with multiple time scale integration and mass preconditioning
- Speeding up Lattice QCD simulations with clover-improved Wilson Fermions
- 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
- Instability in the Molecular Dynamics Step of Hybrid Monte Carlo in Dynamical Fermion Lattice QCD Simulations
- Simulating the Schroedinger functional with two pseudo-fermions
- Instabilities and Non-Reversibility of Molecular Dynamics Trajectories
- Progress in lattice algorithms
- A Numerical Study of the 2-Flavour Schwinger Model with Dynamical Overlap Hypercube Fermions
- Upper and lower Higgs boson mass bounds from a chirally invariant lattice Higgs-Yukawa model
- Reversibility Violation in the Hybrid Monte Carlo Algorithm