Search

We present HPQCD's improved scalar, vector and tensor form factors for $B \to K$ semileptonic decays, using the heavy-HISQ formalism for more accurate normalisation of the weak currents. Working with masses close to the physical $b$ on the finest ensemble and including three ensembles with physical light quarks, we cover the full physical $q^2$ range with good precision. Our uncertainties at $q^2=0$ are a factor of three better than earlier work. We compare Standard Model observables using our form factors to experimental measurements for the rare flavour changing neutral current processes $B \to K \ell^+\ell^-$ and $B \to K \nu\bar{\nu}$ and discuss the significance of the tensions that arise., Comment: 10 pages, 6 figures. Proceedings of the 39th International Symposium on Lattice Field Theory - LATTICE2022

We use HPQCD's recent lattice QCD determination of $B \to K$ scalar, vector and tensor form factors to determine Standard Model differential branching fractions for $B \to K \ell^+\ell^-$, $B\to K \ell_1^+\ell_2^-$ and $B \to K\nu \overline{\nu}$. These form factors are calculated across the full $q^2$ range of the decay and have smaller uncertainties than previous work, particularly at low $q^2$. For $B \to K \ell^+ \ell^-$ we find the Standard Model branching fraction in the $q^2$ region below the squared $J/\psi$ mass to exceed the LHCb results, with tensions as high as $4.2\sigma$ for $B^+\to K^+\mu^+\mu^-$. For the high $q^2$ region we see $2.7\sigma$ tensions. The tensions are much reduced by applying shifts to Wilson coefficients $C_9$ and $C_{10}$ in the effective weak Hamiltonian, moving them away from their Standard Model values consistent with those indicated by other $B$ phenomenology. We also update results for lepton-flavour ratios $R^{\mu}_e$ and $R^{\tau}_{\mu}$ and the `flat term', $F_H^{\ell}$ in the differential branching fraction for $\ell\in\{e,\mu,\tau\}$. Our results for the form-factor-dependent contributions needed for searches for lepton-flavour-violating decays $B\to K\ell^-_1\ell^+_2$ achieve uncertainties of 7%. We also compute the branching fraction $\mathcal{B}(B\to K\nu\bar{\nu})$ with an uncertainty below 10%, for comparison with future experimental results., Comment: 28 pages, 38 figures. Version accepted by Physical Review D. Updated to reflect a change in the scale of alpha_EW

We present the result of lattice QCD calculation of the scalar, vector and tensor form factors for the $B\to K\ell^+\ell^-$ decay, across the full physical range of momentum transfer. We use the highly improved staggered quark (HISQ) formalism for all valence quarks on eight ensembles of gluon field configurations generated by the MILC collaboration. These include four flavours of HISQ quarks in the sea, with three ensembles having the light $u/d$ quarks at physical masses. In the first fully relativistic calculation of these form factors, we use the heavy-HISQ method. This allows us to determine the form factors as a function of heavy quark mass from the $c$ to the $b$, and so we also obtain new results for the $D\to K$ tensor form factor. The advantage of the relativistic formalism is that we can match the lattice weak currents to their continuum counterparts much more accurately than in previous calculations; our scalar and vector currents are renormalised fully nonperturbatively and we use a well-matched intermediate momentum-subtraction scheme for our tensor current. Our scalar and vector $B\to K$ form factors have uncertainties of less than 4% across the entire physical $q^2$ range and the uncertainty in our tensor form factor is less than 7%. Our heavy-HISQ method allows us to map out the dependence on heavy-quark mass of the form factors and we can also see the impact of changing spectator quark mass by comparing to earlier HPQCD results for the same quark weak transition but for heavier mesons., Comment: 30 pages, 22 figures. Version accepted by Physical Review D

An accurate determination of the leading-order hadronic vacuum polarisation (HVP) contribution to the anomalous magnetic moment of the muon is critical to understanding the size and significance of any discrepancy between the Standard Model prediction and experimental results being obtained by the Muon g-2 experiment at Fermilab. The Standard Model prediction is currently based on a data-driven approach to the HVP using experimental results for $\sigma(e^+e^-\rightarrow\,\mathrm{hadrons})$. Lattice QCD aims to provide a result with similar uncertainty from calculated vector-vector correlation functions, but the growth of statistical and systematic errors in the $u/d$ quark correlation functions at large Euclidean time has made this difficult to achieve. We show that restricting the lattice contributions to a one-sided window $0

The precision of lattice QCD calculations has been steadily improving for some time and is now approaching, or has surpassed, the 1% level for multiple quantities. At this level QED effects, i.e. the fact that quarks carry electric as well as color charge, come into play. In this report we will summarise results from the first lattice QCD+QED computations of the properties of ground-state charmonium and bottomonium mesons by the HPQCD Collaboration., Comment: Contribution to the proceedings of the 19th International Conference on Hadron Spectroscopy and Structure (HADRON 2021). Accepted for publication in Suplemento de la Revista Mexicana de F\'isica

As tests of QCD in the bottomonium system, we give the most accurate results to date for the ground-state hyperfine splitting and the $\Upsilon$ leptonic width from full lattice QCD. These quantities are both accurately known from experiment, so can provide a good test of $b$ physics, but previous lattice results have been rather imprecise. We also test the impact on these quantities of the $b$ quark's electric charge. Our results are: $M_{\Upsilon}-M_{\eta_b} = $ 57.5(2.3)(1.0) 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) keV. We also give a new determination of the ratio of the masses for $b$ and $c$ quarks that is completely nonperturbative in lattice QCD and includes the calculation of QED effects for the first time. This gives a result for the $b$ quark mass of $\overline{m}_b(\overline{m}_b,n_f=5) =$ 4.202(21) GeV., Comment: 6 pages, The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology

Semileptonic $D \to{}K \ell \nu$ decays provide one angle of attack to get at the CKM matrix element $V_{cs}$, complementary to the study of leptonic $D_s$ decays. Here, HPQCD present the results of a recently published, improved determination of $V_{cs}$. We discuss a new, precise determination of $D\to K$ scalar and vector form factors from a lattice calculation on eight different $N_f=2+1+1$ MILC gluon field ensembles using the HISQ action, including three with physical light quark masses. When combined with experimental results, we are able to extract $|V_{cs}|=0.9663(80)$ to a sub percent level of precision for the first time. This is achieved using three different methods, which each combine our form factors with different sets of experimental results in different ways, with the results in very good agreement. Our primary method is to use $q^2$-binned data for the differential decay rate, but we also calculate $V_{cs}$ from the total branching fraction and from the value $|V_{cs}|f_+(0)$, which is also quoted by some experiments., Comment: 9 pages, 7 figures. Proceedings of the 38th International Symposium on Lattice Field Theory - LATTICE2021

We provide a 0.8\%-accurate determination of $V_{cs}$ from combining experimental results for the differential rate of $D \rightarrow K$ semileptonic decays with precise form factors that we determine from lattice QCD. This is the first time that $V_{cs}$ has been determined with an accuracy that allows its difference from 1 to be seen. Our lattice QCD calculation uses the Highly Improved Staggered Quark (HISQ) action for all valence quarks on gluon field configurations generated by the MILC collaboration that include the effect of $u$, $d$, $s$ and $c$ HISQ quarks in the sea. We use eight gluon field ensembles with five values of the lattice spacing ranging from 0.15 fm to 0.045 fm and include results with physical $u/d$ quarks for the first time. Our calculated form factors cover the full $q^2$ range of the physical decay process and enable a Standard Model test of the shape of the differential decay rate as well as the determination of $V_{cs}$ from a correlated weighted average over $q^2$ bins. We obtain $|V_{cs}|= 0.9663(53)_{\text{latt}}(39)_{\text{exp}}(19)_{\eta_{EW}}(40)_{\text{EM}}$, where the uncertainties come from lattice QCD, experiment, short-distance electroweak and electromagnetic corrections, respectively. This last uncertainty, neglected for $D \rightarrow K \ell \nu$ hitherto, now needs attention if the uncertainty on $V_{cs}$ is to be reduced further. We also determine $V_{cs}$ values in good agreement using the measured total branching fraction and the rates extrapolated to $q^2=0$. Our form factors enable tests of lepton flavour universality violation. We find the ratio of branching fractions for $D^0 \rightarrow K^-$ with $\mu$ and $e$ in the final state to be $R_{\mu/e}=0.9779(2)_{\text{latt}}(50)_{\mathrm{EM}}$ in the Standard Model, with the uncertainty dominated by that from electromagnetic corrections., Comment: 27 pages, 28 figures. Small changes to text to improve readability. Updated to include the BES average for the D to K branching fraction from 2104.08081

We extend HPQCD's earlier $n_f=4$ lattice-QCD analysis of the ratio of $\overline{\mathrm{MSB}}$ masses of the $b$ and $c$ quark to include results from finer lattices (down to 0.03fm) and a new calculation of QED contributions to the mass ratio. We find that $\overline{m}_b(\mu)/\overline{m}_c(\mu)=4.586(12)$ at renormalization scale $\mu=3$\,GeV. This result is nonperturbative. Combining it with HPQCD's recent lattice QCD$+$QED determination of $\overline{m}_c(3\mathrm{GeV})$ gives a new value for the $b$-quark mass: $\overline{m}_b(3\mathrm{GeV}) = 4.513(26)$GeV. The $b$-mass corresponds to $\overline{m}_b(\overline{m}_b, n_f=5) = 4.202(21)$GeV. These results are the first based on simulations that include QED., Comment: 8 pages, 4 figures

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., Comment: 16 pages, 9 figures Version accepted by PRD: Appendix with tables of fit parameters that allow the dependence on heavy quark mass of the hyperfine splitting and decay constants to be reconstructed added

We present the results of a lattice QCD calculation of the scalar and vector form factors for the unphysical $B_s\to\eta_s$ decay, over the full physical range of $q^2$. This is a useful testing ground both for lattice QCD and for our wider understanding of the behaviour of form factors. Calculations were performed using the highly improved staggered quark (HISQ) action on $N_f = 2 + 1 + 1$ gluon ensembles generated by the MILC Collaboration with an improved gluon action and HISQ sea quarks. We use three lattice spacings and a range of heavy quark masses from that of charm to bottom, all in the HISQ formalism. This permits an extrapolation in the heavy quark mass and lattice spacing to the physical point and nonperturbative renormalisation of the vector matrix element on the lattice. We find results in good agreement with previous work using nonrelativistic QCD $b$ quarks and with reduced errors at low $q^2$, supporting the effectiveness of our heavy HISQ technique as a method for calculating form factors involving heavy quarks. A comparison with results for other decays related by SU(3) flavour symmetry shows that the impact of changing the light daughter quark is substantial but changing the spectator quark has very little effect. We also map out form factor shape parameters as a function of heavy quark mass and compare to heavy quark effective theory expectations for mass scaling at low and high recoil. This work represents an important step in the progression from previous work on heavy-to-heavy decays ($b\to c$) to the numerically more challenging heavy-to-light decays., Comment: published version; 18 pages, 14 figures

Hadron masses are subject to few MeV corrections arising from QED interactions, almost entirely arising from the electric charge of the valence quarks. The QED effects include both self-energy contributions and interactions between the valence quarks/anti-quarks. By combining results from different signs of the valence quark electric charge we are able to isolate the interaction term which is dominated by the Coulomb piece, $\langle \alpha_{\mathrm{QED}}e_{q_1}e_{\overline{q}_2}/r \rangle$, in the nonrelativistic limit. We study this for $D_s$, $\eta_c$ and $J/\psi$ mesons, working in lattice QCD plus quenched QED. We use gluon field configurations that include up, down, strange and charm quarks in the sea at multiple values of the lattice spacing. Our results, including also values for mesons with quarks heavier than charm, can be used to improve phenomenological models for the QED contributions. The QED interaction term carries information about meson structure; we derive effective sizes $\langle 1/r_{\mathrm{eff}} \rangle^{-1}$ for $\eta_c$, $J/\psi$ and $D_s$ of 0.206(8) fm, 0.321(14) fm and 0.307(31) fm respectively., Comment: 8 pages, 3 figures

Lattice QCD calculations of form factors for rare Standard Model processes such as $B \to K \ell^+ \ell^-$ use tensor currents that require renormalisation. These renormalisation factors, $Z_T$, have typically been calculated within perturbation theory and the estimated uncertainties from missing higher order terms are significant. Here we study tensor current renormalisation using lattice implementations of momentum-subtraction schemes. Such schemes are potentially more accurate but have systematic errors from nonperturbative artefacts. To determine and remove these condensate contributions we calculate the ground-state charmonium tensor decay constant, $f_{J/\psi}^T$, which is also of interest in beyond the Standard Model studies. We obtain $f_{J/\psi}^T(\bar{\text{MS}}, 2\ \mathrm{GeV})=0.3927(27)$ GeV, with ratio to the vector decay constant of 0.9569(52), significantly below 1. We also give $Z_T$ factors, converted to the $\bar{\mathrm{MS}}$ scheme, corrected for condensate contamination. This contamination reaches 1.5\% at a renormalisation scale of 2 GeV (in the preferred RI-SMOM scheme) and so must be removed for accurate results., Comment: 12 pages, version accepted for publication in PRD

Two of the elements of the Cabibbo-Kobayashi-Maskawa quark mixing matrix, $|V_{ub}|$ and $|V_{cb}|$, are extracted from semileptonic B decays. The results of the B factories, analysed in the light of the most recent theoretical calculations, remain puzzling, because for both $|V_{ub}|$ and $|V_{cb}|$ the exclusive and inclusive determinations are in clear tension. Further, measurements in the $\tau$ channels at Belle, Babar, and LHCb show discrepancies with the Standard Model predictions, pointing to a possible violation of lepton flavor universality. LHCb and Belle II have the potential to resolve these issues in the next few years. This article summarizes the discussions and results obtained at the MITP workshop held on April 9--13, 2018, in Mainz, Germany, with the goal to develop a medium-term strategy of analyses and calculations aimed at solving the puzzles. Lattice and continuum theorists working together with experimentalists have discussed how to reshape the semileptonic analyses in view of the much higher luminosity expected at Belle II, searching for ways to systematically validate the theoretical predictions in both exclusive and inclusive B decays, and to exploit the rich possibilities at LHCb., Comment: 77 pages

We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $\alpha$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(\alpha^5)$ with negligible numerical uncertainty. The electroweak contribution is suppressed by $(m_\mu/M_W)^2$ and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at $\mathcal{O}(\alpha^2)$ and is due to hadronic vacuum polarization, whereas at $\mathcal{O}(\alpha^3)$ the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first-principle, lattice-QCD approach. The final result reads $a_\mu^\text{SM}=116\,591\,810(43)\times 10^{-11}$ and is smaller than the Brookhaven measurement by 3.7$\sigma$. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. This and the prospects to further reduce the theoretical uncertainty in the near future-which are also discussed here-make this quantity one of the most promising places to look for evidence of new physics., Comment: 196 pages, 103 figures, version published in Phys. Rept., bib files for the citation references are available from: https://muon-gm2-theory.illinois.edu

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$ quark's 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}$., Comment: Added extra discussion on QED setup, some new results to study the effects of strong isospin breaking in the sea (including new Fig. 1) and a fit stability plot for the hyperfine splitting (new Fig. 7). Version accepted for publication in PRD

We present an update on the ongoing calculations by the Fermilab Lattice, HPQCD, and MILC Collaboration of the leading-order (in electromagnetism) hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon. Our project employs ensembles with four flavors of highly improved staggered fermions, physical light-quark masses, and four lattice spacings ranging from $a \approx 0.06$ to 0.15 fm for most of the results thus far., Comment: LATTICE 2019, 7 pages, 7 figures

We examine the renormalisation of flavour-diagonal vector currents in lattice QCD with the aim of understanding and quantifying the systematic errors from nonperturbative artefacts associated with the use of intermediate momentum-subtraction schemes. Our study uses the Highly Improved Staggered Quark (HISQ) action on gluon field configurations that include $n_f=2+1+1$ flavours of sea quarks, but our results have applicability to other quark actions. Renormalisation schemes that make use of the exact lattice vector Ward-Takahashi identity for the conserved current also have renormalisation factors, $Z_V$, for nonconserved vector currents that are free of contamination by nonperturbative condensates. We show this by explicit comparison of two such schemes: that of the vector form factor at zero momentum transfer and the RI-SMOM momentum-subtraction scheme. The two determinations of $Z_V$ differ only by discretisation effects (for any value of momentum-transfer in the RI-SMOM case). The RI$^{\prime}$-MOM scheme, although widely used, does not share this property. We show that $Z_V$ determined in the standard way in this scheme has $\mathcal{O}(1\%)$ nonperturbative contamination that limits its accuracy. Instead we define an RI$^{\prime}$-MOM $Z_V$ from a ratio of local to conserved vector current vertex functions and show that this $Z_V$ is a safe one to use in lattice QCD calculations. We also perform a first study of vector current renormalisation with the inclusion of quenched QED effects on the lattice using the RI-SMOM scheme., Comment: Minor changes to the text and references added. Version accepted by Phys. Rev. D

As the only lattice vector current that does not require renormalisation is the point-split conserved current it is convenient to have a robust, precise and computationally cheap methodology for the calculation of vector current renormalisation factors, $Z_V$. Momentum subtraction schemes, such as RI-SMOM, implemented nonperturbatively on the lattice provide such a method if it can be shown that the systematic errors, e.g. from condensates, are well controlled. We present $Z_V$ calculations for the conserved current in both the RI-SMOM and RI$'$-MOM momentum subtraction schemes as well as local current renormalisation in the RI-SMOM scheme. By performing these calculations at various values of the momentum scale $\mu$ and different lattice spacings we can investigate the presence of power suppressed nonperturbative contributions and compare the results to expectations arising from the Ward-Takahashi identity. Our results show that the RI-SMOM scheme provides a well controlled determination of $Z_V$ but the standard RI$'$-MOM scheme does not. We then present some preliminary uses of these $Z_V$ calculations in charm physics., Comment: Proceedings of Lattice 2019

We calculate the bag parameters for neutral $B$-meson mixing in and beyond the Standard Model, in full four-flavour lattice QCD for the first time. We work on gluon field configurations that include the effect of $u$, $d$, $s$ and $c$ sea quarks with the Highly Improved Staggered Quark (HISQ) action at three values of the lattice spacing and with three $u/d$ quark masses going down to the physical value. The valence $b$ quarks use the improved NRQCD action and the valence light quarks, the HISQ action. Our analysis was blinded. Our results for the bag parameters for all five operators are the most accurate to date. For the Standard Model operator between $B_s$ and $B_d$ mesons we find: $\hat{B}_{B_s}=1.232(53)$, $\hat{B}_{B_d}=1.222(61)$. Combining our results with lattice QCD calculations of the decay constants using HISQ quarks from the Fermilab/MILC collaboration and with experimental values for $B_s$ and $B_d$ oscillation frequencies allows determination of the CKM elements $V_{ts}$ and $V_{td}$. We find $V_{ts} = 0.04189(93)$, $V_{td} = 0.00867(23)$ and $V_{ts}/V_{td} = 0.2071(27)$. Our results agree well (within $2\sigma$) with values determined from CKM unitarity constraints based on tree-level processes (only). Using a ratio to $\Delta M$ in which CKM elements cancel in the Standard Model, we determine the branching fractions ${\text{Br}}(B_s\rightarrow \mu^+\mu^-) = 3.81(18) \times 10^{-9}$ and ${\text{Br}}(B_d\rightarrow \mu^+\mu^-) = 1.031(54) \times 10^{-10}$. We also give results for matrix elements of the operators $R_0$, $R_1$ and $\tilde{R}_1$ that contribute to neutral $B$-meson width differences., Comment: 27 pages, 16 figures. Slight changes to text to clarify arguments. Updated Figure 10. Version accepted by Physical Review D

We present a lattice QCD determination of the $B_s \to D_s \ell\nu$ scalar and vector form factors over the full physical range of momentum transfer. The result is derived from correlation functions computed using the Highly Improved Staggered Quark (HISQ) formalism, on the second generation MILC gluon ensembles accounting for up, down, strange and charm contributions from the sea. We calculate correlation functions for three lattice spacing values and an array of unphysically light $b$-quark masses, and extrapolate to the physical value. Using the HISQ formalism for all quarks means that the lattice current coupling to the $W$ can be renormalized non-perturbatively, giving a result free from perturbative matching errors for the first time. Our results are in agreement with, and more accurate than, previous determinations of these form factors. From the form factors we also determine the ratio of branching fractions that is sensitive to violation of lepton universality: $R(D_s) = \mathcal{B}(B_s\to D_s \tau \nu_{\tau})/\mathcal{B}(B_s\to D_s \ell \nu_{l})$, where $\ell$ is an electron or a muon. We find $R(D_s) = 0.2987(46)$, which is also more accurate than previous lattice QCD results. Combined with a future measurement of $R(D_s)$, this could supply a new test of the Standard Model. We also compare the dependence on heavy quark mass of our form factors to expectations from Heavy Quark Effective Theory., Comment: 22 pages, 17 figures

We present details of a lattice QCD calculation of the $B_s\to D_s^*$ axial form factor at zero recoil using the Highly Improved Staggered Quark (HISQ) formalism on the second generation MILC gluon ensembles that include up, down, strange and charm quarks in the sea. Using the HISQ action for all valence quarks means that the lattice axial vector current that couples to the $W$ can be renormalized fully non-perturbatively, giving a result free of the perturbative matching errors that previous lattice QCD calculations have had. We calculate correlation functions at three values of the lattice spacing, and multiple `$b$'-quark masses, for physical $c$ and $s$. The functional dependence on the $b$-quark mass can be determined and compared to Heavy Quark Effective Theory expectations, and a result for the form factor obtained at the physical value of the $b$-quark mass. We find $\mathcal{F}^{B_s\to D_s^*}(1) = h^s_{A_1}(1) = 0.9020(96)_{\text{stat}}(90)_{\text{sys}}$. This is in agreement with earlier lattice QCD results, which use NRQCD $b$ quarks, with a total uncertainty reduced by more than a factor of two. We discuss implications of this result for the $B\to D^*$ axial form factor at zero recoil and for determinations of $V_{cb}$., Comment: 18 pages, 9 figures

Lattice QCD can provide a direct determination of meson electromagnetic form factors, making predictions for upcoming experiments at Jefferson Lab. The form factors are a reflection of the bound-state nature of the meson and so these calculations give information about how confinement by QCD affects meson internal structure. The region of high squared (space-like) momentum-transfer, $Q^2$, is of particular interest because perturbative QCD predictions take a simple form in that limit that depends on the meson decay constant. We previously showed in\cite{jonnaff} that, up to $Q^2$ of 6 $\mathrm{GeV}^2$, the form factor for a `pseudo-pion' made of strange quarks was significantly larger than the asymptotic perturbative QCD result and showed no sign of heading towards that value at higher $Q^2$. Here we give predictions for real mesons, the $K^+$ and $K^0$, in anticipation of JLAB results for the $K^+$ in the next few years. We also give results for a heavier meson, the $\eta_c$, up to $Q^2$ of 25 $\mathrm{GeV}^2$ for a comparison to perturbative QCD in a higher $Q^2$ regime., Comment: 6 pages, 5 figures, Proceedings of LATTICE2018

We calculate the contribution to the muon anomalous magnetic moment hadronic vacuum polarization from {the} connected diagrams of up and down quarks, omitting electromagnetism. We employ QCD gauge-field configurations with dynamical $u$, $d$, $s$, and $c$ quarks and the physical pion mass, and analyze five ensembles with lattice spacings ranging from $a \approx 0.06$ to~0.15~fm. The up- and down-quark masses in our simulations have equal masses $m_l$. We obtain, in this world where all pions have the mass of the $\pi^0$, $10^{10} a_\mu^{ll}({\rm conn.}) = 637.8\,(8.8)$, in agreement with independent lattice-QCD calculations. We then combine this value with published lattice-QCD results for the connected contributions from strange, charm, and bottom quarks, and an estimate of the uncertainty due to the fact that our calculation does not include strong-isospin breaking, electromagnetism, or contributions from quark-disconnected diagrams. Our final result for the total $\mathcal{O}(\alpha^2)$ hadronic vacuum polarization to the muon's anomalous magnetic moment is~$10^{10}a_\mu^{\rm HVP,LO} = 699(15)_{u,d}(1)_{s,c,b}$, where the errors are from the light-quark and heavy-quark contributions, respectively. Our result agrees with both {\it ab-initio} lattice-QCD calculations and phenomenological determinations from experimental $e^+e^-$-scattering data. It is $1.3\sigma$ below the "no new physics" value of the hadronic-vacuum-polarization contribution inferred from combining the BNL E821 measurement of $a_\mu$ with theoretical calculations of the other contributions., Comment: 19 pages, 12 figures, 6 tables; updated to correct a small mistake in the finite volume correction resulting in small changes to the results, matches published version

Lattice QCD provides several avenues for the high precision determination of quark masses. Using the RI-SMOM scheme applied to lattice calculations with the HISQ action, we obtain mass renormalisation factors that we use to provide strange and charm quark masses with 1% precision. The calculation involves the study of various sources of systematic uncertainty, including an analysis of possible nonperturbative (condensate) contributions. These results allow a comparison of different mass determination methods of comparable precision. In particular we (HPQCD) find good agreement between RI-SMOM and current-current correlator determinations based on the same lattice QCD bare masses, providing a strong test of our understanding of systematic uncertainties., Comment: Proceedings from Lattice 2018

We present progress on an ongoing calculation of the $B_s\to D_s^{(*)} l \nu$ form factors calculated on the $n_f=2+1+1$ MILC ensembles and using the Highly Improved Staggered Quark action for all valence quarks. We perform the calculation at a range of $b$ quark masses (and lattice spacings) so that we can extrapolate to the physical $b$-quark mass., Comment: 6 pages, 3 figures, proceedings of the 36th Annual International Symposium on Lattice Field Theory - LATTICE2018

We derive the Operator Product Expansion whose vacuum expectation value gives the time-moments of the pseudoscalar heavy-light current-current correlator up to and including terms in $\alpha_s^2$ multiplying $\langle\overline{\psi}\psi\rangle/M^3$ and terms in $\alpha_s$ multiplying $\langle \alpha_s G^2 \rangle/M^4$, where $M$ is the heavy-quark mass. Using lattice QCD results for heavy-strange correlators obtained for a variety of heavy quark masses on gluon field configurations including $u$, $d$ and $s$ quarks in the sea at three values of the lattice spacing, we are able to show that the contribution of the strange-quark condensate to the time-moments is very substantial. We use our lattice QCD time-moments and the OPE to determine a value for the condensate, fitting the 4th, 6th, 8th and 10th time-moments simultaneously. Our result, $\langle \overline{s}s \rangle^{\overline{\text{MS}}}(2 \text{GeV}) = -(296(11) \,\mathrm{MeV})^3$, agrees well with HPQCD's earlier, more direct, lattice QCD determination~\cite{McNeile:2012xh}. As well as confirming that the $s$ quark condensate is close in value to the light quark condensate, this demonstrates clearly the consistency of the Operator Product Expansion for fully nonperturbative calculations of matrix elements of short-distance operators in lattice QCD., Comment: 14 pages, 4 figures

We determine the charm and strange quark masses in the $\overline{\text{MS}}$ scheme, using $n_f=2+1+1$ lattice QCD calculations with highly improved staggered quarks (HISQ) and the RI-SMOM intermediate scheme to connect the bare lattice quark masses to continuum renormalisation schemes. Our study covers analysis of systematic uncertainties from this method, including nonperturbative artefacts and the impact of the non-zero physical sea quark masses. We find $m_c^{\overline{\text{MS}}}(3 \text{GeV}) = 0.9896(61)$ GeV and $m_s^{\overline{\text{MS}}}(3 \text{GeV}) = 0.08536(85)$ GeV, where the uncertainties are dominated by the tuning of the bare lattice quark masses. These results are consistent with, and of similar accuracy to, those using the current-current correlator approach coupled to high-order continuum QCD perturbation theory, implemented in the same quark formalism and on the same gauge field configurations. This provides a strong test of the consistency of methods for determining the quark masses to high precision from lattice QCD. We also give updated lattice QCD world averages for $c$ and $s$ quark masses., Comment: 23 pages, 14 figures

We determine the normalisation of scalar and pseudoscalar current operators made from non-relativistic $b$ quarks and Highly Improved Staggered light quarks in lattice Quantum Chromodynamics (QCD) through $\mathcal{O}(\alpha_s)$ and $\Lambda_{\text{QCD}}/m_b$. We use matrix elements of these operators to extract $B$ meson decay constants and form factors, then compare to those obtained using the standard vector and axial-vector operators. This provides a test of systematic errors in the lattice QCD determination of the $B$ meson decay constants and form factors. We provide a new value for the $B$ and $B_s$ meson decay constants from lattice QCD calculations on ensembles that include $u$, $d$, $s$ and $c$ quarks in the sea and those which have the $u/d$ quark mass going down to its physical value. Our results are $f_B=0.196(6)$ GeV, $f_{B_s}=0.236(7)$ GeV and $f_{B_s}/f_B =1.207(7)$, agreeing well with earlier results using the temporal axial current. By combining with these previous results, we provide updated values of $f_B=0.190(4)$ GeV, $f_{B_s}=0.229(5)$ GeV and $f_{B_s}/f_B = 1.206(5)$., Comment: 14 pages, 10 figures

All lattice-QCD calculations of the hadronic-vacuum-polarization contribution to the muon's anomalous magnetic moment to-date have been performed with degenerate up- and down-quark masses. Here we calculate directly the strong-isospin-breaking correction to $a_\mu^{\rm HVP}$ for the first time with physical values of $m_u$ and $m_d$ and dynamical $u$, $d$, $s$, and $c$ quarks, thereby removing this important source of systematic uncertainty. We obtain a relative shift to be applied to lattice-QCD results obtained with degenerate light-quark masses of $\delta a_\mu^{{\rm HVP,} m_u \neq m_d}$= +1.5(7)%, in agreement with estimates from phenomenology and a recent lattice-QCD calculation with unphysically heavy pions., Comment: v2: 6 pages, 2 tables, 2 figures. Additional references and expanded discussion of systematic errors. Version accepted to Physical Review Letters

We compare correlators for pseudoscalar and vector mesons made from valence strange quarks using the clover quark and highly improved staggered quark (HISQ) formalisms in full lattice QCD. We use fully nonperturbative methods to normalise vector and axial vector current operators made from HISQ quarks, clover quarks and from combining HISQ and clover fields. This allows us to test expectations for the renormalisation factors based on perturbative QCD, with implications for the error budget of lattice QCD calculations of the matrix elements of clover-staggered $b$-light weak currents, as well as further HISQ calculations of the hadronic vacuum polarisation. We also compare the approach to the (same) continuum limit in clover and HISQ formalisms for the mass and decay constant of the $\phi$ meson. Our final results for these parameters, using single-meson correlators and neglecting quark-line disconnected diagrams are: $m_{\phi} =$ 1.023(6) GeV and $f_{\phi} = $ 0.238(3) GeV in good agreement with experiment. The results come from calculations in the HISQ formalism using gluon fields that include the effect of $u$, $d$, $s$ and $c$ quarks in the sea with three lattice spacing values and $m_{u/d}$ values going down to the physical point., Comment: 21 pages, 11 figures. Version published by Phys.Rev.D. Minor changes to text. HISQ Z_V value derived in Appendix B for a=0.06fm lattices changed slightly

We give an accurate determination of the vector (electromagnetic) form factor, $F(Q^2)$, for a light pseudoscalar meson up to squared momentum transfer $Q^2$ values of 6 $\mathrm{GeV}^2$ for the first time from full lattice QCD, including $u$, $d$, $s$ and $c$ quarks in the sea at multiple values of the lattice spacing. Our results show good control of lattice discretisation and sea quark mass effects. We study a pseudoscalar meson made of valence $s$ quarks but the qualitative picture obtained applies also to the $\pi$ meson, relevant to upcoming experiments at Jefferson Lab. We find that $Q^2F(Q^2)$ becomes flat in the region between $Q^2$ of 2 $\mathrm{GeV}^2$ and 6 $\mathrm{GeV}^2$, with a value well above that of the asymptotic perturbative QCD expectation, but well below that of the vector-meson dominance pole form appropriate to low $Q^2$ values. Our calculations show that we can reach higher $Q^2$ values in future to shed further light on where the perturbative QCD result emerges., Comment: 8 pages, 6 figures. Version accepted by Physical Review D. Includes results from a third value of the lattice spacing, change to the title and two additional figures in the appendix discussing analysis of systematic and statistical errors

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$., Comment: 14 pages, 10 figures. Discussion of method extended with additional tests and figures added. Typographical errors corrected

The quark-line disconnected diagram is a potentially important ingredient in lattice QCD calculations of the hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon. It is also a notoriously difficult one to evaluate. Here, for the first time, we give an estimate of this contribution based on lattice QCD results that have a statistically significant signal, albeit at one value of the lattice spacing and an unphysically heavy value of the $u/d$ quark mass. We use HPQCD's method of determining the anomalous magnetic moment by reconstructing the Adler function from time-moments of the current-current correlator at zero spatial momentum. Our results lead to a total (including $u$, $d$ and $s$ quarks) quark-line disconnected contribution to $a_{\mu}$ of $-0.15\%$ of the $u/d$ hadronic vacuum polarization contribution with an uncertainty which is 1\% of that contribution., Comment: 9 pages, 3 figures

We present the first calculation of the electromagnetic form factor of the $\pi$ meson at physical light quark masses. We use configurations generated by the MILC collaboration including the effect of $u$, $d$, $s$ and $c$ sea quarks with the Highly Improved Staggered Quark formalism. We work at three values of the lattice spacing on large volumes and with $u$/$d$ quark masses going down to the physical value. We study scalar and vector form factors for a range in space-like $q^2$ from 0.0 to -0.1 $\mathrm{GeV}^2$ and from their shape we extract mean square radii. Our vector form factor agrees well with experiment and we find $\langle r^2 \rangle_V = 0.403(18)(6) \,\mathrm{fm}^2$. For the scalar form factor we include quark-line disconnected contributions which have a significant impact on the radius. We give the first results for SU(3) flavour-singlet and octet scalar mean square radii, obtaining: $\langle r^2 \rangle_S^{\mathrm{singlet}} = 0.506(38)(53) \mathrm{fm}^2$ and $\langle r^2 \rangle_S^{\mathrm{octet}} = 0.431(38)(46) \mathrm{fm}^2$. We discuss the comparison with expectations from chiral perturbation theory., Comment: 13 pages, 11 figures

The exclusive semileptonic decay $B \rightarrow \pi \ell \nu$ is a key process for the determination of the Cabibbo-Kobayashi-Maskawa matrix element $V_{ub}$ from the comparison of experimental rates as a function of $q^2$ with theoretically determined form factors. The sensitivity of the form factors to the $u/d$ quark mass has meant significant systematic uncertainties in lattice QCD calculations at unphysically heavy pion masses. Here we give the first lattice QCD calculations of this process for u/d quark masses going down to their physical values, calculating the $f_0$ form factor at zero recoil to 3\%. We are able to resolve a long-standing controversy by showing that the soft-pion theorem result $f_0(q^2_{max}) = f_B/f_{\pi}$ does hold as $m_{\pi} \rightarrow 0$. We use the Highly Improved Staggered Quark formalism for the light quarks and show that staggered chiral perturbation theory for the $m_{\pi}$ dependence is almost identical to continuum chiral perturbation theory for $f_0$, $f_B$ and $f_{\pi}$. We also give results for other processes such as $B_s \rightarrow K \ell \nu$., Comment: 16 pages, 6 figures

We extend the picture of $B$-meson decay constants obtained in lattice QCD beyond those of the $B$, $B_s$ and $B_c$ to give the first full lattice QCD results for the $B^*$, $B^*_s$ and $B^*_c$. We use improved NonRelativistic QCD for the valence $b$ quark and the Highly Improved Staggered Quark (HISQ) action for the lighter quarks on gluon field configurations that include the effect of $u/d$, $s$ and $c$ quarks in the sea with $u/d$ quark masses going down to physical values. For the ratio of vector to pseudoscalar decay constants, we find $f_{B^*}/f_B$ = 0.941(26), $f_{B^*_s}/f_{B_s}$ = 0.953(23) (both $2\sigma$ less than 1.0) and $f_{B^*_c}/f_{B_c}$ = 0.988(27). Taking correlated uncertainties into account we see clear indications that the ratio increases as the mass of the lighter quark increases. We compare our results to those using the HISQ formalism for all quarks and find good agreement both on decay constant values when the heaviest quark is a $b$ and on the dependence on the mass of the heaviest quark in the region of the $b$. Finally, we give an overview plot of decay constants for gold-plated mesons, the most complete picture of these hadronic parameters to date., Comment: 20 pages, 9 figures. Minor updates to the discussion in several places and some additional references

We present the first lattice QCD calculation of the $B_s$ and $B_d$ mixing parameters with physical light quark masses. We use MILC gluon field configurations that include $u$, $d$, $s$ and $c$ sea quarks at 3 values of the lattice spacing and with 3 values of the $u/d$ quark mass going down to the physical value. We use improved NRQCD for the valence $b$ quarks. Preliminary results show significant improvements over earlier values., Comment: 6 pages, 6 figures, Proceedings of Lattice2014

The charmonium system has several excited states below the energy threshold for decay into $D$ and $\bar{D}$ mesons, which can in principle be studied accurately in lattice QCD. Studies that include many states in the spectrum have typically only been done at one value of the lattice spacing and with relatively heavy light quarks in the sea. Here we give preliminary results for radial and orbital excitation energies for charmonium from a calculation on 2+1+1 MILC configurations at multiple lattice spacings and including physical values for $u/d$ quark masses. We use the HISQ formulation for $c$ to obtain small discretisation errors and smeared operators to improve excited state overlap., Comment: 7 pages, 4 figures. Presented at the 32nd International Symposium on Lattice Field Theory (Lattice 2014), 23-28 June 2014, New York, USA

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., Comment: 21 pages, 15 figures

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$ quark's $\overline{\mathrm{MS}}$ mass: $\alpha_{\overline{\mathrm{MS}}}(M_Z,n_f=5) = 0.11822(74)$ and $m_c(3\mathrm{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., Comment: 16 pages, 9 figures, 6 tables

We describe a new technique to determine the contribution to the anomalous magnetic moment of the muon coming from the hadronic vacuum polarization using lattice QCD. Our method reconstructs the Adler function, using Pad\'{e} approximants, from its derivatives at $q^2=0$ obtained simply and accurately from time-moments of the vector current-current correlator at zero spatial momentum. We test the method using strange quark correlators on large-volume gluon field configurations that include the effect of up and down (at physical masses), strange and charm quarks in the sea at multiple values of the lattice spacing and multiple volumes and show that 1% accuracy is achievable. For the charm quark contributions we use our previously determined moments with up, down and strange quarks in the sea on very fine lattices. We find the (connected) contribution to the anomalous moment from the strange quark vacuum polarization to be $a_\mu^s = 53.41(59) \times 10^{-10}$, and from charm to be $a_\mu^c = 14.42(39)\times 10^{-10}$. These are in good agreement with flavour-separated results from non-lattice methods, given caveats about the comparison. The extension of our method to the light quark contribution and to that from the quark-line disconnected diagram is straightforward., Comment: 7 pages, 4 figures. Updated to version accepted for publication, with additional clarification of the new method developed here

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 \nu$ 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$. }, Comment: 5 pages, 3 figures, Proceedings of Lattice2013

We determine the rate for $D_s^* \rightarrow D_s \gamma$ for the first time from lattice QCD and include the full effect of $u$, $d$ and $s$ sea quarks. The valence quarks are implemented using the Highly Improved Staggered Quark (HISQ) formalism and we normalise the vector current nonperturbatively. We obtain $M(D_s^*)-M(D_s)$ of 148(4) MeV, in good agreement with experiment. The value of the decay constant of the $D_s^*$ we find to be 274(6) MeV, so that $f_{D_s^*}/f_{D_s} = 1.10(2)$. For the radiative decay we find $\Gamma(D_s^* \rightarrow D_s \gamma)$ = 0.066(26) keV. Given the experimental branching fraction for this decay we predict a total width for the $D_s^*$ of 0.070(28) keV, making this the longest lived charged vector meson., Comment: 5 pages, 3 figures

We determine the complete set of axial and vector form factors for the Ds to {\phi}l{\nu} decay from full lattice QCD for the first time. The valence quarks are implemented using the Highly Improved Staggered Quark action and we normalise the appropriate axial and vector currents fully nonperturbatively. The q^2 and angular distributions we obtain for the differential rate agree well with those from the BaBar experiment and, from the total branching fraction, we obtain Vcs = 1.017(63), in good agreement with that from D to Kl{\nu} semileptonic decay. We also find the mass and decay constant of the {\phi} meson in good agreement with experiment, showing that its decay to K{\bar{K}} (which we do not include here) has at most a small effect. We include an Appendix on nonperturbative renormalisation of the complete set of staggered vector and axial vector bilinears needed for this calculation., Comment: 19 pages, 13 figures

We present a calculation of the hyperfine splittings in bottomonium using lattice Nonrelativistic QCD. The calculation includes spin-dependent relativistic corrections through O(v^6), radiative corrections to the leading spin-magnetic coupling and, for the first time, non-perturbative 4-quark interactions which enter at alpha_s^2 v^3. We also include the effect of u,d,s and c quark vacuum polarisation. Our result for the 1S hyperfine splitting is M(Upsilon,1S) - M(eta_b,1S)= 60.0(6.4) MeV. We find the ratio of 2S to 1S hyperfine splittings (M(Upsilon,2S) - M(eta_b,2S))/ (M(Upsilon,1S) - M(eta_b,1S)) = 0.445(28)., Comment: 6 pages, 3 figures. An error was recently found in the coefficients of the 4-quark operators in the lattice NRQCD action that we use here. We have therefore redone this calculation with the corrected coefficients and give the full updated version of the paper here

Calculations of pseudoscalar decay constants of B, Bs, K and pi mesons with physical light quarks are presented. We use HISQ ensembles that include u,d,s and c sea quarks at three lattice spacings. HISQ is used for the valence light quarks and a radiatively improved NRQCD action for the heavy quarks. The key results are f_{B^+}=0.184(4)$ GeV, f_{B_s}=0.224(4) GeV, f_{B_s}/f_{B^+}=1.217(8), f_{K^+}/f_{pi^+}=1.1916(21), f_{K^+}=155.37(34) MeV, giving a significant improvement over previous results that required chiral extrapolation. We also calculate the Wilson flow scale w_0, finding w_0=0.1715(9) fm., Comment: 7 pages, 8 figures, presented at the 31st International Symposium on Lattice Field Theory

We present a new study of the form factors for D -> K semileptonic decay from lattice QCD that allows us to compare the shape of the vector form factor to experiment and, for the first time, to extract V_cs using results from all experimental q^2 bins. The valence quarks are implemented with the Highly Improved Staggered Quark action on MILC configurations that include u, d and s sea quarks. The scalar and vector currents are nonperturbatively normalised and, using phased boundary conditions, we are able to cover the full q^2 range accessible to experiment. Our result is V_cs = 0.963(5)_{expt}(14)_{lattice}. We also demonstrate that the form factors are insensitive to whether the spectator quark is u/d or s, which has implications for other decay channels., Comment: 5 pages, 6 figures

We determine the decay constants of the pi and K mesons on gluon field configurations from the MILC collaboration including u, d, s and c quarks. We use three values of the lattice spacing and u/d quark masses going down to the physical value. We use the w_0 parameter to fix the relative lattice spacing and f_pi to fix the overall scale. This allows us to obtain a value for f{K^+}/f{pi^+} = 1.1916(21). Comparing to the ratio of experimental leptonic decay rates gives |Vus| = 0.22564(28){Br(K^+)}(20){EM}(40){latt}(5){Vud} and the test of unitarity of the first row of the Cabibbo-Kobayashi-Maskawa matrix: |Vud|^2+|Vus|^2+|Vub|^2 - 1 = 0.00009(51)., Comment: 9 pages, 5 figures. Updated to include more discussion of different checks of our fitting procedures and error analysis. Added a figure on autocorrelations. Published version

We present a determination of the b-quark mass accurate through O(\alpha_s^2) in perturbation theory and including partial contributions at O(\alpha_s^3). Nonperturbative input comes from the calculation of the Upsilon and B_s energies in lattice QCD including the effect of u, d and s sea quarks. We use an improved NRQCD action for the b-quark. This is combined with the heavy quark energy shift in NRQCD determined using a mixed approach of high-beta simulation and automated lattice perturbation theory. Comparison with experiment enables the quark mass to be extracted: in the MS bar scheme we find m_b(m_b) = 4.166(43) GeV., Comment: v2 - corrected some typos and an error in the summary plot