BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
N. Ishizuka contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We present results of our trial calculation of the $K \to ππ$ decay amplitudes with the improved Wilson fermion action. Calculations are carried out with $N_f=2+1$ gauge configurations generated with the Iwasaki gauge action and non-perturbatively $O(a)$-improved Wilson fermion action at $a=0.091\,{\rm fm}$, $m_π=280\,{\rm MeV}$ and $m_K=560\,{\rm MeV} (\sim 2 m_π)$ on a $32^3\times 64$ ($La=2.9 {\rm fm}$) lattice.
We perform a lattice QCD study of the $ρ$ meson decay from the $N_f=2+1$ full QCD configurations generated with a renormalization group improved gauge action and a non-perturbatively $O(a)$-improved Wilson fermion action. The resonance parameters, the effective $ρ\toππ$ coupling constant and the resonance mass, are estimated from the $P$-wave scattering phase shift for the isospin I=1 two-pion system. The finite size formulas are employed to calculate the phase shift from the energy on the lattice. Our calculations are carried out at two quark masses, $m_π=410\,{\rm MeV}$ ($m_π/m_ρ=0.46$) and $m_π=300\,{\rm MeV}$ ($m_π/m_ρ=0.35$), on a $32^3\times 64$ ($La=2.9\,{\rm fm}$) lattice at the lattice spacing $a=0.091\,{\rm fm}$. We compare our results at these two quark masses with those given in the previous works using $N_f=2$ full QCD configurations and the experiment.
We perform a lattice QCD study of the $ρ$ meson decay from the $N_f=2+1$ full QCD configurations generated with a renormalization group improved gauge action and a non-perturbatively $O(a)$-improved Wilson fermion action. The resonance parameters, the effective $ρ\toππ$ coupling constant and the resonance mass, are estimated from the $P$-wave scattering phase shift for the isospin I=1 two-pion system. The finite size formulas are employed to calculate the phase shift from the energy on the lattice. Our calculations are carried out at two quark masses, $m_π=410\,{\rm MeV}$ ($m_π/m_ρ=0.46$) and $m_π=300\,{\rm MeV}$ ($m_π/m_ρ=0.35$), on a $32^3\times 64$ ($La=2.9\,{\rm fm}$) lattice at the lattice spacing $a=0.091\,{\rm fm}$. We compare our results at these two quark masses with those given in the previous works using $N_f=2$ full QCD configurations and the experiment.
We present preliminary results on the $ρ$ meson decay width from $N_f=2+1$ full QCD configurations generated by PACS-CS Collaboration. The decay width is estimated from the $P$-wave scattering phase shift for the isospin $I=1$ two-pion system. The finite size formula presented by Lüscher in the center of mass frame and its extension to non-zero total momentum frame by Rummukainen and Gottlieb are employed for the calculations of the phase shift. Our calculations are carried out at $m_π=410\ {\rm MeV}$ ($m_π/m_ρ=0.46$) and $a=0.091\ {\rm fm}$ on a $32^3\times 64$ ($La=2.9 {\rm fm}$) lattice.
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.
I present derivation of Lüscher's finite size formula for the elastic $Nπ$ and the $NN$ scattering system for several angular momenta from the relativistic quantum field theory.
We present an evaluation of the running coupling constant for Nf=2+1 QCD. The Schroedinger functional scheme is used as the intermediate scheme to carry out non-perturbative running from the low energy region, where physical scale is introduced, to deep in the high energy perturbative region, where conversion to the MS-bar scheme is safely performed. Possible systematic errors due to the use of perturbation theory occur only in the conversion from three-flavor to four-flavor running coupling constant near the charm mass threshold, where higher order terms beyond 5th order in the $β$ function may not be negligible. For numerical simulations we adopted Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy region and three lattice spacings to take the continuum limit at each scale. A physical scale is introduced from the previous Nf=2+1 simulation of the CP-PACS/JL-QCD collaboration, which covered the up-down quark mass range heavier than $m_π\sim 500$ MeV.
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.
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.
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.
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.
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.
We apply chiral perturbation theory to the pseudoscalar meson mass and decay constant data obtained in the PACS-CS Project toward 2+1 flavor lattice QCD simulations with the O(a)-improved Wilson quarks. We examine the existence of chiral logarithms in the quark mass range from m_{ud}=47 MeV down to 6 MeV on a (2.8 fm)^3 box with the lattice spacing a=0.09 fm. Several low energy constants are determined. We also discuss the magnitude of finite size effects based on chiral perturbation theory.
We present a lattice QCD calculation of the $ρ$ meson decay width via the $P$-wave scattering phase shift for the I=1 two-pion system. Our calculation uses full QCD gauge configurations for $N_f=2$ flavors generated using a renormalization group improved gauge action and an improved Wilson fermion action on a $12^3\times24$ lattice at $m_π/m_ρ=0.41$ and the lattice spacing $1/a=0.92 {\rm GeV}$. The phase shift calculated with the use of the finite size formula for the two-pion system in the moving frame shows a behavior consistent with the existence of a resonance at a mass close to the vector meson mass obtained in spectroscopy. The decay width estimated from the phase shift is consistent with the experiment, when the quark mass is scaled to the realistic value.
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.
We present preliminary results on the $ρ$ meson decay width estimated from the scattering phase shift of the I=1 two-pion system. The phase shift is calculated by the finite size formula for non-zero total momentum frame (the moving frame) derived by Rummukainen and Gottlieb, using the $N_f=2$ improved Wilson fermion action at $m_π/m_ρ=0.41$ and $L=2.53 {\rm fm}$.
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.
We apply the UV-filtering preconditioner, previously used to improve the Multi-Boson algorithm, to the Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm. The performance test for the algorithm is given for the plaquette gauge action and the $O(a)$-improved Wilson action at $β=5.2, c_{\mathrm{sw}}=2.02, M_π/M_ρ\sim 0.8$ and 0.7 on a $16^3\times 48$ lattice. We find that the UV-filtering reduces the magnitude of the molecular dynamics force from the pseudo fermion by a factor 3 by tuning the UV-filter parameter. Combining with the multi-time scale molecular dynamics integrator we achieve a factor 2 improvement.
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.
We present a study of $B$-parameters for generic $ΔS=2$ four-quark operators in domain wall QCD. Our calculation covers all the $B$-parameters required to study the neutral kaon mixing in the standard model (SM) and beyond it. We evaluate one-loop renormalization factors of the operators employing the plaquette and Iwasaki gauge actions. Numerical simulations are carried out in quenched QCD with both gauge actions on $16^3\times 32\times 16$ and $24^3\times 32\times 16$ at the lattice spacing $1/a\approx 2$GeV. We investigate the relative magnitudes of the non-SM $B$-parameters to the SM one, which are compared with the previous results obtained with the overlap and the clover quark actions.
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.
We calculate the two-pion wave function in the ground state of the I=2 $S$-wave system and find the interaction range between two pions, which allows us to examine the validity of the necessary condition for the finite-volume method for the scattering length proposed by Lüscher. We work in the quenched approximation employing a renormalization group improved gauge action for gluons and an improved Wilson action for quarks at $1/a=1.207(12) {\rm GeV}$ on $16^3 \times 80$, $20^3 \times 80$ and $24^3 \times 80$ lattices. We conclude that the necessary condition is satisfied within the statistical errors for the lattice sizes $L\ge 24$ ($3.92 {\rm fm}$) when the quark mass is in the range that corresponds to $m_π^2 = 0.273-0.736 {\rm GeV}^2$. We obtain the scattering length with a smaller statistical error from the wave function than from the two-pion time correlator.
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.
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.
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.
We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the Iwasaki gauge action and the Luescher-Weisz gauge action. In particular, we test the choice of boundary counter terms and apply a perturbative procedure for removal of lattice artifacts for the simulation results in the extrapolation procedure. We confirm the universality of the step scaling functions at both weak and strong coupling regions. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.
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.
We present results on the light hadron spectrum and quark mass in two-flavor QCD calculated with small sea quark masses down to $m_{PS}/m_{V}=0.35$. The configurations are generated using the RG improved gauge and tadpole-improved clover quark action at $β=1.8$, where $a^{-1} \simeq 1$ GeV. Chiral extrapolations are made using not only polynomials and ChPT in the continuum but also formulae of Wilson chiral perturbation theory (WChPT) including $O(a^2)$ chiral breaking terms. We examine the viability of WChPT and its influence on quark masses.
We present a lattice QCD calculation of phase shift including the chiral and continuum extrapolations in two-flavor QCD. The calculation is carried out for I=2 S-wave $ππ$ scattering. The phase shift is evaluated for two momentum systems, the center of mass and laboratory systems, by using the finite volume method proposed by Lüscher in the center of mass system and its extension to general systems by Rummukainen and Gottlieb. The measurements are made at three different bare couplings $β= 1.80$, 1.95 and 2.10 using a renormalization group improved gauge and a tadpole improved clover fermion action, and employing a set of configurations generated for hadron spectroscopy in our previous work. The illustrative values we obtain for the phase shift in the continuum limit are $δ$(deg.) $= - 3.50(64)$, $ - 9.5(30)$ and $ - 16.9(64)$ for $\sqrt{s}({\rm GeV})$ $=0.4$, $ 0.6$ and $ 0.8$, which are consistent with experiment.
We present a report on a calculation of scattering length for I=2 $S$-wave two-pion system from two-pion wave function. Calculations are made with an RG-improved action for gluons and improved Wilson action for quarks at $a^{-1}=1.207(12) {\rm GeV}$ on $16^3 \times 80$, $20^3 \times 80$ and $24^3 \times 80$ lattices. We investigate the validity of necessary condition for application of Lüscher's formula through the wave function. We find that the condition is satisfied for lattice volumes $L\ge 3.92 {\rm fm}$ for the quark mass range $m_π^2 = 0.273-0.736 {\rm GeV}^2$. We also find that the scattering length can be extracted with a smaller statistical error from the wave function than with a time correlation function used in previous studies.
We extend the study of the light hadron spectrum and the quark mass in two-flavor QCD to smaller sea quark mass, corresponding to $m_{PS}/m_{V}=0.60$--0.35. Numerical simulations are carried out using the RG-improved gauge action and the meanfield-improved clover quark action at $β=1.8$ ($a = 0.2$ fm from $ρ$ meson mass). We observe that the light hadron spectrum for small sea quark mass does not follow the expectation from chiral extrapolations with quadratic functions made from the region of $m_{PS}/m_{V}=0.80$--0.55. Whereas fits with either polynomial or continuum chiral perturbation theory (ChPT) fails, the Wilson ChPT (WChPT) that includes $a^2$ effects associated with explicit chiral symmetry breaking successfully fits the whole data: In particular, WChPT correctly predicts the light quark mass spectrum from simulations for medium heavy quark mass, such as $m_{PS}/m_V \simgt 0.5$. Reanalyzing the previous data %at $m_{PS}/m_{V}=0.80$--0.55 with the use of WChPT, we find the mean up and down quark mass being smaller than the previous result from quadratic chiral extrapolation by approximately 10%, $m_{ud}^{\bar{\rm MS}}(μ=2 {GeV}) = 3.11(17)$ [MeV] in the continuum limit.
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.
We explore application of the domain wall fermion formalism of lattice QCD to calculate the $K\toππ$ decay amplitudes in terms of the $K\toπ$ and $K\to 0$ hadronic matrix elements through relations derived in chiral perturbation theory. Numerical simulations are carried out in quenched QCD using domain-wall fermion action for quarks and an RG-improved gauge action for gluons on a $16^3\times 32\times 16$ and $24^3\times 32\times 16$ lattice at $β=2.6$ corresponding to the lattice spacing $1/a\approx 2$GeV. Quark loop contractions which appear in Penguin diagrams are calculated by the random noise method, and the $ΔI=1/2$ matrix elements which require subtractions with the quark loop contractions are obtained with a statistical accuracy of about 10%. We confirm the chiral properties required of the $K\toπ$ matrix elements. Matching the lattice matrix elements to those in the continuum at $μ=1/a$ using the perturbative renormalization factor to one loop order, and running to the scale $μ=m_c=1.3$ GeV with the renormalization group for $N_f=3$ flavors, we calculate all the matrix elements needed for the decay amplitudes. With these matrix elements, the $ΔI=3/2$ decay amplitude shows a good agreement with experiment in the chiral limit. The $ΔI=1/2$ amplitude, on the other hand, is about 50--60% of the experimental one even after chiral extrapolation. In view ofthe insufficient enhancement of the $ΔI=1/2$ contribution, we employ the experimental values for the real parts of the decay amplitudes in our calculation of $ε'/ε$. We find that the $ΔI=3/2$ contribution is larger than the $ΔI=1/2$ contribution so that $ε'/ε$ is negative and has a magnitude of order $10^{-4}$. Possible reasons for these unsatisfactory results are discussed.
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.
We present a calculation of the scattering phase shift for the I=2 S-wave pion-pion system in the continuum limit with two-flavor full QCD. Calculations are made at three lattice spacings, using the finite volume method of Lüscher in the center of mass frame, and its extension to the laboratory frame.
We report on a calculation of the light hadron spectrum and quark masses in three-flavor dynamical QCD using the non-perturbatively O(a)-improved Wilson quark action and a renormalization-group improved gauge action. Simulations are carried out on a 16^3 \times 32 lattice at β=1.9, where a^{-1} \simeq 2GeV, with 6 ud quark masses corresponding to m_{pi}/m_{rho} \simeq 0.64-0.77 and 2 s quark masses close to the physical value. We observe that the inclusion of dynamical strange quark brings the lattice QCD meson spectrum to good agreement with experiment. Dynamical strange quarks also lead to a reduction of the uds quark masses by about 15%.
We report on a non-perturbative evaluation of the renormalization factors for the vector and axial-vector currents, $Z_V$ and $Z_A$, in the quenched domain-wall QCD (DWQCD) with plaquette and renormalization group improved gauge actions. We take the Dirichlet boundary condition for both gauge and domain-wall fermion fields on the finite box, and introduce the flavor-chiral Ward-Takahashi identities to calculate the renormalization factors. As a test of the method, we numerically confirm the expected relation that $Z_V \simeq Z_A$ in DWQCD. Employing two different box sizes for the numerical simulations at several values of the gauge coupling constant $g^2$ and the domain-wall height $M$, we extrapolate $Z_V$ to the infinite volume to remove $a/L$ errors. We finally give the interpolation formula of $Z_V$ in the infinite volume as a function of $g^2$ and $M$.
We present a fully non-perturbative determination of the $O(a)$ improvement coefficient $c_{\rm SW}$ in three-flavor dynamical QCD for the RG improved as well as the plaquette gauge actions, using the Schrödinger functional scheme. Results are compared with one-loop estimates at weak gauge coupling.
Renormalization constants of vector ($Z_V$) and axial-vector ($Z_A$) currents are determined non-perturbatively in quenched QCD for an RG-improved gauge action and a tadpole-improved clover quark action using the Schrödinger functional method. Meson decay constants $f_ρ$ and $f_π$ show much better scaling when $Z_V$ and $Z_A$ estimated for infinite physical volume are used instead of $Z$-factors from tadpole-improved one-loop perturbation theory.
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.
We report on our study of two-flavor full QCD on anisotropic lattices using $O(a)$-improved Wilson quarks coupled with an RG-improved glue. The bare gauge and quark anisotropies corresponding to the renormalized anisotropy $ξ=a_s/a_t = 2$ are determined as functions of $β$ and $κ$, which covers the region of spatial lattice spacings $a_s\approx 0.28$--0.16 fm and $m_{PS}/m_V\approx 0.6$--0.9. The calibrations of the bare anisotropies are performed with the Wilson loop and the meson dispersion relation at 4 lattice cutoffs and 5--6 quark masses. Using the calibration results we calculate the meson mass spectrum and the Sommer scale $r_0$. We confirm that the values of $r_0$ calculated for the calibration using pseudo scalar and vector meson energy momentum dispersion relation coincide in the continuum limit within errors. This work serves to lay ground toward studies of heavy quark systems and thermodynamics of QCD including the extraction of the equation of state in the continuum limit using Wilson-type quark actions.
We present a detailed study of the charmonium spectrum using anisotropic lattice QCD. We first derive a tree-level improved clover quark action on the anisotropic lattice for arbitrary quark mass. The heavy quark mass dependences of the improvement coefficients, i.e. the ratio of the hopping parameters $ζ=K_t/K_s$ and the clover coefficients $c_{s,t}$, are examined at the tree level. We then compute the charmonium spectrum in the quenched approximation employing $ξ= a_s/a_t = 3$ anisotropic lattices. Simulations are made with the standard anisotropic gauge action and the anisotropic clover quark action at four lattice spacings in the range $a_s$=0.07-0.2 fm. The clover coefficients $c_{s,t}$ are estimated from tree-level tadpole improvement. On the other hand, for the ratio of the hopping parameters $ζ$, we adopt both the tree-level tadpole-improved value and a non-perturbative one. We calculate the spectrum of S- and P-states and their excitations. The results largely depend on the scale input even in the continuum limit, showing a quenching effect. When the lattice spacing is determined from the $1P-1S$ splitting, the deviation from the experimental value is estimated to be $\sim$30% for the S-state hyperfine splitting and $\sim$20% for the P-state fine structure. Our results are consistent with previous results at $ξ= 2$ obtained by Chen when the lattice spacing is determined from the Sommer scale $r_0$. We also address the problem with the hyperfine splitting that different choices of the clover coefficients lead to disagreeing results in the continuum limit.
We report on our study of two-flavor full QCD on anisotropic lattices using $O(a)$-improved Wilson quarks coupled with an RG-improved glue. The bare gauge and quark anisotropies corresponding to the renormalized anisotropy $ξ=a_s/a_t = 2$ are determined as functions of $β$ and $κ$, using the Wilson loop and the meson dispersion relation at several lattice cutoffs and quark masses.
We explore the region of small sea quark masses below $m_{PS}/m_V=0.5$ in two-flavor QCD using a mean-field improved clover quark action and an RG-improved gauge action at $a \simeq 0.2$ fm on $12^3 \times 24$ and $16^3 \times 24$ lattices. We find that instability of the standard BiCGStab algorithm at small quark masses can be mostly removed by the BiCGStab(DS-$L$) algorithm, which employs $L$-th minimal residual polynomials with a dynamical selection of $L$. We also find singular spikes of $ΔH$ in the HMC algorithm at moderate values of $Δτ$. Nature of the spike is studied. We also study finite-size effects and chiral properties of meson masses.
We present results of phase shift for I=2 $S$-wave $ππ$ system with the Wilson fermions in the quenched approximation. The finite size method proposed by Lüscher is employed, and calculations are carried out at $β=5.9$ ($a^{-1}=1.934(16)$ GeV from $m_ρ$) on $24^3 \times 60$, $32^3 \times 60$, and $48^3 \times 60$ lattices.
We present our final results of the charmonium spectrum in quenched QCD on anisotropic lattices. Simulations are made with the plaquette gauge action and a tadpole improved clover quark action employing $ξ= a_s/a_t = 3$. We calculate the spectrum of S- and P-states and their excitation, and study the scaling behavior of mass splittings. Comparison is made with the experiment and previous lattice results. The issue of hyperfine splitting for different choices of the clover coefficients obtained by Klassen is discussed.
We investigate a chiral property of the domain-wall fermion (DWF) system using the four-dimensional hermitian Wilson-Dirac operator $H_W$. A formula expressing the Ward-Takahashi identity quark mass $m_{5q}$ with eigenvalues of this operator is derived, which well explains the $N_5$ dependence of $m_{5q}$ observed in previous numerical simulations. We further discuss the chiral property of DWF in the large volume in terms of the spectra of $H_W$.
We present results for the equation of state for pure SU(3) gauge theory obtained on anisotropic lattices with the anisotropy $ξ\equiv a_s/a_t = 2$. The pressure and energy density are calculated on $N_t / ξ= 4, 5$ and 6 lattices with the integral method. They are found to satisfy the leading $1/N_t^2$ scaling from our coarsest lattice $N_t/ξ=4$. This enables us to carry out well controlled continuum extrapolations. We find that the pressure and energy density agree with those obtained using the isotropic plaquette action, but have smaller and more reliable errors.
We present results of a first study of equation of state in finite-temperature QCD with two flavors of Wilson-type quarks. Simulations are made on lattices with temporal size $N_t=4$ and 6, using an RG-improved action for the gluon sector and a meanfield-improved clover action for the quark sector. The lines of constant physics corresponding to fixed values of the ratio $m_{\rm PS}/m_{\rm V}$ of the pseudo-scalar to vector meson masses at zero temperature are determined, and the beta functions which describe the renormalization-group flow along these lines are calculated. Using these results, the energy density and the pressure are calculated as functions of temperature along the lines of constant physics in the range $m_{\rm PS}/m_{\rm V} = 0.65$--0.95. The quark mass dependence in the equation of state is found to be small for $m_{\rm PS}/m_{\rm V} \simlt 0.8$. Comparison of results for $N_t=4$ and $N_t=6$ lattices show significant scaling violation present in the $N_t=4$ results. At high temperatures the results for $N_t=6$ are quite close to the continuum Stefan-Boltzmann limit, suggesting the possibility of a precise continuum extrapolation of thermodynamic quantities from simulations at $N_t\simgt 6$.
We present a dynamical lattice calculation with 2 flavours for bottomonium states with an additional gluonic excitation. Using improved actions for the quarks and gauge fields at a lattice spacing of $a \approx 0.1$ fm, we find 10.977(61)(62) GeV for the energy of the lowest lying $b\bar bg$-hybrid, where the first error is statistical and the second denotes the systematic uncertainty due to the determination of scale. In a parallel quenched simulation we demonstrate explicitly that vacuum polarisation effects are less than 10% of the splitting with the ground state.
We present preliminary results of scattering length and phase shift for I=2 S-wave $ππ$ system with the Wilson fermions in the quenched approximation. The finite size method presented by Lüscher is employed, and calculations are carried out at $β=5.9$ on a $24^3\times 60$ and $32^3\times 60$ lattice.
We report on a calculation of $B_K$ with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at $β=2.6$ and 2.9 corresponding to $a^{-1}\approx 2$GeV and 3GeV. Effects due to finite fifth dimensional size $N_5$ and finite spatial size $N_σ$ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain $B_K(μ= 2 GeV)= 0.5746(61)(191)$, where the first error is statistical and the second error represents an estimate of scaling violation and ${\cal O}(α^2)$ errors in the renormalization factor added in quadrature, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is consistent, albeit somewhat small, with $B_K(μ= 2 {GeV})= 0.628(42)$ obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for light quark masses are also reported.
We present results of a numerical calculation of lattice QCD with two degenerate flavors of dynamical quarks, identified with up and down quarks, and with a strange quark treated in the quenched approximation. The lattice action and simulation parameters are chosen with a view to carrying out an extrapolation to the continuum limit as well as chiral extrapolations. Gauge configurations are generated with a renormalization-group improved gauge action and a mean field improved clover quark action at three values of $β$ and four sea quark masses. The sizes of lattice are chosen so that the physical spatial size is kept constant. Hadron masses, light quark masses and meson decay constants are measured at five valence quark masses. We also carry out complementary quenched simulations with the same improved actions. The quenched spectrum from this analysis agrees well in the continuum limit with the one of our earlier work using the standard action. We find the two-flavor full QCD meson masses in the continuum limit to be much closer to experimental meson masses than those from quenched QCD. We take these results as manifestations of sea quark effects in two-flavor full QCD. For baryon masses full QCD values for strange baryons are in agreement with experiment, while they differ increasingly with decreasing strange quark content, resulting in a nucleon mass higher than experiment. The pattern suggests finite size effects as a possible origin for this deviation. For light quark masses in the continuum limit we obtain values which are reduced by about 25% compared to the values in quenched QCD. We also present results for decay constants where large scaling violations obstruct a continuum extrapolation. Need for a non-perturbative estimate of renormalization factors is discussed.
Employing the maximum entropy method we extract the spectral functions from meson correlators at four lattice spacings in quenched QCD with the Wilson quark action. We confirm that the masses and decay constants, obtained from the position and the area of peaks, agree well with the results from the conventional exponential fit. For the first excited state, we obtain $m_{π_1} = 660(590)$ MeV, $m_{ρ_1} = 1540(570)$ MeV, and $f_{ρ_1} = 0.085(36)$ in the continuum limit.
We apply the maximum entropy method to extract the spectral functions for pseudoscalar and vector mesons from hadron correlators previously calculated at four different lattice spacings in quenched QCD with the Wilson quark action. We determine masses and decay constants for the ground and excited states of the pseudoscalar and vector channels from position and area of peaks in the spectral functions. We obtain the results, $m_{π_1} = 660(590)$ MeV and $m_{ρ_1} = 1540(570)$ MeV for the masses of the first excited state masses, in the continuum limit of quenched QCD. We also find unphysical states which have infinite mass in the continuum limit, and argue that they are bound states of two doublers of the Wilson quark action. If the interpretation is correct, this is the first time that the state of doublers is identified in lattice QCD numerical simulations.
We present a study of the topological susceptibility in lattice QCD with two degenerate flavors of dynamical quarks. The topological charge is measured on gauge configurations generated with a renormalization group improved gauge action and a mean field improved clover quark action at three values of $β=6/g^2$, corresponding to lattice spacings of $a \approx 0.22$, 0.16 and 0.11 fm, with four sea quark masses at each $β$. The study is supplemented by simulations of pure SU(3) gauge theory with the same gauge action at 5 values of $β$ with lattice spacings 0.09 fm$\simlt a \simlt$0.27 fm. We employ a field theoretic definition of the topological charge together with cooling. For the topological susceptibility in the continuum limit of pure SU(3) gauge theory we obtain $χ_t^{1/4} = 197^{+13}_{-16}$ MeV where the error shows statistical and systematic ones added in quadrature. In full QCD $χ_t$ at heavy sea quark masses is consistent with that of pure SU(3) gauge theory. A decrease of $χ_t$ toward light quark masses, as predicted by the anomalous Ward-Takahashi identity for U(1) chiral symmetry, becomes clearer for smaller lattice spacings. The cross-over in the behavior of $χ_t$ from heavy to light sea quark masses is discussed.
We explore the possibility for an evaluation of non-leptonic $ΔS=1$ $K$ decay amplitudes through the calculation of $K\toπ$ matrix elements using domain-wall QCD. The relation between the physical $K\toππ$ matrix elements and $K\toπ$ matrix elements deduced from chiral perturbation theory is recapitulated. Quenched numerical simulations are performed on an $16^3\times 32\times 16$ lattice at a lattice spacing $a^{-1}\approx 2$GeV for both the standard plaquette gauge action and a renormalization-group improved gauge action, and reasonable signals for $K\toπ$ matrix elements are obtained. Preliminary results are reported on the $K\toππ$ matrix elements, and results from two actions are compared.
We present a calculation of the B and D meson decay constants in lattice QCD with two (Nf=2) flavours of light dynamical quarks, using an O(a)-improved Wilson action for both light and heavy quarks and a renormalization-group improved gauge action. Simulations are made at three values of lattice spacing a=0.22, 0.16, 0.11 fm and four values of sea quark mass in the range m_PS/m_V \~= 0.8-0.6. Our estimate for the continuum values of the decay constants are fBd = 208(10)(11) MeV, fBs = 250(10)(13)(^{+8}_{-0}) MeV, fDd = 225(14)(14) MeV, fDs = 267(13)(17)(^{+10}_{-0}) MeV for Nf=2 where the statistical and systematic errors are separately listed, and the third error for fBs and fDs show uncertainty of determination of strange quark mass. We also carry out a set of quenched simulations using the same action to make a direct examination of sea quark effects. Taking the ratio of results for Nf=2 and Nf=0, we obtain fb^{Nf=2}/fb^{Nf=0} = 1.11(6), fbs^{Nf=2}/fbs^{Nf=0} = 1.14(5), fd^{Nf=2}/\fd^{Nf=0} = 1.03(6), fds^{Nf=2}/\fds^{Nf=0} = 1.07(5). They show a 10-15% increase in the Nf=2 results over those of Nf=0 for the B meson decay constants, while evidence for such a trend is statistically less clear for the D meson decay constants.
We present results for light quark masses from a systematic lattice QCD study with two degenerate flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and a mean-field improved clover quark action for sea quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 and the lattice spacing $a=0.22$--0.11 fm. In the continuum limit we find $m_{ud}^{\bar{MS}}(2 {\rm GeV})=3.44^{+0.14}_{-0.22}$ MeV using the $π$ and $ρ$ meson masses as physical input, and $m_s^{\bar{MS}}(2 {\rm GeV})=88^{+4}_{-6}$ MeV or $90^{+5}_{-11}$ MeV with the $K$ or $ϕ$ meson mass as additional input. The quoted errors represent statistical and systematic combined, the latter including those from continuum and chiral extrapolations, and from renormalization factors. Compared to quenched results, two flavors of dynamical quarks reduce quark masses by about 25%.
Chiral properties of QCD formulated with the domain-wall fermion (DWQCD) are studied using the anomalous quark mass m_{5q} and the spectrum of the 4-dimensional Wilson-Dirac operator. Numerical simulations are made with the standard plaquette gauge action and a renormalization-group improved gauge action. Results are reported on the density of zero eigenvalue obtained with the accumulation method, and a comparison is made with the results for m_{5q}.
We present the final results of the CP-PACS calculation of the light hadron spectrum and quark masses with two flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and a mean-field improved clover quark action for sea quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 and the lattice spacing $a=0.22$--0.11 fm. For the meson spectrum in the continuum limit a clearly improved agreement with experiment is observed compared to the quenched case, demonstrating the importance of sea quark effects. For light quark masses we obtain $m_{ud}^{\bar{MS}}(2GeV)=3.44^{+0.14}_{-0.22}$ MeV and $m_s^{\bar{MS}}(2GeV)=88^{+4}_{-6}$ MeV ($K$-input) and $m_s^{\bar{MS}}(2GeV)=90^{+5}_{-11}$ MeV ($ϕ$-input), which are reduced by about 25% compared to the values in quenched QCD.
We report on a calculation of $B_K$ with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at $β=2.6$ and 2.9 corresponding to $a^{-1}\approx 2$GeV and 3GeV. Effects due to finite fifth dimensional size $N_5$ and finite spatial size $N_σ$ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain $B_K(μ= 2 {GeV})= 0.5746(61)$, where the error is statistical only, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is smaller but consistent with $B_K(μ= 2 {GeV})= 0.628(42)$ obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for strange quark mass are also reported.
We study the finite-temperature phase structure and the transition temperature of QCD with two flavors of dynamical quarks on a lattice with the temporal size $N_t=4$, using a renormalization group improved gauge action and the Wilson quark action improved by the clover term. The region of a parity-broken phase is identified, and the finite-temperature transition line is located on a two-dimensional parameter space of the coupling ($β=6/g^2$) and hopping parameter $K$. Near the chiral transition point, defined as the crossing point of the critical line of the vanishing pion mass and the line of finite-temperature transition, the system exhibits behavior well described by the scaling exponents of the three-dimensional O(4) spin model. This indicates a second-order chiral transition in the continuum limit. The transition temperature in the chiral limit is estimated to be $T_c = 171(4)$ MeV.
We present the results of quenched charmonium spectrum for S- and P-states, obtained by a relativistic heavy quark method on anisotropic lattices. Simulations are carried out using the standard plaquette gauge action and a meanfield-improved clover quark action at $a_t^{-1} = 3$--6 GeV with the renormalized anisotropy fixed to $ξ\equiv a_s/a_t =3$. We study the scaling of our fine and hyperfine mass splittings, and compare with previous results.
We study the effects of two dynamical sea quarks on the spectrum of heavy quarkonia. Within the non-relativistic approach to Lattice QCD we found sizeable changes to the hyperfine splitting, but we could not observe any changes for the fine structure. We also investigated the scaling behaviour of our results for several different lattice spacings.
We investigate effects of action improvement on the light hadron spectrum and the static quark potential in two-flavor QCD for $a^{-1} \approx 1$ GeV and $m_{PS}/m_V = 0.7-0.9$. We compare a renormalization group improved action with the plaquette action for gluons, and the SW-clover action with the Wilson action for quarks. We find a significant improvement in the hadron spectrum by improving the quark action, while the gluon improvement is crucial for a rotationally invariant static potential. We also explore the region of light quark masses corresponding to $m_{PS}/m_V \geq 0.4$ on a 2.7 fm lattice using the improved gauge and quark action. A flattening of the potential is not observed up to 2 fm.
A lattice study of the equation of state for pure SU(3) gauge theory using a renormalization-group (RG) improved action is presented. The energy density and pressure are calculated on a $16^3\times 4$ and a $32^3\times 8$ lattice employing the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%.
We present results for the equation of state for pure SU(3) gauge theory obtained with a renormalization-group (RG) improved action. The energy density and pressure are calculated on a $16^3\times 4$ and a $32^3\times 8$ lattice employing the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%.
We present results on an analysis of the decay constants fB and fBs with two flavours of sea quark. The calculation has been carried out on 3 different bare gauge couplings and 4 sea quark masses at each gauge coupling, with m_pi/m_rho ranging from 0.8 to 0.6. We employ the Fermilab formalism to perform calculations with heavy quarks whose mass is in the range of the b-quark. A comparison with a quenched calculation using the same action is made to elucidate the effects due to the sea quarks.
We report on a study of B mesons on N_f = 2 full QCD configurations using an RG-improved gauge action, NRQCD heavy quark action and tadpole-improved clover light quark action. Results on the heavy-light spectrum and the decay constants from 16^3x32 lattices at a^{-1} ~ 1.5 GeV are presented, and compared with quenched results obtained with the same action combination at matching lattice spacings.
We report on a study of heavy hybrid states using the NRQCD approach on coarse and asymmetric lattices, where we discard vacuum polarisation effects and neglect all spin-correction terms. We find a clear hybrid signal on all our lattices ($a_s= 0.15 ... 0.47$ fm). We have studied in detail the lattice spacing artefacts, finite volume effects and mass dependence. Within the above approximations we predict the hybrid excitation in Charmonium to be 1.323(13) GeV above its ground state. The bottomonium hybrid was found to be 1.542(8) GeV above its ground state.
We present results of our quenched study of the B meson decay constant obtained with a parallel set of simulations with the Wilson and Clover actions at β=5.9, 6.1 and 6.3. Systematic errors associated with the large b-quark mass are analyzed within the Fermilab non-relativistic formalism. As our best estimate in the continuum limit we obtain $f_B=163 \pm 16 MeV$ and $f_{B_s}=175 \pm 18 MeV$ with the Clover action.
The quenched hadron spectrum in the continuum obtained with the Wilson quark action in recent simulations on the CP-PACS is presented. Results for the light quark masses and the QCD scale parameter are reported.
We present first results from a full QCD calculation on the CP-PACS, comparing various actions at $a^{-1} \sim 1 GeV$ and $m_π/m_ρ\approx 0.7$--0.9. We use the plaquette and a renormalization group improved action for the gluons, and the Wilson and the SW-Clover action for quarks. We find that significant improvements in the hadron spectrum results from improving the quarks, while the gluon improvement is required for a rotationally invariant static potential. An ongoing effort towards exploring the chiral limit in full QCD is described.