Related papers: Comparing iterative methods to compute the overlap…
The action of the overlap-Dirac operator on a vector is typically implemented indirectly through a multi-shift conjugate gradient solver. The compute-time required depends upon the condition number, $\kappa$, of the matrix that is used as…
The low-lying eigenvalues of a (sparse) hermitian matrix can be computed with controlled numerical errors by a conjugate gradient (CG) method. This CG algorithm is accelerated by alternating it with exact diagonalisations in the subspace…
I apply a recently developed algorithm for reweighting simulations of lattice QCD from one quark mass to another to simulations performed with overlap fermions in the epsilon regime. I test it by computing the condensate from distributions…
The low-lying eigenvalue spectrum of the QCD Dirac operator in the epsilon-regime is expected to match with that of chiral Random Matrix Theory (ChRMT). We study this correspondence for the case including sea quarks by performing two-flavor…
In this work we introduce a memory-efficient method for computing the action of a Hermitian matrix function on a vector. Our method consists of a rational Lanczos algorithm combined with a basis compression procedure based on rational…
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…
In lattice QCD it is possible, in principle, to determine the parameters in the effective chiral lagrangian (including weak interaction couplings) by performing numerical simulations in the $\epsilon$--regime, i.e. at quark masses where the…
It was recently proposed by the second author to consider lattice formulations of QCD in which complete actions, including the gauge part, are built explicitly from a given Dirac operator D. In a simple example of such theory, the gauge…
Random Matrix Theory has been successfully applied to lattice Quantum Chromodynamics. In particular, a great deal of progress has been made on the understanding, numerically as well as analytically, of the spectral properties of the Wilson…
We analyze how individual eigenvalues of the QCD Dirac operator at nonzero quark chemical potential are distributed in the complex plane. Exact and approximate analytical results for both quenched and unquenched distributions are derived…
We compute the ratio $\Lambda_L/\Lambda_{\bar{MS}}$ between the scale parameter $\Lambda_L$, associated with a lattice formulation of QCD using the overlap-Dirac operator, and $\Lambda_{\bar{MS}}$ of the $\bar{\rm MS}$ renormalization…
Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work…
We evaluate for arbitrary even dimensions the classical continuum limit of the lattice axial anomaly defined by the overlap-Dirac operator. Our calculational scheme is simple and systematic. In particular, a powerful topological argument is…
We compute the one-loop lattice renormalization of the two-quark operators $\bar{\psi} \Gamma \psi$, where $\Gamma$ denotes the generic Dirac matrix, for the lattice formulation of QCD using the overlap-Dirac operator. We also study the…
We consider the task of computing solutions of linear systems that only differ by a shift with the identity matrix as well as linear systems with several different right hand sides. In the past Krylov subspace methods have been developed…
We describe a Lanczos-based algorithm for approximating the product of a rational matrix function with a vector. This algorithm, which we call the Lanczos method for optimal rational matrix function approximation (Lanczos-OR), returns the…
We compare eigenvalue correlations of the Dirac operator with a chemical potential obtained from lattice simulations of quenched QCD with analytic predictions obtained from chiral effective theories in the zero-momentum limit. By comparing…
The overlap hypercube fermion is constructed by inserting a lattice fermion with hypercubic couplings into the overlap formula. One obtains an exact Ginsparg-Wilson fermion, which is more complicated than the standard overlap fermion, but…
A new lattice action is proposed for the overlap Dirac matrix with nonzero chemical potential. It is shown to preserve the full chiral invariance for all values of lattice spacing exactly. It is further demonstrated to arise in the domain…
We investigate the eigenvalues of nearly chiral lattice Dirac operators constructed with five-dimensional implementations. Allowing small violation of the Ginsparg-Wilson relation, the HMC simulation is made much faster while the…