• We calculate in three-flavor lattice QCD the short-distance hadronic matrix elements of all five $\Delta C=2$ four-fermion 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 gauge-field configurations generated with asqtad-improved 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 matrix-element results, we place bounds on the scale of CP-violating new physics in $D^0$~mixing, finding lower limits of about 10--50$\times 10^3$ TeV for couplings of $\mathrm{O}(1)$. To enable our results to be employed in more sophisticated or model-specific phenomenological studies, we provide the correlations among our matrix-element results. For convenience, we also present numerical results in the other commonly-used scheme of Buras, Misiak, and Urban.
  • We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark 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 low-energy 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 heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.
  • We calculate---for the first time in three-flavor lattice QCD---the 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 Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light valence quarks, and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing 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 CKM-unitarity expectations by about 2$\sigma$. These results and others from flavor-changing-neutral 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 gauge-field 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 sea-quark 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 model-independent $z$ expansion to cover the whole kinematically allowed range. We present our final form-factor 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 Standard-Model 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 non-compact 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$ flavor-changing 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 Standard-Model 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 zero-recoil 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 few-MeV accuracy. For the zero-recoil 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+1-flavor asqtad ensembles. We use the Fermilab method for the $b$ quark. We obtain chiral-continuum extrapolations for $E_K$ up to $\sim1.2$ GeV and then extend to the entire kinematic range with the model-independent $z$ expansion.
  • We report on a scale determination with gradient-flow 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 one-loop 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 gauge-field 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 lattice-QCD simulations at the physical point. Our result is the most precise lattice-QCD 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 lowest-lying 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 highly-improved 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 chiral-continuum 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 low-energy 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 leading-order 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 higher-order Wilson low-energy constants.
  • The semileptonic decay channel B -> D tau nu is sensitive to the presence of a scalar current, such as that mediated by a charged-Higgs boson. Recently the BaBar experiment reported the first observation of the exclusive semileptonic decay B -> D tau nu, finding an approximately 2-sigma disagreement with the Standard-Model 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 Standard-Model 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 1-sigma 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 two-Higgs doublet model. Once again, our result is consistent with, but approximately 1-sigma 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 Standard-Model 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 asqtad-improved staggered sea quarks. For the valence light-quarks (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 Wilson-Dirac operator is computed numerically in quenched lattice QCD. We demonstrate that the results given for fixed index of the Wilson-Dirac 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 heavy-light meson correlation functions on the MILC asqtad-improved staggered gauge configurations which include the effects of three light dynamical sea quarks. We simulate with several values of the light valence- and sea-quark 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 heavy-light meson staggered chiral perturbation theory. We renormalize the heavy-light axial current using a mostly nonperturbative method such that only a small correction to unity must be computed in lattice perturbation theory and higher-order 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) flavor-breaking 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 near-zero 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 power-law 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.