Related papers: An iterative method to compute the overlap Dirac o…
In this paper we show how to construct a Dirac operator on a lattice in complete analogy with the continuum. In fact we consider a more general problem, that is, the Dirac operator over an abelian finite group (for which a lattice is a…
This work is concerned with the computation of the action of a matrix function f(A), such as the matrix exponential or the matrix square root, on a vector b. For a general matrix A, this can be done by computing the compression of A onto a…
Ackermann's function can be expressed using an iterative algorithm, which essentially takes the form of a term rewriting system. Although the termination of this algorithm is far from obvious, its equivalence to the traditional recursive…
It is well known that the block Krylov subspace solvers work efficiently for some cases of the solution of differential equations with multiple right-hand sides. In lattice QCD calculation of physical quantities on a given configuration…
We use low lying eigenvectors of the overlap-Dirac operator as a probe of the QCD vacuum. If instantons play a significant role one would expect the low lying eigenmodes of the overlap-Dirac operator to consist mainly of the mixed ``would…
We propose new techniques to implement numerically the overlap-Dirac operator which exploit the physical properties of the underlying theory to avoid nested algorithms. We test these procedures in the two-dimensional Schwinger model and 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…
We propose a fast collocation method based on Krylov subspace iterative solver on general nonuniform grids for the fractional Laplacian problem, in which the fractional operator is presented in a singular integral formulation. The method is…
We discuss the eigenvalue distribution of the overlap Dirac operator in quenched QCD on lattices of size 8^{4}, 10^{4} and 12^{4} at \beta = 5.85 and \beta = 6. We distinguish the topological sectors and study the distributions of the…
New exact upper and lower bounds are derived on the spectrum of the square of the hermitian Wilson Dirac operator. It is hoped that the derivations and the results will be of help in the search for ways to reduce the cost of simulations…
We study the eigenvalues of Dirac operators in QCD with two mass degenerate dynamical fermions. The gauge configurations have been obtained with HMC and the so-called Chirally Improved fermionic action. We compare eigenvalues obtained for…
Representing massless Dirac fermions on a spatial lattice poses a potential challenge known as the Fermion Doubling problem. Addition of a quadratic term to the Dirac Hamiltonian circumvents this problem. We show that the modified…
We study the behavior of the spectrum of the Dirac operator on degenerating families of compact Riemannian surfaces, when the length $t$ of a simple closed geodesic shrinks to zero, under the hypothesis that the spin structure along the…
A new O(N) algorithm based on a recursion method, in which the computational effort is proportional to the number of atoms N, is presented for calculating the inverse of an overlap matrix which is needed in electronic structure calculations…
The authors propose a recycling Krylov subspace method for the solution of a sequence of self-adjoint linear systems. Such problems appear, for example, in the Newton process for solving nonlinear equations. Ritz vectors are automatically…
The low energy eigenmodes of the continuum QCD Dirac operator are extended, but on the lattice, due to discretization effects, the Dirac operator can have localized eigenmodes. These non-physical modes can introduce strong lattice artifacts…
The overlap Dirac operator, which satisfies the Ginsparg-Wilson relation, realizes exact chiral symmetry on the lattice without any unphysical doubler modes. To perform the path integrals, one should, however, note that the overlap fermion…
In this work, we consider Dirac-type operators with a constant delay less than two-fifths of the interval and not less than one-third of the interval. For our considered Dirac-type operators, an incomplete inverse spectral problem is…
A wide variety of (fixed-point) iterative methods for the solution of nonlinear equations (in Hilbert spaces) exists. In many cases, such schemes can be interpreted as iterative local linearization methods, which, as will be shown, can be…
We compare the low eigenvalue spectra of the Overlap Dirac operator on two sets of configurations at $\mu_I/\mu_I^c$ = 0.5 and 1.5 generated with dynamical staggered fermions at these isospin chemical potential on $24^3 \times 6$ lattices.…