R. J. Dowdall

$B \rightarrow π\ell ν$ at zero recoil from lattice QCD with physical $u/d$ quarks

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 ν$ .

B-meson decay constants: a more complete picture from full lattice QCD

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σ$ 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.

Bottomonium hyperfine splittings from lattice NRQCD including radiative and relativistic corrections

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

The size of the pion from full lattice QCD with physical $u$ , $d$ , $s$ and $c$ quarks

We present the first calculation of the electromagnetic form factor of the $π$ 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.

B-meson mixing from full lattice QCD with physical u, d, s and c quarks

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.

High-precision quark masses and QCD coupling from $n_f=4$ lattice QCD

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: $α_{\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.

Strange and charm quark contributions to the anomalous magnetic moment of the muon

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é 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_μ^s = 53.41(59) \times 10^{-10}$ , and from charm to be $a_μ^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.

The $Υ$ and $Υ^{\prime}$ Leptonic Widths, $a_μ^b$ and $m_b$ from full lattice QCD

We determine the decay rate to leptons of the ground-state $Υ$ 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 $Γ(Υ\rightarrow e^+e^-)$ = 1.19(11) keV and $Γ(Υ^{\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 $σ(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_μ^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.

B-meson decay constants from improved lattice NRQCD and physical u, d, s and c sea quarks

We present the first lattice QCD calculation of the decay constants f_B and f_{B_s} with physical light quark masses. We use configurations generated by the MILC collaboration including the effect of u, d, s and c HISQ sea quarks at three lattice spacings and with three u/d quark mass values going down to the physical value. We use improved NRQCD for the valence b quarks. Our results are f_B = 0.186(4) GeV, f_{B_s} = 0.224(5) GeV, f_{B_s}/f_B = 1.205(7) and M_{B_s}-M_B=85(2) MeV, superseding earlier results with NRQCD b quarks. We discuss the implications of our results for the Standard Model rates for B_s -> mu+ mu- and B -> tau nu.

B, Bs, K and pi weak matrix elements with physical light quarks

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.

Bottomonium and B results from full lattice QCD

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 πl ν$ 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_π$ ) accurately and we find good agreement as $M_π \rightarrow 0$ . }

Direct determination of strange and light quark condensates from full lattice QCD

We determine the strange and light quark condensates in full lattice QCD for the first time. This is done by direct calculation of the expectation value of the trace of the quark propagator followed by subtraction of the appropriate perturbative contribution to convert to a value for the condensate in the MS-bar scheme at 2 GeV. We use lattice QCD configurations including u, d, s and c quarks in the sea with u/d quark masses going down to the physical value. We find the ratio of the strange to the light quark condensate to be 1.08(16).

Precision tests of the J/psi from full lattice QCD: mass, leptonic width and radiative decay rate to eta_c

We show results from calculations in full lattice QCD of the mass, leptonic width and radiative decay rate to eta_c of the J/psi meson. These provide few % tests of QCD. Another (1.5%) test comes from comparison of time-moments of the vector charmonium correlator with results derived from the experimental values of R(e+e- to hadrons) in the charm region.

The mass of the b-quark from lattice NRQCD and lattice perturbation theory

We present a determination of the b-quark mass accurate through O(α_s^2) in perturbation theory and including partial contributions at O(α_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.

Vus from pi and K decay constants in full lattice QCD with physical u, d, s and c quarks

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

preprint2012arXiv

Direct determination of the strange and light quark condensates from full lattice QCD

We determine the strange quark condensate from lattice QCD for the first time and compare its value to that of the light quark and chiral condensates. The results come from a direct calculation of the expectation value of the trace of the quark propagator followed by subtraction of the appropriate perturbative contribution, derived here, to convert the non-normal-ordered $m\barψψ$ to the $\bar{MS}$ scheme at a fixed scale. This is then a well-defined physical `nonperturbative' condensate that can be used in the Operator Product Expansion of current-current correlators. The perturbative subtraction is calculated through $\mathcal{O}(α_s)$ and estimates of higher order terms are included through fitting results at multiple lattice spacing values. The gluon field configurations used are `second generation' ensembles from the MILC collaboration that include 2+1+1 flavors of sea quarks implemented with the Highly Improved Staggered Quark action and including $u/d$ sea quarks down to physical masses. Our results are : $<\bar{s}{s}>^{\bar{MS}}(2 \mathrm{GeV})= -(290(15) \mathrm{MeV})^3$ , $<\bar{l}{l}>^{\bar{MS}}(2\, \mathrm{GeV})= -(283(2) \mathrm{MeV})^3$ , where $l$ is a light quark with mass equal to the average of the $u$ and $d$ quarks. The strange to light quark condensate ratio is 1.08(16). The light quark condensate is significantly larger than the chiral condensate in line with expectations from chiral analyses. We discuss the implications of these results for other calculations.

preprint2012arXiv

Precise heavy-light meson masses and hyperfine splittings from lattice QCD including charm quarks in the sea

We present improved results for the B and D meson spectrum from lattice QCD including the effect of u/d,s and c quarks in the sea. For the B mesons the Highly Improved Staggered Quark action is used for the sea and light valence quarks and NonRelativistic QCD for the b quark including O(α_s) radiative corrections to many of the Wilson coefficients for the first time. The D mesons use the Highly Improved Staggered Quark action for both valence quarks on the same sea. We find M_{B_s}-M_B=84(2) MeV, M_{B_s}=5.366(8) GeV, M_{B_c}=6.278(9) GeV, M_{D_s}=1.9697(33) GeV, and M_{D_s}-M_{D}=101(3) MeV. Our results for the B meson hyperfine splittings are M_{B^*}-M_{B}=50(3) MeV, M_{B_s^*}-M_{B_s}=52(3) MeV, in good agreement with existing experimental results. This demonstrates that our perturbative improvement of the NRQCD chromo-magnetic coupling works for both heavyonium and heavy-light mesons. We predict M_{B_c^*}-M_{B_c}=54(3) MeV. We also present first results for the radially excited B_c states as well as the orbitally excited scalar B_c0^* and axial vector B_c1 mesons.

preprint2012arXiv

Precision tests of the J/ψ from full lattice QCD: mass, leptonic width and radiative decay rate to η_c

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

B and bottomonium spectroscopy from lattice NRQCD with charm in the sea

We give results for B, Bs, Bc and bottomonium spectroscopy using NRQCD heavy quarks and HISQ valence and sea quarks. Five MILC ensembles of gluon configurations with three values of the lattice spacing and m(light)/m(strange) values down to 0.1 are used that include 2+1+1 flavours of sea quark. Systematic errors in the NRQCD action are improved through the radiative correction of the coefficients of terms at v^4. Improved results for S-wave and P-wave bottomonium states are discussed as well as a prediction for the full D-wave spectrum. Preliminary results for the ground state B meson masses are also presented.

Fast Fits for Lattice QCD Correlators

We illustrate a technique for fitting lattice QCD correlators to sums of exponentials that is significantly faster than traditional fitting methods --- 10--40 times faster for the realistic examples we present. Our examples are drawn from a recent analysis of the Upsilon spectrum, and another recent analysis of the D -> pi semileptonic form factor. For single correlators, we show how to simplify traditional effective-mass analyses.

Prediction of the bottomonium D-wave spectrum from full lattice QCD

We calculate the full spectrum of D-wave states in the Upsilon system in lattice QCD for the first time, using an improved version of NonRelativistic QCD on coarse and fine "second generation" gluon field configurations from the MILC collaboration that include the effect of up, down, strange and charm quarks in the sea. Taking the 2S-1S splitting to set the lattice spacing, we determine the ^3D_2 - 1S splitting to 2.3%, and find agreement with experiment. Our prediction of the fine structure relative to the ^3D_2 gives the ^3D_3 at 10.181(5) GeV and the ^3D_1 at 10.147(6) GeV. We also discuss the overlap of ^3D_1 operators with ^3S_1 states.