Search

In this second report on our recent numerical simulations of two-flavour QCD, we provide further technical details on the simulations and describe the methods we used to extract the meson masses and decay constants from the generated ensembles of gauge fields. Among the topics covered are the choice of the DD-HMC parameters, the issue of stability, autocorrelations and the statistical error analysis. Extensive data tables are included as well as a short discussion of the quark-mass dependence in partially quenched QCD, supplementing the physics analysis that was presented in the first paper in this series.

Recent conceptual, algorithmic and technical advances allow numerical simulations of lattice QCD with Wilson quarks to be performed at significantly smaller quark masses than was possible before. Here we report on simulations of two-flavour QCD at sea-quark masses from slightly above to approximately 1/4 of the strange-quark mass, on lattices with up to 64x32^3 points and spacings from 0.05 to 0.08 fm. Physical sea-quark effects are clearly seen on these lattices, while the lattice effects appear to be quite small, even without O(a) improvement. A striking result is that the dependence of the pion mass on the sea-quark mass is accurately described by leading-order chiral perturbation theory up to meson masses of about 500 MeV.

We study the spectral gap of the Wilson--Dirac operator in two-flavour lattice QCD as a function of the lattice spacing $a$, the space-time volume $V$ and the current-quark mass $m$. It turns out that the median of the probability distribution of the gap scales proportionally to $m$ and that its width is practically equal to $a/sqrt{V}$. In particular, numerical simulations are safe from accidental zero modes in the large-volume regime of QCD.