
We compare higher moments of baryon numbers measured at the RHIC heavy ion
collision experiments with those by the lattice QCD calculations. We employ the
canonical approach, in which we can access the real chemical potential regions
avoiding the sign problem. In the lattice QCD simulations, we study several
fits of the number density in the pure imaginary chemical potential, and
analyze how these fits affects behaviors at the real chemical potential. In the
energy regions between $\sqrt{s}_{NN}$=19.6 and 200 GeV, the susceptibility
calculated at $T/T_c=0.93$ is consistent with experimental data at $0 \le
\mu_B/T < 1.5$, while the kurtosis shows similar behavior with that of the
experimental data in the small $\mu_B/T$ regions $0 \le \mu_B/T < 0.3$. The
experimental data at $\sqrt{s}_{NN}=$ 11.5 shows quite different behavior. The
lattice result in the deconfinement region,$T/T_c=1.35$, is far from
experimental data.

We study numerically the chromoelectricchromomagnetic asymmetry of the
dimension two gluon condensate as well as the longitudinal gluon propagator at
$T\simeq T_c$ in the Landaugauge $SU(2)$ lattice gauge theory. We show that
substantial correlation between the asymmetry and the Polyakov loop as well as
the correlation between the longitudinal propagator and the Polyakov loop pave
the way to studies of the critical behavior of the asymmetry and the
longitudinal propagator. The respective values of critical exponents and
amplitudes are evaluated.

Using GPGPU techniques and multiprecision calculation we developed the code
to study QCD phase transition line in the canonical approach. The canonical
approach is a powerful tool to investigate sign problem in Lattice QCD. The
central part of the canonical approach is the fugacity expansion of the grand
canonical partition functions. Canonical partition functions $Z_n(T)$ are
coefficients of this expansion. Using various methods we study properties of
$Z_n(T)$. At the last step we perform cubic spline for temperature dependence
of $Z_n(T)$ at fixed $n$ and compute baryon number susceptibility $\chi_B/T^2$
as function of temperature. After that we compute numerically $\partial\chi/
\partial T$ and restore crossover line in QCD phase diagram. We use improved
Wilson fermions and Iwasaki gauge action on the $16^3 \times 4$ lattice with
$m_{\pi}/m_{\rho} = 0.8$ as a sandbox to check the canonical approach. In this
framework we obtain coefficient in parametrization of crossover line
$T_c(\mu_B^2)=T_c\left(c\kappa\, \mu_B^2/T_c^2\right)$ with $\kappa = 0.0453
\pm 0.0099$.

We study the topological structure of $SU(3)$ lattice gluodynamics by cluster
analysis. This methodological study is meant as preparation for full QCD. The
topological charge density is becoming visible in the process of overimproved
gradient flow, which is monitored by means of the the Inverse Participation
Ratio (IPR). The flow is stopped at the moment when calorons dissociate into
dyons due to the overimproved character of the underlying action. This gives
the possibility to simultaneously detect all three dyonic constituents of KvBLL
calorons in the gluonic field. The behaviour of the average Polyakov loop under
(overimproved) gradient flow could be also (as its value) a diagnostics for the
actual phase the configuration is belonging to. Timelike Abelian monopole
currents and specific patterns of the local Polyakov loop are correlated with
the topological clusters.The spectrum of reconstructed cluster charges $Q_{cl}$
corresponds to the phases. It is scattered around $Q_{cl} \approx \pm 1/3$ in
the confined phase, whereas it is $Q_{cl} \approx \pm 0.5 \div 0.7$ for heavy
dyons and $Q_{cl} < 0.3$ for light dyons in the deconfined phase. Heavy dyons
are increasingly suppressed with increasing temperature. The paper is dedicated
to the memory of Michael MuellerPreussker who was a member of our research
group for more than twenty years.

The SU(3) flavour symmetry breaking expansion in up, down and strange quark
masses is extended from hadron masses to meson decay constants. This allows a
determination of the ratio of kaon to pion decay constants in QCD. Furthermore
when using partially quenched valence quarks the expansion is such that SU(2)
isospin breaking effects can also be determined. It is found that the lowest
order SU(3) flavour symmetry breaking expansion (or GellMannOkubo expansion)
works very well. Simulations are performed for 2+1 flavours of clover fermions
at four lattice spacings.

We discuss scale setting in the context of 2+1 dynamical fermion simulations
where we approach the physical point in the quark mass plane keeping the
average quark mass constant. We have simulations at four beta values, and after
determining the paths and lattice spacings, we give an estimation of the
phenomenological values of various Wilson flow scales.

We give a determination of the phenomenological value of the Wilson (or
gradient) flow scales t0 and w0 for 2+1 flavours of dynamical quarks. The
simulations are performed keeping the average quark mass constant, which allows
the approach to the physical point to be made in a controlled manner. O(a)
improved clover fermions are used and together with four lattice spacings this
allows the continuum extrapolation to be taken.

The topological structure of lattice gluodynamics is studied at intermediate
resolution scale in the deconfining phase with the help of a cluster analysis.
UV filtered topological charge densities are determined from a fixed number of
lowlying eigenmodes of the overlap Dirac operator with three types of temporal
boundary conditions applied to the valence quark fields. This method usually
allows to find all three distinguished (anti)dyon constituents in the gauge
field of Kraanvan BaalLeeLu (anti)caloron solutions. The clustering of the
three topological charge densities in Monte Carlo generated configurations is
then used to mark the positions of anticipated (anti)dyons of the corresponding
type. In order to support this interpretation, inside these clusters, we search
also for timelike Abelian monopole currents (defined in the maximally Abelian
gauge) as well as for local holonomies with at least two approximately
degenerated eigenvalues. Our results support the view that light dyonantidyon
pairs  in contrast to the heavy (anti)caloron dyon constituents  contribute
dominantly to thermal YangMills fields in the deconfinement phase. This paper
is dedicated to the memory of Pierre van Baal and Dmitri Igorevich Diakonov who
have influenced our work very much.

The properties of the thermal Abelian colormagnetic monopoles in the
maximally Abelian gauge are studied in the vicinity of the
confinementdeconfinement phase transition in the lattice $SU(3)$ gluodynamics
and lattice QCD. We compute the density and interaction parameters of the
thermal monopoles. We find that the density of the thermal monopoles $\rho(T)$
jumps up near the transition temperature $T_c$. Additionally we present new
results on the percolation transition in $SU(3)$ gluodynamics which is known to
coincide in gluodynamics with the confinementdeconfinement phase transition.

Gluodynamics and twoflavor QCD at nonzero temperature are studied with the
socalled overimproved cooling technique under which caloron solutions may
remain stable. We consider topological configurations either at the first
occuring stable plateau of topological charge or at the first (anti)selfdual
plateau and find the corresponding topological susceptibility at various
temperatures on both sides of the thermal transition or crossover. In pure
gluodynamics the topological susceptibility drops sharply at the deconfinement
temperature while in full QCD it decreases smoothly at temperatures above the
pseudocritical one. The results are close to those calculated by other methods.
We interpret our findings in terms of the (in)stability of calorons with
nontrivial holonomy and their dyon constituents against overimproved cooling.

The critical temperature is computed for N_f = 2 dynamical flavors of
nonperturbatively improved Wilson fermions. The new simulations are performed
on lattices 40^3 x 14 with lattice spacing and pion mass about 0.08 fm and 200
MeV, respectively. We find the deconfinement and chiral phase transitions to
coincide within numerical precision. 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 T_c = 172(3)(6) MeV.

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 present clear evidence of strong effects of Gribov copies in Landau gauge
gluon and ghost propagators computed on the lattice at small momenta by
employing a new approach to Landau gauge fixing and a more effective numerical
algorithm. It is further shown that the new approach substantially decreases
notorious finitevolume effects.

The spatial string tension across a crossover from the low temperature phase
to the high temperature phase is computed in QCD with two flavors of
nonperturbatively improved Wilson fermions at small lattice spacing a \sim
0.12fm. We find that in the low temperature phase spatial string tension agrees
well with zero temperature string tension. Furthermore, it does not show
increasing for temperatures up to T = 1.36 T_{pc}, the highest temperature
considered. Our results agree with some theoretical predictions.

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.

The SU(2) gluon propagator in Landau gauge is studied on the lattice. Our
gauge fixing procedure employs simulated annealing and Z(2)flips. It finds
higher maxima of the gauge functional compared with those obtained with the
standard overrelaxation and leads to systematic deviations of the gluon
propagator in the infrared region. In particular, finitesize effects for
lattice sizes from (1.7 fm)^4 up to (6.5 fm)^4 become weak. The propagator
shows a plateau at p \approx 300 MeV.

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 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.

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.

We study QCD with two flavors of nonperturbatively improved Wilson fermions
at finite temperature on the $16^3 8$ lattice. We determine the transition
temperature at lattice spacings as small as $a \sim 0.12$ fm, and study string
breaking below the finite temperature transition. We find that the static
potential can be fitted by a twostate ansatz, including a string state and a
twomeson state. We investigate the role of Abelian monopoles at finite
temperature.

We study gauge dependence of the recently suggested definition of the singlet
and adjoint potentials in SU(2) lattice gauge theory. We find that in the (time
local) maximal tree axial gauge the singlet potential obtained from the gauge
dependent correlator ${Tr} L(x)L^\dagger(y)$ differs from that computed in the
Coulomb gauge. In the generalized Coulomb gauge we find the range of the
parameter values in which the singlet potential differs from that in the
Coulomb gauge.

We investigate the confining properties of the QCD vacuum with $N_f=2$
flavors of dynamical quarks, and compare the results with the properties of the
quenched theory. We use nonperturbatively $\mathcal{O}(a)$ improved Wilson
fermions to keep cutoff effects small. We focus on color magnetic monopoles.
Among the quantities we study are the monopole density and the monopole
screening length, the static potential and the profile of the color electric
flux tube. We furthermore derive the lowenergy effective monopole action.
Marked differences between the quenched and dynamical vacuum are found.

We study the action and the energy densities of the confining string in the
indirect Z(2) projection of SU(2) lattice gauge theory. We find that the width
of the confining string is proportional to the logarithm of the distance
between the quark and antiquark. In $2q2\bar{q}$ system we observe the effect
of the reconstruction of the flux tube when we change the distance between two
$q\bar{q}$ pairs.

We present results of measurements of various geometrical characteristics of
the monopole clusters in the maximally Abelian projection of SU(2) lattice
gauge theory. We observe scaling for the observables tested. Short clusters
correspond to random walks at small scale but have longrange correlations at
the hadronic scale. On the other hand, the percolating cluster at the hadronic
scale does not correspond to a random walk.