
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 measure correlation functions of the nonperturbatively renormalized
energymomentum tensor in $N_f=2+1$ full QCD at finite temperature by applying
the gradient flow method both to the gauge and quark fields. Our main interest
is to study the conservation law of the energymomentum tensor and to test
whether the linear response relation is properly realized for the entropy
density. By using the linear response relation we calculate the specific heat
from the correlation function. We adopt the nonperturbatively improved Wilson
fermion and Iwasaki gauge action at a fine lattice spacing $=0.07$ fm. In this
paper the temperature is limited to a single value $T=232$ MeV. The $u$, $d$
quark mass is rather heavy with $m_\pi/m_\rho=0.63$ while the $s$ quark mass is
set to approximately its physical value.

We compute the topological charge and its susceptibility in finite
temperature (2+1)flavor QCD on the lattice applying a gradient flow method.
With the Iwasaki gauge action and nonperturbatively $O(a)$improved Wilson
quarks, we perform simulations on a fine lattice
with~$a\simeq0.07\,\mathrm{fm}$ at a heavy $u$, $d$ quark mass with
$m_\pi/m_\rho\simeq0.63$ but approximately physical $s$ quark mass with
$m_{\eta_{ss}}/m_\phi\simeq0.74$. In a temperature range
from~$T\simeq174\,\mathrm{MeV}$ ($N_t=16$) to $697\,\mathrm{MeV}$ ($N_t=4$), we
study two topics on the topological susceptibility. One is a comparison of
gluonic and fermionic definitions of the topological susceptibility. Because
the two definitions are related by chiral WardTakahashi identities, their
equivalence is not trivial for lattice quarks which violate the chiral symmetry
explicitly at finite lattice spacings. The gradient flow method enables us to
compute them without being bothered by the chiral violation. We find a good
agreement between the two definitions with Wilson quarks. The other is a
comparison with a prediction of the dilute instanton gas approximation, which
is relevant in a study of axions as a candidate of the dark matter in the
evolution of the Universe. We find that the topological susceptibility shows a
decrease in $T$ which is consistent with the predicted $\chi_\mathrm{t}(T)
\propto (T/T_{\rm pc})^{8}$ for threeflavor QCD even at low temperature
$T_{\rm pc} < T\le1.5 T_{\rm pc}$.

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.

Recently, we proposed a novel method to define and calculate the
energymomentum tensor (EMT) in lattice gauge theory on the basis of the
YangMills gradient flow [1]. In this proceedings, we summarize the basic idea
and technical steps to obtain the bulk thermodynamic quantities in lattice
gauge theory using this method for the quenched and $(2+1)$flavor QCD. The
revised results of integration measure (trace anomaly) and entropy density of
the quenched QCD with corrected coefficients are shown. Furthermore, we also
show the flow time dependence of the parts of EMT including the dynamical
fermions. This work is based on a jointcollaboration between FlowQCD and WHOT
QCD.

We study the finite density phase transition in the lattice QCD at real
chemical potential. We adopt canonical approach and the canonical partition
function is constructed for Nf=2 QCD. After derivation of the canonical
partition function we calculate observables like the pressure, the quark number
density, its second cumulant and the chiral condensate as a function of the
real chemical potential. We covered a wide range of temperature region starting
from the confining low to the deconfining high temperature. We observe signals
for the deconfinement and the chiral restoration phase transition at real
chemical potential below Tc starting from the confining phase.

The canonical partition function is related to the grand canonical one
through the fugacity expansion and is known to have no sign problem. In this
paper we perform the fugacity expansion by a method of the hopping parameter
expansion in temporal direction for the lattice QCD: winding number expansion.
The canonical partition function is constructed for Nf=2 QCD starting from
gauge configurations at zero chemical potential. After derivation of the
canonical partition function we calculate hadronic observables like chiral
condensate and quark number density and the pressure at the real chemical
potential.

In the previous work, we have shown that the SU(2) chiral symmetry recovered
above the critical temperature gives a strong constraint on the Dirac
eigenvalue spectrum and this constraint is strong enough for a set of anomalous
U(1) chiral symmetry breaking operators to vanish in the thermodynamical and
chiral limits. We use this condition as an input and impose a constraint on the
Landau low energy effective theory of QCD. The only constraint we can set is
that the mass splitting term between the pion and eta meson should vanish. All
the singlet/nonsinglet scalar/pseudoscalar mesons contribute to the effective
theory. We evaluate the renormalization group $\beta$function for the
effective theory using the $\epsilon$expansion at one loop level, but find no
stable infrared fixed point except for the trivial Gaussian one. The chiral
phase transition seems to be of first order.

We reconsider constraints on the eigenvalue density of the Dirac operator in
the chiral symmetric phase of 2 flavor QCD at finite temperature. To avoid
possible ultraviolet(UV) divergences, we work on a lattice, employing the
overlap Dirac operator, which ensures the exact "chiral" symmetry at finite
lattice spacings. Studying multipoint correlation functions in various
channels and taking their thermodynamical limit (and then taking the chiral
limit), we obtain stronger constraints than those found in the previous
studies: both the eigenvalue density at the origin and its first and second
derivatives vanish in the chiral limit of 2 flavor QCD. In addition we show
that the axial U(1) anomaly becomes invisible in susceptibilities of scalar and
pseudo scalar mesons, suggesting that the 2nd order chiral phase transition
with the O(4) scaling is not realized in 2 flavor QCD. Possible lattice
artifacts when nonchiral lattice Dirac operator is employed are briefly
discussed.

We performed a rough estimate of the nonperturbative value of the clover
term coefficient cSW for the APE stout link Wilson fermion. We varied the
number of smearings from Nsmear=1 to 6 and adopted beta values roughly
corresponding to the lattice spacing of 0.1 fm. We used the Schroedinger
functional technique for an evaluation of cSW and found that cSW decreases
monotonically as we increase Nsmear but has a 10% order of deviation from the
tree level value for Nsmear=6.

We reinvestigate constraints on the eigenvalue density of the Dirac operator
in the chiral symmetric phase of 2 flavor QCD at finite temperature, employing
the overlap Dirac operator with the exact chiral symmetry at finite lattice
spacings to avoid possible ultraviolet(UV) divergences. Studying multipoint
correlation functions in various channels in the thermodynamical limit, we
obtain stronger constraints than those found in the previous studies that not
only the eigenvalue density at the origin but also its first and second
derivatives vanish in the chiral limit of 2 flavor QCD. In addition we show
that the axial U(1) anomaly becomes invisible in susceptibilities of scalar and
pseudo scalar mesons.

This study developed handheld and lowcost digital holographic microscopy
(DHM) by adopting an inline type hologram, a webcam, a high power RGB light
emitting diode (LED), and a pinhole. It cost less than 20,000 yen
(approximately 250 US dollars at 80 yen/dollar), and was approximately 120 mm x
80 mm x 55 mm in size. In addition, by adjusting the recordingdistance of a
hologram, the lateral resolution power at the most suitable distance was 17.5
um. Furthermore, this DHM was developed for use in open source libraries, and
is therefore lowcost and can be easily developed by anyone. In this research,
it is the feature to cut down cost and size and to improve the lateral
resolution power further rather than existing reports. This DHM will be a
useful application in fieldwork, education, and so forth.

The renormalization factor and O(a) improvement coefficient of fourquark
operators are calculated perturbatively for the improved Wilson fermion action
with clover term and the Iwasaki gauge action. With an application to the
$K\to\pi\pi$ decay amplitude in mind, the calculation is restricted to the
parity odd operator, for which the operators are multiplicatively renormalized
without any mixing to operators that have different chiral structures.

Renormalization factors of fourquark operators are perturbatively calculated
for the improved Wilson fermion with clover term and the Iwasaki gauge action.
A main application shall be the $K\to\pi\pi$ decay amplitude and the
calculation is restricted to the parity odd operator, for which the operators
are multiplicatively renormalizable without mixing with wrong operators that
have different chiral structures.

We present an evaluation of the quark mass renormalization factor for Nf=2+1
QCD. The Schroedinger functional scheme is employed as the intermediate scheme
to carry out nonperturbative running from the low energy to deep in the high
energy perturbative region. The regularization independent step scaling
function of the quark mass is obtained in the continuum limit. Renormalization
factors for the pseudo scalar density and the axial vector current are also
evaluated for the same action and the bare couplings as two recent large scale
Nf=2+1 simulations; previous work of the CPPACS/JLQCD collaboration, which
covered the updown quark mass range heavier than m_pi=500 MeV and that of
PACSCS collaboration on the physical point using the reweighting technique.

We present an evaluation of the running coupling constant and the quark mass
renormalization factor for $N_f=2+1$ QCD. The Schr\"odinger functional scheme
is used as the intermediate scheme to carry out nonperturbative running from
the low energy region, where physical input is introduced, to deep in the high
energy perturbative region, where conversion to the ${\ovl{\rm MS}}$ scheme is
safely performed.
For numerical simulations we adopted Iwasaki gauge action and
nonperturbatively improved Wilson fermion action with the clover term. Seven
renormalization scales are used to cover from low to high energy region and
three lattice spacings to take the continuum limit at each scale.
Physical inputs are introduced from the previous $N_f=2+1$ simulation of the
CPPACS/JLQCD collaboration, which covered the updown quark mass range
heavier than $m_\pi\sim 500$ MeV, and that of PACSCS collaboration for much
lighter quark masses down to $m_\pi=155$ MeV.

We present a preliminary result of Nf=2+1 QCD running coupling in
Schroedinger functional scheme. We adopted Iwasaki gauge action and
nonperturbatively improved Wilson fermion action with clover term. We use
seven renormalization scales to cover from low energy to high energy
perturbative region and three lattice spacings to take the continuum limit at
each scale.
A scaling behabior of the step scaling function is discussed together with
its renoralization group flow in the continuum. We argue on introduction of the
physical scale through the Sommer scale r0.

We calculate nonperturbative renormalization factors at hadronic scale for
$\Delta S=2$ fourquark operators in quenched domainwall QCD using the
Schr\"{o}dinger functional method. Combining them with the nonperturbative
renormalization group running by the Alpha collaboration, our result yields the
fully nonperturbative renormalization factor, which converts the lattice bare
$B_K$ to the renormalization group invariant (RGI) $\hat{B}_K$. Applying this
to the bare $B_K$ previously obtained by the CPPACS collaboration at
$a^{1}\simeq 2, 3, 4$ GeV, we obtain $\hat{B}_K=0.782(5)(7)$ (equivalent to
$B_K^{\bar{\rm MS}}({\rm NDR}, 2 {\rm GeV}) = 0.565(4)(5)$ by 2loop running)
in the continuum limit, where the first error is statistical and the second is
systematic due to the continuum extrapolation. Except the quenching error, the
total error we have achieved is less than 2%, which is much smaller than the
previous ones. Taking the same procedure, we obtain $m_{u,d}^{\rm
RGI}=5.613(66)$ MeV and $m_s^{\rm RGI}=147.1(17)$ MeV (equivalent to
$m_{u,d}^{\bar{\rm MS}}(2 {\rm GeV})=4.026(48)$ MeV and $m_{s}^{\bar{\rm MS}}(2
{\rm GeV})=105.6(12)$ MeV by 4loop running) in the continuum limit.

We present nonperturbative renormalization factors for $\Delta S=2$
fourquark operators in quenched domainwall QCD using the Schroedinger
functional method. Nonperturbative renormalization factor for $B_K$ is
evaluated at hadronic scale. Combined with the nonperturbative RG running
obtained by the Alpha collaboration, our result yields renormalization factor
which converts lattice bare $B_K$ to the renormalization group invariant one.
We apply the renormalization factor to bare $B_K$ previously obtained by the
CPPACS collaboration with the quenched domainwall QCD(DWQCD). We compare our
result with previous ones obtained by perturbative renormalization factors,
different renormalization schemes or different quark actions. We also show that
chiral symmetry breaking effects in the renormalization factor are numerically
small.

Finite volume renormalization scheme is one of the most fascinating scheme
for nonperturbative renormalization on lattice.
By using the step scaling function one can follow running of renormalized
quantities with reasonable cost.
It has been established the Schroedinger functional is very convenient to
define a field theory in a finite volume for the renormalization scheme.
The Schroedinger functional, which is characterized by a
Dirichlet boundary condition in temporal direction, is well defined and works
well for the YangMills theory and QCD with the Wilson fermion.
However one easily runs into difficulties if one sets the same sort of the
Dirichlet boundary condition for the overlap Dirac operator or the domainwall
fermion.
In this paper we propose an orbifolding projection procedure to impose the
Schroedinger functional Dirichlet boundary condition on the domainwall
fermion.

The Schroedinger functional formalism is given as a field theory in a finite
volume with a Dirichlet boundary condition in temporal direction. When one
tries to construct this formalism with the GinspargWilson fermion including
the overlap Dirac operator and the domainwall fermion one easily runs into
difficulties. The reason is that if the Dirichlet boundary condition is simply
imposed on the Wilson Dirac operator $DW$ inside of the overlap Dirac operator
an exponentially small eigenvalue appears in $DW$, which affects the locality
properties of the operator. In this paper we propose a new procedure to impose
the Schroedinger functional Dirichlet boundary condition on the overlap Dirac
operator using an orbifolding projection.

The N=(2,2) extended super YangMills theory in 2 dimensions is formulated on
the lattice as a dimensional reduction of a 4 dimensional lattice gauge theory.
We use the plaquette action for a bosonic sector and the Wilson or the
overlapDirac operator for a fermion sector. The fermion determinant is real
and, moreover, when the overlapDirac operator is used, semipositive definite.
The flat directions in the target theory become compact and present no subtlety
for a numerical integration along these directions. Any exact supersymmetry
does not exist in our lattice formulation; nevertheless we argue that oneloop
calculable and finite mass counter terms ensure a supersymmetric continuum
limit to all orders of perturbation theory.

In this proceeding we propose a new procedure to impose the Schroedinger
functional Dirichlet boundary condition on the overlap Dirac operator and the
domainwall fermion using an orbifolding projection.
With this procedure the zero mode problem with Dirichlet boundary condition
can easily be avoided.