• The energy-momentum tensor plays an important role in QCD thermodynamics. Its expectation value contains information of the pressure and the energy density as its diagonal part. Further properties like viscosity and specific heat can be extracted from its correlation function. Recently a new method based on the gradient flow was introduced to calculate the energy-momentum tensor on the lattice, and has been successfully applied to quenched QCD. In this paper, we apply the gradient flow method to calculate the energy-momentum tensor in (2+1)-flavor QCD. As the first application of the method with dynamical quarks, we study at a single but fine lattice spacing a=0.07 fm with heavy u and d quarks ($m_\pi/m_\rho=0.63$) and approximately physical s quark. Performing simulations on lattices with Nt=16 to 4, the temperature range of T=174-697 MeV is covered. We find that the results of the pressure and the energy density by the gradient flow method are consistent with the previous results using the T-integration method at T<280 MeV, while the results show disagreement at T>350 MeV (Nt<8), presumably due to the small-Nt lattice artifact of $O((aT)^2)=O(1/N_t^2)$. We also apply the gradient flow method to evaluate the chiral condensate taking advantage of the gradient flow method that renormalized quantities can be directly computed avoiding the difficulty of explicit chiral violation with lattice quarks. We compute the renormalized chiral condensate in the MS-bar scheme at renormalization scale $\mu=2$ GeV with a high precision to study the temperature dependence of the chiral condensate and its disconnected susceptibility. Even with the Wilson-type quark action, we obtain the chiral condensate and its disconnected susceptibility showing a clear signal of pseudocritical temperature at T~190 MeV related to the chiral restoration crossover.
  • We calculate the energy gap (latent heat) and pressure gap between the hot and cold phases of the SU(3) gauge theory at the first order deconfining phase transition point. We perform simulations around the phase transition point with the lattice size in the temporal direction $N_t=6,$ 8 and 12 and extrapolate the results to the continuum limit. The energy density and pressure are evaluated by the derivative method with nonperturabative anisotropy coefficients. We find that the pressure gap vanishes at all values of $N_t$. The spatial volume dependence in the latent heat is found to be small on large lattices. Performing extrapolation to the continuum limit, we obtain $\Delta \epsilon/T^4 = 0.75 \pm 0.17$ and $\Delta (\epsilon -3 p)/T^4 = 0.623 \pm 0.056.$ We also tested a method using the Yang-Mills gradient flow. The preliminary results are consistent with those by the derivative method within the error.
  • We study temperature dependence of the topological susceptibility with the $N_{f}=2+1$ flavors Wilson fermion. We have two major interests in this paper. One is a comparison of gluonic and fermionic definitions of the topological susceptibility. Two definitions are related by the chiral Ward-Takahashi identity but their coincidence is highly non-trivial for the Wilson fermion. By applying the gradient flow both for the gauge and quark fields we find a good agreement of these two measurements. The other is a verification of a prediction of the dilute instanton gas approximation at low temperature region $T_{pc}< T<1.5T_{pc}$, for which we confirm the prediction that the topological susceptibility decays with power $\chi_{t}\propto(T/T_{pc})^{-8}$ for three flavors QCD.
  • The energy-momentum tensor and equation of state are studied in finite-temperature (2+1)-flavor QCD with improved Wilson quarks using the method proposed by Makino and Suzuki based on the gradient flow. We find that the results of the gradient flow are consistent with the previous results using the $T$-integration method at $T \mathrel{\rlap{\raise 0.511ex \hbox{$<$}}{\lower 0.511ex \hbox{$\sim$}}} 280$ MeV ($N_t\mathrel{\rlap{\raise 0.511ex \hbox{$>$}}{\lower 0.511ex \hbox{$\sim$}}}10$), while a disagreement is found at $T \mathrel{\rlap{\raise 0.511ex \hbox{$>$}}{\lower 0.511ex \hbox{$\sim$}}} 350$ MeV ($N_t \mathrel{\rlap{\raise 0.511ex \hbox{$<$}}{\lower 0.511ex \hbox{$\sim$}}} 8$) presumably due to the small-$N_t$ lattice artifact. We also report on the results on the renormalized chiral condensate and its disconnected susceptibility using the method of Hieda and Suzuki. The results show a clear signal of the expected chiral restoration crossover even with Wilson-type quarks which violate the chiral symmetry explicitly.
  • We calculate the energy gap (latent heat) and pressure gap between the hot and cold phases of the SU(3) gauge theory at the first order deconfining phase transition point. We perform simulations around the phase transition point with the lattice size in the temporal direction Nt=6, 8 and 12 and extrapolate the results to the continuum limit. We also investigate the spatial volume dependence. The energy density and pressure are evaluated by the derivative method with non-perturabative anisotropy coefficients. We adopt a multi-point reweighting method to determine the anisotropy coefficients. We confirm that the anisotropy coefficients approach the perturbative values as Nt increases. We find that the pressure gap vanishes at all values of Nt when the non-perturbative anisotropy coefficients are used. The spatial volume dependence in the latent heat is found to be small on large lattices. Performing extrapolation to the continuum limit, we obtain $ \Delta \epsilon/T^4 = 0.75 \pm 0.17 $ and $ \Delta (\epsilon -3 p)/T^4 = 0.623 \pm 0.056.$
  • Chiral phase transition of two flavor QCD at finite quark masses is known to be crossover except near the chiral limit, but it can turn to a first order transition when adding many extra flavors. This property is used to explore the nature of the phase transition of massless two flavor QCD using lattice numerical simulations. The extra heavy flavors being incorporated in the form of the hopping parameter expansion through the reweighting, the number of the extra flavors and their masses appear only in a single parameter, defined by $h$. We determine the critical value of the parameter, at which the first order and the crossover regions are separated, and examine its dependence on the two flavor mass. The lattice calculations are carried out at $N_t$=4, and show that the critical value does not depend on the two flavor mass in the range we have studied ($0.46 \le m_\pi/m_\rho \le 0.66$) and appears to remain finite and positive in the chiral limit, suggesting that the phase transition of massless two flavor QCD is of second order.
  • We perform lattice numerical simulations to study the phase transition of QCD at finite temperature to clarify the nature of the transition of massless two flavor QCD. We investigate QCD with two light and Nf heavy quarks instead of two-flavor QCD, and focus on the light quark mass dependence of the critical heavy mass, below which the transition is of first order. The heavy quarks are incorporated into two flavor configurations in the form of the hopping parameter expansion through the reweighting technique. The nature of the transition is identified by the shape of the constraint effective potential at the critical temperature. Our result indicates that the critical heavy mass remains finite in the chiral limit of the two flavors, suggesting the phase transition of massless two-flavor QCD is of second order.
  • Aiming to understand the phase structure of lattice QCD at nonzero temperature and density, we study the phase transitions of QCD in an extended parameter space, where the number of flavor and quark masses are considered as parameters. Performing simulations of 2 flavor QCD and using the reweighting method, we investigate (2+Nf) flavor QCD at finite density, where two light flavors and Nf massive flavors exist. Calculating probability distribution functions, we determine the critical surface terminating first order phase transitions in the parameter space of the light quark mass, the heavy quark mass and the chemical potential. Through the study of the many flavor system, we discuss the phase structure of QCD at finite density.
  • We discuss the QCD critical point at finite density through the study of many flavor QCD, in which two light flavors and Nf massive flavors exist. Performing simulations of QCD with two flavors of improved Wilson fermions, we calculate probability distribution functions of many flavor QCD at finite temperature and density. The dynamical effects of massive flavors and the chemical potential are added using the reweighting technique. From the shape of the distribution functions, we determine the critical surface separating the first order transition and crossover regions in the parameter space of the light and massive quark masses and the chemical potentials. It is found that the critical massive quark mass becomes larger as the chemical potential increases in (2+Nf) flavor QCD. The indication to the (2+1) flavor QCD is then discussed.
  • 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.