Related papers: Efficiency of different matrix inversion methods a…
We study the efficiency of an improved Multiboson algorithm with two flavours of Wilson fermions in a realistic physical situation ($\beta = 5.60$, $\kappa = 0.156$ on a $16^3 \times 24$ lattice). The performance of this exact algorithm is…
We consider three improvements to seed methods for Hermitian linear systems with multiple right-hand sides: only the Krylov subspace for the first system is used for seeding subsequent right-hand sides, the first right-hand side is solved…
Different formulations of the $4d$ compact lattice QED with staggered fermions (standard Wilson and modified by suppression of lattice artifacts) are investigated by Monte Carlo simulations within the quenched approximation. We show that…
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.…
When lattice QCD is formulated in sectors of fixed quark numbers, the canonical fermion determinants can be expressed explicitly in terms of transfer matrices. This in turn provides a complete factorization of the fermion determinants in…
We apply a recent proposal to speed up the Hybrid-Monte-Carlo simulation of systems with dynamical fermions to two flavour QCD with clover-improvement. The basic idea of our proposal is to split the fermion matrix into two factors with a…
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 briefly describe some of our recent results for the mass spectrum and matrix elements using $O(a)$ improved fermions for quenched QCD. Where possible a comparison is made between improved and Wilson fermions.
Lattice formulations of QCD with Wilson fermions and a chirally twisted quark mass matrix provide an attractive framework for non-perturbative numerical studies. Owing to reparameterization invariance, the limiting continuum theory is just…
We present the results of a precision computation of B_K with Wilson fermions. Simulations are performed at different lattice spacings, enabling continuum limit extrapolations. Two different twisted mass QCD (tmQCD) regularisations are…
We present a comparison of a number of iterative solvers of linear systems of equations for obtaining the fermion propagator in lattice QCD. In particular, we consider chirally invariant overlap and chirally improved Wilson (maximally)…
First results of a recently started simulation of full QCD with two flavours of sea-quarks at a coupling of $\beta = 5.6$ on a $16^3 \times 32$ lattice are presented. Emphasis is laid on the statistical significance that can be achieved by…
Cost estimates for simulations of full QCD with n_f=2 Wilson fermions by hybrid Monte Carlo are presented. The extrapolations are based on the average number of iterations of the iterative solver within the fermionic part of the HMC…
We present an investigation of gauge-invariant smearing for Wilson fermions in quenched lattice QCD on a $24^3 \times 48$ lattice at $\beta = 6.2$. We demonstrate a smearing algorithm that allows a substantial improvement in the…
Heavy-quark effective theory (HQET) is applied to lattice QCD with Wilson fermions at fixed lattice spacing a. This description is possible because heavy-quark symmetries are respected. It is desirable because the ultraviolet cutoff $1/a$…
Within the framework of quenched lattice QCD and using O(a) improved Wilson fermions and non-perturbative renormalisation, a high statistics computation of low moments of the unpolarised nucleon structure functions is given. Particular…
We present a performance comparison of the Kramers equation and the boson algorithms for simulations of QCD with two flavors of dynamical Wilson fermions and gauge group $SU(2)$. Results are obtained on $6^312$, $8^312$ and $16^4$ lattices.…
Lattice QCD calculations require significant computational effort, with the dominant fraction of resources typically spent in the numerical inversion of the Dirac operator. One of the simplest methods to solve such large and sparse linear…
We employ the Complex Langevin method for simulation of complex-valued actions. First, we show how to test for convergence of the method by explicitely computing boundary terms and demonstrate this in a model. Then we investigate the…
In this paper we discuss how the peculiar properties of twisted lattice QCD at maximal twist can be employed to set up a consistent computational scheme in which, despite the explicit breaking of chiral symmetry induced by the presence of…