• Information about the number of Feynman graphs for a given physical process in a given field theory is especially useful for confirming the result of a Feynman graph generator used in an automatic system of perturbative calculations. A method of counting the number of Feynman graphs with weight of symmetry factor was established based on zero-dimensional field theory, and was used in scalar theories and QED. In this article this method is generalized to more complicated models by direct calculation of generating functions on a symbolic calculating system. This method is applied to QCD with and without counter terms, where many higher order are being calculated automatically.
  • We analyze the measured optical conductivity spectra using the density-functional-theory based electronic structure calculation and density-matrix renormalization group calculation of an effective model. We show that, in contrast to a conventional description, the Bose-Einstein condensation of preformed excitons occurs in Ta$_2$NiSe$_5$, despite the fact that a noninteracting band structure is a band-overlap semimetal rather than a small band-gap semiconductor. The system above the transition temperature is therefore not a semimetal, but rather a state of preformed excitons with a finite band gap. A novel insulator state caused by the strong electron-hole attraction is thus established in a real material.
  • 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 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.
  • 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 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.
  • 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.
  • In this paper we show that the phase mixing of continuum Alfv\'{e}n waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfv\'{e}n and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfv\'{e}n waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfv\'{e}n/slow frequency across the magnetic surfaces and is slower than the Alfv\'{e}n/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.
  • 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 show that finite temperature variational cluster approximation (VCA) calculations on an extended Falicov-Kimball model can reproduce angle-resolved photoemission spectroscopy (ARPES) results on Ta2NiSe5 across a semiconductor-to-semiconductor structural phase transition at 325 K. We demonstrate that the characteristic temperature dependence of the flat-top valence band observed by ARPES is reproduced by the VCA calculation on the realistic model for an excitonic insulator only when the strong excitonic fluctuation is taken into account. The present calculations indicate that Ta2NiSe5 falls in the Bose-Einstein condensation regime of the excitonic insulator state.
  • We investigate the triggering mechanisms of plasma eruptions in the solar atmosphere due to interactions between emerging flux and coronal arcade fields by using two-dimensional MHD simulations. We perform parameter surveys with respect to arcade field height, magnetic field strength, and emerging flux location. Our results show that two possible mechanisms exist, and which mechanism is dominant depends mostly on emerging flux location. One mechanism appears when the location of emerging flux is close to the polarity inversion line (PIL) of an arcade field. This mechanism requires reconnection between the emerging flux and the arcade field, as pointed out by previous studies. The other mechanism appears when the location of emerging flux is around the edge of an arcade field. This mechanism does not require reconnection between the emerging flux and the arcade field but does demand reconnection in the arcade field above the PIL. Furthermore, we found that the eruptive condition for this mechanism can be represented by a simple formula.
  • 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.
  • By using the GRACE-Loop system, we calculate the full $\mathcal{O}(\alpha)$ electroweak radiative corrections to the process $e^+e^- \rightarrow e^+e^- \gamma$, which is important for future investigations at the International Linear Collider (ILC). With the GRACE-Loop system, the calculations are checked numerically by three consistency tests: ultraviolet finiteness, infrared finiteness, and gauge-parameter independence. The results show good numerical stability when quadruple precision is used. In the phenomenological results, we find that the electroweak corrections to the total cross section range from $\sim -4\%$ to $\sim -21\%$ when $\sqrt{s}$ varies from $250$ GeV to $1$ TeV. The corrections also significantly affect the differential cross sections, which are a function of the invariant masses and angles and the final-particle energies. Such corrections will play an important role for the high-precision program at the ILC.
  • The full $\mathcal{O}(\alpha)$ electroweak radiative corrections to $t \bar{t} \gamma$ and $e^-e^+\gamma$ productions at the International Linear Collider (ILC) are presented in this paper. The computation is performed with the help of GRACE-Loop system. In the physical results, we discuss on the cross section, electroweak corrections, and the top quark forward-backward asymmetry ($A_{FB}$) which are the function of the center-of-mass energy.
  • Using exact numerical techniques we investigate the nature of excitonic (electron-hole) bound states and the development of exciton coherence in the one-dimensional half-filled extended Falicov-Kimball model. The ground-state phase diagram of the model exhibits, besides band insulator and staggered orbital ordered phases, an excitonic insulator (EI) with power-law correlations. The criticality of the EI state shows up in the von Neumann entropy. The anomalous spectral function and condensation amplitude provide the binding energy and coherence length of the electron-hole pairs which, on their part, point towards a Coulomb interaction driven crossover from BCS-like electron-hole pairing fluctuations to tightly bound excitons. We show that while a mass imbalance between electrons and holes does not affect the location of the BCS-BEC crossover regime it favors staggered orbital ordering to the disadvantage of the EI. Within the BEC regime the quasiparticle dispersion develops a flat valence-band top in accord with the experimental finding for Ta$_2$NiSe$_5$.
  • 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 report on small-cluster exact-diagonalization calculations which prove the formation of electron-hole pairs (excitons) as prerequisite for spontaneous interlayer phase coherence in bilayer systems described by the extended Falicov-Kimball model. Evaluating the anomalous Green's function and momentum distribution function of the pairs, and thereby analyzing the dependence of the exciton binding energy, condensation amplitude, and coherence length on the Coulomb interaction strength, we demonstrate a crossover between a BCS-like electron-hole pairing transition and a Bose-Einstein condensation of tightly bound preformed excitons. We furthermore show that a mass imbalance between electrons and holes tends to suppress the condensation of excitons.
  • We present the full $\mathcal{O}(\alpha)$ electroweak radiative corrections to the process $e^+e^- \rightarrow t \bar{t} \gamma$ at the International Linear Collider (ILC). The computation is performed with the help of the GRACE-Loop system. We present the total cross-section and the top quark forward-backward asymmetry ($A_{FB}$) as a function of the center-of-mass energy and compare them with the process $e^+e^- \rightarrow t \bar{t}$. We find that the value of $A_{FB}$ in $t \bar{t} \gamma$ production is larger than $A_{FB}$ in $t\bar{t}$ production. It is an important result for the measurement of the top quark forward-backward asymmetry at the ILC. Applying a structure function method, we also subtract the QED correction to gain the genuine weak correction in both the $\alpha$ scheme and the $G_{\mu}$ scheme ($\delta_{W}^{G_{\mu}}$). We obtain numerical values for $\delta_{W}^{G_{\mu}}$ which are changing from 2% to -24% when we vary the center-of-mass energy from 360 GeV to 1 TeV.