An exact algorithm for three-flavor QCD with $O(a)$-improved Wilson fermions

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 log 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 give an update on our current project to determine the transition temperature and the order of the deconfinement transition in the chiral limit of two flavour QCD. We use nonperturbatively O(a) improved Wilson fermions of the Sheikholeslami-Wohlert type, employing the efficient deflation accelerated DDHMC algorithm. We start at lattices with N_t>=12 and pion masses below 600 MeV, aiming at chiral and continuum limits with light quarks.

We present a new exact algorithm for estimating all elements of the quark propagator. The advantage of the method is that the exact all-to-all propagator is reproduced in a large but finite number of inversions. The efficacy of the algorithm is tested in Monte Carlo simulations of Wilson quarks in quenched QCD. Applications that are difficult to probe with point propagators are discussed.

An exact, nonlocal, finite step-size algorithm for Monte Carlo simulation of theories with dynamical fermions is proposed. The algorithm is based on obtaining the new configuration U' from the old one U by solving the equation $ M(U') \eta = \omega M(U) \eta$, where $M$ is fermionic operator, $\eta$ is random Gaussian vector, and $\omega$ is random real number close to unity. This algorithm can be used for the acceleration of current simulations in theories with Grassmann var...

We simulate Quantum Chromodynamics (QCD) in four Euclidean dimensions with two (degenerate mass) flavors of dynamical quarks. The Dirac operator we use is the so-called chirally improved Dirac operator. We discuss the algorithm used for the simulation as well as the checks and some results on lattices up to size 8**4 for fermion mass parameters down to 0.1. This is the first attempt to introduce dynamical quarks with the chirally improved Dirac operator.

Considering Ginsparg-Wilson type fermions dynamically in lattice QCD simulations is a challenging task. The hope is to be able to approach smaller pion masses and to eventually reach physical situations. The price to pay is substantially higher computational costs. Here we discuss first results of a dynamical implementation of the so-called Chirally Improved Fermions, a Dirac operator that obeys the Ginsparg-Wilson condition approximately. The simulation is for two species of...

We propose a new method for Hybrid Monte Carlo (HMC) simulations with odd numbers of dynamical fermions on the lattice. It employs a different approach from polynomial or rational HMC. In this method, gamma-five hermiticity of the lattice Dirac operators is crucial and it can be applied to Wilson, domain-wall, and overlap fermions. We compare HMC simulations with two degenerate flavors and (1 + 1) degenerate flavors using optimal domain-wall fermions. The ratio of the efficie...

We present the results of the physical point simulation in 2+1 flavor lattice QCD with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $\beta=1.9$ on a $32^3 \times 64$ lattice. The physical quark masses together with the lattice spacing is determined with $m_\pi$, $m_K$ and $m_\Omega$ as physical inputs. There are two key algorithmic ingredients to make possible the direct simulation at the physical point: One is the mass-precondit...

I report on first results of the ETMC obtained simulating lattice QCD with two degenerate flavours of Wilson twisted mass fermions, for which physical observables are automatically O(a) improved. Recent improvements of the HMC algorithm allow simulations of pseudoscalar masses below 300 MeV on volumes L^3 x 2L with L>2 fm and at values of the lattice spacing a<0.1 fm.

We propose various improvements of finite step-size updating for full QCD on the lattice that might turn finite step-size updating into a viable alternative to the hybrid Monte Carlo algorithm. These improvements are noise reduction of the noisy estimator of the fermion determinant, unbiased inclusion of the hopping parameter expansion and a multi-level Metropolis scheme. First numerical tests are performed for the 2 dimensional Schwinger model with two flavours of Wilson fer...