• ### First lattice calculation of the QED corrections to leptonic decay rates(1711.06537)

Dec. 27, 2017 hep-ph, hep-ex, hep-lat
The leading-order electromagnetic and strong isospin-breaking corrections to the ratio of $K_{\mu 2}$ and $\pi_{\mu 2}$ decay rates are evaluated for the first time on the lattice, following a method recently proposed. The lattice results are obtained using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks. Systematics effects are evaluated and the impact of the quenched QED approximation is estimated. Our result for the correction to the tree-level $K_{\mu 2} / \pi_{\mu 2}$ decay ratio is $-1.22\,(16) \%$ to be compared to the estimate $-1.12\,(21) \%$ based on Chiral Perturbation Theory and adopted by the Particle Data Group.
• ### Electromagnetic corrections to leptonic decay rates of charged pseudoscalar mesons: finite-volume effects(1612.00199)

Dec. 1, 2016 hep-ph, hep-lat
In Carrasco et al. we have recently proposed a method to calculate $O(e^2)$ electromagnetic corrections to leptonic decay widths of pseudoscalar mesons. The method is based on the observation that the infrared divergent contributions (that appear at intermediate stages of the calculation and that cancel in physical quantities thanks to the Bloch-Nordsieck mechanism) are universal, i.e. depend on the charge and the mass of the meson but not on its internal structure. In this talk we perform a detailed analysis of the finite-volume effects associated with our method. In particular we show that also the leading $1/L$ finite-volume effects are universal and perform an analytical calculation of the finite-volume leptonic decay rate for a point-like meson.
• ### Finite-Volume QED Corrections to Decay Amplitudes in Lattice QCD(1611.08497)

Nov. 25, 2016 hep-lat
We demonstrate that the leading and next-to-leading finite-volume effects in the evaluation of leptonic decay widths of pseudoscalar mesons at $O(\alpha)$ are universal, i.e. they are independent of the structure of the meson. This is analogous to a similar result for the spectrum but with some fundamental differences, most notably the presence of infrared divergences in decay amplitudes. The leading non-universal, structure-dependent terms are of $O(1/L^2)$ (compared to the $O(1/L^3)$ leading non-universal corrections in the spectrum). We calculate the universal finite-volume effects, which requires an extension of previously developed techniques to include a dependence on an external three-momentum (in our case, the momentum of the final state lepton). The result can be included in the strategy proposed in Ref.\,\cite{Carrasco:2015xwa} for using lattice simulations to compute the decay widths at $O(\alpha)$, with the remaining finite-volume effects starting at order $O(1/L^2)$. The methods developed in this paper can be generalised to other decay processes, most notably to semileptonic decays, and hence open the possibility of a new era in precision flavour physics.
• ### Electromagnetic corrections to the leptonic decay rates of charged pseudoscalar mesons: lattice results(1610.09668)

Oct. 30, 2016 hep-ph, hep-lat
Electromagnetic effects in the leptonic decay rates $\pi^+ \to \mu^+ \nu$ and $K^+ \to \mu^+ \nu$ are evaluated for the first time on the lattice. Following a method recently proposed in Ref. [1] the emission of virtual photons at leading order in the electromagnetic coupling is evaluated on the lattice and the infrared divergence computed for a point-like meson at finite lattice volume is subtracted. The physical decay rate is then obtained by adding the emission of real and virtual photons regularised with a photon mass. Using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks the feasibility of our approach is demonstrated. Preliminary results for the electromagnetic corrections to charged (neutral) pion and kaon masses as well as to the leptonic decay rates of charged pions and kaons are presented.
• ### 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).
• ### The B*Bpi coupling using relativistic heavy quarks(1506.06413)

April 29, 2016 hep-lat
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.
• ### Prospects for a lattice computation of rare kaon decay amplitudes: $K\to\pi\ell^+\ell^-$ decays(1507.03094)

Jan. 28, 2016 hep-ph, hep-lat
The rare kaon decays $K\to\pi\ell^+\ell^-$ and $K\to\pi\nu\bar{\nu}$ are flavor changing neutral current (FCNC) processes and hence promising channels with which to probe the limits of the standard model and to look for signs of new physics. In this paper we demonstrate the feasibility of lattice calculations of $K\to\pi\ell^+\ell^-$ decay amplitudes for which long-distance contributions are very significant. We show that the dominant finite-volume corrections (those decreasing as powers of the volume) are negligibly small and that, in the four-flavor theory, no new ultraviolet divergences appear as the electromagnetic current $J$ and the effective weak Hamiltonian $H_W$ approach each other. In addition, we demonstrate that one can remove the unphysical terms which grow exponentially with the range of the integration over the time separation between $J$ and $H_W$. We will now proceed to exploratory numerical studies with the aim of motivating further experimental measurements of these decays. Our work extends the earlier study by Isidori, Turchetti and Martinelli which focussed largely on the renormalization of ultraviolet divergences. In a companion paper we discuss the evaluation of the long-distance contributions to $K\to\pi\nu\bar{\nu}$ decays; these contributions are expected to be at the level of a few percent for $K^+$ decays.
• ### 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.
• ### QED Corrections to Hadronic Processes in Lattice QCD(1502.00257)

Feb. 1, 2015 hep-ph, hep-lat
In this paper, for the first time a method is proposed to compute electromagnetic effects in hadronic processes using lattice simulations. The method can be applied, for example, to the leptonic and semileptonic decays of light or heavy pseudoscalar mesons. For these quantities the presence of infrared divergences in intermediate stages of the calculation makes the procedure much more complicated than is the case for the hadronic spectrum, for which calculations already exist. In order to compute the physical widths, diagrams with virtual photons must be combined with those corresponding to the emission of real photons. Only in this way do the infrared divergences cancel as first understood by Bloch and Nordsieck in 1937. We present a detailed analysis of the method for the leptonic decays of a pseudoscalar meson. The implementation of our method, although challenging, is within reach of the present lattice technology.
• ### $K_L-K_S$ mass difference from lattice QCD(1406.0916)

Sept. 23, 2014 hep-lat
We report on the first complete calculation of the $K_L-K_S$ mass difference, $\Delta M_K$, using lattice QCD. The calculation is performed on a 2+1 flavor, domain wall fermion ensemble with a 330MeV pion mass and a 575 MeV kaon mass. We use a quenched charm quark with a 949 MeV mass to implement Glashow-Iliopoulos-Maiani cancellation. For these heavier-than-physical particle masses, we obtain $\Delta M_K =3.19(41)(96)\times 10^{-12}$ MeV, quite similar to the experimental value. Here the first error is statistical and the second is an estimate of the systematic discretization error. An interesting aspect of this calculation is the importance of the disconnected diagrams, a dramatic failure of the OZI rule.
• ### Finite-volume effects in the evaluation of the K_L - K_S mass difference(1401.1362)

Jan. 7, 2014 hep-ph, hep-lat
The RBC and UKQCD collaborations have recently proposed a procedure for computing the K_L-K_S mass difference. A necessary ingredient of this procedure is the calculation of the (non-exponential) finite-volume corrections relating the results obtained on a finite lattice to the physical values. This requires a significant extension of the techniques which were used to obtain the Lellouch-Luscher factor, which contains the finite-volume corrections in the evaluation of non-leptonic kaon decay amplitudes. We review the status of our study of this issue and, although a complete proof is still being developed, suggest the form of these corrections for general volumes and a strategy for taking the infinite-volume limit. The general result reduces to the known corrections in the special case when the volume is tuned so that there is a two-pion state degenerate with the kaon.
• ### The $B^*B\pi$ coupling with relativistic heavy quarks(1311.2251)

Nov. 10, 2013 hep-lat
We report on a calculation of the $B^*B\pi$ coupling in lattice QCD. The strong matrix element $\langle B \pi | B^*\rangle$ is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HM$\chi$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$\chi$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.
• ### 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}$.
• ### $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%.
• ### 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.
• ### The $\eta$ and $\eta^\prime$ mesons from Lattice QCD(1002.2999)

Oct. 20, 2010 hep-ph, hep-lat
The large mass of the ninth pseudoscalar meson, the $\eta^\prime$, is believed to arise from the combined effects of the axial anomaly and the gauge field topology present in QCD. We report a realistic, 2+1 flavor, lattice QCD calculation of the $\eta$ and $\eta^\prime$ masses and mixing which confirms this picture. The physical eigenstates show small octet-singlet mixing with a mixing angle of $\theta = -14.1(2.8)^\circ$. Extrapolation to physical light quark mass gives, with statistical errors only, $m_\eta=573(6)$ MeV and $m_{\eta^\prime}=947(142)$ MeV, consistent with the experimental values of 548 MeV and 958 MeV.
• ### 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.
• ### 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).
• ### Hadronic form factors in Lattice QCD at small and vanishing momentum transfer(hep-lat/0703005)

March 9, 2007 hep-lat
The introduction of partially twisted boundary conditions allows weak and electromagnetic form factors to be evaluated at specified values of the hadronic momenta (and hence momentum transfers) in lattice simulations. We present and demonstrate this technique for the computation of the K->pi semileptonic form factor at zero momentum transfer and for the electromagnetic form factor of the pion at arbitrarily small momentum transfers. These exploratory computations are carried out in full QCD with 3 flavours of sea quarks, but with only two values of m_u=m_d which limits our ability to perform the chiral extrapolations. The results should therefore be viewed primarily as a demonstration of the feasibility of the method. For the K-> pi form factor we compare the new technique to the conventional approach and for the pion form factor we assess our results for very small momentum transfer with the help of chiral perturbation theory.