Related papers: A Parallel SSOR Preconditioner for Lattice QCD
We show how the hopping parameter expansion at order $\kappa^2$ and $\kappa^4$ can be exploited in the simulation of lattice QCD with two flavours of degenerate Wilson fermions. A natural extension of this idea is a "UV-filtering" by using…
The perturbative and nonperturbative renormalisation of quark-antiquark operators in lattice QCD with two flavours of clover fermions is investigated within the research programme of the QCDSF collaboration. Operators with up to three…
Large, sparse linear systems are pervasive in modern science and engineering, and Krylov subspace solvers are an established means of solving them. Yet convergence can be slow for ill-conditioned matrices, so practical deployments usually…
Stencil computations are widely used to simulate the change of state of physical systems across a multidimensional grid over multiple timesteps. The state-of-the-art techniques in this area fall into three groups: cache-aware tiled looping…
We present a method to parallelize the stochastic cutoff (SCO) method, which is a Monte-Carlo method for long-range interacting systems. After interactions are eliminated by the SCO method, we subdivide the lattice into non-interacting…
We compute the quark number susceptibilities in two flavor QCD for staggered fermions by adding the chemical potential as a Lagrange multiplier for the point-split number density term. Since lesser number of quark propagators are required…
We investigate the perturbative and nonperturbative renormalization of composite operators in lattice QCD restricting ourselves to operators that are bilinear in the quark fields (quark-antiquark operators). These include operators which…
An extension of the synchronous parallel kinetic Monte Carlo (pkMC) algorithm developed by Martinez {\it et al} [{\it J.\ Comp.\ Phys.} {\bf 227} (2008) 3804] to discrete lattices is presented. The method solves the master equation…
I review the theoretical foundations, properties as well as the simulation results obtained so far of a variant of the Wilson lattice QCD formulation: Wilson twisted mass lattice QCD. Emphasis is put on the discretization errors and on the…
A new method for computing all elements of the lattice quark propagator is proposed. The method combines the spectral decomposition of the propagator, computing the lowest eigenmodes exactly, with noisy estimators which are 'diluted', i.e.…
It is shown how on-shell O(a) improvement can be implemented non-perturbatively in lattice QCD with Wilson quarks. Improvement conditions are obtained by requiring the PCAC relation to hold exactly in certain matrix elements. These are…
A previously introduced multi-boson technique for the simulation of QCD with dynamical quarks is described and some results of first test runs on a $6^3\times12$ lattice with Wilson quarks and gauge group SU(2) are reported.
We discuss techniques for accelerating the self consistent field (SCF) iteration for solving the Kohn-Sham equations. These techniques are all based on constructing approximations to the inverse of the Jacobian associated with a fixed point…
Scalar and pseudo-scalar resonances decaying to top quarks are common predictions in several scenarios beyond the standard model (SM) and are extensively searched for by LHC experiments. Challenges on the experimental side require…
We consider the renormalisation of lattice QCD operators with one and two covariant derivatives related to the first and second moments of generalised parton distributions and meson distribution amplitudes. Employing the clover fermion…
The efficient solution of moderately large-scale linear systems arising from the KKT conditions in optimal control problems (OCPs) is a critical challenge in robotics. With the stagnation of Moore's law, there is growing interest in…
There has been much recent progress in the understanding and reduction of the computational cost of the Hybrid Monte Carlo algorithm for Lattice QCD as the quark mass parameter is reduced. In this letter we present a new solution to this…
We demonstrate that substantial progress can be achieved in the study of the phase structure of 4-dimensional compact QED by a joint use of hybrid Monte Carlo and multicanonical algorithms, through an efficient parallel implementation. This…
In this talk, we review a QCD factorization based approach to extract parton distribution and correlation functions from lattice QCD calculation of single hadron matrix elements of quark-gluon operators. We argue that although the lattice…
We develop a generalization of low-mode averaging in which the number of low quark modes of the Dirac operator required for a constant variance reduction can be kept independent of the volume by exploiting their local coherence. Typically…