Low-lying Wilson Dirac operator eigenvector mixing in dynamical overlap hybrid Monte Carlo

Current dynamical overlap fermion hybrid Monte Carlo simulations encounter large fermionic forces when there is mixing between eigenvectors of the kernel operator with near zero-eigenvalues. This leads to low acceptance rates when there is a large density of near zero eigenvalues. I present a method where these large forces are eliminated and the large action jumps seen when two eigenvalues approach zero are significantly reduced. This significantly increases the stability of the algorithm, and allows the use of larger integration time steps.

[1]  K. Kanaya,et al.  Dynamical overlap fermion at fixed topology , 2006, hep-lat/0610011.

[2]  Herbert Neuberger.  Exactly massless quarks on the lattice , 1998 .

[3]  P. Weisz,et al.  On-shell improved lattice gauge theories , 1985 .

[4]  Klaus Schilling,et al.  Numerical methods for the QCD overlap operator IV: Hybrid Monte Carlo , 2009, Comput. Phys. Commun..

[5]  H. Matsufuru,et al.  Lattice gauge action suppressing near-zero modes of HW , 2006, hep-lat/0607020.

[6]  V. Rich.  Personal communication , 1989, Nature.

[7]  Z. Fodor,et al.  Dynamical overlap fermions, results with hybrid Monte-Carlo algorithm , 2003, hep-lat/0311010.

[8]  Colin Morningstar,et al.  Analytic smearing of SU(3) link variables in lattice QCD , 2004 .

[9]  Martin Lüscher,et al.  Exact chiral symmetry on the lattice and the Ginsparg-Wilson relation , 1998 .

[10]  S. Aoki,et al.  Two-flavor lattice QCD simulation in the epsilon-regime with exact chiral symmetry , 2007 .

[11]  Current status of Dynamical Overlap project , 2006, hep-lat/0511047.

[12]  S. Duane,et al.  Hybrid Monte Carlo , 1987 .

[13]  R. Narayanan,et al.  Chiral fermions on the lattice , 1993, hep-lat/9308011.

[14]  M. Lüscher.  Chiral gauge theories on the lattice with exact gauge invariance , 1999 .

[15]  P. Weisz,et al.  Computation of the action for on-shell improved lattice gauge theories at weak coupling , 1985 .

[16]  Giampiero Paffuti,et al.  Symanzik's improved lagrangian for lattice gauge theory , 1983 .

[17]  Thomas Lippert,et al.  Topological tunnelling with dynamical overlap fermions , 2008, Comput. Phys. Commun..

[18]  H. V. D. Vorst,et al.  Numerical methods for the QCDd overlap operator. I. Sign-function and error bounds , 2002, hep-lat/0202025.

[19]  Exact 2+1 flavour RHMC simulations , 2005, hep-lat/0409133.

[20]  Computation of the one-loop Symanzik coefficients for the square action , 1997, hep-lat/9701002.

[21]  J. Sexton,et al.  Hamiltonian evolution for the hybrid Monte Carlo algorithm , 1992 .

[22]  Speeding up the hybrid Monte Carlo algorithm for dynamical fermions , 2001, hep-lat/0107019.

[23]  Topology with dynamical overlap fermions , 2006, hep-lat/0510117.

[24]  Physics issues in simulations with dynamical overlap fermions , 2004, hep-lat/0412005.

[25]  S. Schaefer.  Algorithms for dynamical overlap fermions , 2006, hep-lat/0609063.