1+1+1 flavor QCD + QED simulation at the physical point

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.

Physical Point Simulation in 2+1 Flavor Lattice QCD

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 together 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.

2+1 flavor lattice QCD simulation with $O(a)$-improved Wilson quarks

We present simulation details and results for the light hadron spectrum in N f = 2 + 1 lattice QCD with the nonperturbatively O(a)-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at a lattice spacing of 0.09 fm on a (2.9fm)^3 box using the PACS-CS computer. We employ the Luschers domain-decomposed HMC algorithm with several improvements to reduce the degenerate up-down quark mass toward the physical value. So far the resulting pseudoscalar meson mass is ranging from 702MeV down to 156MeV. We discuss on the stability and the efficiency of the algorithm. The light harden spectrum extrapolated at the physical point is compared with the experimental values. We also present the values of the quark masses and the pseudoscalar meson decay constants.

2+1 Flavor Lattice QCD toward the Physical Point

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.

Light hadron spectrum with 2+1 flavor dynamical $O(a)$-improved Wilson quarks

We present preliminary results for the light harden spectrum in $N_f=2+1$ lattice QCD obtained with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at $beta=1.90$ on a $32^3 times 64$ lattice using the PACS-CS computer. We employ Luschers domain-decomposed HMC algorithm to reduce the up-down quark masses toward the physical value. The pseudoscalar meson masses range from 730 MeV down to 210 MeV. We compare the light harden spectrum extrapolated to the physical point with the experimental values.