Search

We present results from the first lattice QCD calculation of the two-pion contributions to the light-quark connected vector-current correlation function obtained from staggered-quark operators. We employ the MILC collaboration's gauge-field ensemble with $2+1+1$ flavors of highly improved staggered sea quarks at a lattice spacing of $a\approx 0.15$ fm with a light sea-quark mass at its physical value. The two-pion contributions allow for a refined determination of the noisy long-distance tail of the vector-current correlation function, which we use to compute the light-quark connected contribution to HVP with improved statistical precision. We compare our results with traditional noise-reduction techniques used in lattice QCD calculations of the light-quark connected HVP, namely the so-called fit and bounding methods. We observe a factor of roughly three improvement in the statistical precision in the determination of the HVP contribution to the muon's anomalous magnetic moment over these approaches. We also lay the group theoretical groundwork for extending this calculation to finer lattice spacings with increased numbers of staggered two-pion taste states.

We report on the ongoing effort of improving the determination of the gradient flow scale on the (2+1+1)-flavor HISQ ensembles generated by the MILC collaboration. We compute the scales $\sqrt{t_0}/a$ and $w_0/a$ with the Wilson and Symanzik flow using three discretizations for the action density: clover, Wilson and tree-level Symanzik improved. For the absolute scale setting, we intend to employ the $\Omega$-baryon mass, but are also using the pion decay constant while the $\Omega$-mass calculations are in progress., Comment: 10 pages, 6 figures, contribution to the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31 - August 4, 2023, Fermilab, USA

A method, known as ``minimal renormalon subtraction'' [Phys. Rev. D 97 (2018) 034503, JHEP 2017 (2017) 62], relates the factorial growth of a perturbative series (in QCD) to the power~$p$ of a power correction $\Lambda^p/Q^p$. ($\Lambda$ is the QCD scale, $Q$ some hard scale.) Here, the derivation is simplified and generalized to any~$p$, more than one such correction, and cases with anomalous dimensions. Strikingly, the well-known factorial growth is seen to emerge already at low or medium orders, as a consequence of constraints on the $Q$ dependence from the renormalization group. The effectiveness of the method is studied with the gluonic energy between a static quark and static antiquark (the ``static energy''). Truncation uncertainties are found to be under control after next-to-leading order, despite the small exponent of the power correction ($p=1$) and associated rapid growth seen in the first four coefficients of the perturbative series., Comment: 25 pp + title page, 8 figures; v2 corrects spelling and grammar typos & conforms with published version

We present a lattice-QCD calculation of the light-quark connected contribution to window observables associated with the leading-order hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, $a_\mu^{\mathrm{HVP,LO}}$. We employ the MILC Collaboration's isospin-symmetric QCD gauge-field ensembles, which contain four flavors of dynamical highly-improved-staggered quarks with four lattice spacings between $a\approx 0.06$-$0.15$~fm and close-to-physical quark masses. We consider several effective-field-theory-based schemes for finite-volume and other lattice corrections and combine the results via Bayesian model averaging to obtain robust estimates of the associated systematic uncertainties. After unblinding, our final results for the intermediate and ``W2'' windows are $a^{ll,{\mathrm W}}_{\mu}(\mathrm{conn.})=206.6(1.0) \times 10^{-10}$ and $a^{ll,\mathrm {W2}}_{\mu}(\mathrm{conn.}) = 100.7(3.2)\times 10^{-10}$, respectively., Comment: v3: Updated to reflect published version, which includes updates to the text in Sections II.B, III.C,D,E. Numerical results unchanged

We present lattice-QCD calculations of the hadronic form factors for the semileptonic decays $D\to\pi\ell\nu$, $D\to K\ell\nu$, and $D_s\to K\ell\nu$. Our calculation uses the highly improved staggered quark (HISQ) action for all valence and sea quarks and includes $N_f=2+1+1$ MILC ensembles with lattice spacings ranging from $a\approx0.12$ fm down to $0.042$ fm. At most lattice spacings, an ensemble with physical-mass light quarks is included. The HISQ action allows all the quarks to be treated with the same relativistic light-quark action, allowing for nonperturbative renormalization using partial conservation of the vector current. We combine our results with experimental measurements of the differential decay rates to determine $|V_{cd}|^{D\to\pi}=0.2238(11)^{\rm Expt}(15)^{\rm QCD}(04)^{\rm EW}(02)^{\rm SIB}[22]^{\rm QED}$ and $|V_{cs}|^{D\to K}=0.9589(23)^{\rm Expt}(40)^{\rm QCD}(15)^{\rm EW}(05)^{\rm SIB}[95]^{\rm QED}$ This result for $|V_{cd}|$ is the most precise to date, with a lattice-QCD error that is, for the first time for the semileptonic extraction, at the same level as the experimental error. Using recent measurements from BES III, we also give the first-ever determination of $|V_{cd}|^{D_s\to K}=0.258(15)^{\rm Expt}(01)^{\rm QCD}[03]^{\rm QED}$ from $D_s\to K \ell\nu$. Our results also furnish new Standard Model calculations of the lepton flavor universality ratios $R^{D\to\pi}=0.98671(17)^{\rm QCD}[500]^{\rm QED}$, $R^{D\to K}=0.97606(16)^{\rm QCD}[500]^{\rm QED}$, and $R^{D_s\to K}=0.98099(10)^{\rm QCD}[500]^{\rm QED}$, which are consistent within $2\sigma$ with experimental measurements. Our extractions of $|V_{cd}|$ and $|V_{cs}|$, when combined with a value for $|V_{cb}|$, provide the most precise test of second-row CKM unitarity, finding agreement with unitarity at the level of one standard deviation., Comment: 92 pages, V2 matches version accepted for publication in PRD. Expanded supplementary material for reconstructing our final results. An implementation of nonlinear shrinkage is also included

This report was prepared on behalf of three Energy Frontier Topical Groups of the Snowmass 2021 Community Planning Exercise. It summarizes the status and implications of studies of strong interactions in high-energy experiments and QCD theory. We emphasize the rich landscape and broad impact of these studies in the decade ahead. Hadronic interactions play a central role in the high-luminosity Large Hadron Collider (LHC) physics program, and strong synergies exist between the (HL-)LHC and planned or proposed experiments at the U.S. Electron-Ion Collider, CERN forward physics experiments, high-intensity facilities, and future TeV-range lepton and hadron colliders. Prospects for precision determinations of the strong coupling and a variety of nonperturbative distribution and fragmentation functions are examined. We also review the potential of envisioned tests of new dynamical regimes of QCD in high-energy and high-density scattering processes with nucleon, ion, and photon initial states. The important role of the high-energy heavy-ion program in studies of nuclear structure and the nuclear medium, and its connections with QCD involving nucleons are summarized. We address ongoing and future theoretical advancements in multi-loop QCD computations, lattice QCD, jet substructure, and event generators. Cross-cutting connections between experimental measurements, theoretical predictions, large-scale data analysis, and high-performance computing are emphasized., Comment: 95 pages (bibliography 30 pages), 28 figures; v.2: minor changes, authors and references added

Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure and spectrum, to serve as a numerical laboratory to reach beyond the Standard Model, or to invent and improve state-of-the-art computational paradigms, the lattice-gauge-theory program is in a prime position to impact the course of developments and enhance discovery potential of a vibrant experimental program in High-Energy Physics over the coming decade. This projection is based on abundant successful results that have emerged using lattice gauge theory over the years: on continued improvement in theoretical frameworks and algorithmic suits; on the forthcoming transition into the exascale era of high-performance computing; and on a skillful, dedicated, and organized community of lattice gauge theorists in the U.S. and worldwide. The prospects of this effort in pushing the frontiers of research in High-Energy Physics have recently been studied within the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021), and the conclusions are summarized in this Topical Report., Comment: 57 pages, 1 figure. Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Topical Group Report for TF05 - Lattice Gauge Theory

Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)., Comment: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

An accurate determination of the leading-order hadronic vacuum polarisation (HVP) contribution to the anomalous magnetic moment of the muon is critical to understanding the size and significance of any discrepancy between the Standard Model prediction and experimental results being obtained by the Muon g-2 experiment at Fermilab. The Standard Model prediction is currently based on a data-driven approach to the HVP using experimental results for $\sigma(e^+e^-\rightarrow\,\mathrm{hadrons})$. Lattice QCD aims to provide a result with similar uncertainty from calculated vector-vector correlation functions, but the growth of statistical and systematic errors in the $u/d$ quark correlation functions at large Euclidean time has made this difficult to achieve. We show that restricting the lattice contributions to a one-sided window $0

We present results for the static energy in ($2+1+1$)-flavor QCD over a wide range of lattice spacings and several quark masses, including the physical quark mass, with ensembles of lattice-gauge-field configurations made available by the MILC Collaboration. We obtain results for the static energy out to distances of nearly $1$~fm, allowing us to perform a simultaneous determination of the scales $r_{1}$ and $r_{0}$, as well as the string tension $\sigma$. For the smallest three lattice spacings we also determine the scale $r_{2}$. Our results for $r_{0}/r_{1}$ and $r_{0}\sqrt{\sigma}$ agree with published ($2+1$)-flavor results. However, our result for $r_{1}/r_{2}$ differs significantly from the value obtained in the ($2+1$)-flavor case, which is most likely due to the effect of the charm quark. We also report results for $r_{0}$, $r_{1}$, and $r_{2}$ in~fm, with the former two being slightly lower than published ($2+1$)-flavor results. We study in detail the effect of the charm quark on the static energy by comparing our results on the finest two lattices with the previously published ($2+1$)-flavor QCD results at similar lattice spacing. We find that for $r > 0.2$~fm our results on the static energy agree with the ($2+1$)-flavor result, implying the decoupling of the charm quark for these distances. For smaller distances, on the other hand, we find that the effect of the dynamical charm quark is noticeable. The lattice results agree well with the two-loop perturbative expression of the static energy incorporating finite charm mass effects. This is the first time that the decoupling of the charm quark is observed and quantitatively analyzed on lattice data of the static energy., Comment: 51 pages, 37 figures, journal version

We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction., Comment: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. This white paper discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators., Comment: 81 pages, contribution to Snowmass 2021

This document provides a comprehensive summary of the state-of-the-art, challenges, and prospects in the experimental and theoretical study of the strong coupling $\alpha_s$. The current status of the seven methods presently used to determine $\alpha_s$ based on: (i) lattice QCD, (ii) hadronic $\tau$ decays, (iii) deep-inelastic scattering and parton distribution functions fits, (iv) electroweak boson decays, hadronic final-states in (v) e+e-, (vi) e-p, and (vii) p-p collisions, and (viii) quarkonia decays and masses, are reviewed. Novel $\alpha_s$ determinations are discussed, as well as the averaging method used to obtain the PDG world-average value at the reference Z boson mass scale, $\alpha_s(m^2_Z)$. Each of the extraction methods proposed provides a "wish list" of experimental and theoretical developments required in order to achieve an ideal permille precision on $\alpha_s(m^2_Z)$ within the next 10 years., Comment: 130 pages, 82 figures. White paper submitted to the Energy Frontier "Proceedings of the US Community Study on the Future of Particle Physics" (Snowmass 2021)

We report on the status of the analysis of the static energy in $2+1+1$-flavor QCD. The static energy is obtained by measuring Wilson line correlators in Coulomb gauge using the HISQ action, yielding the scales $r_{0}/a$, $r_{1}/a$, $r_{2}/a$, their ratios, and the string tension $\sigma r_{i}^{2}$. We put emphasis on the possible effects due to the dynamical charm-quark by comparing the lattice results to continuum results of the static energy with and without a massive flavor at two-loop accuracy. We employ gauge-field ensembles from the HotQCD and MILC Collaborations., Comment: 9 pages, 2 figures, The 38th International Symposium on Lattice Field Theory, LATTICE2021, 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present the first unquenched lattice-QCD calculation of the form factors for the decay $B\rightarrow D^\ast\ell\nu$ at nonzero recoil. Our analysis includes 15 MILC ensembles with $N_f=2+1$ flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from $a\approx 0.15$ fm down to $0.045$ fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence $b$ and $c$ quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element $|V_{cb}|$. We obtain $\left|V_{cb}\right| = (38.40 \pm 0.68_{\textrm{th}} \pm 0.34_{\textrm{exp}} \pm 0.18_{\textrm{EM}})\times 10^{-3}$. The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall $\chi^2\text{/dof} = 126/84$, which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict $R(D^\ast) = 0.265 \pm 0.013$, which confirms the current tension between theory and experiment., Comment: 46 pages, 14 figures. Synthetic data, results and full correlation matrices available in the ancillary files. Version accepted for publication in EPJ C

This work continues our program of lattice-QCD baryon physics using staggered fermions for both the sea and valence quarks. We present a proof-of-concept study that demonstrates, for the first time, how to calculate baryon matrix elements using staggered quarks for the valence sector. We show how to relate the representations of the continuum staggered flavor-taste group $\text{SU}(8)_{FT}$ to those of the discrete lattice symmetry group. The resulting calculations yield the normalization factors relating staggered baryon matrix elements to their physical counterparts. We verify this methodology by calculating the isovector vector and axial-vector charges $g_V$ and $g_A$. We use a single ensemble from the MILC Collaboration with 2+1+1 flavors of sea quark, lattice spacing $a\approx 0.12$ fm, and a pion mass $M_\pi\approx305$ MeV. On this ensemble, we find results consistent with expectations from current conservation and neutron beta decay. Thus, this work demonstrates how highly-improved staggered quarks can be used for precision calculations of baryon properties, and, in particular, the isovector nucleon charges., Comment: 19 pages, 9 figures

In neutrino interactions with nucleons and nuclei, Shallow Inelastic Scattering (SIS) refers to processes, dominated by non-resonant contributions, in the kinematic region where $Q^2$ is small and the invariant mass of the hadronic system, $W$, is above the pion production threshold. The extremely rich science of this complex region, poorly understood both theoretically and experimentally, encompasses the transition from interactions described in terms of hadronic degrees of freedom to interactions with quarks and gluons described by perturbative QCD. Since a large fraction of events in NOvA and DUNE, and in atmospheric neutrino measurements such as IceCube-Upgrade, KM3NeT, Super- and Hyper-Kamiokande, are from this SIS region, there is a definite need to improve our knowledge of this physics. This LoI summarizes the current understandings of the SIS physics and a series of proposals for the path to forward., Comment: 6 pages, Letter of Interest for Snowmass 2021 (https://www.snowmass21.org/docs/files/summaries/NF/SNOWMASS21-NF6_NF1-TF11_TF0-CompF2_CompF0_Katori-094.pdf), NF094

Two of the elements of the Cabibbo-Kobayashi-Maskawa quark mixing matrix, $|V_{ub}|$ and $|V_{cb}|$, are extracted from semileptonic B decays. The results of the B factories, analysed in the light of the most recent theoretical calculations, remain puzzling, because for both $|V_{ub}|$ and $|V_{cb}|$ the exclusive and inclusive determinations are in clear tension. Further, measurements in the $\tau$ channels at Belle, Babar, and LHCb show discrepancies with the Standard Model predictions, pointing to a possible violation of lepton flavor universality. LHCb and Belle II have the potential to resolve these issues in the next few years. This article summarizes the discussions and results obtained at the MITP workshop held on April 9--13, 2018, in Mainz, Germany, with the goal to develop a medium-term strategy of analyses and calculations aimed at solving the puzzles. Lattice and continuum theorists working together with experimentalists have discussed how to reshape the semileptonic analyses in view of the much higher luminosity expected at Belle II, searching for ways to systematically validate the theoretical predictions in both exclusive and inclusive B decays, and to exploit the rich possibilities at LHCb., Comment: 77 pages

We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $\alpha$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(\alpha^5)$ with negligible numerical uncertainty. The electroweak contribution is suppressed by $(m_\mu/M_W)^2$ and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at $\mathcal{O}(\alpha^2)$ and is due to hadronic vacuum polarization, whereas at $\mathcal{O}(\alpha^3)$ the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first-principle, lattice-QCD approach. The final result reads $a_\mu^\text{SM}=116\,591\,810(43)\times 10^{-11}$ and is smaller than the Brookhaven measurement by 3.7$\sigma$. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. This and the prospects to further reduce the theoretical uncertainty in the near future-which are also discussed here-make this quantity one of the most promising places to look for evidence of new physics., Comment: 196 pages, 103 figures, version published in Phys. Rept., bib files for the citation references are available from: https://muon-gm2-theory.illinois.edu

We report updates to an ongoing lattice-QCD calculation of the form factors for the semileptonic decays $B \to \pi \ell \nu$, $B_s \to K \ell \nu$, $B \to \pi \ell^+ \ell^-$, and $B \to K \ell^+ \ell^-$. The tree-level decays $B_{(s)} \to \pi (K) \ell \nu$ enable precise determinations of the CKM matrix element $|V_{ub}|$, while the flavor-changing neutral-current interactions $B \to \pi (K) \ell^+ \ell^-$ are sensitive to contributions from new physics. This work uses MILC's (2+1+1)-flavor HISQ ensembles at approximate lattice spacings between $0.057$ and $0.15$ fm, with physical sea-quark masses on four out of the seven ensembles. The valence sector is comprised of a clover $b$ quark (in the Fermilab interpretation) and HISQ light and $s$ quarks. We present preliminary results for the form factors $f_0$, $f_+$, and $f_T$, including studies of systematic errors., Comment: 7 pages; 3 figures; presented at the 37th International Symposium on Lattice Field Theory, 16-22 June 2019, Wuhan, China

We present nearly final results from our analysis of the form factors for $B\to D^\ast\ell\nu$ decay at nonzero recoil. Our analysis includes 15 MILC asqtad ensembles with $N_f=2+1$ flavors of sea quarks and lattice spacings ranging from $a\approx0.15$ fm down to $0.045$ fm. The valence light quarks employ the asqtad action, whereas the $b$ and $c$ quarks are treated using the Fermilab action. We discuss the impact that our results will have on $\ |V_{cb}\ |$ and $R(D^\ast)$., Comment: 7 pages, 3 figures. Proceedings of the 37th International Symposium on Lattice Field Theory - Lattice 2019, Wuhan (China)

We present an update on the ongoing calculations by the Fermilab Lattice, HPQCD, and MILC Collaboration of the leading-order (in electromagnetism) hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon. Our project employs ensembles with four flavors of highly improved staggered fermions, physical light-quark masses, and four lattice spacings ranging from $a \approx 0.06$ to 0.15 fm for most of the results thus far., Comment: LATTICE 2019, 7 pages, 7 figures

We present the first computation in a program of lattice-QCD baryon physics using staggered fermions for sea and valence quarks. For this initial study, we present a calculation of the nucleon mass, obtaining $964\pm16$ MeV with all sources of statistical and systematic errors controlled and accounted for. This result is the most precise determination to date of the nucleon mass from first principles. We use the highly-improved staggered quark action, which is computationally efficient. Three gluon ensembles are employed, which have approximate lattice spacings $a=0.09$ fm, $0.12$ fm, and $0.15$ fm, each with equal-mass $u$/$d$, $s$, and $c$ quarks in the sea. Further, all ensembles have the light valence and sea $u$/$d$ quarks tuned to reproduce the physical pion mass, avoiding complications from chiral extrapolations or nonunitarity. Our work opens a new avenue for precise calculations of baryon properties, which are both feasible and relevant to experiments in particle and nuclear physics., Comment: 33 pages, 19 figures; published in Physical Review D

The current status of the lattice-QCD calculations of the form factors of the $B\to D^\ast\ell\nu$ semileptonic decay is reviewed. Particular emphasis is given to the most mature calculation at non-zero recoil coming from the Fermilab Lattice and MILC collaborations. Blinded, preliminary results for the form factors are shown, including a preliminary, but detailed error budget. The lattice results seem to favor a large slope at small recoil, in contrast to the latest untagged results coming from the Belle collaboration. A comprehensive comparison between the latest BGL $z$ expansions of Belle, Babar, the lattice and a joint BGL fit including lattice and Belle data is presented, and a roadmap to improve the current calculation is discussed. The current implications for $V_{cb}$ and $R(D^\ast)$ are discussed., Comment: 8 pages, 12 figures. Proceedings of the 2019 Conference on Flavor Physics and CP Violation (FPCP2019). v2 fixes the contour plots of the coefficients of the BGL expansion for BaBar results (b1/b0 and c1/c0)

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in neutrino-oscillation physics, which inevitably entails nucleon and nuclear structure. In addition to discussing pertinent lattice-QCD calculations of nucleon and nuclear matrix elements, the interplay with models of nuclei is discussed. This program of lattice- QCD calculations is relevant to current and upcoming neutrino experiments, becoming increasingly important on the timescale of LBNF/DUNE and HyperK., Comment: 33 pp. One of seven whitepapers from the USQCD Collaboration

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in quark and lepton flavor physics. New data generated at Belle II, LHCb, BES III, NA62, KOTO, and Fermilab E989, combined with precise calculations of the relevant hadronic physics, may reveal what lies beyond the Standard Model. We outline a path toward improvements of the precision of existing lattice-QCD calculations and discuss groundbreaking new methods that allow lattice QCD to access new observables., Comment: USQCD whitepaper

We calculate the contribution to the muon anomalous magnetic moment hadronic vacuum polarization from {the} connected diagrams of up and down quarks, omitting electromagnetism. We employ QCD gauge-field configurations with dynamical $u$, $d$, $s$, and $c$ quarks and the physical pion mass, and analyze five ensembles with lattice spacings ranging from $a \approx 0.06$ to~0.15~fm. The up- and down-quark masses in our simulations have equal masses $m_l$. We obtain, in this world where all pions have the mass of the $\pi^0$, $10^{10} a_\mu^{ll}({\rm conn.}) = 637.8\,(8.8)$, in agreement with independent lattice-QCD calculations. We then combine this value with published lattice-QCD results for the connected contributions from strange, charm, and bottom quarks, and an estimate of the uncertainty due to the fact that our calculation does not include strong-isospin breaking, electromagnetism, or contributions from quark-disconnected diagrams. Our final result for the total $\mathcal{O}(\alpha^2)$ hadronic vacuum polarization to the muon's anomalous magnetic moment is~$10^{10}a_\mu^{\rm HVP,LO} = 699(15)_{u,d}(1)_{s,c,b}$, where the errors are from the light-quark and heavy-quark contributions, respectively. Our result agrees with both {\it ab-initio} lattice-QCD calculations and phenomenological determinations from experimental $e^+e^-$-scattering data. It is $1.3\sigma$ below the "no new physics" value of the hadronic-vacuum-polarization contribution inferred from combining the BNL E821 measurement of $a_\mu$ with theoretical calculations of the other contributions., Comment: 19 pages, 12 figures, 6 tables; updated to correct a small mistake in the finite volume correction resulting in small changes to the results, matches published version

We discuss preliminary results for the vector form factors $f_+^{\{\pi,K\}}$ at zero-momentum transfer for the decays $D\to\pi\ell\nu$ and $D\to K \ell\nu$ using MILC's $N_f = 2+1+1$ HISQ ensembles at four lattice spacings, $a \approx 0.042, 0.06, 0.09$, and 0.12 fm, and various HISQ quark masses down to the (degenerate) physical light quark mass. We use the kinematic constraint $f_+(q^2)= f_0(q^2)$ at $q^2 = 0$ to determine the vector form factor from our study of the scalar current, which yields $f_0(0)$. Results are extrapolated to the continuum physical point in the framework of hard pion/kaon SU(3) heavy-meson-staggered $\chi$PT and Symanzik effective theory. Our calculation improves upon the precision achieved in existing lattice-QCD calculations of the vector form factors at $q^2=0$. We show the values of the CKM matrix elements $|V_{cs}|$ and $|V_{cd}|$ that we would obtain using our preliminary results for the form factors together with recent experimental results, and discuss the implications of these values for the second row CKM unitarity., Comment: 10 pages, 3 figures, proceeding of The 36th Annual International Symposium on Lattice Field Theory

We use lattice QCD to calculate the form factors $f_+(q^2)$ and $f_0(q^2)$ for the semileptonic decay $B_s\to K\ell\nu$. Our calculation uses six MILC asqtad 2+1 flavor gauge-field ensembles with three lattice spacings. At the smallest and largest lattice spacing the light-quark sea mass is set to 1/10 the strange-quark mass. At the intermediate lattice spacing, we use four values for the light-quark sea mass ranging from 1/5 to 1/20 of the strange-quark mass. We use the asqtad improved staggered action for the light valence quarks, and the clover action with the Fermilab interpolation for the heavy valence bottom quark. We use SU(2) hard-kaon heavy-meson rooted staggered chiral perturbation theory to take the chiral-continuum limit. A functional $z$ expansion is used to extend the form factors to the full kinematic range. We present predictions for the differential decay rate for both $B_s\to K\mu\nu$ and $B_s\to K\tau\nu$. We also present results for the forward-backward asymmetry, the lepton polarization asymmetry, ratios of the scalar and vector form factors for the decays $B_s\to K\ell\nu$ and $B_s\to D_s \ell\nu$. Our results, together with future experimental measurements, can be used to determine the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{ub}|$., Comment: 57 pages, 22 figures, 13 tables

We present preliminary blinded results from our analysis of the form factors for $B\rightarrow D^\ast\ell\nu$ decay at non-zero recoil. Our analysis includes 15 MILC asqtad ensembles with $N_f=2+1$ flavors of sea quarks and lattice spacings ranging from $a\approx 0.15$ fm down to $0.045$ fm. The valence light quarks employ the asqtad action, whereas the $b$ and $c$ quarks are treated using the Fermilab action. We discuss the impact that our results will have on $\left|V_{cb}\right|$ and $R(D^\ast)$., Comment: 6 pages, 3 figures, proceedings of the 36th Annual International Symposium on Lattice Field Theory - LATTICE2018

We present high-precision results from lattice QCD for the mass splittings of the low-lying charmonium states. For the valence charm quark, the calculation uses Wilson-clover quarks in the Fermilab interpretation. The gauge-field ensembles are generated in the presence of up, down, and strange sea quarks, based on the improved staggered (asqtad) action, and gluon fields, based on the one-loop, tadpole-improved gauge action. We use five lattice spacings and two values of the light sea quark mass to extrapolate the results to the physical point. An enlarged set of interpolating operators is used for a variational analysis to improve the determination of the energies of the ground states in each channel. We present and implement a continuum extrapolation within the Fermilab interpretation, based on power-counting arguments, and thoroughly discuss all sources of systematic uncertainty. We compare our results for various mass splittings with their experimental values, namely, the 1S hyperfine splitting, the 1P-1S splitting and the P-wave spin-orbit and tensor splittings. Given the uncertainty related to the width of the resonances, we find excellent agreement., Comment: 24 pages, 8 figures

We describe a recent lattice-QCD calculation of the leptonic decay constants of heavy-light pseudoscalar mesons containing charm and bottom quarks and of the masses of the up, down, strange, charm, and bottom quarks. Results for these quantities are of the highest precision to date. Calculations use 24 isospin-symmetric ensembles of gauge-field configurations with six different lattice spacings as small as approximately 0.03 fm and several values of the light quark masses down to physical values of the average up- and down-sea-quark masses. We use the highly-improved staggered quark (HISQ) formulation for valence and sea quarks, including the bottom quark. The analysis employs heavy-quark effective theory (HQET). A novel HQET method is used in the determination of the quark masses., Comment: Talk presented CIPANP2018. 11 pages, LaTeX, 8 pdf figures

Using HISQ $N_f=2+1+1$ MILC ensembles with five different values of the lattice spacing, including four ensembles with physical quark masses, we have performed the most precise computation to date of the $K\to\pi\ell\nu$ vector form factor at zero momentum transfer, $f_+^{K^0\pi^-}(0)=0.9696(15)_\text{stat}(12)_\text{syst}$. This is the first calculation that includes the dominant finite-volume effects, as calculated in chiral perturbation theory at next-to-leading order. Our result for the form factor provides a direct determination of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{us}|=0.22333(44)_{f_+(0)}(42)_\text{exp}$, with a theory error that is, for the first time, at the same level as the experimental error. The uncertainty of the semileptonic determination is now similar to that from leptonic decays and the ratio $f_{K^+}/f_{\pi^+}$, which uses $|V_{ud}|$ as input. Our value of $|V_{us}|$ is in tension at the 2--$2.6\sigma$ level both with the determinations from leptonic decays and with the unitarity of the CKM matrix. In the test of CKM unitarity in the first row, the current limiting factor is the error in $|V_{ud}|$, although a recent determination of the nucleus-independent radiative corrections to superallowed nuclear $\beta$ decays could reduce the $|V_{ud}|^2$ uncertainty nearly to that of $|V_{us}|^2$. Alternative unitarity tests using only kaon decays, for which improvements in the theory and experimental inputs are likely in the next few years, reveal similar tensions. As part of our analysis, we calculated the correction to $f_+^{K\pi}(0)$ due to nonequilibrated topological charge at leading order in chiral perturbation theory, for both the full-QCD and the partially-quenched cases. We also obtain the combination of low-energy constants in the chiral effective Lagrangian $[C_{12}^r+C_{34}^r-(L_5^r)^2](M_\rho)=(2.92\pm0.31)\cdot10^{-6}$., Comment: 42 pages and 12 figures. Expanded discussion of fit methodology. Finite volume error increased, conclusions unchanged. Version accepted by Phys. Rev. D

We present the physics program of the Belle II experiment, located on the intensity frontier SuperKEKB $e^+e^-$ collider. Belle II collected its first collisions in 2018, and is expected to operate for the next decade. It is anticipated to collect 50/ab of collision data over its lifetime. This book is the outcome of a joint effort of Belle II collaborators and theorists through the Belle II theory interface platform (B2TiP), an effort that commenced in 2014. The aim of B2TiP was to elucidate the potential impacts of the Belle II program, which includes a wide scope of physics topics: B physics, charm, tau, quarkonium, electroweak precision measurements and dark sector searches. It is composed of nine working groups (WGs), which are coordinated by teams of theorist and experimentalists conveners: Semileptonic and leptonic B decays, Radiative and Electroweak penguins, phi_1 and phi_2 (time-dependent CP violation) measurements, phi_3 measurements, Charmless hadronic B decay, Charm, Quarkonium(like), tau and low-multiplicity processes, new physics and global fit analyses. This book highlights "golden- and silver-channels", i.e. those that would have the highest potential impact in the field. Theorists scrutinised the role of those measurements and estimated the respective theoretical uncertainties, achievable now as well as prospects for the future. Experimentalists investigated the expected improvements with the large dataset expected from Belle II, taking into account improved performance from the upgraded detector., Comment: 689 pages

We calculate the up-, down-, strange-, charm-, and bottom-quark masses using the MILC highly improved staggered-quark 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 sea-quark masses. We use a new method based on heavy-quark effective theory (HQET) to extract quark masses from heavy-light pseudoscalar meson masses. Combining our analysis with our separate determination of ratios of light-quark 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.130(41)$~MeV, $m_d(2~\text{GeV}) = 4.675(56)$~MeV, $m_s(2~\text{GeV}) = 92.47(69)$~MeV, $m_c(3~\text{GeV}) = 983.7(5.6)$~MeV, and $m_c(m_c) = 1273(10)$~MeV, with four active flavors; and $m_b(m_b) = 4195(14)$~MeV with five active flavors. We also obtain ratios of quark masses $m_c/m_s = 11.783(25)$, $m_b/m_s = 53.94(12)$, and $m_b/m_c = 4.578(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}=555(31)$~MeV in the minimal renormalon-subtracted (MRS) scheme, $\mu_\pi^2 = 0.05(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 [hep-ph]] is the key new aspect of our method., Comment: The published version; 32 pages and 7 figures

We calculate the leptonic decay constants of heavy-light pseudoscalar mesons with charm and bottom quarks in lattice quantum chromodynamics on four-flavor QCD gauge-field configurations with dynamical $u$, $d$, $s$, and $c$ quarks. We analyze over twenty isospin-symmetric ensembles with six lattice spacings down to $a\approx 0.03$~fm and several values of the light-quark mass down to the physical value $\frac{1}{2}(m_u+m_d)$. We employ the highly-improved staggered-quark (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 to-date of the $D$- and $B$-meson decay constants and their ratios, $f_{D^+} = 212.7(0.6)$~MeV, $f_{D_s} = 249.9(0.4)$~MeV, $f_{D_s}/f_{D^+} = 1.1749(16)$, $f_{B^+} = 189.4 (1.4)$~MeV, $f_{B_s} = 230.7(1.3)$~MeV, $f_{B_s}/f_{B^+} = 1.2180(47)$, 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 lattice-QCD calculations, and bring the QCD errors in the Standard-Model predictions for the rare leptonic decays $\overline{\mathcal{B}}(B_s \to \mu^+\mu^-) = 3.64(11) \times 10^{-9}$, $\overline{\mathcal{B}}(B^0 \to \mu^+\mu^-) = 1.00(3) \times 10^{-10}$, and $\overline{\mathcal{B}}(B^0 \to \mu^+\mu^-)/\overline{\mathcal{B}}(B_s \to \mu^+\mu^-) = 0.0273(9)$ to well below other sources of uncertainty. As a byproduct of our analysis, we also update our previously published results for the light-quark-mass ratios and the scale-setting quantities $f_{p4s}$, $M_{p4s}$, and $R_{p4s}$. We obtain the most precise lattice-QCD determination to date of the ratio $f_{K^+}/f_{\pi^+} = 1.1950(^{+16}_{-23})$~MeV., Comment: Errors related to the standard model prediction for the rare leptonic decays are fixed in the abstract and Eqs. (7.44), (7.45), and (8.3)

The study of heavy-light meson masses should provide a way to determine renormalized quark masses and other properties of heavy-light mesons. In the context of lattice QCD, for example, it is possible to calculate hadronic quantities for arbitrary values of the quark masses. In this paper, we address two aspects relating heavy-light meson masses to the quark masses. First, we introduce a definition of the renormalized quark mass that is free of both scale dependence and renormalon ambiguities, and discuss its relation to more familiar definitions of the quark mass. We then show how this definition enters a merger of the descriptions of heavy-light masses in heavy-quark effective theory and in chiral perturbation theory ($\chi$PT). For practical implementations of this merger, we extend the one-loop $\chi$PT corrections to lattice gauge theory with heavy-light mesons composed of staggered fermions for both quarks. Putting everything together, we obtain a practical formula to describe all-staggered heavy-light meson masses in terms of quark masses as well as some lattice artifacts related to staggered fermions. In a companion paper, we use this function to analyze lattice-QCD data and extract quark masses and some matrix elements defined in heavy-quark effective theory., Comment: 22 pp

Using the MILC 2+1 flavor asqtad quark action ensembles, we are calculating the form factors $f_0$ and $f_+$ for the semileptonic $B_s \rightarrow K \ell\nu$ decay. A total of six ensembles with lattice spacing from $\approx0.12$ to 0.06 fm are being used. At the coarsest and finest lattice spacings, the light quark mass $m'_l$ is one-tenth the strange quark mass $m'_s$. At the intermediate lattice spacing, the ratio $m'_l/m'_s$ ranges from 0.05 to 0.2. The valence $b$ quark is treated using the Sheikholeslami-Wohlert Wilson-clover action with the Fermilab interpretation. The other valence quarks use the asqtad action. When combined with (future) measurements from the LHCb and Belle II experiments, these calculations will provide an alternate determination of the CKM matrix element $|V_{ub}|$., Comment: 8 pages, 6 figures, to appear in the Proceedings of Lattice 2017, June 18-24, Granada, Spain

All lattice-QCD calculations of the hadronic-vacuum-polarization contribution to the muon's anomalous magnetic moment to-date have been performed with degenerate up- and down-quark masses. Here we calculate directly the strong-isospin-breaking 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 lattice-QCD results obtained with degenerate light-quark masses of $\delta a_\mu^{{\rm HVP,} m_u \neq m_d}$= +1.5(7)%, in agreement with estimates from phenomenology and a recent lattice-QCD calculation with unphysically heavy pions., Comment: v2: 6 pages, 2 tables, 2 figures. Additional references and expanded discussion of systematic errors. Version accepted to Physical Review Letters

We present preliminary results from our analysis of the form factors for the $\overline{B}\rightarrow D^\ast\ell\overline{\nu}$ decay at non-zero recoil. Our analysis includes 15 MILC asqtad ensembles with $N_f=2+1$ flavors of sea quarks and lattice spacings ranging from $a\approx 0.15$ fm down to $0.045$ fm. The valence light quarks employ the asqtad action, whereas the heavy quarks are treated using the Fermilab action. We conclude with a discussion of future plans and phenomenological implications. When combined with experimental measurements of the decay rate, our calculation will enable a determination of the CKM matrix element $|V_{cb}|$., Comment: 8 pages, 5 figures, Proceedings of the talk presented at the 35th International Symposium on Lattice Field Theory (Lattice2017), Granada, Spain

We report the status of an ongoing lattice-QCD 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 sea-quark masses. We present preliminary results for correlation-function fits., Comment: 8 pages, 3 figures; presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

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., Comment: Published version, 42 pages, 18 figures

The precise measurement of neutrino properties is among the highest priorities in fundamental particle physics, involving many experiments worldwide. Since the experiments rely on the interactions of neutrinos with bound nucleons inside atomic nuclei, the planned advances in the scope and precision of these experiments requires a commensurate effort in the understanding and modeling of the hadronic and nuclear physics of these interactions, which is incorporated as a nuclear model in neutrino event generators. This model is essential to every phase of experimental analyses and its theoretical uncertainties play an important role in interpreting every result. In this White Paper we discuss in detail the impact of neutrino-nucleus interactions, especially the nuclear effects, on the measurement of neutrino properties using the determination of oscillation parameters as a central example. After an Executive Summary and a concise Overview of the issues, we explain how the neutrino event generators work, what can be learned from electron-nucleus interactions and how each underlying physics process - from quasi-elastic to deep inelastic scattering - is understood today. We then emphasize how our understanding must improve to meet the demands of future experiments. With every topic we find that the challenges can be met only with the active support and collaboration among specialists in strong interactions and electroweak physics that include theorists and experimentalists from both the nuclear and high energy physics communities.

We present results for neutral D-meson mixing in 2+1-flavor lattice QCD. We compute the matrix elements for all five operators that contribute to D mixing at short distances, including those that only arise beyond the Standard Model. Our results have an uncertainty similar to those of the ETM collaboration (with 2 and with 2+1+1 flavors). This work shares many features with a recent publication on B mixing and with ongoing work on heavy-light decay constants from the Fermilab Lattice and MILC Collaborations., Comment: 6+1 pp., presented at Lattice 2016