Search

In the early days of QCD, the axial $U(1)$ anomaly was considered as a trigger for the breaking of the $SU(2)_L\times SU(2)_R$ symmetry through topological excitations of gluon fields. However, it has been a challenge for lattice QCD to quantify the effect. In this work, we simulate QCD at high temperatures with chiral fermions. The exact chiral symmetry enables us to separate the contribution from the axial $U(1)$ breaking from others among the susceptibilities in the scalar and pseudoscalar channels. Our result in two-flavor QCD indicates that the chiral susceptibility, which is conventionally used as a probe for $SU(2)_L\times SU(2)_R$ breaking, is actually dominated by the axial $U(1)$ breaking at temperatures $T\ge 165$ MeV., Comment: 9 pages, 3 figures, contribution to LATTICE 2021, published version

We calculate the form factors of the $K \to \pi l \nu$ semileptonic decays in three-flavor lattice QCD, and study their chiral behavior as a function of the momentum transfer and the Nambu-Goldstone boson masses. Chiral symmetry is exactly preserved by using the overlap quark action, which enables us to directly compare the lattice data with chiral perturbation theory (ChPT). We generate gauge ensembles at a lattice spacing of 0.11fm with four pion masses covering 290-540 MeV and a strange quark mass m_s close to its physical value. By using the all-to-all quark propagator, we calculate the vector and scalar form factors with high precision. Their dependence on m_s and the momentum transfer is studied by using the reweighting technique and the twisted boundary conditions for the quark fields. We compare the results for the semileptonic form factors with ChPT at next-to-next-to leading order in detail. While many low-energy constants appear at this order, we make use of our data of the light meson electromagnetic form factors in order to control the chiral extrapolation. We determine the normalization of the form factors as f_+(0) = 0.9636(36)(+57/-35), and observe reasonable agreement of their shape with experiment., Comment: 34 pages, 13 figures

We study the vector and axial-vector current correlators in perturbative and non-perturbative regimes of QCD. The correlators in Euclidean coordinate space are calculated on the lattice using the M\"obius domain-wall fermion formulation at three lattice spacings covering 0.044--0.080~fm. The dynamical quark effects of $2+1$ light flavors are included. The sum $V+A$ and the difference $V-A$ of the vector ($V$) and axial-vector ($A$) current correlators calculated on the lattice after extrapolating to the physical point agree with those converted from the ALEPH experimental data of hadronic $\tau$ decays. The level of the agreement in the $V+A$ channel is about $1.3\sigma$ or smaller in the region of $|x|\ge0.4$~fm, while that in the $V-A$ channel is about $1.8\sigma$ at $|x|=0.74$~fm and smaller at other distances. We also extract the hiral condensate from the short-distance correlators on the lattice using the PCAC relation. Its result extrapolated to the chiral and continuum limit is compatible with other estimates at low energies., Comment: 32 pages, 19 figures

We report on our study of the D meson semileptonic decays in 2+1 flavor lattice QCD. Gauge ensembles are generated at three lattice cutoffs up to 4.5 GeV and with pion masses as low as 300 MeV. We employ the Moebius domain-wall fermion action for both light and charm quarks. We report our preliminary results for the vector and scalar form factors and discuss their dependence on the momentum transfer, quark masses and lattice spacing., Comment: 7 pages, 4 figures, talk presented at the 34th International Symposium on Lattice Field Theory (Lattice 2016), 24-30 July 2016, University of Southampton, UK

We determine the renormalization constants for flavor non-singlet fermion bilinear operators of M\"obius domain-wall fermions. The renormalization condition is imposed on the correlation functions in the coordinate space, such that the non-perturbative lattice calculation reproduces the perturbatively calculated counterpart at short distances. The perturbative expansion is precise as the coefficients are available up to $O(\alpha_s^4)$. We employ $2+1$-flavor lattice ensembles at three lattice spacings in the range 0.044--0.080~fm., Comment: 36 pages, 20 figures

We discuss the fate of the axial U(1) symmetry in 2-flavor QCD at finite temperature, where the non-singlet chiral symmetry is restored. We first summarize the previous theoretical investigation on the relation between the eigenvalue density of the Dirac operator and the axial U(1) symmetry. We show that the eigenvalue density near the origin behaves as $\lambda^\gamma$ with $\gamma > 2$ in the chirally symmetric phase, where $\lambda$ is an eigenvalue. This implies that the axial U(1) symmetry is partially restored, so that the low energy symmetry of the theory becomes SU(2)$\otimes$ SU(2)$\otimes$ Z$_4$. Secondly, we report recent numerical investigations on this issue by lattice QCD simulations with lattice chiral fermions such as Overlap or improved domain-wall fermions. Our preliminary results indicate that the eigenvalue density seems to have a gap at the origin just above $T_c$, the temperature of the chiral symmetry restoration, which implies the axial U(1) symmetry is effectively restored at high temperature. We also point out an importance of the exact lattice chiral symmetry to obtain correct results on this issue., Comment: 10 pages, 5 figures, A talk given at "The 8th International Workshop on Chiral Dynamics, CD2015", 29 June - 03 July, 2015, Pisa, Italy

We present a study of chiral behavior of light meson form factors in QCD with three flavors of overlap quarks. Gauge ensembles are generated at single lattice spacing 0.12 fm with pion masses down to 300 MeV. The pion and kaon electromagnetic form factors and the kaon semileptonic form factors are precisely calculated using the all-to-all quark propagator. We discuss their chiral behavior using the next-to-next-to-leading order chiral perturbation theory., Comment: 7 pages, 4 figures, Proceedings of the 33rd International Symposium on Lattice Field Theory, 14-18 July 2015, Kobe, Japan

We study the $U(1)_A$ anomaly in two-flavor lattice QCD at finite temperature with the M\"obius domain-wall Dirac operator. We generate gauge configurations in the temperature range $(0.9, 1.2) T_c$ on different physical volumes, $L=$ 2--4 fm, and lattice spacings. We measure the difference of the susceptibilities of the flavor non-singlet scalar ($\chi_\delta$) and pseudoscalar ($\chi_\pi$) mesons. They are related by an axial $U(1)$ transformation and the difference vanishes if the axial $U(1)$ symmetry is respected. We identify the source of axial $U(1)$ symmetry breaking at finite temperature in the lowest eigenmodes, for the observable $\chi_\pi - \chi_\delta$. We then reweight the M\"obius domain-wall fermion partition function to that of the overlap-Dirac operator to fully recover chiral symmetry. Our data show a significant discrepancy in the results coming from the M\"obius domain-wall valence quarks, the overlap valence quarks on our DWF configurations and the reweighted ones that have full chiral symmetry. After recovering full chiral symmetry we conclude that the difference $\chi_\pi - \chi_\delta$ shows a suppression in the chiral limit that is compatible with an effective restoration of $U(1)_A$ at $T \gtrsim T_c$ in the scalar meson channels., Comment: 7 pages, 3 figures. Proceedings of the The 33rd International Symposium on Lattice Field Theory, 2015 Kobe. arXiv admin note: text overlap with arXiv:1412.5703

Calculation of the eta-prime meson mass is a notoriously difficult problem, as it requires evaluation of the disconnected diagram which is costly and noisy. In this work, we use a gluonic operator to extract the eta-prime state after smearing the link variables through the Wilson flow. With this choice, one can avoid a large cancellation of pion contribution between the connected and disconnected diagrams. We obtain the eta-prime meson mass on lattices with three different lattice spacings and two physical volumes, which allow us to estimate its continuum and large volume limits., Comment: 7 pages, 3 figures, Proceedings of the 33rd International Symposium on Lattice Field Theory, 14-18 July 2015, Kobe, Japan

We investigate the effects of the violation of the Ginsparg-Wilson (GW) relation in the M\"obius domain-wall fermion formulation on the lattice with finite fifth dimension. Using a decomposion in terms of the eigenmodes of its four-dimensional effective Dirac operator, we isolate the GW-violating terms for various physical quantities including the residual mass and the meson susceptibilities relevant for the effective restoration of the axial U(1) symmetry at finite temperature. Numerical result shows that the GW-violating effect is more significant, or even overwhelming, for the quantities that are dominated by the low-lying eigenmodes., Comment: 23 pages, 9 figures, 2 tables

We study the chiral behavior of the electromagnetic (EM) form factors of pion and kaon in three-flavor lattice QCD. In order to make a direct comparison of the lattice data with chiral perturbation theory (ChPT), we employ the overlap quark action that has exact chiral symmetry. Gauge ensembles are generated at a lattice spacing of 0.11 fm with four pion masses ranging between M_pi \simeq 290 MeV and 540 MeV and with a strange quark mass m_s close to its physical value. We utilize the all-to-all quark propagator technique to calculate the EM form factors with high precision. Their dependence on m_s and on the momentum transfer is studied by using the reweighting technique and the twisted boundary conditions for the quark fields, respectively. A detailed comparison with SU(2) and SU(3) ChPT reveals that the next-to-next-to-leading order terms in the chiral expansion are important to describe the chiral behavior of the form factors in the pion mass range studied in this work. We estimate the relevant low-energy constants and the charge radii, and find reasonable agreement with phenomenological and experimental results., Comment: 59 pages, 34 figures

The flavor-singlet component of the eta prime meson is related to the topological structure of the SU(3) gauge field through the chiral anomaly. We perform a 2+1-flavor lattice QCD calculation and demonstrate that the two-point function of a gluonically defined topological charge density after a short Yang-Mills gradient flow contains the propagation of the eta prime meson, by showing that its mass in the chiral and continuum limit is consistent with the experimental value. The gluonic correlator does not suffer from the contamination of the pion contribution, and the clean signal is obtained at significantly lower numerical cost compared to the conventional method with the quark bilinear operators., Comment: 11pages, 5figures, minor corrections, to appear in PRD

We investigate the axial $U(1)$ symmetry restoration at finite temperature in two flavor QCD. We employ the M\"obius domain-wall formalism that is designed to achieve good chiral symmetry. We show the measurements of a difference of meson susceptibilities, sensitive to the $U(1)_A$ symmetry breaking. The signal is dominated by zero and near-zero modes. By reweighting the measure to that of overlap fermions we find a suppression of the $U(1)_A$ breaking effects above the chiral transition temperature., Comment: 7 pages, 4 figures, contribution to the The 32nd International Symposium on Lattice Field Theory, 23-28 June, 2014, Columbia University New York, NY

Global topological charge decorrelates very slowly or even freezes in fine lattice simulations. On the other hand, its local fluctuations are expected to survive and lead to the correct physical results as long as the volume is large enough. We investigate this issue on recently generated configurations including dynamical domain-wall fermions at lattice spacings a = 0.08 fm and finer. We utilize the Yang-Mills gradient flow to define the topological charge density operator and calculate its long-distance correlation, through which we propose a new method for extracting the topological susceptibility in a sub-volume. This method takes care of the finite volume correction, which reduces the bias caused by the global topological charge. Our lattice data clearly show a shorter auto-correlation time than that of the naive definition using the whole lattice, and are less sensitive to the global topological history. Numerical results show a clear sea-quark mass dependence, which agrees well with the prediction of chiral perturbation theory., Comment: 7 pages, 4 figures, presented at the 32nd International Symposium on Lattice Field Theory (Lattice 2014), June 23-28 2014, New York, USA, references added

We carry out a comparative study among five-dimensional formulations of chirally symmetric fermions about the algorithmic performance, chiral symmetry violation and topological tunneling to find a computationally inexpensive formulation with good chiral symmetry. With our choice of the lattice action, we have launched large-scale simulations on fine lattices aiming at a precision study of light and heavy quark physics. We report on the comparative study, current status of the large-scale simulations, and preliminary results on the residual quark mass and auto-correlation., Comment: 7 pages, talk presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We calculate the kaon semileptonic form factors in lattice QCD with three flavors of dynamical overlap quarks. Gauge ensembles are generated at pion masses as low as 290 MeV and at a strange quark mass near its physical value. We precisely calculate relevant meson correlators using the all-to-all quark propagator. Twisted boundary conditions and the reweighting technique are employed to vary the momentum transfer and the strange quark mass. We discuss the chiral behavior of the form factors by comparing with chiral perturbation theory and experiments., Comment: 7 pages, 8 pdf figures, talk presented at the 30th International Symposium on Lattice Field Theory (Lattice2012), Cairns, Australia, 24-29 June 2012

We consider how to extract the pion form factors in the epsilon regime. Using the correlators with non-zero momenta and taking appropriate ratios of them, we eliminate the dominant finite volume effect from the zero-momentum pion mode. Our preliminary lattice result for the pion charge radius is consistent with the experiment., Comment: 7 pages, 3 figures, Presented at Lattice 2012, June 24-29, 2012, Cairns, Australia. To be published as PoS(Lattice2012)198

We calculate the strange quark content of the nucleon in 2+1-flavor lattice QCD. Chirally symmetric overlap fermion formulation is used to avoid the contamination from up and down quark contents due to an operator mixing between strange and light scalar operators, \bar{s}s and \bar{u}u+\bar{d}d. At a lattice spacing a=0.112(1) fm, we perform calculations at four values of degenerate up and down quark masses, which cover a range of the pion mass M_pi \simeq 300-540 MeV. We employ two different methods: one is a direct method where we calculate the strange quark content by directly inserting the strange scalar operator. The other is an indirect method where the quark content is extracted from a derivative of the nucleon mass in terms of the strange quark mass. With these two methods we obtain consistent results with each other. Our best estimate f_{T_s}=0.009(15)(16) is in good agreement with our previous studies in two-flavor QCD., Comment: 26 pages, 25 eps figures; v2: minor corrections, a reference added, version published in PRD

We report on our calculation of the kaon semileptonic form factors in Nf=2+1 lattice QCD. Chiral symmetry is exactly preserved by using the overlap quark action for a straightforward comparison with chiral perturbation theory (ChPT). We simulate three pion masses down to 290 MeV at a single lattice spacing of 0.11 fm and at a strange quark mass very close to its physical value. The form factors near zero momentum transfer are precisely calculated by using the all-to-all propagator and twisted boundary conditions. We compare the normalizations and slopes of the form factors with ChPT and experiments., Comment: 7 pages, talk presented at the XXIX International Symposium on Lattice Field Theory (Lattice 2011), Squaw Valley, Lake Tahoe, California, July 10-16, 2011

We report on a lattice simulation result for four-dimensional {\cal N}=1 SU(2) super Yang-Mills theory with the dynamical overlap gluino. We study the spectrum of the overlap Dirac operator at three different gluino masses m=0.2, 0.1 and 0.05 with the Iwasaki action on a 8^3 \times 16 lattice. We find that the lowest eigenvalue distributions are in good agreement with the prediction from the random matrix theory. Moreover the mass dependence of the condensate is almost constant, which gives a clean chiral limit. Our results for the gluino condensate in the chiral limit is < \bar{\psi} \psi > r_0^3 = 0.63(12), where r_0 is the Sommer scale., Comment: 7 pages, 3 figures. Talk presented at The XXIX International Symposium on Lattice Field Theory - Lattice 2011, July 10-16, 2011, Squaw Valley, Lake Tahoe, California

A simulation of lattice QCD at (or even below) the physical pion mass is feasible on a small lattice size of \sim 2 fm. The results are, however, subject to large finite volume effects. In order to precisely understand the chiral behavior in a finite volume, we develop a new computational scheme to interpolate the conventional epsilon and p regimes within chiral perturbation theory. In this new scheme, we calculate the two-point function in the pseudoscalar channel, which is described by a set of Bessel functions in an infra-red finite way as in the epsilon regime, while chiral logarithmic effects are kept manifest as in the p regime. The new ChPT formula is compared to our 2+1- flavor lattice QCD data near the physical up and down quark mass, mud \sim 3 MeV on an L \sim 1.8 fm lattice. We extract the pion mass = 99(4) MeV, from which we attempt a chiral "interpolation" of the observables to the physical point., Comment: 7 pages, 2 figures. Talk presented at the XXIX International Symposium on Lattice Field Theory (Lattice 2011), July 10-16, 2011, Squaw Valley, Lake Tahoe, California

We calculate the strange quark content of the nucleon directly from its disconnected three-point function in N_f=2+1 QCD. Chiral symmetry is crucial to avoid a possibly large contamination due to operator mixing, and is exactly preserved by employing the overlap quark action. We also use the all-to-all quark propagator and the low-mode averaging technique in order to accurately calculate the relevant nucleon correlator. Our preliminary result extrapolated to the physical point is f_{T_s} = 0.013(12)(16). This is in good agreement with our previous estimate in N_f=2 QCD as well as those from our indirect calculations using the Feynman-Hellmann theorem., Comment: 7pages, talk presented at The XXVIII International Symposium on Lattice Field Theory, Lattice2010

We report on our calculation of pion and kaon form factors in three-flavor QCD using the overlap quark action. Gauge ensembles are generated on a 16^3 \times 48 lattice at a lattice spacing of 0.11 fm with pion masses down to 310 MeV. Connected and disconnected meson correaltors are calculated using the all-to-all quark propagator. We present our preliminary analysis on the chiral behavior of the electromagnetic and scalar form factors as well as a comparison of the shape of the K to pi form factors with experiment., Comment: 7pages, talk presented at The XXVIII International Symposium on Lattice Field Theory, Lattice2010

Strange quark content of the nucleon is calculated in dynamical lattice QCD employing the overlap fermion formulation. For this quantity, exact chiral symmetry guaranteed by the Ginsparg-Wilson relation is crucial to avoid large contamination due to a possible operator mixing with $\bar{u}u+\bar{d}d$. Gauge configurations are generated with two dynamical flavors on a 16^3 x 32 lattice at a lattice spacing a \simeq 0.12fm. We directly calculate the relevant three-point function on the lattice including a disconnected strange quark loop utilizing the techniques of all-to-all quark propagator and low-mode averaging. Our result f_{T_s} = 0.032(8)(22), is in good agreement with our previous indirect estimate using the Feynman-Hellmann theorem., Comment: 31 pages, 22 figures; version published in PRD

We present some preliminary results of the project on finite temperature QCD with overlap fermions at KEK. We performed a series of simulations to assess the effects of fixing the topological sector at finite temperature and we will show the first calculations of topological susceptibility and meson masses for quenched and full QCD., Comment: 7 pages. 5 figures. Contribution to the proceedings of The XXVIII International Symposium on Lattice Field Theory, Lattice2010, June 14-19, 2010, Villasimius, Italy

We perform a precise calculation of the chiral condensate in QCD using lattice QCD with 2+1 flavors of dynamical overlap quarks. Up and down quark masses cover a range between 3 and 100 MeV on a 16^3x48 lattice at a lattice spacing around 0.11 fm. At the lightest sea quark mass, the finite volume system on the lattice is in the epsilon-regime. By matching the low-lying eigenvalue spectrum of the Dirac operator with the prediction of chiral perturbation theory at the next-to-leading order, we determine the chiral condensate in 2+1-flavor QCD with strange quark mass fixed at its physical value as Sigma (MS-bar at 2 GeV) = [242(04)(^+19_-18}) MeV}]^3, where the errors are statistical and systematic, respectively., Comment: 4 pages, 3 figures, errors in table 1 and fig.3 corrected. Published in PRL

We calculate the nucleon strange quark content directly from disconnected three-point functions. Numerical simulations are carried out in two-flavor QCD using the overlap quark action with up and down quark masses down to a fifth of the physical strange quark mass. To improve the statistical accuracy, we calculate the nucleon two-point functions with the low-mode averaging technique, whereas the all-to-all quark propagator is used for the disconnected quark loop. We obtain the y parameter, which is the ratio of the strange and light quark contents, y = 0.024(45) at the physical point. This is in a good agreement with our earlier calculation from the nucleon spectrum through the Feynman-Hellmann theorem., Comment: 7 pages, talk presented at The XXVII International Symposium on Lattice Field Theory (Lattice2009), July 26-31 2009, Peking University, Beijing, China

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 calculate mesonic two-point functions in the epsilon-regime of two-flavor QCD on the lattice with exact chiral symmetry. We use gauge configurations of size 16^3 32 at the lattice spacing a \sim 0.11 fm generated with dynamical overlap fermions. The sea quark mass is fixed at \sim 3 MeV and the valence quark mass is varied in the range 1-4 MeV, both of which are in the epsilon-regime. We find a good consistency with the expectations from the next-to-leading order calculation in the epsilon-expansion of (partially quenched) chiral perturbation theory. From a fit we obtain the pion decay constant F=87.3(5.6) MeV and the chiral condensate Sigma^{MS}=[239.8(4.0) MeV ]^3 up to next-to-next-to-leading order contributions., Comment: 20 pages, 12 figures, final version to appear in PRD

We calculate the meson correlators in the $\epsilon$-regime of two-flavor QCD. On a $16^3\times 32$ lattice with $a\sim 0.11$ fm, the lattice simulations are performed with the dynamical overlap fermions. We reduce the sea quark mass down to $\sim$ 3 MeV and the valence quark masses are taken in the range 1-4 MeV. The meson correlators in various channels are compared with the predictions of (partially quenched) chiral perturbation theory (ChPT). Including the NLO order of the $\epsilon$-expansion, we extract the leading-order low energy constants of ChPT, the pion decay constant $F$ and the chiral condensate $\Sigma$, as $F=87.3(5.5)$ MeV and $\Sigma^{\bar{\mathrm{MS}}}=[237.8(4.0){MeV}]^3$., Comment: 7pages, 4figures. Presented at Lattice 2007, Regensburg, Germany, 30 Jul - 4 Aug 2007, minor corrections in figures

We report on the status of the dynamical overlap QCD simulation project by the JLQCD collaboration. After completing two-flavor QCD simulation on a 16^3x32 lattice at lattice spacing a 0.12 fm, we started a series of runs with 2+1 flavors. In this report, we describe an outline of our algorithms, parameter choices, and some early physics results of this second phase of our project., Comment: 7 pages, 4 figures, presented at the XXV International Symposium on Lattice Field Theory, July 30-4 August 2007, Regensburg, Germany

We report on our calculation of the pion electromagnetic form factor with two-flavors of dynamical overlap quarks. Gauge configurations are generated using the Iwasaki gauge action on a 16^3 \times 32 lattice at the lattice spacing of 0.12fm with sea quark masses down to m_s/6, where m_s is the physical strange quark mass. We describe our setup to measure the form factor through all-to-all quark propagators and present preliminary results., Comment: 7 pages, 5 figures, talk presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4, 2007, Regensburg, Germany

We calculate the electromagnetic contribution to the pion mass difference, $\Delta m^2_\pi=m^2_{\pi^+}-m^2_{\pi^0}$, in the chiral limit through the $VV-AA$ type vacuum polarization using Das-Guralnik-Mathur-Low-Young (DGMLY) sum rule. The calculation is made with two-flavors of dynamical overlap fermions on a $16^3\times 32$ lattice at $a\sim$0.12 fm. The exact chiral symmetry of the overlap fermion is essential to control the systematic error in the difference $VV-AA$. We obtain $\Delta m_\pi^2 = 1024(100) {\rm MeV^2}$ combining the lattice data with the perturbative contribution in the high momentum region evaluated by the operator product expansion. By analyzing the momentum dependence of the vacuum polarization, we also obtain pion decay constant $f_\pi$ and the low-energy constants $L_{10}^r$ in the chiral limit., Comment: 7 pages, 2 figures, talk (Hadron spectroscopy) presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4, 2007, Regensburg, Germany

We report on a calculation of $B_K$ with two-flavor dynamical overlap fermions on a $16^3 \times 32$ lattice at $a\sim 0.12$ fm. The results are compared with the PQChPT prediction of quark mass dependence. The systematic errors due to finite volume effects and fixing topology are discussed., Comment: 7 pages, 5 figures, talk presented at Lattice 2007 (Weak decays and matrix elements)

We perform lattice simulations of two-flavor QCD using Neuberger's overlap fermion, with which the exact chiral symmetry is realized at finite lattice spacings. The epsilon-regime is reached by decreasing the light quark mass down to 3 MeV on a 16^3 32 lattice with a lattice spacing \sim 0.11 fm. We find a good agreement of the low-lying Dirac eigenvalue spectrum with the analytical predictions of the chiral random matrix theory, which reduces to the chiral perturbation theory in the epsilon-regime. The chiral condensate is extracted as \Sigma(2 GeV) = (251(7)(11) MeV)^3, where the errors are statistical and an estimate of the higher order effects in the epsilon-expansion., Comment: 10pages, 4figures

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 the two-flavor QCD simulation in the epsilon-regime using the overlap fermion formulation. Sea quark mass is reduced to ~ 2 MeV on a 16^3x32 lattice with the lattice spacing ~ 0.11fm. Topological charge is fixed at Q=0. We compare the eigenvalue distribution of the overlap-Dirac operator with the prediction of the chiral random matrix theory. Preliminary results on meson correlators are also reported., Comment: 7 pages, talk presented at Lattice 2006 (chiral symmetry); to appear in Proceedings of Science, PoS(LAT2006)050

We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the numerical costs. Our main target is a $N_f=2$ simulation with overlap fermion on a $16^3\times 32$ lattice at lattice spacing around 0.12 fm., Comment: 7 pages, 8 figures, talk given at Lattice 2006 (Algorithms, Machines, and Networks), to appear in Proceedings Of Science, POS(LAT2006)031

We launched a project to perform dymanical fermion simulations using the overlap fermion formulation for sea quarks. In order to avoid the appearace of near-zero modes of the hermitian Wilson-Dirac operator $H_W$, we introduce a pair of extra Wilson fermions with a large negative mass term. Crossing of the topological boundary is then strictly prohibited, and the topological charge is conserved during simulations. It makes the simulations substantially faster compared to the algorithms which allow the topology change. We discuss on the finite volume effects due to the fixed global topology., Comment: 7 pages, 4 figures, talk presented at Lattice 2006 (chiral symmetry); to appear in Proceedings of Science, PoS(LAT2006)052

We perform a systematic study of low-lying eigenmodes of $H_w$ with various gauge actions to find the optimal choice for dynamical overlap fermion simulations, with which one may achieve lower numerical cost for HMC and better locality property of the overlap kernel. For this purpose, our study is made with emphasis on the distribution of low-lying eigenvalues and the mobility edge with and without dynamical overlap fermions., Comment: 7 pages, 7 figures, talk presented at Lattice 2006 (chiral symmetry). References added

We propose a lattice action including unphysical Wilson fermions with a negative mass m_0 of the order of the inverse lattice spacing. With this action, the exact zero mode of the hermitian Wilson-Dirac operator H_W(m_0) cannot appear and near-zero modes are strongly suppressed. By measuring the spectral density rho(lambda_W), we find a gap near lambda_W=0 on the configurations generated with the standard and improved gauge actions. This gap provides a necessary condition for the proof of the exponential locality of the overlap-Dirac operator by Hernandez, Jansen, and Luescher. Since the number of near-zero modes is small, the numerical cost to calculate the matrix sign function of H_W(m_0) is significantly reduced, and the simulation including dynamical overlap fermions becomes feasible. We also introduce a pair of twisted mass pseudo-fermions to cancel the unwanted higher mode effects of the Wilson fermions. The gauge coupling renormalization due to the additional fields is then minimized. The topological charge measured through the index of the overlap-Dirac operator is conserved during continuous evolutions of gauge field variables., Comment: 22 pages, 6 figures; a section added, the version to appear in Phys. Rev. D

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)

Three-flavor QCD simulation with the $O(a)$-improved Wilson fermion action is made employing an exact fermion algorithm developed for odd number of quark flavors. For the plaquette gauge action, an unexpected first-order phase transition is found in the strong coupling regime ($\beta\lesssim$ 5.0) at relatively heavy quark masses ($m_{\mathrm{PS}}/m_{\mathrm{V}}\sim$ 0.74--0.87). Strong metastability persists on a large lattice of size $12^3\times 32$, which indicates that the transition has a bulk nature. The phase gap becomes smaller toward weaker couplings and vanishes at $\beta\simeq 5.0$, which corresponds to a lattice spacing $a\simeq$ 0.1 fm. The phase transition is not found if the improved gauge actions are employed. Our results imply that realistic simulations of QCD with three flavors of dynamical Wilson-type fermions at lattice spacings in the range $a=$ 0.1--0.2 fm require use of improved gauge actions. Possible origins of the phase transition is discussed., Comment: REVTeX4, 9 pages, 11 figures

We present an unquenched lattice calculation for the $B^0-\bar{B}^0$ transition amplitude. The calculation, carried out at an inverse lattice spacing $1/a$ = 2.22(4) GeV, incorporates two flavors of dynamical quarks described by the fully $O(a)$-improved Wilson fermion action and heavy quarks described by NRQCD. A particular attention is paid to the uncertainty that arises from the chiral extrapolation, especially the effect of pion loops, for light quarks, which we find could be sizable for the leptonic decay of the $B_d$ meson, whereas it is small for the $B_s$ meson and the $B$ parameters. We obtain $f_{B_d}=191(10)(^{+12}_{-22})$ MeV, $f_{B_s}/f_{B_d}=1.13(3)(^{+13}_{-\ 2})$, $B_{B_d}(m_b)=0.836(27)(^{+56}_{-62})$, $B_{B_s}/B_{B_d}=1.017(16)(^{+56}_{-17})$, $\xi=1.14(3)(^{+13}_{- 2})$, where the first error is statistical, and the second is systematic, including uncertainties due to chiral extrapolation, finite lattice spacing, heavy quark expansion and perturbative operator matching., Comment: 4 pages, 3 figures; version to appear in PRL

We present a lattice QCD calculation of the heavy quark expansion parameters $\mu_{\pi}^2$ and $\mu_G^2$ for heavy-light mesons and heavy-light-light baryons. The calculation is carried out on a 20$^3\times$48 lattice at $\beta$ = 6.0 in the quenched approximation, using the lattice NRQCD action for heavy quarks. We obtain the parameters $\mu_{\pi}^2$ and $\mu_G^2$ in two different methods: a direct calculation of the matrix elements and an indirect calculation through the mass spectrum, and confirm that the both methods give consistent results. We also discuss an application to the lifetime ratios., Comment: REVTeX4, 45 pages, 37 eps figures

We present a high statistics study of the light hadron spectrum and quark masses in QCD with two flavors of dynamical quarks. Numerical simulations are carried out using the plaquette gauge action and the O(a)-improved Wilson quark action at \beta=5.2, where the lattice spacing is found to be a=0.0887(11)fm from \rho meson mass, on a 20^3\times 48 lattice. At each of five sea quark masses corresponding to m_{PS}/m_{V} \simeq 0.8-0.6, we generate 12000 trajectories using the symmetrically preconditioned Hybrid Monte Carlo algorithm. Finite spatial volume effects are investigated employing 12^3 \times 48, 16^3 \times 48 lattices. We also perform a set of simulations in quenched QCD with the same lattice actions at a similar lattice spacing to those for the full QCD runs. In the meson sector we find clear evidence of sea quark effects. The J parameter increases for lighter sea quark masses, and the full QCD meson masses are systematically closer to experiment than in quenched QCD. Careful finite-size studies are made to ascertain that these are not due to finite-size effects. Evidence of sea quark effects is less clear in the baryon sector due to larger finite-size effects. We also calculate light quark masses and find m_{ud}^{MS}(2GeV) =3.223(+0.046/-0.069)MeV and m_s^{MS}(2GeV)=84.5(+12.0/-1.7)MeV which are about 20% smaller than in quenched QCD., Comment: 52 pages, LaTeX2e with RevTeX4; a new figure (Fig.5) added, Fig.20 and typos corrected, the final version to appear in PRD

We test the one-loop chiral perturbation theory formula on unquenched lattice data of pseudoscalar meson decay constants. The chiral extrapolation including the effect of the chiral logarithm is attempted and its uncertainty is discussed., Comment: 3 pages, 4 figures, talk presented at Lattice2002(spectrum)

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., Comment: 3 pages, 2 figures, Lattice2002(algor)

We present the final results of the JLQCD calculation of the light hadron spectrum and quark masses with two flavors of dynamical quarks using the plaquette gauge action and fully $O(a)$-improved Wilson quark action at $\beta=5.2$. We observe that sea quark effects lead to a closer agreement of the strange meson and baryon masses with experiment and a reduction of quark masses by about 25., Comment: 3 pages, 5 figures, Lattice2002 (spectrum)