No Arabic abstract
We calculate the mass difference between the $Upsilon$ and $eta_b$ and the $Upsilon$ leptonic width from lattice QCD using the Highly Improved Staggered Quark formalism for the $b$ quark and including $u$, $d$, $s$ and $c$ quarks in the sea. We have results for lattices with lattice spacing as low as 0.03 fm and multiple heavy quark masses, enabling us to map out the heavy quark mass dependence and determine values at the $b$ quark mass. Our results are: $M_{Upsilon} -M_{eta_b} = 57.5(2.3)(1.0) ,mathrm{MeV}$ (where the second uncertainty comes from neglect of quark-line disconnected correlation functions) and decay constants, $f_{eta_b}=724(12)$ MeV and $f_{Upsilon} =677.2(9.7)$ MeV, giving $Gamma(Upsilon rightarrow e^+e^-) = 1.292(37)(3) ,mathrm{keV}$. The hyperfine splitting and leptonic width are both in good agreement with experiment, and provide the most accurate lattice QCD results to date for these quantities by some margin. At the same time results for the time moments of the vector-vector correlation function can be compared to values for the $b$ quark contribution to $sigma(e^+e^- rightarrow mathrm{hadrons})$ determined from experiment. Moments 4--10 provide a 2% test of QCD and yield a $b$ quark contribution to the anomalous magnetic moment of the muon of 0.300(15)$times 10^{-10}$. Our results, covering a range of heavy quark masses, may also be useful to constrain QCD-like composite theories for beyond the Standard Model physics.
We determine the decay rate to leptons of the ground-state $Upsilon$ meson and its first radial excitation in lattice QCD for the first time. We use radiatively-improved NRQCD for the $b$ quarks and include $u$, $d$, $s$ and $c$ quarks in the sea with $u/d$ masses down to their physical values. We find $Gamma(Upsilon rightarrow e^+e^-)$ = 1.19(11) keV and $Gamma(Upsilon^{prime} rightarrow e^+e^-)$ = 0.69(9) keV, both in good agreement with experiment. The decay constants we obtain are included in a summary plot of meson decay constants from lattice QCD given in the Conclusions. We also test time-moments of the vector current-current correlator against values determined from the $b$ quark contribution to $sigma(e^+e^- rightarrow mathrm{hadrons})$ and calculate the $b$-quark piece of the hadronic vacuum polarisation contribution to the anomalous magnetic moment of the muon, $a_{mu}^b = 0.271(37) times 10^{-10}$. Finally we determine the $b$-quark mass, obtaining in the $overline{MS}$ scheme, $overline{m}_b(overline{m}_b, n_f=5)$ = 4.196(23) GeV, the most accurate result from lattice QCD to date.
We have performed the first $n_f = 2+1+1$ lattice QCD computations of the properties (masses and decay constants) of ground-state charmonium mesons. Our calculation uses the HISQ action to generate quark-line connected two-point correlation functions on MILC gluon field configurations that include $u/d$ quark masses going down to the physical point, tuning the $c$ quark mass from $M_{J/psi}$ and including the effect of the $c$ quarks electric charge through quenched QED. We obtain $M_{J/psi}-M_{eta_c}$ (connected) = 120.3(1.1) MeV and interpret the difference with experiment as the impact on $M_{eta_c}$ of its decay to gluons, missing from the lattice calculation. This allows us to determine $Delta M_{eta_c}^{mathrm{annihiln}}$ =+7.3(1.2) MeV, giving its value for the first time. Our result of $f_{J/psi}=$ 0.4104(17) GeV, gives $Gamma(J/psi rightarrow e^+e^-)$=5.637(49) keV, in agreement with, but now more accurate than experiment. At the same time we have improved the determination of the $c$ quark mass, including the impact of quenched QED to give $overline{m}_c(3,mathrm{GeV})$ = 0.9841(51) GeV. We have also used the time-moments of the vector charmonium current-current correlators to improve the lattice QCD result for the $c$ quark HVP contribution to the anomalous magnetic moment of the muon. We obtain $a_{mu}^c = 14.638(47) times 10^{-10}$, which is 2.5$sigma$ higher than the value derived using moments extracted from some sets of experimental data on $R(e^+e^- rightarrow mathrm{hadrons})$. This value for $a_{mu}^c$ includes our determination of the effect of QED on this quantity, $delta a_{mu}^c = 0.0313(28) times 10^{-10}$.
We have developed two methods for handling $b$ quarks in lattice QCD. One uses NRQCD (now improved to include radiative corrections) and the other uses Highly Improved Staggered Quarks (HISQ), extrapolating to the $b$ quark from lighter masses and using multiple lattice spacings to control discretisation errors. Comparison of results for the two different methods gives confidence in estimates of lattice QCD systematic errors, since they are very different in these two cases. Here we show results for heavyonium hyperfine splittings and vector current-current correlator moments using HISQ quarks, to add to earlier results testing the heavy HISQ method with pseudoscalar mesons. We also show the form factor for $B rightarrow pi l u$ decay at zero recoil using NRQCD $b$ quarks and $u/d$ quarks with physical masses. This allows us to test the soft pion theorem relation ($f_0(q^2_{max})=f_B/f_{pi}$) accurately and we find good agreement as $M_{pi} rightarrow 0$. }
We show results for the Upsilon spectrum calculated in lattice QCD including for the first time vacuum polarization effects for light u and d quarks as well as s quarks. We use gluon field configurations generated by the MILC collaboration. The calculations compare the results for a variety of u and d quark masses, as well as making a comparison to quenched results (in which quark vacuum polarisation is ignored) and results with only u and d quarks. The b quarks in the Upsilon are treated in lattice Nonrelativistic QCD through NLO in an expansion in the velocity of the b quark. We concentrate on accurate results for orbital and radial splittings where we see clear agreement with experiment once u, d and s quark vacuum polarisation effects are included. This now allows a consistent determination of the parameters of QCD. We demonstrate this consistency through the agreement of the Upsilon and B spectrum using the same lattice bare b quark mass. A one-loop matching to continuum QCD gives a value for the b quark mass in full lattice QCD for the first time. We obtain m_b^{bar{MS}}(m_b^{bar{MS}}) = 4.4(3) GeV. We are able to give physical results for the heavy quark potential parameters, r_0 = 0.469(7) fm and r_1 = 0.321(5) fm. Results for the fine structure in the spectrum and the Upsilon leptonic width are also presented. We predict the Upsilon - eta_b splitting to be 61(14) MeV, the Upsilon^{prime} - eta_b^{prime} splitting as 30(19) MeV and the splitting between the h_b and the spin-average of the chi_b states to be less than 6 MeV. Improvements to these calculations that will be made in the near future are discussed.
We determine the contribution to the anomalous magnetic moment of the muon from the $alpha^2_{mathrm{QED}}$ hadronic vacuum polarization diagram using full lattice QCD and including $u/d$ quarks with physical masses for the first time. We use gluon field configurations that include $u$, $d$, $s$ and $c$ quarks in the sea at multiple values of the lattice spacing, multiple $u/d$ masses and multiple volumes that allow us to include an analysis of finite-volume effects. We obtain a result for $a_{mu}^{mathrm{HVP,LO}}$ of $667(6)(12)$, where the first error is from the lattice calculation and the second includes systematic errors from missing QED and isospin-breaking effects and from quark-line disconnected diagrams. Our result implies a discrepancy between the experimental determination of $a_{mu}$ and the Standard Model of 3$sigma$.