Related papers: Deflation Methods in Fermion Inverters
The computation of eigenvalues of real symmetric tridiagonal matrices frequently proceeds by a sequence of QR steps with shifts. We introduce simple shift strategies, functions sigma satisfying natural conditions, taking each n x n matrix T…
We propose a framework to study the properties of the Lefschetz thimbles decomposition for lattice fermion models approaching the thermodynamic limit. The proposed set of algorithms includes the Schur complement solver and the exact…
Achieving an accurate description of fermionic systems typically requires considerably many more orbitals than fermions. Previous resource analyses of quantum chemistry simulation often failed to exploit this low fermionic number…
We have computed the decay constants for the $B$ and $D$ mesons, using quenched lattices at $\beta=6.3$, by interpolating between the static approximation of Eichten and the conventional (``heavy'' Wilson fermion) method. A more careful…
Tensor network methods have progressed from variational techniques based on matrix-product states able to compute properties of one-dimensional condensed-matter lattice models into methods rooted in more elaborate states such as projected…
Representing the time-evolution operator as a tensor network constitutes a key ingredient in several algorithms for studying quantum lattice systems at finite temperature or in a non-equilibrium setting. For a Hamiltonian composed of…
We investigate fermion--anti-fermion production in 1+1 dimensional QED using real-time lattice techniques. In this non-perturbative approach the full quantum dynamics of fermions is included while the gauge field dynamics can be accurately…
We have technically improved the non-perturbative renormalization method, proposed by Martinelli et al., by using quark momentum sources and sinks. Composite two-fermion operators up to three derivatives have been measured for Wilson…
The improvement of simulations of QCD with dynamical Wilson fermions by combining the Hybrid Monte Carlo algorithm with parallel tempering is studied on $10^4$ and $12^4$ lattices. As an indicator for decorrelation the topological charge is…
We develop the formalism for the evaluation of density-density correlators in lattice QCD that includes techniques for the computation of the all-to-all propagators involved. A novel technique in this context is the implementation of the…
Lattice computations in the Hamiltonian formulation have so far mainly focused on staggered fermions. In these proceedings, we study Wilson fermions in the Hamiltonian formulation and propose a new method to determine the resulting mass…
An algorithm for the numerical inversion of large matrices, the biconjugate gradient algorithm (BGA), is investigated in view of its use for Monte Carlo simulations of fermionic field theories. It is compared with the usual conjugate…
We summarize our recent investigations of lattice QCD with dynamical overlap fermions. We sketch algorithmic issues and our approach to solving them. We show our measurement of the topological susceptibility. We describe a computation of…
We report the findings of our extensive study of the spectra of flavoured mesons in lattice gauge theories with symplectic gauge group and fermion matter content treated in the quenched approximation. For the $Sp(4)$, $Sp(6)$, and $Sp(8)$…
Lattice QCD should allow a derivation of the $\Delta I=1/2$ rule from first principles, but numerical calculations to date have been plagued by a variety of problems. After a brief review of these problems, we present several new methods…
In this paper we present deflation and augmentation techniques that have been designed to accelerate the convergence of Krylov subspace methods for the solution of linear systems of equations. We review numerical approaches both for linear…
This paper proposes a novel matrix rank-one decomposition for quaternion Hermitian matrices, which admits a stronger property than the previous results in (sturm2003cones,huang2007complex,ai2011new). The enhanced property can be used to…
We study a model of quantum mechanical fermions with matrix-like index structure (with indices $N$ and $L$) and quartic interactions, recently introduced by Anninos and Silva. We compute the partition function exactly with $q$-deformed…
The fermion determinant is a highly non-local object and its logarithm is an extensive quantity. For these reasons it is widely believed that the determinant cannot be treated in acceptance steps of gauge link configurations that differ in…
Hadronic matrix elements of operators relevant to nucleon decay in grand unified theories are calculated numerically using lattice QCD. In this context, the domain-wall fermion formulation, combined with non-perturbative renormalization, is…