Related papers: Accelerating Dynamical Fermion Computations using …

Representing massless Dirac fermions on a spatial lattice poses a potential challenge known as the Fermion Doubling problem. Addition of a quadratic term to the Dirac Hamiltonian circumvents this problem. We show that the modified…

We suggest and implement a new Monte Carlo strategy for correlated models involving fermions strongly coupled to classical degrees of freedom, with accurate handling of quenched disorder as well. Current methods iteratively diagonalise the…

We consider the Riemann manifold Hamiltonian Monte Carlo (RMHMC) method for solving statistical inverse problems governed by partial differential equations (PDEs). The power of the RMHMC method is that it exploits the geometric structure…

An exact, nonlocal, finite step-size algorithm for Monte Carlo simulation of theories with dynamical fermions is proposed. The algorithm is based on obtaining the new configuration U' from the old one U by solving the equation $ M(U') \eta…

We show that in practical simulations of lattice QCD with two dynamical light fermion species the PHMC algorithm samples configuration space differently from the commonly used HMC algorithm.

We present an exact version of the local bosonic algorithm for the simulation of dynamical quarks in lattice QCD. This version is based on a non-hermitian polynomial approximation of the inverse of the quark matrix. A Metropolis test…

Using a dual representation of lattice fermion models that is based on spin-charge transformation and fermionisation of the original description, I derive an algorithm for diagrammatic Monte Carlo simulation of strongly correlated systems.…

We compare the performance of the Kramers Equation Monte Carlo (KMC) Algorithm with that of the Hybrid Monte Carlo (HMC) algorithm for numerical simulations with dynamical Kogut-Susskind fermions. Using the lattice Gross-Neveu model in 2…

We present a quantum algorithm to compute the logarithm of the determinant of the fermion matrix, assuming access to a classical lattice gauge field configuration. The algorithm uses the quantum eigenvalue transform, and quantum mean…

We describe how to incorporate symmetries of the Hamiltonian into auxiliary-field quantum Monte Carlo calculations (AFQMC). Focusing on the case of Abelian symmetries, we show that the computational cost of most steps of an AFQMC…

We present a method for direct hybrid Monte Carlo simulation of graphene on the hexagonal lattice. We compare the results of the simulation with exact results for a unit hexagonal cell system, where the Hamiltonian can be solved…

I present a summary of recent algorithmic developments for lattice field theories. In particular I give a pedagogical introduction to the new Multicanonical algorithm, and discuss the relation between the Hybrid Overrelaxation and Hybrid…

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,…

We introduce a lattice fermion-Higgs model with one component `reduced staggered' fermions. In order to use the fermion field as efficiently as possible we couple the two {\em staggered} flavors to the O(4) Higgs field leading to a model…

Metal hydrides are promising candidates for hydrogen storage applications. From a materials discovery perspective, an accurate, efficient computational workflow is urgently required that can rapidly analyze/predict thermodynamic properties…

A method based on the Monte Carlo inversion of the Dirac operator on the lattice provides low noise results for the correlations entering the definition of the heavy meson decay constant in the static limit. The method is complementary to…

We discuss the adaptation of the Hybrid Monte Carlo algorithm to overlap fermions. We derive a method which can be used to account for the delta function in the fermionic force caused by the differential of the sign function. We discuss the…

We enable the automatic construction of Hybrid Monte Carlo (HMC) forces in lattice gauge theory by performing reverse-mode automatic differentiation at the level of optimized LLVM intermediate representation, making the approach applicable…

We introduce a semistochastic implementation of the power method to compute, for very large matrices, the dominant eigenvalue and expectation values involving the corresponding eigenvector. The method is semistochastic in that the matrix…

We introduce a recent symplectic integration scheme derived for solving physically motivated systems with non-separable Hamiltonians. We show its relevance to Riemannian manifold Hamiltonian Monte Carlo (RMHMC) and provide an alternative to…