Source author record

C. Kelly

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 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.

preprint2018arXiv

Combinatorial Quantum Gravity: Emergence of Geometric Space from Random Graphs

We review and extend the recently proposed model of combinatorial quantum gravity. Contrary to previous discrete approaches, this model is defined on (regular) random graphs and is driven by a purely combinatorial version of Ricci curvature, the Ollivier curvature, defined on generic metric spaces equipped with a Markov chain. It dispenses thus of notions such as simplicial complexes and Regge calculus and is ideally suited to extend quantum gravity to combinatorial structures which have a priori nothing to do with geometry. Indeed, our results show that geometry and general relativity emerge from random structures in a second-order phase transition due to the condensation of cycles on random graphs, a critical point that defines quantum gravity non-perturbatively according to asymptotic safety. In combinatorial quantum gravity the entropy area law emerges naturally as a consequence of infinite-dimensional critical behaviour on networks rather than on lattices. We propose thus that the entropy area law is a signature of the random graph nature of space-(time) on the smallest scales.

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).

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.

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.

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

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.

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

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

preprint2012arXiv

Opening the Rome-Southampton window for operator mixing matrices

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.

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.

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.

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%.

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.

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.

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.