We study various improved staggered quark Dirac operators on quenched gluon backgrounds in lattice QCD. We find a clear separation of the spectrum of eigenvalues into high chirality, would-be zero modes and others, in accordance with the Index Theorem. We find the expected clustering of the non-zero modes into quartets as we approach the continuum limit. The predictions of random matrix theory for the epsilon regime are well reproduced. We conclude that improved staggered quarks near the continuum limit respond correctly to QCD topology.
We investigate numerically the spectral flow introduced by Adams for the staggered Dirac operator on realistic gauge configurations. We study both the unimproved and the HISQ Dirac operators. We compare the spectral flow index with the index obtained by identifying low-lying modes of large chirality.
The MILC collaborations simulations with improved staggered quarks are being extended with runs at a lattice spacing of 0.06 fm with quark masses down to one tenth the strange quark mass. We give a brief introduction to these new simulations and the determination of the lattice spacing. Then we combine these new runs with older results to study the masses of the nucleon and the Omega minus in the continuum and chiral limits.
It is well established that lattice artifacts can be suppressed substantially by the use of SU(3)-projected smeared links in the fermion action. An example is the Highly Improved Staggered Quark action where the ASQ-like effective links are constructed from reunitarized Fat7 links. A general procedure is presented for computing the derivative of the fermion action with respect to the base links (fermion force) - a key component in dynamical simulations using molecular dynamics evolution. The method is iterative and can be applied to actions with arbitrary levels of smearing and reunitarization. The cost of calculating the fermion force is determined for the ASQ action and the HISQ action. Test results show that calculating the HISQ force is about two times more expensive than the ASQ force.
We report on a calculation of $B_c$ ground state and radial excitation energies, obtained from heavy-charm highly improved staggered quark (HISQ) correlators computed on MILC gauge ensembles, with lattice spacings down to $a=0.044$ fm. Using HISQ valence quarks on progressively finer lattices allows us to simulate up to the $b$-quark mass. In particular we focus on the $B_c(2S)$ energy, which we compare with O(alpha_s)-improved non-relativistic QCD results computed on the same ensembles and recent experimental results from ATLAS.
We present the first computation in a program of lattice-QCD baryon physics using staggered fermions for sea and valence quarks. For this initial study, we present a calculation of the nucleon mass, obtaining $964pm16$ MeV with all sources of statistical and systematic errors controlled and accounted for. This result is the most precise determination to date of the nucleon mass from first principles. We use the highly-improved staggered quark action, which is computationally efficient. Three gluon ensembles are employed, which have approximate lattice spacings $a=0.09$ fm, $0.12$ fm, and $0.15$ fm, each with equal-mass $u$/$d$, $s$, and $c$ quarks in the sea. Further, all ensembles have the light valence and sea $u$/$d$ quarks tuned to reproduce the physical pion mass, avoiding complications from chiral extrapolations or nonunitarity. Our work opens a new avenue for precise calculations of baryon properties, which are both feasible and relevant to experiments in particle and nuclear physics.