
We study the doublewinding Wilson loops in the SU(N) YangMills theory on
the lattice. We discuss how the area law falloff of the doublewinding Wilson
loop average is modified by changing the enclosing contours C1 and C2 for
various values of the number of color N. By using the strong coupling
expansion, we evaluate the doublewinding Wilson loop average in the lattice
SU(N) YangMills theory. Moreover, we compute the doublewinding Wilson loop
average by lattice Monte Carlo simulations for SU(2) and SU(3). We further
discuss the results from the viewpoint of the NonAbelian Stokes theorem in the
higher representations.

Studying SU(3) gauge theories with increasing number of light fermions is
relevant both for understanding the strong dynamics of QCD and for constructing
strongly interacting extensions of the Standard Model (e.g. UV completions of
composite Higgs models). In order to contrast these manyflavors strongly
interacting theories with QCD, we study the flavorsinglet spectrum as an
interesting probe. In fact, some composite Higgs models require the Higgs boson
to be the lightest flavorsinglet scalar in the spectrum of a strongly
interacting new sector with a well defined hierarchy with the rest of the
states. Moreover, introducing many light flavors at fixed number of colors can
influence the dynamics of the lightest flavorsinglet pseudoscalar. We present
the ongoing study of these flavorsinglet channels using multiple
interpolating operators on highstatistics ensembles generated by the LatKMI
collaboration and we compare results with available data obtained by the
Lattice Strong Dynamics collaboration. For the theory with 8 flavors, the two
collaborations have generated configurations that complement each others with
the aim to tackle the massless limit using the largest possible volumes.

We examine how the average of doublewinding Wilson loops depends on the
number of color $N$ in the $SU(N)$ YangMills theory. In the case where the two
loops $C_1$ and $C_2$ are identical, we derive the exact operator relation
which relates the doublewinding Wilson loop operator in the fundamental
representation to that in the higher dimensional representations depending on
$N$. By taking the average of the relation, we find that the
differenceofareas law for the area law falloff recently claimed for $N=2$ is
excluded for $N \geq 3$, provided that the string tension obeys the Casimir
scaling for the higher representations. In the case where the two loops are
distinct, we argue that the area law follows a novel law $(N  3)A_1/(N1)+A_2$
with $A_1$ and $A_2 (A_1<A_2)$ being the minimal areas spanned respectively by
the loops $C_1$ and $C_2$, which is neither sumofareas ($A_1+A_2$) nor
differenceofareas ($A_2  A_1$) law when ($N\geq3$). Indeed, this behavior
can be confirmed in the twodimensional $SU(N)$ YangMills theory exactly.

Based on the highly improved staggered quark action, we perform lattice
simulations of $N_f=8$ QCD and confirm our previous observation of a
flavorsinglet scalar meson (denoted as $\sigma$) as light as the pion and
various "walking signals" through lowlying spectra, with higher statistics,
smaller fermion masses $m_f$, and larger volumes. We measure $M_\pi$, $F_\pi$,
$M_\rho$, $M_{a_0}$, $M_{a_1}$, $M_{b_1}$, $M_N$, $M_\sigma$, $F_\sigma$,
$\langle \bar{\psi} \psi\rangle$ (both directly and through the GMOR relation),
and the string tension. The data are consistent with the spontaneously broken
phase of the chiral symmetry, in agreement with the previous results: ratios of
the quantities to $M_\pi$ monotonically increase in the smaller $m_f$ region
towards the chiral limit similarly to $N_f=4$ QCD, in sharp contrast to
$N_f=12$ QCD where the ratios become flattened. The hyperscaling relation holds
with roughly a universal value of the anomalous dimension, $\gamma_m \simeq 1$,
with a notable exception of $M_\pi$ with $\gamma_m \simeq 0.6$ as in the
previous results. This is a salient feature ("walking signal") of $N_f=8$,
unlike either $N_f=4$ which has no hyperscaling relation at all, or $N_f=12$
QCD which exhibits universal hyperscaling. We further confirm the previous
observation of the light $\sigma$ with mass comparable to the pion in the
studied $m_f$ region. In a chiral limit extrapolation of the $\sigma$ mass
using the dilaton chiral perturbation theory and also using the simple linear
fit, we find the value consistent with the 125 GeV Higgs boson within errors.
Our results suggest that the theory could be a good candidate for walking
technicolor model, having anomalous dimension $\gamma_m \simeq 1$ and a light
flavorsinglet scalar meson as a technidilaton, which can be identified with
the 125 GeV composite Higgs in $N_f=8$ onefamily model.

The dual superconductivity is a promising mechanism for quark confinement. We
have presented a new formulation of the YangMills theory on the lattice that
enables us to change the original nonAbelian gauge field into the new field
variables such that one of them called the restricted field gives the dominant
contribution to quark confinement in the gauge independent way. We have pointed
out that the SU(3) YangMills theory has another reformulation using new field
variables (minimal option), in addition to the way adopted by Cho, Faddeev and
Niemi (maximal option). In the past lattice conferences, we have shown the
numerical evidences that support the nonAbelian dual superconductivity using
the minimal option for the SU(3) YangMills theory. This result should be
compared with Abelian dual superconductivity obtained in the maximal option
which is a gauge invariant extension of the conventional Abelian projection
method in the maximal Abelian gauge.
In this talk, we focus on discriminating between two reformulations, i.e.,
maximal and minimal options of the $SU(3)$ YangMills theory from the viewpoint
of dual superconductivity for quark confinement. We investigate the
confinement/deconfinement phase transitions at finite temperature in both
options, which are compared with the original YangMills theory. For this
purpose, we measure the distribution of Polyakovloops and the Polyakovloop
average, the correlation function of the Polyakov loops and the distribution of
the chromoelectric flux connecting a quark and antiquark in both confinement
and deconfinement phases.

SU(3) gauge theory with eight massless flavours is believed to be walking,
while the corresponding twelve and fourflavour appear IRconformal and
confining respectively. Looking at the simulations performed by the LatKMI
collaboration of these theories, we use the topological susceptibility as an
additional probe of the IR dynamics. By drawing a comparison with SU(3) pure
gauge theory, we see a dynamical quenching effect emerge at larger number of
flavours, which is suggestive of emerging nearconformal and conformal
behaviour.

We have pointed out that the $SU(3)$ YangMills theory has a new way of
reformulation using new field variables (minimal option), in addition to the
conventional option adopted by Cho, Faddeev and Niemi (maximal option). The
reformulation enables us to change the original nonAbelian gauge field into
the new field variables such that one of them called the restricted field gives
the dominant contribution to quark confinement in the gaugeindependent way. In
the minimal option, especially, the restricted field is nonAbelian $U(2)$ and
involves the nonAbelian magnetic monopole. In the preceding lattice
conferences, we have accumulated the numerical evidences for the nonAbelian
magneticmonopole dominance in addition to the restricted nonAbelian field
dominance for quark confinement supporting the nonAbelian dual
superconductivity using the minimal option for the SU(3) YangMills theory.
This should be compared with the maximal option which is a gauge invarient
version of the Abelian projection in the maximal Abelian gauge: the restricted
field is Abelian $U(1) \times U(1)$ and involves only the Abelian magnetic
monopole, just like the Abelian projection.
In this talk, we focus on discriminating between two reformulations, i.e.,
maximal and minimal options of $SU(3)$ YangMills theory for quark confinement
from the viewpoint of dual superconductivity. For this purpose, we measure the
distribution of the chromoelectric flux connecting a quark and an antiquark and
the induced magneticmonopole current around the flux tube.

As a part of the project studying large $N_f$ QCD, the LatKMI Collaboration
has been investigating the SU(3) gauge theory with four fundamental fermions
(fourflavor QCD). The main purpose of studying fourflavor QCD is to provide a
qualitative comparison to $N_f= 8$, $12$, $16$ QCD; however, a quantitative
comparison to realworld QCD is also interesting. To make such comparisons more
meaningful, it is desirable to use the same kind of lattice action
consistently, so that qualitative difference of different theories are less
affected by artifacts of lattice discretization. Here, we adopt the
highlyimproved staggered quark action with the treelevel Symanzik gauge
action (HISQ/tree), which is exactly the same as the setup for our simulations
for $SU(3)$ gauge theories with $N_f=8$, $12$ and $16$ fundamental
fermions~\cite{Aoki:2013xza, Aoki:2012eq, Aoki:2014oma}. In the next section,
we show the fermion mass dependence of $F_\pi$, $\langle\bar{\psi}\psi\rangle$,
$M_\pi$, $M_\rho$, $M_N$ and their chiral extrapolations. In section 3,
preliminary results of the measurement of the mass of the flavorsinglet scalar
bound state will be reported.

In the preceeding works, we have given a nonAbelian dual superconductivity
picture for quark confinement, and demonstrated the numerical evidences on the
lattice. In this talk, we discuss the confinement and deconfinement phase
transition at finite temperature in view of the dual superconductivity. We
investigate chromomagnetic monopole currents induced by chromoelectric flux in
both confinement and deconfinement phase by the numerical simulations on a
lattice at finite temperature, and discuss the role of the chromomagnetic
monopole in the confinement/deconfinement phase transition.

We give an analytical derivation of the confinement/deconfinement phase
transition at finite temperature in the $SU(N)$ YangMills theory in the
$D$dimensional space time for $D>2$. For this purpose, we use a novel
reformulation of the YangMills theory which allows the gaugeinvariant gluonic
mass term, and calculate analytically the effective potential of the Polyakov
loop average concretely for the $SU(2)$ and $SU(3)$ YangMills theories by
including the gaugeinvariant dynamical gluonic mass $M$. For $D=4$, we give an
estimate on the transition temperature $T_d$ as the ratio $T_d/M$ to the mass
$M$ which has been measured on the lattice at zero temperature and is
calculable also at finite temperature. We show that the order of the phase
transition at $T_d$ is the second order for $SU(2)$ and weakly first order for
$SU(3)$ YangMills theory. We elucidate what is the mechanism for quark
confinement and deconfinement at finite temperature and why the phase
transition occurs at a certain temperature. These initial results are obtained
easily based on the analytical calculations of the "oneloop type" in the first
approximation. We discuss also how these results are improved to eliminate the
artifacts obtained for some thermodynamic observables.

In the search for a composite Higgs boson in walking technicolor models, many
flavor QCD, in particular with $N_f=8$, is an attractive candidate, and has
been found to have a composite flavorsinglet scalar as light as the pion.
Based on lattice simulations of this theory with the HISQ action, we will
present our preliminary results on the scalar decay constant using the
fermionic bilinear operator, and on the mass of the lightest baryon state which
could be a dark matter candidate. Combining these two results, implications for
dark matter direct detection are also discussed.

LatKMI Collaboration discusses the topological insights in manyflavor QCD on
the lattice. We explore walking/conformal/confining phase in $N_\mathrm{f}$ =
4, 8 and 12 (in particular $N_\mathrm{f}$ = 8) lattice QCD via the topological
charge and susceptibility, eigenvalues and anomalous dimension.

Motivated by recent progress on many flavor QCD on a lattice, we investigate
conformal/walking dynamics by using SchwingerDyson (SD) equation within an
improved ladder approximation for twoloop running coupling. By numerically
solving the SD equation, we obtain a pole mass $m_{p}$, pion decay constant
$f_{\pi}$, and investigate the chiral symmetry breaking and mass anomalous
dimension $\gamma_{m}$ in the presence of IR cutoffs $\Lambda_{\mathrm{IR}}$.
We find that the chiral symmetry breaking is suppressed \ if IR cutoff
$\Lambda_{\mathrm{IR}}$ becomes larger than the critical \ value near the
dynamical mass ($\Lambda_{\mathrm{IR}}$ $\simeq m_{D}$) In the conformal phase
the $\gamma_{m}$ is strongly suppressed by IR cutoffs for $\Lambda
_{\mathrm{IR}}$ $\simeq m_{p}$. We, then, obtain finite size hyperscaling (FSS)
relation by adapting a linearized approximation for the SD equation, and
examine the $\gamma_{m}$ The results offer valuable insight and suggestion for
analyses in lattice gauge theories.

First, we give a gaugeindependent definition of chromomagnetic monopoles in
$SU(N)$ YangMills theory which is derived through a nonAbelian Stokes theorem
for the Wilson loop operator. Then we discuss how such magnetic monopoles can
give a nontrivial contribution to the Wilson loop operator for understanding
the area law of the Wilson loop average. Next, we discuss how the magnetic
monopole condensation picture are compatible with the vortex condensation
picture as another promising scenario for quark confinement. We analyze the
profile function of the magnetic flux tube as the nonAbelian vortex solution
of $U(N)$ gaugeHiggs model, which is to be compared with numerical simulations
of the $SU(N)$ YangMills theory on a lattice. This analysis gives an estimate
of the string tension based on the vortex condensation picture, and possible
interactions between two nonAbelian vortices.

The purpose of this paper is to review the recent progress in understanding
quark confinement. The emphasis of this review is placed on how to obtain a
manifestly gaugeindependent picture for quark confinement supporting the dual
superconductivity in the YangMills theory, which should be compared with the
Abelian projection proposed by 't Hooft. The basic tools are novel
reformulations of the YangMills theory based on change of variables extending
the decomposition of the $SU(N)$ YangMills field due to Cho, DuanGe and
FaddeevNiemi, together with the combined use of extended versions of the
DiakonovPetrov version of the nonAbelian Stokes theorem for the $SU(N)$
Wilson loop operator. Moreover, we give the lattice gauge theoretical versions
of the reformulation of the YangMills theory which enables us to perform the
numerical simulations on the lattice. In fact, we present some numerical
evidences for supporting the dual superconductivity for quark confinement. The
numerical simulations include the derivation of the linear potential for static
interquark potential, i.e., nonvanishing string tension, in which the
"Abelian" dominance and magnetic monopole dominance are established,
confirmation of the dual Meissner effect by measuring the chromoelectric flux
tube between quarkantiquark pair, the induced magneticmonopole current, and
the type of dual superconductivity, etc. In addition, we give a direct
connection between the topological configuration of the YangMills field such
as instantons/merons and the magnetic monopole.

The spectrum of twelveflavor QCD has been studied in details by the LatKMI
collaboration. In this proceeding we present our updated results for the
spectrum obtained with the HISQ action at two lattice spacings, several volumes
and fermion masses. In particular, we emphasize the existence of a
flavorsinglet scalar state parametrically light with respect to the rest of
the spectrum, first reported in our paper. This feature is expected to be
present for theories in the conformal window, but the lattice calculation of
such a state is difficult and requires noisereduction techniques together with
large statistics, in order to evaluate disconnected diagrams. Being able to
provide a robust observed connection between a light flavorsinglet scalar and
(near)conformality is an important step towards observing a light composite
Higgs boson in walking technicolor theories on the lattice. We also show
updated results for the mass anomalous dimension $\gamma_m$ obtained from
various spectral quantities, including the string tension, under the assumption
that the theory is inside the conformal window.

We have proposed the nonAbelian dual superconductivity in SU(3) YangMills
theory for the mechanism of quark confinement,and we presented the numerical
evidences in preceding lattice conferences by using the proposed gauge link
decomposition to extract magnetic monopole in the gauge invariant way. In this
talk, we focus on the dual Meissner effects in view of the magnetic monopole in
SU(3) YangMills theory. We measure the chromoelectric and chromomagnetic flux
due to a pair of quark and antiquark source at finite temperature. Then, we
measure the correlation function of Polyakov loops and Polyakov loop average at
various temperatures, and investigate chromomagnetic monopole current induced
by chromomagnetic flux in both confinement and deconfinement phase. We will
discuss the role of the chromoelectric monopole in confinement/deconfinement
phase transition.

We propose the reformulations of the $SU(N)$ YangMills theory toward quark
confinement and mass gap. In fact, we have given a new framework for
reformulating the $SU(N)$ YangMills theory using new field variables. This
includes the preceding works given by Cho, Faddeev and Niemi, as a special case
called the maximal option in our reformulations. The advantage of our
reformulations is that the original nonAbelian gauge field variables can be
changed into the new field variables such that one of them called the
restricted field gives the dominant contribution to quark confinement in the
gaugeindependent way. Our reformulations can be combined with the $SU(N)$
extension of the DiakonovPetrov version of the nonAbelian Stokes theorem for
the Wilson loop operator to give a gaugeinvariant definition for the magnetic
monopole in the $SU(N)$ YangMills theory without the scalar field. In the
socalled minimal option, especially, the restricted field is nonAbelian and
involves the nonAbelian magnetic monopole with the stability group $U(N1)$.
This suggests the nonAbelian dual superconductivity picture for quark
confinement. This should be compared with the maximal option: the restricted
field is Abelian and involves only the Abelian magnetic monopoles with the
stability group $U(1)^{N1}$, just like the Abelian projection. We give some
applications of this reformulation, e.g., the stability for the homogeneous
chromomagnetic condensation of the Savvidy type, the large N treatment for
deriving the dimensional transmutation and understanding the mass gap, and also
the numerical simulations on a lattice which are given by Dr. Shibata in a
subsequent talk.

In oder to investigate quark confinement, we give a new reformulation of the
$SU(N)$ YangMills theory on a lattice and present the results of the numerical
simulations of the $SU(3)$ YangMills theory on a lattice. The numerical
simulations include the derivation of the linear potential for static
interquark potential, i.e., nonvanishing string tension, in which the
"Abelian" dominance and magnetic monopole dominance are established,
confirmation of the dual Meissner effect by measuring the chromoelectric flux
tube between quarkantiquark pair, the induced magneticmonopole current, and
the type of dual superconductivity, etc.

In the $SU(2)$ YangMills theory on the fourdimensional Euclidean lattice,
we confirm the gaugeindependent "Abelian" dominance (or the restricted field
dominance) and gaugeindependent magneticmonopole dominance in the string
tension of the linear potential extracted from the Wilson loop in the
fundamental representation. The dual Meissner effect is observed by
demonstrating the squeezing of the chromoelectric field flux connecting a pair
of quark and antiquark. In addition, the circular magneticmonopole current is
induced around the chromoelectric flux. The type of the dual superconductivity
is also determined by fitting the result with the dual GinzburgLandau model.
Thus the dual superconductor picture for quark confinement is supported in a
gaugeindependent manner. These results are obtained based on a reformulation
of the lattice YangMills theory based on the change of variables a la
ChoDuanGeFaddeevNiemi combined with a nonAbelian Stokes theorem for the
Wilson loop operator. We give a new procedure (called the reduction) for
obtaining the color direction field which plays the central role in this
reformulation.

We present the first observation of a flavorsinglet scalar meson as light as
the pion in $N_f=8$ QCD on the lattice, using the Highly Improved Staggered
Quark action. Such a light scalar meson can be regarded as a composite Higgs
with mass 125 GeV. In accord with our previous lattice results showing that the
theory exhibits walking behavior, the light scalar may be a technidilaton, a
pseudo NambuGoldstone boson of the approximate scale symmetry in walking
technicolor.

The dual superconductivity is a promising mechanism for quark confinement. We
have proposed the nonAbelian dual superconductivity picture for SU(3)
YangMills theory, and showed the restricted field dominance (called
conventionally Abelian dominance), and nonAbelian magnetic monopole dominance
in the string tension. We have further demonstrated by measuring the
chromoelectric flux that the nonAbelian dual Meissner effect exists and
determined that the dual superconductivity for SU(3) case is of type I, which
is in sharp contrast to the SU(2) case: the border of type I and type II.
In this talk, we focus on the confinement/deconfinement phase transition and
the nonAbelian dual superconductivity at a finite temperature: We measure the
Polyakov loop average and correlator and investigate the restricted field
dominance in the Polyakov loop. Then, we measure the chromoelectric flux
between a pair of static quark and antiquark created by a pair of Polyakov
loops, and investigate the nonAbelian dual Meissner effect and its relevance
to the phase transition.

The dual superconductivity is a promising mechanism for quark confinement. We
proposed the nonAbelian dual superconductivity picture for SU(3) YangMills
theory, and demonstrated the restricted field dominance (called conventionally
"Abelian" dominance), and nonAbelian magnetic monopole dominance in the string
tension. In the last conference, we have demonstrated by measuring the
chromoelectric flux that the nonAbelian dual Meissner effect exists and
determined that the dual superconductivity for SU(3) case is of type I, which
is in sharp contrast to the SU(2) case: the border of type I and type II.
In this talk, we focus on the confinement/deconfinemen phase transition and
the nonAbelian dual superconductivity at finite temperature: We measure the
chromoelectric flux between a pair of static quark and antiquark at finite
temperature, and investigate its relevance to the phase transition and the
nonAbelian dual Meissner effect.

We report the calculation of the flavorsinglet scalar in the SU(3) gauge
theory with the degenerate twelve fermions in the fundamental representation
using a HISQtype action at a fixed $\beta$. In order to reduce the large
statistical error coming from the vacuumsubtracted disconnected correlator, we
employ a noise reduction method and a large number of configurations. We
observe that the flavorsinglet scalar is lighter than the pion in this theory
from the calculations with the fermion bilinear and gluonic operators. This
peculiar feature is considered to be due to the infrared conformality of this
theory, and it is a promissing signal for a walking technicolor, where a light
composite Higgs boson is expected to emerge by approximate conformal dynamics.

We measure glueball masses and the string tension in twelveflavour QCD,
aiming at comparing the emerging gluonic spectrum to the mesonic one. When
approaching the critical surface at zero quark mass, the hierarchy of masses in
the different sectors of the spectrum gives a new handle to determine the
existence of an infrared fixed point. We describe the details of our gluonic
measurements and the results obtained on a large number of gauge configurations
generated with the HISQ action. In particular, we focus on the scalar glueball
and its mixing with a flavoursinglet fermionic state, which is lighter than
the pseudoscalar (wouldbe pion) state. The results are interesting in view of
a light composite Higgs boson in walking technicolor theories.