Related papers: Algorithms for Domain Wall Fermions
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…
Normalized hypercubic smearing improves the behavior of dynamical Wilson-clover fermions, but has the unwanted side effect that it can occasionally produce spikes in the fermion force. These spikes originate in the chain rule connecting the…
Close to the chiral limit, many calculations in numerical lattice QCD can potentially be accelerated using low-mode deflation techniques. In this paper it is shown that the recently introduced domain-decomposed deflation subspaces can be…
Hamiltonian Monte Carlo is a widely used algorithm for sampling from posterior distributions of complex Bayesian models. It can efficiently explore high-dimensional parameter spaces guided by simulated Hamiltonian flows. However, the…
We compute non-perturbatively the renormalization coefficients of scalar and pseudoscalar operators, local vector and axial currents, conserved vector and axial currents, and $O^{\Delta S=2}_{LL}$ over a wide range of energy scales using a…
We present a universal quantum Monte Carlo algorithm for simulating arbitrary high-spin (spin greater than 1/2) Hamiltonians, based on the recently developed permutation matrix representation (PMR) framework. Our approach extends a…
Hamiltonian Monte Carlo (HMC) algorithms which combine numerical approximation of Hamiltonian dynamics on finite intervals with stochastic refreshment and Metropolis correction are popular sampling schemes, but it is known that they may…
We construct a new holographic description of QCD using domain wall fermions. The construction consists of probe D7 branes in a D5 brane geometry describing quarks on a 4+1d defect in a 5+1d gauge theory. We then compactify one dimension of…
We present the Adaptive Aggregation-based Domain Decomposition Multigrid method extended to the twisted mass fermion discretization action. We show comparisons of results as a function of tuning the parameters that enter the twisted mass…
Variational optimization of neural-network representations of quantum states has been successfully applied to solve interacting fermionic problems. Despite rapid developments, significant scalability challenges arise when considering…
We present a variant of the HMC algorithm with mass preconditioning (Hasenbusch acceleration) and multiple time scale integration. We have tested this variant for standard Wilson fermions at beta=5.6 and at pion masses ranging from 380 MeV…
We introduce a variant of the Hybrid Monte Carlo (HMC) algorithm to address large-deviation statistics in stochastic hydrodynamics. Based on the path-integral approach to stochastic (partial) differential equations, our HMC algorithm…
In the study of QCD dynamics, C* boundary conditions are physically relevant in certain cases. In this paper we study the implementation of these boundary conditions in the lattice formulation of full QCD with staggered fermions. In…
We derive the exactly conserved vector, and almost conserved axial currents for rational approximations to the overlap operator with a general Mobius kernel. The approach maintains manifest Hermiticity, and allows matrix elements of the…
In recent years, the variational Monte Carlo (VMC) calculations of projected entangled pair states (PEPS) has emerged as a competitive method for computing the ground states of many-body quantum systems. This method is particularly…
We construct domain wall fermions with a staggered kernel and investigate their spectral and chiral properties numerically in the Schwinger model. In some relevant cases we see an improvement of chirality by more than an order of magnitude…
I describe a method that places the fermion fields and the gauge fields on different lattice spacings during the Hybrid Monte Carlo generation of Ginsparg-Wilson dynamical ensembles. The idea is motivated by Wilson's formulation of the…
We present a general technique for addressing sign problems that arise in Monte Carlo simulations of field theories. This method deforms the domain of the path integral to a manifold in complex field space that maximizes the average sign…
Over the past years GPUs have been successfully applied to the task of inverting the fermion matrix in lattice QCD calculations. Even strong scaling to capability-level supercomputers, corresponding to O(100) GPUs or more has been achieved.…
For big data analysis, high computational cost for Bayesian methods often limits their applications in practice. In recent years, there have been many attempts to improve computational efficiency of Bayesian inference. Here we propose an…