
The energymomentum 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 energymomentum tensor on the
lattice, and has been successfully applied to quenched QCD.
In this paper, we apply the gradient flow method to calculate the
energymomentum 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=174697 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 Tintegration method at T<280 MeV, while the
results show disagreement at T>350 MeV (Nt<8), presumably due to the smallNt
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 MSbar
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 Wilsontype 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 YangMills 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 WardTakahashi
identity but their coincidence is highly nontrivial 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 energymomentum tensor and equation of state are studied in
finitetemperature (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
Wilsontype 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 nonperturabative anisotropy coefficients. We adopt a multipoint
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
nonperturbative 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
twoflavor 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 twoflavor 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 MonteCarlo 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.