• ### Calculation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment(1801.07224)

Jan. 22, 2018 hep-ph, hep-lat
We present a first-principles lattice QCD+QED calculation at physical pion mass of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. The total contribution of up, down, strange, and charm quarks including QED and strong isospin breaking effects is found to be $a_\mu^{\rm HVP~LO}=715.4(16.3)(9.2) \times 10^{-10}$, where the first error is statistical and the second is systematic. By supplementing lattice data for very short and long distances with experimental R-ratio data using the compilation of Ref. [1], we significantly improve the precision of our calculation and find $a_\mu^{\rm HVP~LO} = 692.5(1.4)(0.5)(0.7)(2.1) \times 10^{-10}$ with lattice statistical, lattice systematic, R-ratio statistical, and R-ratio systematic errors given separately. This is the currently most precise determination of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. In addition, we present the first lattice calculation of the light-quark QED correction at physical pion mass.
• ### Neutral kaon mixing beyond the Standard Model with nf=2+1 chiral fermions part II: Non Perturbative Renormalisation of the $\Delta F=2$ four-quark operators(1708.03552)

Aug. 11, 2017 hep-ph, hep-lat
We compute the renormalisation factors (Z-matrices) of the $\Delta 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.
• ### Domain wall QCD with physical quark masses(1411.7017)

May 16, 2016 hep-lat
We present results for several light hadronic quantities ($f_\pi$, $f_K$, $B_K$, $m_{ud}$, $m_s$, $t_0^{1/2}$, $w_0$) obtained from simulations of 2+1 flavor domain wall lattice QCD with large physical volumes and nearly-physical pion masses at two lattice spacings. We perform a short, O(3)%, extrapolation in pion mass to the physical values by combining our new data in a simultaneous chiral/continuum `global fit' with a number of other ensembles with heavier pion masses. We use the physical values of $m_\pi$, $m_K$ and $m_\Omega$ to determine the two quark masses and the scale - all other quantities are outputs from our simulations. We obtain results with sub-percent statistical errors and negligible chiral and finite-volume systematics for these light hadronic quantities, including: $f_\pi$ = 130.2(9) MeV; $f_K$ = 155.5(8) MeV; the average up/down quark mass and strange quark mass in the $\bar {\rm MS}$ scheme at 3 GeV, 2.997(49) and 81.64(1.17) MeV respectively; and the neutral kaon mixing parameter, $B_K$, in the RGI scheme, 0.750(15) and the $\bar{\rm MS}$ scheme at 3 GeV, 0.530(11).
• ### Lattice calculation of the leading strange quark-connected contribution to the muon $g-2$(1602.01767)

May 3, 2016 hep-ph, hep-lat
We present results for the leading hadronic contribution to the muon anomalous magnetic moment due to strange quark-connected vacuum polarisation effects. Simulations were performed using RBC--UKQCD's $N_f=2+1$ domain wall fermion ensembles with physical light sea quark masses at two lattice spacings. We consider a large number of analysis scenarios in order to obtain solid estimates for residual systematic effects. Our final result in the continuum limit is $a_\mu^{(2)\,{\rm had},\,s}=53.1(9)\left(^{+1}_{-3}\right)\times10^{-10}$.
• ### Calculation of the hadronic vacuum polarization disconnected contribution to the muon anomalous magnetic moment(1512.09054)

Dec. 30, 2015 hep-ph, hep-lat
We report the first lattice QCD calculation of the hadronic vacuum polarization disconnected contribution to the muon anomalous magnetic moment at physical pion mass. The calculation uses a refined noise-reduction technique which enabled the control of statistical uncertainties at the desired level with modest computational effort. Measurements were performed on the $48^3 \times 96$ physical-pion-mass lattice generated by the RBC and UKQCD collaborations. We find $a_\mu^{\rm HVP~(LO)~DISC} = -9.6(3.3)(2.3)\times 10^{-10}$, where the first error is statistical and the second systematic.
• ### The Low Energy Constants of $SU(2)$ Partially Quenched Chiral Perturbation Theory from $N_{f}=2+1$ Domain Wall QCD(1511.01950)

Nov. 5, 2015 hep-ph, hep-lat
We have performed fits of the pseudoscalar masses and decay constants, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at next-to leading order (NLO) and next-to-next-to leading order (NNLO). We report values for 9 NLO and 8 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the size of successive terms in the chiral expansion and use our large set of low energy constants to make predictions for mass splittings due to QCD isospin breaking effects and the S-wave $\pi \pi$ scattering lengths. We conclude that, for the range of pseudoscalar masses explored in this work, $115~\mathrm{MeV} \lesssim m_{\rm PS} \lesssim 430~\mathrm{MeV}$, the NNLO $SU(2)$ expansion is quite robust and can fit lattice data with percent-scale accuracy.
• ### Conserved currents for Mobius Domain Wall Fermions(1411.5728)

Nov. 20, 2014 hep-lat
We derive the exactly conserved vector, and almost conserved axial currents for rational approximations to the overlap operator with a general Mobius kernel. The approach maintains manifest Hermiticity, and allows matrix elements of the currents to be constructed at no extra cost after solution of the usual 5d system of equations, similar to the original approach of Furman and Shamir for domain wall Fermions.
• ### A Combined NNLO Lattice-Continuum Determination of $L_{10}^r$(1403.6729)

April 23, 2014 hep-ph, hep-ex, hep-lat
The renormalized next-to-leading-order (NLO) chiral low-energy constant, $L_{10}^r$, is determined in a complete next-to-next-to-leading-order (NNLO) analysis, using a combination of lattice and continuum data for the flavor $ud$ $V-A$ correlator and results from a recent chiral sum-rule analysis of the flavor-breaking combination of $ud$ and $us$ $V-A$ correlator differences. The analysis also fixes two combinations of NNLO low-energy constants, the determination of which is crucial to the precision achieved for $L_{10}^r$. Using the results of the flavor-breaking chiral $V-A$ sum rule obtained with current versions of the strange hadronic $\tau$ branching fractions as input, we find $L_{10}^r(m_\rho )\, =\, -0.00346(32)$. This result represents the first NNLO determination of $L_{10}^r$ having all inputs under full theoretical and/or experimental control, and the best current precision for this quantity.
• ### Lattice Input on the Inclusive $\tau$ Decay $V_{us}$ Puzzle(1312.1716)

Dec. 5, 2013 hep-ph, hep-ex, hep-lat
Recent analyses of flavor-breaking hadronic-$\tau$-decay-based sum rules produce values of $\vert V_{us}\vert$ $\sim 3\sigma$ low compared to 3-family unitarity expectations. An unresolved systematic issue is the significant variation in $\vert V_{us}\vert$ produced by different prescriptions for treating the slowly converging $D=2$ OPE series. We investigate the reliability of these prescriptions using lattice data for various flavor-breaking correlators and show the fixed-scale prescription is clearly preferred. Preliminary updates of the conventional $\tau$-based, and related mixed $\tau$-electroproduction-data-based, sum rule analyses incorporating B-factory results for low-multiplicity strange $\tau$ decay mode distributions are then performed. Use of the preferred FOPT $D=2$ OPE prescription is shown to significantly reduce the discrepancy between 3-family unitarity expectations and the sum rule results.
• ### $L_{10}^r$ From a Combined NNLO Lattice, Continuum Analysis of the Light Quark V-A Correlator(1311.0397)

Nov. 2, 2013 hep-ph, hep-lat
A combination of lattice and continuum data for the light-quark V-A correlator, supplemented by results from a chiral sum-rule analysis of the flavor-breaking flavor $ud$-$us$ V-A correlator difference, is shown to make possible a high-precision NNLO determination of the renormalized NLO chiral low-energy constant $L_{10}^r$. Key to this determination is the ability to simultaneously fix the two combinations of NNLO low-energy constants also entering the analysis. With current versions of the strange hadronic $\tau$ branching fractions required as input to the flavor-breaking V-A sum rule, we find $L_{10}^r(m_\rho ) = -0.00346(29)$. This represents both the best current precision for $L_{10}^r$, and the first NNLO determination having all errors under full control.
• ### New results from the lattice on the theoretical inputs to the hadronic tau determination of V_us(1301.4930)

Jan. 21, 2013 hep-ph, hep-ex, hep-lat
Recent sum rule determinations of |V_us|, employing flavor-breaking combinations of hadronic tau decay data, are significantly lower than either expectations based on 3-family unitarity or determinations from K_ell3 and Gamma[K_mu2]/Gamma[pi_mu2]. We use lattice data to investigate the accuracy/reliability of the OPE representation of the flavor-breaking correlator combination entering the tau decay analyses. The behavior of an alternate correlator combination, constructed to reduce problems associated with the slow convergence of the D = 2 OPE series, and entering an alternate sum rule requiring both electroproduction cross-section and hadronic tau decay data, is also investigated. Preliminary updates of both analyses, with the lessons learned from the lattice data in mind, are also presented.
• ### Some continuum physics results from the lattice V-A correlator(1301.2565)

Jan. 21, 2013 hep-ph, hep-lat
We present preliminary results on extractions of the chiral LECs L_10 and C_87 and constraints on the excited pseudoscalar state pi(1300) and pi(1800) decay constants obtained from an analysis of lattice data for the flavor ud light quark V-A correlator. A comparison of the results for the correlator to the corresponding mildly-model-dependent continuum results (based primarily on experimental hadronic tau decay data) is also given
• ### Lattice determination of the $K \to (\pi\pi)_{I=2}$ Decay Amplitude $A_2$(1206.5142)

June 22, 2012 hep-ph, hep-lat
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(\pi\pi)_{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 $\epsilon^\prime/\epsilon$ 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 $\epsilon^\prime/\epsilon$ is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.25 \pm 0.44_{\textrm{stat}} \pm 1.19_{\textrm{syst}}) \times 10^{-4}$.
• ### Opening the Rome-Southampton window for operator mixing matrices(1109.1223)

Jan. 24, 2012 hep-lat
We show that the running of operators which mix under renormalization can be computed fully non-perturbatively as a product of continuum step scaling matrices. These step scaling matrices are obtained by taking the "ratio" of Z matrices computed at different energies in an RI-MOM type scheme for which twisted boundary conditions are an essential ingredient. Our method allows us to relax the bounds of the Rome-Southampton window. We also explain why such a method is important in view of the light quark physics program of the RBC-UKQCD collaborations. To illustrate our method, using n_f=2+1 domain-wall fermions, we compute the non-perturbative running matrix of four-quark operators needed in K->pipi decay and neutral kaon mixing. Our results are then compared to perturbation theory.
• ### $K$ to $\pi\pi$ Decay amplitudes from Lattice QCD(1106.2714)

June 14, 2011 hep-ph, hep-lat
We report a direct lattice calculation of the $K$ to $\pi\pi$ decay matrix elements for both the $\Delta 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\pi\pi$ 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 $\Delta =1/2$ rule and testing the standard model of CP violation in the kaon system.
• ### Continuum Limit of $B_K$ from 2+1 Flavor Domain Wall QCD(1012.4178)

Dec. 19, 2010 hep-ph, hep-lat
We determine the neutral kaon mixing matrix element $B_K$ in the continuum limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action at two different lattice spacings. These lattice fermions have near exact chiral symmetry and therefore avoid artificial lattice operator mixing. We introduce a significant improvement to the conventional NPR method in which the bare matrix elements are renormalized non-perturbatively in the RI-MOM scheme and are then converted into the MSbar scheme using continuum perturbation theory. In addition to RI-MOM, we introduce and implement four non-exceptional intermediate momentum schemes that suppress infrared non-perturbative uncertainties in the renormalization procedure. We compute the conversion factors relating the matrix elements in this family of RI-SMOM schemes and MSbar at one-loop order. Comparison of the results obtained using these different intermediate schemes allows for a more reliable estimate of the unknown higher-order contributions and hence for a correspondingly more robust estimate of the systematic error. We also apply a recently proposed approach in which twisted boundary conditions are used to control the Symanzik expansion for off-shell vertex functions leading to a better control of the renormalization in the continuum limit. We control chiral extrapolation errors by considering both the NLO SU(2) chiral effective theory, and an analytic mass expansion. We obtain $B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_\chi(2)_{FV}(11)_{NPR}$. This corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_\chi(3)_{FV}(15)_{NPR}$. Adding all sources of error in quadrature we obtain $\hat{B}_K = 0.749(27)_{combined}$, with an overall combined error of 3.6%.
• ### Lattice Results for Low Moments of Light Meson Distribution Amplitudes(1011.5906)

Nov. 25, 2010 hep-lat
As part of the UKQCD and RBC collaborations' N_f=2+1 domain-wall fermion phenomenology programme, we calculate the first two moments of the light-cone distribution amplitudes of the pseudoscalar mesons pion and kaon and the (longitudinally-polarised) vector mesons rho, K-star and phi. We obtain the desired quantities with good precision and are able to discern the expected quark-mass dependence of SU(3)-flavour breaking effects. An important ingredient of the calculation is the nonperturbative renormalisation of lattice operators using the RI'/MOM technique.
• ### Light Meson Distribution Amplitudes(1011.2883)

Nov. 12, 2010 hep-lat
We calculated the first two moments of the light-cone distribution amplitudes for the pseudoscalar mesons ($\pi$ and $K$) and the longitudinally polarised vector mesons ($\rho$, $K^*$ and $\phi$) as part of the UKQCD and RBC collaborations' $N_f=2+1$ domain-wall fermion phenomenology programme. These quantities were obtained with a good precision and, in particular, the expected effects of $SU(3)$-flavour symmetry breaking were observed. Operators were renormalised non-perturbatively and extrapolations to the physical point were made, guided by leading order chiral perturbation theory. The main results presented are for two volumes, $16^3\times 32$ and $24^3\times 64$, with a common lattice spacing. Preliminary results for a lattice with a finer lattice spacing, $32^3\times64$, are discussed and a first look is taken at the use of twisted boundary conditions to extract distribution amplitudes.
• ### K->pi form factors with reduced model dependence(1004.0886)

April 6, 2010 hep-lat
• ### Non-perturbative renormalization of quark bilinear operators and B_K using domain wall fermions(0712.1061)

Dec. 6, 2007 hep-lat
We present a calculation of the renormalization coefficients of the quark bilinear operators and the K-Kbar mixing parameter B_K. The coefficients relating the bare lattice operators to those in the RI/MOM scheme are computed non-perturbatively and then matched perturbatively to the MSbar scheme. The coefficients are calculated on the RBC/UKQCD 2+1 flavor dynamical lattice configurations. Specifically we use a 16^3 x 32 lattice volume, the Iwasaki gauge action at beta=2.13 and domain wall fermions with L_s=16.
• ### 2+1 flavor domain wall QCD on a (2 fm)^3 lattice: light meson spectroscopy with Ls = 16(hep-lat/0701013)

Oct. 23, 2007 hep-lat
We present results for light meson masses and pseudoscalar decay constants from the first of a series of lattice calculations with 2+1 dynamical flavors of domain wall fermions and the Iwasaki gauge action. The work reported here was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice, which amounts to a spatial volume of (2 fm)^3 in physical units. The number of sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these simulations. Three values of input light sea quark masses, m_l^{sea} \approx 0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the physical light quark limit, whilst the heavier sea quark mass was fixed to approximately the physical strange quark mass m_s. The exact rational hybrid Monte Carlo algorithm was used to evaluate the fractional powers of the fermion determinants in the ensemble generation. We have found that f_\pi = 127(4) MeV, f_K = 157(5) MeV and f_K/f_\pi = 1.24(2), where the errors are statistical only, which are in good agreement with the experimental values.
• ### Lattice Results for Vector Meson Couplings and Parton Distribution Amplitudes(0710.0869)

Oct. 3, 2007 hep-lat
We present results for the couplings of light vector mesons to vector and tensor currents and on the low moments of the light-cone distribution amplitudes of the pion and kaon. The calculations are performed on the RBC and UKQCD collaborations' ensembles generated with the Iwasaki gauge action and with 2+1 flavours of domain wall fermions. The (preliminary) results for the ratios of the couplings of the vector meson to the vector and tensor currents (f_V and f^T_V respectively) in the MSbar scheme at 2GeV are: f_\rho^T/f_\rho=0.681(20); f_{K^\ast}^T/f_{K^\ast}=0.712(11) and f_\phi^T/f_\phi=0.751(9). For the first moment of the kaon's distribution amplitude we find (in the same scheme and at the same scale) <xi>_K=0.029(2) and for the second moment <xi^2>_\pi=0.28(3) and <xi^2>_K=0.27(2).