Related papers: The PHMC algorithm for simulations of dynamical fe…

Some algorithms for the numerically exact treatment of fermion determinants are summarised. This is not supposed to be a review, rather a concise handbook. The audience is expected to have a basic understanding of how to put fermions on a…

This paper describes a state-of-the-art parallel Lattice QCD Monte Carlo code for staggered fermions, purposely designed to be portable across different computer architectures, including GPUs and commodity CPUs. Portability is achieved…

Hamiltonian Monte Carlo (HMC) is a popular Markov chain Monte Carlo (MCMC) algorithm that generates proposals for a Metropolis-Hastings algorithm by simulating the dynamics of a Hamiltonian system. However, HMC is sensitive to large time…

Close to the continuum the lattice spacing affects the smallest eigenvalues of the Wilson Dirac operator in a very specific manner determined by the way in which the discretization breaks chiral symmetry. These effects can be computed…

We propose a method for multi-scale hybrid simulations of molecular dynamics (MD) and computational fluid dynamics (CFD). In the method, usual lattice-mesh based simulations are applied for CFD level, but each lattice is associated with a…

Our predictions for particle physics processes are realized in a chain of complex simulators. They allow us to generate high-fidelity simulated data, but they are not well-suited for inference on the theory parameters with observed data. We…

One of the open challenges in quantum computing is to find meaningful and practical methods to leverage quantum computation to accelerate classical machine learning workflows. A ubiquitous problem in machine learning workflows is sampling…

It is well-known that molecular dynamics integrators, which are used for lattice quantum chromodynamics (QCD), suffer from instabilities and possess a rather low order of the accuracy. Hence, it is highly desirable to construct a new class…

A numerical method is presented for reproducing fermionic quantum gas microscope experiments in equilibrium. By employing nested componentwise direct sampling of fermion pseudo-density matrices, as they arise naturally in determinantal…

The main purpose of this paper is to facilitate the communication between the Analytic, Probabilistic and Algorithmic communities. We present a proof of convergence of the Hamiltonian (Hybrid) Monte Carlo algorithm from the point of view of…

It has recently been demonstrated in quenched lattice simulations that the distribution of the low-lying eigenvalues of the QCD Dirac operator is universal and described by random-matrix theory. We present first evidence that this…

In the epsilon-regime of QCD the main features of the spectrum of the low-lying eigenvalues of the (euclidean) Dirac operator are expected to be described by a certain universality class of random matrix models. In particular, the latter…

Particle MCMC involves using a particle filter within an MCMC algorithm. For inference of a model which involves an unobserved stochastic process, the standard implementation uses the particle filter to propose new values for the stochastic…

An increasingly large number of HPC systems rely on heterogeneous architectures combining traditional multi-core CPUs with power efficient accelerators. Designing efficient applications for these systems has been troublesome in the past as…

We compare analytic predictions of non-Hermitian chiral random matrix theory with the complex Dirac operator eigenvalue spectrum of two-colour lattice gauge theory with dynamical fermions at nonzero chemical potential. The Dirac eigenvalues…

I extend to QCD an efficient method for lattice gauge theory with dynamical fermions. Once the eigenvalues of the Dirac operator and the density of states of pure gluonic configurations at a set of plaquette energies (proportional to the…

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…

Two cluster algorithms, based on constructing and flipping loops, are presented for worldline quantum Monte Carlo simulations of fermions and are tested on the one-dimensional repulsive Hubbard model. We call these algorithms the loop-flip…

We show how the integrators used for the molecular dynamics step of the Hybrid Monte Carlo algorithm can be further improved. These integrators not only approximately conserve some Hamiltonian $H$ but conserve exactly a nearby shadow…

Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not…