Strong-Isospin-Breaking Correction to the Muon Anomalous Magnetic Moment from Lattice QCD at the Physical Point

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_{μ}^{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 δa_{μ}^{HVP,m_{u}≠m_{d}}=+1.5(7)%, in agreement with estimates from phenomenology.

Bs → Kℓv form factors with 2+1 flavors

Using the MILC 2+1 flavor asqtad quark action ensembles, we are calculating the form factors f0 and f+ for the semileptonic Bs → Kℓv decay. A total of six ensembles with lattice spacing from ≈ 0.12 to 0.06 fm are being used. At the coarsest and finest lattice spacings, the light quark mass m’l is one-tenth the strange quark mass m’s. At the intermediate lattice spacing, the ratio m’l/m’s ranges from 0.05 to 0.2. The valence b quark is treated using the Sheikholeslami-Wohlert Wilson-clover action with the Fermilab interpretation. The other valence quarks use the asqtad action. When combined with (future) measurements from the LHCb and Belle II experiments, these calculations will provide an alternate determination of the CKM matrix element |Vub|.

B→πll Form Factors for New Physics Searches from Lattice QCD

The rare decay B→πℓ^{+}ℓ^{-} arises from b→d flavor-changing neutral currents and could be sensitive to physics beyond the standard model. Here, we present the first ab initio QCD calculation of the B→π tensor form factor f_{T}. Together with the vector and scalar form factors f_{+} and f_{0} from our companion work [J. A. Bailey et al., Phys. Rev. D 92, 014024 (2015)], these parametrize the hadronic contribution to B→π semileptonic decays in any extension of the standard model. We obtain the total branching ratio BR(B^{+}→π^{+}μ^{+}μ^{-})=20.4(2.1)×10^{-9} in the standard model, which is the most precise theoretical determination to date, and agrees with the recent measurement from the LHCb experiment [R. Aaij et al., J. High Energy Phys. 12 (2012) 125].

Determination of |V(us)|| from a lattice QCD calculation of the K → πℓν semileptonic form factor with physical quark masses

We calculate the kaon semileptonic form factor f+(0) from lattice QCD, working, for the first time, at the physical light-quark masses. We use gauge configurations generated by the MILC Collaboration with Nf = 2 + 1 + 1 flavors of sea quarks, which incorporate the effects of dynamical charm quarks as well as those of up, down, and strange. We employ data at three lattice spacings to extrapolate to the continuum limit. Our result, f+(0) = 0.9704(32), where the error is the total statistical plus systematic uncertainty added in quadrature, is the most precise determination to date. Combining our result with the latest experimental measurements of K semileptonic decays, one obtains the Cabibbo-Kobayashi-Maskawa matrix element |V(us)| = 0.22290(74)(52), where the first error is from f+(0) and the second one is from experiment. In the first-row test of Cabibbo-Kobayashi-Maskawa unitarity, the error stemming from |V(us)| is now comparable to that from |V(ud)|.

Refining new-physics searches in B→Dτν with lattice QCD

The semileptonic decay channel B→Dτν is sensitive to the presence of a scalar current, such as that mediated by a charged-Higgs boson. Recently, the BABAR experiment reported the first observation of the exclusive semileptonic decay B→Dτ(-)ν, finding an approximately 2σ disagreement with the standard-model prediction for the ratio R(D)=BR(B→Dτν)/BR(B→Dℓν), where ℓ = e,μ. We compute this ratio of branching fractions using hadronic form factors computed in unquenched lattice QCD and obtain R(D)=0.316(12)(7), where the errors are statistical and total systematic, respectively. This result is the first standard-model calculation of R(D) from ab initio full QCD. Its error is smaller than that of previous estimates, primarily due to the reduced uncertainty in the scalar form factor f(0)(q(2)). Our determination of R(D) is approximately 1σ higher than previous estimates and, thus, reduces the tension with experiment. We also compute R(D) in models with electrically charged scalar exchange, such as the type-II two-Higgs-doublet model. Once again, our result is consistent with, but approximately 1σ higher than, previous estimates for phenomenologically relevant values of the scalar coupling in the type-II model. As a by-product of our calculation, we also present the standard-model prediction for the longitudinal-polarization ratio P(L)(D)=0.325(4)(3).

Charmed-meson decay constants in three-flavor lattice QCD

We present the first lattice QCD calculation with realistic sea quark content of the D+-meson decay constant f(D+). We use the MILC Collaboration's publicly available ensembles of lattice gauge fields, which have a quark sea with two flavors (up and down) much lighter than a third (strange). We obtain f(D+)=201+/-3+/-17 MeV, where the errors are statistical and a combination of systematic errors. We also obtain f(Ds)=249+/-3+/-16 MeV for the Ds meson.

Semileptonic decays of d mesons in three-flavor lattice QCD

We present the first three-flavor lattice QCD calculations for D-->pilnu and D-->Klnu semileptonic decays. Simulations are carried out using ensembles of unquenched gauge fields generated by the MILC Collaboration. With an improved staggered action for light quarks, we are able to simulate at light quark masses down to 1/8 of the strange mass. Consequently, the systematic error from the chiral extrapolation is much smaller than in previous calculations with Wilson-type light quarks. Our results for the form factors at q(2)=0 are f(D-->pi)(+)(0)=0.64(3)(6) and f(D-->K)(+)(0)=0.73(3)(7), where the first error is statistical and the second is systematic, added in quadrature. Combining our results with experimental branching ratios, we obtain the Cabibbo-Kobayashi-Maskawa matrix elements |V(cd)|=0.239(10)(24)(20) and |V(cs)|=0.969(39)(94)(24), where the last errors are from experimental uncertainties.

High-precision lattice QCD confronts experiment

The recently developed Symanzik-improved staggered-quark discretization allows unquenched lattice-QCD simulations with much smaller (and more realistic) quark masses than previously possible. To test this formalism, we compare experiment with a variety of nonperturbative calculations in QCD drawn from a restricted set of "gold-plated" quantities. We find agreement to within statistical and systematic errors of 3% or less. We discuss the implications for phenomenology and, in particular, for heavy-quark physics.