BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

C. T. Sachrajda

C. T. Sachrajda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phhep-thnucl-th

Catalog footprint

What is connected

43works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2022arXiv

Lattice determination of $I= 0$ and 2 $ππ$ scattering phase shifts with a physical pion mass

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Virtual Photon Emission in Leptonic Decays of Charged Pseudoscalar Mesons

We study the radiative leptonic decays $P\to\ellν_\ell\,\ell^{\prime\,+}\ell^{\prime\,-}$, where $P$ is a pseudoscalar meson and $\ell$ and $\ell^\prime$ are charged leptons. In such decays the emitted photon is off-shell and, in addition to the "point-like" contribution in which the virtual photon is emitted either from the lepton or the meson treated as a point-like particle, four structure-dependent (SD) form factors contribute to the amplitude. We present a strategy for the extraction of the SD form factors and implement it in an exploratory lattice computation of the decay rates for the four channels of kaon decays ($\ell,\ell^\prime=e,μ$). It is the SD form factors which describe the interaction between the virtual photon and the internal hadronic structure of the decaying meson, and in our procedure we separate the SD and point-like contributions to the amplitudes. We demonstrate that the form factors can be extracted with good precision and, in spite of the unphysical quark masses used in our simulation ($m_π\simeq 320\,$MeV and $m_K\simeq 530\,$MeV), the results for the decay rates are in reasonable semiquantitative agreement with experimental data (for the channels where these exist). Following this preparatory work, the emphasis of our future work will be on obtaining results at physical quark masses and on the control of the systematic uncertainties associated with discretisation and finite-volume errors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

First lattice calculation of radiative leptonic decay rates of pseudoscalar mesons

We present a non-perturbative lattice calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar ν_\ell γ$, where $P$ is a pseudoscalar meson and $\ell$ is a charged lepton. Together with the non-perturbative determination of the corrections to the processes $P\to \ell \bar ν_\ell$ due to the exchange of a virtual photon, this allows accurate predictions at $O(α_{em})$ to be made for leptonic decay rates for pseudoscalar mesons ranging from the pion to the $D_s$ meson. We are able to separate unambiguously and non-pertubatively the point-like contribution, from the structure-dependent, infrared-safe, terms in the amplitude. The fully non-perturbative $O(a)$ improved calculation of the inclusive leptonic decay rates will lead to the determination of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements also at $O(α_{em})$. Prospects for a precise evaluation of leptonic decay rates with emission of a hard photon are also very interesting, especially for the decays of heavy $D$ and $B$ mesons for which currently only model-dependent predictions are available to compare with existing experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Domain wall QCD with physical quark masses

We present results for several light hadronic quantities ($f_π$, $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_π$, $m_K$ and $m_Ω$ 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_π$ = 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).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electromagnetic corrections to the leptonic decay rates of charged pseudoscalar mesons: lattice results

Electromagnetic effects in the leptonic decay rates $π^+ \to μ^+ ν$ and $K^+ \to μ^+ ν$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Erratum: Standard-model prediction for direct CP violation in $K\toππ$ decay

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Prospects for a lattice computation of rare kaon decay amplitudes: $K\toπ\ell^+\ell^-$ decays

The rare kaon decays $K\toπ\ell^+\ell^-$ and $K\toπν\barν$ 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π\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πν\barν$ decays; these contributions are expected to be at the level of a few percent for $K^+$ decays.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Review of lattice results concerning low-energy particle physics

We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Standard-model prediction for direct CP violation in $K\toππ$ decay

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

$K \rightarrow ππ$ $ΔI=3/2$ decay amplitude in the continuum limit

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

QED Corrections to Hadronic Processes in Lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Low Energy Constants of $SU(2)$ Partially Quenched Chiral Perturbation Theory from $N_{f}=2+1$ Domain Wall QCD

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 $ππ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

$K_L-K_S$ mass difference from lattice QCD

We report on the first complete calculation of the $K_L-K_S$ mass difference, $Δ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 $Δ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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Domain Wall QCD with Near-Physical Pions

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_π=127(3)_{stat}(3)_{sys} MeV and f_K = 152(3)_{stat}(2)_{sys} MeV respectively (f_K/f_π= 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}.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Finite-volume effects in the evaluation of the K_L - K_S mass difference

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Determination of the $A_2$ amplitude of $K \rightarrow ππ$ decays

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 $ΔI=1/2$ rule.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kaon Mixing Beyond the Standard Model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Long distance contribution to the $K_L-K_S$ mass difference

We develop and demonstrate techniques needed to compute the long distance contribution to the $K_{L}$-$K_{S}$ mass difference, $ΔM_K$, in lattice QCD and carry out a first, exploratory calculation of this fundamental quantity. The calculation is performed on 2+1 flavor, domain wall fermion, $16^3\times32$ configurations with a 421 MeV pion mass and an inverse lattice spacing $1/a=1.73$ GeV. We include only current-current operators and drop all disconnected and double penguin diagrams. The short distance part of the mass difference in a 2+1 flavor calculation contains a quadratic divergence cut off by the lattice spacing. Here, this quadratic divergence is eliminated through the GIM mechanism by introducing a valence charm quark. The inclusion of the charm quark makes the complete calculation accessible to lattice methods provided the discretization errors associated with the charm quark can be controlled. The long distance effects are discussed for each parity channel separately. While we can see a clear signal in the parity odd channel, the signal to noise ratio in the parity even channel is exponentially decreasing as the separation between the two weak operators increases. We obtain a mass difference $ΔM_K$ which ranges from $6.58(30)\times 10^{-12}$ MeV to $11.89(81)\times 10^{-12}$ MeV for kaon masses varying from 563 MeV to 839 MeV. Extensions of these methods are proposed which promise accurate results for both $ΔM_K$ and $ε_K$, including long distance effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The $B^*Bπ$ coupling with relativistic heavy quarks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice determination of the $K \to (ππ)_{I=2}$ Decay Amplitude $A_2$

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

$K$ to $ππ$ Decay amplitudes from Lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Continuum Limit Physics from 2+1 Flavor Domain Wall QCD

We present physical results obtained from simulations using 2+1 flavors of domain wall quarks and the Iwasaki gauge action at two values of the lattice spacing $a$, ($a^{-1}$=\,1.73\,(3)\,GeV and $a^{-1}$=\,2.28\,(3)\,GeV). On the coarser lattice, with $24^3\times 64\times 16$ points, the analysis of ref.[1] is extended to approximately twice the number of configurations. The ensembles on the finer $32^3\times 64\times 16$ lattice are new. We explain how we use lattice data obtained at several values of the lattice spacing and for a range of quark masses in combined continuum-chiral fits in order to obtain results in the continuum limit and at physical quark masses. We implement this procedure at two lattice spacings, with unitary pion masses in the approximate range 290--420\,MeV (225--420\,MeV for partially quenched pions). We use the masses of the $π$ and $K$ mesons and the $Ω$ baryon to determine the physical quark masses and the values of the lattice spacing. While our data are consistent with the predictions of NLO SU(2) chiral perturbation theory, they are also consistent with a simple analytic ansatz leading to an inherent uncertainty in how best to perform the chiral extrapolation that we are reluctant to reduce with model-dependent assumptions about higher order corrections. Our main results include $f_π=124(2)_{\rm stat}(5)_{\rm syst}$\,MeV, $f_K/f_π=1.204(7)(25)$ where $f_K$ is the kaon decay constant, $m_s^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(96.2\pm 2.7)$\,MeV and $m_{ud}^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(3.59\pm 0.21)$\,MeV\, ($m_s/m_{ud}=26.8\pm 1.4$) where $m_s$ and $m_{ud}$ are the mass of the strange-quark and the average of the up and down quark masses respectively, $[Σ^{\msbar}(2 {\rm GeV})]^{1/3} = 256(6)\; {\rm MeV}$, where $Σ$ is the chiral condensate, the Sommer scale $r_0=0.487(9)$\,fm and $r_1=0.333(9)$\,fm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The $K\to(ππ)_{I=2}$ Decay Amplitude from Lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

$K\to(ππ)_{I=2}$ decays and twisted boundary conditions

We propose a new method to evaluate the Lellouch-Lüscher factor which relates the $ΔI=3/2$ $K\toππ$ matrix elements computed on a finite lattice to the physical (infinite-volume) decay amplitudes. The method relies on the use of partially twisted boundary conditions, which allow the s-wave $ππ$ phase shift to be computed as an almost continuous function of the centre-of-mass relative momentum and hence for its derivative to be evaluated. We successfully demonstrate the feasibility of the technique in an exploratory computation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Continuum Limit of $B_K$ from 2+1 Flavor Domain Wall QCD

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)_χ(2)_{FV}(11)_{NPR}$. This corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_χ(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%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

K->pi form factors with reduced model dependence

Using partially twisted boundary conditions we compute the K->pi semi-leptonic form factors in the range of momentum transfers 0 <~ q^2 <= q^2_{max}=(mK-mpi)^2 in lattice QCD with N_f=2+1 dynamical flavours. In this way we are able to determine f+(0) without any interpolation in the momentum transfer, thus eliminating one source of systematic error. This study confirms our earlier phenomenological ansatz for the strange quark mass dependence of the scalar form factor. We identify and estimate potentially significant NNLO effects in the chiral expansion that guides the extrapolation of the data to the physical point. Our main result is f+(0) = 0.9599(34)(^{+31}_{-43})(14)$, where the first error is statistical, the second error is due to the uncertainties in the chiral extrapolation of the lattice data and the last error is an estimate of potential discretisation effects.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Light Meson Distribution Amplitudes

We calculated the first two moments of the light-cone distribution amplitudes for the pseudoscalar mesons ($π$ and $K$) and the longitudinally polarised vector mesons ($ρ$, $K^*$ and $ϕ$) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The $η$ and $η^\prime$ mesons from Lattice QCD

The large mass of the ninth pseudoscalar meson, the $η^\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 $η$ and $η^\prime$ masses and mixing which confirms this picture. The physical eigenstates show small octet-singlet mixing with a mixing angle of $θ= -14.1(2.8)^\circ$. Extrapolation to physical light quark mass gives, with statistical errors only, $m_η=573(6)$ MeV and $m_{η^\prime}=947(142)$ MeV, consistent with the experimental values of 548 MeV and 958 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Continuum results for light hadrons from 2+1 flavor DWF ensembles

From simultaneous fits to the data from the 2+1 flavor DWF ensembles generated by the RBC and UKQCD collaborations at two different lattice spacings, we present preliminary continuum results for light hadrons. We focus on light pseudoscalar decay constants and quark masses. Several approaches to the calculation of the lattice spacing are discussed and the errors associated with the chiral extrapolation are explained. We make use of reweighting in the dynamical strange quark mass such that our ensembles have a self-consistently determined strange quark mass.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Kl3 and pion form factors using partially twisted boundary conditions

We compute the Kl3 and pion form factors using partially twisted boundary conditions. The twists are chosen so that the Kl3 form factors are calculated directly at zero momentum transfer (q^2=0), removing the need for a q^2 interpolation, while the pion form factor is determined at values of q^2 close to q^2=0. The simulations are performed on an ensemble of the RBC/UKQCD collaboration's gauge configurations with Domain Wall Fermions and the Iwaski gauge action with an inverse lattice spacing of 1.73(3) GeV. Simulating at a single pion mass of 330 MeV, we find the pion charge radius to be < r^2>_{330 MeV}=0.354(31) fm^2 which, using NLO SU(2) chiral perturbation theory, translates to a value of <r_π^2>=0.418(31) fm^2 for a physical pion. For the value of the Kl3 form factor, f_{Kπ}^+(q^2), determined directly at q^2=0, we find a value of f_{Kπ}^+(0)=0.9742(41) at this particular quark mass, which agrees well with our earlier result (0.9774(35)) obtained using the standard, indirect method.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Matrix elements of DI=3/2 K-->pi pi decays

We present a numerical computation of matrix elements of DI=3/2 K-->pi pi decays by using Wilson fermions. In order to extrapolate to the physical point we work at unphysical kinematics and we resort to Chiral Perturbation Theory at the next-to-leading order. In particular we explain the case of the electroweak penguins O_{7,8} which can contribute significantly in the theoretical prediction of epsilon'/epsilon. The study is done at beta=6.0 on a 24^3x64 lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Kaon Weak Matrix Elements with Wilson Fermions

We present results of several numerical studies with Wilson fermions relevant for kaon physics. We compute the B_K parameter by using two different methods and extrapolate to the continuum limit. Our preliminary result is B_K(2 GeV)=0.66(7). Delta I=3/2 K->pi pi matrix elements are obtained by using the next-to-leading order expressions derived in chiral perturbation theory in which the low energy constants are determined by the lattice results computed at unphysical kinematics. From the simulation at beta=6.0 our (preliminary) results read: <pi pi|O_7(2 GeV)|K>_{I=2}=0.14(1)(1) GeV^3 and <pi pi|O_8(2GeV)|K>_{I=2}=0.69(6)(6) GeV^3.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Extraction of K-->pipi Matrix Elements with Wilson Fermions

We present the status of a lattice calculation for the K-->pipi matrix elements of the (delta S=1) effective weak Hamiltonian, directly with two pion in the final state. We study the energy shift of two pion in a finite volume both in the I=0 and I=2 channels. We explain a method to avoid the Goldstone pole contamination in the computation of renormalization constants for (delta I=3/2) operators. Finally we show some preliminary results for the matrix elements of (delta I=1/2) operators. Our quenched simulation is done at beta=6.0, with Wilson fermions, on a (24^3 X 64) lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Finite Volume Effects in Weak Hadronic Decays

In this talk we discuss finite-volume computations of two-body hadronic decays below the inelastic threshold (e.g. $K\toππ$ decays). In particular we show how the relation between finite-volume matrix elements and physical amplitudes, recently derived by Lellouch and Lüscher, can be extended to all elastic states under the inelastic threshold. We also provide a derivation of the Lüscher quantization condition directly in quantum field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

QCD Factorization in B -> pi K, pi pi Decays and Extraction of Wolfenstein Parameters

In the heavy-quark limit, the hadronic matrix elements entering nonleptonic $B$-meson decays into two light mesons can be calculated from first principles including ``nonfactorizable'' strong-interaction corrections. The $B\toπK,ππ$ decay amplitudes are computed including electroweak penguin contributions, SU(3) violation in the light-cone distribution amplitudes, and an estimate of power corrections from chirally-enhanced terms and annihilation graphs. The results are then used to reduce the theoretical uncertainties in determinations of the weak phases $γ$ and $α$. In that way, new constraints in the $(\barρ,\barη)$ plane are derived. Predictions for the $B\toπK, ππ$ branching ratios and CP asymmetries are also presented. A good global fit to the (in part preliminary) experimental data on the branching fractions is obtained without taking recourse to phenomenological models.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Heavy Quark Physics From Lattice QCD

We review the application of lattice QCD to the phenomenology of b- and c-quarks. After a short discussion of the lattice techniques used to evaluate hadronic matrix elements and the corresponding systematic uncertainties, we summarise results for leptonic decay constants, B--Bbar mixing, semileptonic and rare radiative decays. A discussion of the determination of heavy quark effective theory parameters is followed by an explanation of the difficulty in applying lattice methods to exclusive nonleptonic decays.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Renormalons and the Heavy Quark Effective Theory

We propose a non-perturbative method for defining the higher dimensional operators which appear in the Heavy Quark Effective Theory (HQET), such that their matrix elements are free of renormalon singularities, and diverge at most logarithmically with the ultra-violet cut-off. Matrix elements of these operators can be computed numerically in lattice simulations of the HQET. We illustrate our procedures by presenting physical definitions of the binding energy ($\lb$) and of the kinetic energy (-$λ_1/2m_Q$) of the heavy quark in a hadron. This allows us to define a ``subtracted pole mass", whose inverse can be used as the expansion parameter in applications of the HQET.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Invisible Renormalon

We study the structure of renormalons in the Heavy Quark Effective Theory, by expanding the heavy quark propagator in powers of $1/m_Q$. We demonstrate that the way in which renormalons appear depends on the regularisation scheme used to define the effective theory. In order to investigate the relation between ultraviolet renormalons and power divergences of matrix elements of higher-dimensional operators in the heavy quark expansion, we perform calculations in dimensional regularisation and in three different cut-off regularisation schemes. In the case of the kinetic energy operator, we find that the leading ultraviolet renormalon which corresponds to a quadratic divergence, is absent in all but one (the lattice) regularisation scheme. The nature of this ``invisible renormalon'' remains unclear.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Quenched Hadrons using Wilson and O(a)-Improved Fermion Actions at beta=6.2

We present the first study of the light hadron spectrum and decay constants for quenched QCD using an O(a)-improved nearest-neighbour Wilson fermion action at β=6.2. We compare the results with those obtained using the standard Wilson fermion action, on the same set of 18 gauge field configurations of a 24^3 times 48 lattice. For pseudoscalar meson masses in the range 330-800 MeV, we find no significant difference between the results for the two actions. The scales obtained from the string tension and mesonic sector are consistent, but differ from that derived from baryon masses. The ratio of the pseudoscalar decay constant to the vector meson mass is roughly independent of quark mass as observed experimentally, and in approximate agreement with the measured value.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

The Hyperfine Splitting in Charmonium: Lattice Computations Using the Wilson and Clover Fermion Actions

We compute the hyperfine splitting $m_{J/ψ}-m_{η_c}$ on the lattice, using both the Wilson and $O(a)$-improved (clover) actions for quenched quarks. The computations are performed on a $24^3\times48$ lattice at $β= 6.2$, using the same set of 18 gluon configurations for both fermion actions. We find that the splitting is 1.83\err{13}{15} times larger with the clover action than with the Wilson action, demonstrating the sensitivity of the spin-splitting to the magnetic moment term which is present in the clover action. However, even with the clover action the result is less than half of the physical mass-splitting. We also compute the decay constants $f_{η_c}$ and $f^{-1}_{J/ψ}$, both of which are considerably larger when computed using the clover action than with the Wilson action. For example for the ratio $f^{-1}_{J/ψ}/f^{-1}_ρ$ we find 0.32\err{1}{2} with the Wilson action and $0.48\pm 3$ with the clover action (the physical value is 0.44(2)).

0 citations0 reviewsNo code linkedOpen access