Related papers: Parallel tempering in full QCD with Wilson fermion…
The hybrid Monte Carlo (HMC) algorithm is a ubiquitous method in computational physics with applications ranging from condensed matter to lattice QCD and beyond. However, HMC simulations often suffer from long autocorrelation times,…
We construct new Ginsparg-Wilson fermions for QCD by inserting an approximately chiral Dirac operator - which involves ingredients of a perfect action - into the overlap formula. This accelerates the convergence of the overlap Dirac…
We present a massively parallel quantum Monte Carlo based implementation of real-space dynamical mean-field theory for general inhomogeneous correlated fermionic lattice systems. As a first application, we study magnetic order in a binary…
We propose a modification of the Hybrid Monte Carlo (HMC) algorithm that overcomes the topological freezing of a two-dimensional $U(1)$ gauge theory with and without fermion content. This algorithm includes reversible jumps between…
We present recent results of SESAM's large scale lattice simulation of QCD with two dynamical flavours of Wilson fermions. The emphasis is on future prospects in the extraction of flavour-singlet matrix elements, i.e the pi-nucleon sigma…
We present an exact dynamical QCD simulation algorithm for the $O(a)$-improved Wilson fermion with odd number of flavors. Our algorithm is an extension of the non-Hermitian polynomials HMC algorithm proposed by Takaishi and de Forcrand…
Current dynamical overlap fermion hybrid Monte Carlo simulations encounter large fermionic forces when there is mixing between near zero-eigenvectors of the kernel operator. This leads to low acceptance rates when there is a large density…
In this short note, we show how the parallel adaptive Wang-Landau (PAWL) algorithm of Bornn et al. (2013) can be used to automate and improve simulated tempering algorithms. While Wang-Landau and other stochastic approximation methods have…
Current trends in parallel processors call for the design of efficient massively parallel algorithms for scientific computing. Parallel algorithms for Monte Carlo simulations of thermodynamic ensembles of particles have received little…
Stabilized Wilson fermions are a reformulation of Wilson clover fermions that incorporates several numerical stabilizing techniques, but also a local change of the fermion action - the original clover term being replaced with an…
In the QCD phase diagram, the dependence of the pseudo-critical temperature, $T_{\rm{pc}}$, on the baryon chemical potential, $\mu_B$, is of fundamental interest. The variation of $T_{\rm{pc}}$ with $\mu_B$ is normally captured by $\kappa$,…
Modern problems in astronomical Bayesian inference require efficient methods for sampling from complex, high-dimensional, often multi-modal probability distributions. Most popular methods, such as Markov chain Monte Carlo sampling, perform…
We describe recent results concerning the behavior of lattice QCD with light dynamical Wilson and Staggered quarks. We show that it is possible to reach regions of parameter space with light pions $m_\pi\approx 0.2/a$ using Wilson fermions.…
We present an update on our previous determination of the Lambda parameter in QCD. The main emphasis is on results for two flavours of light dynamical quarks, where we can now almost double the amount of data used, including values at…
Multimodal structures in the sampling density (e.g. two competing phases) can be a serious problem for traditional Markov Chain Monte Carlo (MCMC), because correct sampling of the different structures can only be guaranteed for infinite…
I present one-loop perturbative calculations of matching coefficients between matrix elements in continuum regulated QCD and lattice QCD with overlap fermions, with emphasis a recently-proposed variant discretization of the overlap. These…
We propose a new sampler that integrates the protocol of parallel tempering with the Nos\'e-Hoover (NH) dynamics. The proposed method can efficiently draw representative samples from complex posterior distributions with multiple isolated…
A status report is presented on the large-volume simulations in the Schroedinger functional with two flavours of O(a) improved Wilson quarks performed by the ALPHA collaboration. The physics goal is to set the scale for the computation of…
We have developed a technique to accelerate the acquisition of effectively uncorrelated configurations for off-lattice models of dense polymer melts which makes use of both parallel tempering and large scale Monte Carlo moves. The method is…
We present an exact pseudofermion action for hybrid Monte Carlo simulation (HMC) of one-flavor domain-wall fermion (DWF), with the effective 4-dimensional Dirac operator equal to the optimal rational approximation of the overlap-Dirac…