Comparing the R algorithm and RHMC for staggered fermions

M. A. Clark; B. Joó; A. D. Kennedy

doi:10.1016/S0920-5632(03)01748-1

Comparing the R algorithm and RHMC for staggered fermions

M. A. Clark, B. Joó, A. D. Kennedy

Research output: Contribution to journal › Article › peer-review

17 Scopus citations

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(δτ²) errors, but show that the relevant expansion parameter is (δτ/m)², 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 δτ ≈ m.

Original language	English
Pages (from-to)	1015-1017
Number of pages	3
Journal	Nuclear Physics B - Proceedings Supplements
Volume	119
DOIs	https://doi.org/10.1016/S0920-5632(03)01748-1
State	Published - May 2003
Externally published	Yes

Funding

Our thanks go to Robert, Edwards, Ivan HorvLth, Stefan Sint, and Tlrs Wenger for discussions and for providing the RHMC’ code in the SZIN soft,ware system[4]. Work supported by the Europeau (,‘omnlunity’s Human potential programme under HPRN-CT-2000-00145 Hadrons/Lat,ticeQ(‘D.

Funders	Funder number
Europeau	HPRN-CT-2000-00145 Hadrons/Lat

Access to Document

10.1016/S0920-5632(03)01748-1

Cite this

@article{c75fd04894d84440aec5d799479a372b,

title = "Comparing the R algorithm and RHMC for staggered fermions",

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(δτ2) errors, but show that the relevant expansion parameter is (δτ/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 δτ ≈ m.",

author = "Clark, {M. A.} and B. Jo{\'o} and Kennedy, {A. D.}",

year = "2003",

month = may,

doi = "10.1016/S0920-5632(03)01748-1",

language = "English",

volume = "119",

pages = "1015--1017",

journal = "Nuclear Physics B - Proceedings Supplements",

issn = "0920-5632",

publisher = "Elsevier",

}

TY - JOUR

T1 - Comparing the R algorithm and RHMC for staggered fermions

AU - Clark, M. A.

AU - Joó, B.

AU - Kennedy, A. D.

PY - 2003/5

Y1 - 2003/5

N2 - 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(δτ2) errors, but show that the relevant expansion parameter is (δτ/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 δτ ≈ m.

AB - 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(δτ2) errors, but show that the relevant expansion parameter is (δτ/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 δτ ≈ m.

UR - http://www.scopus.com/inward/record.url?scp=0037947768&partnerID=8YFLogxK

U2 - 10.1016/S0920-5632(03)01748-1

DO - 10.1016/S0920-5632(03)01748-1

M3 - Article

AN - SCOPUS:0037947768

SN - 0920-5632

VL - 119

SP - 1015

EP - 1017

JO - Nuclear Physics B - Proceedings Supplements

JF - Nuclear Physics B - Proceedings Supplements

ER -

Comparing the R algorithm and RHMC for staggered fermions

Abstract

Funding

Access to Document

Other files and links

Fingerprint

Cite this