We present a new lattice QCD analysis of heavy-quark pseudoscalar-pseudoscalar correlators, using gluon configurations from the MILC collaboration that include vacuum polarization from $u$, $d$, $s$ and $c$~quarks($n_f=4$). We extract new values for
the QCD coupling and for the $c$ quarks $overline{mathrm{MS}}$ mass: $alpha_{overline{mathrm{MS}}}(M_Z,n_f=5) = 0.11822(74)$ and $m_c(3mathrm{GeV}, n_f=4) = 0.9851(63)$GeV. These agree well with our earlier simulations using $n_f=3$ sea quarks, vindicating the perturbative treatment of $c$ quarks in that analysis. We also obtain a new nonperturbative result for the ratio of $c$~and $s$~quark masses: $m_c/m_s=11.652(65)$. This ratio implies $m_s(2,mathrm{GeV}, n_f=3)=93.6(8)$MeV when it is combined with our new~$c$~mass. Combining $m_c/m_s$ with our earlier $m_b/m_c$ gives $m_b/m_s=52.55(55)$, which is several standard deviations (but only 4%) away from the Georgi-Jarlskop prediction from certain GUTs. Finally we obtain an $n_f=4$ estimate for $m_b/m_c=4.528(54)$ which agrees well with our earlier $n_f=3$ result. The new ratio implies~$m_b(m_b,n_f=5)=4.162(48)$GeV.
We demonstrate that the O(a) taste mixing exhibited in standard textbook presentations of staggered quarks is an artifact of the particular definition of the flavor fields in those presentations, and has nothing to do with the underlying precision of
staggered-quark actions, despite continuing comments to the contrary in the current literature. To illustrate this point we introduce a new coordinate-space definition of the flavor fields that suppresses the O(a) term by two additional powers of a. In fact there are no errors at all from this mechanism. The only source of taste mixing comes from the exchange of highly-virtual gluons and enters in O(a^2). We review the idiosyncrasies of Symanzik improvement for naive/staggerd-quark actions, and show how these results follow from that program.
We illustrate a technique for fitting lattice QCD correlators to sums of exponentials that is significantly faster than traditional fitting methods --- 10--40 times faster for the realistic examples we present. Our examples are drawn from a recent an
alysis of the Upsilon spectrum, and another recent analysis of the D -> pi semileptonic form factor. For single correlators, we show how to simplify traditional effective-mass analyses.
We extend our earlier lattice-QCD analysis of heavy-quark correlators to smaller lattice spacings and larger masses to obtain new values for the c mass and QCD coupling, and, for the first time, values for the b mass: m_c(3GeV,n_f=4)=0.986(6)GeV, alp
ha_msb(M_Z,n_f=5)=0.1183(7), and m_b(10GeV,n_f=5)=3.617(25)GeV. These are among the most accurate determinations by any method. We check our results using a nonperturbative determination of the mass ratio m_b(mu,n_f)/m_c(mu,n_f); the two methods agree to within our 1% errors and taken together imply m_b/m_c=4.51(4). We also update our previous analysis of alpha_msb from Wilson loops to account for revised values for r_1 and r_1/a, finding a new value alpha_msb(M_Z,n_f=5)=0.1184(6); and we update our recent values for light-quark masses from the ratio m_c/m_s. Finally, in the Appendix, we derive a procedure for simplifying and accelerating complicated least-squares fits.
We use lattice QCD simulations, with MILC configurations (including vacuum polarization from u, d, and s quarks), to update our previous determinations of the QCD coupling constant. Our new analysis uses results from 6 different lattice spacings and
12 different combinations of sea-quark masses to significantly reduce our previous errors. We also correct for finite-lattice-spacing errors in the scale setting, and for nonperturbative chiral corrections to the 22 short-distance quantities from which we extract the coupling. Our final result is alpha_V(7.5GeV,nf=3) = 0.2120(28), which is equivalent to alpha_msbar(M_Z,n_f=5)= 0.1183(8). We compare this with our previous result, which differs by one standard deviation.
We use lattice QCD simulations, with MILC gluon configurations and HISQ c-quark propagators, to make very precise determinations of moments of charm-quark pseudoscalar, vector and axial-vector correlators. These moments are combined with new four-loo
p results from continuum perturbation theory to obtain several new determinations of the MSbar mass of the charm quark and of the MSbar coupling. We find m_c(3GeV)=0.986(10)GeV, or, equivalently, m_c(m_c)=1.268(9)GeV, both for n_f=4 flavors; and alpha_msb(3GeV,n_f=4)=0.251(6), or, equivalently, alpha_msb(M_Z,n_f=5)=0.1174(12). The new mass agrees well with results from continuum analyses of the vector correlator using experimental data for e+e- annihilation (instead of using lattice QCD simulations). These lattice and continuum results are the most accurate determinations to date of this mass. Ours is also one of the most accurate determinations of the QCD coupling by any method.
