Search

An accurate calculation of the leading-order hadronic vacuum polarisation (LOHVP) contribution to the anomalous magnetic moment of the muon ($a_\mu$) is key to determining whether a discrepancy, suggesting new physics, exists between the Standard Model and experimental results. This calculation can be expressed as an integral over Euclidean time of a current-current correlator $G(t)$, where $G(t)$ can be calculated using lattice QCD or, with dispersion relations, from experimental data for $e^+e^-\to\mbox{hadrons}$. The BMW/DMZ collaboration recently presented a hybrid approach in which $G(t)$ is calculated using lattice QCD for most of the contributing $t$ range, but using experimental data for the largest $t$ (lowest energy) region. Here we study the advantages of varying the position $t=t_1$ separating lattice QCD from data-driven contributions. The total LOHVP contribution should be independent of $t_1$, providing both a test of the experimental input and the robustness of the hybrid approach. We use this criterion and a correlated fit to show that Fermilab/HPQCD/MILC lattice QCD results from 2019 strongly favour the CMD-3 cross-section data for $e^+e^-\to\pi^+\pi^-$ over a combination of earlier experimental results for this channel. Further, the resulting total LOHVP contribution obtained is consistent with the result obtained by BMW/DMZ, and supports the scenario in which there is no significant discrepancy between the experimental value for $a_\mu$ and that expected in the Standard Model. We then discuss how improved lattice results in this hybrid approach could provide a more accurate total LOHVP across a wider range of $t_1$ values with an uncertainty that is smaller than that from either lattice QCD or data-driven approaches on their own., Comment: 9 pages, 5 figures. Minor changes to text, added references. Version accepted by Phys. Rev. D

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

We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction., Comment: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

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

The current status of the lattice-QCD calculations of the form factors of the $B\to D^\ast\ell\nu$ semileptonic decay is reviewed. Particular emphasis is given to the most mature calculation at non-zero recoil coming from the Fermilab Lattice and MILC collaborations. Blinded, preliminary results for the form factors are shown, including a preliminary, but detailed error budget. The lattice results seem to favor a large slope at small recoil, in contrast to the latest untagged results coming from the Belle collaboration. A comprehensive comparison between the latest BGL $z$ expansions of Belle, Babar, the lattice and a joint BGL fit including lattice and Belle data is presented, and a roadmap to improve the current calculation is discussed. The current implications for $V_{cb}$ and $R(D^\ast)$ are discussed., Comment: 8 pages, 12 figures. Proceedings of the 2019 Conference on Flavor Physics and CP Violation (FPCP2019). v2 fixes the contour plots of the coefficients of the BGL expansion for BaBar results (b1/b0 and c1/c0)

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

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 present a lattice calculation of the electromagnetic (EM) effects on the masses of light pseudoscalar mesons. The simulations employ 2+1 dynamical flavors of asqtad QCD quarks, and quenched photons. Lattice spacings vary from $\approx 0.12$ fm to $\approx 0.045$ fm. We compute the quantity $\epsilon$, which parameterizes the corrections to Dashen's theorem for the $K^+$-$K^0$ EM mass splitting, as well as $\epsilon_{K^0}$, which parameterizes the EM contribution to the mass of the $K^0$ itself. An extension of the nonperturbative EM renormalization scheme introduced by the BMW group is used in separating EM effects from isospin-violating quark mass effects. We correct for leading finite-volume effects in our realization of lattice electrodynamics in chiral perturbation theory, and remaining finite-volume errors are relatively small. While electroquenched effects are under control for $\epsilon$, they are estimated only qualitatively for $\epsilon_{K^0}$, and constitute one of the largest sources of uncertainty for that quantity. We find $\epsilon = 0.78(1)_{\rm stat}({}^{+\phantom{1}8}_{-11})_{\rm syst}$ and $\epsilon_{K^0}=0.035(3)_{\rm stat}(20)_{\rm syst}$. We then use these results on 2+1+1 flavor pure QCD HISQ ensembles and find $m_u/m_d = 0.4529(48)_{\rm stat}( {}_{-\phantom{1}67}^{+150})_{\rm syst}$., Comment: Version published in Phys. Rev. D. Compared to v1, more discussion of nonperturbative EM renormalization scheme, of statistical errors (with 3 added figures), and of choice QED_{TL} in finite volume. 87 pages, 35 figures

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 present preliminary results from our analysis of the form factors for the $\overline{B}\rightarrow D^\ast\ell\overline{\nu}$ decay at non-zero recoil. Our analysis includes 15 MILC asqtad ensembles with $N_f=2+1$ flavors of sea quarks and lattice spacings ranging from $a\approx 0.15$ fm down to $0.045$ fm. The valence light quarks employ the asqtad action, whereas the heavy quarks are treated using the Fermilab action. We conclude with a discussion of future plans and phenomenological implications. When combined with experimental measurements of the decay rate, our calculation will enable a determination of the CKM matrix element $|V_{cb}|$., Comment: 8 pages, 5 figures, Proceedings of the talk presented at the 35th International Symposium on Lattice Field Theory (Lattice2017), Granada, Spain

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

The MILC Collaboration has completed production running of electromagnetic effects on light mesons using asqtad improved staggered quarks. In these calculations, we use quenched photons in the noncompact formalism. We study four lattice spacings from $\approx\!0.12\:$fm to $\approx\!0.045\:$fm. To study finite-volume effects, we used six spatial lattice sizes $L/a=12$, 16, 20, 28, 40, and 48, at $a\!\approx\!0.12\:$fm. We update our preliminary values for the correction to Dashen's theorem ($\epsilon$) and the quark-mass ratio $m_u/m_d$., Comment: 7 pages, 3 figures, to appear in the Proceedings of Lattice 2015, The 33rd International Symposium on Lattice Field Theory, held 14 -18 July 2015, at Kobe International Conference Center, Kobe, Japan. Will be published by Proceedings of Science, pos.sissa.it, PoS(LATTICE2015)

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

For some time, the MILC Collaboration has been studying electromagnetic effects on light mesons. These calculations use fully dynamical QCD, but only quenched photons, which suffices to NLO in XPT. That is, the sea quarks are electrically neutral, while the valence quarks carry charge. For the photons we use the non-compact formalism. We have new results with lattice spacing as small as 0.045 fm and a large range of volumes. We consider how well chiral perturbation theory describes these results and the implications for light quark masses, Comment: Comments: 6 pages, 4 figures. Proceedings of the XXVI IUPAP Conference on Computational Physics (CCP2014), held at Boston University

We study the exclusive semileptonic $B$-meson decays $B\to K(\pi)\ell^+\ell^-$, $B\to K(\pi)\nu\bar\nu$, and $B\to\pi\tau\nu$, computing observables in the Standard model using the recent lattice-QCD results for the underlying form factors from the Fermilab Lattice and MILC Collaborations. These processes provide theoretically clean windows into physics beyond the Standard Model because the hadronic uncertainties are now under good control for suitably binned observables. For example, the resulting partially integrated branching fractions for $B\to\pi\mu^+\mu^-$ and $B\to K\mu^+\mu^-$ outside the charmonium resonance region are 1-2$\sigma$ higher than the LHCb Collaboration's recent measurements, where the theoretical and experimental errors are commensurate. The combined tension is 1.7$\sigma$. Combining the Standard-Model rates with LHCb's measurements yields values for the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{td}|=7.45{(69)}\times10^{-3}$, $|V_{ts}|=35.7(1.5)\times10^{-3}$, and $|V_{td}/V_{ts}|=0.201{(20)}$, which are compatible with the values obtained from neutral $B_{(s)}$-meson oscillations and have competitive uncertainties. Alternatively, taking the CKM matrix elements from unitarity, we constrain new-physics contributions at the electroweak scale. The constraints on the Wilson coefficients ${\rm Re}(C_9)$ and ${\rm Re}(C_{10})$ from $B\to\pi\mu^+\mu^-$ and $B\to K\mu^+\mu^-$ are competitive with those from $B\to K^* \mu^+\mu^-$, and display a 2.0$\sigma$ tension with the Standard Model. Our predictions for $B\to K(\pi)\nu\bar\nu$ and $B\to\pi\tau\nu$ are close to the current experimental limits., Comment: V3: Typos in Eq. (5.13) and text corrected. Reference added. Version published in Phys. Rev. D

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 report on a calculation of the B*Bpi coupling in lattice QCD. The strong matrix element for a B* to Bpi transition is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HMChPT) for B mesons. We carry out our calculation directly at the b-quark mass using a non-perturbatively tuned clover action that controls discretization effects of order pa and (ma)^n for all n. Our analysis is performed on RBC/UKQCD gauge configurations using domain-wall fermions and the Iwasaki gauge action at two lattice spacings of ainverse = 1.729(25) GeV, ainverse = 2.281(28) GeV, and unitary pion masses down to 290 MeV. We achieve good statistical precision and control all systematic uncertainties, giving a final result for the HMChPT coupling g_b = 0.56(3)stat(7)sys in the continuum and at the physical light-quark masses. This is the first calculation performed directly at the physical b-quark mass and lies in the region one would expect from carrying out an interpolation between previous results at the charm mass and at the static point., Comment: Postprint: version accepted for publication

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 report on a scale determination with gradient-flow techniques on the $N_f=2+1+1$ highly improved staggered quark ensembles generated by the MILC Collaboration. The ensembles include four lattice spacings, ranging from approximately 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ and their tree-level improvements, $\sqrt{t_{0,{\rm imp}}}$ and $w_{0,{\rm imp}}$, are computed on each ensemble using Symanzik flow and the cloverleaf definition of the energy density $E$. Using a combination of continuum chiral-perturbation theory and a Taylor-series ansatz for the lattice-spacing and strong-coupling dependence, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. We determine the scales $\sqrt{t_0} = 0.1416({}_{-5}^{+8})$ fm and $w_0 = 0.1714({}_{-12}^{+15})$ fm, where the errors are sums, in quadrature, of statistical and all systematic errors. The precision of $w_0$ and $\sqrt{t_0}$ is comparable to or more precise than the best previous estimates, respectively. We then find the continuum mass dependence of $\sqrt{t_0}$ and $w_0$, which will be useful for estimating the scales of new ensembles. We also estimate the integrated autocorrelation length of $\langle E(t) \rangle$. For long flow times, the autocorrelation length of $\langle E \rangle$ appears to be comparable to that of the topological charge., Comment: 50 pages, 12 pdf figures, 10 tables, v2: Figure 11 and Table 9 updated to include additional results by other groups, v3: official PRD release with revised charm mass dependence, autocorrelation sections, and other miscellaneous changes, v4: fixed one typo

We calculate the $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses as light as $M_\pi \approx 290$ MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) "hard-pion" chiral perturbation theory. We provide complete error budgets for the form factors $f_+(q^2)$ and $f_0(q^2)$ at three momenta that span the $q^2$ range accessible in our numerical simulations. We extrapolate these results to $q^2 = 0$ using a model-independent $z$-parametrization and present our final form factors as the $z$-coefficients and the matrix of correlations between them. Our results agree with other lattice determinations using staggered light quarks and provide important independent cross-checks. Both $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ decays enable a determination of the CKM matrix element $|V_{ub}|$. To illustrate this, we perform a combined $z$-fit of our numerical $B\to\pi\ell\nu$ form-factor data with the experimental branching-fraction measurements leaving the relative normalization as a free parameter; we obtain $|V_{ub}| = 3.61(32) \times 10^{-3}$, where the error includes statistical and systematic uncertainties. This approach can be applied to $B_s\to K \ell\nu$ decay to determine $|V_{ub}|$ once the process has been measured experimentally. Finally, in anticipation of future measurements, we make predictions for $B \to \pi\ell\nu$ and $B_s\to K \ell\nu$ Standard-Model differential branching fractions and forward-backward asymmetries., Comment: V3: References and discussions added. Large tables move to appendix. 39 pages, 22 tables, 23 figures. Version published in Physical Review D

We report on the status of our calculation of the hadronic matrix elements for neutral $B$-meson mixing with asqtad sea and valence light quarks and using the Wilson clover action with the Fermilab interpretation for the $b$ quark. We calculate the matrix elements of all five local operators that contribute to neutral $B$-meson mixing both in and beyond the Standard Model. We use MILC ensembles with $N_f=2+1$ dynamical flavors at four different lattice spacings in the range $a \approx 0.045$--$0.12$~fm, and with light sea-quark masses as low as 0.05 times the physical strange quark mass. We perform a combined chiral-continuum extrapolation including the so-called wrong-spin contributions in simultaneous fits to the matrix elements of the five operators. We present a complete systematic error budget and conclude with an outlook for obtaining final results from this analysis., Comment: 7 pages, 4 figures, presented at the 32nd International Symposium on Lattice Field Theory (Lattice 2014), 23-28 June, 2014, Columbia University, New York, NY

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 report on a preliminary scale determination with gradient-flow techniques on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The ensembles include four lattice spacings, ranging from 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ are computed using Symanzik flow and the cloverleaf definition of $\langle E \rangle$ on each ensemble. Then both scales and the meson masses $aM_\pi$ and $aM_K$ are adjusted for mistunings in the charm mass. Using a combination of continuum chiral perturbation theory and a Taylor series ansatz in the lattice spacing, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. Our preliminary results are $\sqrt{t_0} = 0.1422(7)$fm and $w_0 = 0.1732(10)$fm. We also find the continuum mass-dependence of $w_0$., Comment: 7 pages, 3 pdf figures, 2 tables, presented at Lattice 2014, Columbia University, June 23-28, 2014

We report on the MILC Collaboration calculation of electromagnetic effects on light pseudoscalar mesons. The simulations employ asqtad staggered dynamical quarks in QCD plus quenched photons, with lattice spacings varying from 0.12 to 0.06 fm. Finite volume corrections for the MILC realization of lattice electrodynamics have been calculated in chiral perturbation theory and applied to the lattice data. These corrections differ from those calculated by Hayakawa and Uno because our treatment of zero modes differs from theirs. Updated results for the corrections to "Dashen's theorem" are presented., Comment: 7 pages, 2 figures. Presented at Lattice 2014, Columbia University, June 23-28, 2014

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 present recent results on the QCD equation of state with 2+1+1 flavors of highly improved staggered quarks (HISQ). We focus on three sets of ensembles with temporal extent 6, 8 and 10, that reach up to temperatures of 967, 725 and 580 MeV, respectively. The strange and charm quark masses are tuned to the physical values and the light quarks mass is set to one fifth of the strange. This corresponds to a Goldstone pion of about 300 MeV., Comment: 7 pages, 3 figures, talk presented at the XXXI International Symposium on Lattice Field Theory, July 29 - August 3, 2013, Mainz, Germany

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

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