• We study spatial isovector meson correlators in $N_f=2$ QCD with dynamical domain-wall fermions on $32^3\times 8$ lattices at temperatures $T=220-380$ MeV. We measure the correlators of spin-one ($J=1$) operators including vector, axial-vector, tensor and axial-tensor. Restoration of chiral $U(1)_A$ and $SU(2)_L \times SU(2)_R$ symmetries of QCD implies degeneracies in vector--axial-vector ($SU(2)_L \times SU(2)_R$) and tensor--axial-tensor ($U(1)_A$) pairs, which are indeed observed at temperatures above $T_c$. Moreover, we observe an approximate degeneracy of all $J=1$ correlators with increasing temperature. This approximate degeneracy suggests emergent $SU(2)_{CS}$ and $SU(4)$ symmeries at high temperatures, that mix left- and right-handed quarks.
  • We compute the topological susceptibility $\chi_t$ of lattice QCD with $2+1$ dynamical quark flavors described by the M\"obius domain wall fermion. Violation of chiral symmetry as measured by the residual mass is kept at $\sim$1 MeV or smaller. We measure the fluctuation of the topological charge density in a `slab' sub-volume of the simulated lattice using the method proposed by Bietenholz {\it et al.} The quark mass dependence of $\chi_t$ is consistent with the prediction of chiral perturbation theory, from which the chiral condensate is extracted as $\Sigma^{\overline{\rm MS}} (\mbox{2GeV}) = [274(13)(29)\mbox{MeV}]^3$, where the first error is statistical and the second one is systematic. Combining the results for the pion mass $M_\pi$ and decay constant $F_\pi$, we obtain $\chi_t = 0.229(03)(13)M_\pi^2F_\pi^2$ at the physical point.
  • We study the vector and axial-vector current correlators in perturbative and non-perturbative regimes of QCD. The correlators in Euclidean coordinate space are calculated on the lattice using the M\"obius domain-wall fermion formulation at three lattice spacings covering 0.044--0.080~fm. The dynamical quark effects of $2+1$ light flavors are included. The sum $V+A$ and the difference $V-A$ of the vector ($V$) and axial-vector ($A$) current correlators calculated on the lattice after extrapolating to the physical point agree with those converted from the ALEPH experimental data of hadronic $\tau$ decays. The level of the agreement in the $V+A$ channel is about $1.3\sigma$ or smaller in the region of $|x|\ge0.4$~fm, while that in the $V-A$ channel is about $1.8\sigma$ at $|x|=0.74$~fm and smaller at other distances. We also extract the hiral condensate from the short-distance correlators on the lattice using the PCAC relation. Its result extrapolated to the chiral and continuum limit is compatible with other estimates at low energies.
  • We calculate the form factors of the $K \to \pi l \nu$ semileptonic decays in three-flavor lattice QCD, and study their chiral behavior as a function of the momentum transfer and the Nambu-Goldstone boson masses. Chiral symmetry is exactly preserved by using the overlap quark action, which enables us to directly compare the lattice data with chiral perturbation theory (ChPT). We generate gauge ensembles at a lattice spacing of 0.11fm with four pion masses covering 290-540 MeV and a strange quark mass m_s close to its physical value. By using the all-to-all quark propagator, we calculate the vector and scalar form factors with high precision. Their dependence on m_s and the momentum transfer is studied by using the reweighting technique and the twisted boundary conditions for the quark fields. We compare the results for the semileptonic form factors with ChPT at next-to-next-to leading order in detail. While many low-energy constants appear at this order, we make use of our data of the light meson electromagnetic form factors in order to control the chiral extrapolation. We determine the normalization of the form factors as f_+(0) = 0.9636(36)(+57/-35), and observe reasonable agreement of their shape with experiment.
  • We study the axial U(1) symmetry at finite temperature in two-flavor lattice QCD. Employing the Mobius domain-wall fermions, we generate gauge configurations slightly above the critical temperature Tc with different lattice sizes L = 2-4 fm. Our action allows frequent topology tunneling while keeping good chiral symmetry close enough to that of overlap fermions. This allows us to recover full chiral symmetry by an overlap/domain-wall reweighting. Above the phase transition, a strong suppression of the low-lying modes is observed in both of overlap and domain-wall Dirac spectra. We, however, find a sizable violation of the Ginsparg-Wilson relation in the Mobius domain-wall Dirac eigenmodes, which dominates the signals of the axial U(1) symmetry breaking near the chiral limit. We also find that the use of overlap fermion only in the valence sector is dangerous since it suffers from the artifacts due to partial quenching. Reweighting the Mobius domain-wall fermion determinant to that of the overlap fermion, we observe the axial U(1) breaking to vanish in the chiral limit, which is stable against the changes of the lattice volume and lattice spacing.
  • We determine the renormalization constants for flavor non-singlet fermion bilinear operators of M\"obius domain-wall fermions. The renormalization condition is imposed on the correlation functions in the coordinate space, such that the non-perturbative lattice calculation reproduces the perturbatively calculated counterpart at short distances. The perturbative expansion is precise as the coefficients are available up to $O(\alpha_s^4)$. We employ $2+1$-flavor lattice ensembles at three lattice spacings in the range 0.044--0.080~fm.
  • We present a study of chiral behavior of light meson form factors in QCD with three flavors of overlap quarks. Gauge ensembles are generated at single lattice spacing 0.12 fm with pion masses down to 300 MeV. The pion and kaon electromagnetic form factors and the kaon semileptonic form factors are precisely calculated using the all-to-all quark propagator. We discuss their chiral behavior using the next-to-next-to-leading order chiral perturbation theory.
  • We calculate the spectral function of the QCD Dirac operator using the four-dimensional effective operator constructed from the Mobius domain-wall implementation. We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number. The spectrum in the entire eigenvalue range is obtained with a single set of measurements. Results on 2+1-flavor ensembles with Mobius domain-wall sea quarks at lattice spacing ~ 0.08 fm are shown.
  • We report results of masses and decay constants of light and charmed pseudo-scalar mesons using lattice QCD with M\"obius domain-wall fermions. Using this formulation we are able to compute pseudo-scalar decay constants through the pseudo-scalar density operator as well as with the axial-vector current. Results are shown from several lattice spacings and pion masses between 230 MeV and 500 MeV. We present an analysis of these results at different quark masses to show the chiral properties of the light mesons masses and decay constants.
  • We investigate the correlators of flavor non-singlet bilinear operators calculated on the lattice at short distances. In the continuum theory, non-perturbative effects are encoded in the form of the operator product expansion (OPE). We test the prediction of OPE by comparing lattice results with those in the continuum theory. We also determine the renormalization factors of quark currents.
  • We investigate the phase structure of the SU(3) gauge theory with $N_f=8$ by numerical simulations employing the massless Domain-Wall fermions.Our aim is to study directly the massless quark region, since it is the most important region to clarify the properties of conformal theories. When the number of flavor is within the conformal window, it is claimed recently with Wilson quarks that there is the conformal region at the small quark mass region in the parameter space in addition to the confining phase and the deconfining phase. We study the properties of the conformal region investing the spatial Polyakov loops and the temporal meson propagators. Our data imply that there is the conformal region, and a phase transition between the confining phase and the conformal region takes place. These results are consistent with the claim that the conformal window is between $7$ and $16$. Progress reports on other related studies are also presented.
  • The flavor-singlet component of the eta prime meson is related to the topological structure of the SU(3) gauge field through the chiral anomaly. We perform a 2+1-flavor lattice QCD calculation and demonstrate that the two-point function of a gluonically defined topological charge density after a short Yang-Mills gradient flow contains the propagation of the eta prime meson, by showing that its mass in the chiral and continuum limit is consistent with the experimental value. The gluonic correlator does not suffer from the contamination of the pion contribution, and the clean signal is obtained at significantly lower numerical cost compared to the conventional method with the quark bilinear operators.
  • Calculation of the eta-prime meson mass is a notoriously difficult problem, as it requires evaluation of the disconnected diagram which is costly and noisy. In this work, we use a gluonic operator to extract the eta-prime state after smearing the link variables through the Wilson flow. With this choice, one can avoid a large cancellation of pion contribution between the connected and disconnected diagrams. We obtain the eta-prime meson mass on lattices with three different lattice spacings and two physical volumes, which allow us to estimate its continuum and large volume limits.
  • We study the chiral behavior of the electromagnetic (EM) form factors of pion and kaon in three-flavor lattice QCD. In order to make a direct comparison of the lattice data with chiral perturbation theory (ChPT), we employ the overlap quark action that has exact chiral symmetry. Gauge ensembles are generated at a lattice spacing of 0.11 fm with four pion masses ranging between M_pi \simeq 290 MeV and 540 MeV and with a strange quark mass m_s close to its physical value. We utilize the all-to-all quark propagator technique to calculate the EM form factors with high precision. Their dependence on m_s and on the momentum transfer is studied by using the reweighting technique and the twisted boundary conditions for the quark fields, respectively. A detailed comparison with SU(2) and SU(3) ChPT reveals that the next-to-next-to-leading order terms in the chiral expansion are important to describe the chiral behavior of the form factors in the pion mass range studied in this work. We estimate the relevant low-energy constants and the charge radii, and find reasonable agreement with phenomenological and experimental results.
  • We study the non-perturbative determination of the renormalization constants of flavor non-singlet quark bilinear operators on the lattice. The renormalization condition is imposed on correlation functions of bilinear operators in the coordinate space. The results are converted to the value at 2 GeV in the $\rm\overline{MS}$ scheme by a perturbative matching. The calculation is carried out on gauge configurations generated with the Mobius domain-wall fermions at two lattice spacings $a^{-1} = 2.4$ GeV and $a^{-1} = 3.6$ GeV.
  • Global topological charge decorrelates very slowly or even freezes in fine lattice simulations. On the other hand, its local fluctuations are expected to survive and lead to the correct physical results as long as the volume is large enough. We investigate this issue on recently generated configurations including dynamical domain-wall fermions at lattice spacings a = 0.08 fm and finer. We utilize the Yang-Mills gradient flow to define the topological charge density operator and calculate its long-distance correlation, through which we propose a new method for extracting the topological susceptibility in a sub-volume. This method takes care of the finite volume correction, which reduces the bias caused by the global topological charge. Our lattice data clearly show a shorter auto-correlation time than that of the naive definition using the whole lattice, and are less sensitive to the global topological history. Numerical results show a clear sea-quark mass dependence, which agrees well with the prediction of chiral perturbation theory.
  • We carry out a comparative study among five-dimensional formulations of chirally symmetric fermions about the algorithmic performance, chiral symmetry violation and topological tunneling to find a computationally inexpensive formulation with good chiral symmetry. With our choice of the lattice action, we have launched large-scale simulations on fine lattices aiming at a precision study of light and heavy quark physics. We report on the comparative study, current status of the large-scale simulations, and preliminary results on the residual quark mass and auto-correlation.
  • We investigate the eigenvalues of nearly chiral lattice Dirac operators constructed with five-dimensional implementations. Allowing small violation of the Ginsparg-Wilson relation, the HMC simulation is made much faster while the eigenvalues are not significantly affected. We discuss the possibility of reweighting the gauge configurations generated with domain-wall fermions to those of exactly chiral lattice fermions.
  • We calculate the kaon semileptonic form factors in lattice QCD with three flavors of dynamical overlap quarks. Gauge ensembles are generated at pion masses as low as 290 MeV and at a strange quark mass near its physical value. We precisely calculate relevant meson correlators using the all-to-all quark propagator. Twisted boundary conditions and the reweighting technique are employed to vary the momentum transfer and the strange quark mass. We discuss the chiral behavior of the form factors by comparing with chiral perturbation theory and experiments.
  • We report on our calculation of the kaon semileptonic form factors in Nf=2+1 lattice QCD. Chiral symmetry is exactly preserved by using the overlap quark action for a straightforward comparison with chiral perturbation theory (ChPT). We simulate three pion masses down to 290 MeV at a single lattice spacing of 0.11 fm and at a strange quark mass very close to its physical value. The form factors near zero momentum transfer are precisely calculated by using the all-to-all propagator and twisted boundary conditions. We compare the normalizations and slopes of the form factors with ChPT and experiments.
  • We report on our calculation of pion and kaon form factors in three-flavor QCD using the overlap quark action. Gauge ensembles are generated on a 16^3 \times 48 lattice at a lattice spacing of 0.11 fm with pion masses down to 310 MeV. Connected and disconnected meson correaltors are calculated using the all-to-all quark propagator. We present our preliminary analysis on the chiral behavior of the electromagnetic and scalar form factors as well as a comparison of the shape of the K to pi form factors with experiment.
  • We perform a detailed study of the phase diagram of the lattice Higgs SU(2) model with fixed Higgs field length. Consistently with previsions based on the Fradkin Shenker theorem we find a first order transition line with an endpoint whose position we determined. The diagram also shows cross-over lines: the cross-over corresponding to the pure SU(2) bulk is also present at nonzero coupling with the Higgs field and merges with the one that continues the line of first order transition beyond the critical endpoint. At high temperature the first order line becomes a crossover, whose position moves by varying the temperature.
  • The Higgs SU(2) model with fixed Higgs length is usually believed to have two different phases at high gauge coupling (\beta), separated by a line of first order transitions but not distinuguished by any typical symmetry associated with a local order parameter, as first proved by Fradkin and Shenker. We show that in regions of the parameter space where it is usually supposed to be a first order phase transition only a smooth crossover is in fact present.
  • In QCD with adjoint fermions (aQCD) the deconfining transition takes place at a lower temperature than the chiral transition. We study the two transitions by use of the Polyakov Loop, the monopole order parameter and the chiral condensate. The deconfining transition is first order, the chiral is a crossover. The order parameters for confinement are not affected by the chiral transition. We conclude that the degrees of freedom relevant to confinement are different from those describing chiral symmetry.
  • QCD with fermions in the adjoint representation (aQCD) is a model for which a deconfinement and a chiral phase transition take place at different temperatures. In this work, we present a study of the deconfinement transition in the dual superconductor picture based on the evaluation of an operator which carries magnetic charge. The expectation value of this operator signals monopole condensation and is an order parameter for deconfinement as in the case of fermions in the fundamental representation. We find a sharp first order deconfinement transition. We also study the effects of the chiral transition on the monopole order parameter and find them negligible.