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.

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.

The precise measurement of neutrino properties is among the highest priorities in fundamental particle physics, involving many experiments worldwide. Since the experiments rely on the interactions of neutrinos with bound nucleons inside atomic nuclei, the planned advances in the scope and precision of these experiments requires a commensurate effort in the understanding and modeling of the hadronic and nuclear physics of these interactions, which is incorporated as a nuclear model in neutrino event generators. This model is essential to every phase of experimental analyses and its theoretical uncertainties play an important role in interpreting every result. In this White Paper we discuss in detail the impact of neutrino-nucleus interactions, especially the nuclear effects, on the measurement of neutrino properties using the determination of oscillation parameters as a central example. After an Executive Summary and a concise Overview of the issues, we explain how the neutrino event generators work, what can be learned from electron-nucleus interactions and how each underlying physics process - from quasi-elastic to deep inelastic scattering - is understood today. We then emphasize how our understanding must improve to meet the demands of future experiments. With every topic we find that the challenges can be met only with the active support and collaboration among specialists in strong interactions and electroweak physics that include theorists and experimentalists from both the nuclear and high energy physics communities.

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.

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.

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.

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.

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.

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

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.

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.

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.

We update the lattice calculation of the $B\to\pi$ semileptonic form factors, which have important applications to the CKM matrix element $|V_{ub}|$ and the $B\to\pi\ell^+\ell^-$ rare decay. We use MILC asqtad ensembles with $N_f=2+1$ sea quarks and over a range of lattice spacings $a \approx 0.045$--$0.12$ fm. We perform a combined chiral and continuum extrapolation of our lattice data using SU(2) staggered chiral perturbation theory in the hard pion limit. To extend the results for the form factors to the full kinematic range, we take a functional approach to parameterize the form factors using the Bourrely-Caprini-Lellouch formalism in a model-independent way. Our analysis is still blinded with an unknown off-set factor which will be disclosed when we present the final results.

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.

By combining experimentally measured partial branching fractions for the semileptonic decays D to K l nu and D to pi l nu with lattice calculations of the form factors f_+^{D to K}(q^2) and f_+^{D to pi}(q^2), one can extract the CKM matrix elements |V_{cs}| and |V_{cd}|. We are calculating the form factors by using Fermilab charm and asqtad staggered light and strange quarks on 2+1 flavor asqtad staggered ensembles generated by the MILC Collaboration. We vary the light valence quark masses from 0.4m_s to 0.05m_s (m_s is the strange sea-quark mass), and the lattice spacings, from about 0.12 fm to about 0.045 fm. We extrapolate to the physical light-quark mass and the continuum limit using heavy-light meson staggered chiral perturbation theory about the SU(2) and SU(3) limits, compare the resulting (preliminary) form factor shapes with experiment, and discuss our errors.

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.

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

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.

We give an update on simulation results for the decay constants f_B, f_{B_s}, f_D and f_{D_s}. These decay constants are important for precision tests of the standard model, in particular entering as inputs to the global CKM unitarity triangle fit. The results presented here make use of the MILC (2+1)-flavor asqtad ensembles, with heavy quarks incorporated using the clover action with the Fermilab method. Partially quenched, staggered chiral perturbation theory is used to extract the decay constants at the physical point. In addition, we give error projections for a new analysis in progress, based on an extended data set.

We calculate the leptonic decay constants of B_{(s)} and D_{(s)} mesons in lattice QCD using staggered light quarks and Fermilab bottom and charm quarks. We compute the heavy-light meson correlation functions on the MILC asqtad-improved staggered gauge configurations which include the effects of three light dynamical sea quarks. We simulate with several values of the light valence- and sea-quark masses (down to ~m_s/10) and at three lattice spacings (a ~ 0.15, 0.12, and 0.09 fm) and extrapolate to the physical up and down quark masses and the continuum using expressions derived in heavy-light meson staggered chiral perturbation theory. We renormalize the heavy-light axial current using a mostly nonperturbative method such that only a small correction to unity must be computed in lattice perturbation theory and higher-order terms are expected to be small. We obtain f_{B^+} = 196.9(8.9) MeV, f_{B_s} = 242.0(9.5) MeV, f_{D^+} = 218.9(11.3) MeV, f_{D_s} = 260.1(10.8) MeV, and the SU(3) flavor-breaking ratios f_{B_s}/f_{B} = 1.229(26) and f_{D_s}/f_{D} = 1.188(25), where the numbers in parentheses are the total statistical and systematic uncertainties added in quadrature.

The measured partial widths of the semileptonic decays D to K l nu and D to pi l nu can be combined with the form factors calculated on the lattice to extract the CKM matrix elements |V_cs| and |V_cd|. The lattice calculations can be checked by comparing the form factor shapes from the lattice and experiment. We have generated a sizable data set by using heavy clover quarks with the Fermilab interpretation for charm and asqtad staggered light quarks on 2+1 flavor MILC ensembles with lattice spacings of approximately 0.12, 0.09, 0.06, and 0.045 fm. Preliminary fits to staggered chiral perturbation theory suggest that we can reduce the uncertainties in the form factors at zero invariant mass to below 5%.

We give a brief overview and progress report on our lattice QCD calculation of neutral B mixing hadronic matrix elements needed for Standard Model and Beyond the Standard Model physics. Reference [1] contains more details and results.

We report on the status of the Fermilab-MILC calculation of the form factor f_+^{K pi}(0), needed to extract the CKM matrix element |V_{us}| from experimental data on K semileptonic decays. The HISQ formulation is used in the simulations for the valence quarks, while the sea quarks are simulated with the asqtad action (MILC N_f=2+1 configurations). We discuss the general methodology of the calculation, including the use of twisted boundary conditions to get values of the momentum transfer close to zero and the different techniques applied for the correlators fits. We present initial results for lattice spacings a=0.12fm and a=0.09fm, and several choices of the light quark masses.

Recently, the Fermilab heavy-quark action was extended to include dimension-six and -seven operators in order to reduce the discretization errors. In this talk, we present results of the first numerical simulations with this action (the OK action), where we study the masses of the quarkonium and heavy-light systems. We calculate combinations of masses designed to test improvement and compare results obtained with the OK action to their counterparts obtained with the clover action. Our preliminary results show a clear improvement.

One of the major challenges of particle physics has been to gain an in-depth understanding of the role of quark flavor and measurements and theoretical interpretations of their results have advanced tremendously: apart from masses and quantum numbers of flavor particles, there now exist detailed measurements of the characteristics of their interactions allowing stringent tests of Standard Model predictions. Among the most interesting phenomena of flavor physics is the violation of the CP symmetry that has been subtle and difficult to explore. Till early 1990s observations of CP violation were confined to neutral $K$ mesons, but since then a large number of CP-violating processes have been studied in detail in neutral $B$ mesons. In parallel, measurements of the couplings of the heavy quarks and the dynamics for their decays in large samples of $K, D$, and $B$ mesons have been greatly improved in accuracy and the results are being used as probes in the search for deviations from the Standard Model. In the near future, there will be a transition from the current to a new generation of experiments, thus a review of the status of quark flavor physics is timely. This report summarizes the results of the current generation of experiments that is about to be completed and it confronts these results with the theoretical understanding of the field.