Search

We investigate the finite size effect on the vector meson and the baryon sectors using a subset of the "PACS10" configurations which are generated, keeping the space-time volumes over (10 fm$)^4$ in 2+1 flavor QCD at the physical point. Comparing the results on (5.5 fm$)^4$ and (10.9 fm$)^4$ lattices the ground states of octet baryons , which are stable on the lattice, show no finite size effect within less than 0.5% level of statistical errors. For those of vector mesons, which are unstable on the lattice, we observe that the effective masses are well below the experimental resonance levels both on (5.5 fm$)^4$ and (10.9 fm$)^4$ lattices. For the decuplet baryon sector we have found that the time dependence of the effective mass looks quite similar to that for the vector meson sector including the $\Omega$ baryon channel. We discuss its origin due to a possible mixing with the nearby multihadron states. Since the $\Xi$ baryon mass can be determined with the smallest ambiguity among the vector meson and the baryon masses, we use it together with the pion and kaon masses as the physical inputs to determine the physical point., Comment: 9 pages, 8 figures: published version

We present our result for the $K\to\pi\pi$ decay amplitudes for both the $\Delta I=1/2$ and $3/2$ processes with the improved Wilson fermion action. In order to realize the physical kinematics, where the pions in the final state have finite momenta, we consider the decay process $K({\bf p}) \to \pi({\bf p}) + \pi({\bf 0})$ in the nonzero momentum frame with momentum ${\bf p}=(0,0,2\pi/L)$ on the lattice. 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}$ ($1/a=2.176\,{\rm GeV}$), $m_\pi=260\,{\rm MeV}$, and $m_K=570\,{\rm MeV}$ on a $48^3\times 64$ ($La=4.4\,{\rm fm}$) lattice. For these parameters the energy of the $K$ meson is set at that of two-pion in the final state. We obtain ${\rm Re}A_2 = 2.431(19) \times10^{ -8}\,{\rm GeV}$, ${\rm Re}A_0 = 51(28) \times10^{ -8}\,{\rm GeV}$, and $\epsilon'/\epsilon = 1.9(5.7) \times10^{-3}$ for a matching scale $q^* =1/a$ where the errors are statistical. The dependence on the matching scale $q^*$ of these values is weak. The systematic error arising from the renormalization factors is expected to be around $1.3\%$ for ${\rm Re}A_2$ and $11 \%$ for ${\rm Re}A_0$. Prospects toward calculations with the physical quark mass are discussed., Comment: LaTeX2e, 19 pages, 5 eps figures, uses revtex4 and graphicx. arXiv admin note: substantial text overlap with arXiv:1505.05289. Published in PRD

We investigate the finite size effect on pseudoscalar meson masses and decay constants using a subset of the "PACS10" configurations which are generated keeping the space-time volumes over (10 fm$)^4$ in 2+1 flavor QCD at the physical point. We have tried two kinds of analyses, fixing $\kappa$ values or measured axial Ward identity quark masses. Comparing the results on (5.4 fm$)^4$ and (10.8 fm$)^4$ lattices, we have found a sizable finite size effect on the pseudoscalar meson sector in the former analysis: a 2.1(8)%, 4.8(1.6)%, and 0.36(31)% finite size effect on $m_\pi$, $m_{\rm ud}$, and $f_\pi$, respectively, on the (5.4 fm$)^4$ lattice. For the latter analysis, the finite size effect on the pseudoscalar meson decay constants is 0.66(33)% for $f_\pi$, 0.26(13)% for $f_K$, and 0.40(32)% for $f_K/f_\pi$. These values with two-sigma error bars are consistent with the predictions from the full one-loop SU(3) chiral perturbation theory, which are 0.20% for $f_\pi$, 0.08% for $f_K$, and 0.13% for $f_K/f_\pi$. The finite size effect on the pseudoscalar meson masses is hardly detected under the current statistical precision., Comment: 7 pages, 10 figures, published version

We investigate the nature of the chiral phase transition in the massless two-flavor QCD using the renormalization group improved gauge action and the Wilson quark action on $32^3\times 16$, $24^3\times 12$, and $16^3\times 8$ lattices. We calculate the spacial and temporal propagators of the iso-triplet mesons in the pseudo-scalar ($PS$), scalar ($S$), vector ($V$) and axial-vector ($AV$) channels on the lattice of three sizes. We first verify that the RG scaling is excellently satisfied for all cases. This is consistent with the claim that the chiral phase transition is second order. Then we compare the spacial and temporal effective masses between the axial partners, i.e. $PS$ vs $S$ and $V$ vs $AV$, on each of the three size lattices. We find the effective masses of all of six cases for the axial partners agree remarkably. This is consistent with the claim that at least $Z_4$ subgroup of the $U_A(1)$ symmetry in addition to the $SU_A(2)$ symmetry is recovered at the chiral phase transition point., Comment: 4 figures, 1 table. arXiv admin note: text overlap with arXiv:1704.03134

We investigate the nature of the chiral phase transition in the massless two-flavor QCD using the renormalization group improved gauge action and the Wilson quark action on $32^3\times 16$, $24^3\times 12$, and $16^3\times 8$ lattices. Based on the renormalization group equation, we derive the scaling relation for the effective masses of mesons at the chiral phase transition point. If the chiral phase transition is second order, the effective masses as a function of the rescaled time/space do not depend on the lattice size and show the universal behavior. We find that our numerical simulations on the three sizes of lattices are excellently on the scaling curves, which is consistent with the second order phase transition., Comment: 9 pages, 4 figures, 2 tables

We generate $2+1$ flavor QCD configurations near the physical point on a $96^4$ lattice employing the 6-APE stout smeared Wilson clover action with a nonperturbative $c_{\rm SW}$ and the Iwasaki gauge action at $\beta=1.82$. The physical point is estimated based on the chiral perturbation theory using several data points generated by the reweighting technique from the simulation point, wherer $m_\pi$,$m_K$ and $m_\Omega$ are used as physical inputs. The physics results include the quark masses, the hadron spectrum, the pseudoscalar meson decay constants and nucleon sigma terms, using the nonperturbative renormalization factors evaluated with the Schrodinger functional method., Comment: 7 pages, 10 figures. Proceedings of the 33rd International Symposium on Lattice Field Theory, July 14-18, 2015, Kobe, Japan

We present the quark mass and axial current renormalization factors for the RG-improved Iwasaki gauge action and three flavors of the stout smeared $O(a)$-improved Wilson quark action. We employ $\alpha=0.1$ and $n_{\mathrm{step}}=6$ for the stout link smearing parameters and all links in the quark action are replaced with the smeared links. Using the Schr\"{o}dinger functional scheme we evaluate the renormalization factors at $\beta=1.82$ where large scale simulations are being carried out., Comment: 7 pages, 2 figures, PoS.cls, PoS(LATTICE 2015)271

We investigate the phase structure of the SU(3) gauge theory with $N_f=8$ by numerical simulations employing the massless Domain-Wall fermions.Our aim is to study directly the massless quark region, since it is the most important region to clarify the properties of conformal theories. When the number of flavor is within the conformal window, it is claimed recently with Wilson quarks that there is the conformal region at the small quark mass region in the parameter space in addition to the confining phase and the deconfining phase. We study the properties of the conformal region investing the spatial Polyakov loops and the temporal meson propagators. Our data imply that there is the conformal region, and a phase transition between the confining phase and the conformal region takes place. These results are consistent with the claim that the conformal window is between $7$ and $16$. Progress reports on other related studies are also presented., Comment: 7 pages, 12 figures, talk presented at Lattice 2015, PoS (LATTICE2015)

We investigate the properties of quarks and gluons above the chiral phase transition temperature $T_c,$ using the RG improved gauge action and the Wilson quark action with two degenerate quarks mainly on a $32^3\times 16$ lattice. In the one-loop perturbation theory, the thermal ensemble is dominated by the gauge configurations with effectively $Z(3)$ center twisted boundary conditions, making the thermal expectation value of the spatial Polyakov loop take a non-trivial $Z(3)$ center. This is in agreement with our lattice simulation of high temperature QCD. We further observe that the temporal propagator of massless quarks at extremely high temperature $\beta=100.0 \, (T \simeq10^{58} T_c)$ remarkably agrees with the temporal propagator of free quarks with the $Z(3)$ twisted boundary condition for $t/L_t \geq 0.2$, but differs from that with the $Z(3)$ trivial boundary condition. As we increase the mass of quarks $m_q$, we find that the thermal ensemble continues to be dominated by the $Z(3)$ twisted gauge field configurations as long as $m_q \le 3.0 \, T$ and above that the $Z(3)$ trivial configurations come in. The transition is essentially identical to what we found in the departure from the conformal region in the zero-temperature many-flavor conformal QCD on a finite lattice by increasing the mass of quarks. We argue that the behavior is consistent with the renormalization group analysis at finite temperature., Comment: 16 pages, 9 figures; 4 tables, an appendix added

We present our result for the $K\to\pi\pi$ 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 with 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^3\times 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., Comment: LaTeX2e, 31 pages, 12 eps figures, uses revtex4 and graphicx, published in PRD

We present results for the $K\to\pi\pi$ 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^3\times 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., Comment: 7 pages, 3 figures, Contribution to the 32nd International Symposium on Lattice Field Theory (Lattice 2014), 23-28 June 2014, Columbia University, New York, NY, USA

We present results of our trial calculation of the $K \to \pi\pi$ decay amplitudes with the improved Wilson fermion action. Calculations are carried out with $N_f=2+1$ gauge configurations generated 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=560\,{\rm MeV} (\sim 2 m_\pi)$ on a $32^3\times 64$ ($La=2.9 {\rm fm}$) lattice., Comment: 7 pages, 3 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany (Added references in reference table of v1.)

We investigate SU(3) gauge theories in four dimensions with Nf fundamental fermions, on a lattice using the Wilson fermion. Clarifying the vacuum structure in terms of Polyakov loops in spatial directions and properties of temporal propagators using a new method "local analysis", we conjecture that the "conformal region" exists together with the confining region and the deconfining region in the phase structure parametrized by beta and K, both in the cases of the large Nf QCD within the conformal window (referred as Conformal QCD) with an IR cutoff and small Nf QCD at T/Tc>1 with Tc being the chiral transition temperature (referred as High Temperature QCD). Our numerical simulation on a lattice of the size 16^3 x 64 shows the following evidence of the conjecture. In the conformal region we find the vacuum is the nontrivial Z(3) twisted vacuum modified by non-perturbative effects and temporal propagators of meson behave at large t as a power-law corrected Yukawa-type decaying form. The transition from the conformal region to the deconfining region or the confining region is a sharp transition between different vacua and therefore it suggests a first order transition both in Conformal QCD and in High Temperature QCD. Within our fixed lattice simulation, we find that there is a precise correspondence between Conformal QCD and High Temperature QCD in the temporal propagators under the change of the parameters Nf and T/Tc respectively. In particular, we find the correspondence between Conformal QCD with Nf = 7 and High Temperature QCD with Nf=2 at T ~ 2 Tc being in close relation to a meson unparticle model. From this we estimate the anomalous mass dimension gamma* = 1.2 (1) for Nf=7. We also show that the asymptotic state in the limit T/Tc --> infty is a free quark state in the Z(3) twisted vacuum., Comment: 48 pages, 47 figures and 7tables; two subsections added on "in relation to mass deformed CFT" and "finite size scaling"; a version to be published in PRD

We give a new perspective on the properties of quarks and gluons at finite temperature T in N_f = 2 ~ 6 QCD. We point out the existence of an IR fixed point for the gauge coupling constant at T>T_c (T_c is the chiral phase transition temperature). Based on this observation we predict theoretically and verify numerically that the correlation functions of a meson G(t) at T/T_c > 1 decay with a power-law corrected Yukawa-type decaying form, G(t)=c exp(-m t)/t^alpha in the "conformal region" defined by m < c Lambda_IR, where Lambda_IR is the IR cutoff, m is the characteristic scale of the spectrum in the meson cannel and c is a constant of order 1. The decaying form is the characteristics of conformal theories with an IR cutoff. We discuss in detail how the resulting hyper scaling relation of physical observables may modify the existing argument about the order of the chiral phase transition in the N_f=2 case., Comment: 7 pages, 5figures, 1 table; revised to make the argument more clear with new figures and a table

We give a new perspective on the dynamics of conformal theories realized in the SU(N) gauge theory, when the number of flavors N_f is within the conformal window. Motivated by the RG argument on conformal theories with a finite IR cutoff \Lambda_{IR}, we conjecture that the propagator of a meson G_H(t) on a lattice behaves at large t as a power-law corrected Yukawa-type decaying form G_H(t) = c_H \exp{(-m_H t)}/t^{\alpha_H} instead of the exponentially decaying form c_H\exp{(-m_H t)}, in the small quark mass region where m_H \le c \Lambda_{IR}: m_H is the mass of the ground state hadron in the channel H and c is a constant of order 1. The transition between the "conformal region" and the "confining region" is a first order transition. Our numerical results verify the predictions for the N_f=7 case and the N_f=16 case in the SU(3) gauge theory with the fundamental representation., Comment: 4 pages, 4 figures; minor modifications as published

We investigate the charmed baryon mass spectrum using the relativistic heavy quark action on 2+1 flavor PACS-CS configurations previously generated on $32^3 \times 64$ lattice. The dynamical up-down and strange quark masses are tuned to their physical values, reweighted from those employed in the configuration generation. At the physical point, the inverse lattice spacing determined from the $\Omega$ baryon mass gives $a^{-1}=2.194(10)$ GeV, and thus the spatial extent becomes $L = 32 a = 2.88(1)$ fm. Our results for the charmed baryon masses are consistent with experimental values, except for the mass of $\Xi_{cc}$, which has been measured by only one experimental group so far and has not been confirmed yet by others. In addition, we report values of other doubly and triply charmed baryon masses, which have never been measured experimentally., Comment: 12 pages, 9 figures. arXiv admin note: substantial text overlap with arXiv:1212.0073

We study the algorithmic optimization and performance tuning of the Lattice QCD clover-fermion solver for the K computer. We implement the L\"uscher's SAP preconditioner with sub-blocking in which the lattice block in a node is further divided to several sub-blocks to extract enough parallelism for the 8-core CPU SPARC64$^{\mathrm{TM}}$ VIIIfx of the K computer. To achieve a better convergence property we use the symmetric successive over-relaxation (SSOR) iteration with {\it locally-lexicographical} ordering for the sub-blocks in obtaining the block inverse. The SAP preconditioner is included in the single precision BiCGStab solver of the nested BiCGStab solver. The single precision part of the computational kernel are solely written with the SIMD oriented intrinsics to achieve the best performance of the \SPARC on the K computer. We benchmark the single precision BiCGStab solver on the three lattice sizes: $12^3\times 24$, $24^3\times 48$ and $48^3\times 96$, with fixing the local lattice size in a node at $6^3\times 12$. We observe an ideal weak-scaling performance from 16 nodes to 4096 nodes. The performance of a computational kernel exceeds 50% efficiency, and the single precision BiCGstab has $\sim26% susutained efficiency., Comment: 7 pages, 4 figures, Proceedings of The 30th International Symposium on Lattice Field Theory, June 24-29, 2012, Cairns, Australia

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^3\times 64$ lattice. We evaluate the values of the up, down and strange quark masses individually with non-perturbative QCD renormalization., Comment: 6 pages, 10 figres

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 $\rho\to\pi\pi$ 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^3\times 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., Comment: LaTeX2e, 24 pages, 9 eps figures, uses revtex4 and graphicx, published in PRD

We investigate the charm quark system using the relativistic heavy quark action on 2+1 flavor PACS-CS configurations previously generated on $32^3 \times 64$ lattice. The dynamical up-down and strange quark masses are set to the physical values by using the technique of reweighting to shift the quark hopping parameters from the values employed in the configuration generation. At the physical point, the lattice spacing equals $a^{-1}=2.194(10)$ GeV and the spatial extent $L=2.88(1)$ fm. The charm quark mass is determined by the spin-averaged mass of the 1S charmonium state, from which we obtain $m_{\rm charm}^{\msbar}(\mu = m_{\rm charm}^{\msbar}) = 1.260(1)(6)(35)$ GeV, where the errors are due to our statistics, scale determination and renormalization factor. An additional systematic error from the heavy quark is of order $\alpha_s^2 f(m_Q a)(a \Lambda_{QCD})$, which is estimated to be a percent level if the factor $f(m_Q a)$ analytic in $m_Q a$ is of order unity. Our results for the charmed and charmed-strange meson decay constants are $f_D=226(6)(1)(5)$ MeV, $f_{D_s}=257(2)(1)(5)$ MeV, again up to the heavy quark errors of order $\alpha_s^2 f(m_Q a)(a \Lambda_{QCD})$. Combined with the CLEO values for the leptonic decay widths, these values yield $|V_{cd}| = 0.205(6)(1)(5)(9)$, $|V_{cs}| = 1.00(1)(1)(3)(3)$, where the last error is on account of the experimental uncertainty of the decay widths., Comment: 16 pages, 12 figures

We present preliminary results on the $\rho$ meson decay width from $N_f=2+1$ full QCD configurations generated by PACS-CS Collaboration. The decay width is estimated from the $P$-wave scattering phase shift for the isospin $I=1$ two-pion system. The finite size formula presented by L\"uscher in the center of mass frame and its extension to non-zero total momentum frame by Rummukainen and Gottlieb are employed for the calculations of the phase shift. Our calculations are carried out at $m_\pi=410\ {\rm MeV}$ ($m_\pi/m_\rho=0.46$) and $a=0.091\ {\rm fm}$ on a $32^3\times 64$ ($La=2.9 {\rm fm}$) lattice., Comment: 7 pages, 4 figures, talk given at the XXVIII International Symposium on Lattice Filed Theory, June 14-19,2010, Villasimius, Sardinia Italy

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 bare 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_\pi\sim 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., Comment: 26 pages, 17 Postscript figures. Two tables are updated

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., Comment: 20 pages, 17 figures, version to appear in Phys. Rev. D

We present the International Lattice Data Grid (ILDG), a loosely federated grid of grids for sharing data from Lattice Quantum Chromodynamics (LQCD) simulations. The ILDG comprises of metadata, file format and web-service standards, which can be used to wrap regional data-grid interfaces, allowing seamless access to catalogues and data in a diverse set of collaborating regional grids. We discuss the technological underpinnings of the ILDG, primarily the metadata and the middleware, and offer a critique of its various aspects with the hindsight of the design work and the first full year of production., Comment: 13 pages

We present an evaluation of the running coupling constant for Nf=2+1 QCD. The Schroedinger functional scheme is used as the intermediate scheme to carry out non-perturbative running from the low energy region, where physical scale is introduced, to deep in the high energy perturbative region, where conversion to the MS-bar scheme is safely performed. Possible systematic errors due to the use of perturbation theory occur only in the conversion from three-flavor to four-flavor running coupling constant near the charm mass threshold, where higher order terms beyond 5th order in the $\beta$ function may not be negligible. For numerical simulations we adopted 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 region and three lattice spacings to take the continuum limit at each scale. A physical scale is introduced from the previous Nf=2+1 simulation of the CP-PACS/JL-QCD collaboration, which covered the up-down quark mass range heavier than $m_\pi\sim 500$ MeV., Comment: 20 pages, 11 Postscript figures; corrected for publication

We investigate the quark mass dependence of baryon masses in 2+1 flavor lattice QCD using SU(3) heavy baryon chiral perturbation theory up to one-loop order. The baryon mass data used for the analyses are obtained for the degenerate up-down quark mass of 3 MeV to 24 MeV and two choices of the strange quark mass around the physical value. We find that the SU(3) chiral expansion fails to describe both the octet and the decuplet baryon data if phenomenological values are employed for the meson-baryon couplings. The SU(2) case is also examined for the nucleon. We observe that higher order terms are controlled only around the physical point. We also evaluate finite size effects using SU(3) heavy baryon chiralperturbation theory, finding small values of order 1% even at the physical point., Comment: 26 pages, 10 tables, 58 figures

We study heavy-heavy and heavy-light quark systems for charm with a relativistic heavy quark action in 2+1 flavor lattice QCD. Configurations are generated by the PACS-CS Collaboration at the lattice spacing is $a=0.09$ fm with the lattice size of $32^3\times 64$ employing the $O(a)$-improved Wilson quark action and the Iwasaki gauge action. We present preliminary results for the charmonium spectrum and the $D$ and $D_s$ meson decay constants evaluated at 3.5 MeV$< m_{\rm ud}<$ 12 MeV with $m_{\rm s}$ around the physical value. We investigate the dynamical quark mass dependences of the hyperfine and the orbital splittings. The decay constants are compared with the recent experimental values., Comment: 7 Pages, 8 figures, talk presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

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 Luscher's 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., Comment: 7 pages, 15 figures, 1 table talk presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

We investigate the quark mass dependence of meson and baryon masses obtained from 2+1 flavor dynamical quark simulations performed by the PACS-CS Collaboration. With the use of SU(2) and SU(3) chiral perturbation theories up to NLO, we examine the chiral behavior of the pseudoscalar meson masses and the decay constants in terms of the degenerate up-down quark mass ranging form 3 MeV to 24 MeV and two choices of the strange quark mass around the physical value. We discuss the convergence of the SU(2) and SU(3) chiral expansions and present the results for the low energy constants. We find that the SU(3) expansion is not convergent at NLO for the physical strange quark mass. The chiral behavior of the nucleon mass is also discussed based on the SU(2) heavy baryon chiral perturbation theory up to NNLO. Our results show that the expansion is well behaved only up to m_pi^2 ~ 0.2 GeV^2., Comment: 7 Pages, 8 figures, talk presented at the XXVI International Symposium on Lattice Field Theory, July 14-19, 2008, Williamsburg, Virginia, USA

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^3\times 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., Comment: 43 pages, 48 figures

We perform numerical simulations of lattice QCD with two flavors of dynamical overlap quarks, which have exact chiral symmetry on the lattice. While this fermion discretization is computationally demanding, we demonstrate the feasibility to simulate reasonably large and fine lattices by a careful choice of the lattice action and algorithmic improvements. Our production runs are carried out on a 16^3 \times 32 lattice at a single lattice spacing around 0.12 fm. We explore the sea quark mass region down to m_s/6, where m_s is the physical strange quark mass, for a good control of the chiral extrapolation in future calculations of physical observables. We describe in detail our setup and algorithmic properties of the production simulations and present results for the static quark potential to fix the lattice scale and the locality of the overlap operator., Comment: 38 pages, 28 figures; v2: minor changes in text, references added, version published in PRD

We apply chiral perturbation theory to the pseudoscalar meson mass and decay constant data obtained in the PACS-CS Project toward 2+1 flavor lattice QCD simulations with the O(a)-improved Wilson quarks. We examine the existence of chiral logarithms in the quark mass range from m_{ud}=47 MeV down to 6 MeV on a (2.8 fm)^3 box with the lattice spacing a=0.09 fm. Several low energy constants are determined. We also discuss the magnitude of finite size effects based on chiral perturbation theory., Comment: 7 pages, 6 figures, talk presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg, Germany

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 L\"uscher's 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., Comment: 7 pages, 11 figures, talk presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4, 2007, Regensburg, Germany

QCDml is an XML-based markup language designed for sharing QCD configurations and ensembles world-wide via the International Lattice Data Grid (ILDG). Based on the latest release, we present key ingredients of the QCDml in order to provide some starting points for colleagues in this community to markup valuable configurations and submit them to the ILDG., Comment: 7 pages, 3 figures, poster presented at Lattice2007, Regensburg, Germany (Algorithms and machines)

We present a lattice QCD calculation of the $\rho$ meson decay width via the $P$-wave scattering phase shift for the I=1 two-pion system. Our calculation uses full QCD gauge configurations for $N_f=2$ flavors generated using a renormalization group improved gauge action and an improved Wilson fermion action on a $12^3\times24$ lattice at $m_\pi/m_\rho=0.41$ and the lattice spacing $1/a=0.92 {\rm GeV}$. The phase shift calculated with the use of the finite size formula for the two-pion system in the moving frame shows a behavior consistent with the existence of a resonance at a mass close to the vector meson mass obtained in spectroscopy. The decay width estimated from the phase shift is consistent with the experiment, when the quark mass is scaled to the realistic value., Comment: LaTeX2e, 16 pages, 6 eps figures, uses revtex4 and graphicx

We calculate the light meson spectrum and the light quark masses by lattice QCD simulation, treating all light quarks dynamically and employing the Iwasaki gluon action and the nonperturbatively O(a)-improved Wilson quark action. The calculations are made at the squared lattice spacings at an equal distance a^2~0.005, 0.01 and 0.015 fm^2, and the continuum limit is taken assuming an O(a^2) discretization error. The light meson spectrum is consistent with experiment. The up, down and strange quark masses in the \bar{MS} scheme at 2 GeV are \bar{m}=(m_{u}+m_{d})/2=3.55^{+0.65}_{-0.28} MeV and m_s=90.1^{+17.2}_{-6.1} MeV where the error includes statistical and all systematic errors added in quadrature. These values contain the previous estimates obtained with the dynamical u and d quarks within the error., Comment: 4 pages, 3 figures, revtex4; v2: contents partly modified, published version

We perform a numerical test of a relativistic heavy quark(RHQ) action, recently proposed by Tsukuba group, in quenched lattice QCD at $a\simeq 0.1$ fm. With the use of the improvement parameters previously determined at one-loop level for the RHQ action, we investigate a restoration of rotational symmetry for heavy-heavy and heavy-light meson systems around the charm quark mass. We focused on two quantities, the meson dispersion relation and the pseudo-scalar meson decay constants. It is shown that the RHQ action significantly reduces the discretization errors due to the charm quark mass. We also calculate the S-state hyperfine splittings for the charmonium and charmed-strange mesons and the $D_s$ meson decay constant. The remaining discretization errors in the physical quantities are discussed., Comment: 21 pages, 16 figures. A reference and a comment added, a major modification in appendix, several minor changes in the abstract and the main text. Errors in affiliation are corrected. Version appeared in JHEP

We discuss a possibility that the Neutron Electric Dipole Moment (NEDM) can be calculated in lattice QCD simulations in the presence of the CP violating $\theta$ term. In this paper we measure the energy difference between spin-up and spin-down states of the neutron in the presence of an uniform and static external electric field. We first test this method in quenched QCD with the RG improved gauge action on a $16^3\times 32$ lattice at $a^{-1}\simeq$ 2 GeV, employing two different lattice fermion formulations, the domain-wall fermion and the clover fermion for quarks, at relatively heavy quark mass $(m_{PS}/m_V \simeq 0.85)$. We obtain non-zero values of NEDM from calculations with both fermion formulations. We next consider some systematic uncertainties of our method for NEDM, using $24^3\times 32$ lattice at the same lattice spacing only with the clover fermion. We finally investigate the quark mass dependence of NEDM and observe a non-vanishing behavior of NEDM toward the chiral limit. We interpret this behavior as a manifestation of the pathology in the quenched approximation., Comment: LaTeX2e, 51 pages, 43 figures, uses revtex4 and graphicx, References and comments added, typos corrected, accepted by PRD

We present a study of $B$-parameters for generic $\Delta S=2$ four-quark operators in domain wall QCD. Our calculation covers all the $B$-parameters required to study the neutral kaon mixing in the standard model (SM) and beyond it. We evaluate one-loop renormalization factors of the operators employing the plaquette and Iwasaki gauge actions. Numerical simulations are carried out in quenched QCD with both gauge actions on $16^3\times 32\times 16$ and $24^3\times 32\times 16$ at the lattice spacing $1/a\approx 2$GeV. We investigate the relative magnitudes of the non-SM $B$-parameters to the SM one, which are compared with the previous results obtained with the overlap and the clover quark actions., Comment: 7 pages, talk presented at Lattice2006 (Electroweak Decays and Mixing), PoS format

We report on a study of 2+1 flavor lattice QCD with the $O(a)$-improved Wilson quarks on a $16^3\times 32$ lattice at the lattice spacing $1/a\approx 2$GeV employing Luescher's domain-decomposed HMC(LDDHMC) algorithm. This is dedicated to a preliminary study for the PACS-CS project which plans to complete the Wilson-clover $N_f=2+1$ program lowering the up-down quark masses close to the physical values as much as possible. We focus on three issues: (i) how light quark masses we can reach with LDDHMC, (ii) efficiency of the algorithm compared with the conventional HMC, (iii) parameter choice for the production runs on PACS-CS., Comment: 7 pages, talk presented at Lattice2006 (Algorithm), PoS format

We present a summary of results of the joint CP-PACS and JLQCD project toward a 2+1 flavor full QCD simulation with the O(a)-improved Wilson quark formalism and the Iwasaki gauge action. Configurations were generated during 2002-2005 at three lattice spacings, a~0.076, 0.100 and 0.122 fm, keeping the physical volume constant at (2.0fm)^3. Up and down quark masses are taken in the range m_{PS}/m_V~0.6-0.78. We have completed the analysis for the light meson spectrum and quark masses in the continuum limit using the full configuration set. The predicted meson masses reproduce experimental values in the continuum limit at a 1% level. The average up and down, and strange quark masses turn out to be m_{ud}^{\bar{MS}}(\mu=2 GeV)=3.50(14)({}^{+26}_{-15}) MeV and m_s^{\bar{MS}}(\mu=2 GeV)=91.8(3.9)({}^{+6.8}_{-4.1}) MeV. We discuss our future strategy toward definitive results on hadron spectroscopy with the Wilson-clover formalism., Comment: 7 pages, 4 figures, talk presented at Lattice2006 (Hadron Spectroscopy), PoS format (v3: typos corrected)

We apply the UV-filtering preconditioner, previously used to improve the Multi-Boson algorithm, to the Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm. The performance test for the algorithm is given for the plaquette gauge action and the $O(a)$-improved Wilson action at $\beta=5.2, c_{\mathrm{sw}}=2.02, M_{\pi}/M_{\rho}\sim 0.8$ and 0.7 on a $16^3\times 48$ lattice. We find that the UV-filtering reduces the magnitude of the molecular dynamics force from the pseudo fermion by a factor 3 by tuning the UV-filter parameter. Combining with the multi-time scale molecular dynamics integrator we achieve a factor 2 improvement., Comment: 7 pages, 1 figure, talk presented at Lattice2006 (Algorithm), PoS format; v2 references added

We present an overview of our project of large-scale simulations with dynamical overlap fermions. The first production run in two-flavor QCD is on-going using the Iwasaki gauge action on a 16^3 \times 32 lattice at the lattice spacing of 0.12fm with six sea quark masses down to m_{s,phys}/6, where m_{s,phys} is the physical strange quark mass. We briefly introduce our choice of the lattice action and simulation algorithm, and describe the present status of the production run. Preliminary results on the light meson masses and the static quark potential are also reported., Comment: 7 pages, 8 figures, talk presented at Lattice 2006 (Chiral Symmetry)

We report on a direct lattice calculation of the neutron EDM(NEDM) using the external electric field method in both quenched and full QCD. In quenched QCD, we use a $24^3\times 32$ lattice at $\beta=2.6$ with the Iwasaki gauge action and the clover fermion action to examine the viability of this method. In particular we investigate possible effects of violation of the periodic boundary condition of the external electric field on the NEDM signal. We also study the quark mass dependence of NEDM in quenched QCD, and observe that NEDM seems to remain non-zero toward the chiral limit because of the quenched artifact. In 2-flavor full QCD we employ configurations generated by the CP-PACS collaboration on a $24^3\times 48$ lattice at $\beta=2.1$ with the same gluon and quark actions as in the quenched case. Since the number of configurations is limited, we employ 8 different source points per one configuration and take an average over them. Our preliminary result at three quark masses ($m_{PS}/m_V\simeq 0.81, 0.76, 0.69$) indicates that non-zero value for NEDM can be obtained in full QCD. Statistical errors, however, are still too large to show the theoretically expected behaviour for NEDM in full QCD that it vanishes in the chiral limit., Comment: 7 pages, 15 figures, a talk given at the XXIV International Symposium on Lattice Field Theory (Lattice 2006, Hadronic Interactions and Structure), Tucson, Arizona, USA, July 23-28, 2006

Masses of the eta and eta-prime mesons are estimated in Nf=2+1 lattice QCD with the non-perturbatively O(a) improved Wilson quark action and the Iwasaki RG-improved gluon action, using CP-PACS/JLQCD configurations on a 16^3 x 32 lattice at beta=1.83 (lattice spacing is 0.122 fm). We apply a stochastic noise estimator technique combined with smearing method to evaluate correlators among flavor SU(2) singlet pseudoscalar operators and strange pseudoscalar operators for 10 combinations of up/down and strange quark masses. The correlator matrix is then diagonalized to identify signals for mass eigenstates. Masses of the ground state and the first excited state extrapolated to the physical point are m_eta= 0.545(16) GeV and m_eta-prime= 0.871(46) GeV, being close to the experimental values of the eta and eta-prime masses., Comment: 7 pages, 7 figures, talk presented at Lattice2006 (Hadron Spectroscopy)

We present preliminary results on the $\rho$ meson decay width estimated from the scattering phase shift of the I=1 two-pion system. The phase shift is calculated by the finite size formula for non-zero total momentum frame (the moving frame) derived by Rummukainen and Gottlieb, using the $N_f=2$ improved Wilson fermion action at $m_\pi/m_\rho=0.41$ and $L=2.53 {\rm fm}$., Comment: 7 pages, 6 figures, Talk at Lattice 2006(spectroscopy), to be published in PoS(LAT2006)110

We present a calculation of pion electromagnetic and scalar form factors in two-flavor QCD with the non-perturbatively O(a)-improved Wilson fermion. Chiral extrapolation of the corresponding charge radius is discussed based on the chiral perturbation theory., Comment: 6 pages, 4 figures, presented at Lattice 2005 (Weak matrix elements), Dublin, July 25-30, 2005, to appear in Proceedings of Science

We present a calculation of the kaon form factors in two-flavor QCD with the non-perturbatively $O(a)$-improved Wilson quark action. In order to achieve a few percent accuracy in the study of SU(3) breaking effects, we use a set of double ratios of the matrix elements, with which the bulk of the statistical fluctuation and the multiplicative renormalization factors cancel., Comment: 6 pages, 7 figures, talk presented at Lattice2005 (Weak matrix elements)

We describe our plan to develop a large-scale cluster system with a peak speed of 14.3Tflops for lattice QCD at the Center for Computational Sciences, University of Tsukuba, as a successor to the current 0.6Tflops CP-PACS computer. The system consist of 2560 nodes connected by a 16x16x10 three-dimensional hyper crossbar network. Each node has a single low-voltage 2.8GHz Xeon processor and 2GBytes of memory with 6.4GBytes/sec bandwidth, and 160 GBytes of disk in RAID1 mode. The network link in each of the three directions is made of dual Gigabit Ethernet with the peak throughput of 250MByte/sec. Hence each node has an aggregate network bandwidth of 750MByte/sec. The system will run under Linux and SCore, and an extension of the PM driver is developed for the network. The system will be developed jointly with Hitachi Limited. The installation is scheduled in the first quarter of Japanese Fiscal 2006 (April-June 2006) and the start of operation is expected in July 2006., Comment: 6 pages, 2 figures, talk presented at Lattice2005 (Algorithm and Machines)

We carry out a feasibility study toward a lattice QCD calculation of the neutron electric dipole moment (NEDM) in the presence of the $\theta$ term using two different approaches. In the first method, we calculate the CP-odd electromagnetic form factor $F_3$, which becomes the NEDM in the zero momentum transfer limit. At the first order in $\theta$, we derive a formula connecting the lattice three-point function to the CP-odd electromagnetic form factor. In the second method we directly extract the NEDM from the energy difference between spin-up and spin-down neutron states in the presence of a constant electric field, without expanding a small but non-zero $\theta$. We test both approaches numerically, employing the domain-wall quark action with the RG improved gauge action in quenched QCD at $a^{-1}\simeq 2$ GeV on a $16^3\times 32\times 16$ lattice, and further applying the second method to the clover quark action at a similar lattice spacing and nucleon mass. We obtain good signals from both approaches. In particular the second method works well with both fermion formulations., Comment: 6 pages, 7 eps-figures, talk presented at Lattice 2005 (Chiral fermions), Trinity College, Dublin, Ireland, 25-30 Jul 2005