Search

We compute the standard Euclidean window of the hadronic vacuum polarization using multiple independent blinded analyses. We improve the continuum and infinite-volume extrapolations of the dominant quark-connected light-quark isospin-symmetric contribution and address additional sub-leading systematic effects from sea-charm quarks and residual chiral-symmetry breaking from first principles. We find $a_\mu^{\rm W} = 235.56(65)(50) \times 10^{-10}$, which is in $3.8\sigma$ tension with the recently published dispersive result of Colangelo et al., $a_\mu^{\rm W} = 229.4(1.4) \times 10^{-10}$, and in agreement with other recent lattice determinations. We also provide a result for the standard short-distance window. The results reported here are unchanged compared to our presentation at the Edinburgh workshop of the g-2 Theory Initiative in 2022., Comment: 24 pages, 15 figures

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., Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

Phase shifts for $s$-wave $����$ scattering in both the $I=0$ and $I=2$ channels are determined from a lattice QCD calculation performed on 741 gauge configurations obeying G-parity boundary conditions with a physical pion mass and lattice size of $32^3\times 64$. These results support our recent study of direct CP violation in $K\to����$ decay \cite{Abbott:2020hxn}, improving our earlier 2015 calculation \cite{Bai:2015nea}. The phase shifts are determined for both stationary and moving $����$ systems, at three ($I=0$) and four ($I=2$) different total momenta. We implement several $����$ interpolating operators including a scalar bilinear "$��$" operator and paired single-pion bilinear operators with the constituent pions carrying various relative momenta. Several techniques, including correlated fitting and a bootstrap determination of p-values have been used to refine the results and a comparison with the generalized eigenvalue problem (GEVP) method is given. A detailed systematic error analysis is performed which allows phase shift results to be presented at a fixed energy., v3: Add a subsection "Higher partial wave correction", and correct the unit of scattering length. 88 pages and 14 figures v2: 1). Add reference 29 as an example of pipi scattering calculation above 4mpi threshold. 2). Modify the wording on page 3 for the footage. 3). Correct the \sigma operator on page 17

RBC/UKQCD is preparing a calculation of leptonic decay rates including isospin breaking corrections using a perturbative approach to include NLO contributions from QED effects. We present preliminary numerical results for a contribution to the leptonic pion decay rate and report on exploratory studies of computational techniques based on all-to-all propagators., Comment: Proceedings

We compute the renormalisation factors (Z-matrices) of the $��F=2$ four-quark operators needed for Beyond the Standard Model (BSM) kaon mixing. We work with nf=2+1 flavours of Domain-Wall fermions whose chiral-flavour properties are essential to maintain a continuum-like mixing pattern. We introduce new RI-SMOM renormalisation schemes, which we argue are better behaved compared to the commonly-used corresponding RI-MOM one. We find that, once converted to MS, the Z-factors computed through these RI-SMOM schemes are in good agreement but differ significantly from the ones computed through the RI-MOM scheme. The RI-SMOM Z-factors presented here have been used to compute the BSM neutral kaon mixing matrix elements in the companion paper [1]. We argue that the renormalisation procedure is responsible for the discrepancies observed by different collaborations, we will investigate and elucidate the origin of these differences throughout this work.

In this document we address an error discovered in the ensemble generation for our calculation of the $I=0$ $K\to����$ amplitude (Phys. Rev. Lett. 115, 212001 (2015), arXiv:1505.07863) whereby the same random numbers were used for the two independent quark flavors, resulting in small but measurable correlations between gauge observables separated by 12 units in the y-direction. We conclude that the effects of this error are negligible compared to the overall errors on our calculation., 2 pages, 1 figure

We report on a calculation of the B*Bpi coupling in lattice QCD. The strong matrix element for a B* to Bpi transition is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HMChPT) for B mesons. We carry out our calculation directly at the b-quark mass using a non-perturbatively tuned clover action that controls discretization effects of order pa and (ma)^n for all n. Our analysis is performed on RBC/UKQCD gauge configurations using domain-wall fermions and the Iwasaki gauge action at two lattice spacings of ainverse = 1.729(25) GeV, ainverse = 2.281(28) GeV, and unitary pion masses down to 290 MeV. We achieve good statistical precision and control all systematic uncertainties, giving a final result for the HMChPT coupling g_b = 0.56(3)stat(7)sys in the continuum and at the physical light-quark masses. This is the first calculation performed directly at the physical b-quark mass and lies in the region one would expect from carrying out an interpolation between previous results at the charm mass and at the static point., Comment: Postprint: version accepted for publication

We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first value is in approximate agreement with the experimental result: Re$(A_0) = 3.3201(18) \times 10^{-7}$ GeV while the second can be used to compute the direct CP violating ratio Re$(\varepsilon'/\varepsilon)=1.38(5.15)(4.59)\times 10^{-4}$, which is $2.1��$ below the experimental value $16.6(2.3)\times 10^{-4}$. The real part of $A_0$ is CP conserving and serves as a test of our method while the result for Re$(\varepsilon'/\varepsilon)$ provides a new test of the standard-model theory of CP violation, one which can be made more accurate with increasing computer capability., 9 pages, 3 figures. Updated to match published version

We calculate the $B \to��\ell��$ and $B_s \to K \ell��$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses as light as $M_��\approx 290$ MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) "hard-pion" chiral perturbation theory. We provide complete error budgets for the form factors $f_+(q^2)$ and $f_0(q^2)$ at three momenta that span the $q^2$ range accessible in our numerical simulations. We extrapolate these results to $q^2 = 0$ using a model-independent $z$-parametrization and present our final form factors as the $z$-coefficients and the matrix of correlations between them. Our results agree with other lattice determinations using staggered light quarks and provide important independent cross-checks. Both $B \to��\ell��$ and $B_s \to K \ell��$ decays enable a determination of the CKM matrix element $|V_{ub}|$. To illustrate this, we perform a combined $z$-fit of our numerical $B\to��\ell��$ form-factor data with the experimental branching-fraction measurements leaving the relative normalization as a free parameter; we obtain $|V_{ub}| = 3.61(32) \times 10^{-3}$, where the error includes statistical and systematic uncertainties. This approach can be applied to $B_s\to K \ell��$ decay to determine $|V_{ub}|$ once the process has been measured experimentally. Finally, in anticipation of future measurements, we make predictions for $B \to ��\ell��$ and $B_s\to K \ell��$ Standard-Model differential branching fractions and forward-backward asymmetries., V3: References and discussions added. Large tables move to appendix. 39 pages, 22 tables, 23 figures. Version published in Physical Review D

We present new results for the amplitude $A_2$ for a kaon to decay into two pions with isospin $I=2$: Re$A_2 = 1.50(4)_\mathrm{stat}(14)_\mathrm{syst}\times 10^{-8}$ GeV; Im$A_2 = -6.99(20)_\mathrm{stat}(84)_\mathrm{syst}\times 10^{-13}$ GeV. These results were obtained from two ensembles generated at physical quark masses (in the isospin limit) with inverse lattice spacings $a^{-1}=1.728(4)$ GeV and $2.358(7)$ GeV. We are therefore able to perform a continuum extrapolation and hence largely to remove the dominant systematic uncertainty from our earlier results, that due to lattice artefacts. The only previous lattice computation of $K\to����$ decays at physical kinematics was performed using an ensemble at a single, rather coarse, value of the lattice spacing ($a^{-1}\simeq 1.37(1)$ GeV). We confirm the observation that there is a significant cancellation between the two dominant contributions to Re$A_2$ which we suggest is an important ingredient in understanding the $��I=1/2$ rule, Re$A_0$/Re$A_2\simeq 22.5$, where the subscript denotes the total isospin of the two-pion final state. Our result for $A_2$ implies that the electroweak penguin contribution to $��^\prime/��$ is Re($��^\prime/��)_\textrm{EWP}=-(6.6\pm 1.0)\times 10^{-4}$., 46 pages, 12 figures

We calculate the B-meson decay constants f_B, f_Bs, and their ratio in unquenched lattice QCD using domain-wall light quarks and relativistic b-quarks. We use gauge-field ensembles generated by the RBC and UKQCD collaborations using the domain-wall fermion action and Iwasaki gauge action with three flavors of light dynamical quarks. We analyze data at two lattice spacings of a ~ 0.11, 0.086 fm with unitary pion masses as light as M_pi ~ 290 MeV; this enables us to control the extrapolation to the physical light-quark masses and continuum. For the b-quarks we use the anisotropic clover action with the relativistic heavy-quark interpretation, such that discretization errors from the heavy-quark action are of the same size as from the light-quark sector. We renormalize the lattice heavy-light axial-vector current using a mostly nonperturbative method in which we compute the bulk of the matching factor nonperturbatively, with a small correction, that is close to unity, in lattice perturbation theory. We also improve the lattice heavy-light current through O(alpha_s a). We extrapolate our results to the physical light-quark masses and continuum using SU(2) heavy-meson chiral perturbation theory, and provide a complete systematic error budget. We obtain f_B0 = 199.5(12.6) MeV, f_B+ = 195.6(14.9) MeV, f_Bs = 235.4(12.2) MeV, f_Bs/f_B0 = 1.197(50), and f_Bs/f_B+ = 1.223(71), where the errors are statistical and total systematic added in quadrature. These results are in good agreement with other published results and provide an important independent cross check of other three-flavor determinations of $B$-meson decay constants using staggered light quarks., 25 pages, 13 tables, 12 figures. Version 2 revised for PRD. Sign in Eq. (27) corrected. Chiral-continuum extrapolations of fB and fBs/fB changed. Results for fB and fBs/fB, plots, and text updated accordingly. Updated references

We report on a calculation of the $B^*B��$ coupling in lattice QCD. The strong matrix element $\langle B ��| B^*\rangle$ is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HM$��$PT) for $B$-mesons. We carry out our calculation directly at the $b$-quark mass using a non-perturbatively tuned clover action that controls discretisation effects of order $|\vec{p}a|$ and $(ma)^n$ for all $n$. Our analysis is performed on RBC/UKQCD gauge configurations using domain wall fermions and the Iwasaki gauge action at two lattice spacings of $a^{-1}=1.73(3)$ GeV, $a^{-1}=2.28(3)$ GeV, and unitary pion masses down to 290 MeV. We achieve good statistical precision and control all systematic uncertainties, giving a final result for the HM$��$PT coupling $g_b = 0.569(48)_{stat}(59)_{sys}$ in the continuum and at the physical light-quark masses. This is the first calculation performed directly at the physical $b$-quark mass and lies in the region one would expect from carrying out an interpolation between previous results at the charm mass and at the static point., 7 pages, 2 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We describe the computation of the amplitude A_2 for a kaon to decay into two pions with isospin I=2. The results presented in the letter Phys.Rev.Lett. 108 (2012) 141601 from an analysis of 63 gluon configurations are updated to 146 configurations giving Re$A_2=1.381(46)_{\textrm{stat}}(258)_{\textrm{syst}} 10^{-8}$ GeV and Im$A_2=-6.54(46)_{\textrm{stat}}(120)_{\textrm{syst}}10^{-13}$ GeV. Re$A_2$ is in good agreement with the experimental result, whereas the value of Im$A_2$ was hitherto unknown. We are also working towards a direct computation of the $K\to(����)_{I=0}$ amplitude $A_0$ but, within the standard model, our result for Im$A_2$ can be combined with the experimental results for Re$A_0$, Re$A_2$ and $��^\prime/��$ to give Im$A_0/$Re$A_0= -1.61(28)\times 10^{-4}$ . Our result for Im\,$A_2$ implies that the electroweak penguin (EWP) contribution to $��^\prime/��$ is Re$(��^\prime/��)_{\mathrm{EWP}} = -(6.25 \pm 0.44_{\textrm{stat}} \pm 1.19_{\textrm{syst}}) \times 10^{-4}$., 59 pages, 11 figures

There has been much speculation as to the origin of the ��I = 1/2 rule (Re A_0/Re A_2 \simeq 22.5). We find that the two dominant contributions to the ��I=3/2, K \to ���� correlation functions have opposite signs leading to a significant cancellation. This partial cancellation occurs in our computation of Re A_2 with physical quark masses and kinematics (where we reproduce the experimental value of A_2) and also for heavier pions at threshold. For Re A_0, although we do not have results at physical kinematics, we do have results for pions at zero-momentum with m_�� \simeq 420 MeV (Re A_0/Re A_2=9.1(2.1)) and m_�� \simeq 330 MeV (Re A_0/Re A_2=12.0(1.7)). The contributions which partially cancel in Re A_2 are also the largest ones in Re A_0, but now they have the same sign and so enhance this amplitude. The emerging explanation of the ��I=1/2 rule is a combination of the perturbative running to scales of O(2 GeV), a relative suppression of Re A_2 through the cancellation of the two dominant contributions and the corresponding enhancement of Re A_0. QCD and EWP penguin operators make only very small contributions at such scales., 5 pages, 4 figures. v2 minor revisions to coincide w/ published version

We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two ��-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}$\,A_2=-(6.83 \pm 0.51_{\textrm{stat}} \pm 1.30_{\textrm{syst}})\,10^{-13}\,{\rm GeV}$. Moreover combining our result for Im\,$A_2$ with experimental values of Re\,$A_2$, Re\,$A_0$ and $��^\prime/��$, we obtain the following value for the unknown ratio Im\,$A_0$/Re\,$A_0$ within the Standard Model: $\mathrm{Im}\,A_0/\mathrm{Re}\,A_0=-1.63(19)_{\mathrm{stat}}(20)_{\mathrm{syst}}\times10^{-4}$. One consequence of these results is that the contribution from Im\,$A_2$ to the direct CP violation parameter $��^{\prime}$ (the so-called Electroweak Penguin, EWP, contribution) is Re$(��^\prime/��)_{\mathrm{EWP}} = -(6.52 \pm 0.49_{\textrm{stat}} \pm 1.24_{\textrm{syst}}) \times 10^{-4}$. We explain why this calculation of $A_2$ represents a major milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding of CP-violation in kaon decays., 5 pages, 1 figure

We report a direct lattice calculation of the $K$ to $����$ decay matrix elements for both the $��I=1/2$ and 3/2 amplitudes $A_0$ and $A_2$ on 2+1 flavor, domain wall fermion, $16^3\times32\times16$ lattices. This is a complete calculation in which all contractions for the required ten, four-quark operators are evaluated, including the disconnected graphs in which no quark line connects the initial kaon and final two-pion states. These lattice operators are non-perturbatively renormalized using the Rome-Southampton method and the quadratic divergences are studied and removed. This is an important but notoriously difficult calculation, requiring high statistics on a large volume. In this paper we take a major step towards the computation of the physical $K\to����$ amplitudes by performing a complete calculation at unphysical kinematics with pions of mass 422\,MeV at rest in the kaon rest frame. With this simplification we are able to resolve Re$(A_0)$ from zero for the first time, with a 25% statistical error and can develop and evaluate methods for computing the complete, complex amplitude $A_0$, a calculation central to understanding the $��=1/2$ rule and testing the standard model of CP violation in the kaon system., 40 pages, 12 figures

We are very grateful to the Fermilab Directorate and the Fermilab Theoretical Physics Department for their financial and logistical support of the first workshop of the Muon g -2 Theory Initiative (held near Fermilab in June 2017) [123], which was crucial for its success, and indeed for the successful start of the Initiative. Financial support for this workshop was also provided by the Fermilab Distinguished Scholars program, the Universities Research Association through a URA Visiting Scholar award, the Riken Brookhaven Research Center, and the Japan Society for the Promotion of Science under Grant No. KAKEHNHI-17H02906. We thank Shoji Hashimoto, Toru Iijima, Takashi Kaneko, and Shohei Nishida for hosting the HVP workshop at KEK [124] and the KEK Theory Center and the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics for providing logistical and financial support. The HLbL workshop at the University of Connecticut [125] was hosted by the University of Connecticut Physics Department. We also gratefully acknowledge support for the second plenary workshop in Mainz [126] from the Deutsche Forschungsgemeinschaft via the Cluster of Excellence "Precision Physics, Fundamental Interactions and Structure of Matter'' (PRISMA), the Collaborative Research Centre "The low-energy frontier of the Standard Model'' (SFB 1044), as well as the Helmholtz Institute Mainz. And finally, we thank the Institute for Nuclear Theory at the University of Washington for hosting the third plenary workshop [127] and for its kind hospitality and stimulating research environment. This workshop was supported in part by the U.S. Department of Energy, Office of Science, under Award Nos. DE-FG02-00ER41132, DE-SC0020106, and by the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics. This review benefited from discussions with O. Cata, N. Christ, L.Y. Dai, H. Davoudiasl, S. Fayer, S. Ganguly, A. Gasparian, S. Hashimoto, T. Iijima, K. Kampf, D. Kawall, I. Larin, Z. Pagel, M. Petschlies, A. Rebhan, K. Schilcher, K. Shimomura, E. Shintani, D. Steffen, S. Tracz, C. Tu, and T. Yamazaki. The work in this paper was supported by CNRS, by Conacyt (Ciencia Basica 2015) under Grant No. 250628, by CONACyT-Mexico under Grant No. CB2014-22117, by Coordinacion de la Investigacion Cientifica (CIC-UMSNH) under Grant No. 4.10, by Danmarks Frie Forskningsfond under Grant No. 8021-00122B, by Deutsche Forschungsgemeinschaft Collaborative Research Centers CRC 1044, CRC 1044 -204404729, CRC 110, and under Grant No. HI 2048/1-1, Prisma Cluster for Excellence PRISMA+ EXC2118/1, STO/876/6-1, by European Research Council under the European Union's Horizon 2020 research and innovation programme under Grant Agreement Nos. 668679, 757646, 771971-SIMDAMA, 813942, by the European Union H2020-MSCA-COFUND2016 under Grant No. 754510, by the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Grant Agreement No. 843134, by European Union EuroPLEx Grant H2020-MSCA-ITN-2018-813942, by European Union STRONG 2020 project under Grant Agreement No. 824093, by the Excellence Initiative of Aix-Marseille University -A*MIDEX, a French "Investissements d'Avenir" program, through the Chaire d'Excellence program and the OCEVU Laboratoire d'Excellence (ANR-11-LABX-0060), by the Fermilab Distinguished Scholars program, by Fondo SEP-Cinvestav under Grant No. 142, Fundacao para a Ciencia e a Tecnologia under Grant No. SFRH/BPD/109443/2015, by Generalitat de Catalunya under Grant No. 2017SGR1069, by the Helmholtz Association (German Federal Ministry of Education and Research), by the Helmholtz-Institut Mainz, by the Istituto Nazionale di Fisica Nucleare (INFN), by the Isaac Newton Trust, by the Japan Society for the Promotion of Science under Grant Nos. KAKENHI-15H05742, 16K05317, 16K05323, 16K05338, 17H01133, 17H02906, 18H05226, 19K21872, 20K03926, 20K03960, by Junta de Andalucia under Grant No. A-FQM-467-UGR18, by KEK, by Ministerio de Ciencia e Innovacion under Grant No. CICYTFEDER-FPA2017-86989-P, by Ministerio de Industria, Economia y Competitividad under Grant Nos. FPA2016-78220-C3-3-P, FPA2017-86989-P, PGC2018-094857-B-I00, SEV-2016-0588, SEV-2016-0597, by Laboratoires d'Excellence FIRST-TF grants, a French "Investissements d'Avenir'' program, by the Ministry of Science and Higher Education of the Russian Federation under Grant Agreement No. 14.W03.31.0026, by the National Research Foundation of South Africa, by the Natural Sciences and Engineering Research Council of Canada, by the Portuguese Science Foundation (FCT) Investigator Grant IF/00898/2015, by the Romanian Ministry of Education and Research under Grant No. PN19060101, by the Russian Science Foundation under Grant No. RSF 18-12-00128, by the Secretaria d'Universitats i Recerca del Departament d'Economia i Coneixement de la Generalitat de Catalunya under Grant No. 2017 SGR 1069, by the Swedish Research Council under Grant Nos. 2016-05996, 2019-03779, by the Swiss National Science Foundation under Grant Nos. PP00P2_176884, PCEFP2_181117, by The Leverhulme Trust under Grant No. ECF-2019-223, by the UK Science and Technology Facilities Council (STFC) under Grant Nos. ST/N504130/1, ST/P000290/1, ST/P00055X/1, ST/P000630/1, ST/P000711/1, ST/P000746/1, ST/S000879/1, ST/S000925/1, by the U.S. Department of Energy, Office of Science, Office of High Energy Physics under Award Nos. DE-SC0009919, DE-SC0009998, DE-SC0010005, DE-SC0010120, DE-SC0010339, DE-SC0012391, DE-SC0012704, DE-SC0013682, DE-SC0013895, DE-SC0015655, by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics under Award Nos. DE-AC02-05CH11231, DE-FG02-00ER41132, DE-FG02-97ER41020, by the U.S.-Japan Science and Technology Cooperation Program in High Energy Physics, "Incubation Platform for Intensity Frontier", by the U.S. National Science Foundation under Grant Nos. NSF-PHY-1316222, PHY14-14614, PHY17-19626, and PHY19-23131, and by the U.S. National Institute of Standards and Technology (NIST) Precision Measurement Grant Program under Award No. 60NANB16D271. This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics., 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 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 with negligible numerical uncertainty. The electroweak contribution is suppressed by 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 and is due to hadronic vacuum polarization, whereas at 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 and is smaller than the Brookhaven measurement by 3.7. 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., Fermilab Directorate, Fermilab Theoretical Physics Department, Fermilab Distinguished Scholars program, Universities Research Association through a URA Visiting Scholar award, Riken Brookhaven Research Center, Ministry of Education, Culture, Sports, Science and Technology, Japan (MEXT) Japan Society for the Promotion of Science KAKEHNHI-17H02906 KAKENHI-15H05742 16K05317 16K05323 16K05338 17H01133 17H02906 18H05226 19K21872 20K03926 20K03960, High Energy Accelerator Research Organization (KEK), U.S.-Japan Science and Technology Cooperation Program in High Energy Physics, Deutsche Forschungsgemeinschaft via the Cluster of Excellence "Precision Physics, Fundamental Interactions and Structure of Matter'' (PRISMA), Collaborative Research Centre "The low-energy frontier of the Standard Model'' SFB 1044, Helmholtz Institute Mainz, United States Department of Energy (DOE) DE-FG02-00ER41132 DE-SC0020106, Centre National de la Recherche Scientifique (CNRS), Conacyt (Ciencia Basica 2015) 250628, Consejo Nacional de Ciencia y Tecnologia (CONACyT) CB2014-22117, Coordinacion de la Investigacion Cientifica (CIC-UMSNH) 4.10, Danmarks Frie Forskningsfond 8021-00122B, German Research Foundation (DFG) CRC 1044 CRC 1044 -204404729 CRC 110 HI 2048/1-1, Prisma Cluster for Excellence PRISMA+ EXC2118/1 STO/876/6-1, European Research Council (ERC) 668679 757646 771971-SIMDAMA 813942, European Union (EU) 754510, European Union (EU) 843134, European Union EuroPLEx Grant H2020-MSCA-ITN-2018-813942, European Union STRONG 2020 project 824093, French National Research Agency (ANR), OCEVU Laboratoire d'Excellence ANR-11-LABX-0060, Fondo SEP-Cinvestav 142, Portuguese Foundation for Science and Technology SFRH/BPD/109443/2015, Generalitat de Catalunya 2017SGR1069, Helmholtz Association (German Federal Ministry of Education and Research), Istituto Nazionale di Fisica Nucleare (INFN), Isaac Newton Trust, Junta de Andalucia A-FQM-467-UGR18, Instituto de Salud Carlos III Spanish Government CICYTFEDER-FPA2017-86989-P, Ministerio de Industria, Economia y Competitividad FPA2016-78220-C3-3-P FPA2017-86989-P PGC2018-094857-B-I00 SEV-2016-0588 SEV-2016-0597, Laboratoires d'Excellence FIRST-TF grants, Ministry of Science and Higher Education of the Russian Federation 14.W03.31.0026, National Research Foundation - South Africa, Natural Sciences and Engineering Research Council of Canada (NSERC) CGIAR, Portuguese Science Foundation (FCT) Investigator Grant IF/00898/2015, Romanian Ministry of Education and Research PN19060101, Russian Science Foundation (RSF) RSF 18-12-00128, Generalitat de Catalunya 2017 SGR 1069, Swedish Research Council 2016-05996 2019-03779, Swiss National Science Foundation (SNSF) PP00P2_176884 PCEFP2_181117, Leverhulme Trust ECF-2019-223, UK Research & Innovation (UKRI) Science & Technology Facilities Council (STFC) ST/N504130/1 ST/P000290/1 ST/P00055X/1 ST/P000630/1 ST/P000711/1 ST/P000746/1 ST/S000879/1 ST/S000925/1, United States Department of Energy (DOE) DE-SC0009919 DE-SC0009998 DE-SC0010005 DE-SC0010120 DE-SC0010339 DE-SC0012391 DE-SC0012704 DE-SC0013682 DE-SC0013895 DE-SC0015655, United States Department of Energy (DOE) DE-FG02-00ER41132 DE-AC02-05CH11231 DE-FG02-97ER41020 DE-AC02-07CH11359, U.S.-Japan Science and Technology Cooperation Program in High Energy Physics, "Incubation Platform for Intensity Frontier", National Science Foundation (NSF) NSF-PHY-1316222 PHY14-14614 PHY17-19626 PHY19-23131, U.S. National Institute of Standards and Technology (NIST) Precision Measurement Grant 60NANB16D271

The Hadronic Vacuum Polarization (HVP) is a dominant contribution to the theoretical uncertainty of the muon anomalous magnetic moment. The uncertainty in a lattice QCD calculation of the connected light-quark contribution to the HVP is dominated by the long-distance region of the vector correlation function. Explicit studies of the exclusive channels of the HVP diagram make it possible to reconstruct the long-distance behavior of the correlation function. This removes most of the statistical uncertainty of the correlation function. In these proceedings, preliminary results of an exclusive study of the isospin symmetric connected-only vector-vector correlation function using a hybrid of distillation and A2A techniques are presented. The computation is performed on 2+1 flavor M\"obius Domain Wall Fermion ensembles with physical pion mass. Reconstruction of the long-distance correlation function will enable lattice-only calculations of the HVP to achieve precision similar to estimates of the HVP from the R-ratio method on the timescale of the new experimental measurements of the muon anomalous magnetic moment., Comment: 7 pages, 4 figures; Presented at the 37th International Symposium on Lattice Field Theory, 16-22 June 2019, Wuhan, China

We propose a method to compute the Wilson coefficients of the weak effective Hamiltonian to all orders in the strong coupling constant using Lattice QCD simulations. We perform our calculations adopting an unphysically light weak boson mass of around $2~\mathrm{GeV}$. We demonstrate that systematic errors for the Wilson coefficients $C_1$ and $C_2$, related to the current-current four-quark operators, can be controlled and present a path towards precise determinations in subsequent works., 16 pages, 9 figures