Related papers: Optimizing the domain wall fermion Dirac operator …
Lattice simulations of Quantum Chromodynamics (QCD) enable one to calculate the low-energy properties of the strong interaction among quarks and gluons based on the first principle. The most time-consuming part of the numerical simulations…
Critical slowing down for the Krylov Dirac solver presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. We propose a new multi-grid approach for chiral fermions, applicable to both…
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…
It is now a noticeable trend in High Performance Computing that the systems are becoming more and more heterogeneous. Compute nodes with a host CPU are being equipped with accelerators, the latter being a GPU or FPGA cards or both. In many…
We investigate the eigenvalues of nearly chiral lattice Dirac operators constructed with five-dimensional implementations. Allowing small violation of the Ginsparg-Wilson relation, the HMC simulation is made much faster while the…
In QCD chiral symmetry is explicitly broken by quark masses, the effect of which can be described reliably by chiral perturbation theory. Effects of explicit chiral symmetry breaking by the lattice regularisation of the Dirac operator,…
We perform hybrid Monte-Carlo simulation of $N_f=2+1+1 $ lattice QCD with domain-wall quarks at the physical point. The simulation is carried out on the $ L^3 \times T = 64^3 \times 64 $ lattice with lattice spacing $a \sim 0.064 $ fm ($ L…
A modern Fortran implementation of three Dirac operators (Wilson, Brillouin, Susskind) in lattice QCD is presented, based on OpenMP shared-memory parallelization and SIMD pragmas. The main idea is to apply a Dirac operator to $N_v$ vectors…
Solving discretized versions of the Dirac equation represents a large share of execution time in lattice Quantum Chromodynamics (QCD) simulations. Many high-performance computing (HPC) clusters use graphics processing units (GPUs) to offer…
The intersection of two ferromagnetic domain walls placed on the surface of topological insulators provides a one-way beam splitter for domain-wall Dirac fermions. Based on an analytic expression for a static two-soliton magnetic texture we…
Markov Chain Monte Carlo simulations of lattice Quantum Chromodynamics (QCD) are the only known tool to investigate non-perturbatively the theory of the strong interaction and are required to perform precision tests of the Standard Model of…
We introduce BriskStream, an in-memory data stream processing system (DSPSs) specifically designed for modern shared-memory multicore architectures. BriskStream's key contribution is an execution plan optimization paradigm, namely RLAS,…
Recent developments have shown that a lot can be gained for QCD simulations from GPU hardware. This can be exploited especially in the case of Ginsparg-Wilson fermions when the com putational costs are particularly high. In this work, we…
This paper reviews the most popular methods which are used in lattice QCD to compute the determinant of the lattice Dirac operator: Gaussian integral representation and noisy methods. Both of them lead naturally to matrix function problems.…
In this tutorial, we pedagogically review recent developments in the field of non-interacting fermionic phases of matter, focussing on the low energy description of higher-order topological insulators in terms of the Dirac equation. Our aim…
I review the lattice formulations of vector-like gauge theories (e.g. QCD) with domain-wall/overlap fermions, and discuss how to optimize the chiral symmetry for any finite $ N_s $ (sites in the fifth dimension). In this formulation, quark…
The spin-momentum locking at the Dirac surface state of a topological insulator (TI) offers a distinct possibility of a highly efficient charge-to-spin current (C-S) conversion compared with spin Hall effects in conventional paramagnetic…
We develop a GPU-accelerated hybrid quantum Monte Carlo (QMC) algorithm to solve the fundamental yet difficult problem of $U(1)$ gauge field coupled to fermions, which gives rise to a $U(1)$ Dirac spin liquid state under the description of…
We introduce PVSC-DTM (Parallel Vectorized Stencil Code for Dirac and Topological Materials), a library and code generator based on a domain-specific language tailored to implement the specific stencil-like algorithms that can describe…
A simple minded approach to implement three discretizations of the Dirac operator (staggered, Wilson, Brillouin) on two architectures (KNL and core i7) is presented. The idea is to use a high-level compiler along with OpenMP parallelization…