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

K. -I. Ishikawa

K. -I. Ishikawa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phhep-thnucl-thphysics.comp-ph

Catalog footprint

What is connected

36works
5topics
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

36 published item(s)

preprint2016arXiv

Non-trivial Center Dominance in High Temperature QCD

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 $β=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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2+1 flavor QCD simulation on a $96^4$ lattice

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 $β=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_π$,$m_K$ and $m_Ω$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Calculation of $K \to ππ$ decay amplitudes with improved Wilson fermion action in lattice QCD

We present our result for the $K\toππ$ decay amplitudes for both the $Δ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_π$. 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_π=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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

IR fixed points in $SU(3)$ gauge Theories

We propose a novel RG method to specify the location of the IR fixed point in lattice gauge theories and apply it to the $SU(3)$ gauge theories with $N_f$ fundamental fermions. It is based on the scaling behavior of the propagator through the RG analysis with a finite IR cut-off, which we cannot remove in the conformal field theories in sharp contrast with the confining theories. The method also enables us to estimate the anomalous mass dimension in the continuum limit at the IR fixed point. We perform the program for $N_f=16, 12, 8 $ and $N_f=7$ and indeed identify the location of the IR fixed points in all cases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mass and Axial current renormalization in the Schrödinger functional scheme for the RG-improved gauge and the stout smeared $O(a)$-improved Wilson quark actions

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 $α=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ödinger functional scheme we evaluate the renormalization factors at $β=1.82$ where large scale simulations are being carried out.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conformal Behavior in QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conformal Theories with an IR cutoff

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 Λ_{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^{α_H} instead of the exponentially decaying form c_H\exp{(-m_H t)}, in the small quark mass region where m_H \le c Λ_{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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global Structure of Conformal Theories in the SU(3) Gauge Theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Charmed baryons at the physical point in 2+1 flavor lattice QCD

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 $Ω$ 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 $Ξ_{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.

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

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

We present the results of 1+1+1 flavor QCD+QED simulation at the physical point, in which the dynamical quark effects in QED and the up-down quark mass difference are incorporated by the reweighting technique. The physical quark masses together with the lattice spacing are determined with $m_{π^+}$, $m_{K^+}$, $m_{K^0}$ and $m_{Ω^-}$ 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 $β=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.

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
preprint2012arXiv

Multi-block/multi-core SSOR preconditioner for the QCD quark solver for K computer

We study the algorithmic optimization and performance tuning of the Lattice QCD clover-fermion solver for the K computer. We implement the Lüscher'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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charm quark system at the physical point of 2+1 flavor lattice QCD

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}(μ= 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 $α_s^2 f(m_Q a)(a Λ_{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 $α_s^2 f(m_Q a)(a Λ_{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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modified Block BiCGSTAB for Lattice QCD

We present results for application of block BiCGSTAB algorithm modified by the QR decomposition and the SAP preconditioner to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD on a $32^3 \times 64$ lattice at almost physical quark masses. The QR decomposition improves convergence behaviors in the block BiCGSTAB algorithm suppressing deviation between true residual and recursive one. The SAP preconditioner applied to the domain-decomposed lattice helps us minimize communication overhead. We find remarkable cost reduction thanks to cache tuning and reduction of number of iterations.

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

Non-perturbative renormalization of quark mass in Nf=2+1 QCD with the Schroedinger functional scheme

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_π\sim 500$ MeV and that of PACS-CS collaboration for much lighter quark masses down to $m_π=155$ MeV. The quark mass renormalization factor is used to renormalize bare PCAC masses in these simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Physical Point Simulation in 2+1 Flavor Lattice QCD

We present the results of the physical point simulation in 2+1 flavor lattice QCD with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $β=1.9$ on a $32^3 \times 64$ lattice. The physical quark masses together with the lattice spacing is determined with $m_π$, $m_K$ and $m_Ω$ 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{ö}dinger functional method. The results are compared with the previous ones obtained by the chiral extrapolation method.

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

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

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

We present simulation details and results for the light hadron spectrum in N f = 2 + 1 lattice QCD with the nonperturbatively O(a)-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at a lattice spacing of 0.09 fm on a (2.9fm)^3 box using the PACS-CS computer. We employ the 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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

2+1 Flavor Lattice QCD toward the Physical Point

We present the first results of the PACS-CS project which aims to simulate 2+1 flavor lattice QCD on the physical point with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Numerical simulations are carried out at the lattice spacing of $a=0.0907(13)$fm on a $32^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_π$, $m_K$ and $m_Ω$ 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_π=134.0(4.2)$MeV, $f_K=159.4(3.1)$MeV and $f_K/f_π=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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Charm quark system in 2+1 flavor lattice QCD using the PACS-CS configurations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

SU(2) and SU(3) chiral perturbation theory analyses on meson and baryon masses in 2+1 flavor lattice QCD

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.

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

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

We present preliminary results for the light harden spectrum in $N_f=2+1$ lattice QCD obtained with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at $β=1.90$ on a $32^3 \times 64$ lattice using the PACS-CS computer. We employ Lüscher'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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

2+1 Flavor Lattice QCD with Luescher's Domain-Decomposed HMC Algorithm

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.

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
preprint2005arXiv

Light hadron spectrum and quark masses in 2+1 flavor QCD

CP-PACS and JLQCD collaborations are carrying out a joint project of the 2+1 flavor full QCD simulation. Gauge configurations are generated for the non-perturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action using PHMC algorithm at three lattice spacings, $a\sim 0.076$, 0.010 and 0.122 fm, with a fixed physical volume $(2.0 fm)^3$. We present analysis for the light meson spectrum and quark masses in the continuum limit, which are determined using data obtained from the simulations at the two coarser lattices. Our simulations reproduce experimental values of meson masses. The ud and strange quark masses turn out to be $m_{ud}^{\bar{MS}}(μ=2 GeV)=3.34(23) MeV$ and $m_s^{\bar{MS}}(μ=2 GeV)=86.7(5.9) MeV$. We also show preliminary results at our finest lattice spacing for which simulations are still being continued.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The PACS-CS Project

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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Non-perturbative renormalization of meson decay constants in quenched QCD for a renormalization group improved gauge action

Renormalization constants ($Z$-factors) of vector and axial-vector currents are determined non-perturbatively in quenched QCD for a renormalization group improved gauge action and a tadpole improved clover quark action using the Schrödinger functional method. Non-perturbative values of $Z$-factors turn out to be smaller than one-loop perturbative values by $O(15%)$ at lattice spacing of $a^{-1}\approx$ 1 GeV. The pseudoscalar and vector meson decay constants calculated with the non-perturbative $Z$-factors show a much better scaling behavior compared to previous results obtained with tadpole improved one-loop $Z$-factors. In particular, the non-perturbative $Z$-factors normalized at infinite physical volume show that scaling violation of the decay constants are within about 10% up to the lattice spacing $a^{-1}\sim 1$ GeV. The continuum estimates obtained from data in the range $a^{-1}=$ 1 -- 2 GeV agree with those determined from finer lattices ($a^{-1}\sim 2-4$ GeV) with the standard action.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Flavor Singlet Meson Mass in the Continuum Limit in Two-Flavor Lattice QCD

We present results for the mass of the eta-prime meson in the continuum limit for two-flavor lattice QCD, calculated on the CP-PACS computer, using a renormalization-group improved gauge action, and Sheikholeslami and Wohlert's fermion action with tadpole-improved csw. Correlation functions are measured at three values of the coupling constant beta corresponding to the lattice spacing a approx. 0.22, 0.16, 0.11 fm and for four values of the quark mass parameter kappa corresponding to mpi over mrho approx. 0.8, 0.75, 0.7 and 0.6. For each beta, kappa pair, 400-800 gauge configurations are used. The two-loop diagrams are evaluated using a noisy source method. We calculate eta-prime propagators using local sources, and find that excited state contributions are much reduced by smearing. A full analysis for the smeared propagators gives metaprime=0.960(87)+0.036-0.248 GeV, in the continuum limit, where the second error represents the systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations.

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