Related papers: Partial-Global Stochastic Metropolis Update for Dy…
We consider the problem of inference for nonlinear, multivariate diffusion processes, satisfying It\^o stochastic differential equations (SDEs), using data at discrete times that may be incomplete and subject to measurement error. Our…
Tempering is used to change the quark mass while remaining in equilibrium between the trajectories of a standard hybrid Monte Carlo simulation of four flavours of staggered fermions. The algorithm is faster for small enough quark masses,…
We propose an exact simulation algorithm for lattice QCD with dynamical Kogut-Susskind fermion in which the N_f-flavor fermion operator is defined as the N_f/4-th root of the Kogut-Susskind (KS) fermion operator. The algorithm is an…
We report on studies of simple matrix elements from simulations with two flavors of sea quarks, both staggered and Wilson. We show the decay constants of vector and pseudoscalar mesons. The effects of sea quarks are small. These simulations…
We propose an efficient Markov Chain Monte Carlo method for sampling equilibrium distributions for stochastic lattice models, capable of handling correctly long and short-range particle interactions. The proposed method is a Metropolis-type…
We describe a method to study the rate at which modes decorrelate in numerical simulations. We study the XY model updated with the Metropolis and Wolff dynamics respectively and compute the rate at which each eigenvector of the dynamics…
We propose a modification of the Hybrid-Monte-Carlo algorithm that allows for a larger step-size of the integration scheme at constant acceptance rate. The key ingredient is the splitting of the pseudo-fermion action into two parts. We test…
Accurately calculating the mass of flavor-singlet meson states from numerical lattice simulations is an important milestone for lattice QCD. Careful measurement of the full pseudoscalar flavor-singlet propagator is also a crucial step in…
It is well established that lattice artifacts can be suppressed substantially by the use of SU(3)-projected smeared links in the fermion action. An example is the Highly Improved Staggered Quark action where the ASQ-like effective links are…
I discuss the behaviour of algorithms for dynamical fermions as the sea-quark mass decreases. I focus on the Hybrid-Monte-Carlo (HMC) algorithm applied to two degenerate flavours of Wilson fermions. First, I briefly review the performance…
We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the…
Molecular dynamics algorithms are subject to some amount of error dependent on the size of the time step that is used. This error can be corrected by periodically updating the system with a Metropolis criteria, where the integration step is…
Based on recent work by Adams, I construct a lattice fermion operator that fully lifts the staggered flavor degeneracy. The resulting operator is of Wilson type but smaller by a factor of 4, better conditioned and contains 3 instead of 15…
We present results from simulations of two flavors of dynamical overlap fermions on 8^4 lattices at three values of the sea quark mass and a lattice spacing of about 0.16 fm. We measure the topological susceptibility and the chiral…
State-of-the-art simulations of discrete gauge theories are based on Markov chains with local changes in the field space, which however at very fine lattice spacings are notoriously difficult due to separated topological sectors of the…
We calculate the one-loop matching coefficients between continuum and lattice four-fermion operators for lattice operators constructed using staggered fermions and improved by the use of fattened links. In particular, we consider hypercubic…
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.…
In this paper, we study the asymptotic variance of sample path averages for inhomogeneous Markov chains that evolve alternatingly according to two different $\pi$-reversible Markov transition kernels $P$ and $Q$. More specifically, our main…
We propose a modification of the Hybrid-Monte-Carlo algorithm that allows for a larger step-size of the integration scheme at constant acceptance rate. The key ingredient is that the pseudo-fermion action is split into two parts. We test…
A study of principle is conducted on the inclusion of the fermionic determinant as a Metropolis acceptance correction. It is carried out in the 2-D Schwinger model to prepare later applications to the Schr"odinger functional. A mixed…