Related papers: Speeding up the HMC algorithm: Some new results
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…
We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a…
We investigate the performance of the hybrid Monte Carlo algorithm, the standard algorithm used for lattice QCD simulations involving fermions, in updating non-trivial global topological structures. We find that the hybrid Monte Carlo…
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 study the performance of an automated hybrid Monte Carlo (HMC) approach for conditional simulation of a recently proposed, single-parameter Gibbs Markov random field (Gibbs MRF). The MRF is based on a modified version of the planar…
We discuss how the peculiar properties of maximally twisted Wilson fermions can be exploited to set up a consistent LQCD computational scheme in which the CP-conserving matrix elements of the $\Delta S =1,2$ effective weak Hamiltonian can…
In this work, we introduce three algorithmic improvements to reduce the cost and improve the scaling of orbital space variational Monte Carlo (VMC). First, we show that by appropriately screening the one- and two-electron integrals of the…
Self-learning Monte Carlo method (SLMC), using a trained effective model to guide Monte Carlo sampling processes, is a powerful general-purpose numerical method recently introduced to speed up simulations in (quantum) many-body systems. In…
Perturbative expansions of several small Wilson loops are computed through next-to-next-to-leading order in unquenched lattice QCD, from Monte Carlo simulations at weak couplings. This approach provides a much simpler alternative to…
We study the efficiency of an improved Multiboson algorithm with two flavours of Wilson fermions in a realistic physical situation ($\beta = 5.60$, $\kappa = 0.156$ on a $16^3 \times 24$ lattice). The performance of this exact algorithm is…
Recently, the Fermilab heavy-quark action was extended to include dimension-six and -seven operators in order to reduce the discretization errors. In this talk, we present results of the first numerical simulations with this action (the OK…
The use of mass preconditioning or Hasenbusch filtering in modern Hybrid Monte Carlo simulations is common. At light quark masses, multiple filters (three or more) are typically used to reduce the cost of generating dynamical gauge fields;…
We address a long standing issue and determine the decorrelation efficiency of the Hybrid Monte Carlo algorithm (HMC), for full QCD with Wilson fermions, with respect to vacuum topology. On the basis of five state-of-the art QCD vacuum…
We discuss our implementation of dynamical Ginsparg-Wilson type fermions using a stout-smeared chirally improved Dirac operator. Such operators have been studied extensively in quenched calculations within the Bern-Graz-Regensburg (BGR)…
I derive the equation of motion in molecular dynamics for doing full lattice QCD simulations with clover quarks. The even-odd preconditioning technique, expected to significantly reduce the computational effort, is further developed for the…
In this paper we study asymptotic properties of different data-augmentation-type Markov chain Monte Carlo algorithms sampling from mixture models comprising discrete as well as continuous random variables. Of particular interest to us is…
We present an exact local bosonic algorithm for the simulation of dynamical fermions in lattice QCD. It is based on a non-hermitian polynomial approximation of the inverse of the quark matrix and a global Metropolis accept/reject correction…
I describe a generalization of the hybrid Monte Carlo (HMC) algorithm in which the molecular dynamics (MD) steps utilize Nambu generalized Hamiltonian dynamics. Characterized by multiple Hamiltonian functions, this formalism allows me to…
We report on the determination of the improvement coefficient $c_{\rm A}$ for the non-singlet axial-vector current $A^{a}_{\mu}(x)$ in the framework of stabilised Wilson-Clover fermions within three flavour lattice QCD. This is done by…
We discuss a mixed-action approach in which sea quarks are regularised using non-perturbatively ${\rm O}(a)$ improved Wilson fermions, while a fully-twisted tmQCD action is used for valence quarks. In this setup, automatic ${\rm O}(a)$…