
All latticeQCD calculations of the hadronicvacuumpolarization contribution
to the muon's anomalous magnetic moment todate have been performed with
degenerate up and downquark masses. Here we calculate directly the
strongisospinbreaking 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 latticeQCD results obtained with degenerate
lightquark masses of $\delta a_\mu^{{\rm HVP,} m_u \neq m_d}$= +1.5(7)%, in
agreement with estimates from phenomenology and a recent latticeQCD
calculation with unphysically heavy pions.

We calculate in threeflavor lattice QCD the shortdistance hadronic matrix
elements of all five $\Delta C=2$ fourfermion 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 gaugefield configurations generated with
asqtadimproved 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
matrixelement results, we place bounds on the scale of CPviolating new
physics in $D^0$~mixing, finding lower limits of about 1050$\times 10^3$ TeV
for couplings of $\mathrm{O}(1)$. To enable our results to be employed in more
sophisticated or modelspecific phenomenological studies, we provide the
correlations among our matrixelement results. For convenience, we also present
numerical results in the other commonlyused 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 neutrinonucleus
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
electronnucleus interactions and how each underlying physics process  from
quasielastic 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 calculatefor the first time in threeflavor lattice QCDthe 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 CabibboKobayashiMaskawa (CKM) unitarity triangle. Our
calculation employs gaugefield ensembles generated by the MILC Collaboration
with four lattice spacings and pion masses close to the physical value. We use
the asqtadimproved staggered action for the light valence quarks, and the
Fermilab method for the bottom quark. We use heavylight meson chiral
perturbation theory modified to include latticespacing 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 CKMunitarity expectations by about
2$\sigma$. These results and others from flavorchangingneutral 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 gaugefield 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 seaquark 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 modelindependent $z$ expansion to cover the whole kinematically
allowed range. We present our final formfactor 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 StandardModel 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 latticeQCD 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 12$\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
StandardModel rates with LHCb's measurements yields values for the
CabibboKobayashiMaskawa (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
newphysics 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 latticeQCD 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+1flavor lattice configurations at four lattice spacings
and lightquark masses down to 1/20 of the physical strangequark mass. We
extrapolate the lattice form factors to the continuum using staggered chiral
perturbation theory in the hardpion and SU(2) limits. We employ a
modelindependent $z$ parameterization to extrapolate our lattice form factors
from largerecoil momentum to the full kinematic range. We introduce a new
functional method to propagate information from the chiralcontinuum
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 latticeQCD 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 formfactor 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
formfactor 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 preciselydetermined than from our latticeQCD 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$ flavorchanging
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 StandardModel prediction.

We present the first unquenched latticeQCD 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 gaugefield configurations generated with 2+1 flavors of
asqtadimproved 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 seaquark 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 asqtadimproved
staggered fermions. We extrapolate our results to the physical point using
rooted staggered heavylight meson chiral perturbation theory. We then
parameterize the form factors and extend them to the full kinematic range using
modelindependent 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
formfactor 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 zerorecoil 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 fewMeV accuracy. For the
zerorecoil 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+1flavor asqtad ensembles. We
use the Fermilab method for the $b$ quark. We obtain chiralcontinuum
extrapolations for $E_K$ up to $\sim1.2$ GeV and then extend to the entire
kinematic range with the modelindependent $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
BourrelyCapriniLellouch formalism in a modelindependent way. Our analysis is
still blinded with an unknown offset factor which will be disclosed when we
present the final results.

We present an update of the FermilabMILC 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 SheikoleslamiWohlert
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
seaquark mass), and the lattice spacings, from about 0.12 fm to about 0.045
fm. We extrapolate to the physical lightquark mass and the continuum limit
using heavylight 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 highlyimproved 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 chiralcontinuum 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 chargedHiggs boson. Recently the
BaBar experiment reported the first observation of the exclusive semileptonic
decay B > D tau nu, finding an approximately 2sigma disagreement with the
StandardModel 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 StandardModel 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 1sigma 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 twoHiggs doublet model. Once
again, our result is consistent with, but approximately 1sigma 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 StandardModel 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 asqtadimproved
staggered sea quarks. For the valence lightquarks (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 heavylight meson correlation functions on the MILC
asqtadimproved staggered gauge configurations which include the effects of
three light dynamical sea quarks. We simulate with several values of the light
valence and seaquark 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 heavylight meson
staggered chiral perturbation theory. We renormalize the heavylight axial
current using a mostly nonperturbative method such that only a small correction
to unity must be computed in lattice perturbation theory and higherorder 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) flavorbreaking 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 FermilabMILC 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 heavyquark action was extended to include
dimensionsix 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 heavylight 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 indepth
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 CPviolating 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.