Improved lattice actions for Kogut-Susskind quarks have been shown to improve rotational symmetry and flavor symmetry. In this work we find improved scaling behavior of the rho and nucleon masses expressed in units of a length scale obtained from the static quark potential, and better behavior of the Dirac operator in instanton backgrounds.
We present a study for the pion decay constant $f_pi$ in the quenched approximation to lattice QCD with the Kogut-Susskind (KS) quark action, with the emphasis given to the renormalization problems. Numerical simulations are carried out at the couplings $beta = 6.0$ and 6.2 on $32^3times 64$ and $48^3times 64$ lattices, respectively. The pion decay constant is evaluated for all KS flavors via gauge invariant and non-invariant axial vector currents with the renormalization constants calculated by both non-perturbative method and perturbation theory. We obtain $f_pi = 89(6)$ MeV in the continuum limit as the best value using the partially conserved axial vector current, which requires no renormalization. From a study for the other KS flavors we find that the results obtained with the non-perturbative renormalization constants are well convergent among the KS flavors in the continuum limit, confirming restoration of $rm SU(4)_A$ flavor symmetry, while perturbative renormalization still leaves an apparent flavor breaking effect even in the continuum limit.
Report is made of a systematic scaling study of the finite-temperature chiral phase transition of two-flavor QCD with the Kogut-Susskind quark action based on simulations on $L^3times4$ ($L$=8, 12 and 16) lattices at the quark mass of $m_q=0.075, 0.0375, 0.02$ and 0.01. Our finite-size data show that a phase transition is absent for $m_qgeq 0.02$, and quite likely also at $m_q=0.01$. The scaling behavior of susceptibilities as a function of $m_q$ is consistent with a second-order transition at $m_q=0$. However, the exponents deviate from the O(2) or O(4) values theoretically expected.
Recently, the Fermilab heavy-quark action was extended to include dimension-six and -seven operators in order to reduce the discretization errors. In this talk, we present results of the first numerical simulations with this action (the OK action), where we study the masses of the quarkonium and heavy-light systems. We calculate combinations of masses designed to test improvement and compare results obtained with the OK action to their counterparts obtained with the clover action. Our preliminary results show a clear improvement.
We present a non-perturbative calculation for the pion decay constant with quenched Kogut-Susskind quarks. Numerical simulations are carried out at $beta = 6.0$ and 6.2 with various operators extending over all flavors. The renormalization correction is applied for each flavor by computing non-perturbative renormalization constants, and it is compared with a perturbative calculation. We also study the behavior of $f_pi$ in the continuum limits for both non-perturbative and perturbative calculations. The results in the continuum limit is also discussed.
We present a polynomial Hybrid Monte Carlo (PHMC) algorithm as an exact simulation algorithm with dynamical Kogut-Susskind fermions. The algorithm uses a Hermitian polynomial approximation for the fractional power of the KS fermion matrix. The systematic error from the polynomial approximation is removed by the Kennedy-Kuti noisy Metropolis test so that the algorithm becomes exact at a finite molecular dynamics step size. We performed numerical tests with $N_f$$=$2 case on several lattice sizes. We found that the PHMC algorithm works on a moderately large lattice of $16^4$ at $beta$$=$5.7, $m$$=$0.02 ($m_{mathrm{PS}}/m_{mathrm{V}}$$sim$0.69) with a reasonable computational time.