Search

We calculate the form factors for the $B \to D^*\ellν_\ell$ decay in 2+1 flavor lattice QCD. For all quark flavors, we employ the Möbius domain-wall action, which preserves chiral symmetry to a good precision. Our gauge ensembles are generated at three lattice cutoffs $a^{-1} \sim 2.5$, 3.6 and 4.5 GeV with pion masses as low as $M_π\sim 230$ MeV. The physical lattice size $L$ satisfies the condition $M_πL \geq 4$ to control finite volume effects (FVEs), while we simulate a smaller size at the smallest $M_π$ to directly examine FVEs. The bottom quark masses are chosen in a range from the physical charm quark mass to $0.7 a^{-1}$ to control discretization effects. We extrapolate the form factors to the continuum limit and physical quark masses based on heavy meson chiral perturbation theory at next-to-leading order. Then the recoil parameter dependence is parametrized using a model independent form leading to our estimate of the decay rate ratio between the tau ($\ell = τ$) and light lepton ($\ell = e,μ$) channels $R(D^*) = 0.252(22)$ in the Standard Model. A simultaneous fit with recent data from the Belle experiment yields $|V_{cb}| = 39.19(90)\times 10^{-3}$, which is consistent with previous exclusive determinations, and shows good consistency in the kinematical distribution of the differential decay rate between the lattice and experimental data., 37 pages, 13 figures

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., 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

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 $��$ and $��_b$ and the $��$ 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_�� -M_{��_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_{��_b}=724(12)$ MeV and $f_�� =677.2(9.7)$ MeV, giving $��(��\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 $��(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., 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 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)_{��_{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 ��$ 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 $��$ and $e$ in the final state to be $R_{��/e}=0.9779(2)_{\text{latt}}(50)_{\mathrm{EM}}$ in the Standard Model, with the uncertainty dominated by that from electromagnetic corrections., 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 calculate the $B \to D^{(*)} \ell ��$ form factors in 2+1 flavor relativistic lattice QCD by employing the M��bius domain-wall action for all quark flavors. Our simulations are carried out at lattice cut-offs $a^{-1} \sim 2.5$, 3.6 and 4.5 GeV with the bottom quark masses up to 0.7 $a^{-1}$ to control discretization effects. We extrapolate the form factors to the continuum limit and physical quark masses, and discuss systematic uncertainties of the form factors., 9 pages, 4 figures, talk presented at the 38th International Symposium on Lattice Field Theory (LATTICE2021), 26th-30th July, 2021, Zoom/Gather@MIT, USA

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 ��$ 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., 9 pages, 7 figures. Proceedings of the 38th International Symposium on Lattice Field Theory - LATTICE2021

We analyze the inclusive semileptonic $B \to X_c \ell��$ structure functions in 2+1-flavor lattice QCD. The M��bius domain-wall fermion action is used for light, strange, charm and bottom quarks. The structure function receives contributions from various exclusive modes, including the dominant S-wave states $D^{(*)}_s$ as well as the P-wave states $D_s^{**}$. We can identify them in the lattice data, from which we put some constraints on the $B_s \to D_s^{**}\ell��$ form factors., 7 pages, 3 figures, talk presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China

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/��}$ and including the effect of the $c$ quark's electric charge through quenched QED. We obtain $M_{J/��}-M_{��_c}$ (connected) = 120.3(1.1) MeV and interpret the difference with experiment as the impact on $M_{��_c}$ of its decay to gluons, missing from the lattice calculation. This allows us to determine $��M_{��_c}^{\mathrm{annihiln}}$ =+7.3(1.2) MeV, giving its value for the first time. Our result of $f_{J/��}=$ 0.4104(17) GeV, gives $��(J/��\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_��^c = 14.638(47) \times 10^{-10}$, which is 2.5$��$ 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_��^c$ includes our determination of the effect of QED on this quantity, $��a_��^c = 0.0313(28) \times 10^{-10}$., 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 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}$., 18 pages, 9 figures

We report on our calculation of the B \to D^(*) \ell ��form factors in 2+1 flavor lattice QCD. The M��bius domain-wall action is employed for light, strange, charm and bottom quarks. At lattice cutoffs 1/a \sim 2.4, 3.6 and 4.5 GeV, we simulate bottom quark masses up to 0.7/a to control discretization errors. The pion mass is as low as 230 MeV. We extrapolate the form factors to the continuum limit and physical quark masses, and make a comparison with recent phenomenological analyses., 7 pages, 4 figures, talk presented at the 37th International Symposium on Lattice Field Theory (Lattice2019), 16-22 June 2019, Wuhan, China

We present a lattice QCD determination of the $B_s \to D_s \ell��$ 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 ����_��)/\mathcal{B}(B_s\to D_s \ell ��_{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., 22 pages, 17 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 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 present the first calculation of the pion electromagnetic form factor at physical light quark masses. This form factor parameterises the deviations from the behaviour of a point-like particle when a photon hits the pion. These deviations result from the internal structure of the pion and can thus be calculated in QCD. We use three sets (different lattice spacings) of $n_f = 2+1+1$ lattice configurations generated by the MILC collaboration. The Highly Improved Staggered Quark formalism (HISQ) is used for all of the sea and valence quarks. Using lattice configurations with $u$/$d$ quark masses very close to the physical value is a big advantage, as we avoid the chiral extrapolation. We study the shape of the vector ($f_+$) form factor in the $q^2$ range from $0$ to $-0.15$~GeV$^2$ and extract the mean square radius, $\langle r^2_v\rangle$. The shape of the vector form factor and the resulting radius is compared with experiment. We also discuss the scalar form factor and radius extracted from that, which is not directly accessible to experiment. We have also calculated the contributions from the disconnected diagrams to the scalar form factor at small $q^2$ and discuss their impact on the scalar radius $\langle r^2_s\rangle$., Comment: Talk presented at the 33rd International Symposium on Lattice Field Theory (Lattice 2015), 14-18 July 2015, Kobe, Japan

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 determine the contribution to the anomalous magnetic moment of the muon from the $��^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_��^{\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_��$ and the Standard Model of 3$��$., 14 pages, 10 figures. Discussion of method extended with additional tests and figures added. Typographical errors corrected

The exclusive semileptonic decay $B \rightarrow ��\ell ��$ 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_��$ does hold as $m_�� \rightarrow 0$. We use the Highly Improved Staggered Quark formalism for the light quarks and show that staggered chiral perturbation theory for the $m_��$ dependence is almost identical to continuum chiral perturbation theory for $f_0$, $f_B$ and $f_��$. We also give results for other processes such as $B_s \rightarrow K \ell ��$., 16 pages, 6 figures

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 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 ��l�� 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�� semileptonic decay. We also find the mass and decay constant of the �� 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., 19 pages, 13 figures

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., 5 pages, 6 figures

We calculate the J/�� mass, leptonic width and radiative decay rate to ����_c from lattice QCD including u, d and s quarks in the sea for the first time. We use the Highly Improved Staggered Quark formalism and nonperturbatively normalised vector currents for the leptonic and radiative decay rates. Our results are: M_{J/��} -M_{��_c} = 116.5(3.2) MeV; ��(J/�� to e^+e^-) = 5.48(16) keV; ��(J/�� to ����_c) = 2.49(19) keV. The first two are in good agreement with experiment, with ��(J/�� to e^+e^-) providing a test of a decay matrix element in QCD, independent of CKM uncertainties, to 2%. At the same time results for the time moments of the correlation function can be compared to values from the charm contribution to ��(e^+e^- to hadrons), giving a 1.5% test of QCD. Our results show that an improved experimental error would enable a similarly strong test from ��(J/�� to ����_c)., 20 pages, 16 figures. Added some text plus a figure and two tables comparing the time moments of the charmonium vector correlator to results extracted from experiment for the charm contribution to e^+e^- -> hadrons. This allows a 1.5% test of QCD. Further minor amendments to the text. Version accepted for publication by Physical Review D

This is a follow-up to our earlier work on the energies and radial distributions of heavy-light mesons. The heavy quark is taken to be static (infinitely heavy) and the light quark has a mass about that of the strange quark. We now concentrate on the energies of the excited states with higher angular momentum and with a radial node. A new improvement is the use of hypercubic blocking in the time direction. The calculation is carried out with dynamical fermions on a 16 cubed times 32 lattice with a lattice spacing approximately 0.1 fm generated using a non-perturbatively improved clover action. In nature the closest equivalent of this heavy-light system is the B_s meson, which allows us to compare our lattice calculations to experimental results (where available) or to give a prediction where the excited states, particularly P-wave states, should lie. We pay special attention to the spin-orbit splitting, to see which one of the states (for a given angular momentum L) has the lower energy. An attempt is made to understand these results in terms of the Dirac equation., Comment: 35 pages. v3: Data from two new lattices added. New results in several chapters

This is a follow-up to our earlier work for the energies and the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a mass about that of the strange quark. We study the radial distributions of higher angular momentum states, namely P- and D-wave states, using a "fuzzy" static quark. A new improvement is the use of hypercubic blocking in the time direction, which effectively constrains the heavy quark to move within a 2a hypercube (a is the lattice spacing). The calculation is carried out with dynamical fermions on a 16^3 times 32 lattice with a lattice spacing approximately 0.10 fm generated using the non-perturbatively improved clover action. The configurations were generated by the UKQCD Collaboration using lattice action parameters beta = 5.2, c_SW = 2.0171 and kappa = 0.1350. In nature the closest equivalent of this heavy-light system is the B_s meson. Attempts are now being made to understand these results in terms of the Dirac equation., Presented at XXIV International Symposium on Lattice Field Theory (Lattice 2006, Tucson, Arizona)

This is a follow-up to our earlier work for the energies and the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a mass about that of the strange quark. We now study excited states of these mesons with higher angular momentum and with radial nodes. The calculation is carried out with dynamical fermions on a 16^3 times 32 lattice with a lattice spacing approximately 0.10 fm. The lattice configurations were generated by the UKQCD Collaboration. Attempts are now being made to understand these results in terms of the Dirac equation. In nature the closest equivalent of this heavy-light system is the B_s meson, which allows us to compare our lattice calculations to experimental results (where available) or give a prediction where the P-wave states should lie. We pay particular attention to the spin-orbit splitting, to see which one of the states (for a given angular momentum L) has the lower energy., Comment: Based on a talk presented at the Final Euridice meeting in Kazimierz, Poland, 24-27 August 2006. v2: Table I and Figure 2 updated

This is a follow-up to our earlier work [Phys. Rev. D 65, 014512 (2002); Eur. Phys. J. C 28, 79 (2003); Phys. Rev. D 69, 094505 (2004)] for the energies and the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a mass about that of the strange quark. We study the radial distributions of higher angular momentum states, namely P- and D-wave states. In nature the closest equivalent of this heavy-light system is the B_s meson. The calculation is carried out with dynamical fermions on a 16^3 times 32 lattice with a lattice spacing of about 0.10 fm generated with the non-perturbatively improved clover action. It is shown that several features of the energies and radial distributions are in qualitative agreement with what one expects from a simple one-body Dirac equation interpretation., 6 pages, poster presented at Lattice 2005 (Heavy quarks)

This is a follow-up to our earlier work for the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a mass about that of the strange quark. The calculation is again carried out with dynamical fermions on a 16^3x24 lattice with a lattice spacing of about 0.14 fm. It is shown that several features of the S- and P-wave distributions are in qualitative agreement with what one expects from a simple one-body Dirac equation interpretation., 5 pages, 2 figures, Quark Confinement and the Hadron Spectrum VI, Sardinia, Italy, September, 2004

We determine the energies of the excited states of a heavy-light meson $Q\bar{q}$, with a static heavy quark and light quark with mass approximately that of the strange quark from both quenched lattices and with dynamical fermions. We are able to explore the energies of orbital excitations up to L=3, the spin-orbit splitting up to L=2 and the first radial excitation. These $b \bar{s}$ mesons will be very narrow if their mass is less than 5775 MeV -- the $BK$ threshold. We investigate this in detail and present evidence that the scalar meson (L=1) will be very narrow and that as many as 6 $b \bar{s}$ excited states will have energies close to the $BK$ threshold and should also be relatively narrow., Comment: 17 pages, 6 ps figures

The energies of different angular momentum states of a heavy-light meson were measured on a lattice in Michael and Peisa, Phys. Rev. D 58, 34506 (1998). We have now repeated this study using several different lattices, quenched and unquenched, that have different physical lattice sizes, clover coefficients and quark-gluon couplings. The heavy quark is taken to be infinitely heavy, whereas the light quark mass is approximately that of the strange quark. By interpolating and extrapolating in the light quark mass we can thus compare the lattice results with B and B_s meson experiments. Most interesting is the lowest P-wave B_s state, since it is possible that it lies below the $BK$ threshold and hence is very narrow. Unfortunately, there are no experimental results on P-wave B or B_s mesons available at present. In addition to the energy spectrum, we measured earlier also vector (charge) and scalar (matter) radial distributions of the light quark in the S-wave states of a heavy-light meson on a lattice (Phys. Rev. D 65, 014512 (2002) and Eur. Phys. J. C 28, 79 (2003)). Now we are extending the study of radial distributions to P-wave states., Comment: Poster presented at Lattice2004(heavy), Fermilab, June 21-26, 2004

The energies of different angular momentum states of a heavy-light meson were measured on a lattice in PRD 58, 34506 (1998) . We have now continued this study using several different lattices, quenched and unquenched, that have different physical lattice sizes, clover coefficients, hopping parameters and quark-gluon couplings. The heavy quark is taken to be infinitely heavy, whereas the light quark mass is approximately that of the strange quark. By interpolating in the heavy and light quark masses we can thus compare the lattice results with the B_s meson. Most interesting is the lowest P-wave B_s state, since it is possible that it lies below the BK threshold and hence is very narrow. Unfortunately, there are no experimental results on P-wave B_s mesons available at present. In addition to the energy spectrum, we measured earlier also vector (charge) and scalar (matter) radial distributions of the light quark in the S-wave states of a heavy-light meson on a lattice (PRD 65, 014512 (2002) and Eur. Phys. J. C 28, 79 (2003)). Now we are extending the study of radial distributions to P-wave states., Comment: 5 pages. A talk presented at the Workshop on QCD in Extreme Environments

A knowledge of the radial distributions of quarks inside hadrons could lead to a better understanding of the QCD description of these hadrons and possibly suggest forms for phenomenological models. As a step in this direction, in an earlier work, the charge (vector) and matter (scalar) radial distributions of heavy-light mesons were measured in the quenched approximation on a 16^3x24 lattice with a lattice spacing of 'a' approx. 0.17 fm, and a hopping parameter corresponding to a light quark mass about that of the strange quark. Here several improvements are now made: 1) The configurations are generated using dynamical fermions with a approx 0.14 fm; 2) Many more gauge configurations areincluded; 3) The distributions at many off-axis, in addition to on-axis, points are measured; 4) The data analysis is much more complete. In particular, distributions involving excited states are extracted. The exponential decay of the charge and matter distributions can be described by mesons of mass 0.9 +- 0.1 and 1.5 +- 0.1 GeV respectively -- values that are consistent with those of vector and scalar $q\bar{q}$-states calculated directly with the same lattice parameters., 35 pages, 6 figures, Latex

For a heavy-light meson with a static heavy quark, we can explore the light quark distribution. The charge and matter radial distributions of these heavy-light mesons are measured on a 16^3 * 24 lattice at beta=5.7 and a hopping parameter corresponding to a light quark mass about that of the strange quark. Both distributions can be well fitted up to 4 lattice spacings (r approx 0.7 fm) with the exponential form w_i^2(r), where w_i(r)=A exp(-r/r_i). For the charge(c) and matter(m) distributions r_c approx 0.32(2) fm and r_m approx 0.24(2) fm. We also discuss the normalisation of the total charge and matter integrated over all space, finding 1.30(5) and 0.4(1) respectively., 31 pages including 7 ps figures

In nuclear many-body problems the short-range correlation between two nucleons is well described by the corresponding correlation in the {two}-body problem. Therefore, as a first step in any attempt at an analogous description of many-quark systems, it is necessary to know the two-quark correlation. With this in mind, we study the light quark distribution in a heavy-light meson with a static heavy quark. The charge and matter radial distributions of these heavy-light mesons are measured on a lattice with a light quark mass about that of the strange quark. Both distributions can be well fitted upto r approx 0.7 fm with the exponential form w_i^2(r), where w_i(r)=A exp(-r/r_i). For the charge(c) and matter(m) distributions r_c approx 0.32(2) fm and r_m \approx 0.24(2) fm. We also discuss the normalisation of the total charge (defined to be unity in the continuum limit) and matter integrated over all space, finding 1.30(5) and 0.4(1) respectively for a lattice spacing approx 0.17 fm., Comment: 8 pages, 3 ps figures

A brief review is given of attempts to understand the energies of four-quark systems calculated on a lattice in terms of nuclear-physics-inspired many-body techniques involving interquark potentials. Results are given for the next stage of this study where the wavefunctions of heavy-light mesons are also calculated on a lattice., 3 pages, one figure. Contribution to Confinement IV held in Vienna, 2000

We describe a new technique (presented in arXiv:1403.1778) to determine the contribution to the anomalous magnetic moment (g-2) of the muon coming from the hadronic vacuum polarisation using lattice QCD. Our method uses Pad\'{e} approximants to reconstruct the Adler function from its derivatives at $q^2=0$. These are 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 calculated on MILC Collaboration's $n_f$ = 2+1+1 HISQ ensembles at multiple values of the lattice spacing, multiple volumes and multiple light sea quark masses (including physical pion mass configurations)., Comment: 7 pages, 5 figures, a talk presented at 32nd International Symposium on Lattice Field Theory - LATTICE 2014, June 23 - June 28, 2014, NY, USA

We are very grateful to the Fermilab Directorate and the Fermilab Theoretical Physics Department for their financial and logistical support of the first workshop of the Muon g -2 Theory Initiative (held near Fermilab in June 2017) [123], which was crucial for its success, and indeed for the successful start of the Initiative. Financial support for this workshop was also provided by the Fermilab Distinguished Scholars program, the Universities Research Association through a URA Visiting Scholar award, the Riken Brookhaven Research Center, and the Japan Society for the Promotion of Science under Grant No. KAKEHNHI-17H02906. We thank Shoji Hashimoto, Toru Iijima, Takashi Kaneko, and Shohei Nishida for hosting the HVP workshop at KEK [124] and the KEK Theory Center and the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics for providing logistical and financial support. The HLbL workshop at the University of Connecticut [125] was hosted by the University of Connecticut Physics Department. We also gratefully acknowledge support for the second plenary workshop in Mainz [126] from the Deutsche Forschungsgemeinschaft via the Cluster of Excellence "Precision Physics, Fundamental Interactions and Structure of Matter'' (PRISMA), the Collaborative Research Centre "The low-energy frontier of the Standard Model'' (SFB 1044), as well as the Helmholtz Institute Mainz. And finally, we thank the Institute for Nuclear Theory at the University of Washington for hosting the third plenary workshop [127] and for its kind hospitality and stimulating research environment. This workshop was supported in part by the U.S. Department of Energy, Office of Science, under Award Nos. DE-FG02-00ER41132, DE-SC0020106, and by the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics. This review benefited from discussions with O. Cata, N. Christ, L.Y. Dai, H. Davoudiasl, S. Fayer, S. Ganguly, A. Gasparian, S. Hashimoto, T. Iijima, K. Kampf, D. Kawall, I. Larin, Z. Pagel, M. Petschlies, A. Rebhan, K. Schilcher, K. Shimomura, E. Shintani, D. Steffen, S. Tracz, C. Tu, and T. Yamazaki. The work in this paper was supported by CNRS, by Conacyt (Ciencia Basica 2015) under Grant No. 250628, by CONACyT-Mexico under Grant No. CB2014-22117, by Coordinacion de la Investigacion Cientifica (CIC-UMSNH) under Grant No. 4.10, by Danmarks Frie Forskningsfond under Grant No. 8021-00122B, by Deutsche Forschungsgemeinschaft Collaborative Research Centers CRC 1044, CRC 1044 -204404729, CRC 110, and under Grant No. HI 2048/1-1, Prisma Cluster for Excellence PRISMA+ EXC2118/1, STO/876/6-1, by European Research Council under the European Union's Horizon 2020 research and innovation programme under Grant Agreement Nos. 668679, 757646, 771971-SIMDAMA, 813942, by the European Union H2020-MSCA-COFUND2016 under Grant No. 754510, by the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Grant Agreement No. 843134, by European Union EuroPLEx Grant H2020-MSCA-ITN-2018-813942, by European Union STRONG 2020 project under Grant Agreement No. 824093, by the Excellence Initiative of Aix-Marseille University -A*MIDEX, a French "Investissements d'Avenir" program, through the Chaire d'Excellence program and the OCEVU Laboratoire d'Excellence (ANR-11-LABX-0060), by the Fermilab Distinguished Scholars program, by Fondo SEP-Cinvestav under Grant No. 142, Fundacao para a Ciencia e a Tecnologia under Grant No. SFRH/BPD/109443/2015, by Generalitat de Catalunya under Grant No. 2017SGR1069, by the Helmholtz Association (German Federal Ministry of Education and Research), by the Helmholtz-Institut Mainz, by the Istituto Nazionale di Fisica Nucleare (INFN), by the Isaac Newton Trust, by the Japan Society for the Promotion of Science under Grant Nos. KAKENHI-15H05742, 16K05317, 16K05323, 16K05338, 17H01133, 17H02906, 18H05226, 19K21872, 20K03926, 20K03960, by Junta de Andalucia under Grant No. A-FQM-467-UGR18, by KEK, by Ministerio de Ciencia e Innovacion under Grant No. CICYTFEDER-FPA2017-86989-P, by Ministerio de Industria, Economia y Competitividad under Grant Nos. FPA2016-78220-C3-3-P, FPA2017-86989-P, PGC2018-094857-B-I00, SEV-2016-0588, SEV-2016-0597, by Laboratoires d'Excellence FIRST-TF grants, a French "Investissements d'Avenir'' program, by the Ministry of Science and Higher Education of the Russian Federation under Grant Agreement No. 14.W03.31.0026, by the National Research Foundation of South Africa, by the Natural Sciences and Engineering Research Council of Canada, by the Portuguese Science Foundation (FCT) Investigator Grant IF/00898/2015, by the Romanian Ministry of Education and Research under Grant No. PN19060101, by the Russian Science Foundation under Grant No. RSF 18-12-00128, by the Secretaria d'Universitats i Recerca del Departament d'Economia i Coneixement de la Generalitat de Catalunya under Grant No. 2017 SGR 1069, by the Swedish Research Council under Grant Nos. 2016-05996, 2019-03779, by the Swiss National Science Foundation under Grant Nos. PP00P2_176884, PCEFP2_181117, by The Leverhulme Trust under Grant No. ECF-2019-223, by the UK Science and Technology Facilities Council (STFC) under Grant Nos. ST/N504130/1, ST/P000290/1, ST/P00055X/1, ST/P000630/1, ST/P000711/1, ST/P000746/1, ST/S000879/1, ST/S000925/1, by the U.S. Department of Energy, Office of Science, Office of High Energy Physics under Award Nos. DE-SC0009919, DE-SC0009998, DE-SC0010005, DE-SC0010120, DE-SC0010339, DE-SC0012391, DE-SC0012704, DE-SC0013682, DE-SC0013895, DE-SC0015655, by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics under Award Nos. DE-AC02-05CH11231, DE-FG02-00ER41132, DE-FG02-97ER41020, by the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics, "Incubation Platform for Intensity Frontier", by the U.S. National Science Foundation under Grant Nos. NSF-PHY-1316222, PHY14-14614, PHY17-19626, and PHY19-23131, and by the U.S. National Institute of Standards and Technology (NIST) Precision Measurement Grant Program under Award No. 60NANB16D271. This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics., 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 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 with negligible numerical uncertainty. The electroweak contribution is suppressed by 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 and is due to hadronic vacuum polarization, whereas at 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 and is smaller than the Brookhaven measurement by 3.7. 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., Fermilab Directorate, Fermilab Theoretical Physics Department, Fermilab Distinguished Scholars program, Universities Research Association through a URA Visiting Scholar award, Riken Brookhaven Research Center, Ministry of Education, Culture, Sports, Science and Technology, Japan (MEXT) Japan Society for the Promotion of Science KAKEHNHI-17H02906 KAKENHI-15H05742 16K05317 16K05323 16K05338 17H01133 17H02906 18H05226 19K21872 20K03926 20K03960, High Energy Accelerator Research Organization (KEK), U.S.-Japan Science and Technology Cooperation Program in High Energy Physics, Deutsche Forschungsgemeinschaft via the Cluster of Excellence "Precision Physics, Fundamental Interactions and Structure of Matter'' (PRISMA), Collaborative Research Centre "The low-energy frontier of the Standard Model'' SFB 1044, Helmholtz Institute Mainz, United States Department of Energy (DOE) DE-FG02-00ER41132 DE-SC0020106, Centre National de la Recherche Scientifique (CNRS), Conacyt (Ciencia Basica 2015) 250628, Consejo Nacional de Ciencia y Tecnologia (CONACyT) CB2014-22117, Coordinacion de la Investigacion Cientifica (CIC-UMSNH) 4.10, Danmarks Frie Forskningsfond 8021-00122B, German Research Foundation (DFG) CRC 1044 CRC 1044 -204404729 CRC 110 HI 2048/1-1, Prisma Cluster for Excellence PRISMA+ EXC2118/1 STO/876/6-1, European Research Council (ERC) 668679 757646 771971-SIMDAMA 813942, European Union (EU) 754510, European Union (EU) 843134, European Union EuroPLEx Grant H2020-MSCA-ITN-2018-813942, European Union STRONG 2020 project 824093, French National Research Agency (ANR), OCEVU Laboratoire d'Excellence ANR-11-LABX-0060, Fondo SEP-Cinvestav 142, Portuguese Foundation for Science and Technology SFRH/BPD/109443/2015, Generalitat de Catalunya 2017SGR1069, Helmholtz Association (German Federal Ministry of Education and Research), Istituto Nazionale di Fisica Nucleare (INFN), Isaac Newton Trust, Junta de Andalucia A-FQM-467-UGR18, Instituto de Salud Carlos III Spanish Government CICYTFEDER-FPA2017-86989-P, Ministerio de Industria, Economia y Competitividad FPA2016-78220-C3-3-P FPA2017-86989-P PGC2018-094857-B-I00 SEV-2016-0588 SEV-2016-0597, Laboratoires d'Excellence FIRST-TF grants, Ministry of Science and Higher Education of the Russian Federation 14.W03.31.0026, National Research Foundation - South Africa, Natural Sciences and Engineering Research Council of Canada (NSERC) CGIAR, Portuguese Science Foundation (FCT) Investigator Grant IF/00898/2015, Romanian Ministry of Education and Research PN19060101, Russian Science Foundation (RSF) RSF 18-12-00128, Generalitat de Catalunya 2017 SGR 1069, Swedish Research Council 2016-05996 2019-03779, Swiss National Science Foundation (SNSF) PP00P2_176884 PCEFP2_181117, Leverhulme Trust ECF-2019-223, UK Research & Innovation (UKRI) Science & Technology Facilities Council (STFC) ST/N504130/1 ST/P000290/1 ST/P00055X/1 ST/P000630/1 ST/P000711/1 ST/P000746/1 ST/S000879/1 ST/S000925/1, United States Department of Energy (DOE) DE-SC0009919 DE-SC0009998 DE-SC0010005 DE-SC0010120 DE-SC0010339 DE-SC0012391 DE-SC0012704 DE-SC0013682 DE-SC0013895 DE-SC0015655, United States Department of Energy (DOE) DE-FG02-00ER41132 DE-AC02-05CH11231 DE-FG02-97ER41020 DE-AC02-07CH11359, U.S.-Japan Science and Technology Cooperation Program in High Energy Physics, "Incubation Platform for Intensity Frontier", National Science Foundation (NSF) NSF-PHY-1316222 PHY14-14614 PHY17-19626 PHY19-23131, U.S. National Institute of Standards and Technology (NIST) Precision Measurement Grant 60NANB16D271