Related papers: Speeding up the HMC algorithm: Some new results
We introduce a dynamical fermion algorithm which is based on the hybrid Monte Carlo (HMC) algorithm, but without pseudofermions. The molecular dynamics steps in HMC are retained except the derivatives with respect to the gauge fields are…
We discuss how the integrators used for the Hybrid Monte Carlo (HMC) algorithm not only approximately conserve some Hamiltonian $H$ but exactly conserve a nearby shadow Hamiltonian (\tilde H), and how the difference $\Delta H \equiv \tilde…
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…
We present a practical strategy to optimize a set of Hybrid Monte Carlo parameters in simulations of QCD and QCD-like theories. We specialize to the case of mass-preconditioning, with multiple time-step Omelyan integrators. Starting from…
Variational inference lies at the core of many state-of-the-art algorithms. To improve the approximation of the posterior beyond parametric families, it was proposed to include MCMC steps into the variational lower bound. In this work we…
We present results for pseudoscalar decay constants of heavy-light mesons using both quenched and N_f = 2 dynamical fermion configurations. A variety of fermion actions is investigated : Wilson, nonperturbative clover, and fat-link clover.…
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.…
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…
Considering Ginsparg-Wilson type fermions dynamically in lattice QCD simulations is a challenging task. The hope is to be able to approach smaller pion masses and to eventually reach physical situations. The price to pay is substantially…
We test the reweighting of the quark determinant of O(a) improved Wilson fermions in the domain-decomposed hybrid Monte-Carlo algorithm. Specifically, we implement a reweighting in a twisted-mass parameter proposed by Palombi and L\"uscher…
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 present simulation results for the 2-flavour Schwinger model with dynamical Ginsparg-Wilson fermions. Our Dirac operator is constructed by inserting an approximately chiral hypercube operator into the overlap formula, which yields the…
The improved Wilson quark action - the clover action - is constructed to have smaller discretization errors than the normal Wilson quark action. We test this in a quenched spectroscopy computation on 6 lattice ensembles with spacings from…
Results from two studies of full QCD with two flavours of dynamical Wilson fermions are presented. At beta=5.6, the region 0.83 > m_pi/m_rho > 0.56 at m_pia > 0.23 L^{-1} is explored. The SESAM collaboration has generated ensembles of about…
Stabilized Wilson fermions are a reformulation of Wilson clover fermions that incorporates several numerical stabilizing techniques, but also a local change of the fermion action - the original clover term being replaced with an…
We show that it is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with Wilson fermions by taking arithmetic averages of correlators computed in theories regularized…
In this paper we discuss how the peculiar properties of twisted lattice QCD at maximal twist can be employed to set up a consistent computational scheme in which, despite the explicit breaking of chiral symmetry induced by the presence of…
Simulations of odd flavors QCD can be 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…
We discuss how dynamical fermion computations may be made yet cheaper by using symplectic integrators that conserve energy much more accurately without decreasing the integration step size. We first explain why symplectic integrators…
Within the HMC algorithm, we discuss how, by using the shadow Hamiltonian and the Poisson brackets, one can achieve a simple factorization in the dependence of the Hamiltonian violations upon either the algorithmic parameters or the…