Related papers: Algorithm Shootout: R versus RHMC
We discuss the lattice formulation of gauge theories with fermions in arbitrary representations of the color group, and present in detail the implementation of the HMC/RHMC algorithm for simulating dynamical fermions. We discuss the…
In QCD chiral symmetry is explicitly broken by quark masses, the effect of which can be described reliably by chiral perturbation theory. Effects of explicit chiral symmetry breaking by the lattice regularisation of the Dirac operator,…
This talk presents results of a comparitive study of iterative algorithms like minimal residue ($MR$) and conjugate gradient ($CG$, $BiCG\gamma_5$, and \bicgstab) used for inverting the Dirac matrix $M$. The tests were done on the…
We examine a (3+1)-dimensional model of staggered lattice fermions with a four-fermion interaction and Z(2) chiral symmetry using the Hamiltonian formulation. This model cannot be simulated with standard fermion algorithms because those…
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…
There has been much recent progress in the understanding and reduction of the computational cost of the Hybrid Monte Carlo algorithm for Lattice QCD as the quark mass parameter is reduced. In this letter we present a new solution to this…
We test iterative algorithms, MR, QMR$\gamma_5$ and BiCG$\gamma_5$, to compare their efficiency in matrix inversion with multi-quarks (shifted matrices) within one iteration process. Our results on the 8^3 x 12 and 16^3 x 24 show that MR…
A cluster algorithm is constructed and applied to study the chiral limit of the strongly coupled lattice Schwinger model involving staggered fermions. The algorithm is based on a novel loop representation of the model. Finite size scaling…
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…
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…
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 develop an improved lattice action for heavy quarks based on Brillouin-type fermions, that have excellent energy-momentum dispersion relation. The leading discretization errors of $O(a)$ and $O(a^2)$ are eliminated at tree-level. We…
We investigate the validity of the square rooting procedure of the staggered determinant in the context of the Schwinger model. We find some evidence that at fixed physical quark mass the square root of the staggered determinant becomes…
We discuss the lattice formulation of gauge theories with fermions in arbitrary representations of the color group, and present the implementation of the RHMC algorithm for simulating dynamical Wilson fermions. A first dataset is presented…
Astrochemical simulations are a powerful tool for revealing chemical evolution in the interstellar medium. Astrochemical calculations require efficient processing of large matrices for the chemical networks. The large chemical reaction…
Some results of test runs on a $6^3\times 12$ lattice with Wilson quarks and gauge group SU(2) for a previously proposed fermion algorithm by A. Slavnov are presented.
The development of improved algorithms for QCD on the lattice has enabled us to do calculations at small quark masses and get control over the chiral extrapolation. Also finer lattices have become possible, however, a severe slowing down…
Ridge regression (RR) is an important machine learning technique which introduces a regularization hyperparameter $\alpha$ to ordinary multiple linear regression for analyzing data suffering from multicollinearity. In this paper, we present…
This paper presents a formulation of lattice fermions applicable to all quark masses, large and small. We incorporate interactions from previous light-fermion and heavy-fermion methods, and thus ensure a smooth connection to these limiting…
Models of new physics induce K_0-K_0bar mixing through operators having Dirac structures other than the "left-left" form of the Standard Model. To carry out the chiral-continuum extrapolation of results from numerical simulations, one needs…