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

The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on lattices of size up to 32x24x24x24, with lattice spacings a~0.08 fm and at bare current-quark masses as low as 21 MeV.

We present preliminary results of the PACS-CS project which simulates 2+1 flavor lattice QCD toward the physical point with the nonperturbatively O(a)-improved Wilson quark action and the Iwasaki gauge action. Calculations are carried out at beta=1.9 on a 32^3x64 lattice with the use of the domain-decomposed HMC algorithm to reduce the up-down quark mass. The resulting pseudoscalar meson masses range from 730 MeV down to 210 MeV. We discuss the physical results including the ...

We present an algorithmic study for the simulation of two massless flavors of O(a) improved Wilson quarks with Schroedinger functional boundary conditions. The algorithm used is Hybrid Monte Carlo with two pseudo-fermion fields as proposed by M. Hasenbusch. A gain in CPU cost of a factor two is reached when compared to one pseudo-fermion field due to the larger possible step-size. This study is integrated in the ALPHA project for the computation of the running of the renormal...

We present an exact local bosonic algorithm for the simulation of dynamical fermions in lattice QCD. We show that this algorithm is a viable alternative to the Hybrid Monte Carlo algorithm.

This is the write-up of three lectures on algorithms for dynamical fermions that were given at the ILFTN workshop 'Perspectives in Lattice QCD' in Nara during November 2005. The first lecture is on the fundamentals of Markov Chain Monte Carlo methods and introduces the Hybrid Monte Carlo (HMC) algorithm and symplectic integrators; the second lecture covers topics in approximation theory and thereby introduces the Rational Hybrid Monte Carlo (RHMC) algorithm and ways of evadin...

We report on our simulation of finite temperature lattice QCD with two flavors of ${\cal O}(a)$ Symanzik-improved fermions and ${\cal O}(a^2)$ Symanzik-improved glue. Our thermodynamic simulations were performed on an $8^3 \times 4$ lattice, and we have performed complementary zero temperature simulations on an $8^3 \times 16$ lattice. We compare our results to those from simulations with two flavors of Wilson fermions and discuss the improvement resulting from use of the imp...

As algorithmic developments have driven down the cost of simulating degenerate light quark flavors the relative cost of simulating single quark flavors with the Rational Hybrid Monte Carlo (RHMC) algorithm has become more expensive. TWQCD has proposed an exact one-flavor algorithm (EOFA) that allows for HMC simulations of a single quark flavor without taking a square root of the fermion determinant. We have independently implemented EOFA in the Columbia Physics System (CPS) a...

We present details of our investigation of the Parallel Tempering algorithm. We consider the application of action matching technology to the selection of parameters. We then present a simple model of the autocorrelations for a particular parallel tempered system. Finally we present numerical results from applying the algorithm to lattice QCD with O(a)-improved dynamical fermions with twin ensemble systems

Lattice QCD calculations including the effects of one or more non-degenerate sea quark flavors are conventionally performed using the Rational Hybrid Monte Carlo (RHMC) algorithm, which computes the square root of the determinant of $\mathscr{D}^{\dagger} \mathscr{D}$, where $\mathscr{D}$ is the Dirac operator. The special case of two degenerate quark flavors with the same mass is described directly by the determinant of $\mathscr{D}^{\dagger} \mathscr{D}$ --- in particular, ...

The standard hybrid Monte Carlo algorithm is known to simulate even flavors QCD only. Simulations of odd flavors QCD, however, can be also performed in the framework of the hybrid Monte Carlo algorithm where the inverse of the fermion matrix is approximated by a polynomial. In this exploratory study we perform three flavors QCD simulations. We make a comparison of the hybrid Monte Carlo algorithm and the R-algorithm which also simulates odd flavors systems but has step-size e...