Related papers: Comparing the R algorithm and RHMC for staggered f…
We describe a Fourier Accelerated Hybrid Monte Carlo algorithm suitable for dynamical fermion simulations of non-gauge models. We test the algorithm in supersymmetric quantum mechanics viewed as a one-dimensional Euclidean lattice field…
Simulating from the multivariate truncated normal distribution (MTN) is required in various statistical applications yet remains challenging in high dimensions. Currently available algorithms and their implementations often fail when the…
The hybrid Monte Carlo (HMC) algorithm is arguably the most efficient sampling method for general probability distributions of continuous variables. Together with exact Fourier acceleration (EFA) the HMC becomes equivalent to direct…
We develop a GPU-accelerated hybrid quantum Monte Carlo (QMC) algorithm to solve the fundamental yet difficult problem of $U(1)$ gauge field coupled to fermions, which gives rise to a $U(1)$ Dirac spin liquid state under the description of…
For thermodynamics studies it is desirable to simulate two degenerate flavors and retain at least a remnant of the chiral symmetry. Staggered fermions can achieve this at the cost of rooting the determinant. Rooting can be avoided using…
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…
We compare the Hybrid Monte Carlo (HMC) and the Kramers equation algorithms for simulations of QCD with two flavors of dynamical Wilson fermions and gauge group $SU(2)$. The results for the performance of both algorithms are obtained on…
Tuning the durations of the Hamiltonian flow in Hamiltonian Monte Carlo (also called Hybrid Monte Carlo) (HMC) involves a tradeoff between computational cost and sampling quality, which is typically challenging to resolve in a satisfactory…
We present a polynomial Hybrid Monte Carlo (PHMC) algorithm as an exact simulation algorithm with dynamical Kogut-Susskind fermions. The algorithm uses a Hermitian polynomial approximation for the fractional power of the KS fermion matrix.…
We report on our implementation of the RHMC algorithm for the simulation of lattice QCD with two staggered flavors on Graphics Processing Units, using the NVIDIA CUDA programming language. The main feature of our code is that the GPU is not…
We simulate two dynamical, mass degenerate light quarks on 16^3x32 lattices with a spatial extent of 2.4 fm using the Chirally Improved Dirac operator. The simulation method, the implementation of the action and signals of equilibration are…
We discuss the construction of a chiral random matrix model for staggered fermions. This model includes $O(a^2)$ corrections to the continuum limit of staggered fermions and is related to the zero momentum limit of the Lee-Sharpe Lagrangian…
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…
Mass preconditioned HMC and DD-HMC are among the most popular algorithms to simulate Wilson fermions. We present a comparison of the performance of the two algorithms for realistic quark masses and lattice sizes. In particular, we use the…
We present a set of related Hybrid Monte Carlo methods to simulate an arbitrary number of dynamical overlap fermions. Each fermion is represented by a chiral pseudo-fermion field. The new algorithm reduces critical slowing down in the…
We continue the investigation on the applications of the Kramers equation to the numerical simulation of field theoretic models. In a previous paper we have described the theory and proposed various algorithms. Here, we compare the simplest…
Complete spectra of the staggered Dirac operator $\Dirac$ are determined in quenched four-dimensional $SU(2)$ gauge fields, and also in the presence of dynamical fermions. Periodic as well as antiperiodic boundary conditions are used. An…
The dual-fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations,…
The main purpose of this paper is to facilitate the communication between the Analytic, Probabilistic and Algorithmic communities. We present a proof of convergence of the Hamiltonian (Hybrid) Monte Carlo algorithm from the point of view of…
We present an exact pseudofermion action for hybrid Monte Carlo simulation (HMC) of one-flavor domain-wall fermion (DWF), with the effective 4-dimensional Dirac operator equal to the optimal rational approximation of the overlap-Dirac…