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 $mu_q$ up to the fourth order at $mu_q=0$. We adopt several improvement techniques in the evaluation. We study thermodynamic quantities and quark number susceptibilities at finite $mu_q$ using these derivatives. We find enhancement of the quark number susceptibility at finite $mu_q$, in accordance with previous observations using staggered-type quarks. This suggests the existence of a nearby critical point.
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_tsimgt 6$.
We study the curvature of the chiral transition/crossover line between the low-temperature hadronic phase and the high-temperature quark-gluon-plasma phase at low densities, performing simulations of two-flavor QCD with improved Wilson quarks. After confirming that the chiral order parameter defined by a Ward-Takahashi identity is consistent with the scaling of the O(4) universality class at zero chemical potential, we extend the scaling analysis to finite chemical potential to determine the curvature of the chiral transition/crossover line at low densities assuming the O(4) universality. To convert the curvature in lattice units to that of the $T_c(mu_B)$ in physical units, we evaluate the lattice scale applying a gradient flow method. We find $kappa=0.0006(1)$ in the chiral limit, which is much smaller than that obtained in (2+1)-flavor QCD with improved staggered quarks.
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_pi/m_rho 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$.
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.
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.