Methods for pseudoscalar flavor-singlet mesons with staggered fermions

We present the first $2+1$ flavor lattice QCD calculations of pseudoscalar flavor-singlet propagators using improved staggered fermions. We explore the relevant techniques and discuss prospects for the larger scale studies now in progress. The disconnected correlator is shown to have a highly non-Gaussian distribution and reliable estimates of the errors require care.