Related papers: Numerical Methods for the QCD Overlap Operator:III…
We investigate a number of algorithms that calculate the quark propagators for the overlap-Dirac fermion operator. The QCD simulations were performed at beta = 5.9 with a lattice volume of 16**3*32.
We review the numerical analysis' understanding of Krylov subspace methods for solving (non-hermitian) systems of equations and discuss its implications for lattice gauge theory computations using the example of the Wilson fermion matrix.…
A parallelizable SSOR preconditioning scheme for Krylov subspace iterative solvers in lattice QCD applications involving Wilson fermions is presented. In actual Hybrid Monte Carlo and quark propagator calculations it helps to reduce the…
The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace…
We present a method for the numerical calculation of derivatives of functions of general complex matrices. The method can be used in combination with any algorithm that evaluates or approximates the desired matrix function, in particular…
Improvements of various methods to compute the sign function of the hermitian Wilson-Dirac matrix within the overlap operator are presented. An optimal partial fraction expansion (PFE) based on a theorem of Zolotarev is given. Benchmarks…
We compute fermionic observables relevant to the study of chiral symmetry in quenched QCD using the Overlap-Dirac operator for a wide range of the fermion mass. We use analytical results to disentangle the contribution from exact zero modes…
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…
We derive an augmented Krylov subspace method with subspace recycling for computing a sequence of matrix function applications on a set of vectors. The matrix is either fixed or changes as the sequence progresses. We assume consecutive…
Rational approximations of the matrix sign function lead to multishift methods. For non-Hermitian matrices long recurrences can cause storage problems, which can be circumvented with restarts. Together with deflation we obtain efficient…
A recently developed numerical method for the calculation of derivatives of functions of general complex matrices, which can also be combined with implicit matrix function approximations such as Krylov-Ritz type algorithms, is presented. An…
In order to develop fast inversion algorithms we have used overlap solvers in two dimensions. Lattice QED theory with U(1) group symmetry in two dimensional space-times dimensions has always been a testing ground for algorithms. By the…
We present numerical simulation of QCD with two dynamical quark flavors described by the overlap fermion action on a $16^3\times 32\times (0.12 {\rm fm})^4$ lattice. We calculate pseudo-scalar masses and decay constants and investigate…
We briefly review the overlap formalism for chiral gauge theories, the overlap Dirac operator for massless fermions and its connection to domain wall fermions. We describe properties of the overlap Dirac operator, and methods to implement…
I show how to avoid a two level nested conjugate gradient procedure in the context of Hybrid Monte Carlo with the overlap fermionic action. The resulting procedure is quite similar to Hybrid Monte Carlo with domain wall fermions, but is…
We present an update of the light meson spectrum with $N_f$=2+1 overlap fermions on a $16^3\times 48$ lattice at five different up and down quark masses and two strange quark masses. Based on our experience with the previous simulation with…
We present a modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD. A larger number of test vectors than that used in conventional multigrid is generated by the smoother. This set of test…
We calculate the scalar and tensor charges of the nucleon in 2+1-flavor lattice QCD, for which the systematics of the renormalization of the disconnected diagram is well controlled. Numerical simulations are performed at a single lattice…
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.
We discuss properties of thermal Quantum Chromodynamics obtained by means of lattice simulations with overlap fermions. This fermion discretisation preserves chiral symmetry at finite lattice spacing. We present details of the formulation…