Related papers: Polynomial Filtered HMC -- an algorithm for lattic…
I apply a recently developed algorithm for reweighting simulations of lattice QCD from one quark mass to another to simulations performed with overlap fermions in the epsilon regime. I test it by computing the condensate from distributions…
The pseudofermion action of the Hybrid Monte Carlo (HMC) algorithm for dynamical fermions is modified to directly incorporate Incomplete LU (ILU) factorisation. This reduces the stochastic noise and allows a larger molecular dynamics…
We simulate Quantum Chromodynamics (QCD) in four Euclidean dimensions with two (degenerate mass) flavors of dynamical quarks. The Dirac operator we use is the so-called chirally improved Dirac operator. We discuss the algorithm used for the…
We give details of our precise determination of the light quark masses m_{ud}=(m_u+m_d)/2 and m_s in 2+1 flavor QCD, with simulated pion masses down to 120 MeV, at five lattice spacings, and in large volumes. The details concern the action…
We discuss the lattice formulation of gauge theories with fermions in arbitrary representations of the color group, and present in detail the implementation of the HMC/RHMC algorithm for simulating dynamical fermions. We discuss the…
We report on the simulation of QCD with light dynamical quarks using the two-step multi-boson (TSMB) algorithm. In an exploratory study with two flavours of quarks at lattice spacing about 0.27 fm and with quark mass down to one sixth of…
We investigate a hybrid numerical algorithm aimed at the large-scale cosmological N-body simulation for the on-going and the future high precious sky surveys. It makes use of a truncated Fast Multiple Method (FMM) for short-range gravity,…
The numerical computation of many hadronic correlation functions is exceedingly difficult due to the exponentially decreasing signal-to-noise ratio with the distance between source and sink. Multilevel integration methods, using independent…
We consider how to accelerate fermionic molecular dynamics algorithms by introducing n pseudofermion fields coupled with the nth root of the fermionic kernel. This reduces the maximum pseudofermionic force, and thus allows a larger…
The determination of quark masses from lattice QCD simulations requires a non-perturbative renormalization procedure and subsequent scale evolution to high energies, where a conversion to the commonly used MS-bar scheme can be safely…
We present the results of an effort to accelerate a Rational Hybrid Monte Carlo (RHMC) program for lattice quantum chromodynamics (QCD) simulation for 2 flavours of staggered fermions on multiple Kepler K20X GPUs distributed on different…
In this paper, we study first-order methods on a large variety of low-rank matrix optimization problems, whose solutions only live in a low dimensional eigenspace. Traditional first-order methods depend on the eigenvalue decomposition at…
Accurate measurements of K, D and B meson mixing amplitudes provide stringent constraints in the Unitary Triangle analysis, as well as useful bounds on New Physics scales. Lattice QCD provides a non perturbative tool to compute the hadronic…
Hamiltonian Monte Carlo (HMC) is a state-of-the-art Markov chain Monte Carlo sampling algorithm for drawing samples from smooth probability densities over continuous spaces. We study the variant most widely used in practice, Metropolized…
We consider a lattice-inspired random matrix model for the QCD chiral phase transition at finite chemical potential. Useful features of the usual RMM for QCD at finite chemical potential are reobtained, some being brought closer to their…
I report the status of the dynamical overlap fermion project by the JLQCD and TWQCD collaborations. So far, the simulations have been completed with two flavors of overlap sea quarks in a wide range of sea quark mass corresponding the…
Understanding the behavior of interacting fermions is of fundamental interest in many fields ranging from condensed matter to high energy physics. Developing numerically efficient and accurate simulation methods is an indispensable part of…
We propose a simple iterative algorithm to construct the optimal multi-configuration approximation of an $N$-fermion wave function. That is, $M\geq N $ single-particle orbitals are sought iteratively so that the projection of the given wave…
We describe a new Hybrid Monte Carlo (HMC) algorithm for dynamical overlap fermions, which improves the rate of topological index changes by adding an additional (intensive) term to the action for the molecular dynamics part of the…
We introduce efficient numerical methods for generic HJM equations of interest rate theory by means of high-order weak approximation schemes. These schemes allow for QMC implementations due to the relatively low dimensional integration…