An exact Polynomial Hybrid Monte Carlo algorithm for dynamical Kogut-Susskind fermions

We present an exact dynamical QCD simulation algorithm for the $O(a)$-improved Wilson fermion with odd number of flavors. Our algorithm is an extension of the non-Hermitian polynomials HMC algorithm proposed by Takaishi and de Forcrand previously. In our algorithm, the systematic errors caused by the polynomial approximation of the inverse of Dirac operator is removed by a noisy-Metropolis test. For one flavor quark it is achieved by taking the square root of the correction m...

Since its first description fifty years ago, the Metropolis Monte Carlo method has been used in a variety of different ways for the simulation of continuum quantum many-body systems. This paper will consider some of the generalizations of the Metropolis algorithm employed in quantum Monte Carlo: Variational Monte Carlo, dynamical methods for projector monte carlo ({\it i.e.} diffusion Monte Carlo with rejection), multilevel sampling in path integral Monte Carlo, the sampling ...

There has been much recent progress in the understanding and reduction of the computational cost of the Hybrid Monte Carlo algorithm for Lattice QCD as the quark mass parameter is reduced. In this letter we present a new solution to this problem, where we represent the fermionic determinant using (n) pseudofermion fields, each with an (\nth) root kernel. We implement this within the framework of the Rational Hybrid Monte Carlo (RHMC) algorithm. We compare this algorithm with ...

Three topics concerning fermion simulation algorithms are discussed: 1.) A performance comparison of the multiboson technique to simulate dynamical fermions and the Kramers equation algorithm, 2.) the question of reversibility in the Hybrid Monte Carlo algorithm and 3.) the implementation of Symanzik's improvement program for dynamical Wilson fermions.

We performed dynamical simulations with HYP smeared staggered fermions using the recently proposed partial-global stochastic Metropolis algorithm with fermion matrix reduction and determinant breakup improvements. In this paper we discuss our choice of the action parameters and study the autocorrelation time both with four and two fermionic flavors at different quark mass values on approximately 10 fm^4 lattices. We find that the update is especially efficient with two flavor...

We offer a new proposal for the Monte Carlo treatment of many-fermion systems in continuous space. It is based upon Diffusion Monte Carlo with significant modifications: correlated pairs of random walkers that carry opposite signs; different functions ``guide'' walkers of different signs; the Gaussians used for members of a pair are correlated; walkers can cancel so as to conserve their expected future contributions. We report results for free-fermion systems and a fermion fl...

We present an exact Monte Carlo algorithm designed to sample theories where the energy is a sum of many couplings of decreasing strength. The algorithm avoids the computation of almost all non-leading terms. Its use is illustrated by simulating SU(2) lattice gauge theory with a 5-loop improved action. A new approach for dynamical fermion simulations is proposed.

We extend the dual algorithm recently described for pure, non-abelian Yang-Mills on the lattice to the case of lattice fermions coupled to Yang-Mills, by constructing an ergodic Metropolis algorithm for dynamic fermions that is local, exact, and built from gauge-invariant boson-fermion coupled configurations. For concreteness, we present in detail the case of three dimensions, for the group SU(2) and staggered fermions, however the algorithm readily generalizes with regard to...

A software library is presented for the polynomial expansion method (PEM) of the density of states (DOS) developed by two of the authors (N.F. and Y. M.). The library provides all necessary functions for the use of the PEM and its truncated version (TPEM) in a model independent way. The PEM/TPEM replaces the exact diagonalization of the one electron sector in models for fermions coupled to classical fields. The computational cost of the algorithm is O(N) -- with N the number ...

We present a detailed design of a (P)HMC simulation algorithm for N_f = 2 + 1 + 1 maximally twisted Wilson quark flavours. The algorithm retains even/odd and mass-shift preconditionings combined with multiple Molecular Dynamics time scales for both the light mass degenerate, u and d, quarks and the heavy mass non-degenerate, s and c, quarks. Various non-standard aspects of the algorithm are discussed, among which those connected to the use of a polynomial approximation for th...