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

R. D. Mawhinney

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

ResearcherUnclaimed source record
hep-lathep-phnucl-th

Catalog footprint

What is connected

28works
3topics
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

28 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
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

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

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

NLO and NNLO Low Energy Constants for $SU(2)$ Chiral Perturbation Theory

We have performed global fits of $f_π$ and $m_π$, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at 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 convergence of the expansion and use our large set of low energy constants to make predictions for the pion mass splitting 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
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

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

The QCD chiral transition, $\ua$ symmetry and the Dirac spectrum using domain wall fermions

We report on a study of the finite-temperature QCD transition region for temperatures between 139 and 196 MeV, with a pion mass of 200 MeV and two space-time volumes: $24^3\times8$ and $32^3\times8$, where the larger volume varies in linear size between 5.6 fm (at T=139 MeV) and 4.0 fm (at T=195 MeV). These results are compared with the results of an earlier calculation using the same action and quark masses but a smaller, $16^3\times8$ volume. The chiral domain wall fermion formulation with a combined Iwasaki and dislocation suppressing determinant ratio gauge action are used. This lattice action accurately reproduces the $\sua$ and $\ua$ symmetries of the continuum. Results are reported for the chiral condensates, connected and disconnected susceptibilities and the Dirac eigenvalue spectrum. We find a pseudo-critical temperature, $T_c$, of approximately 165 MeV consistent with previous results and strong finite volume dependence below $T_c$. Clear evidence is seen for $\ua$ symmetry breaking above $T_c$ which is quantitatively explained by the measured density of near-zero modes in accordance with the dilute instanton gas approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD phase transition with physical-mass, chiral quarks

We report on the first lattice calculation of the QCD phase transition using chiral fermions at physical values of the quark masses. This calculation uses 2+1 quark flavors, spatial volumes between (4 fm$)^3$ and (11 fm$)^3$ and temperatures between 139 and 196 MeV . Each temperature was calculated using a single lattice spacing corresponding to a temporal Euclidean extent of $N_t=8$. The disconnected chiral susceptibility, $χ_{\rm disc}$ shows a pronounced peak whose position and height depend sensitively on the quark mass. We find no metastability in the region of the peak and a peak height which does not change when a 5 fm spatial extent is increased to 10 fm. Each result is strong evidence that the QCD ``phase transition'' is not first order but a continuous cross-over for $m_π=135$ MeV. The peak location determines a pseudo-critical temperature $T_c = 155(1)(8)$ MeV. Chiral $SU(2)_L\times SU(2)_R$ symmetry is fully restored above 164 MeV, but anomalous $U(1)_A$ symmetry breaking is non-zero above $T_c$ and vanishes as $T$ is increased to 196 MeV.

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
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
preprint2012arXiv

The chiral transition and U(1)_A symmetry restoration from lattice QCD using Domain Wall Fermions

We present results on both the restoration of the spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)_A symmetry in finite temperature QCD at zero chemical potential using lattice QCD. We employ domain wall fermions on lattices with fixed temporal extent N_τ= 8 and spatial extent N_σ= 16 in a temperature range of T = 139 - 195 MeV, corresponding to lattice spacings of a \approx 0.12 - 0.18 fm. In these calculations, we include two degenerate light quarks and a strange quark at fixed pion mass m_π= 200 MeV. The strange quark mass is set near its physical value. We also present results from a second set of finite temperature gauge configurations at the same volume and temporal extent with slightly heavier pion mass. To study chiral symmetry restoration, we calculate the chiral condensate, the disconnected chiral susceptibility, and susceptibilities in several meson channels of different quantum numbers. To study U(1)_A restoration, we calculate spatial correlators in the scalar and pseudo-scalar channels, as well as the corresponding susceptibilities. Furthermore, we also show results for the eigenvalue spectrum of the Dirac operator as a function of temperature, which can be connected to both U(1)_A and chiral symmetry restoration via Banks-Casher relations.

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

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.

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

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

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

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

Baryon Number, Strangeness and Electric Charge Fluctuations in QCD at High Temperature

We analyze baryon number, strangeness and electric charge fluctuations as well as their correlations in QCD at high temperature. We present results obtained from lattice calculations performed with an improved staggered fermion action (p4-action) at two values of the lattice cut-off with almost physical up and down quark masses and a physical value for the strange quark mass. We compare these results, with an ideal quark gas at high temperature and a hadron resonance gas model at low temperature. We find that fluctuations and correlations are well described by the former already for temperatures about 1.5 times the transition temperature. At low temperature qualitative features of the lattice results are quite well described by a hadron resonance gas model. Higher order cumulants, which become increasingly sensitive to the light pions, however show deviations from a resonance gas in the vicinity of the transition temperature.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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''.

0 citations0 reviewsNo code linkedOpen access