Kramers equation algorithm for simulations of QCD with two flavors of Wilson fermions and gauge group SU(2)

Abstract

We compare the Hybrid Monte Carlo (HMC) and the Kramers equation algorithms for simulations of QCD with two flavors of dynamical Wilson fermions and gauge group $SU(2)$. The results for the performance of both algorithms are obtained on $6^{3}12$, $12^4$ and $16^4$ lattices at a pion to $\rho$ meson mass ratio of $m_{\pi }/m_{\rho }\approx 0.9$. We find that the Kramers equation algorithm gives an equally good performance as the HMC algorithm. We demonstrate that the classical equations of motion used in these algorithms lack reversibility in practical simulations and behave like those of a chaotic dynamical system with a Liapunov exponent $\nu \approx 0.75$.

20pages , uuencoded postscript file

References in corpus (2)

• A new simulation algorithm for lattice QCD with dynamical quarks
• Study of a new simulation algorithm for dynamical quarks on the APE-100 parallel computer

Cited by in corpus (29)

• Lattice QCD without topology barriers
• A Polynomial Hybrid Monte Carlo Algorithm
• Stochastic locality and master-field simulations of very large lattices
• Master-field simulations of O(a)-improved lattice QCD: Algorithms, stability and exactness
• Implementation of Symanzik's Improvement Program for Simulations of Dynamical Wilson Fermions in Lattice QCD
• The PHMC algorithm for simulations of dynamical fermions: I -- description and properties
• Spectral function at high temperature in the classical approximation
• Comparative Benchmarks of full QCD Algorithms
• Simulating the Schroedinger functional with two pseudo-fermions
• Instabilities and Non-Reversibility of Molecular Dynamics Trajectories
• Improvements of Lüscher's local bosonic fermion algorithm
• Non-renormalizability of the HMC algorithm
• Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms
• Investigation of New Methods for Numerical Stochastic Perturbation Theory in ${\phi }^4$ Theory
• Recent Developments in Fermion Simulation Algorithms
• Improved Pseudofermion Approach for All-Point Propagators
• Performance Tests of the Kramers Equation and Boson Algorithms for Simulations of QCD
• Critical Dynamics of the Hybrid Monte Carlo Algorithm
• Study of Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms
• Ordering monomial factors of polynomials in the product representation
• Improvements of the local bosonic algorithm
• Unquenched QCD with Light Quarks
• Reversibility Violation in the Hybrid Monte Carlo Algorithm
• Kramers Equation Algorithm with Kogut-Susskind Fermions on Lattice
• Dynamic relaxation of SU(2) lattice gauge theory in (3+1) dimensions
• Hubbard Model with Lüscher fermions
• Can we do better than Hybrid Monte Carlo in Lattice QCD?
• Non-Reversibility of Molecular Dynamics Trajectories
• Progress in generating gauge ensembles with Stabilized Wilson Fermions