Researcher profile

C. T. Sachrajda

C. T. Sachrajda contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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

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.

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&#39; $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.

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.

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&#39;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.

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&#39;/epsilon. The study is done at beta=6.0 on a 24^3x64 lattice.

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.

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&#39;&#39; 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.

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.

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