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

The exchange antisymmetry between identical fermions gives rise to the well known fermion sign problem, in the form of large cancellation between positive and negative contribution to the partition function, making any simulation methods which directly sample this partition function exponentially difficult to converge. In this work, we employ path integral molecular dynamics (PIMD) and build upon the recently discovered fictitious particle model to investigate the fermion sig...

An improved algorithm is proposed for Monte Carlo methods to study fermion systems interacting with adiabatical fields. To obtain a weight for each Monte Carlo sample with a fixed configuration of adiabatical fields, a series expansion using Chebyshev polynomials is applied. By introducing truncations of matrix operations in a systematic and controlled way, it is shown that the cpu time is reduced from O(N^3) to O(N) where N is the system size. Benchmark results show that the...

We discuss recent algorithmic improvements in simulating finite temperature QCD on a lattice. In particular, the Rational Hybrid Monte Carlo(RHMC) algorithm is employed to generate lattice configurations for 2+1 flavor QCD. Unlike the Hybrid R Algorithm, RHMC is reversible, admitting a Metropolis accept/reject step that eliminates the $\mathcal{O}(\delta t^2)$ errors inherent in the R Algorithm. We also employ several algorithmic speed-ups, including multiple time scales, the...

We consider recent progress in algorithms for generating gauge field configurations that include the dynamical effects of light fermions. We survey what has been achieved in recent state-of-the-art computations, and examine the trade-offs between performance and control of systematic errors. We briefly review the use of polynomial and rational approximations in Hybrid Monte Carlo algorithms, and some of the theory of on-shell chiral fermions on the lattice. This provides a th...

We develop the self-learning Monte Carlo (SLMC) method, a general-purpose numerical method recently introduced to simulate many-body systems, for studying interacting fermion systems. Our method uses a highly-efficient update algorithm, which we design and dub "cumulative update", to generate new candidate configurations in the Markov chain based on a self-learned bosonic effective model. From general analysis and numerical study of the double exchange model as an example, we...

The advantages of using Multi-Step corrections for simulations of lattice gauge theories with dynamical fermions will be discussed. This technique is suited for algorithms based on the Multi-Boson representation of the dynamical fermions as well as for the Hybrid Monte-Carlo (HMC) algorithm and variants of the latter, like the Polynomial-HMC. Especially the latter has the power to deal with an odd number of fermion fields--an essential feature necessary for realistic QCD-simu...

This review gives an overview on the research of algorithms for dynamical fermions used in large scale lattice QCD simulations. First a short overview on the state-of-the-art of ensemble generation at the physical point is given. Followed by an overview on necessary steps towards simulation of large lattices with the Hybrid Monte Carlo algorithm. Here, the status of iterative solvers and tuning procedures for numerical integrators within the molecular dynamics are discuss...

In this paper we propose new algorithm to reduce autocorrelation in Markov chain Monte-Carlo algorithms for euclidean field theories on the lattice. Our proposing algorithm is the Hybrid Monte-Carlo algorithm (HMC) with restricted Boltzmann machine. We examine the validity of the algorithm by employing the phi-fourth theory in three dimension. We observe reduction of the autocorrelation both in symmetric and broken phase as well. Our proposing algorithm provides consistent ce...

In the present paper we examine the effects of noise on Monte Carlo algorithms, a problem raised previously by Kennedy and Kuti (Phys. Rev. Lett. {\bf 54}, 2473 (1985)). We show that the effects of introducing unbiased noise into the acceptance/rejection phase of the conventional Metropolis approach are surprisingly modest, and, to a significant degree, largely controllable. We present model condensed phase numerical applications to support these conclusions.

The past few years have seen considerable progress in algorithmic development for the generation of gauge fields including the effects of dynamical fermions. The Rational Hybrid Monte Carlo (RHMC) algorithm, where Hybrid Monte Carlo is performed using a rational approximation in place the usual inverse quark matrix kernel is one of these developments. This algorithm has been found to be extremely beneficial in many areas of lattice QCD (chiral fermions, finite temperature, Wi...