Index theorem and random matrix theory for improved staggered quarks

We study various improved staggered quark Dirac operators on quenched gluon backgrounds in lattice QCD generated using a Symanzik-improved gluon action. We find a clear separation of the spectrum of eigenvalues into would-be zero modes and others. The number of would-be zero modes depends on the topological charge as expected from the Index Theorem, and their chirality expectation value is large. The remaining modes have low chirality and show clear signs of clustering into quartets and approaching the random matrix theory predictions for all topological charge sectors. We conclude that improvement of the fermionic and gauge actions moves the staggered quarks closer to the continuum limit where they respond correctly to QCD topology.