Search

We present the first unquenched lattice-QCD calculation of the form factors for the decay $B\rightarrow D^\ast\ell\nu$ at nonzero recoil. Our analysis includes 15 MILC ensembles with $N_f=2+1$ flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from $a\approx 0.15$ fm down to $0.045$ fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence $b$ and $c$ quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element $|V_{cb}|$. We obtain $\left|V_{cb}\right| = (38.40 \pm 0.68_{\textrm{th}} \pm 0.34_{\textrm{exp}} \pm 0.18_{\textrm{EM}})\times 10^{-3}$. The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall $\chi^2\text{/dof} = 126/84$, which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict $R(D^\ast) = 0.265 \pm 0.013$, which confirms the current tension between theory and experiment., Comment: 46 pages, 14 figures. Synthetic data, results and full correlation matrices available in the ancillary files. Version accepted for publication in EPJ C

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 report updates to an ongoing lattice-QCD calculation of the form factors for the semileptonic decays $B \to \pi \ell \nu$, $B_s \to K \ell \nu$, $B \to \pi \ell^+ \ell^-$, and $B \to K \ell^+ \ell^-$. The tree-level decays $B_{(s)} \to \pi (K) \ell \nu$ enable precise determinations of the CKM matrix element $|V_{ub}|$, while the flavor-changing neutral-current interactions $B \to \pi (K) \ell^+ \ell^-$ are sensitive to contributions from new physics. This work uses MILC's (2+1+1)-flavor HISQ ensembles at approximate lattice spacings between $0.057$ and $0.15$ fm, with physical sea-quark masses on four out of the seven ensembles. The valence sector is comprised of a clover $b$ quark (in the Fermilab interpretation) and HISQ light and $s$ quarks. We present preliminary results for the form factors $f_0$, $f_+$, and $f_T$, including studies of systematic errors., Comment: 7 pages; 3 figures; presented at the 37th International Symposium on Lattice Field Theory, 16-22 June 2019, Wuhan, China

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

We discuss preliminary results for the vector form factors $f_+^{\{\pi,K\}}$ at zero-momentum transfer for the decays $D\to\pi\ell\nu$ and $D\to K \ell\nu$ using MILC's $N_f = 2+1+1$ HISQ ensembles at four lattice spacings, $a \approx 0.042, 0.06, 0.09$, and 0.12 fm, and various HISQ quark masses down to the (degenerate) physical light quark mass. We use the kinematic constraint $f_+(q^2)= f_0(q^2)$ at $q^2 = 0$ to determine the vector form factor from our study of the scalar current, which yields $f_0(0)$. Results are extrapolated to the continuum physical point in the framework of hard pion/kaon SU(3) heavy-meson-staggered $\chi$PT and Symanzik effective theory. Our calculation improves upon the precision achieved in existing lattice-QCD calculations of the vector form factors at $q^2=0$. We show the values of the CKM matrix elements $|V_{cs}|$ and $|V_{cd}|$ that we would obtain using our preliminary results for the form factors together with recent experimental results, and discuss the implications of these values for the second row CKM unitarity., Comment: 10 pages, 3 figures, proceeding of The 36th Annual International Symposium on Lattice Field Theory

We use lattice QCD to calculate the form factors $f_+(q^2)$ and $f_0(q^2)$ for the semileptonic decay $B_s\to K\ell\nu$. Our calculation uses six MILC asqtad 2+1 flavor gauge-field ensembles with three lattice spacings. At the smallest and largest lattice spacing the light-quark sea mass is set to 1/10 the strange-quark mass. At the intermediate lattice spacing, we use four values for the light-quark sea mass ranging from 1/5 to 1/20 of the strange-quark mass. We use the asqtad improved staggered action for the light valence quarks, and the clover action with the Fermilab interpolation for the heavy valence bottom quark. We use SU(2) hard-kaon heavy-meson rooted staggered chiral perturbation theory to take the chiral-continuum limit. A functional $z$ expansion is used to extend the form factors to the full kinematic range. We present predictions for the differential decay rate for both $B_s\to K\mu\nu$ and $B_s\to K\tau\nu$. We also present results for the forward-backward asymmetry, the lepton polarization asymmetry, ratios of the scalar and vector form factors for the decays $B_s\to K\ell\nu$ and $B_s\to D_s \ell\nu$. Our results, together with future experimental measurements, can be used to determine the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{ub}|$., Comment: 57 pages, 22 figures, 13 tables

Motivated by the success of the flavour physics programme carried out over the last decade at the Large Hadron Collider (LHC), we characterize in detail the physics potential of its High-Luminosity and High-Energy upgrades in this domain of physics. We document the extraordinary breadth of the HL/HE-LHC programme enabled by a putative Upgrade II of the dedicated flavour physics experiment LHCb and the evolution of the established flavour physics role of the ATLAS and CMS general purpose experiments. We connect the dedicated flavour physics programme to studies of the top quark, Higgs boson, and direct high-$p_T$ searches for new particles and force carriers. We discuss the complementarity of their discovery potential for physics beyond the Standard Model, affirming the necessity to fully exploit the LHC's flavour physics potential throughout its upgrade eras., Comment: Report from Working Group 4 on the Physics of the HL-LHC, and Perspectives at the HE-LHC, 292 pages

We describe a recent lattice-QCD calculation of the leptonic decay constants of heavy-light pseudoscalar mesons containing charm and bottom quarks and of the masses of the up, down, strange, charm, and bottom quarks. Results for these quantities are of the highest precision to date. Calculations use 24 isospin-symmetric ensembles of gauge-field configurations with six different lattice spacings as small as approximately 0.03 fm and several values of the light quark masses down to physical values of the average up- and down-sea-quark masses. We use the highly-improved staggered quark (HISQ) formulation for valence and sea quarks, including the bottom quark. The analysis employs heavy-quark effective theory (HQET). A novel HQET method is used in the determination of the quark masses., Comment: Talk presented CIPANP2018. 11 pages, LaTeX, 8 pdf figures

Using HISQ $N_f=2+1+1$ MILC ensembles with five different values of the lattice spacing, including four ensembles with physical quark masses, we have performed the most precise computation to date of the $K\to\pi\ell\nu$ vector form factor at zero momentum transfer, $f_+^{K^0\pi^-}(0)=0.9696(15)_\text{stat}(12)_\text{syst}$. This is the first calculation that includes the dominant finite-volume effects, as calculated in chiral perturbation theory at next-to-leading order. Our result for the form factor provides a direct determination of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{us}|=0.22333(44)_{f_+(0)}(42)_\text{exp}$, with a theory error that is, for the first time, at the same level as the experimental error. The uncertainty of the semileptonic determination is now similar to that from leptonic decays and the ratio $f_{K^+}/f_{\pi^+}$, which uses $|V_{ud}|$ as input. Our value of $|V_{us}|$ is in tension at the 2--$2.6\sigma$ level both with the determinations from leptonic decays and with the unitarity of the CKM matrix. In the test of CKM unitarity in the first row, the current limiting factor is the error in $|V_{ud}|$, although a recent determination of the nucleus-independent radiative corrections to superallowed nuclear $\beta$ decays could reduce the $|V_{ud}|^2$ uncertainty nearly to that of $|V_{us}|^2$. Alternative unitarity tests using only kaon decays, for which improvements in the theory and experimental inputs are likely in the next few years, reveal similar tensions. As part of our analysis, we calculated the correction to $f_+^{K\pi}(0)$ due to nonequilibrated topological charge at leading order in chiral perturbation theory, for both the full-QCD and the partially-quenched cases. We also obtain the combination of low-energy constants in the chiral effective Lagrangian $[C_{12}^r+C_{34}^r-(L_5^r)^2](M_\rho)=(2.92\pm0.31)\cdot10^{-6}$., Comment: 42 pages and 12 figures. Expanded discussion of fit methodology. Finite volume error increased, conclusions unchanged. Version accepted by Phys. Rev. D

We calculate the up-, down-, strange-, charm-, and bottom-quark masses using the MILC highly improved staggered-quark ensembles with four flavors of dynamical quarks. We use ensembles at six lattice spacings ranging from $a\approx0.15$~fm to $0.03$~fm and with both physical and unphysical values of the two light and the strange sea-quark masses. We use a new method based on heavy-quark effective theory (HQET) to extract quark masses from heavy-light pseudoscalar meson masses. Combining our analysis with our separate determination of ratios of light-quark masses we present masses of the up, down, strange, charm, and bottom quarks. Our results for the $\overline{\text{MS}}$-renormalized masses are $m_u(2~\text{GeV}) = 2.130(41)$~MeV, $m_d(2~\text{GeV}) = 4.675(56)$~MeV, $m_s(2~\text{GeV}) = 92.47(69)$~MeV, $m_c(3~\text{GeV}) = 983.7(5.6)$~MeV, and $m_c(m_c) = 1273(10)$~MeV, with four active flavors; and $m_b(m_b) = 4195(14)$~MeV with five active flavors. We also obtain ratios of quark masses $m_c/m_s = 11.783(25)$, $m_b/m_s = 53.94(12)$, and $m_b/m_c = 4.578(8)$. The result for $m_c$ matches the precision of the most precise calculation to date, and the other masses and all quoted ratios are the most precise to date. Moreover, these results are the first with a perturbative accuracy of $\alpha_s^4$. As byproducts of our method, we obtain the matrix elements of HQET operators with dimension 4 and 5: $\overline{\Lambda}_\text{MRS}=555(31)$~MeV in the minimal renormalon-subtracted (MRS) scheme, $\mu_\pi^2 = 0.05(22)~\text{GeV}^2$, and $\mu_G^2(m_b)=0.38(2)~\text{GeV}^2$. The MRS scheme [Phys. Rev. D97, 034503 (2018), arXiv:1712.04983 [hep-ph]] is the key new aspect of our method., Comment: The published version; 32 pages and 7 figures

We calculate the leptonic decay constants of heavy-light pseudoscalar mesons with charm and bottom quarks in lattice quantum chromodynamics on four-flavor QCD gauge-field configurations with dynamical $u$, $d$, $s$, and $c$ quarks. We analyze over twenty isospin-symmetric ensembles with six lattice spacings down to $a\approx 0.03$~fm and several values of the light-quark mass down to the physical value $\frac{1}{2}(m_u+m_d)$. We employ the highly-improved staggered-quark (HISQ) action for the sea and valence quarks; on the finest lattice spacings, discretization errors are sufficiently small that we can calculate the $B$-meson decay constants with the HISQ action for the first time directly at the physical $b$-quark mass. We obtain the most precise determinations to-date of the $D$- and $B$-meson decay constants and their ratios, $f_{D^+} = 212.7(0.6)$~MeV, $f_{D_s} = 249.9(0.4)$~MeV, $f_{D_s}/f_{D^+} = 1.1749(16)$, $f_{B^+} = 189.4 (1.4)$~MeV, $f_{B_s} = 230.7(1.3)$~MeV, $f_{B_s}/f_{B^+} = 1.2180(47)$, where the errors include statistical and all systematic uncertainties. Our results for the $B$-meson decay constants are three times more precise than the previous best lattice-QCD calculations, and bring the QCD errors in the Standard-Model predictions for the rare leptonic decays $\overline{\mathcal{B}}(B_s \to \mu^+\mu^-) = 3.64(11) \times 10^{-9}$, $\overline{\mathcal{B}}(B^0 \to \mu^+\mu^-) = 1.00(3) \times 10^{-10}$, and $\overline{\mathcal{B}}(B^0 \to \mu^+\mu^-)/\overline{\mathcal{B}}(B_s \to \mu^+\mu^-) = 0.0273(9)$ to well below other sources of uncertainty. As a byproduct of our analysis, we also update our previously published results for the light-quark-mass ratios and the scale-setting quantities $f_{p4s}$, $M_{p4s}$, and $R_{p4s}$. We obtain the most precise lattice-QCD determination to date of the ratio $f_{K^+}/f_{\pi^+} = 1.1950(^{+16}_{-23})$~MeV., Comment: Errors related to the standard model prediction for the rare leptonic decays are fixed in the abstract and Eqs. (7.44), (7.45), and (8.3)

Using the MILC 2+1 flavor asqtad quark action ensembles, we are calculating the form factors $f_0$ and $f_+$ for the semileptonic $B_s \rightarrow K \ell\nu$ decay. A total of six ensembles with lattice spacing from $\approx0.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 $|V_{ub}|$., Comment: 8 pages, 6 figures, to appear in the Proceedings of Lattice 2017, June 18-24, Granada, Spain

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 report the status of an ongoing lattice-QCD calculation of form factors for exclusive semileptonic decays of $B$ mesons with both charged currents ($B\to\pi\ell\nu$, $B_s\to K\ell\nu$) and neutral currents ($B\to\pi\ell^+\ell^-$, $B\to K\ell^+\ell^-$). The results are important for constraining or revealing physics beyond the Standard Model. This work uses MILC's (2+1+1)-flavor ensembles with the HISQ action for the sea and light valence quarks and the clover action in the Fermilab interpretation for the $b$ quark. Simulations are carried out at three lattice spacings down to $0.088$ fm, with both physical and unphysical sea-quark masses. We present preliminary results for correlation-function fits., Comment: 8 pages, 3 figures; presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

We calculate in three-flavor lattice QCD the short-distance hadronic matrix elements of all five $\Delta C=2$ four-fermion operators that contribute to neutral $D$-meson mixing both in and beyond the Standard Model. We use the MILC Collaboration's $N_f = 2+1$ lattice gauge-field configurations generated with asqtad-improved staggered sea quarks. We also employ the asqtad action for the valence light quarks and use the clover action with the Fermilab interpretation for the charm quark. We analyze a large set of ensembles with pions as light as $M_\pi \approx 180$ MeV and lattice spacings as fine as $a\approx 0.045$ fm, thereby enabling good control over the extrapolation to the physical pion mass and continuum limit. We obtain for the matrix elements in the $\overline{\text{MS}}$-NDR scheme using the choice of evanescent operators proposed by Beneke \emph{et al.}, evaluated at 3 GeV, $\langle D^0|\mathcal{O}_i|\bar{D}^0 \rangle = \{0.0805(55)(16), -0.1561(70)(31), 0.0464(31)(9), 0.2747(129)(55), 0.1035(71)(21)\}~\text{GeV}^4$ ($i=1$--5). The errors shown are from statistics and lattice systematics, and the omission of charmed sea quarks, respectively. To illustrate the utility of our matrix-element results, we place bounds on the scale of CP-violating new physics in $D^0$~mixing, finding lower limits of about 10--50$\times 10^3$ TeV for couplings of $\mathrm{O}(1)$. To enable our results to be employed in more sophisticated or model-specific phenomenological studies, we provide the correlations among our matrix-element results. For convenience, we also present numerical results in the other commonly-used scheme of Buras, Misiak, and Urban., Comment: Published version, 42 pages, 18 figures

We present results for neutral D-meson mixing in 2+1-flavor lattice QCD. We compute the matrix elements for all five operators that contribute to D mixing at short distances, including those that only arise beyond the Standard Model. Our results have an uncertainty similar to those of the ETM collaboration (with 2 and with 2+1+1 flavors). This work shares many features with a recent publication on B mixing and with ongoing work on heavy-light decay constants from the Fermilab Lattice and MILC Collaborations., Comment: 6+1 pp., presented at Lattice 2016

We present a progress report on our calculation of the decay constants $f_B$ and $f_{B_s}$ from lattice-QCD simulations with highly-improved staggered quarks. Simulations are carried out with several heavy valence-quark masses on $(2+1+1)$-flavor ensembles that include charm sea quarks. We include data at six lattice spacings and several light sea-quark masses, including an approximately physical-mass ensemble at all but the smallest lattice spacing, 0.03 fm. This range of parameters provides excellent control of the continuum extrapolation to zero lattice spacing and of heavy-quark discretization errors. Finally, using the heavy-quark effective theory expansion we present a method of extracting from the same correlation functions the charm- and bottom-quark masses as well as some low-energy constants appearing in the heavy-quark expansion., Comment: 7 pages, 3 figures, Lattice 2016

We discuss the reduction of errors in the calculation of the form factor $f_+^{K \pi}(0)$ with HISQ fermions on the $N_f=2+1+1$ MILC configurations from increased statistics on some key ensembles, new data on ensembles with lattice spacings down to 0.042 fm and the study of finite-volume effects within staggered ChPT. We also study the implications for the unitarity of the CKM matrix in the first row and for current tensions with leptonic determinations of $\vert V_{us}\vert$., Comment: 7 pages, 2 figures, to appear in the Proceedings of Lattice 2016, The 34th International Symposium on Lattice Field Theory, held 24-30 July 2016, at the University of Southampton, UK. v2: references corrected

We calculate---for the first time in three-flavor lattice QCD---the hadronic matrix elements of all five local operators that contribute to neutral $B^0$- and $B_s$-meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral $B$-meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio $\xi = 1.206(18)(6)$, where the second error stems from the omission of charm sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from $B$ mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light valence quarks, and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral $B$-meson oscillation frequencies to determine the CKM matrix elements $|V_{td}| = 8.00(34)(8) \times 10^{-3}$, $|V_{ts}| = 39.0(1.2)(0.4) \times 10^{-3}$, and $|V_{td}/V_{ts}| = 0.2052(31)(10)$, which differ from CKM-unitarity expectations by about 2$\sigma$. These results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels., Comment: 75 pp, 17 figs. Ver 2 fixes typos; corrects mistakes resulting in slight changes to results, correlation matrices; updates decay constants to agree with recent PDG update; corrects uncertainties for tree-level CKM matrix elements used in comparison, slightly reducing tensions; includes additional analyses that support mostly-nonperturbative matching; expands discussion of isospin-breaking effects

We give a progress report on a project aimed at a high-precision calculation of the decay constants $f_B$ and $f_{B_s}$ from simulations with HISQ heavy and light valence and sea quarks. Calculations are carried out with several heavy valence-quark masses on ensembles with 2+1+1 flavors of HISQ sea quarks at five lattice spacings and several light sea-quark mass ratios $m_{ud}/m_s$, including approximately physical sea-quark masses. This range of parameters provides excellent control of the continuum limit and of heavy-quark discretization errors. We present a preliminary error budget with projected uncertainties of 2.2~MeV and 1.5~MeV for $f_B$ and $f_{B_s}$, respectively., Comment: 7 pages, 6 figures, Lattice 2015

We compute the form factors for the $B \to Kl^+l^-$ semileptonic decay process in lattice QCD using gauge-field ensembles with 2+1 flavors of sea quark, generated by the MILC Collaboration. The ensembles span lattice spacings from 0.12 to 0.045 fm and have multiple sea-quark masses to help control the chiral extrapolation. The asqtad improved staggered action is used for the light valence and sea quarks, and the clover action with the Fermilab interpretation is used for the heavy $b$ quark. We present results for the form factors $f_+(q^2)$, $f_0(q^2)$, and $f_T(q^2)$, where $q^2$ is the momentum transfer, together with a comprehensive examination of systematic errors. Lattice QCD determines the form factors for a limited range of $q^2$, and we use the model-independent $z$ expansion to cover the whole kinematically allowed range. We present our final form-factor results as coefficients of the $z$ expansion and the correlations between them, where the errors on the coefficients include statistical and all systematic uncertainties. We use this complete description of the form factors to test QCD predictions of the form factors at high and low $q^2$. We also compare a Standard-Model calculation of the branching ratio for $B \to Kl^+l^-$ with experimental data., Comment: V2: Fig.7 added. Typos text corrected. Reference added. Version published in Phys. Rev. D

The rare decay $B\to\pi\ell^+\ell^-$ arises from $b\to 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\to\pi$ 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 parameterize the hadronic contribution to $B\to\pi$ semileptonic decays in any extension of the Standard Model. We obtain the total branching ratio ${\text{BR}}(B^+\to\pi^+\mu^+\mu^-)=20.4(2.1)\times10^{-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.$, JHEP 1212, 125 (2012)]. Note added: after this paper was submitted for publication, LHCb announced a new measurement of the differential decay rate for this process [T. Tekampe, talk at DPF 2015], which we now compare to the shape and normalization of the Standard-Model prediction., Comment: V3: Corrected errors in results for Standard-Model differential and total decay rates in abstract, Fig. 3, Table IV, and outlook. Added new preliminary LHCb data to Fig. 3 and brief discussion after outlook. Replaced outdated correlation matrix in Table III with correct final version. Other minor wording changes and references added. 7 pages, 4 tables, 3 figures

We present a lattice-QCD calculation of the $B\to\pi\ell\nu$ semileptonic form factors and a new determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad 2+1-flavor lattice configurations at four lattice spacings and light-quark masses down to 1/20 of the physical strange-quark mass. We extrapolate the lattice form factors to the continuum using staggered chiral perturbation theory in the hard-pion and SU(2) limits. We employ a model-independent $z$ parameterization to extrapolate our lattice form factors from large-recoil momentum to the full kinematic range. We introduce a new functional method to propagate information from the chiral-continuum extrapolation to the $z$ expansion. We present our results together with a complete systematic error budget, including a covariance matrix to enable the combination of our form factors with other lattice-QCD and experimental results. To obtain $|V_{ub}|$, we simultaneously fit the experimental data for the $B\to\pi\ell\nu$ differential decay rate obtained by the BaBar and Belle collaborations together with our lattice form-factor results. We find $|V_{ub}|=(3.72\pm 0.16)\times 10^{-3}$ where the error is from the combined fit to lattice plus experiments and includes all sources of uncertainty. Our form-factor results bring the QCD error on $|V_{ub}|$ to the same level as the experimental error. We also provide results for the $B\to\pi\ell\nu$ vector and scalar form factors obtained from the combined lattice and experiment fit, which are more precisely-determined than from our lattice-QCD calculation alone. These results can be used in other phenomenological applications and to test other approaches to QCD., Comment: 63 pages, 48 figures; v2: minor changes in Sec. IV, Table X, modified Fig.14,16, results unchanged

We present the first unquenched lattice-QCD calculation of the hadronic form factors for the exclusive decay $\overline{B} \rightarrow D \ell \overline{\nu}$ at nonzero recoil. We carry out numerical simulations on fourteen ensembles of gauge-field configurations generated with 2+1 flavors of asqtad-improved staggered sea quarks. The ensembles encompass a wide range of lattice spacings (approximately 0.045 to 0.12 fm) and ratios of light (up and down) to strange sea-quark masses ranging from 0.05 to 0.4. For the $b$ and $c$ valence quarks we use improved Wilson fermions with the Fermilab interpretation, while for the light valence quarks we use asqtad-improved staggered fermions. We extrapolate our results to the physical point using rooted staggered heavy-light meson chiral perturbation theory. We then parameterize the form factors and extend them to the full kinematic range using model-independent functions based on analyticity and unitarity. We present our final results for $f_+(q^2)$ and $f_0(q^2)$, including statistical and systematic errors, as coefficients of a series in the variable $z$ and the covariance matrix between these coefficients. We then fit the lattice form-factor data jointly with the experimentally measured differential decay rate from BaBar to determine the CKM matrix element, $|V_{cb}|=(39.6 \pm 1.7_{\rm QCD+exp} \pm 0.2_{\rm QED})\times 10^{-3}$. As a byproduct of the joint fit we obtain the form factors with improved precision at large recoil. Finally, we use them to update our calculation of the ratio $R(D)$ in the Standard Model, which yields $R(D) = 0.299(11)$., Comment: 47 pages, 32 figures

We present an update on our calculation of the short-distance $D^0$-meson mixing hadronic matrix elements. The analysis is performed on the MILC collaboration's $N_f=2+1$ asqtad configurations. We use asqtad light valence quarks and the Sheikoleslami-Wohlert action with the Fermilab interpretation for the valence charm quark. SU(3), partially quenched, rooted, staggered heavy-meson chiral perturbation theory is used to extrapolate to the chiral-continuum limit. Systematic errors arising from the chiral-continuum extrapolation, heavy-quark discretization, and quark-mass uncertainties are folded into the statistical errors from the chiral-continuum fits with methods of Bayesian inference. A preliminary error budget for all five operators is presented., Comment: 7 pages, 1 figure, LATTICE2014 proceedings

We compute the $B\to\pi\ell\nu$ semileptonic form factors and update the determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad ensembles with $N_f=2+1$ sea quarks at four different lattice spacings in the range $a \approx 0.045$~fm to $0.12$~fm. The lattice form factors are extrapolated to the continuum limit using SU(2) staggered chiral perturbation theory in the hard pion limit, followed by an extrapolation in $q^2$ to the full kinematic range using a functional $z$-parameterization. The extrapolation is combined with the experimental measurements of the partial branching fraction to extract $|V_{ub}|$. Our preliminary result is $|V_{ub}|=(3.72\pm 0.14)\times 10^{-3}$, where the error reflects both the lattice and experimental uncertainties, which are now on par with each other., Comment: 7 pages, 4 figures; talk presented at Lattice 2014. the 32nd International Symposium on Lattice Field Theory, 23-28 June, 2014, Columbia University New York, NY

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. We also obtain $f_{K^+}/f_{\pi^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively., Comment: 7 pages, 1 figure. Proceedings of the 32nd International Symposium on Lattice Field Theory; 23-28 June, 2014, Columbia University, New York

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD using the experimentally determined value of $f_{\pi^+}$ for normalization. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors---up, down, strange, and charm---and with both physical and unphysical values of the light sea-quark masses. The use of physical pions removes the need for a chiral extrapolation, thereby eliminating a significant source of uncertainty in previous calculations. Four different lattice spacings ranging from $a\approx 0.06$ fm to $0.15$ fm are included in the analysis to control the extrapolation to the continuum limit. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. The errors on our results for the charm decay constants and their ratio are approximately two to four times smaller than those of the most precise previous lattice calculations. We also obtain $f_{K^+}/f_{\pi^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively., Comment: v2: minor clarifications and additions; version published in Phys. Rev. D. (73 pages, 26 figures.)

We present a study of the $D$ and $B$ leptonic decay constants on the MILC $N_f=2+1$ asqtad gauge ensembles using asqtad-improved staggered light quarks and clover heavy quarks in the Fermilab interpretation. Our previous analysis \cite{Bazavov:2011aa} computed the decay constants at lattice spacings $a \approx 0.14, 0.11$ and $0.083$ fm. We have extended the simulations to finer $a \approx 0.058$ and $0.043$ fm lattice spacings, and have also increased statistics; this allows us to address many important sources of uncertainty. Technical advances include a two-step two-point fit procedure, better tuning of the heavy quark masses and a better determination of the axial-vector current matching. The present analysis remains blinded, so here we focus on the improvements and their predicted impact on the error budget compared to the prior analysis., Comment: LATTICE 2013; added missing .bbl file

We compute the zero-recoil form factor for the semileptonic decay $\bar{B}^0\to D^{*+}\ell^-\bar{\nu}$ (and modes related by isospin and charge conjugation) using lattice QCD with three flavors of sea quarks. We use an improved staggered action for the light valence and sea quarks (the MILC \asqtad\ configurations), and the Fermilab action for the heavy quarks. Our calculations incorporate higher statistics, finer lattice spacings, and lighter quark masses than our 2008 work. As a byproduct of tuning the new data set, we obtain the $D_s$ and $B_s$ hyperfine splittings with few-MeV accuracy. For the zero-recoil form factor, we obtain $\mathcal{F}(1)=0.906(4)(12)$, where the first error is statistical and the second is the sum in quadrature of all systematic errors. With the latest HFAG average of experimental results and a cautious treatment of QED effects, we find $|V_{cb}| = (39.04 \pm 0.49_\text{expt} \pm 0.53_\text{QCD} \pm 0.19_\text{QED})\times10^{-3}$. The QCD error is now commensurate with the experimental error., Comment: 53 pages, 12 figures; expanded discussion of correlator fits, typos corrected, conforms to version published in PRD

We calculate the form factors for the semileptonic decays $B_s\to K\ell\nu$ and $B\to K\ell\ell$ with lattice QCD. We work at several lattice spacings and a range of light quark masses, using the MILC 2+1-flavor asqtad ensembles. We use the Fermilab method for the $b$ quark. We obtain chiral-continuum extrapolations for $E_K$ up to $\sim1.2$ GeV and then extend to the entire kinematic range with the model-independent $z$ expansion., Comment: 7. pp, 6 figs. Presented at the 31st International Symposium on Lattice Field Theory - LATTICE 2013, July 29 - August 3, 2013, Mainz, Germany

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 $N_f=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}|$., Comment: 6 pages, 2 figures; version published in PRL

We update our determinations of $f_{D^+}$, $f_{D_s}$, $f_K$, and quark mass ratios from simulations with four flavors of HISQ dynamical quarks. The availability of ensembles with light quarks near their physical mass means that we can extract physical results with only small corrections for valence- and sea-quark mass mistunings instead of a chiral extrapolation. The adjusted valence-quark masses and lattice spacings may be determined from an ensemble-by-ensemble analysis, and the results for the quark mass ratios then extrapolated to the continuum limit. Our central values of the charmed meson decay constants, however, come from an alternative analysis, which uses staggered chiral perturbation theory for the heavy-light mesons, and allows us to incorporate data at unphysical quark masses where statistical errors are often smaller. A jackknife analysis propagated through all of these steps takes account of the correlations among all the quantities used in the analysis. Systematic errors from the finite spatial size and EM effects are estimated by varying the parameters in the analysis, and systematic errors from the assumptions in the continuum extrapolation are estimated from the spread of values from different extrapolations., Comment: presented at Lattice 2013, Mainz, Germany, July 29 - August 3, 2013. 14 pages, 7 figures

The determination of $\vert V_{us}\vert$ from kaon semileptonic decays requires the value of the form factor $f_+(q^2=0)$, which can be calculated precisely on the lattice. We provide the one-loop partially quenched staggered chiral perturbation theory expressions that may be employed to analyze staggered simulations of $f_+(q^2)$ with three light flavors. We consider both the case of a mixed action, where the valence and sea sectors have different staggered actions, and the standard case where these actions are the same. The momentum transfer $q^2$ of the form factor is allowed to have an arbitrary value. We give results for the generic situation where the $u$, $d$, and $s$ quark masses are all different, $N_f=1+1+1$, and for the isospin limit, $N_f=2+1$. The expression we obtain for $f_+(q^2)$ is independent of the mass of the (valence) spectator quark. In the limit of vanishing lattice spacing, our results reduce to the one-loop continuum partially quenched expression for $f_+(q^2)$, which has not previously been reported in the literature for the $N_f=1+1+1$ case. Our expressions have already been used in staggered lattice analyses of $f_+(0)$, and should prove useful in future calculations as well., Comment: 33 pages, 5 figures; v2: some referencing changed

We present results for the form factor $f_+^{K \pi}(0)$, needed to extract the CKM matrix element $|V_{us}|$ from experimental data on semileptonic $K$ decays, on the HISQ $N_f=2+1+1$ MILC configurations. The HISQ action is also used for the valence sector. The data set used for our final result includes three different values of the lattice spacing and data at the physical light quark masses. We discuss the error budget and how this calculation improves on our previous determination of $f_+^{K \pi}(0)$ on the asqtad $N_f=2+1$ MILC configurations., Comment: 7 pages, 1 figure, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany; v2: minor changes

We present the status of our calculation of hadronic matrix elements for $D$- and $B$-meson mixing. We use a large set of the MILC collaboration's $N_f=2+1$ asqtad ensembles, which includes lattice spacings in the range $a\approx0.12$-0.045 fm, and up/down to strange quark mass ratios as low as 0.05. The asqtad action is also employed for the light valence quarks. For the heavy quarks we use the Sheikholeslami-Wohlert action with the Fermilab interpretation. Our calculation covers the complete set of five local operators needed to describe $B$-meson mixing in the Standard Model and Beyond. In the charm sector, our calculation of local mixing matrix elements may be used to constrain new physics models. We present final correlator fit results on the full data set for the $B$-meson mixing project and preliminary fit results for the $D$-meson mixing project., Comment: 7 pages, 4 figures, Lattice 2013 proceedings

We present an update of the Fermilab-MILC Collaboration's calculation of hadronic matrix elements for B^0-\bar{B^0} mixing. This work is a more extended analysis than our recent publication of the SU(3)-breaking ratio xi [arXiv:1205.7013]. We use the asqtad staggered action for light valence quarks in combination with the Fermilab interpretation of the Sheikoleslami-Wohlert action for heavy quarks. The calculations use MILC's 2+1 flavor asqtad ensembles. Ensembles include four lattice spacings from approximately 0.125 fm to 0.045 fm and up/down to strange quark mass ratios as low as 0.05. Our calculation covers the complete set of five operators needed to describe B mixing in the Standard Model and beyond. In addition to an update including a fuller set of analyzed data, we comment on the form of the staggered ChPT extrapolation function., Comment: 7 pages, 2 figures; Proceedings of the 30th International Symposium on Lattice Field Theory, June 24-29, 2012, Cairns, Australia

Using staggered fermions and twisted boundary conditions, we calculate the K meson semileptonic decay vector form factor at zero momentum transfer. The HISQ formulation is used for the valence quarks, while the sea quarks are simulated with the asqtad action (MILC N_f=2+1 configurations). For the chiral and continuum extrapolation we use two-loop continuum CHPT, supplemented by partially quenched staggered CHPT at one loop. Our result is f_+^{K\pi}(0) = 0.9667+-0.0023+-0.0033, where the first error is statistical and the second is the sum in quadrature of the systematic uncertainties. This result is the first N_f=2+1 calculation with two lattice spacings and a controlled continuum extrapolation. It is also the most precise result to date for the vector form factor and, although the central value is larger than previous unquenched lattice calculations, it is compatible with them within errors. Combining our value for f_+^{K\pi}(0) with the latest experimental measurements of K semileptonic decays, we obtain |V_{us}| = 0.2238+-0.0009+-0.0005, where the first error is from f_+^{K\pi}(0) and the second one is experimental. As a byproduct of our calculation, we obtain the combination of low-energy constants [C_{12}^r+C_{34}^r-(L_5^r)^2](M_\rho) = (3.62+-1.00)x10^{-6}., Comment: 28 pages, 6 figures

We present our methods to fit the two point correlators for light, strange, and charmed pseudoscalar meson physics with the highly improved staggered quark (HISQ) action. We make use of the least-squares fit including the full covariance matrix of the correlators and including Gaussian constraints on some parameters. We fit the correlators on a variety of the HISQ ensembles. The lattice spacing ranges from 0.15 fm down to 0.06 fm. The light sea quark mass ranges from 0.2 times the strange quark mass down to the physical light quark mass. The HISQ ensembles also include lattices with different volumes and with unphysical values of the strange quark mass. We use the results from this work to obtain our preliminary results of $f_D$, $f_{D_s}$, $f_{D_s}/f_{D}$, and ratios of quark masses presented in another talk [1]., Comment: Proceedings of Lattice 2012 Int. Symposium

We report on the status of our kaon semileptonic form factor calculations using the highly-improved staggered quark (HISQ) formulation to simulate the valence fermions. We present results for the form factor f_+^{K \pi}(0) on the asqtad N_f=2+1 MILC configurations, discuss the chiral-continuum extrapolation, and give a preliminary estimate of the total error. We also present a more preliminary set of results for the same form factor but with the sea quarks also simulated with the HISQ action; these results include data at the physical light quark masses. The improvements that we expect to achieve with the use of the HISQ configurations and simulations at the physical quark masses are briefly discussed., Comment: Proceedings of Lattice 2012 International Symposium, 7 pages, 4 figures

We present preliminary results for light, strange and charmed pseudoscalar meson physics from simulations using four flavors of dynamical quarks with the highly improved staggered quark (HISQ) action. These simulations include lattice spacings ranging from 0.15 to 0.06 fm, and sea-quark masses both above and at their physical value. The major results are charm meson decay constants f_D, f_{D_s} and f_{D_s}/f_D and ratios of quark masses. This talk will focus on our procedures for finding the decay constants on each ensemble, the continuum extrapolation, and estimates of systematic error., Comment: Proceedings of Lattice 2012 Int. Symposium

The semileptonic decay channel B -> D tau nu 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 tau nu, finding an approximately 2-sigma disagreement with the Standard-Model prediction for the ratio R(D)=BR(B->D tau nu)/BR(B->D l nu), where l=e,mu. 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-sigma 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-sigma higher than, previous estimates for phenomenologically relevant values of the scalar coupling in the type II model. As a byproduct of our calculation, we also present the Standard-Model prediction for the longitudinal polarization ratio P_L (D)= 0.325(4)(3)., Comment: 6 pages, 3 figures, 1 table. References and text added. Fig. 3 R(D) in 2HDM II corrected and conclusions modified. Standard-Model R(D) unchanged

We study SU(3)-breaking effects in the neutral B_d-\bar B_d and B_s-\bar B_s systems with unquenched N_f=2+1 lattice QCD. We calculate the relevant matrix elements on the MILC collaboration's gauge configurations with asqtad-improved staggered sea quarks. For the valence light-quarks (u, d, and s) we use the asqtad action, while for b quarks we use the Fermilab action. We obtain \xi=f_{B_s}\sqrt{B_{B_s}}/f_{B_d}\sqrt{B_{B_d}}=1.268+-0.063. We also present results for the ratio of bag parameters B_{B_s}/B_{B_d} and the ratio of CKM matrix elements |V_{td}|/|V_{ts}|. Although we focus on the calculation of \xi, the strategy and techniques described here will be employed in future extended studies of the B mixing parameters \Delta M_{d,s} and \Delta\Gamma_{d,s} in the Standard Model and beyond., Comment: 36 pages, 7 figures

We calculate form-factor ratios between the semileptonic decays \bar{B}->D^+\ell^-\bar{\nu} and \bar{B}_s->D_s^+\ell^-\bar{\nu} with lattice QCD. These ratios are a key theoretical input in a new strategy to determine the fragmentation fractions of the neutral B decays, which are needed for measurements of BR(B^0_s-> \mu^+\mu^-). We use the MILC ensembles of gauge configurations with 2+1 flavors of sea quarks at two lattice spacings of approximately 0.12 fm and 0.09 fm. We use the model-independent z-parametrization to extrapolate our simulation results at small recoil toward maximum recoil. Our results for the form-factor ratios are $f_0^{(s)}(M^2_\pi)/f_0^{(d)}(M^2_K) =1.046(44)_{stat.}(15)_{syst.}$ and $f_0^{(s)}(M^2_\pi)/f_0^{(d)}(M^2_\pi)=1.054(47)_{stat.}(17)_{syst.}$. In contrast to a QCD sum-rule calculation, no significant departure from U-spin (d<->s) symmetry is observed., Comment: 30 pages, 11 figures. Fig. 1 updated. Table II added. Conforms with version published in Physical Review D, except typos fixed, as in the PRD Erratum, in Table V (previously Table IV in arXiv v1). Results unchanged

We report on the status of our lattice-QCD calculation of the hadronic contribution to $B_d^0$ and $B^0_s$ mixing, with 2+1 flavors of dynamical sea quarks. Preliminary results for hadronic mixing matrix elements are given for a basis of five four-quark, dimension-six, $\Delta B=2$ mixing operators that spans the space of all possible hadronic mixing contributions in the Standard Model and beyond. At the intermediate stage of analysis reported on in this work, our errors are competitive with published Standard Model matrix element results. For beyond the Standard Model matrix elements, this is the first unquenched calculation and the first new lattice-QCD calculation in ten years., Comment: 13 pages, 6 figures. Proceedings of the XXIX International Symposium on Lattice Field Theory - Lattice 2011, July 10-16, 2011, Squaw Valley, Lake Tahoe, California. Ver 2: New Figs. 2-4 and Table 4 corrects scripting error. Table 4 now includes results for BJU as well as BBGLN choice of evanescent operators