Related papers: Decorrelating Topology with HMC
We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong…
Traditional Markov Chain Monte Carlo methods suffer from low acceptance rate, slow mixing and low efficiency in high dimensions. Hamiltonian Monte Carlo resolves this issue by avoiding the random walk. Hamiltonian Monte Carlo (HMC) is a…
We derive non-asymptotic quantitative bounds for convergence to equilibrium of the exact preconditioned Hamiltonian Monte Carlo algorithm (pHMC) on a Hilbert space. As a consequence, explicit and dimension-free bounds for pHMC applied to…
We present results from a 2+1 flavor DWF calculation at 1/a = 3 GeV and discuss strategies for similar calculations at finer lattice spacings which will target charm physics. At weak coupling the autocorrelation time of the global…
We present the preliminary tests on two modifications of the Hybrid Monte Carlo (HMC) algorithm. Both algorithms are designed to travel much farther in the Hamiltonian phase space for each trajectory and reduce the autocorrelations among…
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 describe an HMC algorithm for dynamical overlap fermions which makes use of their good chiral properties. We test the algorithm in the Schwinger model. Topological sectors are readily changed even in the massless case.
Master-field simulations offer an approach to lattice QCD in which calculations are performed on a small number of large-volume gauge-field configurations. The latter is advantageous for simulations in which the global topological charge is…
The performance of the two-step multiboson (TSMB) algorithm is investigated in comparison with the hybrid Monte Carlo (HMC) method for compact lattice QED with standard Wilson fermions both in the Coulomb and confinement phases. The…
The topological susceptibility is an important quantity in QCD, which can be computed using lattice methods. However, at a fine lattice spacing, or when using high quality chirally symmetric quarks, algorithms which proceed in small update…
An algorithm for separating the high- and low-frequency molecular dynamics modes in Hybrid Monte Carlo simulations of gauge theories with dynamical fermions is presented. The separation is based on splitting the pseudo-fermion action into…
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…
We address a long standing problem regarding topology in lattice simulations of QCD with unimproved Wilson fermions. Earlier attempt with unimproved Wilson fermions at \beta =5.6 to verify the suppression of topological susceptibility with…
Global topological charge decorrelates very slowly or even freezes in fine lattice simulations. On the other hand, its local fluctuations are expected to survive and lead to the correct physical results as long as the volume is large…
We compare the Hybrid Monte Carlo (HMC) and the Kramers equation algorithms for simulations of QCD with two flavors of dynamical Wilson fermions and gauge group $SU(2)$. The results for the performance of both algorithms are obtained on…
We study the performance of QCD simulations with dynamical Wilson fermions by combining the Hybrid Monte Carlo algorithm with parallel tempering on $10^4$ and $12^4$ lattices. In order to compare tempered with standard simulations,…
We present a variant of the HMC algorithm with mass preconditioning (Hasenbusch acceleration) and multiple time scale integration. We have tested this variant for standard Wilson fermions at beta=5.6 and at pion masses ranging from 380 MeV…
The past few years have seen considerable progress in algorithmic development for the generation of gauge fields including the effects of dynamical fermions. The Rational Hybrid Monte Carlo (RHMC) algorithm, where Hybrid Monte Carlo is…
Hamiltonian Monte Carlo (HMC) is a Markov chain algorithm for sampling from a high-dimensional distribution with density $e^{-f(x)}$, given access to the gradient of $f$. A particular case of interest is that of a $d$-dimensional Gaussian…
The topological susceptibility is one of the few physical quantities that directly measure the properties of the QCD vacuum. Chiral perturbation theory predicts that in the small quark mass limit the topological susceptibility depends…