Comparing the R algorithm and RHMC for staggered fermions

被引:11
作者
Clark, MA [1 ]
Joó, B [1 ]
Kennedy, AD [1 ]
机构
[1] Univ Edinburgh, Sch Phys, Edinburgh EH9 3JZ, Midlothian, Scotland
关键词
D O I
10.1016/S0920-5632(03)01748-1
中图分类号
O412 [相对论、场论]; O572.2 [粒子物理学];
学科分类号
摘要
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(deltatau(2)) errors, but show that the relevant expansion parameter is (deltaxi/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 deltatau approximate to m.
引用
收藏
页码:1015 / 1017
页数:3
相关论文
共 2 条
  • [1] HYBRID-MOLECULAR-DYNAMICS ALGORITHMS FOR THE NUMERICAL-SIMULATION OF QUANTUM CHROMODYNAMICS
    GOTTLIEB, S
    LIU, W
    TOUSSAINT, D
    RENKEN, RL
    SUGAR, RL
    [J]. PHYSICAL REVIEW D, 1987, 35 (08): : 2531 - 2542
  • [2] A new exact method for dynamical fermion computations with non-local actions
    Horváth, I
    Kennedy, AD
    Sint, S
    [J]. NUCLEAR PHYSICS B-PROCEEDINGS SUPPLEMENTS, 1999, 73 : 834 - 836