Comparing the R algorithm and RHMC for staggered fermions

DOI10.1016/S0920-5632(03)01748-1MaRDI QIDQ1410028zbMATH OpenFDO

Authors Michael A. Clark, Bálint Joó, A. D. Kennedy

Publication date 13 October 2003

Published in Nuclear Physics. B. Proceedings Supplements (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/hep-lat/0209035

zbMATH Keywords

lattice QCD simulationsprobability distribution

Mathematics Subject Classification ID

Quantum field theory on lattices (81T25) Strong interaction, including quantum chromodynamics (81V05)

Abstract: The R algorithm is widely used for simulating two flavours of dynamical staggered fermions. We give a simple proof that the algorithm converges to the desired probability distribution to within O(dt^2) errors, but show that the relevant expansion parameter is (dt/m)^2, m being the quark mass. The Rational Hybrid Monte Carlo (RHMC) algorithm provides an exact (i.e., has no step size errors) alternative for simulating the square root of the staggered Dirac operator. We propose using it to test the validity of the R algorithm for simulations carried out with dt m.

Recommendations

• A new exact method for dynamical fermion computations with non-local actions
• An exact Polynomial Hybrid Monte Carlo algorithm for dynamical Kogut-Susskind fermions
• Numerical methods for the QCD overlap operator IV: Hybrid Monte Carlo
• Speeding up the hybrid Monte Carlo algorithm for dynamical fermions
• A new way to deal with fermions in Monte Carlo simulations

Cited in

(2)

• Chiral limit of staggered fermions at strong couplings: A loop representation
• A new exact method for dynamical fermion computations with non-local actions

This page was built for publication: Comparing the R algorithm and RHMC for staggered fermions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1410028)