Source author record

S. Ejiri

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

68works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

68 published item(s)

preprint2015arXiv

Multipoint reweighting method and beta functions for the calculation of QCD equation of state

We study a reweighting method aiming at numerical studies of QCD at finite density, in which the conventional Monte-Carlo method cannot be applied directly. One of the most important problems in the reweighting method is the overlap problem. To solve it, we propose to perform simulations at several simulation points and combine their results in the data analyses. In this report, we introduce this multipoint reweighting method and test if the method works well by measuring histograms of physical quantities. Using this method, we calculate the meson masses as continuous functions of the gauge coupling beta and the hopping parameters kappa in QCD at zero density. We then determine lines of constant physics in the (beta, kappa) space and evaluate the derivatives of the lattice spacing with respect to beta and kappa along the lines of constant physics (inverse of the beta functions), which are needed in a calculation of the equation of state.

preprint2015arXiv

Multipoint reweighting method and its applications to lattice QCD

The reweighting method is widely used in numerical studies of QCD, in particular, for the cases in which the conventional Monte-Carlo method cannot be applied directly, e.g., finite density QCD. However, the application range of the reweighing method is restricted due to several problems. One of the most severe problems here is the overlap problem. To solve it, we examine a multipoint reweighting method in which simulations at several simulation points are combined in the data analyses. We systematically study the applicability and limitation of the multipoint reweighting method in two-flavor QCD at zero density. Measuring histograms of physical quantities at a series of simulation points, we apply the multipoint reweighting method to calculate the meson masses as continuous functions of the gauge coupling $β$ and the hopping parameters $κ$. We then determine lines of constant physics and beta functions, which are needed in a calculation of the equation of state at finite temperature.

preprint2015arXiv

Towards the QCD equation of state at the physical point using Wilson fermion

We study the (2+1)-flavor QCD at nonzero temperatures using nonperturbatively improved Wilson quarks of the physical masses by the fixed scale approach. We perform physical point simulations at finite temperatures with the coupling parameters which were adopted by the PACS-CS Collaboration in their studies using the reweighting technique. Zero temperature values are obtained on the PACS-CS configurations which are open to the public on the ILDG/JLDG. Finite temperature configurations are generated with the RHMC algorithm. The lattice sizes are $32^3 \times N_t$ with $N_t=14$, 13, $\cdots$, 4 which correspond to $T \approx 160$--550 MeV. We present results of some basic observables at these temperatures and the status of our calculation of the equation of state.

preprint2014arXiv

Histograms in heavy-quark QCD at finite temperature and density

We study the phase structure of lattice QCD with heavy quarks at finite temperature and density by a histogram method. We determine the location of the critical point at which the first-order deconfining transition in the heavy-quark limit turns into a crossover at intermediate quark masses through a change of the shape of the histogram under variation of coupling parameters. We estimate the effect of the complex phase factor which causes the sign problem at finite density, and show that, in heavy-quark QCD, the effect is small around the critical point. We determine the critical surface in 2+1 flavor QCD in the heavy-quark region at all values of the chemical potential mu including mu=infty.

preprint2013arXiv

Scaling properties of the chiral phase transition in the low density region of two-flavor QCD with improved Wilson fermions

We study scaling behavior of a chiral order parameter in the low density region, performing a simulation of two-flavor QCD with improved Wilson quarks. The scaling behavior of the chiral order parameter defined by a Ward-Takahashi identity agrees with the scaling function of the three-dimensional O(4) spin model at zero chemical potential. We extend the scaling study to finite density QCD. Applying the reweighting method and calculating derivatives of the chiral order parameter with respect to the chemical potential, the scaling properties of the chiral phase transition are discussed in the low density region. We moreover calculate the curvature of the phase boundary of the chiral phase transition in the temperature and chemical potential plane assuming the O(4) scaling relation.

preprint2012arXiv

Application of fixed scale approach to static quark free energies in quenched and 2+1 flavor lattice QCD with improved Wilson quark action

Free energies between static quarks and Debye screening masses in the quark-gluon plasma are studied on the basis of Polyakov-line correlations in lattice simulations of 2+1 flavors QCD with the renormalization-group improved gluon action and the $O(a)$-improved Wilson quark action. We perform simulations at $m_{\rm PS}/m_{\rm V} = 0.63$ (0.74) for light (strange) flavors with lattice sizes of $32^3 \times N_t$ with $N_t=4$--12. We adopt the fixed-scale approach, where temperature can be varied without changing the spatial volume and renormalization factor. We find that, at short distance, the free energies of static quarks in color-singlet channel converge to the static-quark potential evaluated from the Wilson-loop at zero-temperature, in accordance with the expected insensitivity of short distance physics to the temperature. At long distance, the free energies of static quarks approach to twice the single-quark free energies, implying that the interaction between static quarks is fully screened. The screening properties can be well described by the screened Coulomb form with appropriate Casimir factor at high temperature. We also discuss a limitation of the fixed-scale approach at high temperature.

preprint2012arXiv

Equation of state in 2+1 flavor QCD with improved Wilson quarks by the fixed scale approach

We study the equation of state in 2+1 flavor QCD with nonperturbatively improved Wilson quarks coupled with the RG-improved Iwasaki glue. We apply the $T$-integration method to nonperturbatively calculate the equation of state by the fixed-scale approach. With the fixed-scale approach, we can purely vary the temperature on a line of constant physics without changing the system size and renormalization constants. Unlike the conventional fixed-$N_t$ approach, it is easy to keep scaling violations small at low temperature in the fixed scale approach. We study 2+1 flavor QCD at light quark mass corresponding to $m_π/m_ρ\simeq 0.63$, while the strange quark mass is chosen around the physical point. Although the light quark masses are heavier than the physical values yet, our equation of state is roughly consistent with recent results with highly improved staggered quarks at large $N_t$.

preprint2012arXiv

Finite density QCD phase transition in the heavy quark region

We extend our previous study of the QCD phase structure in the heavy quark region to non-zero chemical potentials. To identify the critical point where the first order deconfining transition terminates, we study an effective potential defined by the probability distribution function of the plaquette and the Polyakov loop. The reweighting technique is shown to be powerful in evaluating the effective potential in a wide range of the plaquette and Polyakov loop expectation values. We adopt the cumulant expansion to overcome the sign problem in the calculation of complex phase of the quark determinant. We find that the method provides us with an intuitive and powerful way to study the phase structure. We estimate the location of the critical point at finite chemical potential in the heavy quark region.

preprint2012arXiv

Phase structure of finite temperature QCD in the heavy quark region

We study the quark mass dependence of the finite temperature QCD phase transition in the heavy quark region using an effective potential defined through the probability distribution function of the average plaquette. Performing a simulation of SU(3) pure gauge theory, we first confirm that the distribution function has two peaks indicating that the phase transition is of first order in the heavy quark limit, while the first order transition turns into a crossover as the quark mass decreases from infinity, where the mass dependence of the distribution function is evaluated by the reweighting method combined with the hopping parameter expansion. We determine the endpoint of the first order transition region for N_f=1, 2, 3 and 2+1 cases. The quark mass dependence of the latent heat is also evaluated in the first order transition region.

preprint2012arXiv

Probability distribution functions in the finite density lattice QCD

We study the phase structure of QCD at high temperature and density by lattice QCD simulations adopting a histogram method. We try to solve the problems which arise in the numerical study of the finite density QCD, focusing on the probability distribution function (histogram). As a first step, we investigate the quark mass dependence and the chemical potential dependence of the probability distribution function as a function of the Polyakov loop when all quark masses are sufficiently large, and study the properties of the distribution function. The effect from the complex phase of the quark determinant is estimated explicitly. The shape of the distribution function changes with the quark mass and the chemical potential. Through the shape of the distribution, the critical surface which separates the first order transition and crossover regions in the heavy quark region is determined for the 2+1-flavor case.

preprint2011arXiv

Charmonium spectral functions with the variational method in zero and finite temperature lattice QCD

We propose a method to evaluate spectral functions on the lattice based on a variational method. On a lattice with a finite spatial extent, spectral functions consist of discrete spectra only. Adopting a variational method, we calculate the locations and the heights of spectral functions at low-lying discrete spectra. We first test the method in the case of analytically solvable free Wilson quarks at zero and finite temperatures and confirm that the method well reproduces the analytic results for low-lying spectra. We find that we can systematically improve the results by increasing the number of trial states. We then apply the method to calculate the charmonium spectral functions for S and P-wave states at zero-temperature in quenched QCD and compare the results with those obtained using the conventional maximum entropy method (MEM). The results for the ground state are consistent with the location and the area of the first peak in spectral functions from the MEM, while the variational method leads to a mass which is closer to the experimental value for the first excited state. We also investigate the temperature dependence of the spectral functions for S-wave states below and above $T_c$. We obtain no clear evidences for dissociation of $J/ψ$ and $η_c$ up to 1.4$T_c$.

preprint2011arXiv

Histogram method in finite density QCD with phase quenched simulations

We propose a new approach to finite density QCD based on a histogram method with phase quenched simulations at finite chemical potential. Integrating numerically the derivatives of the logarithm of the quark determinant with respect to the chemical potential, we calculate the reweighting factor and the complex phase of the quark determinant. The complex phase is handled with a cumulant expansion to avoid the sign problem. We examine the applicability of this method.

preprint2011arXiv

Scaling behavior of chiral phase transition in two-flavor QCD with improved Wilson quarks at finite density

We study scaling behavior of a chiral order parameter performing a simulation of two-flavor QCD with improved Wilson quarks. It has been shown that the scaling behavior of the chiral order parameter defined by a Ward-Takahashi identity agrees with the scaling function of the three-dimensional O(4) spin model at zero chemical potential. We extend the scaling study to finite density QCD. Calculating derivatives of the chiral order parameter with respect to the chemical potential in two-flavor QCD, the scaling property of chiral phase transition is discussed in the low density region. We moreover calculate the curvature of the phase boundary of the chirl phase transition in the temperature and chemical potential plane assuming the O(4) scaling relation.

preprint2010arXiv

An application of the variational analysis to calculate the meson spectral functions

We present a new method to calculate meson spectral functions (SPFs) on the lattice based on a variational method. Because, on a finite volume lattice, the meson SPFs have discrete spectra only, a suitable way to extract such discrete signals is needed. Using a variational method, we can calculate several discrete quantities such as the position and the area of spectral peaks for low-lying states. Moreover data accuracy can be improved by increasing the number of basis functions. In this report, we first confirm our method in the free quark case and show that our method works well. Then, we apply the method to a quenched lattice QCD simulation and calculate the charmonium SPFs for S and P-waves at zero temperature. Our results for the ground state are well consistent with the position and the area of the lowest peaks of charmonium SPFs calculated by the conventional maximum entropy method. For first excited states, the signals may be reliablly extracted with our method because the charmonium mass converges to a value close to the experimental one when the number of basis functions is increased. We also investigate the SPFs for S-wave charmonia at below and above $T_c$. Our results suggest that $J/ψ$ and $η_c$ may survive up to 1.4$T_c$.

preprint2010arXiv

Electric and Magnetic Screening Masses at Finite Temperature from Generalized Polyakov-Line Correlations in Two-flavor Lattice QCD

Screenings of the quark-gluon plasma in electric and magnetic sectors are studied on the basis of generalized Polyakov-line correlation functions in lattice QCD simulations with two flavors of improved Wilson quarks. Using the Euclidean-time reflection ($\R$) and the charge conjugation ($\Ca$), electric and magnetic screening masses are extracted in a gauge invariant manner. Long distance behavior of the standard Polyakov-line correlation in the quark-gluon plasma is found to be dictated by the magnetic screening. Also, ratio of the two screening masses agrees with that obtained from the dimensionally-reduced effective field theory and the ${\cal N}=4$ supersymmetric Yang-Mills theory.

preprint2010arXiv

Equation of State and Heavy-Quark Free Energy at Finite Temperature and Density in Two Flavor Lattice QCD with Wilson Quark Action

We study the equation of state at finite temperature and density in two-flavor QCD with the RG-improved gluon action and the clover-improved Wilson quark action on a $ 16^3 \times 4$ lattice. Along the lines of constant physics at $m_{\rm PS}/m_{\rm V} = 0.65$ and 0.80, we compute the second and forth derivatives of the grand canonical partition function with respect to the quark chemical potential $μ_q = (μ_u+μ_d)/2$ and the isospin chemical potential $μ_I = (μ_u-μ_d)/2$ at vanishing chemical potentials, and study the behaviors of thermodynamic quantities at finite $μ_q$ using these derivatives for the case $μ_I=0$. In particular, we study density fluctuations at none-zero temperature and density by calculating the quark number and isospin susceptibilities and their derivatives with respect to $μ_q$. To suppress statistical fluctuations, we also examine new techniques applicable at low densities. We find a large enhancement in the fluctuation of quark number when the density increased near the pseudo-critical temperature, suggesting a critical point at finite $μ_q$ terminating the first order transition line between hadronic and quark gluon plasma phases. This result agrees with the previous results using staggered-type quark actions qualitatively. Furthermore, we study heavy-quark free energies and Debye screening masses at finite density by measuring the first and second derivatives of these quantities for various color channels of heavy quark-quark and quark-anti-quark pairs. The results suggest that, to the leading order of $μ_q$, the interaction between two quarks becomes stronger at finite densities, while that between quark and anti-quark becomes weaker.

preprint2010arXiv

Equation of State for physical quark masses

We calculate the QCD equation of state for temperatures corresponding to the transition region with physical mass values for two degenerate light quark flavors and a strange quark using an improved staggered fermion action (p4-action) on lattices with temporal extent N_tau=8. We compare our results with previous calculations performed at twice larger values of the light quark masses as well as with results obtained from a resonance gas model calculation. We also discuss the deconfining and chiral aspects of the QCD transition in terms of renormalized Polyakov loop, strangeness fluctuations and subtracted chiral condensate. We show that compared to the calculations performed at twice larger value of the light quark mass the transition region shifts by about 5 MeV toward smaller temperatures

preprint2010arXiv

QCD thermodynamics at zero and finite densities with improved Wilson quarks

The WHOT-QCD Collaboration is pushing forward lattice studies of QCD at finite temperatures and densities using improved Wilson quarks. We first present results on QCD at zero and finite densities with two flavors of degenerate quarks (N_F=2 QCD) adopting the conventional fixed-Nt approach. We then report on the status of a study of N_F=2+1 QCD adopting a fixed-scale approach armed with the T-integration method which we have developed.

preprint2010arXiv

The order of the deconfinement phase transition in a heavy quark mass region

We study the quark mass dependence of the QCD phase transition by an effective potential defined through the distribution function of observables. As a test of the method, we study the first order deconfinement phase transition in the heavy quark mass limit and its fate at lighter quark masses. We confirm that the distribution function for the plaquette has two peaks indicating that the phase transition is of first order in the heavy quark limit. We then study the quark mass dependence of the distribution function by a reweighting method combined with the hopping parameter expansion. We find that the first order transition turns into a crossover as the quark mass decreases. We determine the critical point for the cases of $N_f$=1, 2, 3 and 2+1. We find that the probability distribution function provides us with a powerful tool to study the order of transitions.

preprint2009arXiv

Heavy-quark free energy at finite temperature with 2+1 flavors of improved Wilson quarks in fixed scale approach

The free energy between a static quark and an antiquark is studied by using the color-singlet Polyakov-line correlation at finite temperature. We perform simulations on $32^3 \times 12$, 10, 8, 6, 4 lattices in the high temperature phase with the RG-improved gluon action and 2+1 flavors of the clover-improved Wilson quark action. Since the simulations are based on the fixed scale approach that the temperature can be varied without changing the spatial volume and renormalization factor, it is possible to investigate temperature dependence of the heavy-quark free energy without any adjustment of the overall constant. We find that, the heavy-quark free energies at short distance converge to the heavy-quark potential evaluated from the Wilson-loop operator at zero temperature, in accordance with the expected insensitivity of short distance physics to the temperature. At long distance, the heavy-quark free energies approach to twice the single-quark free energies, implying that the interaction between heavy quarks is screened. The Debye screening mass obtained from the long range behavior of the heavy-quark free energy is compared with results of the thermal perturbation theory and those of $N_f=2$ and $N_f=0$ lattice simulations.

preprint2009arXiv

On the magnetic equation of state in (2+1)-flavor QCD

A first study of critical behavior in the vicinity of the chiral phase transition of (2+1)-flavor QCD is presented. We analyze the quark mass and volume dependence of the chiral condensate and chiral susceptibilities in QCD with two degenerate light quark masses and a strange quark. The strange quark mass (m_s) is chosen close to its physical value; the two degenerate light quark masses (m_l) are varied in a wide range 1/80 \le m_l/m_s \le 2/5, where the smallest light quark mass value corresponds to a pseudo-scalar Goldstone mass of about 75 MeV. All calculations are performed with staggered fermions on lattices with temporal extent Nt=4. We show that numerical results are consistent with O(N) scaling in the chiral limit. We find that in the region of physical light quark mass values, m_l/m_s \simeq 1/20, the temperature and quark mass dependence of the chiral condensate is already dominated by universal properties of QCD that are encoded in the scaling function for the chiral order parameter, the magnetic equation of state. We also provide evidence for the influence of thermal fluctuations of Goldstone modes on the chiral condensate at finite temperature. At temperatures below, but close to the chiral phase transition at vanishing quark mass, this leads to a characteristic dependence of the light quark chiral condensate on the square root of the light quark mass.

preprint2009arXiv

Towards the equation of state in 2+1 flavor QCD with improved Wilson quarks in the fixed scale approach

We report on the status of our study towards the equation of state in 2+1 flavor QCD with improved Wilson quarks. To reduce the computational cost which is quite demanding for Wilson-type quarks, we adopt the fixed scale approach, i.e. the temperature T is varied by N_t at fixed lattice spacing. Since the conventional integral method to obtain the pressure is inapplicable at a fixed scale, we adopt the "T-integral method", to calculate the pressure non-perturbatively. Reduction of the computational cost of T=0 simulations thus achieved is indispensable to study EOS in QCD with dynamical quarks.

preprint2008arXiv

Equation of state at finite density in two-flavor QCD with improved Wilson quarks

We study the equation of state in two-flavor QCD at finite temperature and density. Simulations are made with the RG-improved gluon action and the clover-improved Wilson quark action. Along the lines of constant physics for $m_{\rm PS}/m_{\rm V} = 0.65$ and 0.80, we compute the derivatives of the quark determinant with respect to the quark chemical potential $μ_q$ up to the fourth order at $μ_q=0$. We adopt several improvement techniques in the evaluation. We study thermodynamic quantities and quark number susceptibilities at finite $μ_q$ using these derivatives. We find enhancement of the quark number susceptibility at finite $μ_q$, in accordance with previous observations using staggered-type quarks. This suggests the existence of a nearby critical point.

preprint2008arXiv

Magnetic and electric screening masses from Polyakov-loop correlations

Screening properties of the quark gluon plasma are studied from Polyakov-loop correlation in lattice QCD simulations with two flavors of improved Wilson quarks at temperatures $T/\Tpc \simeq 1$--4 where $\Tpc$ is the pseudocritical temperature. Using the Euclidean-time reflection symmetry and the charge conjugation symmetry, we introduce various types of Polyakov-loop correlation functions and extract screening masses in magnetic and electric sectors. We find that the temperature dependence of the screening masses are well described by the weak coupling expansion. We also find that a ratio of the screening masses in the electric sector to the magnetic sector shows qualitative agreement with a prediction from the dimensionally-reduced effective field theory and the N=4 supersymmetric Yang-Mills theory at $1.3 < T/\Tpc < 3$.

preprint2008arXiv

Thermodynamics of SU(3) gauge theory at fixed lattice spacing

We study thermodynamics of SU(3) gauge theory at fixed scales on the lattice, where we vary temperature by changing the temporal lattice size N_t=(Ta_t)^{-1}. In the fixed scale approach, finite temperature simulations are performed on common lattice spacings and spatial volumes. Consequently, we can isolate thermal effects in observables from other uncertainties, such as lattice artifact, renormalization factor, and spatial volume effect. Furthermore, in the EOS calculations, the fixed scale approach is able to reduce computational costs for zero temperature subtraction and parameter search to find lines of constant physics, which are demanding in full QCD simulations. As a test of the approach, we study the thermodynamics of the SU(3) gauge theory on isotropic and anisotropic lattices. In addition to the equation of state, we calculate the critical temperature and the static quark free energy at a fixed scale.

preprint2007arXiv

Heavy-Quark Free Energy, Debye Mass, and Spatial String Tension at Finite Temperature in Two Flavor Lattice QCD with Wilson Quark Action

We study Polyakov loop correlations and spatial Wilson loop at finite Temperature in two-flavor QCD simulations with the RG-improved gluon action and the clover-improved Wilson quark action on a $ 16^3 \times 4$ lattice. From the line of constant physics at $m_{\rm PS}/m_{\rm V}=0.65$ and 0.80, we extract the heavy-quark free energies, the effective running coupling $g_{\rm eff}(T)$ and the Debye screening mass $m_D(T)$ for various color channels of heavy quark--quark and quark--anti-quark pairs above the critical temperature. The free energies are well approximated by the screened Coulomb form with the appropriate Casimir factors at high temperature. The magnitude and the temperature dependence of the Debye mass are compared to those of the next-to-leading order thermal perturbation theory and to a phenomenological formula in terms of $g_{\rm eff}(T)$. We make a comparison between our results with the Wilson quark action and the previous results with the staggered quark action. The spatial string tension is also studied in the high temperature phase and is compared to the next-to-next-leading order prediction in an effective theory with dimensional reduction.

preprint2005arXiv

Heavy quark free energies and screening at finite temperature and density

We study the free energies of heavy quarks calculated from Polyakov loop correlation functions in full 2-flavour QCD using the p4-improved staggered fermion action. A small but finite Baryon number density is included via Taylor expansion of the fermion determinant in the Baryo-chemical potential mu. For temperatures above Tc we extract Debye screening masses from the large distance behaviour of the free energies and compare their mu-dependence to perturbative results.

preprint2005arXiv

The QCD equation of state for two flavours at non-zero chemical potential

We present results of a simulation of 2 flavour QCD on a $16^3\times4$ lattice using p4-improved staggered fermions with bare quark mass $m/T=0.4$. Derivatives of the thermodynamic grand canonical partition function $Z(V,T,μ_u,μ_d)$ with respect to chemical potentials $μ_{u,d}$ for different quark flavours are calculated up to sixth order, enabling estimates of the pressure and the quark number density as well as the chiral condensate and various susceptibilities as functions of $μ_{u,d}$ via Taylor series expansion. Results are compared to high temperature perturbation theory as well as a hadron resonance gas model. We also analyze baryon as well as isospin fluctuations and discuss the relation to the chiral critical point in the QCD phase diagram. We moreover discuss the dependence of the heavy quark free energy on the chemical potential.

preprint2003arXiv

Calculation of Non-Leptonic Kaon Decay Amplitudes from $K\toπ$ Matrix Elements in Quenched Domain-Wall QCD

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.

preprint2003arXiv

The Equation of State for Two Flavor QCD at Non-zero Chemical Potential

We present results of a simulation of QCD on a 4x16^3 lattice with 2 continuum flavors of p4-improved staggered fermion with mass m/T=0.4. Derivatives of the thermodynamic grand potential with respect to quark chemical potential mu_q up to fourth order are calculated, enabling estimates of the pressure, quark number density and associated susceptibilities as functions of mu_q via Taylor series expansion. Discretisation effects associated with various staggered fermion formulations are discussed in some detail. In addition it is possible to estimate the radius of convergence of the expansion as a function of temperature. We also discuss the calculation of energy and entropy densities which are defined via mixed derivatives of the thermodynamic grand potential with respect to the bare couplings and quark masses.

preprint2002arXiv

Charmonium Spectrum from Quenched Anisotropic Lattice QCD

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.

preprint2002arXiv

Light Hadron Spectrum and Quark Masses from Quenched Lattice QCD

We present details of simulations for the light hadron spectrum in quenched QCD carried out on the CP-PACS parallel computer. Simulations are made with the Wilson quark action and the plaquette gauge action on 32^3x56 - 64^3x112 lattices at four lattice spacings (a \approx 0.1-0.05 fm) and the spatial extent of 3 fm. Hadronic observables are calculated at five quark masses (m_{PS}/m_V \approx 0.75 - 0.4), assuming the u and d quarks being degenerate but treating the s quark separately. We find that the presence of quenched chiral singularities is supported from an analysis of the pseudoscalar meson data. We take m_π, m_ρand m_K (or m_ϕ) as input. After chiral and continuum extrapolations, the agreement of the calculated mass spectrum with experiment is at a 10% level. In comparison with the statistical accuracy of 1-3% and systematic errors of at most 1.7% we have achieved, this demonstrates a failure of the quenched approximation for the hadron spectrum: the meson hyperfine splitting is too small, and the octet masses and the decuplet mass splittings are both smaller than experiment. Light quark masses are calculated using two definitions: the conventional one and the one based on the axial-vector Ward identity. The two results converge toward the continuum limit, yielding m_{ud}=4.29(14)^{+0.51}_{-0.79} MeV. The s quark mass depends on the strange hadron mass chosen for input: m_s = 113.8(2.3)^{+5.8}_{-2.9} MeV from m_K and m_s = 142.3(5.8)^{+22.0}_{-0} MeV from m_ϕ, indicating again a failure of the quenched approximation. We obtain Λ_{\bar{MS}}^{(0)}= 219.5(5.4) MeV. An O(10%) deviation from experiment is observed in the pseudoscalar meson decay constants.

preprint2001arXiv

Charmonium spectrum from quenched QCD on anisotropic 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.

preprint2001arXiv

Chiral property of domain-wall fermion from eigenvalues of 4D Wilson-Dirac Operator

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

preprint2001arXiv

Equation of state for pure SU(3) gauge theory on anisotropic lattices

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.

preprint2001arXiv

Equation of state in finite-temperature QCD with two flavors of improved Wilson quarks

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

preprint2001arXiv

Hybrid Quarkonia with Dynamical Sea Quarks

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.

preprint2001arXiv

Kaon B parameter from quenched domain-wall 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)(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.

preprint2001arXiv

Light Hadron Spectroscopy with Two Flavors of Dynamical Quarks on the Lattice

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.

preprint2001arXiv

Topological Susceptibility in Lattice QCD with Two Flavors of Dynamical Quarks

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.

preprint2000arXiv

Calculation of $K\to π$ matrix elements in quenched domain-wall QCD

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.

preprint2000arXiv

Decay constants of B and D mesons from improved relativistic lattice QCD with two flavours of sea quarks

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.

preprint2000arXiv

Dynamical Quark Effects on Light Quark Masses

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

preprint2000arXiv

Eigenvalues of the hermitian Wilson-Dirac operator and chiral properties of the domain-wall fermion

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

preprint2000arXiv

Full QCD Light Hadron Spectrum and Quark Masses: Final Results from CP-PACS

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.

preprint2000arXiv

Kaon B parameter from quenched domain-wall 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.

preprint2000arXiv

Phase structure and critical temperature of two-flavor QCD with a renormalization group improved gauge action and clover improved Wilson quark action

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.

preprint2000arXiv

Quenched charmonium spectrum on anisotropic lattices

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.

preprint1999arXiv

Equation of state for pure SU(3) gauge theory with renormalization group improved action

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

preprint1999arXiv

Equation of state for pure SU(3) gauge theory with RG improved action

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

preprint1999arXiv

Heavy Quarkonia from Anisotropic and Isotropic Lattices

We report on recent results for the spectrum of heavy quarkonia. Using coarse and anisotropic lattices we achieved an unprecedented control over statistical and systematic errors for higher excited states such as exotic hybrid states. In a parallel study on isotropic lattices we also investigate the effect of two dynamical flavours on the spin structure of charmonium and bottomonium for several symmetric lattices.

preprint1999arXiv

Heavy-light decay constants from clover heavy quark action in QCD with two flavors of dynamical quarks

We present results on an analysis of the decay constants f_B and f_Bs 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 detailed comparison with a quenched calculation using the same action is made to elucidate the effects due to the sea quarks.

preprint1999arXiv

Heavy-light decay constants from relativistic Nf=2,0 QCD

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.

preprint1999arXiv

Heavy-light spectrum and decay constant from NRQCD with two flavors of dynamical 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.

preprint1999arXiv

Hybrid Quarkonia on Asymmetric Lattices

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.

preprint1999arXiv

Light hadron spectrum and quark masses in QCD with two flavors of dynamical quarks

We present updated results of the CP-PACS calculation of the light hadron spectrum in $N_{\rm f}=2$ full QCD. Simulations are made with an RG-improved gauge action and a tadpole-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.09 fm. A comparison of the full QCD spectrum with new quenched results, obtained with the same improved action, shows clearly the existence of sea quark effects in vector meson masses. Results for light quark masses in $N_{\rm f}=2$ QCD are also presented.

preprint1999arXiv

Quenched Light Hadron Spectrum

We present results of a large-scale simulation for the flavor non-singlet light hadron spectrum in quenched lattice QCD with the Wilson quark action. Hadron masses are calculated at four values of lattice spacing in the range $a \approx$ 0.1 - 0.05 fm on lattices with a physical extent of 3 fm at five quark masses corresponding to $m_π/m_ρ\approx 0.75$ - 0.4. The calculated spectrum in the continuum limit shows a systematic deviation from experiment, though the magnitude of deviation is contained within 11%. Results for decay constants and light quark masses are also reported.

preprint1999arXiv

Quenched QCD with domain-wall fermions on coarse lattices

We investigate the existence of chiral zero modes at a^{-1} \simeq 1 GeV in quenched domain-wall QCD. Simulations are carried out for the plaquette and an RG-improved gauge actions on a 12^3x24xN_s lattice with N_s=10-50. We find that the pion mass in the chiral limit remains non-vanishing as N_s\to\infty for both gauge actions. Possible origins of this non-vanishing pion mass are discussed.

preprint1998arXiv

Full QCD light hadron spectrum from the CP-PACS

We report on an on-going two-flavor full QCD study on CP-PACS using an RG-improved gauge action and a tadpole-improved SW quark action. Runs are made for three lattice spacings $a^{-1}\approx 0.9$, 1.3, and 2.5 GeV on $12^3\times24$, $16^3\times32$, and $24^3\times48$ lattices. Four sea quark masses having $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 are simulated, for each of which hadron masses are evaluated for valence quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.5. Results for hadron and light quark masses are presented and compared with those obtained in quenched QCD.

preprint1998arXiv

Quenched Light Hadron Spectrum with the Wilson Quark Action: Final Results from CP-PACS

We report the final results of the CP-PACS calculation for the quenched light hadron spectrum with the Wilson quark action. Our data support the presence of quenched chiral singularities, and this motivates us to use mass formulae based on quenched chiral perturbation theory in order to extrapolate hadron masses to the physical point. Hadron masses and decay constants in the continuum limit show unambiguous systematic deviations from experiment. We also report the results for light quark masses.

preprint1998arXiv

The static quark potential in full QCD

We report results on the static quark potential in two-flavor full QCD. The calculation is performed for three values of lattice spacing $a^{-1}\approx 0.9, 1.3$ and 2.5 GeV on $12^3{\times}24, 16^3{\times}32$ and $24^3{\times}48$ lattices respectively, at sea quark masses corresponding to $m_π/m_ρ\approx 0.8-0.6$. An RG-improved gauge action and a tadpole-improved SW clover quark action are employed. We discuss scaling of $m_ρ/\sqrtσ$ and effects of dynamical quarks on the potential.