Related papers: The PHMC algorithm for simulations of dynamical fe…
We investigate the effects of low-lying fermion eigenmodes on the QCD partition function in the $\epsilon$-regime. The fermion determinant is approximated by a truncated product of low-lying eigenvalues of the overlap-Dirac operator. With…
In this lecture we review recent lattice QCD studies of the statistical properties of the eigenvalues of the QCD Dirac operator. We find that the fluctuations of the smallest Dirac eigenvalues are described by chiral Random Matrix Theories…
Quantum field theories underlie all of our understanding of the fundamental forces of nature. The are relatively few first principles approaches to the study of quantum field theories [such as quantum chromodynamics (QCD) relevant to the…
In a series of recent papers we have shown how the dynamical behavior of certain classical systems can be analyzed using operators evolving according to Heisenberg-like equations of motions. In particular, we have shown that raising and…
We discuss simulations with different lattice Dirac operators for sea and valence quarks. A goal of such a "mixed" action approach is to probe deeper the chiral regime of QCD by enabling simulations with light valence quarks. This is…
Geometric numerical integration has recently been exploited to design symplectic accelerated optimization algorithms by simulating the Lagrangian and Hamiltonian systems from the variational framework introduced in Wibisono et al. In this…
We propose a modification of the Hybrid-Monte-Carlo algorithm that allows for a larger step-size of the integration scheme at constant acceptance rate. The key ingredient is the splitting of the pseudo-fermion action into two parts. We test…
We introduce and compare three different Monte Carlo determinantal algorithms that allow one to compute dynamical quantities, such as the self-energy, of fermionic systems in their thermodynamic limit. We show that the most efficient…
We develop a GPU-accelerated hybrid quantum Monte Carlo (QMC) algorithm to solve the fundamental yet difficult problem of $U(1)$ gauge field coupled to fermions, which gives rise to a $U(1)$ Dirac spin liquid state under the description of…
Traditional heterogeneous parallel algorithms, designed for heterogeneous clusters of workstations, are based on the assumption that the absolute speed of the processors does not depend on the size of the computational task. This assumption…
We compute the condensate in QCD with two flavors of dynamical fermions using numerical simulation. The simulations use overlap fermions, and the condensate is extracted by fitting the distribution of low lying eigenvalues of the Dirac…
We present a simple derivation of a Feynman-Kac type formula to study fermionic systems. In this approach the real time or the imaginary time dynamics is expressed in terms of the evolution of a collection of Poisson processes. A computer…
Presented is a quantum computing representation of Dirac particle dynamics. The approach employs an operator splitting method that is an analytically closed-form product decomposition of the unitary evolution operator. This allows the Dirac…
We consider dynamically constrained Monte-Carlo dynamics and show that this leads to the generation of long ranged effective interactions. This allows us to construct a local algorithm for the simulation of charged systems without ever…
The use of APE smearing or other blocking techniques in lattice fermion actions can provide many advantages. There are many variants of these fat link actions in lattice QCD currently, such as FLIC fermions. The FLIC fermion formalism makes…
We present a new approach to simulating Hamiltonian dynamics based on implementing linear combinations of unitary operations rather than products of unitary operations. The resulting algorithm has superior performance to existing simulation…
We introduce a new algorithm which we call the {Rational Hybrid Monte Carlo} Algorithm (RHMC). This method uses a rational approximation to the fermionic kernel together with a noisy Kennedy-Kuti acceptance step to give an efficient…
We propose to stabilise HMC simulations of lattice QCD with very light Wilson quarks by splitting the quark determinant into two factors and by treating the factor that includes the contribution of the low modes of the Dirac operator as a…
In this paper, we develop the self-learning Monte-Carlo (SLMC) algorithm for non-abelian gauge theory with dynamical fermions in four dimensions to resolve the autocorrelation problem in lattice QCD. We perform simulations with the…
We have used Fortran 90 to implement lattice QCD. We have designed a set of machine independent modules that define fields (gauge, fermions, scalars, etc...) and overloaded operators for all possible operations between fields, matrices and…