Related papers: Partial-Global Stochastic Metropolis Update for Dy…
The Metropolis-Hastings algorithm allows one to sample asymptotically from any probability distribution $\pi$. There has been recently much work devoted to the development of variants of the MH update which can handle scenarios where such…
The two-flavour Schwinger model is used to test the local boson action algorithm of M Luescher. The autocorrelation time is found to rise linearly with the number of auxiliary boson fields. An extension to the algorithm is proposed which…
In an experimental study of single enzyme reactions, it has been proposed that the rate constants of the enzymatic reactions fluctuate randomly, according to a given distribution. To quantify the uncertainty arising from random rate…
We study the finite temperature phase transition of four flavor staggered fermions with hypercubic fat link actions on N_t=4 and N_t=6 temporal lattices. Our fat links are constructed with hypercubic blocking (HYP) and therefore are very…
We propose an adaptive Metropolis-Hastings algorithm in which sampled data are used to update the proposal distribution. We use the samples found by the algorithm at a particular step to form the information-theoretically optimal mean-field…
We calculate one-loop matching factors for bilinear operators composed of improved staggered fermions. We compare the results for different improvement schemes used in the recent literature, all of which involve the use of smeared links.…
Staggered fermion shift symmetries correspond to translations of the fermion field within the unit cell of a hypercubic lattice. They satisfy an algebra and in four Euclidean dimensions can be related to a discrete subgroup of an $SU(4)$…
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…
We report on recent progress in employing the Highly Improved Staggered Quark (HISQ) action introduced by the HPQCD/UKQCD collaboration in simulations with dynamical fermions. The HISQ action is an order $a^2$ Symanzik-improved action with…
Particle Metropolis-Hastings enables Bayesian parameter inference in general nonlinear state space models (SSMs). However, in many implementations a random walk proposal is used and this can result in poor mixing if not tuned correctly…
We present a new multiple-try Metropolis-Hastings algorithm designed to be especially beneficial when a tailored proposal distribution is available. The algorithm is based on a given acyclic graph $G$, where one of the nodes in $G$, $k$…
We describe a new Hybrid Monte Carlo (HMC) algorithm for dynamical overlap fermions, which improves the rate of topological index changes by adding an additional (intensive) term to the action for the molecular dynamics part of the…
Stochastic hydrodynamics provides a dynamical framework for the evolution of fluctuations in heavy-ion collisions, but poses significant challenges in numerical simulations. We present an algorithm for the simulation of non-relativistic…
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 simulate Quantum Chromodynamics in four Euclidean dimensions with two (degenerate mass) flavors of dynamical quarks. The Dirac operator is the so-called chirally improved operator that has been studied so far in quenched calculations. We…
Staggered fermions with 4 tastes are expected to describe 4-flavor QCD in the continuum limit, therefore at finite lattice spacing the staggered determinant should be equivalent to an SU(4) flavor-symmetric system up to lattice artifacts.…
We study the impact of hypercubic (HYP) smearing on the size of taste breaking for staggered fermions, comparing to unimproved and to asqtad-improved staggered fermions. As in previous studies, we find a substantial reduction in…
The R algorithm is widely used for simulating two flavours of dynamical staggered fermions. We give a simple proof that the algorithm converges to the desired probability distribution to within O(dt^2) errors, but show that the relevant…
We propose various improvements of finite step-size updating for full QCD on the lattice that might turn finite step-size updating into a viable alternative to the hybrid Monte Carlo algorithm. These improvements are noise reduction of the…
We present a new action for highly improved staggered fermions. We show that perturbative calculations for the new action are well-behaved where those of the conventional staggered action are badly behaved. We discuss the effects of the new…