Related papers: The role of multigrid algorithms for LQCD
We apply the topology conserving gauge action proposed by Luescher to the four-dimensional lattice QCD simulation in the quenched approximation. With this gauge action the topological charge is stabilized along the hybrid Monte Carlo…
Rather than regarding the restriction of current lattice QCD simulations to quark masses that are 5--10 times larger than those observed, we note that this presents a wonderful opportunity to deepen our understanding of QCD. Just as it has…
Random matrix theory and chiral Lagrangians offer a convenient tool for the exact calculation of microscopic spectral correlators of the Dirac operator in a well-defined finite-volume scaling regime.
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…
The QCD partition function for the Wilson Dirac operator, $D_W$, at nonzero lattice spacing $a$ can be expressed in terms of a chiral Lagrangian as a systematic expansion in the quark mass, the momentum and $a^2$. Starting from this chiral…
Some algorithmic details of our $N_f=2+1$ QCD mixed action simulations with overlap valence and improved Wilson sea quarks are presented.
The standard QCD action is improved by the addition of irrelevant operators built with chiral composites. An effective Lagrangian is derived in terms of auxiliary fields, which has the form of the phenomenological chiral Lagrangians. Our…
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…
The concept of a Dirac algebroid, which is a linear almost Dirac structure on a vector bundle, was designed to generate phase equations for mechanical systems with linear nonholonomic constraints. We apply it to systems with magnetic-like…
Many iterative parallel-in-time algorithms have been shown to be highly efficient for diffusion-dominated partial differential equations (PDEs), but are inefficient or even divergent when applied to advection-dominated PDEs. We consider the…
I make a number of observations concerning discretization errors in twisted-mass lattice QCD that can be deduced by applying chiral perturbation theory including lattice artifacts. (1) The line along which the PCAC quark mass vanishes in…
Lattice QCD simulations are computationally expensive, with the solution of the Dirac equation being the major computational bottleneck of many calculations. We introduce a novel gauge-equivariant neural-network architecture for…
The microscopic spectrum of the QCD Dirac operator is shown to obey random matrix model statistics in the bulk region of the spectrum close to the origin using finite-volume partition functions.
Graphical Processing Units (GPUs) are more and more frequently used for lattice QCD calculations. Lattice studies often require computing the quark propagators for several masses. These systems can be solved using multi-shift inverters but…
Multigrid methods were invented for the solution of discretized partial differential equations in ordered systems. The slowness of traditional algorithms is overcome by updates on various length scales. In this article we discuss…
We present a parameter-free Wilson-type lattice Dirac operator with an 81-point stencil for the covariant derivative and the Laplacian which attempts to minimize the breaking of rotational symmetry near the boundary of the Brillouin zone.…
Some exact relations for the spectral density $\rho(\lambda)$ of the Euclidean Dirac operator in $QCD$ are derived. They follow directly from the chiral symmetry of the $QCD$ lagrangian with massless quarks. New results are obtained both in…
We compute the chiral condensate in 2+1-flavor QCD through the spectrum of low-lying eigenmodes of Dirac operator. The number of eigenvalues of the Dirac operator is evaluated using a stochastic method with an eigenvalue filtering technique…
We present a progress report on a new class of multigrid solver algorithm suitable for the solution of 5d chiral fermions such as Domain Wall fermions and the Continued Fraction overlap. Unlike HDCG \cite{Boyle:2014rwa}, the algorithm works…
We describe recent results obtained as part of the High Energy Monte Carlo Grand Challenge (HEMCGC) project concerning the behaviour of lattice QCD with light dynamical Wilson quarks. We show that it is possible to reach regions of…