• We highlight the progress, current status, and open challenges of QCD-driven 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 strongly-coupled, complex systems in particle and condensed-matter 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 non-compact $3\hm + 1$--dimensional Abelian gauge field and $2\hm + 1$--dimensional Kogut-Susskind fermions. It was found that defects shift the phase transition point semimetal-insulator towards higher values of a substrate permittivity.
  • We study some of the local CP-odd and magnetic properties of the non-Abelian 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 low-dimensional defects with d=2..3, while a sequence of cooling steps gives rise to four-dimensional instantons and hence a four-dimensional 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 non-Abelian vacuum induced by strong external magnetic field. We perform calculations in the quenched SU(3) lattice gauge theory with tadpole-improved Luscher-Weisz 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 non-zero 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 CP-odd 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 Yang-Mills 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 near-zero 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 (hadron-scale) 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 CP-odd generation of the electric current of quarks directed along the magnetic field]; (4) a CP-odd 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 non-perturbatively 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 low-lying eigenmodes of the kinetic operator for scalar particles, in the color adjoint representation of Yang-Mills 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 low-lying eigenmodes, and find a possible connection to multi-critical matrix models of order m=1.
  • We study correlations between low-lying modes of the overlap Dirac operator and vacuum defects, center vortices and three-dimensional volumes, in lattice SU(2) gluodynamics. The low-lying 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 near-zero 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 non-perturbatively 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 dimension-2 operators both in the MA and the LA gauges. This determination of the coherence length is supported by theoretical and numerical observation that the dimension-2 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 Faddeev-Popov operator are all extended rather than localized (as required for confinement) despite the similarity of the kinetic and Faddeev-Popov operators.
  • As a probe of the Yang-Mills 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 high-lying 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 gauge-Higgs theory. Localization is also absent in the spectrum of the Coulomb gauge Faddeev-Popov operator, as required in Coulomb gauge confinement scenarios.
  • We study the distribution of color electric flux of the three-quark system in quenched and full QCD (with N_f = 2 flavors of dynamical quarks) at zero and finite temperature. To reduce ultra-violet 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 Y--shape 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 Coulomb--like.
  • We present results obtained in QCD with two flavors of non-perturbatively 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 three-dimensional 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 three-dimensional 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 Frohlich-Marchetti 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.