We calculate the spectral function of the QCD Dirac operator using the four-dimensional effective operator constructed from the Mobius domain-wall implementation. We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number. The spectrum in the entire eigenvalue range is obtained with a single set of measurements. Results on 2+1-flavor ensembles with Mobius domain-wall sea quarks at lattice spacing ~ 0.08 fm are shown.
We calculate the charmonium correlators on the lattice with $n_f = 2+ 1$ Moebius domain wall fermion, and extract the charm quark mass and the strong coupling constant. Time moments are defined by current-current correlators, which have been calculated in the continuum theory by perturbation theory. We extract the charm quark mass by matching the lattice results with the corresponding perturbative QCD calculations, using the recently generated ensembles by the JLQCD collaboration at lattice spacings $a = 0.083, 0.055$, and $0.044$ fm.
We compute the topological susceptibility $chi_t$ of lattice QCD with $2+1$ dynamical quark flavors described by the Mobius domain wall fermion. Violation of chiral symmetry as measured by the residual mass is kept at $sim$1 MeV or smaller. We measure the fluctuation of the topological charge density in a `slab sub-volume of the simulated lattice using the method proposed by Bietenholz {it et al.} The quark mass dependence of $chi_t$ is consistent with the prediction of chiral perturbation theory, from which the chiral condensate is extracted as $Sigma^{overline{rm MS}} (mbox{2GeV}) = [274(13)(29)mbox{MeV}]^3$, where the first error is statistical and the second one is systematic. Combining the results for the pion mass $M_pi$ and decay constant $F_pi$, we obtain $chi_t = 0.229(03)(13)M_pi^2F_pi^2$ at the physical point.
We investigate a chiral property of the domain-wall fermion (DWF) system using the four-dimensional hermitian Wilson-Dirac operator $H_W$. A formula expressing the Ward-Takahashi identity quark mass $m_{5q}$ with eigenvalues of this operator is derived, which well explains the $N_5$ dependence of $m_{5q}$ observed in previous numerical simulations. We further discuss the chiral property of DWF in the large volume in terms of the spectra of $H_W$.
Chiral properties of QCD formulated with the domain-wall fermion (DWQCD) are studied using the anomalous quark mass m_{5q} and the spectrum of the 4-dimensional Wilson-Dirac operator. Numerical simulations are made with the standard plaquette gauge action and a renormalization-group improved gauge action. Results are reported on the density of zero eigenvalue obtained with the accumulation method, and a comparison is made with the results for m_{5q}.
Lattice QCD calculations including the effects of one or more non-degenerate sea quark flavors are conventionally performed using the Rational Hybrid Monte Carlo (RHMC) algorithm, which computes the square root of the determinant of $mathscr{D}^{dagger} mathscr{D}$, where $mathscr{D}$ is the Dirac operator. The special case of two degenerate quark flavors with the same mass is described directly by the determinant of $mathscr{D}^{dagger} mathscr{D}$ --- in particular, no square root is necessary --- enabling a variety of algorithmic developments, which have driven down the cost of simulating the light (up and down) quarks in the isospin-symmetric limit of equal masses. As a result, the relative cost of single quark flavors --- such as the strange or charm --- computed with RHMC has become more expensive. This problem is even more severe in the context of our measurements of the $Delta I = 1/2$ $K rightarrow pi pi$ matrix elements on lattice ensembles with $G$-parity boundary conditions, since $G$-parity is associated with a doubling of the number of quark flavors described by $mathscr{D}$, and thus RHMC is needed for the isospin-symmetric light quarks as well. In this paper we report on our implementation of the exact one flavor algorithm (EOFA) introduced by the TWQCD collaboration for simulations including single flavors of domain wall quarks. We have developed a new preconditioner for the EOFA Dirac equation, which both reduces the cost of solving the Dirac equation and allows us to re-use the bulk of our existing high-performance code. Coupling these improvements with careful tuning of our integrator, the time per accepted trajectory in the production of our 2+1 flavor $G$-parity ensembles with physical pion and kaon masses has been decreased by a factor of 4.2.