We present our result for the $Ktopipi$ decay amplitudes for both the $Delta I=1/2$ and $3/2$ processes with the improved Wilson fermion action. Expanding on the earlier works by Bernard {it et al.} and by Donini {it et al.}, we show that mixings wit
h four-fermion operators with wrong chirality are absent even for the Wilson fermion action for the parity odd process in both channels due to CPS symmetry. Therefore, after subtraction of an effect from the lower dimensional operator, a calculation of the decay amplitudes is possible without complications from operators with wrong chirality, as for the case with chirally symmetric lattice actions. As a first step to verify the possibility of calculations with the Wilson fermion action, we consider the decay amplitudes at an unphysical quark mass $m_K sim 2 m_pi$. Our calculations are carried out with $N_f=2+1$ gauge configurations generated with the Iwasaki gauge action and nonperturbatively $O(a)$-improved Wilson fermion action at $a=0.091,{rm fm}$, $m_pi=280,{rm MeV}$, and $m_K=580,{rm MeV}$ on a $32^3times 64$ ($La=2.9,{rm fm}$) lattice. For the quark loops in the penguin and disconnected contributions in the $I=0$ channel, the combined hopping parameter expansion and truncated solver method work very well for variance reduction. We obtain, for the first time with a Wilson-type fermion action, that ${rm Re}A_0 = 60(36) times10^{ -8},{rm GeV}$ and ${rm Im}A_0 =-67(56) times10^{-12},{rm GeV}$ for a matching scale $q^* =1/a$. The dependence on the matching scale $q^*$ for these values is weak.
We present results for the $Ktopipi$ decay amplitudes for both the $Delta I=1/2$ and $3/2$ channels. This calculation is carried out on 480 gauge configurations in $N_f=2+1$ QCD generated over 12,000 trajectories with the Iwasaki gauge action and non
-perturbatively $O(a)$-improved Wilson fermion action at $a=0.091,{rm fm}$, $m_pi=280,{rm MeV}$ and $m_K=580,{rm MeV}$ on a $32^3times 64$ ($La=2.9,{rm fm}$) lattice. For the quark loops in the Penguin and disconnected contributions in the $I=0$ channel, the combined hopping parameter expansion and truncated solver techniques work very well for variance reduction. We obtain, for the first time with a Wilson-type fermion action, that ${rm Re}A_0 = 60(36) times10^{ -8},{rm GeV}$ and ${rm Im}A_0 =-67(56) times10^{-12},{rm GeV}$ for a matching scale $q^* =1/a$. The dependence on the matching scale is weak.
We present the results of 1+1+1 flavor QCD+QED simulation at the physical point, in which the dynamical quark effects in QED and the up-down quark mass difference are incorporated by the reweighting technique. The physical quark masses together with
the lattice spacing are determined with $m_{pi^+}$, $m_{K^+}$, $m_{K^0}$ and $m_{Omega^-}$ as physical inputs. Calculations are carried out using a set of 2+1 flavor QCD configurations near the physical point generated by the non-perturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $beta=1.9$ on a $32^3times 64$ lattice. We evaluate the values of the up, down and strange quark masses individually with non-perturbative QCD renormalization.
We perform a lattice QCD study of the $rho$ meson decay from the $N_f=2+1$ full QCD configurations generated with a renormalization group improved gauge action and a non-perturbatively $O(a)$-improved Wilson fermion action. The resonance parameters,
the effective $rhotopipi$ coupling constant and the resonance mass, are estimated from the $P$-wave scattering phase shift for the isospin I=1 two-pion system. The finite size formulas are employed to calculate the phase shift from the energy on the lattice. Our calculations are carried out at two quark masses, $m_pi=410,{rm MeV}$ ($m_pi/m_rho=0.46$) and $m_pi=300,{rm MeV}$ ($m_pi/m_rho=0.35$), on a $32^3times 64$ ($La=2.9,{rm fm}$) lattice at the lattice spacing $a=0.091,{rm fm}$. We compare our results at these two quark masses with those given in the previous works using $N_f=2$ full QCD configurations and the experiment.
We present an evaluation of the quark mass renormalization factor for Nf=2+1 QCD. The Schroedinger functional scheme is employed as the intermediate scheme to carry out non-perturbative running from the low energy region, where renormalization of bar
e mass is performed on the lattice, to deep in the high energy perturbative region, where the conversion to the renormalization group invariant mass or the MS-bar scheme is safely carried out. For numerical simulations we adopted the Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy regions and three lattice spacings to take the continuum limit at each scale. The regularization independent step scaling function of the quark mass for the Nf=2+1 QCD is obtained in the continuum limit. Renormalization factors for the pseudo scalar density and the axial vector current are also evaluated for the same action and the bare couplings as two recent large scale Nf=2+1 simulations; previous work of the CP-PACS/JLQCD collaboration, which covered the up-down quark mass range heavier than $m_pisim 500$ MeV and that of PACS-CS collaboration for much lighter quark masses down to $m_pi=155$ MeV. The quark mass renormalization factor is used to renormalize bare PCAC masses in these simulations.
We present the results of the physical point simulation in 2+1 flavor lattice QCD with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $beta=1.9$ on a $32^3 times 64$ lattice. The physical quark masses toge
ther with the lattice spacing is determined with $m_pi$, $m_K$ and $m_Omega$ as physical inputs. There are two key algorithmic ingredients to make possible the direct simulation at the physical point: One is the mass-preconditioned domain-decomposed HMC algorithm to reduce the computational cost. The other is the reweighting technique to adjust the hopping parameters exactly to the physical point. The physics results include the hadron spectrum, the quark masses and the pseudoscalar meson decay constants. The renormalization factors are nonperturbatively evaluated with the Schr{o}dinger functional method. The results are compared with the previous ones obtained by the chiral extrapolation method.
I present derivation of Luschers finite size formula for the elastic $Npi$ and the $NN$ scattering system for several angular momenta from the relativistic quantum field theory.
We present the first results of the PACS-CS project which aims to simulate 2+1 flavor lattice QCD on the physical point with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Numerical simulations are carried out
at the lattice spacing of $a=0.0907(13)$fm on a $32^3times 64$ lattice with the use of the DDHMC algorithm to reduce the up-down quark mass. Further algorithmic improvements make possible the simulation whose ud quark mass is as light as the physical value. The resulting PS meson masses range from 702MeV down to 156MeV, which clearly exhibit the presence of chiral logarithms. An analysis of the PS meson sector with SU(3) ChPT reveals that the NLO corrections are large at the physical strange quark mass. In order to estimate the physical ud quark mass, we employ the SU(2) chiral analysis expanding the strange quark contributions analytically around the physical strange quark mass. The SU(2) LECs ${bar l}_3$ and ${bar l}_4$ are comparable with the recent estimates by other lattice QCD calculations. We determine the physical point together with the lattice spacing employing $m_pi$, $m_K$ and $m_Omega$ as input. The hadron spectrum extrapolated to the physical point shows an agreement with the experimental values at a few % level of statistical errors, albeit there remain possible cutoff effects. We also find that our results of $f_pi=134.0(4.2)$MeV, $f_K=159.4(3.1)$MeV and $f_K/f_pi=1.189(20)$ with the perturbative renormalization factors are compatible with the experimental values. For the physical quark masses we obtain $m_{rm ud}^msbar=2.527(47)$MeV and $m_{rm s}^msbar=72.72(78)$MeV extracted from the axial-vector Ward-Takahashi identity with the perturbative renormalization factors.
