Researcher profile

R. D. Mawhinney

R. D. Mawhinney 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.

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.

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

Meson screening masses from lattice QCD with two light and the strange quark

We present results for screening masses of mesons built from light and strange quarks in the temperature range of approximately between 140 MeV to 800 MeV. The lattice computations were performed with 2+1 dynamical light and strange flavors of improved (p4) staggered fermions along a line of constant physics defined by a pion mass of about 220 MeV and a kaon mass of 500 MeV. The lattices had temporal extents Nt = 4, 6 and 8 and aspect ratios of Ns / Nt \geq 4. At least up to a temperature of 140 MeV the pseudo-scalar screening mass remains almost equal to the corresponding zero temperature pseudo-scalar (pole) mass. At temperatures around 3Tc (Tc being the transition temperature) the continuum extrapolated pseudo-scalar screening mass approaches very close to the free continuum result of 2 πT from below. On the other hand, at high temperatures the vector screening mass turns out to be larger than the free continuum value of 2 πT. The pseudo-scalar and the vector screening masses do not become degenerate even for a temperature as high as 4Tc. Using these mesonic spatial correlation functions we have also investigated the restoration of chiral symmetry and the effective restoration of the axial symmetry. We have found that the vector and the axial-vector screening correlators become degenerate, indicating chiral symmetry restoration, at a temperature which is consistent with the QCD transition temperature obtained in previous studies. On the other hand, the pseudo-scalar and the scalar screening correlators become degenerate only at temperatures larger than 1.3Tc, indicating that the effective restoration of the axial symmetry takes place at a temperature larger than the QCD transition temperature.

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

Equation of State for physical quark masses

We calculate the QCD equation of state for temperatures corresponding to the transition region with physical mass values for two degenerate light quark flavors and a strange quark using an improved staggered fermion action (p4-action) on lattices with temporal extent N_tau=8. We compare our results with previous calculations performed at twice larger values of the light quark masses as well as with results obtained from a resonance gas model calculation. We also discuss the deconfining and chiral aspects of the QCD transition in terms of renormalized Polyakov loop, strangeness fluctuations and subtracted chiral condensate. We show that compared to the calculations performed at twice larger value of the light quark mass the transition region shifts by about 5 MeV toward smaller temperatures

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

The finite temperature QCD using 2+1 flavors of domain wall fermions at N_t = 8

We study the region of the QCD phase transition using 2+1 flavors of domain wall fermions (DWF) and a $16^3 \times 8$ lattice volume with a fifth dimension of $L_s = 32$. The disconnected light quark chiral susceptibility, quark number susceptibility and the Polyakov loop suggest a chiral and deconfining crossover transition lying between 155 and 185 MeV for our choice of quark mass and lattice spacing. In this region the lattice scale deduced from the Sommer parameter $r_0$ is $a^{-1} \approx 1.3$ GeV, the pion mass is $\approx 300$ MeV and the kaon mass is approximately physical. The peak in the chiral susceptibility implies a pseudo critical temperature $T_c = 171(10)(17)$ MeV where the first error is associated with determining the peak location and the second with our unphysical light quark mass and non-zero lattice spacing. The effects of residual chiral symmetry breaking on the chiral condensate and disconnected chiral susceptibility are studied using several values of the valence $L_s$.

preprint2007arXiv

2+1 flavor domain wall QCD on a (2 fm)^3 lattice: light meson spectroscopy with Ls = 16

We present results for light meson masses and pseudoscalar decay constants from the first of a series of lattice calculations with 2+1 dynamical flavors of domain wall fermions and the Iwasaki gauge action. The work reported here was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice, which amounts to a spatial volume of (2 fm)^3 in physical units. The number of sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these simulations. Three values of input light sea quark masses, m_l^{sea} \approx 0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the physical light quark limit, whilst the heavier sea quark mass was fixed to approximately the physical strange quark mass m_s. The exact rational hybrid Monte Carlo algorithm was used to evaluate the fractional powers of the fermion determinants in the ensemble generation. We have found that f_π= 127(4) MeV, f_K = 157(5) MeV and f_K/f_π= 1.24(2), where the errors are statistical only, which are in good agreement with the experimental values.

preprint2006arXiv

The Kaon B-parameter from Quenched Domain-Wall QCD

We present numerical results for the kaon B-parameter, B_K, determined in the quenched approximation of lattice QCD. Our simulations are performed using domain-wall fermions and the renormalization group improved, DBW2 gauge action which combine to give quarks with good chiral symmetry at finite lattice spacing. Operators are renormalized non-perturbatively using the RI/MOM scheme. We study scaling by performing the simulation on two different lattices with a^{-1} = 1.982(30) and 2.914(54) GeV. We combine this quenched scaling study with an earlier calculation of B_K using two flavors of dynamical, domain-wall quarks at a single lattice spacing to obtain B_K(MS,NDR,mu=2GeV)=0.563(21)(39)(30), were the first error is statistical, the second systematic (without quenching errors) and the third estimates the error due to quenching.

preprint2003arXiv

Hardware and software status of QCDOC

QCDOC is a massively parallel supercomputer whose processing nodes are based on an application-specific integrated circuit (ASIC). This ASIC was custom-designed so that crucial lattice QCD kernels achieve an overall sustained performance of 50% on machines with several 10,000 nodes. This strong scalability, together with low power consumption and a price/performance ratio of $1 per sustained MFlops, enable QCDOC to attack the most demanding lattice QCD problems. The first ASICs became available in June of 2003, and the testing performed so far has shown all systems functioning according to specification. We review the hardware and software status of QCDOC and present performance figures obtained in real hardware as well as in simulation.

preprint2000arXiv

QCDOC: A 10-teraflops scale computer for lattice QCD

The architecture of a new class of computers, optimized for lattice QCD calculations, is described. An individual node is based on a single integrated circuit containing a PowerPC 32-bit integer processor with a 1 Gflops 64-bit IEEE floating point unit, 4 Mbyte of memory, 8 Gbit/sec nearest-neighbor communications and additional control and diagnostic circuitry. The machine's name, QCDOC, derives from ``QCD On a Chip''.