• Observation of neutrinoless double beta decay, a lepton number violating process that has been proposed to clarify the nature of neutrino masses, has spawned an enormous world-wide experimental effort. Relating nuclear decay rates to high-energy, beyond the Standard Model (BSM) physics requires detailed knowledge of non-perturbative QCD effects. Using lattice QCD and taking advantage of effective field theory methods, we compute the model-independent leading-order matrix elements of short-range operators, which arise due to heavy BSM mediators, that contribute to this decay. Contributions from short-range operators may prove to be equally important to or even more important than those from long-range Majorana neutrino exchange.
  • 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.
  • 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).
  • We compute the hadronic matrix elements of the four-quark operators needed for the study of $K^0-{\bar K^0}$ mixing beyond the Standard Model. We have used $n_f=2+1$ flavours of domain wall fermion (DWF) at two values of the lattice spacing ($a\approx0.08$ and $a\approx0.11\text{fm}$) and with lightest unitary pion mass of $\approx300\,\text{MeV}$. Renormalisation is performed non-perturbatively and the impact of different intermediate momentum schemes is investigated.
  • 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.
  • We study the heavy quark mass dependence of selected observables constructed from heavy-light meson correlation functions in small-volume two-flavour lattice QCD after taking the continuum limit. The light quark mass is tuned to zero, whereas the range of available heavy quark masses $m_h$ covers a region extending from around the charm to beyond the bottom quark mass scale. This allows entering the asymptotic mass-scaling regime as $1/m_h \to 0$ and performing well-controlled extrapolations to the infinite-mass limit. Our results are then compared to predictions obtained in the static limit of continuum Heavy Quark Effective Theory (HQET), in order to verify non-perturbatively that HQET is an effective theory of QCD. While in general we observe a nice agreement at the few-% level, we find it to be less convincing for the small-volume pseudoscalar decay constant when perturbative matching is involved.
  • 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\pi\pi$ 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 $\Delta 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 $\epsilon^\prime/\epsilon$ is Re($\epsilon^\prime/\epsilon)_\textrm{EWP}=-(6.6\pm 1.0)\times 10^{-4}$.
  • 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.
  • We present physical results for a variety of light hadronic quantities obtained via a combined analysis of three 2+1 flavour domain wall fermion ensemble sets. For two of our ensemble sets we used the Iwasaki gauge action with beta=2.13 (a^-1=1.75(4) GeV) and beta=2.25 (a^-1=2.31(4) GeV) and lattice sizes of 24^3 x 64 and 32^3 x 64 respectively, with unitary pion masses in the range 293(5)-417(10) MeV. The extent L_s for the 5^th dimension of the domain wall fermion formulation is L_s=16 in these ensembles. In this analysis we include a third ensemble set that makes use of the novel Iwasaki+DSDR (Dislocation Suppressing Determinant Ratio) gauge action at beta = 1.75 (a^-1=1.37(1) GeV) with a lattice size of 32^3 x 64 and L_s=32 to reach down to partially-quenched pion masses as low as 143(1) MeV and a unitary pion mass of 171(1) MeV, while retaining good chiral symmetry and topological tunneling. We demonstrate a significant improvement in our control over the chiral extrapolation, resulting in much improved continuum predictions for the above quantities. The main results of this analysis include the pion and kaon decay constants, f_\pi=127(3)_{stat}(3)_{sys} MeV and f_K = 152(3)_{stat}(2)_{sys} MeV respectively (f_K/f_\pi = 1.199(12)_{stat}(14)_{sys}); the average up/down quark mass and the strange-quark mass in the MSbar-scheme at 3 GeV, m_{ud}(MSbar, 3 GeV) = 3.05(8)_{stat}(6)_{sys} MeV and m_s(MSbar, 3 GeV) = 83.5(1.7)_{stat}(1.1)_{sys}; the neutral kaon mixing parameter in the MSbar-scheme at 3 GeV, B_K(MSbar,3 GeV) = 0.535(8)_{stat}(13)_{sys}, and in the RGI scheme, \hat B_K = 0.758(11)_{stat}(19)_{sys}; and the Sommer scales r_1 = 0.323(8)_{stat}(4)_{sys} fm and r_0 = 0.480(10)_{stat}(4)_{sys} (r_1/r_0 = 0.673(11)_{stat}(3)_{sys}). We also obtain values for the SU(2) ChPT effective couplings, \bar{l_3} = 2.91(23)_{stat}(7)_{sys}$ and \bar{l_4} = 3.99(16)_{stat}(9)_{sys}.
  • We present a computation of B-meson decay constants from lattice QCD simulations within the framework of Heavy Quark Effective Theory for the b-quark. The next-to-leading order corrections in the HQET expansion are included non-perturbatively. Based on Nf=2 gauge field ensembles, covering three lattice spacings a (0.08-0.05)fm and pion masses down to 190MeV, a variational method for extracting hadronic matrix elements is used to keep systematic errors under control. In addition we perform a careful autocorrelation analysis in the extrapolation to the continuum and to the physical pion mass limits. Our final results read fB=186(13)MeV, fBs=224(14)MeV and fBs/fB=1.203(65). A comparison with other results in the literature does not reveal a dependence on the number of dynamical quarks, and effects from truncating HQET appear to be negligible.
  • We report our final estimate of the b-quark mass from $N_f=2$ lattice QCD simulations using Heavy Quark Effective Theory non-perturbatively matched to QCD at $O(1/m_h)$. Treating systematic and statistical errors in a conservative manner, we obtain $\overline{m}_{\rm b}^{\overline{\rm MS}}(2 {\rm GeV})=4.88(15)$ GeV after an extrapolation to the physical point.
  • 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.
  • We review the status of recent calculations by the RBC-UKQCD collaboration of the complex amplitude $A_2$, corresponding to the decay of a kaon to a two pion state with total isospin 2. In particular, we present preliminary results from two new ensembles: $48^3 \times 96$ with $a^{-1}=1.73$ GeV and $64^3 \times 128$ with $a^{-1}=2.3$ GeV, both at physical kinematics. Both ensembles were generated Iwasaki gauge action and domain wall fermion action with 2+1 flavours. These results, in comparison to our earlier ones on a $32^3$ DSDR lattice with $a^{-1}=1.36$ GeV, enable us to significantly reduce the discretization errors. The partial cancellation between the two dominant contractions contributing to Re($A_2$) has been confirmed and we believe that this cancellation is a major contribution to the $\Delta I=1/2$ rule.
  • 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.
  • We report on an ongoing calculation of hadronic matrix elements needed to parameterize K-Kbar mixing in generic BSM scenarios, using domain wall fermions (DWF) at two lattice spacings. Recent work by the SWME collaboration shows a significant disagreement with our previous results for two of these quantities. Since the origin of this disagreement is unknown, it is important to reduce the various uncertainties. In this work, we are using N_f=2+1 DWF with Iwasaki gauge action at inverse lattice spacings of 2.31 and 1.75 GeV, with multiple unitary pions on each ensemble, the lightest being 290 and 330 MeV on the finer and coarser of the two ensembles respectively. This extends previous work by the addition of a second lattice spacing (a^{-1}\approx 1.75 GeV). Renormalization is carried out non-perturbatively in the RI/MOM scheme and converted perturbatively to MSbar.
  • We report the final results of the ALPHA collaboration for some B-physics observables: $f_B$, $f_{B_s}$ and $m_b$. We employ CLS configurations with 2 flavors of $O(a)$ improved Wilson fermions in the sea and pion masses ranging down to 190 MeV. The b-quark is treated in HQET to order $1/m_b$. The renormalization, the matching and the improvement were performed non-perturbatively, and three lattice spacings reaching $a=0.048$ fm are used in the continuum extrapolation.
  • There has been much speculation as to the origin of the \Delta I = 1/2 rule (Re A_0/Re A_2 \simeq 22.5). We find that the two dominant contributions to the \Delta I=3/2, K \to \pi \pi{} 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_\pi{} \simeq 420 MeV (Re A_0/Re A_2=9.1(2.1)) and m_\pi{} \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 \Delta 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.
  • 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.
  • 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
  • The 2012 PDG reports a tension at the level of $3 \sigma$ between two exclusive determinations of $|V_{ub}|$. They are obtained by combining the experimental branching ratios of $B \to \tau \nu$ and $B \to \pi l \nu$ (respectively) with a theoretical computation of the hadronic matrix elements $\fB$ and the $B \to \pi$ form factor $f_+(q^2)$. To understand the tension, improved precision and a careful analysis of the systematics involved are necessary. We report the results of the ALPHA collaboration for $\fB$ from the lattice with 2 flavors of $O(a)$ improved Wilson fermions. We employ HQET, including $1/m_b$ corrections, with pion masses ranging down to $\approx$ 190 MeV. Renormalization and matching were performed non-perturbatively, and three lattice spacings reaching $a^{-1}\approx 4.1$ GeV are used in the continuum extrapolation. We also present progress towards a computation of $f_+(q^2)$, to directly compare two independent exclusive determinations of $|V_{ub}|$ with each other and with inclusive determinations. Additionally, we report on preliminary results for $\fBq{s}$, needed for the analysis of $B_s \to \mu^+\mu^-$.}
  • We present our analysis of B physics quantities using non-perturbatively matched Heavy Quark Effective Theory (HQET) in Nf= 2 lattice QCD on the CLS ensembles. Using all-to-all propagators, HYP-smeared static quarks, and the Generalized Eigenvalue Problem (GEVP) approach with a conservative plateau selection procedure, we are able to systematically control all sources of error. With significantly increased statistics compared to last year, our preliminary results are mb (mb) = 4.22(10)(4)z GeV for the MS b-quark mass, and fB = 193(9)stat (4)_\chi MeV and fBs = 219(12)stat MeV for the B-meson decay constants.
  • We compute the hadronic matrix elements of the four-quark operators needed for the study of neutral kaon mixing beyond the Standard Model (SM). We use nf=2+1 flavours of domain-wall fermions (DWF) which exhibit good chiral-flavour symmetry. The renormalization is performed non-perturbatively through the RI-MOM scheme and our results are converted perturbatively to MSbar. The computation is performed on a single lattice spacing a=0.086 fm with a lightest unitary pion mass of 290 MeV. The various systematic errors, including the discretisation effects, are estimated and discussed. Our results confirm a previous quenched study, where large ratios of non-SM to SM matrix elements were obtained.
  • 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}$.
  • 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.
  • We compute the renormalization factors of four-quark operators needed for the study of $K\to\pi\pi$ decay in the $\Delta I=3/2$ channel. We evaluate the Z-factors at a low energy scale ($\mu_0=1.145 \GeV$) using four different non-exceptional RI-SMOM schemes on a large, coarse lattice ($a\sim 0.14\fm$) on which the bare matrix elements are also computed. Then we compute the universal, non-perturbative, scale evolution matrix of these renormalization factors between $\mu_0$ and $3\GeV$. We give the numerical results for the different steps of the computation in two different non-exceptional lattice schemes, and the connection to $\msbar$ at $3\GeV$ is made using one-loop perturbation theory.