
The STEREO experiment is a very short baseline reactor antineutrino
experiment aiming at testing the hypothesis of light sterile neutrinos as an
explanation of the deficit of the observed neutrino interaction rate with
respect to the predicted rate, known as the Reactor Antineutrino Anomaly. The
detector center is located 10 m away from the compact, highly $^{235}$U
enriched core of the research nuclear reactor of the Institut Laue Langevin in
Grenoble, France. This paper describes the STEREO site, the detector components
and associated shielding designed to suppress the external sources of
background which were characterized on site. It reports the performances in
terms of detector response and energy reconstruction.

We study the effects of failure to equilibrate the squared topological charge
$Q^2$ on lattice calculations of pseudoscalar masses and decay constants. The
analysis is based on chiral perturbation theory calculations of the dependence
of these quantities on the QCD vacuum angle $\theta$. For the lightlight
partially quenched case, we rederive the known chiral perturbation theory
results of Aoki and Fukaya, but using the nonperturbativelyvalid chiral theory
worked out by Golterman, Sharpe and Singleton, and by Sharpe and Shoresh. We
then extend these calculations to heavylight mesons. Results when staggered
tasteviolations are important are also presented. The derived $Q^2$ dependence
is compared to that of simulations using the MILC collaboration's ensembles of
lattices with four flavors of HISQ dynamical quarks. We find agreement, albeit
with large statistical errors. These results can be used to correct for the
leading effects of unequilibrated $Q^2$, or to make estimates of the systematic
error coming from the failure to equilibrate $Q^2$. In an appendix, we show
that the partially quenched chiral theory may be extended beyond a lower bound
on valence masses discovered by Sharpe and Shoresh. Subtleties occurring when a
seaquark mass vanishes are discussed in another appendix.

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 calculate the up, down, strange, charm, and bottomquark masses using
the MILC highly improved staggeredquark ensembles with four flavors of
dynamical quarks. We use ensembles at six lattice spacings ranging from
$a\approx0.15~$fm to $0.03~$fm and with both physical and unphysical values of
the two light and the strange seaquark masses. We use a new method based on
heavyquark effective theory (HQET) to extract quark masses from heavylight
pseudoscalar meson masses. Combining our analysis with our separate
determination of ratios of lightquark masses we present masses of the up,
down, strange, charm, and bottom quarks. Our results for the
$\overline{\text{MS}}$renormalized masses are $m_u(2~\text{GeV}) =
2.118(38)~$MeV, $m_d(2~\text{GeV}) = 4.690(54)~$MeV, $m_s(2~\text{GeV}) =
92.52(69)~$MeV, $m_c(3~\text{GeV}) = 984.3(5.6)~$MeV, and $m_c(m_c) =
1273(10)~$MeV, with four active flavors; and $m_b(m_b) = 4197(14)~$MeV with
five active flavors. We also obtain ratios of quark masses $m_c/m_s =
11.784(22)$, $m_b/m_s = 53.93(12)$, and $m_b/m_c = 4.577(8)$. The result for
$m_c$ matches the precision of the most precise calculation to date, and the
other masses and all quoted ratios are the most precise to date. Moreover,
these results are the first with a perturbative accuracy of $\alpha_s^4$. As
byproducts of our method, we obtain the matrix elements of HQET operators with
dimension 4 and 5: $\overline{\Lambda}_\text{MRS}=552(30)~$MeV in the minimal
renormalonsubtracted (MRS) scheme, $\mu_\pi^2 = 0.06(22)~\text{GeV}^2$, and
$\mu_G^2(m_b)=0.38(2)~\text{GeV}^2$. The MRS scheme [Phys. Rev. D97, 034503
(2018), arXiv:1712.04983 [hepph]] is the key new aspect of our method.

We calculate the leptonic decay constants of heavylight pseudoscalar mesons
with charm and bottom quarks in lattice quantum chromodynamics on fourflavor
QCD gaugefield configurations with dynamical $u$, $d$, $s$, and $c$ quarks. We
analyze over twenty isospinsymmetric ensembles with six lattice spacings down
to $a\approx 0.03$~fm and several values of the lightquark mass down to the
physical value $\frac{1}{2}(m_u+m_d)$. We employ the highlyimproved
staggeredquark (HISQ) action for the sea and valence quarks; on the finest
lattice spacings, discretization errors are sufficiently small that we can
calculate the $B$meson decay constants with the HISQ action for the first time
directly at the physical $b$quark mass. We obtain the most precise
determinations todate of the $D$ and $B$meson decay constants and their
ratios, $f_{D^+} = 212.6 (0.5)$~MeV, $f_{D_s} = 249.8(0.4)$~MeV,
$f_{D_s}/f_{D^+} = 1.1749(11)$, $f_{B^+} = 189.4(1.4)$~MeV, $f_{B_s} =
230.7(1.2)$~MeV, $f_{B_s}/f_{B^+} = 1.2180(49)$, where the errors include
statistical and all systematic uncertainties. Our results for the $B$meson
decay constants are three times more precise than the previous best latticeQCD
calculations, and bring the QCD errors in the StandardModel predictions for
the rare leptonic decays $\overline{\mathcal{B}}(B_s \to \mu^+\mu^) = 3.65(11)
\times 10^{9}$, $\overline{\mathcal{B}}(B^0 \to \mu^+\mu^) = 1.00(3) \times
10^{11}$, and $\overline{\mathcal{B}}(B^0 \to
\mu^+\mu^)/\overline{\mathcal{B}}(B_s \to \mu^+\mu^) = 0.00264(7)$ to well
below other sources of uncertainty. As a byproduct of our analysis, we also
update our previously published results for the lightquarkmass ratios and the
scalesetting quantities $f_{p4s}$, $M_{p4s}$, and $R_{p4s}$. We obtain the
most precise latticeQCD determination to date of the ratio $f_{K^+}/f_{\pi^+}
= 1.1950(^{+15}_{22})$~MeV.

We report the status of an ongoing latticeQCD calculation of form factors
for exclusive semileptonic decays of $B$ mesons with both charged currents
($B\to\pi\ell\nu$, $B_s\to K\ell\nu$) and neutral currents
($B\to\pi\ell^+\ell^$, $B\to K\ell^+\ell^$). The results are important for
constraining or revealing physics beyond the Standard Model. This work uses
MILC's (2+1+1)flavor ensembles with the HISQ action for the sea and light
valence quarks and the clover action in the Fermilab interpretation for the $b$
quark. Simulations are carried out at three lattice spacings down to $0.088$
fm, with both physical and unphysical seaquark masses. We present preliminary
results for correlationfunction fits.

We present a progress report on our calculation of the decay constants $f_B$
and $f_{B_s}$ from latticeQCD simulations with highlyimproved staggered
quarks. Simulations are carried out with several heavy valencequark masses on
$(2+1+1)$flavor ensembles that include charm sea quarks. We include data at
six lattice spacings and several light seaquark masses, including an
approximately physicalmass ensemble at all but the smallest lattice spacing,
0.03 fm. This range of parameters provides excellent control of the continuum
extrapolation to zero lattice spacing and of heavyquark discretization errors.
Finally, using the heavyquark effective theory expansion we present a method
of extracting from the same correlation functions the charm and bottomquark
masses as well as some lowenergy constants appearing in the heavyquark
expansion.

We discuss the reduction of errors in the calculation of the form factor
$f_+^{K \pi}(0)$ with HISQ fermions on the $N_f=2+1+1$ MILC configurations from
increased statistics on some key ensembles, new data on ensembles with lattice
spacings down to 0.042 fm and the study of finitevolume effects within
staggered ChPT. We also study the implications for the unitarity of the CKM
matrix in the first row and for current tensions with leptonic determinations
of $\vert V_{us}\vert$.

We report a combined experimental and theoretical study of the candidate
typeII Weyl semimetal MoTe2. Using laserbased angleresolved photoemission we
resolve multiple distinct Fermi arcs on the inequivalent top and bottom (001)
surfaces. All surface states observed experimentally are reproduced by an
electronic structure calculation for the experimental crystal structure that
predicts a topological Weyl semimetal state with 8 typeII Weyl points. We
further use systematic electronic structure calculations simulating different
Weyl point arrangements to discuss the robustness of the identified Weyl
semimetal state and the topological character of Fermi arcs in MoTe2.

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 report on a scale determination with gradientflow techniques on the
$N_f=2+1+1$ highly improved staggered quark ensembles generated by the MILC
Collaboration. The ensembles include four lattice spacings, ranging from
approximately 0.15 to 0.06 fm, and both physical and unphysical values of the
quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ and their treelevel
improvements, $\sqrt{t_{0,{\rm imp}}}$ and $w_{0,{\rm imp}}$, are computed on
each ensemble using Symanzik flow and the cloverleaf definition of the energy
density $E$. Using a combination of continuum chiralperturbation theory and a
Taylorseries ansatz for the latticespacing and strongcoupling dependence,
the results are simultaneously extrapolated to the continuum and interpolated
to physical quark masses. We determine the scales $\sqrt{t_0} =
0.1416({}_{5}^{+8})$ fm and $w_0 = 0.1714({}_{12}^{+15})$ fm, where the
errors are sums, in quadrature, of statistical and all systematic errors. The
precision of $w_0$ and $\sqrt{t_0}$ is comparable to or more precise than the
best previous estimates, respectively. We then find the continuum mass
dependence of $\sqrt{t_0}$ and $w_0$, which will be useful for estimating the
scales of new ensembles. We also estimate the integrated autocorrelation length
of $\langle E(t) \rangle$. For long flow times, the autocorrelation length of
$\langle E \rangle$ appears to be comparable to that of the topological charge.

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 give a progress report on a project aimed at a highprecision calculation
of the decay constants $f_B$ and $f_{B_s}$ from simulations with HISQ heavy and
light valence and sea quarks. Calculations are carried out with several heavy
valencequark masses on ensembles with 2+1+1 flavors of HISQ sea quarks at five
lattice spacings and several light seaquark mass ratios $m_{ud}/m_s$,
including approximately physical seaquark masses. This range of parameters
provides excellent control of the continuum limit and of heavyquark
discretization errors. We present a preliminary error budget with projected
uncertainties of 2.2~MeV and 1.5~MeV for $f_B$ and $f_{B_s}$, respectively.

I describe some of the many connections between lattice QCD and effective
field theories, focusing in particular on chiral effective theory, and, to a
lesser extent, Symanzik effective theory. I first discuss the ways in which
effective theories have enabled and supported lattice QCD calculations.
Particular attention is paid to the inclusion of discretization errors, for a
variety of lattice QCD actions, into chiral effective theory. Several other
examples of the usefulness of chiral perturbation theory, including the
encoding of partial quenching and of twisted boundary conditions, are also
described. In the second part of the talk, I turn to results from lattice QCD
for the low energy constants of the two and threeflavor chiral theories. I
concentrate here on mesonic quantities, but the dependence of the nucleon mass
on the pion mass is also discussed. Finally I describe some recent preliminary
lattice QCD calculations by the MILC Collaboration relating to the threeflavor
chiral limit.

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 an update on our calculation of the shortdistance $D^0$meson
mixing hadronic matrix elements. The analysis is performed on the MILC
collaboration's $N_f=2+1$ asqtad configurations. We use asqtad light valence
quarks and the SheikoleslamiWohlert action with the Fermilab interpretation
for the valence charm quark. SU(3), partially quenched, rooted, staggered
heavymeson chiral perturbation theory is used to extrapolate to the
chiralcontinuum limit. Systematic errors arising from the chiralcontinuum
extrapolation, heavyquark discretization, and quarkmass uncertainties are
folded into the statistical errors from the chiralcontinuum fits with methods
of Bayesian inference. A preliminary error budget for all five operators is
presented.

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 leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$,
and the quarkmass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD
using the experimentally determined value of $f_{\pi^+}$ for normalization. We
use the MILC highly improved staggered quark (HISQ) ensembles with four
dynamical quark flavorsup, down, strange, and charmand with both physical
and unphysical values of the light seaquark masses. The use of physical pions
removes the need for a chiral extrapolation, thereby eliminating a significant
source of uncertainty in previous calculations. Four different lattice spacings
ranging from $a\approx 0.06$ fm to $0.15$ fm are included in the analysis to
control the extrapolation to the continuum limit. Our primary results are
$f_{D^+} = 212.6(0.4)({}^{+1.0}_{1.2})\ \mathrm{MeV}$, $f_{D_s} =
249.0(0.3)({}^{+1.1}_{1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} =
1.1712(10)({}^{+29}_{32})$, where the errors are statistical and total
systematic, respectively. The errors on our results for the charm decay
constants and their ratio are approximately two to four times smaller than
those of the most precise previous lattice calculations. We also obtain
$f_{K^+}/f_{\pi^+} = 1.1956(10)({}^{+26}_{18})$, updating our previous result,
and determine the quarkmass ratios $m_s/m_l = 27.35(5)({}^{+10}_{7})$ and
$m_c/m_s = 11.747(19)({}^{+59}_{43})$. When combined with experimental
measurements of the decay rates, our results lead to precise determinations of
the CKM matrix elements $V_{us} = 0.22487(51) (29)(20)(5)$,
$V_{cd}=0.217(1) (5)(1)$ and $V_{cs}= 1.010(5)(18)(6)$, where the errors
are from this calculation of the decay constants, the uncertainty in the
experimental decay rates, structuredependent electromagnetic corrections, and,
in the case of $V_{us}$, the uncertainty in $V_{ud}$, respectively.

We report on a preliminary scale determination with gradientflow techniques
on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration.
The ensembles include four lattice spacings, ranging from 0.15 to 0.06 fm, and
both physical and unphysical values of the quark masses. The scales
$\sqrt{t_0}/a$ and $w_0/a$ are computed using Symanzik flow and the cloverleaf
definition of $\langle E \rangle$ on each ensemble. Then both scales and the
meson masses $aM_\pi$ and $aM_K$ are adjusted for mistunings in the charm mass.
Using a combination of continuum chiral perturbation theory and a Taylor series
ansatz in the lattice spacing, the results are simultaneously extrapolated to
the continuum and interpolated to physical quark masses. Our preliminary
results are $\sqrt{t_0} = 0.1422(7)$fm and $w_0 = 0.1732(10)$fm. We also find
the continuum massdependence of $w_0$.

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$,
and the quarkmass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD. We
use the MILC highly improved staggered quark (HISQ) ensembles with four
dynamical quark flavors. Our primary results are $f_{D^+} =
212.6(0.4)({}^{+1.0}_{1.2})\ \mathrm{MeV}$, $f_{D_s} =
249.0(0.3)({}^{+1.1}_{1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} =
1.1712(10)({}^{+29}_{32})$, where the errors are statistical and total
systematic, respectively. We also obtain $f_{K^+}/f_{\pi^+} =
1.1956(10)({}^{+26}_{18})$, updating our previous result, and determine the
quarkmass ratios $m_s/m_l = 27.35(5)({}^{+10}_{7})$ and $m_c/m_s =
11.747(19)({}^{+59}_{43})$. When combined with experimental measurements of
the decay rates, our results lead to precise determinations of the CKM matrix
elements $V_{us} = 0.22487(51) (29)(20)(5)$, $V_{cd}=0.217(1) (5)(1)$ and
$V_{cs}= 1.010(5)(18)(6)$, where the errors are from this calculation of the
decay constants, the uncertainty in the experimental decay rates,
structuredependent electromagnetic corrections, and, in the case of
$V_{us}$, the uncertainty in $V_{ud}$, respectively.

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

We calculate the kaon semileptonic form factor $f_+(0)$ from lattice QCD,
working, for the first time, at the physical lightquark masses. We use gauge
configurations generated by the MILC collaboration with $N_f=2+1+1$ flavors of
sea quarks, which incorporate the effects of dynamical charm quarks as well as
those of up, down, and strange. We employ data at three lattice spacings to
extrapolate to the continuum limit. Our result, $f_+(0) = 0.9704(32)$, where
the error is the total statistical plus systematic uncertainty added in
quadrature, is the most precise determination to date. Combining our result
with the latest experimental measurements of $K$ semileptonic decays, one
obtains the CabibboKobayashiMaskawa matrix element
$V_{us}=0.22290(74)(52)$, where the first error is from $f_+(0)$ and the
second one is from experiment. In the firstrow test of
CabibboKobayashiMaskawa unitarity, the error stemming from $V_{us}$ is now
comparable to that from $V_{ud}$.

We present results for the form factor $f_+^{K \pi}(0)$, needed to extract
the CKM matrix element $V_{us}$ from experimental data on semileptonic $K$
decays, on the HISQ $N_f=2+1+1$ MILC configurations. The HISQ action is also
used for the valence sector. The data set used for our final result includes
three different values of the lattice spacing and data at the physical light
quark masses. We discuss the error budget and how this calculation improves on
our previous determination of $f_+^{K \pi}(0)$ on the asqtad $N_f=2+1$ MILC
configurations.

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.