
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 review lattice results related to pion, kaon, D and Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark masses,
the form factor f+(0), arising in the semileptonic K > pi transition at zero
momentum transfer, as well as the decay constant ratio fK/fpi and its
consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe
the results obtained on the lattice for some of the lowenergy constants of
SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the
determination of the BK parameter of neutral kaon mixing as well as the
additional four B parameters that arise in theories of physics beyond the
Standard Model. The latter quantities are an addition compared to the previous
review. For the heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle fit. Finally, we review the status of lattice
determinations of the strong coupling constant alpha_s.

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.

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

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 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 report on a scale determination with gradientflow techniques on the $N_f
= 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The lattice
scale $w_0/a$, originally proposed by the BMW collaboration, is computed using
Symanzik flow at four lattice spacings ranging from 0.15 to 0.06 fm. With a
Taylor series ansatz, the results are simultaneously extrapolated to the
continuum and interpolated to physical quark masses. We give a preliminary
determination of the scale $w_0$ in physical units, along with associated
systematic errors, and compare with results from other groups. We also present
a first estimate of autocorrelation lengths as a function of flowtime for these
ensembles.

We present results from our simulations of quantum chromodynamics (QCD) with
four flavors of quarks: u, d, s, and c. These simulations are performed with a
oneloop Symanzik improved gauge action, and the highly improved staggered
quark (HISQ) action. We are generating gauge configurations with four values of
the lattice spacing ranging from 0.06 fm to 0.15 fm, and three values of the
light quark mass, including the value for which the Goldstone pion mass is
equal to the physical pion mass. We discuss simulation algorithms, scale
setting, taste symmetry breaking, and the autocorrelations of various
quantities. We also present results for the topological susceptibility which
demonstrate the improvement of the HISQ configurations relative to those
generated earlier with the asqtad improved staggered action.

A calculation of the ratio of leptonic decay constants f_{K^+}/f_{\pi^+}
makes possible a precise determination of the ratio of CKM matrix elements
V_{us}/V_{ud} in the Standard Model, and places a stringent constraint on
the scale of new physics that would lead to deviations from unitarity in the
first row of the CKM matrix. We compute f_{K^+}/f_{\pi^+} numerically in
unquenched lattice QCD using gaugefield ensembles recently generated that
include four flavors of dynamical quarks: up, down, strange, and charm. We
analyze data at four lattice spacings a ~ 0.06, 0.09, 0.12, and 0.15 fm with
simulated pion masses down to the physical value 135 MeV. We obtain
f_{K^+}/f_{\pi^+} = 1.1947(26)(37), where the errors are statistical and total
systematic, respectively. This is our first physics result from our N_f = 2+1+1
ensembles, and the first calculation of f_{K^+}/f_{\pi^+} from latticeQCD
simulations at the physical point. Our result is the most precise latticeQCD
determination of f_{K^+}/f_{\pi^+}, with an error comparable to the current
world average. When combined with experimental measurements of the leptonic
branching fractions, it leads to a precise determination of V_{us}/V_{ud} =
0.2309(9)(4) where the errors are theoretical and experimental, respectively.

We calculate pseudoscalar masses on gauge configurations containing the
effects of 2+1 flavors of dynamical asqtad quarks and quenched
electromagnetism. The lattice spacings vary from 0.12 to 0.06 fm. The masses
are fit with staggered chiral perturbation theory including NLO electromagnetic
terms. We attempt to extract the fit parameters for the electromagnetic
contributions, while taking into account the finite volume effects, and
extrapolate them to the physical limit.

We investigate $U_A(1)$ breaking above $T_c$ in terms of the Dirac spectrum
on configurations with (2+1)flavors, using the HISQ action. The strange quark
mass is at its physical value. We use several light quark masses corresponding
to the Goldstone pion masses in the range of about 115  230 MeV on lattices
of size 32$^3 \times$8 and 48$^3 \times$8. We calculate the 100 lowestlying
Dirac eigenvalues at temperatures below and above $T_c$. We investigate the
volume dependence of the Dirac eigenvalue density to determine whether there is
a gap around zero, which can appear if $U_A(1)$ symmetry is restored in the
chiral symmetric phase. We also investigate the quark mass dependence of the
Dirac eigenvalue density at zero and check whether there is a linear behavior
that would signal the $U_A(1)$ breaking above $T_c$.

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.

We show how the leading physical and Wilson lowenergy constants associated
with Wilson fermions in lattice gauge theory can be determined individually by
using spectral information of the Wilson Dirac operator with fixed index at
finite volume. The methods are demonstrated in simulations with leadingorder
improved Wilson fermions. In addition to the expected suppression of the
leading term in Wilson chiral perturbation theory we observe a substantial
reduction also of the higherorder Wilson lowenergy constants.

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 analyze topological charge contributions from classical SU(2) center
vortices with shapes of planes and spheres using different topological charge
definitions, namely the center vortex picture of topological charge, a discrete
version of F\~{F} in the plaquette and hypercube definitions and the lattice
index theorem. For the latter the zeromodes of the Dirac operator in the
fundamental and adjoint representations using both the overlap and asqtad
staggered fermion formulations are investigated. We find several problems for
the individual definitions and discuss the discrepancies between the different
topological charge definitions. Our results show that the interpretation of
topological charge in the background of center vortices is rather subtle.

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 present results on the chiral and deconfinement properties of the QCD
transition at finite temperature. Calculations are performed with 2+1 flavors
of quarks using the p4, asqtad and HISQ/tree actions. Lattices with temporal
extent N_tau=6, 8 and 12 are used to understand and control discretization
errors and to reliably extrapolate estimates obtained at finite lattice
spacings to the continuum limit. The chiral transition temperature is defined
in terms of the phase transition in a theory with two massless flavors and
analyzed using O(N) scaling fits to the chiral condensate and susceptibility.
We find consistent estimates from the HISQ/tree and asqtad actions and our main
result is T_c=154 +/ 9 MeV.

The microscopic spectral density of the Hermitian WilsonDirac operator is
computed numerically in quenched lattice QCD. We demonstrate that the results
given for fixed index of the WilsonDirac operator can be matched by the
predictions from Wilson chiral perturbation theory. We test successfully the
finite volume and the mass scaling predicted by Wilson chiral perturbation
theory at fixed lattice spacing.

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.

We report on the behavior of the eigenvalue distribution of the Dirac
operator in (2+1)flavor QCD at finite temperature, using the HISQ action. We
calculate the eigenvalue density at several values of the temperature close to
the pseudocritical temperature. For this study we use gauge field
configurations generated on lattices of size $32^3 \times 8$ with two light
quark masses corresponding to pion masses of about 160 and 115 MeV. We find
that the eigenvalue density below $T_c$ receives large contributions from
nearzero modes which become smaller as the temperature increases or the light
quark mass decreases. Moreover we find no clear evidence for a gap in the
eigenvalue density up to 1.1$T_c$. We also analyze the eigenvalue density near
$T_c$ where it appears to show a powerlaw behavior consistent with what is
expected in the critical region near the second order chiral symmetry restoring
phase transition in the massless limit.

We present an update of our calculations of the decay constants of the D,
D_s, B, and B_s mesons in unquenched 2+1 flavor QCD. We use the MILC library of
improved staggered gauge ensembles at lattice spacings 0.09, 0.12, and 0.15 fm,
clover heavy quarks with the Fermilab normalizations, and improved staggered
light valence quarks.