Related papers: The RHMC Algorithm for 2 Flavours of Dynamical Sta…
We introduce and compare three different Monte Carlo determinantal algorithms that allow one to compute dynamical quantities, such as the self-energy, of fermionic systems in their thermodynamic limit. We show that the most efficient…
Hamiltonian Monte Carlo (HMC) improves the computational efficiency of the Metropolis algorithm by reducing its random walk behavior. Riemannian Manifold HMC (RMHMC) further improves HMC's performance by exploiting the geometric properties…
We report on our study of the Riemannian Manifold HMC (RMHMC) algorithm with the mass term for the gauge momenta replaced by rational functions of the gauge covariant Laplace operator. A comparison of HMC and RMHMC on a 2+1+1 flavor…
One of the open challenges in quantum computing is to find meaningful and practical methods to leverage quantum computation to accelerate classical machine learning workflows. A ubiquitous problem in machine learning workflows is sampling…
Staggered fermions with smeared links can have greatly improved chiral properties. In a recent paper we introduced a simple and effective method to simulate four flavors of staggered smeared link fermions. In this work we extend the four…
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…
Smearing the gauge links of dynamical configurations removes small scale unphysical vacuum fluctuations and thus improves the chiral properties of lattice fermions. Recently we proposed the hypercubic smearing (HYP) that improves the flavor…
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…
We present first, exploratory results of a hybrid Monte-Carlo algorithm for dynamical, n_f=2, four-dimensional QCD with overlap fermions. As expected, the computational requirements are typically two orders of magnitude larger for the…
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…
Chiral logarithms in $m_\pi^2$ are calculated at one loop, taking into account the leading contributions to flavor symmetry breaking due to staggered fermions. I treat both the full QCD case (2+1 light dynamical flavors) and the quenched…
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…
We introduce a lattice fermion-Higgs model with one component `reduced staggered' fermions. In order to use the fermion field as efficiently as possible we couple the two {\em staggered} flavors to the O(4) Higgs field leading to a model…
We present a calculation of the charmonium spectrum with three flavors of dynamical staggered quarks from gauge configurations that were generated by the MILC collaboration. We use the Fermilab action for the valence charm quarks. Our…
We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong…
We propose a hybrid Monte Carlo (HMC) technique applicable to high-dimensional multivariate normal distributions that effectively samples along chaotic trajectories. The method is predicated on the freedom of choice of the HMC momentum…
Randomized quasi-Monte Carlo (RQMC) sampling can bring orders of magnitude reduction in variance compared to plain Monte Carlo (MC) sampling. The extent of the efficiency gain varies from problem to problem and can be hard to predict. This…
We propose improved estimators to compute the reweighting factors which are needed for lattice QCD calculations that rely on twisted-mass reweighting for the light quark contribution and the Rational Hybrid Monte Carlo (RHMC) algorithm for…
We discuss hybrid Monte Carlo algorithms for odd-flavor lattice QCD simulations. The algorithms include a polynomial approximation which enables us to simulate odd-flavor QCD in the framework of the hybrid Monte Carlo algorithm. In order to…
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…