
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.

We develop methods to calculate the electroweak gauge boson contribution to
the effective Higgs potential in the context of composite Higgs models, using
lattice gauge theory. The calculation is analogous to that of the
electromagnetic mass splitting of the pion multiplet in QCD. We discuss
technical details of carrying out this calculation, including modeling of the
momentum and fermionmass dependence of the underlying currentcurrent
correlation function; direct integration of the correlation function over
momentum; and fits based on the minimalhadron approximation. We show results
of a numerical study using valence overlap fermions, carried out in an SU(4)
gauge theory with two flavors of Dirac fermions in the twoindex antisymmetric
representation.

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 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 study an SU(3) gauge theory with Nf=8 degenerate flavors of light fermions
in the fundamental representation. Using the domain wall fermion formulation,
we investigate the light hadron spectrum, chiral condensate and electroweak S
parameter. We consider a range of light fermion masses on two lattice volumes
at a single gauge coupling chosen so that IR scales approximately match those
from our previous studies of the two and sixflavor systems. Our results for
the Nf=8 spectrum suggest spontaneous chiral symmetry breaking, though fits to
the fermion mass dependence of spectral quantities do not strongly disfavor the
hypothesis of massdeformed infrared conformality. Compared to Nf=2 we observe
a significant enhancement of the chiral condensate relative to the symmetry
breaking scale F, similar to the situation for Nf=6. The reduction of the S
parameter, related to parity doubling in the vector and axialvector channels,
is also comparable to our sixflavor results.

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.

This is the report of the Computing Frontier working group on Lattice Field
Theory prepared for the proceedings of the 2013 Community Summer Study
("Snowmass"). We present the future computing needs and plans of the U.S.
lattice gauge theory community and argue that continued support of the U.S.
(and worldwide) latticeQCD effort is essential to fully capitalize on the
enormous investment in the highenergy physics experimental program. We first
summarize the dramatic progress of numerical latticeQCD simulations in the
past decade, with some emphasis on calculations carried out under the auspices
of the U.S. LatticeQCD Collaboration, and describe a broad program of
latticeQCD calculations that will be relevant for future experiments at the
intensity and energy frontiers. We then present details of the computational
hardware and software resources needed to undertake these calculations.

Composite dark matter candidates, which can arise from new stronglycoupled
sectors, are wellmotivated and phenomenologically interesting, particularly in
the context of asymmetric generation of the relic density. In this work, we
employ lattice calculations to study the electromagnetic form factors of
electroweakneutral darkmatter baryons for a threecolor, QCDlike theory with
Nf = 2 and 6 degenerate fermions in the fundamental representation. We
calculate the (connected) charge radius and anomalous magnetic moment, both of
which can play a significant role for direct detection of composite dark
matter. We find minimal Nf dependence in these quantities. We generate
massdependent crosssections for dark matternucleon interactions and use them
in conjunction with experimental results from XENON100, excluding dark matter
candidates of this type with masses below 10 TeV.

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.