We develop an implementation for a recently proposed noisy Monte Carlo approach to the simulation of lattice QCD with dynamical fermions by incorporating the full fermion determinant directly. Our algorithm uses a quenched gauge field update with a shifted gauge coupling to minimize fluctuations in the trace of the logarithm of the Wilson-Dirac matrix. The details of tuning the gauge coupling shift as well as results for the distribution of noisy estimators in our implementation are given. We present data for some basic observables from the noisy method, as well as the acceptance rate information and discuss potential autocorrelation and sign violation effects. Both the results and the efficiency of the algorithm are compared against those of the hybrid Monte Carlo algorithm.
|Journal||Physical Review D - Particles, Fields, Gravitation and Cosmology|
|State||Published - 2003|
ASJC Scopus subject areas
- Nuclear and High Energy Physics
- Physics and Astronomy (miscellaneous)