Related papers: The PHMC algorithm for simulations of dynamical fe…
We present a hybrid Path Integral Monte Carlo (hPIMC) algorithm to calculate real-time quantum thermal correlation functions and demonstrate its application to open quantum systems. The hPIMC algorithm leverages the successes of classical…
For Hermitian positive definite linear systems and eigenvalue problems, the eigCG algorithm is a memory efficient algorithm that solves the linear system and simultaneously computes some of its eigenvalues. The algorithm is based on the…
We consider how to accelerate fermionic molecular dynamics algorithms by introducing n pseudofermion fields coupled with the nth root of the fermionic kernel. This reduces the maximum pseudofermionic force, and thus allows a larger…
We describe a random matrix model suitable for the simulation of the eigenvalues of the Dirac operator on the lattice for Wilson fermions. We compare the obtained global eigenvalue spectrum for various values of the hopping parameter \kappa…
We apply the Hybrid Monte Carlo method to the simulation of overlap fermions. We give the fermionic force for the molecular dynamics update. We present early results on a small dynamical chiral ensemble.
The Hypercube operator determines a variant of the approximate, truncated perfect fermion action. In this pilot study we are going to report on first experiences in dynamical QCD simulations with the Hypercube fermions. We apply this…
The extreme computational costs of calculating the sign of the Wilson matrix within the overlap operator have so far prevented four dimensional dynamical overlap simulations on realistic lattice sizes, because the computational power…
We present a quantum algorithm for the simulation of molecular systems that is asymptotically more efficient than all previous algorithms in the literature in terms of the main problem parameters. As in previous work [Babbush et al., New…
We report on our progress in using the overlap-Dirac fermion operator in simulations of lattice QCD. We have investigated the Lanczos based method of Borici, as well as various rational approximations, to calculate the step function in the…
In quantum chromodynamics (QCD) at nonzero chemical potential, the eigenvalues of the Dirac operator are scattered in the complex plane. Can the fluctuation properties of the Dirac spectrum be described by universal predictions of…
We define a sparse hermitian lattice Dirac matrix, $H$, coupling $2n+1$ Dirac fermions. When $2n$ fermions are integrated out the induced action for the last fermion is a rational approximation to the hermitian overlap Dirac operator. We…
A numerical simulation algorithm for lattice QCD is described, in which the short- and long-distance effects of the sea quarks are treated separately. The algorithm can be regarded, to some extent, as an implementation at the quantum level…
Previous investigations have shown that the canonical approach to simulating QCD at finite density is promising. The algorithm we used in our earlier work employs an exact calculation of the fermionic determinant which limits the size of…
We show that quantum number preserving Ans\"{a}tze for variational optimization in quantum chemistry find an elegant mapping to ultracold fermions in optical superlattices. Using native Hubbard dynamics, trial ground states of molecular…
I explore computer simulations of the dynamics of small multi-fermion lattice systems. The method is more general, but I concentrate on Hubbard type models where the fermions hop between a small number of connected sites. I use the natural…
In this work, we introduce a simple modification of the Monte Carlo algorithm, which we call step Monte Carlo (sMC). The sMC approach allows to simulate processes far from equilibrium and obtain information about the dynamic properties of…
The derivation of the Feynman rules for lattice perturbation theory from actions and operators is complicated, especially for highly improved actions such as HISQ. This task is, however, both important and particularly suitable for…
Dynamics of a charged particle in the canonical coordinates is a Hamiltonian system, and the well-known symplectic algorithm has been regarded as the de facto method for numerical integration of Hamiltonian systems due to its long-term…
Chiral perturbation theory for eigenvalue distributions, and equivalently random matrix theory, has recently been extended to include lattice effects for Wilson fermions. We test the predictions by comparison to eigenvalue distributions of…
The linked cell list algorithm is an essential part of molecular simulation software, both molecular dynamics and Monte Carlo. Though it scales linearly with the number of particles, there has been a constant interest in increasing its…