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

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.