
The main aspects of chiral symmetry in QCD are presented. The necessity of
its spontaneous breakdown is explained. Some lowenergy theorems are reviewed.
The role of chiral effective Lagrangians in the formulation and realization of
chiral perturbation theory is emphasized. The consequences of the presence of
anomalies are sketched.

The motivations of the 1/N expansion method in quantum field theory are
explained. The method is first illustrated with the O(N) model of scalar
fields. A second example is considered with the twodimensional GrossNeveu
model of fermion fields with global U(N) and discrete chiral symmetries. The
case of QCD is briefly sketched.

The main aspects of a gaugeinvariant approach to the description of quark
dynamics in the nonperturbative regime of QCD are first reviewed. In
particular, the role of the parallel transport operation in constructing
gaugeinvariant Green's functions is presented, and the relevance of Wilson
loops for the representation of the interaction is emphasized. Recent
developments, based on the use of polygonal lines for the parallel transport
operation, are then presented. An integrodifferential equation is obtained for
the quark Green's function defined with a phase factor along a single, straight
line segment. It is solved exactly and analytically in the case of
twodimensional QCD in the large $N_c$ limit. The solution displays the
dynamical mass generation phenomenon for quarks, with an infinite number of
branchcut singularities that are stronger than simple poles.

The phenomenon of dynamical quark mass generation is studied in QCD within
the framework of a gauge invariant formalism. An exact relationship is
established between the equation satisfied by the scalar part of the twopoint
gauge invariant quark Green's function and the quarkantiquark bound state
equation in the chiral limit. A possible nontrivial solution of the former
yields a massless pseudoscalar solution of the bound state equation with
vanishing total momentum. The result is also corroborated by the corresponding
WardTakahashi identity. The problem is explicitly solved in twodimensional
QCD in the large$N_c$ limit.

We highlight the progress, current status, and open challenges of QCDdriven
physics, in theory and in experiment. We discuss how the strong interaction is
intimately connected to a broad sweep of physical problems, in settings ranging
from astrophysics and cosmology to stronglycoupled, complex systems in
particle and condensedmatter physics, as well as to searches for physics
beyond the Standard Model. We also discuss how success in describing the strong
interaction impacts other fields, and, in turn, how such subjects can impact
studies of the strong interaction. In the course of the work we offer a
perspective on the many research streams which flow into and out of QCD, as
well as a vision for future developments.

Polygonal lines are used for the paths of the gluon field phase factors
entering in the definition of gauge invariant quark Green's functions. This
allows classification of the Green's functions according to the number of
segments the polygonal lines contain. Functional relations are established
between Green's functions with polygonal lines with different numbers of
segments. An integrodifferential equation is obtained for the quark twopoint
Green's function with a path along a single straight line segment where the
kernels are represented by a series of Wilson loop averages along polygonal
contours. The equation is exactly and analytically solved in the case of
twodimensional QCD in the large$N_c$ limit. The solution displays generation
of an infinite number of dynamical quark masses accompanied with branch point
singularities that are stronger than simple poles. An approximation scheme,
based on the counting of functional derivatives of Wilson loops, is proposed
for the resolution of the equation in four dimensions.

Using gauge invariant quark Green's functions, defined with pathordered
gluon field phase factors along polygonal lines, and functional relations among
them, two compatible bound state equations of the Dirac type are established
for quarkantiquark systems, each relative to the quark or to the antiquark of
the system. The kernels of the bound state equations are defined through a
series of Wilson loop averages along closed polygonal contours and their
functional derivatives on them. A sufficient criterion for spontaneous chiral
symmetry breaking is derived, relating the Goldstone boson wave function in the
zero total momentum limit with the scalar part of the gauge invariant quark
twopoint Green's function.

The properties of the gauge invariant twopoint quark Green's function are
studied in the largeNc limit of twodimensional QCD. The analysis is done by
means of an exact integrodifferential equation. The Green's function is found
infrared finite, with singularities in the momentum squared variable
represented by an infinite number of threshold type branch points with a power
3/2, starting at positive mass squared values, with cuts lying on the positive
real axis. The expression of the Green's function is analytically determined.

The gauge invariant quark Green's function, defined with a pathordered phase
factor along a straight line, is studied in twodimensional QCD in the
largeN_c limit by means of an exact integrodifferential equation. It is found
to be infrared finite with singularities represented by an infinite number of
threshold type branch points with a power of 3/2, starting at positive mass
squared values. The Green's function is analytically determined.

Using an exact integrodifferential equation we study the properties of the
gauge invariant quark Green's function, defined with a pathordered gluon field
phase factor along a straight line, in twodimensional QCD in the largeN_c
limit. The Green's function is found to be infrared finite with singularities
represented by an infinite number of threshold type branch points with a power
equal to 3/2, starting at positive mass squared values. The solution is
analytically determined.

The gauge invariant quark Green's function with a pathordered phase factor
along a straightline is studied in twodimensional QCD in the largeNc limit
by means of an exact integrodifferential equation. Its spectral functions are
analytically determined. They are infrared finite and lie on the positive real
axis of the complex plane of the momentum squared variable, corresponding to
momenta in the forward light cone. Their singularities are represented by an
infinite number of threshold type branch points with powerlaw 3/2, starting
at positive mass values, characterized by an integer number n and increasing
with n. The analytic expression of the Green's function for all momenta is
presented. The appearance of strong threshold singularities is suggestive of
the fact that quarks could not be observed as asymptotic states.

We consider gauge invariant quark twopoint Green's functions in which the
gluonic phase factor follows a skewpolygonal line. Using a particular
representation for the quark propagator in the presence of an external gluon
field, functional relations between Green's functions with different numbers of
segments of the polygonal lines are established. An integral equation is
obtained for the Green's function having a phase factor along a single straight
line. The related kernels involve Wilson loops with skewpolygonal contours and
with functional derivatives along the sides of the contours.

Gauge invariant quark twopoint Green's functions defined with pathordered
gluon field phase factors along skewpolygonal lines joining the quark to the
antiquark are considered. Functional relations between Green's functions with
different numbers of path segments are established. An integral equation is
obtained for the Green's function defined with a phase factor along a single
straight line. The equation implicates an infinite series of twopoint Green's
functions, having an increasing number of path segments; the related kernels
involve Wilson loops with contours corresponding to the skewpolygonal lines of
the accompanying Green's function and with functional derivatives along the
sides of the contours. The series can be viewed as an expansion in terms of the
global number of the functional derivatives of the Wilson loops. The
lowestorder kernel, which involves a Wilson loop with two functional
derivatives, provides the framework for an approximate resolution of the
equation.

The energy and decay width of the piK atom are evaluated in the framework of
the quasipotentialconstraint theory approach. The main electromagnetic and
isospin symmetry breaking corrections to the lowestorder formulas for the
energy shift from the Coulomb binding energy and for the decay width are
calculated. They are estimated to be of the order of a few per cent. We display
formulas to extract the strong interaction Swave piK scattering lengths from
future experimental data concerning the piK atom.

The quarkantiquark gauge invariant Green function is studied through its
dependence on Wilson loops. The latter are saturated, in the large Nc limit and
for large contours, by minimal surfaces. A covariant bound state equation is
derived which in the centerofmass frame and at equaltimes takes the form of
a BreitSalpeter type equation. The largedistance interaction potentials
reduce in the static case to a confining linear vector potential. In general,
the interaction potentials involve contributions having the structure of flux
tube like terms.

Wilson loop averages are evaluated for large contours and in the large N
limit by means of minimal surfaces. This allows the study of the
quarkantiquark gauge invariant Green function through its dependence on Wilson
loops. A covariant bound state equation is derived which in the centerofmass
frame and at equaltimes takes the form of a BreitSalpeter type equation. The
interaction potentials reduce in the static case to a confining linear vector
potential. For moving quarks, flux tube like contributions are present. The
nonrelativistic limit is considered.

Results obtained in recent years in the strange quark sector of chiral
perturbation theory are reviewed and the theoretical relevance of this sector
for probing the phase structure of QCD at zero temperature with respect to the
variation of the number of massless quarks is emphasized.

The largedistance dynamics in quarkonium systems is investigated, in the
large N limit, through the saturation of Wilson loop averages by minimal
surfaces. Using a representation for the quark propagator in the presence of
the external gluon field based on the use of pathordered phase factors, a
covariant threedimensional bound state equation of the BreitSalpeter type is
derived, in which the interaction potentials are provided by the
energymomentum vector of the straight segment joining the quark to the
antiquark and carrying a constant linear energy density, equal to the string
tension. The interaction potentials are confining and reduce to the linear
vector potential in the static case and receive, for moving quarks,
contributions from the moments of inertia of the straight segment. The
selfenergy parts of the quark propagators induce spontaneous breakdown of
chiral symmetry with a mechanism identical to that of the exchange of one
Coulombgluon. The nonrelativistic and ultrarelativistic properties of the
bound state spectrum are studied.

The pionium lifetime is calculated in the framework of the
quasipotentialconstraint theory approach, including the sizable
electromagnetic corrections. The framework of generalized chiral perturbation
theory allows then an analysis of the lifetime value as a function of the pion
Swave scattering lengths with isospins I=0,2, the latter being dependent on
the quark condensate value.

Pionium lifetime corrections to the nonrelativistic formula are calculated in
the framework of the quasipotentialconstraint theory approach. The calculation
extends an earlier evaluation, made in the scheme of standard chiral
perturbation theory, to the scheme of generalized chiral perturbation theory,
in which the quark condensate is left as a free parameter. The pionium lifetime
is calculated as a function of the combination $(a_0^0a_0^2)$ of the $\pi\pi$
Swave scattering lengths with isospin I=0,2.

The relationship between the pionium lifetime and the $\pi\pi$ scattering
lengths is established, including the sizable electromagnetic corrections. The
bound state formalism that is used is that of constraint theory which provides
a covariant threedimensional reduction of the BetheSalpeter equation. The
framework of generalized chiral perturbation theory allows then an analysis of
the lifetime value as a function of the $\pi\pi$ scattering lengths, the latter
being dependent on the quark condensate value.

Pionium lifetime corrections are evaluated in the frameworks of constrained
BetheSalpeter equation and chiral perturbation theory. Corrections of order
$O(\alpha)$ are calculated with respect to the conventional lowestorder
formula, in which the strong interaction amplitude has been calculated to
twoloop order with charged pion masses. The total correction is found to be of
the order of $(6.1\pm 2.0)%$.

The pionium decay width is evaluated in the framework of chiral perturbation
theory and the relativistic bound state formalism of constraint theory.
Corrections of order O(\alpha) are calculated with respect to the conventional
lowestorder formula, in which the strong interaction amplitude has been
calculated to twoloop order with charged pion masses. Strong interaction
corrections from secondorder perturbation theory of the bound state wave
equation are found to be of the order of 0.4%. Electromagnetic radiative
corrections, due to pionphoton interactions, are estimated to be of the order
of 0.1%. Electromagnetic mass shift insertions in internal propagators produce
a correction of the order of 0.3%. The correction due to the passage from the
strong interaction scattering amplitude evaluated with the mass parameter fixed
at the charged pion mass to the amplitude evaluated with the mass parameter
fixed at the neutral pion mass is found to be of the order of 6.4%.

The behavior of the twoparticle Green's function in QED is analyzed in the
limit when one of the particles becomes infinitely massive. It is found that
the dependences of the Green's function on the relative times of the ingoing
and outgoing particles factorize and that the bound state spectrum is the same
as that of the Dirac equation with the static potential created by the heavy
particle. The BetheSalpeter wave function is also determined in terms of the
Dirac wave function. The present result excludes the existence, in the above
limit, of abnormal solutions due to relative time excitations as predicted by
the BetheSalpeter equation in the ladder approximation.

Using connection with quantum field theory, the infinitesimal covariant
abelian gauge transformation laws of relativistic twoparticle constraint
theory wave functions and potentials are established and weak invariance of the
corresponding wave equations shown. Because of the threedimensional projection
operation, these transformation laws are interaction dependent. Simplifications
occur for local potentials, which result, in each formal order of perturbation
theory, from the infrared leading effects of multiphoton exchange diagrams. In
this case, the finite gauge transformation can explicitly be represented, with
a suitable approximation and up to a multiplicative factor, by a momentum
dependent unitary operator that acts in $x$space as a local dilatation
operator. The latter is utilized to reconstruct from the Feynman gauge the
potentials in other linear covariant gauges. The resulting effective potential
of the final PauliSchr\"odinger type eigenvalue equation has the gauge
invariant attractive singularity $\alpha^2/r^2$, leading to a gauge invariant
critical coupling constant $\alpha_c =1/2$.