Search

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., The XXIX Internaitonal Symposium on Lattice Field Theory, Lattice2011 July 11-16, 2011 The Village at Squaw Valley, US

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., The XXIX Internaitonal Symposium on Lattice Field Theory, Lattice2011 July 10-16, 2011 Squaw Lake Tahoe, California

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 Nf = 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., The XXVIII Internaitonal Symposium on Lattice Field Theory, Lattice2010 July 14-19, 2010 Villasimius, Italy

We present the status of our study on the equation of state in 2+1 flavor QCD with nonperturbatively improved Wilson quarks coupled with the RG improved glue. We apply the Tintegration method to non-perturbatively calculate the equation of state by the fixed-scale approach., The XXVIII Internaitonal Symposium on Lattice Field Theory, Lattice2010 July 14-19, 2010 Villasimius, Italy

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 threedimensional 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., The XXVIII Internaitonal Symposium on Lattice Field Theory, Lattice2010 July 14-19, 2010 Villasimius, Sardinia Italy

We report on the status of our study towards the equation of state in 2+1 flavor QCD with improvedWilson 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 Nt 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., The XXVII Internaitonal Symposium on Lattice Field Theory - LAT2009 July 26-31, 2009 Peking University, Beijing, China

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 323 £12, 10, 8, 6, 4 lattices in the high temperature phase with the RG-improved gluon action and 2+1 flavors of the clover-improvedWilson 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 Nf = 2 and Nf = 0 lattice simulations., The XXVII Internaitonal Symposium on Lattice Field Theory, July 26-31, 2009 Peking University, Beijing, China

Screening properties of the quark gluon plasma are studied from the Polyakov-line correlations in lattice QCD simulations with two flavors of improvedWilson quarks at temperatures T=Tpc '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-line correlations and extract screening masses in magnetic and electric sectors. We 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 theN =4 supersymmetric Yang-Mills theory at 1:3 < T=Tpc < 3., The XXVI Internaitonal Symposium on Lattice Field Theory July 14-19 2008 Williamsburg, Virginia, USA

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 mPS/mV = 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., The XXVI Internaitonal Symposium on Lattice Field Theory July 14-19 2008 Williamsburg, Virginia, USA

We present results for mesonic propagators in temporal and spatial direction and for topological properties at T below and above the deconfining transition in quenched QCD. We use anisotropic lattices and Wilson fermions., Comment: 6 pages, 7 figures, Talk given at 16th International Symposium on Lattice Field Theory (LATTICE 98(hightemp)) , Boulder, CO, 13-18 Jul 1998. (Replaced: Fig.4 corrected, further minor modifications in legends and text.)

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 $\beta=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., Comment: 3 pages, 3 figures, Talk presented at Lattice2004(spectrum), Fermilab, June 21-26, 2004

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., 7 pages, 12 figures, proceedings of the 30th International Symposium on Lattice Field Theory, June 24-29, 2012, Cairns, Australia

We study thermodynamic properties of 2+1 flavor QCD with improved Wilson quarks coupled with the RG improved Iwasaki glue, using the fixed scale approach. We present the results for the equation of state, renormalized Polyakov loop, and chiral condensate., Comment: 7 pages, 7 figures, proceedings of the 30th International Symposium on Lattice Field Theory, June 24-29, 2012, Cairns, Australia

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., Comment: 7 pages, 7 figures; XXIX International Symposium on Lattice Field Theory - Lattice 2011

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., Comment: 13 pages, 17 figures, 3 tables

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., Comment: 7 pages, 6 figures, poster presented at the XXVIII International Symposium on Lattice Field Theory, June 14-19 2010, Villasimius, Italy

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., Comment: 7 pages, 11 figures

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., 7 pages, 7 figure, Talk presented at the XXVIII International Symposium on Lattice Field Theory, Lattice2010, Villasimius, Italy, June 2010

We present the status of our study on the equation of state in 2+1 flavor QCD with non-perturbatively improved Wilson quarks coupled with the RG improved glue. We apply the T-integration method to non-perturbatively calculate the equation of state by the fixed-scale approach., Comment: 7 pages, 7 figures, talk presented at the XXVIII International Symposium on Lattice Field Theory, June 14-19 2010, Villasimius, Italy

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., To appear in the proceedigns of 27th International Symposium on Lattice Field Theory (Lattice 2009), Beijing, China, 25-31 July 2009

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$., 7 pages, 5 figures, Talk given at the XXVI International Symposium on Lattice Field Theory (LATTICE 2008), July 14-19, 2008, Williamsburg, Virginia, USA

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\"odinger 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., Comment: 19 pages, 18 eps figures. Corrected address

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 $\beta=1.8$ ($a = 0.2$ fm from $\rho$ 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}}(\mu=2 {GeV}) = 3.11(17)$ [MeV] in the continuum limit., Comment: 33 pages

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 $\xi=a_s/a_t = 2$ are determined as functions of $\beta$ and $\kappa$, 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., Comment: 16 pages, 23 figures, Version accepted for publication in Physical Review D

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 $\Delta H$ in the HMC algorithm at moderate values of $\Delta\tau$. Nature of the spike is studied. We also study finite-size effects and chiral properties of meson masses., Comment: 3 pages, 4 figures, Lattice2002(spectrum)

We study relativistic charmonium on very fine quenched lattices (beta=6.4 and 6.6). We concentrate on the calculation of the hyperfine splitting between eta_c and J/psi, aiming for a controlled continuum extrapolation of this quantity. Results for the eta_c and J/psi wave functions are also presented., Lattice2001(heavyquark), 3 pages

We study the response of the heavy-light meson mass to the chemical potential ($\partial{m}/\partial{\mu}$) at high temperature and at zero chemical potential. Preliminary results show that $\partial{m}/\partial{\mu}|_{\mu=0}$ is negative in the confinement phase and is positive in the deconfinement phase. This may suggest the possibility for the suppression of $J/\psi$ without the quark gluon plasma formation. We study the response of meson masses to the chemical potential ($\partial{m}/\partial{\mu}$) at high temperature and at zero chemical potential on $N_f=2$ lattice with staggered fermions. Preliminary results for the meson composed of different quarks show that $\partial{m}/\partial{\mu}|_{\mu=0}$ is negative in the confinement phase and positive in the deconfinement phase.

Responses to chemical potential of the pseudoscalar meson screening mass and the chiral condensate in lattice QCD are investigated. On a $16 \times 8^2 \times \4$ lattice with two flavors of staggered quarks the first and second responses below and above $T_c$ are evaluated. Different behavior in the low and the high temperature phases are observed, which may be explained as a consequence of the chiral symmetry breaking and restoration., Comment: Lattice2001(hightemp), 3pages, 3figures

We study the response of meson masses to the chemical potential ($\partial{m}/\partial{\mu}$) at high temperature and at zero chemical potential on $N_f=2$ lattice with staggered fermions. Preliminary results for the meson composed of different quarks show that $\partial{m}/\partial{\mu}|_{\mu=0}$ is negative in the confinement phase and positive in the deconfinement phase., Comment: Lattice 2000 (Finite Density), 4 pages, 2 figures

We analyze temporal and spatial meson correlators in quenched lattice QCD at T>0. Above T_c we find different masses and (spatial) screening masses, signals of plasma formation, and indication of persisting mesonic excitations., 4 pages, 5 figures, revised formulation

We study the equation of state in 2+1 flavor QCD with nonperturbatively improved Wilson quarks coupled with the renormalization group-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-Nt 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ρ≃0.63, while the strange quark mass is chosen around the physical point. Although the light quark masses are still heavier than the physical values, our equation of state is roughly consistent with recent results with highly improved staggered quarks at large Nt.

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 (C), 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, the ratio of the two screening masses agrees with that obtained from the dimensionally-reduced effective field theory and the N=4 supersymmetric Yang-Mills theory.