
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.

The results of the simulations by Monte Carlo method of graphene with
structural defects are presented.
The calculations are performed within an effective quantum field theory with
noncompact $3\hm + 1$dimensional Abelian gauge field and $2\hm +
1$dimensional KogutSusskind fermions. It was found that defects shift the
phase transition point semimetalinsulator towards higher values of a substrate
permittivity.

We study some of the local CPodd and magnetic properties of the nonAbelian
vacuum with use of overlap fermions within the quenched lattice gauge theory.
Among these properties are the following: inhomogeneous spatial distribution of
the topological charge density (chirality for massless fermions) in SU(2)
gluodynamics (for uncooled gauge configurations the chirality is localized on
lowdimensional defects with d=2..3, while a sequence of cooling steps gives
rise to fourdimensional instantons and hence a fourdimensional structure of
the chirality distribution); finite local fluctuations of the chirality growing
with the strength of an external magnetic field; magnetization and
susceptibility of the QCD vacuum in SU(3) theory; magnetic catalysis of the
chiral symmetry breaking, and the electric conductivity of the QCD vacuum in
strong magnetic fields.

We study some properties of the nonAbelian vacuum induced by strong external
magnetic field. We perform calculations in the quenched SU(3) lattice gauge
theory with tadpoleimproved LuscherWeisz action and chirally invariant
lattice Dirac operator. The following results are obtained: The chiral symmetry
breaking is enhanced by the magnetic field. The chiral condensate depends on
the strength of the applied field as a power function with exponent n = 1.6 +/
0.2. There is a paramagnetic polarization of the vacuum. The corresponding
susceptibility and other magnetic properties are calculated and compared with
the theoretical estimations. There are nonzero local fluctuations of the
chirality and electromagnetic current, which grow with the magnetic field
strength. These fluctuations can be a manifestation of the Chiral Magnetic
Effect (CME).

We study the correlator of two vector currents in quenched $SU\lr{2}$ lattice
gauge theory with a chirally invariant lattice Dirac operator with a constant
external magnetic field. It is found that in the confinement phase the
correlator of the components of the current parallel to the magnetic field
decays much slower than in the absence of a magnetic field, while for other
components the correlation length slightly decreases. We apply the maximal
entropy method to extract the corresponding spectral function. In the limit of
zero frequency this spectral function yields the electric conductivity of the
quenched theory. We find that in the confinement phase the external magnetic
field induces nonzero electric conductivity along the direction of the field,
transforming the system from an insulator into an anisotropic conductor. In the
deconfinement phase the conductivity does not exhibit any sizable dependence on
the magnetic field.

The critical temperature and the nature of the QCD finite temperature phase
transition are determined for Nf=2 dynamical flavors of nonperturbatively
improved Wilson fermions. The calculations are performed on large lattices with
temporal extents Nt=12, 10 and 8, and lattice spacings down to a = 0.075 fm. We
find the deconfinement and chiral phasetransitions to take place at the same
temperature. Our results are in broad agreement with a second order phase
transition in the chiral limit. The critical temperature at the physical quark
mass is found to be Tc = 174(3)(6) MeV.

We show numerically that quarks develop an electric dipole moment in the
direction of a sufficiently intense magnetic field due to local fluctuations of
topological charge. This anomalous CPodd effect is a spin analogue of the
Chiral Magnetic Effect in QCD.

A possible experimental observation of the chiral magnetic effect in heavy
ion collisions at RHIC was recently reported by the STAR Collaboration. We
study signatures of this effect in SU(2) lattice gluodynamics with the chirally
invariant Dirac operator. We find that at zero temperature the local
fluctuations of an electric current of quarks and chirality fluctuations
increase with external Abelian magnetic field. The external magnetic field
leads to spatial separation of the quark's electric charges. The separation
increases with the strength of the magnetic field. As temperature gets higher
the dependence of these quantities on the strength of the magnetic field
becomes weaker. In the deconfinement phase the local fluctuations of the chiral
density and of the spatial components of the quarks electric current are large
and are almost independent on the external magnetic field. The local
fluctuations of the electric charge density decrease with the strength of the
magnetic field in this phase.

The chiral magnetization properties of cold and hot vacua are studied using
quenched simulations in lattice YangMills theory. In weak external magnetic
fields the magnetization is proportional to the first power of the magnetic
field. We evaluate numerically the coefficient of the proportionality (the
chiral susceptibility) using nearzero eigenmodes of overlap fermions. We found
that the product of the chiral susceptibility and the chiral condensate equals
to 46(3) MeV. This value is very close to the phenomenological value of 50 MeV.
In strong fields the magnetization is a nonlinear function of the applied
magnetic field. We find that the nonlinear features of the magnetization are
well described by an inverse tangent function. The magnetization is weakly
sensitive to temperature in the confinement phase.

Vacuum of Quantum Chromodynamics in very strong (hadronscale) magnetic
fields exhibits many interesting nonperturbative effects. Some of these effects
can be studied with the help of lattice simulations in quenched QCD. We review
our recent results demonstrating that very strong external magnetic fields lead
to (1) the enhancement of the chiral symmetry breaking [the quark condensate
rises with the increase of the external magnetic field]; (2) the chiral
magnetization of the QCD vacuum [spins of the quarks turn parallel to the
external field]; (3) the chiral magnetic effect [a CPodd generation of the
electric current of quarks directed along the magnetic field]; (4) a CPodd
generation of the anomalous quark electric dipole moment along the axis of
magnetic field. The first three effects were already predicted theoretically,
and subsequently observed numerically in our simulations, while the fourth
effect is a new result.

The chiral magnetic effect is the generation of electric current of quarks
along external magnetic field in the background of topologically nontrivial
gluon fields. There is a recent evidence that this effect is observed by the
STAR Collaboration in heavy ion collisions at RHIC. In our paper we study
qualitative signatures of the chiral magnetic effect using quenched lattice
simulations. We find indications that the electric current is indeed enhanced
in the direction of the magnetic field both in equilibrium configurations of
the quantum gluon fields and in a smooth gluon background with nonzero
topological charge. In the confinement phase the magnetic field enhances the
local fluctuations of both the electric charge and chiral charge densities. In
the deconfinement phase the effects of the magnetic field become smaller,
possibly due to thermal screening. Using a simple model of a fireball we obtain
a good agreement between our data and experimental results of the STAR
Collaboration.

We present results of finite temperature QCD from N_f=2 flavors of
nonperturbatively improved Wilson fermions on lattices with N_t=8, 10 and 12.
The transition temperature T_c is determined for the first time at lattice
spacings as low as a = 0.08 fm, albeit at pion masses 1.3 < r_0 m_\pi. We
furthermore compute the screening masses in various color channels at T>T_c.

We study the spectrum of lowlying eigenmodes of the kinetic operator for
scalar particles, in the color adjoint representation of YangMills theory. The
kinetic operator is the covariant Laplacian, plus a constant which serves to
renormalize mass. In the pure gauge theory, our data indicates that the
interval between the lowest eigenvalue and the mobility edge tends to infinity
in the continuum limit. On these grounds, it is suggested that the perturbative
expression for the scalar propagator may be misleading even at distance scales
that are small compared to the confinement scale. We also measure the density
of lowlying eigenmodes, and find a possible connection to multicritical
matrix models of order m=1.

We study correlations between lowlying modes of the overlap Dirac operator
and vacuum defects, center vortices and threedimensional volumes, in lattice
SU(2) gluodynamics. The lowlying modes are apparently sensitive to topological
properties of the underlying gluon field configurations while the vacuum
defects are crucial for the confinement. We find distinct positive correlation
in both cases. In case of vortices the correlation is stronger.

After fixing the Maximal Abelian gauge in SU(2) lattice gauge theory we
decompose the nonabelian gauge field into the so called monopole field and the
modified nonabelian field with monopoles removed. We then calculate respective
static potentials and find that the potential due to the modified nonabelian
field is nonconfining while, as is well known, the monopole field potential is
linear. Furthermore, we show that the sum of these potentials approximates the
nonabelian static potential with 5% or higher precision at all distances
considered. We conclude that at large distances the monopole field potential
describes the classical energy of the hadronic string while the modified
nonabelian field potential describes the string fluctuations. Similar
decomposition was observed to work for the adjoint static potential. A check
was also made of the center projection in the direct center gauge. Two static
potentials, determined by projected $Z_2$ and by modified nonabelian field
without $Z_2$ component were calculated. It was found that their sum is a
substantially worse approximation of the SU(2) static potential than that found
in the monopole case. It is further demonstrated that similar decomposition can
be made for the flux tube action/energy density.

We consider properties of zero and nearzero modes for overlap fermion
operator in SU(2) lattice gluodynamics. The density of the states is of the
order of Lambda(QCD) while the localization volume of the modes tends to zero
in physical units with the lattice spacing tending to zero. The situation
changes drastically when we study "vortex removed" configurations.

We review briefly recent results of lattice simulations on 3d domains in the
vacuum state of SU(2) gluodynamics. The defects are defined as unification of
all the negative links in central projection under condition that the total
number of negative links is minimized. In the continuum limit, negative links
correspond, generally speaking to singular fields. The data indicate that total
volume of the defects scales in physical units. We consider also correlator of
negative links. The correlator scales in physical units as well, within the
error bars. A new observation reported here is a strong anisotropy of the
correlator.

We present results obtained in QCD with two flavors of nonperturbatively
improved Wilson fermions at finite temperature on $16^3 \times 8$ and $24^3
\times 10$ lattices. We determine the transition temperature in the range of
quark masses $0.6<m_\pi/m_\rho<0.8$ at lattice spacing a$\approx$0.1 fm and
extrapolate the transition temperature to the continuum and to the chiral
limits. We also discuss the order of phase transition.

The type of the vacuum is studied numerically in the maximally Abelian (MA)
gauge and in the Landau (LA) gauge of SU(2) gluodynamics. The type of the
vacuum is determined by a ratio between the dual coherence and the dual
penetration lengths. The dual penetration length is determined from
correlations between Wilson loops and electric fields in both gauges. The dual
coherence length is found from correlations between Wilson loops and
dimension2 operators both in the MA and the LA gauges. This determination of
the coherence length is supported by theoretical and numerical observation that
the dimension2 gluon operators in the studied gauges have a strong correlation
with the monopole current determined in the MA gauge. We find numerically that
the dual penetration lengths and the dual coherence lengths in the LA and the
MA gauges are almost the same. Therefore we conclude, that in both gauges the
type of the vacuum in the confinement phase is near to the border between the
type 1 and the type 2 dual superconductors.

We study numerically the eigenmode spectrum of the covariant lattice
Laplacian, in the fundamental SU(2) color group representation. It is found
that eigenmodes at the lower and upper ends of the spectrum are localized, and
that the localization volume scales. In contrast, the eigenmodes of the lattice
FaddeevPopov operator are all extended rather than localized (as required for
confinement) despite the similarity of the kinetic and FaddeevPopov operators.

As a probe of the YangMills vacuum, we study numerically the eigenmode
spectrum of the covariant lattice Laplacian operator. We find that the
eigenmodes at the low and high ends of the spectrum are localized in finite
regions whose volume is insensitive to the lattice volume. We also find that
the vacuum is seen very differently by localized modes of the covariant
Laplacian in different representations of the gauge group. In the fundamental
representation, the data suggests that the localization volume is finite in
physical units set by the string tension, and localization disappears when
center vortices are removed. In the adjoint and j=3/2 representations the low
and highlying modes are far more localized, and the localization volume
appears to scale to zero, in physical units, in the continuum limit. The
adjoint Laplacian is insensitive to vortex removal, but we find that
exponential localization is absent for adjoint eigenmodes in the Higgs phase of
a gaugeHiggs theory. Localization is also absent in the spectrum of the
Coulomb gauge FaddeevPopov operator, as required in Coulomb gauge confinement
scenarios.

We study the distribution of color electric flux of the threequark system in
quenched and full QCD (with N_f = 2 flavors of dynamical quarks) at zero and
finite temperature. To reduce ultraviolet fluctuations, the calculations are
done in the abelian projected theory fixed to the maximally abelian gauge. In
the confined phase we find clear evidence for a Yshape flux tube surrounded
and formed by the solenoidal monopole current, in accordance with the dual
superconductor picture of confinement. In the deconfined, high temperature
phase monopoles cease to condense, and the distribution of the color electric
field becomes Coulomblike.

We present results obtained in QCD with two flavors of nonperturbatively
improved Wilson fermions at finite temperature on $16^3 \times 8$ and $24^3
\times 10$ lattices. We determine the transition temperature in the range of
quark masses $0.6<m_\pi/m_\rho<0.8$ at lattice spacing a$\approx$0.1 fm and
extrapolate the transition temperature to the continuum and to the chiral
limits.

Performing lattice simulations of the four dimensional SU(2)gluodynamics we
find evidence for existence of threedimensional domains whose total volume
scales in physical units. Technically, the domains are defined in terms of the
minimal density of negative links in Z(2) projection of gauge fields. The
volume can be viewed also as the minimal volume bound by the center vortices.
We argue that the threedimensional domains are closely related to confinement.

We study the temperature dependence of the monopole condensate in different
Abelian projections of the SU(2) lattice gauge theory. Using the
FrohlichMarchetti monopole creation operator we show numerically that the
monopole condensate depends on the choice of the Abelian projection. Contrary
to the claims in the current literature we observe that in the Abelian Polyakov
gauge and in the field strength gauge the monopole condensate does not vanish
at the critical temperature and thus is not an order parameter.