
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 zerodimensional 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
densityfunctionaltheory based electronic structure calculation and
densitymatrix renormalization group calculation of an effective model. We show
that, in contrast to a conventional description, the BoseEinstein condensation
of preformed excitons occurs in Ta$_2$NiSe$_5$, despite the fact that a
noninteracting band structure is a bandoverlap semimetal rather than a small
bandgap 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 electronhole
attraction is thus established in a real material.

We study the vector and axialvector current correlators in perturbative and
nonperturbative regimes of QCD. The correlators in Euclidean coordinate space
are calculated on the lattice using the M\"obius domainwall fermion
formulation at three lattice spacings covering 0.0440.080~fm. The dynamical
quark effects of $2+1$ light flavors are included. The sum $V+A$ and the
difference $VA$ of the vector ($V$) and axialvector ($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 $VA$ channel is about
$1.8\sigma$ at $x=0.74$~fm and smaller at other distances. We also extract
the hiral condensate from the shortdistance 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
threeflavor lattice QCD, and study their chiral behavior as a function of the
momentum transfer and the NambuGoldstone 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 290540 MeV and a strange quark mass m_s close to its physical value.
By using the alltoall 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 nexttonextto leading order in detail. While many
lowenergy 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 twoflavor lattice
QCD. Employing the Mobius domainwall fermions, we generate gauge
configurations slightly above the critical temperature Tc with different
lattice sizes L = 24 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/domainwall
reweighting. Above the phase transition, a strong suppression of the lowlying
modes is observed in both of overlap and domainwall Dirac spectra. We,
however, find a sizable violation of the GinspargWilson relation in the Mobius
domainwall 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 domainwall 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 Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark 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 lowenergy 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 heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle fit. Finally, we review the status of lattice
determinations of the strong coupling constant alpha_s.

We determine the renormalization constants for flavor nonsinglet fermion
bilinear operators of M\"obius domainwall fermions. The renormalization
condition is imposed on the correlation functions in the coordinate space, such
that the nonperturbative 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.0440.080~fm.

We calculate the spectral function of the QCD Dirac operator using the
fourdimensional effective operator constructed from the Mobius domainwall
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+1flavor
ensembles with Mobius domainwall sea quarks at lattice spacing ~ 0.08 fm are
shown.

We report results of masses and decay constants of light and charmed
pseudoscalar mesons using lattice QCD with M\"obius domainwall fermions.
Using this formulation we are able to compute pseudoscalar decay constants
through the pseudoscalar density operator as well as with the axialvector
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 flavorsinglet 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+1flavor lattice QCD calculation and demonstrate that the twopoint
function of a gluonically defined topological charge density after a short
YangMills 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+1flavor 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 allmodeaveraging technique to calculate the
relevant nucleon correlation functions, and the disconnected quark loop is
efficiently calculated by using the alltoall 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 etaprime 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 etaprime 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 etaprime 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 threeflavor 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 alltoall 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 nexttonexttoleading 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
lowenergy 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 crossfield 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 crossfield 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 crossfield 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 domainwall fermions at lattice spacings a = 0.08 fm and
finer. We utilize the YangMills gradient flow to define the topological charge
density operator and calculate its longdistance correlation, through which we
propose a new method for extracting the topological susceptibility in a
subvolume. 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 autocorrelation 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 seaquark 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 FalicovKimball model can reproduce angleresolved
photoemission spectroscopy (ARPES) results on Ta2NiSe5 across a
semiconductortosemiconductor structural phase transition at 325 K. We
demonstrate that the characteristic temperature dependence of the flattop
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 BoseEinstein 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 twodimensional 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 threeflavor 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 socalled 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 nonzero momenta
and taking appropriate ratios of the two and three point functions, we can
eliminate the contribution from the zeromomentum pion mode. Then the remaining
finite volume effect is a small perturbation from the nonzero modes. Our
lattice data agree with this theoretical prediction and the extracted pion
charge radius is consistent with the experiment.

By using the GRACELoop 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 GRACELoop system, the calculations are checked
numerically by three consistency tests: ultraviolet finiteness, infrared
finiteness, and gaugeparameter 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 finalparticle
energies. Such corrections will play an important role for the highprecision
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 GRACELoop system. In the physical results, we discuss on the cross
section, electroweak corrections, and the top quark forwardbackward asymmetry
($A_{FB}$) which are the function of the centerofmass energy.

Using exact numerical techniques we investigate the nature of excitonic
(electronhole) bound states and the development of exciton coherence in the
onedimensional halffilled extended FalicovKimball model. The groundstate
phase diagram of the model exhibits, besides band insulator and staggered
orbital ordered phases, an excitonic insulator (EI) with powerlaw
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 electronhole pairs which, on their
part, point towards a Coulomb interaction driven crossover from BCSlike
electronhole pairing fluctuations to tightly bound excitons. We show that
while a mass imbalance between electrons and holes does not affect the location
of the BCSBEC crossover regime it favors staggered orbital ordering to the
disadvantage of the EI. Within the BEC regime the quasiparticle dispersion
develops a flat valenceband top in accord with the experimental finding for
Ta$_2$NiSe$_5$.

We investigate the eigenvalues of nearly chiral lattice Dirac operators
constructed with fivedimensional implementations. Allowing small violation of
the GinspargWilson 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 domainwall fermions to
those of exactly chiral lattice fermions.

We report on smallcluster exactdiagonalization calculations which prove the
formation of electronhole pairs (excitons) as prerequisite for spontaneous
interlayer phase coherence in bilayer systems described by the extended
FalicovKimball 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 BCSlike
electronhole pairing transition and a BoseEinstein 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
GRACELoop system. We present the total crosssection and the top quark
forwardbackward asymmetry ($A_{FB}$) as a function of the centerofmass
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 forwardbackward 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 centerofmass energy from
360 GeV to 1 TeV.