• 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 analyse the performance of a gas time projection chamber (TPC) as a high-performance gamma-ray telescope and polarimeter in the e$^+$e$^-$ pair creation regime. We use data collected at a gamma-ray beam of known polarisation. The TPC provides two orthogonal projections $(x,z)$ and $(y,z)$ of the tracks induced by each conversion in the gas volume. We use a simple vertex finder in which vertices and pseudo-tracks exiting from them are identified. We study the various contributions to the single-photon angular resolution using Monte Carlo simulations and compare them with the experimental data and find that they are in excellent agreement. The distribution of the azimutal angle of pair conversions shows a bias due to the non-cylindrical-symmetric structure of the detector. This bias would average out for a long duration exposure on a space mission, but for this pencil-beam characterisation we have ensured its accurate simulation by a double systematics control scheme, data taking with the detector rotated at several angles with respect to the beam polarisation direction and systematics control with a non-polarised beam. We measure, for the first time, the polarisation asymmetry of a linearly polarised gamma-ray beam in the low energy pair creation regime. This sub-GeV energy range is critical for cosmic sources as their spectra are power laws which fall quickly as a function of increasing energy. This work could pave the way to extending polarised gamma-ray astronomy beyond the MeV energy regime.
  • 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 report on our study of the D meson semileptonic decays in 2+1 flavor lattice QCD. Gauge ensembles are generated at three lattice cutoffs up to 4.5 GeV and with pion masses as low as 300 MeV. We employ the Moebius domain-wall fermion action for both light and charm quarks. We report our preliminary results for the vector and scalar form factors and discuss their dependence on the momentum transfer, quark masses and lattice spacing.
  • 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 review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.
  • Current $\gamma$-ray telescopes suffer from a gap in sensitivity in the energy range between 100keV and 100MeV, and no polarisation measurement has ever been done on cosmic sources above 1MeV. Past and present e$^+$e$^-$ pair telescopes are limited at lower energies by the multiple scattering of electrons in passive tungsten converter plates. This results in low angular resolution, and, consequently, a drop in sensitivity to point sources below 1GeV. The polarisation information, which is carried by the azimuthal angle of the conversion plane, is lost for the same reasons. HARPO (Hermetic ARgon POlarimeter) is an R\&D program to characterise the operation of a gaseous detector (a Time Projection Chamber or TPC) as a high angular-resolution and sensitivity telescope and polarimeter for $\gamma$ rays from cosmic sources. It represents a first step towards a future space instrument in the MeV-GeV range. We built and characterised a 30cm cubic demonstrator [SPIE 91441M], and put it in a polarised $\gamma$-ray beam at the NewSUBARU accelerator in Japan. Data were taken at photon energies from 1.74MeV to 74MeV, and with different polarisation configurations. We describe the experimental setup in beam. We then describe the software we developed to reconstruct the photon conversion events, with special focus on low energies. We also describe the thorough simulation of the detector used to compare results. Finally we will present the performance of the detector as extracted from this analysis and preliminary measurements of the polarisation asymmetry. This beam-test qualification of a gas TPC prototype in a $\gamma$-ray beam could open the way to high-performance $\gamma$-ray astronomy and polarimetry in the MeV-GeV energy range in the near future.
  • 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.
  • Current {\gamma}-ray telescopes based on photon conversions to electron-positron pairs, such as Fermi, use tungsten converters. They suffer of limited angular resolution at low energies, and their sensitivity drops below 1 GeV. The low multiple scattering in a gaseous detector gives access to higher angular resolution in the MeV-GeV range, and to the linear polarisation of the photons through the azimuthal angle of the electron-positron pair. HARPO is an R&D program to characterise the operation of a TPC (Time Projection Chamber) as a high angular-resolution and sensitivity telescope and polarimeter for {\gamma} rays from cosmic sources. It represents a first step towards a future space instrument. A 30 cm cubic TPC demonstrator was built, and filled with 2 bar argon-based gas. It was put in a polarised {\gamma}-ray beam at the NewSUBARU accelerator in Japan in November 2014. Data were taken at different photon energies from 1.7 MeV to 74 MeV, and with different polarisation configurations. The electronics setup is described, with an emphasis on the trigger system. The event reconstruction algorithm is quickly described, and preliminary measurements of the polarisation of 11 MeVphotons are shown.
  • 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.
  • 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.
  • We report on our calculation of the nucleon axial and tensor charges in 2+1-flavor QCD with dynamical overlap quarks. Gauge ensembles are generated at a single lattice spacing 0.12 fm and at a strange quark mass close to its physical value. We employ the all-mode-averaging technique to calculate the relevant nucleon correlation functions, and the disconnected quark loop is efficiently calculated by using the all-to-all quark propagator. We present our preliminary results for the isoscalar and isovector charges obtained at pion masses $m_\pi$ = 450 and 540 MeV.
  • 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 determine the strong coupling constant $\alpha_s$ from a lattice calculation of vacuum polarization functions (VPF) in three-flavor QCD with dynamical overlap fermions. Fitting lattice data of VPF to the continuum perturbative formula including the operator product expansion, we extract the QCD scale parameter $\Lambda_{\overline{MS}}^{(3)}$. At the $Z$ boson mass scale, we obtain $\alpha_s^{(5)}(M_Z)=0.1181(3)(^{+14}_{-12})$, where the first error is statistical and the second is our estimate of various systematic uncertainties.
  • We calculate the electromagnetic pion form factor in lattice QCD with 2+1 flavors of the dynamical overlap quarks. Up and down quark masses are set below their physical values so that the system is in the so-called epsilon regime with the small size of our lattice ~ 1.8 fm. The finite volume corrections are generally expected to be ~ 100% in the epsilon regime. We, however, find a way to automatically cancel the dominant part of them. Inserting non-zero momenta and taking appropriate ratios of the two and three point functions, we can eliminate the contribution from the zero-momentum pion mode. Then the remaining finite volume effect is a small perturbation from the non-zero modes. Our lattice data agree with this theoretical prediction and the extracted pion charge radius is consistent with the experiment.
  • 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.