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 measure correlation functions of the nonperturbatively renormalized energy-momentum 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 energy-momentum 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 Ward-Takahashi 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 three-flavor 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 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.

Recently, we proposed a novel method to define and calculate the energy-momentum tensor (EMT) in lattice gauge theory on the basis of the Yang-Mills 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 joint-collaboration 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/non-singlet scalar/pseudo-scalar 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 infra-red 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 ultra-violet(UV) divergences, we work on a lattice, employing the overlap Dirac operator, which ensures the exact "chiral" symmetry at finite lattice spacings. Studying multi-point 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 non-chiral lattice Dirac operator is employed are briefly discussed.

We performed a rough estimate of the non-perturbative 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 ultra-violet(UV) divergences. Studying multi-point 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 low-cost digital holographic microscopy (DHM) by adopting an in-line 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 recording-distance 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 low-cost 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 four-quark 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 four-quark 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 non-perturbative 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 CP-PACS/JLQCD collaboration, which covered the up-down quark mass range heavier than m_pi=500 MeV and that of PACS-CS 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 non-perturbative 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 non-perturbatively 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 CP-PACS/JL-QCD collaboration, which covered the up-down quark mass range heavier than $m_\pi\sim 500$ MeV, and that of PACS-CS 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 non-perturbatively 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 non-perturbative renormalization factors at hadronic scale for $\Delta S=2$ four-quark operators in quenched domain-wall QCD using the Schr\"{o}dinger functional method. Combining them with the non-perturbative renormalization group running by the Alpha collaboration, our result yields the fully non-perturbative 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 CP-PACS 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 2-loop 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 4-loop running) in the continuum limit.

We present non-perturbative renormalization factors for $\Delta S=2$ four-quark operators in quenched domain-wall QCD using the Schroedinger functional method. Non-perturbative renormalization factor for $B_K$ is evaluated at hadronic scale. Combined with the non-perturbative 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 CP-PACS collaboration with the quenched domain-wall 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 non-perturbative 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 Yang-Mills 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 domain-wall fermion. In this paper we propose an orbifolding projection procedure to impose the Schroedinger functional Dirichlet boundary condition on the domain-wall 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 Ginsparg-Wilson fermion including the overlap Dirac operator and the domain-wall 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 Yang-Mills 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 overlap-Dirac operator for a fermion sector. The fermion determinant is real and, moreover, when the overlap-Dirac operator is used, semi-positive 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 one-loop 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 domain-wall fermion using an orbifolding projection. With this procedure the zero mode problem with Dirichlet boundary condition can easily be avoided.