No Arabic abstract
We describe details of 1/a ~ 2.2Gev, L ~ 3 fm dynamical domain wall fermion simulations which will allow us to do a more systematic continuum extrapolation in combination with existing simu- lations. Details of the simulations such as algorithm choices and machine performance, as well as results of basic measurements are presented. These configurations are presently being generated on the QCDOC machine at Edinburgh and the DOE QCDOC machine at Brookhaven as part of a joint project with LHPC.
We present renormalization constants of overlap quark bilinear operators on 2+1-flavor domain wall fermion configurations. This setup is being used by the chiQCD collaboration in calculations of physical quantities such as strangeness in the nucleon and the strange and charm quark masses. The scale independent renormalization constant for the axial vector current is computed using the Ward Identity. The renormalization constants for scalar, pseudoscalar and vector current are calculated in the RI-MOM scheme. Results in the MS-bar scheme are also given. The step scaling function of quark masses in the RI-MOM scheme is computed as well. The analysis uses, in total, six different ensembles of three sea quarks each on two lattices with sizes 24^3x64 and 32^3x64 at spacings a=(1.73 GeV)^{-1} and (2.28 GeV)^{-1}, respectively.
Nucleon-structure calculations of isovector vector- and axialvector-current form factors, transversity and scalar charge, and quark momentum and helicity fractions are reported from two recent 2+1-flavor dynamical domain-wall fermions lattice-QCD ensembles generated jointly by the RIKEN-BNL-Columbia and UKQCD Collaborations with Iwasaki $times$ dislocation-suppressing-determinatn-ratio gauge action at inverse lattice spacing of 1.378(7) GeV and pion mass values of 249.4(3) and 172.3(3) MeV.
The overlap fermion propagator is calculated on 2+1 flavor domain wall fermion gauge configurations on 16^3 x 32, 24^3 x 64 and 32^3 x 64 lattices. With HYP smearing and low eigenmode deflation, it is shown that the inversion of the overlap operator can be expedited by ~ 20 times for the 16^3 x 32 lattice and ~ 80 times for the 32^3 x 64 lattice. Through the study of hyperfine splitting, we found that the O(m^2a^2) error is small and these dynamical fermion lattices can adequately accommodate quark mass up to the charm quark. The low energy constant Delta_{mix} which characterizes the discretization error of the pion made up of a pair of sea and valence quarks in this mixed action approach is calculated via the scalar correlator with periodic and anti-periodic boundary conditions. It is found to be small which shifts a 300 MeV pion mass by ~ 10 to 19 MeV on these sets of lattices. We have studied the signal-to-noise issue of the noise source for the meson and baryon. It is found that the many-to-all meson and baryon correlators with Z_3 grid source and low eigenmode substitution is efficient in reducing errors for the correlators of both mesons and baryons. With 64-point Z_3 grid source and low-mode substitution, it can reduce the statistical errors of the light quark (m_{pi} ~ 200 - 300 MeV) meson and nucleon correlators by a factor of ~ 3-4 as compared to the point source. The Z_3 grid source itself can reduce the errors of the charmonium correlators by a factor of ~ 3.
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^3times 64times 16$ points, the analysis of ref.[1] is extended to approximately twice the number of configurations. The ensembles on the finer $32^3times 64times 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 $pi$ and $K$ mesons and the $Omega$ 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_pi=124(2)_{rm stat}(5)_{rm syst}$,MeV, $f_K/f_pi=1.204(7)(25)$ where $f_K$ is the kaon decay constant, $m_s^{bar{textrm{MS}}}(2,textrm{GeV})=(96.2pm 2.7)$,MeV and $m_{ud}^{bar{textrm{MS}}}(2,textrm{GeV})=(3.59pm 0.21)$,MeV, ($m_s/m_{ud}=26.8pm 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, $[Sigma^{msbar}(2 {rm GeV})]^{1/3} = 256(6); {rm MeV}$, where $Sigma$ is the chiral condensate, the Sommer scale $r_0=0.487(9)$,fm and $r_1=0.333(9)$,fm.
We report our numerical lattice QCD calculations of the isovector nucleon form factors for the vector and axialvector currents: the vector, induced tensor, axialvector, and induced pseudoscalar form factors. The calculation is carried out with the gauge configurations generated with N_f=2+1 dynamical domain wall fermions and Iwasaki gauge actions at beta = 2.13, corresponding to a cutoff 1/a = 1.73 GeV, and a spatial volume of (2.7 fm)^3. The up and down quark masses are varied so the pion mass lies between 0.33 and 0.67 GeV while the strange quark mass is about 12% heavier than the physical one. We calculate the form factors in the range of momentum transfers, 0.2 < q^2 < 0.75 GeV^2. The vector and induced tensor form factors are well described by the conventional dipole forms and result in significant underestimation of the Dirac and Pauli mean-squared radii and the anomalous magnetic moment compared to the respective experimental values. We show that the axialvector form factor is significantly affected by the finite spatial volume of the lattice. In particular in the axial charge, g_A/g_V, the finite volume effect scales with a single dimensionless quantity, m_pi L, the product of the calculated pion mass and the spatial lattice extent. Our results indicate that for this quantity, m_pi L > 6 is required to ensure that finite volume effects are below 1%.