• 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 chromoelectric-chromomagnetic asymmetry of the dimension two gluon condensate as well as the longitudinal gluon propagator at $T\simeq T_c$ in the Landau-gauge $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 multi-precision 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 Mueller-Preussker 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 Gell-Mann-Okubo 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 low-lying 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 Kraan-van Baal-Lee-Lu (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 time-like 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 dyon-antidyon pairs - in contrast to the heavy (anti)caloron dyon constituents - contribute dominantly to thermal Yang-Mills 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 color-magnetic monopoles in the maximally Abelian gauge are studied in the vicinity of the confinement-deconfinement 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 confinement-deconfinement phase transition.
  • Gluodynamics and two-flavor QCD at non-zero temperature are studied with the so-called 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 non-trivial 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 finite-volume 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 non-perturbatively 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 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.
  • 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, finite-size 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 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.
  • 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.
  • We study QCD with two flavors of non-perturbatively 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 two-state ansatz, including a string state and a two-meson 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 non-perturbatively $\mathcal{O}(a)$ improved Wilson fermions to keep cut-off 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 low-energy 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 long-range correlations at the hadronic scale. On the other hand, the percolating cluster at the hadronic scale does not correspond to a random walk.