BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

N. Yamada

N. Yamada appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phgr-qcastro-ph.IMhep-exmath.AP

Catalog footprint

What is connected

58works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

58 published item(s)

preprint2016arXiv

Topological susceptibility at high temperature on the lattice

QCD topological susceptibility at high temperature, $χ_t(T)$, provides an important input for the estimate of the axion abundance in the present Universe. While the model independent determination of $χ_t(T)$ should be possible from the first principles using lattice QCD, existing methods fail at high temperature, since not only the probability that non-trivial topological sectors appear in the configuration generation process but also the local topological fluctuations get strongly suppressed. We propose a novel method to calculate the temperature dependence of topological susceptibility at high temperature. A feasibility test is performed on a small lattice in the quenched approximation, and the results are compared with the prediction of the dilute instanton gas approximation. It is found that the method works well especially at very high temperature and the result is consistent with the instanton calculus down to $T\sim 2\, T_c$ within the statistical uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Lagrangian Approach to Weakly Coupled Hamilton-Jacobi Systems

We study a class of weakly coupled Hamilton-Jacobi systems with a specific aim to perform a qualitative analysis in the spirit of weak KAM theory. Our main achievement is the definition of a family of related action functionals containing the Lagrangians obtained by duality from the Hamiltonians of the system. We use them to characterize, by means of a suitable estimate, all the subsolutions of the system, and to explicitly represent some subsolutions enjoying an additional maximality property. A crucial step for our analysis is to put the problem in a suitable random frame. Only some basic knowledge of measure theory is required, and the presentation is accessible to readers without background in probability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong coupling constant from vacuum polarization functions in three-flavor lattice QCD with dynamical overlap fermions

We determine the strong coupling constant $α_s$ from a lattice calculation of vacuum polarization functions (VPF) in three-flavor QCD with dynamical overlap fermions. Fitting lattice data of VPF to the continuum perturbative formula including the operator product expansion, we extract the QCD scale parameter $Λ_{\overline{MS}}^{(3)}$. At the $Z$ boson mass scale, we obtain $α_s^{(5)}(M_Z)=0.1181(3)(^{+14}_{-12})$, where the first error is statistical and the second is our estimate of various systematic uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lattice Study on quantum-mechanical dynamics of two-color QCD with six light flavors

We investigate the chiral properties of SU(2) gauge theory with six flavors, i.e. six light Dirac fermions in the fundamental representations by lattice simulation, and point out that the spontaneous breakdown of chiral symmetry does not occur in this system. The quark mass dependence of the mesonic spectrum provides an evidence for such a possibility. The decay constant tends to be increased by the finite size effect, which is opposite to the behavior predicted by chiral perturbation theory and indicates that the long distance dynamics in the six-flavor theory could be different from the theory with chiral symmetry breaking. The subtracted chiral condensate, whose utility is demonstrated by the simulation of two-flavor theory, is shown to vanish in the chiral limit within the precision of available data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Running coupling constant and mass anomalous dimension of six-flavor SU(2) gauge theory

In the exploration of viable models of dynamical electroweak symmetry breaking, it is essential to locate the lower end of the conformal window and know the mass anomalous dimensions there for a variety of gauge theories. We calculate, with the Schrödinger functional scheme, the running coupling constant and the mass anomalous dimension of SU(2) gauge theory with six massless Dirac fermions in the fundamental representation. The calculations are performed on $6^4$ - $24^4$ lattices over a wide range of lattice bare couplings to take the continuum limit. The discretization errors for both quantities are removed perturbatively. We find that the running slows down and comes to a stop at $0.06 \lesssim 1/g^2 \lesssim 0.15$ where the mass anomalous dimension is estimated to be $0.26 \lesssim γ^*_m \lesssim 0.74$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice study on two-color QCD with six flavors of dynamical quarks

We study the dynamics of SU(2) gauge theory with NF=6 Dirac fermions by means of lattice simulation to investigate if they are appropriate to realization of electroweak symmetry breaking. The discrete analogue of beta function for the running coupling constant defined under the Schroedinger functional boundary condition are computed on the lattices up to linear size of L/a=24 and preclude the existence of infrared fixed point below 7.6. Gluonic observables such as heavy quark potential, string tension, Polyakov loop suggest that the target system is in the confining phase even in the massless quark limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chiral interpolation in a finite volume

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Determination of the chiral condensate from QCD Dirac spectrum on the lattice

We calculate the chiral condensate of QCD with 2, 2+1 and 3 flavors of sea quarks. Lattice QCD simulations are performed employing dynamical overlap fermions with up and down quark masses covering a range between 3 and 100 MeV. On L ~ 1.8-1.9 fm lattices at a lattice spacing ~ 0.11 fm, we calculate the eigenvalue spectrum of the overlap-Dirac operator. By matching the lattice data with the analytical prediction from chiral perturbation theory at the next-to-leading order, the chiral condensate in the massless limit of up and down quarks is determined.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-precision CTE measurement of hybrid C/SiC composite for cryogenic space telescopes

This paper presents highly precise measurements of thermal expansion of a "hybrid" carbon-fiber reinforced silicon carbide composite, HB-Cesic\textregistered - a trademark of ECM, in the temperature region of \sim310-10K. Whilst C/SiC composites have been considered to be promising for the mirrors and other structures of space-borne cryogenic telescopes, the anisotropic thermal expansion has been a potential disadvantage of this material. HB-Cesic\textregistered is a newly developed composite using a mixture of different types of chopped, short carbon-fiber, in which one of the important aims of the development was to reduce the anisotropy. The measurements indicate that the anisotropy was much reduced down to 4% as a result of hybridization. The thermal expansion data obtained are presented as functions of temperature using eighth-order polynomials separately for the horizontal (XY-) and vertical (Z-) directions of the fabrication process. The average CTEs and their dispersion (1σ) in the range 293-10K derived from the data for the XY- and Z-directions were 0.805$\pm$0.003\times10$^{-6}$ K$^{-1}$ and 0.837\pm0.001\times10$^{-6}$ K$^{-1}$, respectively. The absolute accuracy and the reproducibility of the present measurements are suggested to be better than 0.01\times10$^{-6}$ K$^{-1}$ and 0.001\times(10)^{-6} K^{-1}, respectively. The residual anisotropy of the thermal expansion was consistent with our previous speculation regarding carbon-fiber, in which the residual anisotropy tended to lie mainly in the horizontal plane.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Running coupling constant of ten-flavor QCD with the Schrödinger functional method

Walking technicolor theory attempts to realize electroweak symmetry breaking as the spontaneous chiral symmetry breakdown caused by the gauge dynamics with slowly varying gauge coupling constant and large mass anomalous dimension. Many-flavor QCD is one of the candidates owning these features. We focus on the SU(3) gauge theory with ten flavors of massless fermions in the fundamental representation, and compute the gauge coupling constant in the Schrödinger functional scheme. Numerical simulation is performed with $O(a)$-unimproved lattice action, and the continuum limit is taken in linear in lattice spacing. We observe evidence that this theory possesses an infrared fixed point.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two-photon decay of $π^0$ from two-flavor lattice QCD

We study the correction to the radiative $π^0$ decay width due to finite light quark mass. Using lattice QCD with the overlap fermion formulation, we calculate the three-point function of the form $<PV_μV_ν>$ in the (Euclidean) momentum space, which corresponds to the $π^0\rightarrow γ^*γ^*$ amplitude. To fit the lattice data, we use two different modifications of vector meson dominance (VMD) ansatz. One is a combined form of VMD with the next-to-leading order (NLO) chiral perturbation theory (ChPT), and the other is a resummed form of pion-loop diagrams. We extract one of the low energy constants in NLO ChPT, and estimate $π^0\rightarrow γγ$ decay width including finite mass correction.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Determination of the chiral condensate from 2+1-flavor lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electromagnetic mass splittings of the low lying hadrons and quark masses from 2+1 flavor lattice QCD+QED

Results are presented for the electromagnetic mass splittings of the low lying hadrons. These are used to determine the non-degenerate light quark masses. It is found that m_u=2.24(10)(34), m_d=4.65(15)(32), and $m_s=97.6(2.9)(5.5)$ MeV (MSbar scheme, 2 GeV scale). The first error is statistical and the second systematic. We find the lowest order electromagnetic splitting (m_pi+-m_pi0)_QED=3.38(23) MeV, the splittings including next-to-leading order, (m_pi+-m_pi0)_QED=4.50(23) MeV, (m_K+-m_K0)_QED=1.87(10) MeV, and the m_u != m_d contribution to the kaon mass difference, (m_K+-m_K0)_(m_u-m_d)=-5.840(96) MeV. All errors are statistical only, and the next-to-leading order pion splitting is only approximate; it does not contain all next-to-leading order contributions. We also computed the proton-neutron mass difference, including for the first time, QED interactions in a realistic 2+1 flavor calculation. We find $(m_p-m_n)_{\rm QED}=0.383(68)$ MeV, (m_p-m_n)_(m_u-m_d)=-2.51(14) MeV, and the total m_p-m_n=-2.13(16)(70) MeV, where the first error is statistical, and the second, part of the systematic error. We use domain wall fermions and the Iwasaki gauge action (gauge coupling beta=2.13). We use two lattice sizes, 16^3 and 24^3, to address finite volume effects. Non-compact QED is treated in the quenched approximation. We present new results for the electromagnetic low energy constants in SU(3) and SU(2) partially-quenched chiral perturbation theory to the next-to-leading order, obtained from fits to our data. Detailed analysis of systematic errors in our results and methods for improving them are discussed. Finally, new analytic results for SU(2)_L x SU(2)_R-plus-kaon chiral perturbation theory, including the one-loop logs proportional to alpha_em*m, are given.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Improving many flavor QCD simulations using multiple GPUs

We accelerate many-flavor lattice QCD simulations using multiple GPUs. Multiple pseudo-fermion fields are introduced additively and independently for each flavor in the many-flavor HMC algorithm. Using the independence of each pseudo-fermion field and the blocking technique for the quark solver, we can assign the solver task to each GPU card. In this report we present the blocking technique for the many-flavor dynamical QCD simulations. We investigate the effect of the blocking and the acceleration with the multiple GPUs for the Schrödinger functional simulations with Wilson SU(3) plaquette gauge action and $N_f=10$ Wilson fermions. Five pseudo-fermion fields are introduced and the quark solver task is distributed in the ratio of 2:3 to two GPUs. We expect a 40% timing reduction from the single GPU case and have observed a 34% timing reduction in the test simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nucleon strange quark content in 2+1-flavor QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Study of the running coupling constant in 10-flavor QCD with the Schrödinger functional method

The electroweak gauge symmetry is allowed to be spontaneously broken by the strongly interacting vector-like gauge dynamics. When the gauge coupling of a theory runs slowly in a wide range of energy scale, the theory is a candidate for walking technicolor. This may open up the possibility that the origin of all masses may be traced back to the gauge theory. We use the \SF method to see whether the gauge coupling of 10-flavor QCD "walks" or not. Preliminary result is reported.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

$π^0$ to two-photon decay in lattice QCD

We calculate the neutral pion ($π^0$) to off-shell two photon ($γ^*γ^*$) transition form factor in lattice QCD. The transition form factor can be extracted from the three-point function of the form (axial-vector)-(vector)-(vector) as a function of off-shell two-photon momentum. Since the axial-anomaly plays an important role in the $π^0\toγγ$ decay process, we employ the overlap fermion, which preserves the exact chiral symmetery on the lattice. After extrapolating to the chiral and the vanishing photon momentum limit with a fit function based on vector meson dominance (VMD) model, we find that the Adler-Bell-Jackiw anomaly is correctly reproduced.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Calculation of nucleon strange quark content with dynamical overlap quarks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Chiral properties of light mesons with $N_f=2+1$ overlap fermions

We present an update of the light meson spectrum with $N_f$=2+1 overlap fermions on a $16^3\times 48$ lattice at five different up and down quark masses and two strange quark masses. Based on our experience with the previous simulation with $N_f=2$, we carry out the chiral extrapolation with the prediction of the chiral perturbation theory at the next-to-next-to leading order. We also check the consistency of our analysis by using alternative chiral extrapolation with a reduced theory in which the strange quark mass is integrated out.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Flavor-singlet mesons in N_f=2+1 QCD with dynamical overlap quarks

We report on our study of flavor-singlet mesons in three-flavor QCD with dynamical overlap quarks. Gauge ensembles are generated on a 16^3 x 48 lattice at a lattice spacing of 0.10 fm with the strange quark masses around its physical value m_s,phys and up and down quark masses down to m_s,phys/5. Connected and disconnected meson correlators are calculated using the all-to-all quark propagator. We present our preliminary results on the spectrum of flavor-singlet pseudoscalar and vector mesons.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-perturbative renormalization of bilinear operators with dynamical overlap fermions

Using the non-perturbative renormalization technique, we calculate the renormalization factors for quark bilinear operators made of overlap fermions on the lattice. The background gauge field is generated by the JLQCD and TWQCD collaborations including dynamical effects of two or 2+1 flavors of light quarks on a 16$^3\times$32 or 16$^3\times$48 lattice at lattice spacing around 0.1 fm. By reducing the quark mass close to the chiral limit, where the finite volume system enters the so-called $ε$-regime, the unwanted effect of spontaneous chiral symmetry breaking on the renormalization factors is suppressed. On the lattices in the conventional $p$-regime, this effect is precisely subtracted by separately calculating the contributions from the chiral condensate.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nucleon sigma term and strange quark content in 2+1-flavor QCD with dynamical overlap fermions

We study the sigma term and the strange quark content of nucleon in 2+1-flavor QCD with dynamical overlap fermions. We analyze the lattice data of nucleon mass taken at two different strange quark masses with five values of up and down quark masses each. Using the reweighting technique, we study the strange quark mass dependence of the nucleon and extract the strange quark content.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Pion form factors from two-flavor lattice QCD with exact chiral symmetry

We calculate pion vector and scalar form factors in two-flavor lattice QCD and study the chiral behavior of the vector and scalar radii <r^2>_{V,S}. Numerical simulations are carried out on a 16^3 x 32 lattice at a lattice spacing of 0.12 fm with quark masses down to \sim m_s/6, where m_s is the physical strange quark mass. Chiral symmetry, which is essential for a direct comparison with chiral perturbation theory (ChPT), is exactly preserved in our calculation at finite lattice spacing by employing the overlap quark action. We utilize the so-called all-to-all quark propagator in order to calculate the scalar form factor including the contributions of disconnected diagrams and to improve statistical accuracy of the form factors. A detailed comparison with ChPT reveals that the next-to-next-to-leading-order contributions to the radii are essential to describe their chiral behavior in the region of quark mass from m_s/6 to m_s/2. Chiral extrapolation based on two-loop ChPT yields <r^2>_V=0.409(23)(37)fm and <r^2>_S=0.617(79)(66)fm, which are consistent with phenomenological analysis. We also present our estimates of relevant low-energy constants.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Study of the running coupling constant in 10-flavor QCD with the Schroedinger functional method

The electroweak gauge symmetry is allowed to be spontaneously broken by the strongly interacting vector-like gauge dynamics. When the gauge coupling of a theory runs slowly in a wide range of energy scale, the theory is extremely interesting. This may open up the possibility that the origin of all masses may be traced back to the gauge theory. We use the SF method to determine the scale dependence of the gauge coupling of 10-flavor QCD. Preliminary results are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Calculation of the nucleon sigma term and strange quark content with two flavors of dynamical overlap fermions

We present a calculation of the nucleon sigma term on two-flavor QCD configurations with dynamical overlap fermions. We analyse the lattice data for the nucleon mass using the baryon chiral perturbation theory. Using partially quenched data sets, we extract the connected and disconnected contributions to the nucleon sigma term separately. Chiral symmetry on the lattice simplifies the determination of the disconnected contribution. We find that the strange quark content, which determines the neutralino dark matter reaction rate with nucleon through the Higgs boson exchange, is much smaller than the previous lattice results.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Light meson spectrum with $N_f=2+1$ dynamical overlap fermions

We report on a numerical simulation with 2+1 dynamical flavors of overlap fermions. We calculate pseudo-scalar masses and decay constants on a $16^3\times 48 \times (0.11 {\rm fm})^4$ lattice at five different up and down quark masses and two strange quark masses. The lightest pion mass corresponds to $\approx 310$ MeV. We also study the validity of the chiral perturbation theory using the results of the numerical simulation with two dynamical flavors and conclude that the one-loop formulae cannot be directly applied in the strange quark mass region. We therefore extrapolate our 2+1-flavor results to the chiral limit by fitting the data to the two-loop formulae of the chiral perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Light quark masses from unquenched lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Pion vector and scalar form factors with dynamical overlap quarks

We calculate the pion vector and scalar form factors in two-flavor QCD. Gauge configurations are generated with dynamical overlap quarks on a 16^3 x 32 lattice at a lattice spacing of 0.12 fm with sea quark masses down to a sixth of the physical strange quark mass. Contributions of disconnected diagrams to the scalar form factor is calculated employing the all-to-all quark propagators. We present a detailed comparison of the vector and scalar radii with chiral perturbation theory to two loops.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Topological susceptibility in (2+1)-flavor lattice QCD with overlap fermion

We determine the topological susceptibility χ_t in the topologically-trivial sector generated by lattice simulations of N_f = 2+1 QCD with overlap Dirac fermion, on a 16^3 x 48 lattice with lattice spacing ~ 0.11 fm, for five sea quark masses m_q ranging from m_s/6 to m_s (where m_s is the physical strange quark mass). The χ_t is extracted from the plateau (at large time separation) of the 2-point and 4-point time-correlation functions of the flavor-singlet pseudoscalar meson η', which arises from the finite size effect due to fixed topology. In the small m_q regime, our result of χ_t agrees with the chiral effective theory. Using the formula χ_t = Σ(m_u^{-1} + m_d^{-1} + m_s^{-1})^{-1} by Leutwyler-Smilga, we obtain the chiral condensate Σ^{MSbar}(2 GeV) = [249(4)(2) MeV]^3.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

2+1 flavor light hadron spectrum and quark masses with the O(a) improved Wilson-clover quark formalism

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}}(μ=2 GeV)=3.50(14)({}^{+26}_{-15}) MeV and m_s^{\bar{MS}}(μ=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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Two-flavor lattice QCD in the epsilon-regime and chiral Random Matrix Theory

The low-lying eigenvalue spectrum of the QCD Dirac operator in the epsilon-regime is expected to match with that of chiral Random Matrix Theory (ChRMT). We study this correspondence for the case including sea quarks by performing two-flavor QCD simulations on the lattice. Using the overlap fermion formulation, which preserves exact chiral symmetry at finite lattice spacings, we push the sea quark mass down to \sim 3 MeV on a 16^3\times 32 lattice at a lattice spacing a \simeq 0.11 fm. We compare the low-lying eigenvalue distributions and find a good agreement with the analytical predictions of ChRMT. By matching the lowest-lying eigenvalue we extract the chiral condensate, Σ(2 GeV)[MSbar] = [251(7)(11) MeV]^3, where errors represent statistical and higher order effects in the epsilon expansion. We also calculate the eigenvalue distributions on the lattices with heavier sea quarks at two lattice spacings. Although the epsilon expansion is not applied for those sea quarks, we find a reasonable agreement of the Dirac operator spectrum with ChRMT. The value of Sigma, after extrapolating to the chiral limit, is consistent with the estimate in the epsilon-regime.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

An estimate of the eta and eta-prime meson masses in Nf=2+1 lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Nonperturbative $O(a)$ improvement of the Wilson quark action with the RG-improved gauge action using the Schrödinger functional method

We perform a nonperturbative determination of the $O(a)$-improvement coefficient $c_{\rm SW}$ and the critical hopping parameter $κ_c$ for $N_f$=3, 2, 0 flavor QCD with the RG-improved gauge action using the Schrödinger functional method. In order to interpolate $c_{\rm SW}$ and $κ_c$ as a function of the bare coupling, a wide range of $β$ from the weak coupling region to the moderately strong coupling points used in large-scale simulations is studied. Corrections at finite lattice size of $O(a/L)$ turned out to be large for the RG-improved gauge action, and hence we make the determination at a size fixed in physical units using a modified improvement condition. This enables us to avoid $O(a)$ scaling violations which would remain in physical observables if $c_{\rm SW}$ determined for a fixed lattice size $L/a$ is used in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Kaon B-parameter from Quenched Domain-Wall QCD

We present numerical results for the kaon B-parameter, B_K, determined in the quenched approximation of lattice QCD. Our simulations are performed using domain-wall fermions and the renormalization group improved, DBW2 gauge action which combine to give quarks with good chiral symmetry at finite lattice spacing. Operators are renormalized non-perturbatively using the RI/MOM scheme. We study scaling by performing the simulation on two different lattices with a^{-1} = 1.982(30) and 2.914(54) GeV. We combine this quenched scaling study with an earlier calculation of B_K using two flavors of dynamical, domain-wall quarks at a single lattice spacing to obtain B_K(MS,NDR,mu=2GeV)=0.563(21)(39)(30), were the first error is statistical, the second systematic (without quenching errors) and the third estimates the error due to quenching.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A canonical formalism of f(T)-type gravity in terms of Lie derivatives

A canonical formalism of f(R)-type gravity is proposed, resolving the problem in the formalism of Buchbinder and Lyakhovich(BL). The new coordinates corresponding to the time derivatives of the metric are taken to be its Lie derivatives which is the same as in BL. The momenta canonically conjugate to them and Hamiltonian density are defined similarly to the formalism of Ostrogradski. It is shown that our method surely resolves the problem of BL.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Non-perturbative O(a)-improvement of Wilson quark action in three-flavor QCD with plaquette gauge action

We perform a non-perturbative determination of the O(a)-improvement coefficient c_SW for the Wilson quark action in three-flavor QCD with the plaquette gauge action. Numerical simulations are carried out in a range of β=12.0-5.2 on a single lattice size of 8^3x16 employing the Schrödinger functional setup of lattice QCD. As our main result, we obtain an interpolation formula for c_SW and the critical hopping parameter K_c as a function of the bare coupling. This enables us to remove O(a) scaling violation from physical observables in future numerical simulation in the wide range of β. Our analysis with a perturbatively modified improvement condition for c_SW suggests that finite volume effects in c_SW are not large on the 8^3x16 lattice. We investigate N_f dependence of c_SW by additional simulations for N_f=4, 2 and 0 at β=9.6. As a preparatory step for this study, we also determine c_SW in two-flavor QCD at β=5.2. At this β, several groups carried out large-scale calculations of the hadron spectrum, while no systematic determination of c_SW has been performed.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Bulk first-order phase transition in three-flavor lattice QCD with $O(a)$-improved Wilson fermion action at zero temperature

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 ($β\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 $β\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.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

An Exact Algorithm for Any-flavor Lattice QCD with Kogut-Susskind Fermion

We propose an exact simulation algorithm for lattice QCD with dynamical Kogut-Susskind fermion in which the N_f-flavor fermion operator is defined as the N_f/4-th root of the Kogut-Susskind (KS) fermion operator. The algorithm is an extension of the Polynomial Hybrid Monte Carlo (PHMC) algorithm to KS fermions. The fractional power of the KS fermion operator is approximated with a Hermitian Chebyshev polynomial, with which we can construct an algorithm for any number of flavors. The error which arises from the approximation is corrected by the Kennedy-Kuti noisy Metropolis test. Numerical simulations are performed for the two-flavor case for several lattice parameters in order to confirm the validity and the practical feasibility of the algorithm. In particular tests on a 16^4 lattice with a quark mass corresponding to m_{PS}/m_V ~ 0.68 are successfully accomplished. We conclude that our algorithm provides an attractive exact method for dynamical QCD simulations with KS fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Heavy quark expansion parameters from lattice NRQCD

We present a lattice QCD calculation of the heavy quark expansion parameters $μ_π^2$ and $μ_G^2$ for heavy-light mesons and heavy-light-light baryons. The calculation is carried out on a 20$^3\times$48 lattice at $β$ = 6.0 in the quenched approximation, using the lattice NRQCD action for heavy quarks. We obtain the parameters $μ_π^2$ and $μ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Study of finite volume effects in the non-perturbative determination of $\csw$ with the SF method in full three-flavor lattice QCD

The non-perturbative $\csw$ determined by the Schrödinger functional (SF) method with the RG-improved gauge action in dynamical $N_f=3$ QCD shows a finite volume effect when the numerical simulations are carried out at a constant lattice size $L/a$. We remove the unwanted finite volume effect by keeping physical lattice extent $L$ at a constant. The details of the method and the result obtained for non-perturbative $\csw$ with a constant $L$ are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

An exact Polynomial Hybrid Monte Carlo algorithm for dynamical Kogut-Susskind fermions

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 $β$$=$5.7, $m$$=$0.02 ($m_{\mathrm{PS}}/m_{\mathrm{V}}$$\sim$0.69) with a reasonable computational time.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

I=2 Pion Scattering Length with the Wilson Fermion

The calculation of the I=2 pion scattering length in quenched lattice QCD is revisited. The calculation is carried out with the Wilson fermion action employing Lüscher's finite size scaling method at $β=5.9$, 6.1, and 6.3 corresponding to the range of lattice spacing $a\simeq 0.12 - 0.07$ fm. We obtain in the continuum limit $a_0m_π= -0.0410(69)$, which is consistent with the prediction of chiral perturbation theory $a_0m_π=-0.0444(10)$.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Polynomial Hybrid Monte Carlo algorithm for lattice QCD with an odd number of flavors

We present a polynomial hybrid Monte Carlo (PHMC) algorithm for lattice QCD with odd numbers of flavors of O(a)-improved Wilson quark action. The algorithm makes use of the non-Hermitian Chebyshev polynomial to approximate the inverse square root of the fermion matrix required for an odd number of flavors. The systematic error from the polynomial approximation is removed by a noisy Metropolis test for which a new method is developed. Investigating the property of our PHMC algorithm in the N_f=2 QCD case, we find that it is as efficient as the conventional HMC algorithm for a moderately large lattice size (16^3 times 48) with intermediate quark masses (m_{PS}/m_V ~ 0.7-0.8). We test our odd-flavor algorithm through extensive simulations of two-flavor QCD treated as an N_f = 1+1 system, and comparing the results with those of the established algorithms for N_f=2 QCD. These tests establish that our PHMC algorithm works on a moderately large lattice size with intermediate quark masses (16^3 times 48, m_{PS}/m_V ~ 0.7-0.8). Finally we experiment with the (2+1)-flavor QCD simulation on small lattices (4^3 times 8 and 8^3 times 16), and confirm the agreement of our results with those obtained with the R algorithm and extrapolated to a zero molecular dynamics step size.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

One-loop renormalization of heavy-light currents

We calculate the mass dependent renormalization factors of heavy-light bilinears at one-loop order of perturbation theory, when the heavy quark is treated with the Fermilab formalism. We present numerical results for the Wilson and Sheikholeslami-Wohlert actions, with and without tree-level rotation. We find that in both cases our results smoothly interpolate from the static limit to the massless limit. We also calculate the mass dependent Brodsky-Lepage-Mackenzie scale $q^*$, with and without tadpole-improvement.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

$B \to πl\barν$ Form Factors with NRQCD Heavy Quark and Clover Light Quark Actions

We report results on semileptonic $B\toπl\barν$ decay form factors near $q^2_{\rm max}$ using NRQCD heavy quark and clover light quark actions and currents improved through $O(αa)$. An inconsistency with the soft pion relation $f^0(q^2_{\rm max})= f_B/f_π$ found in a previous work is confirmed, and a possible solution with nonperturbative renormalization is discussed. We find that $f^+(q^2)$ is well described by the $B^*$ pole near $q^2_{\rm max}$, and its $1/M_B$ scaling is also consistent with the prediction of the pole dominance model.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

$O(α_{s}a)$ matching coefficients for axial vector current and $ΔB$$=$2 operator

We present a calculation of the perturbative matching coefficients including mixing with higher dimensional operators for the temporal component of the heavy-light axial current, $A_{4}$, and the $ΔB=2$ operator, $O_S$. For $A_{4}^{\scriptsize static, NRQCD}$, calculations with various RG-improved gauge actions are peformed. Matching coefficients with NRQCD and heavy-clover actions are also compared.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Non-perturbative determination of quark masses in quenched lattice QCD with the Kogut-Susskind fermion action

We report results of quark masses in quenched lattice QCD with the Kogut-Susskind fermion action, employing the Reguralization Independent scheme (RI) of Martinelli et al. to non-perturbatively evaluate the renormalization factor relating the bare quark mass on the lattice to that in the continuum. Calculations are carried out at β=6.0, 6.2, and 6.4, from which we find $m^{\bar{MS}}_{ud} (2 GeV)= 4.23(29) MeV$ for the average up and down quark mass and, with the $ϕ$ meson mass as input, $m^{\bar{MS}}_{s} (2 GeV)= 129(12) MeV$ for the strange mass in the continuum limit. These values are about 20% larger than those obtained with the one-loop perturbative renormalization factor.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Nucleon decay matrix elements with the Wilson quark action: an update

We present preliminary results of a new lattice computation of hadronic matrix elements of baryon number violating operators which appear in the low-energy effective Lagrangian of (SUSY-)Grand Unified Theories. The contribution of irrelevant form factor which has caused an underestimate of the matrix elements in previous studies is subtracted in this calculation. Our results are 2$\sim$4 times larger than the most conservative values often employed in phenomenological analyses of nucleon decay with specific GUT models.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Pion decay constant in quenched QCD with Kogut-Susskind quarks

We present a non-perturbative calculation for the pion decay constant with quenched Kogut-Susskind quarks. Numerical simulations are carried out at $β= 6.0$ and 6.2 with various operators extending over all flavors. The renormalization correction is applied for each flavor by computing non-perturbative renormalization constants, and it is compared with a perturbative calculation. We also study the behavior of $f_π$ in the continuum limits for both non-perturbative and perturbative calculations. The results in the continuum limit is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Non-perturbative renormalization factors of bilinear quark operators for Kogut-Susskind fermions and light quark masses in quenched QCD

Light quark masses are computed for Kogut-Susskind fermions by evaluating non-perturbatively the renormalization factor for bilinear quark operators. Calculations are carried out in the quenched approximation at β=6.0, 6.2, and 6.4. For the average up and down quark mass we find $m_{\bar MS}(2 GeV)= 4.15(27) MeV$ in the continuum limit, which is significantly larger than $3.51(20) MeV$ ($q^*=1/a$) or $3.40(21) MeV$ ($q^*=π/a$) obtained with the one-loop perturbative renormalization factor.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A study of $O(1/m_{Q}^{2})$ corrections for $f_{B}$ with lattice NRQCD

We investigate higher order effects in the nonrelativistic expansion of lattice QCD on the heavy-light meson decay constants and some other quantities in order to understand the truncation error of NRQCD. While our numerical results have large $O(a)$ and $O(α_s)$ errors due to the use of Wilson light quark action and the tree-level matching, we find that the truncation error of higher order relativistic corrections are adequately small around the mass of the $b$ quark. We also present a perturbatively matched results through 1-loop level without operator mixing effects.

0 citations0 reviewsNo code linkedOpen access