
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 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 study the nonperturbative formulation of quantum gravity defined via
Euclidean dynamical triangulations (EDT) in an attempt to make contact with
Weinberg's asymptotic safety scenario. We find that a finetuning is necessary
in order to recover semiclassical behavior. Such a finetuning is generally
associated with the breaking of a target symmetry by the lattice regulator; in
this case we argue that the target symmetry is the general coordinate
invariance of the theory. After introducing and finetuning a nontrivial local
measure term, we find no barrier to taking a continuum limit, and we find
evidence that fourdimensional, semiclassical geometries are recovered at long
distance scales in the continuum limit. We also find that the spectral
dimension at short distance scales is consistent with 3/2, a value that could
resolve the tension between asymptotic safety and the holographic entropy
scaling of black holes. We argue that the number of relevant couplings in the
continuum theory is one, once symmetry breaking by the lattice regulator is
accounted for. Such a theory is maximally predictive, with no adjustable
parameters. The cosmological constant in Planck units is the only relevant
parameter, which serves to set the lattice scale. The cosmological constant in
Planck units is of order 1 in the ultraviolet and undergoes renormalization
group running to small values in the infrared. If these findings hold up under
further scrutiny, the lattice may provide a nonperturbative definition of a
renormalizable quantum field theory of general relativity with no adjustable
parameters and a cosmological constant that is naturally small in the infrared.

We study a formulation of lattice gravity defined via Euclidean dynamical
triangulations (EDT). After finetuning a nontrivial local measure term we
find evidence that fourdimensional, semiclassical geometries are recovered at
long distance scales in the continuum limit. Furthermore, we find that the
spectral dimension at short distance scales is consistent with 3/2, a value
that is also observed in the causal dynamical triangulation (CDT) approach to
quantum gravity.

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

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 report on the status of our calculation of the hadronic matrix elements
for neutral $B$meson mixing with asqtad sea and valence light quarks and using
the Wilson clover action with the Fermilab interpretation for the $b$ quark. We
calculate the matrix elements of all five local operators that contribute to
neutral $B$meson mixing both in and beyond the Standard Model. We use MILC
ensembles with $N_f=2+1$ dynamical flavors at four different lattice spacings
in the range $a \approx 0.045$$0.12$~fm, and with light seaquark masses as
low as 0.05 times the physical strange quark mass. We perform a combined
chiralcontinuum extrapolation including the socalled wrongspin contributions
in simultaneous fits to the matrix elements of the five operators. We present a
complete systematic error budget and conclude with an outlook for obtaining
final results from this analysis.

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.

We present a study of the $D$ and $B$ leptonic decay constants on the MILC
$N_f=2+1$ asqtad gauge ensembles using asqtadimproved staggered light quarks
and clover heavy quarks in the Fermilab interpretation. Our previous analysis
\cite{Bazavov:2011aa} computed the decay constants at lattice spacings $a
\approx 0.14, 0.11$ and $0.083$ fm. We have extended the simulations to finer
$a \approx 0.058$ and $0.043$ fm lattice spacings, and have also increased
statistics; this allows us to address many important sources of uncertainty.
Technical advances include a twostep twopoint fit procedure, better tuning of
the heavy quark masses and a better determination of the axialvector current
matching. The present analysis remains blinded, so here we focus on the
improvements and their predicted impact on the error budget compared to the
prior analysis.

We present recent results on the QCD equation of state with 2+1+1 flavors of
highly improved staggered quarks (HISQ). We focus on three sets of ensembles
with temporal extent 6, 8 and 10, that reach up to temperatures of 967, 725 and
580 MeV, respectively. The strange and charm quark masses are tuned to the
physical values and the light quarks mass is set to one fifth of the strange.
This corresponds to a Goldstone pion of about 300 MeV.