
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.

In this paper we study the confinement/deconfinement transition in lattice
$SU(2)$ QCD at finite quark density and zero temperature. The simulations are
performed on an $32^4$ lattice with rooted staggered fermions at a lattice
spacing $a = 0.044 \mathrm{~fm}$. This small lattice spacing allowed us to
reach very large baryon density (up to quark chemical potential $\mu_q > 2000
\mathrm{~MeV}$) avoiding strong lattice artifacts. In the region $\mu_q\sim
1000 \mathrm{~MeV}$ we observe for the first time the confinement/deconfinement
transition which manifests itself in rising of the Polyakov loop and vanishing
of the string tension $\sigma$. After the deconfinement is achieved at $\mu_q >
1000 \mathrm{~MeV}$, we observe a monotonous decrease of the spatial string
tension $\sigma_s$ which ends up with $\sigma_s$ vanishing at $\mu_q > 2000
\mathrm{~MeV}$. From this observation we draw the conclusion that the
confinement/deconfinement transition at finite density and zero temperature is
quite different from that at finite temperature and zero density. Our results
indicate that in very dense matter the quarkgluon plasma is in essence a
weakly interacting gas of quarks and gluons without a magnetic screening mass
in the system, sharply different from a quarkgluon plasma at large
temperature.

I report on the first application of a novel, generalized Bayesian
reconstruction (BR) method for spectral functions to the characterization of
QCD constituents. These spectral functions find applications in offshell
kinetics of the quarkgluon plasma and in calculations of transport
coefficients. The new BR method is applied to Euclidean propagator data,
obtained in Landau gauge on lattices with $N_f=2+1+1$ dynamical flavors by the
"twisted mass at finite temperature" (tmfT) collaboration. The deployed
reconstruction method is designed for spectral functions that can exhibit
positivity violation (opposed to that of hadronic bound states). The
transversal and longitudinal gluon spectral functions show a robust structure
composed of quasiparticle peak and a negative trough. Characteristic
differences between the hadronic and the plasma phase and between the two
channels become visible. We obtain the temperature dependence of the
transversal and longitudinal gluon masses.

Substantial experimental and theoretical efforts worldwide are devoted to
explore the phase diagram of strongly interacting matter. At LHC and top RHIC
energies, QCD matter is studied at very high temperatures and nearly vanishing
netbaryon densities. There is evidence that a QuarkGluonPlasma (QGP) was
created at experiments at RHIC and LHC. The transition from the QGP back to the
hadron gas is found to be a smooth cross over. For larger netbaryon densities
and lower temperatures, it is expected that the QCD phase diagram exhibits a
rich structure, such as a firstorder phase transition between hadronic and
partonic matter which terminates in a critical point, or exotic phases like
quarkyonic matter. The discovery of these landmarks would be a breakthrough in
our understanding of the strong interaction and is therefore in the focus of
various highenergy heavyion research programs. The Compressed Baryonic Matter
(CBM) experiment at FAIR will play a unique role in the exploration of the QCD
phase diagram in the region of high netbaryon densities, because it is
designed to run at unprecedented interaction rates. Highrate operation is the
key prerequisite for highprecision measurements of multidifferential
observables and of rare diagnostic probes which are sensitive to the dense
phase of the nuclear fireball. The goal of the CBM experiment at SIS100
(sqrt(s_NN) = 2.7  4.9 GeV) is to discover fundamental properties of QCD
matter: the phase structure at large baryonchemical potentials (mu_B > 500
MeV), effects of chiral symmetry, and the equationofstate at high density as
it is expected to occur in the core of neutron stars. In this article, we
review the motivation for and the physics programme of CBM, including
activities before the start of data taking in 2022, in the context of the
worldwide efforts to explore highdensity QCD matter.

We apply the liquid droplet model to describe the clustering phenomenon in
SU(2) gluodynamics, especially, in the vicinity of the deconfinement phase
transition. In particular, we analyze the size distributions of clusters formed
by the Polyakov loops of the same sign. Within such an approach this phase
transition can be considered as the transition between two types of liquids
where one of the liquids (the largest droplet of a certain Polyakov loop sign)
experiences a condensation, while the other one (the next to largest droplet of
opposite Polyakov loop sign) evaporates. The clusters of smaller sizes form two
accompanying gases, and their size distributions are described by the liquid
droplet parameterization. By fitting the lattice data we have extracted the
value of Fisher exponent $\tau =$ 1.806 $\pm$ 0.008. Also we found that the
temperature dependences of the surface tension of both gaseous clusters are
entirely different below and above the phase transition and, hence, they can
serve as an order parameter. The critical exponents of the surface tension
coefficient in the vicinity of the phase transition are found. Our analysis
shows that the temperature dependence of the surface tension coefficient above
the critical temperature has a $T^2$ behavior in one gas of clusters and $T^4$
in the other one.

The liquid droplet formula is applied to an analysis of the properties of
geometrical (anti)clusters formed in SU(2) gluodynamics by the Polyakov loops
of the same sign. Using this approach, we explain the phase transition in SU(2)
gluodynamics as a transition between two liquids during which one of the liquid
droplets (the largest cluster of a certain Polyakov loop sign) experiences a
condensation, while the droplet of another liquid (the next to the largest
cluster of the opposite sign of Polyakov loop) evaporates. The clusters of
smaller sizes form two accompanying gases, which behave oppositely to their
liquids. The liquid droplet formula is used to analyze the size distributions
of the gaseous (anti)clusters. The fit of these distributions allows us to
extract the temperature dependence of surface tension and the value of Fisher
topological exponent $\tau$ for both kinds of gaseous clusters. It is shown
that the surface tension coeficient of gaseous (anti)clusters can serve as an
order parameter of the deconfinement phase transition in SU(2) gluodynamics.
The Fisher topological exponent $\tau$ of clusters and anticlusters is found to
have the same value $1.806 \pm 0.008$. This value disagrees with the famous
Fisher droplet model, but it agrees well with an exactly solvable model of the
nuclear liquidgas phase transition. This finding may evidence for the fact
that the SU(2) gluodynamics and this exactly solvable model of nuclear
liquidgas phase transition are in the same universality class.

In this paper we carry out a lowtemperature scan of the phase diagram of
dense twocolor QCD with $N_f=2$ quarks. The study is conducted using lattice
simulation with rooted staggered quarks. At small chemical potential we observe
the hadronic phase, where the theory is in a confining state, chiral symmetry
is broken, the baryon density is zero and there is no diquark condensate. At
the critical point $\mu = m_{\pi}/2$ we observe the expected second order
transition to BoseEinstein condensation of scalar diquarks. In this phase the
system is still in confinement in conjunction with nonzero baryon density, but
the chiral symmetry is restored in the chiral limit. We have also found that in
the first two phases the system is well described by chiral perturbation
theory. For larger values of the chemical potential the system turns into
another phase, where the relevant degrees of freedom are fermions residing
inside the Fermi sphere, and the diquark condensation takes place on the Fermi
surface. In this phase the system is still in confinement, chiral symmetry is
restored and the system is very similar to the quarkyonic state predicted by
SU($N_c$) theory at large $N_c$.

We study lattice QCD with $N_f=2$ Wilson fermions at nonzero imaginary
chemical potential and nonzero temperature. We relate the Roberge  Weiss phase
transition to the properties of dyons which are constituents of the KvBLL
calorons. We present numerical evidence that the characteristic features of the
spectral gap of the overlap Dirac operator as function of an angle modifying
the boundary condition are determined by the $Z_3$ sector of the respective
imaginary chemical potential. We then demonstrate that dyon excitations in
thermal configurations could be responsible (in line with perturbative
excitations) for these phenomena.

We study the topological structure of QCD by cluster analysis. The fermionic
topological charge density is constructed from low lying modes of overlap Dirac
operator for three types of temporal boundary conditions for fermion field.
This gives the possibility to mark all three dyon constitutents of KvBLL
caloron in gluonic fields. The gluonic topological charge density is appearing
in the process of overimproved gradient flow process stopped at the moment when
it maximally matches the fermionic topological charge density. This corresponds
to the smearing of gluonic fields up to the scale set by dyon size. The
timelike Abelian monopoles and specific KvBLL pattern of Polyakov line
correlate with topological clusters.

In this paper we report on lattice simulations of SU(3)QCD with nonzero
chiral chemical potential. We focus on the influence of the chiral chemical
potential on the confinement/deconfinement phase transition and the
breaking/restoration of chiral symmetry. The simulation is carried out with
dynamical Wilson fermions. We find that the critical temperature rises as the
chiral chemical potential grows.

Here we present several remarkable irregularities at chemical freezeout
which are found using an advanced version of the hadron resonance gas model.
The most prominent of them are the sharp peak of the trace anomaly existing at
chemical freezeout at the center of mass energy 4.9 GeV and two sets of highly
correlated quasiplateaus in the collision energy dependence of the entropy per
baryon, total pion number per baryon, and thermal pion number per baryon which
we found at the center of mass energies 3.84.9 GeV and 7.610 GeV. The low
energy set of quasiplateaus was predicted a long time ago. On the basis of the
generalized shockadiabat model we demonstrate that the low energy correlated
quasiplateaus give evidence for the anomalous thermodynamic properties inside
the quarkgluonhadron mixed phase. It is also shown that the trace anomaly
sharp peak at chemical freezeout corresponds to the trace anomaly peak at the
boundary between the mixed phase and quark gluon plasma. We argue that the high
energy correlated quasiplateaus may correspond to a second phase transition
and discuss its possible origin and location. Besides we suggest two new
observables which may serve as clear signals of these phase transformations.

We present an ongoing project aimed at determining the thermodynamic Equation
of State (EoS) of quarkgluon matter from lattice QCD with two generations of
dynamical quarks. We employ the Wilson twisted mass implementation for the
fermionic fields and the improved Iwasaki gauge action. Relying on $T=0$ data
obtained by the ETM Collaboration the strange and charm quark masses are fixed
at their physical values, while the pion mass takes four values in the range
from 470 MeV down to 210 MeV. The temperature is varied within a fixedlattice
scale approach. The values for the pseudocritical temperature are obtained from
various observables. For the EoS we show preliminary results for the pure
gluonic contribution obtained at the pion mass value 370 MeV, where we can
compare with previously obtained results with $N_f=2$ degenerate light
flavours.

We present results for the topological susceptibility at nonzero temperature
obtained from lattice QCD with four dynamical quark flavours. We apply
different smoothing methods, including gradient Wilson flow and overimproved
cooling, before calculating the susceptibility. It is shown that the considered
smoothing techniques basically agree among each other, and that there are
simple scaling relations between flow time and the number of cooling/smearing
steps. The topological susceptibility exhibits a surprisingly slow decrease at
high temperature.

We study finite (physical) volume and scaling violation effects of the Landau
gauge ghost propagator as well as of the running coupling $\alpha_s(p)$ in the
$SU(2)$ lattice gauge theory. We consider lattices with physical linear sizes
between $aL \simeq 3$ fm and $7$ fm and values of lattice spacing between
$a=0.2$ fm and $0.07$ fm. To fix the gauge we apply an efficient gauge fixing
method aimed at finding extrema as close as possible to the global maximum of
the gauge functional. We find finite volume effects to be small for the lattice
size $aL \simeq 3$ fm at momenta $p\, \ge \, 0.6$ GeV. For the same lattice
size we study extrapolations to the continuum limit of the ghost dressing
function as well as for the running coupling with momenta chosen between $p =
0.41$ GeV and $3.2$ GeV. We present fit formulae for the continuum limit of
both observables in this momentum range. Our results testify in favor of the
decoupling behavior in the infrared limit.

The phase diagram of twocolor QCD with nonzero chiral chemical potential is
studied by means of lattice simulation. We focus on the influence of a chiral
chemical potential on the confinement/deconfinement phase transition and the
breaking/restoration of chiral symmetry. The simulation is carried out with
dynamical staggered fermions without rooting. The dependences of the Polyakov
loop, the chiral condensate and the corresponding susceptibilities on the
chiral chemical potential and the temperature are presented. The critical
temperature is observed to increase with increasing chiral chemical potential.

We report on lattice QCD results for the thermodynamic equation of state of
quarkgluon matter obtained with Nf=2 degenerate quark flavors. For the fermion
field discretization we are using the Wilsontwisted mass prescription.
Simulations have been carried out at three values of the bare quark masses
corresponding to pion masses of 360, 430 and 640 MeV. We highlight the
importance of a good control of the lattice cutoff dependence of the trace
anomaly which we have studied at several values of the inverse temperature 1/T
= a Nt with a timelike lattice extent up to Nt=12. We contrast our results
with those of other groups obtained for Nf=0 and Nf=2+1. At low temperature we
also confront them with hadron resonance gas model predictions for the trace
anomaly.

The phase diagram of twocolor QCD with a chiral chemical potential is
studied on the lattice. The focus is on the confinement/deconfinement phase
transition and the breaking/restoration of chiral symmetry. The simulations are
carried out with dynamical staggered fermions without rooting. The dependence
of the Polyakov loop, the chiral condensate and the corresponding
susceptibilities on the chiral chemical potential and the temperature are
presented.

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.

Twocolor lattice QCD with N_f=4 staggered fermion degrees of freedom (no
rooting trick is applied) with equal electric charge q is studied in a
homogeneous magnetic background field B and at nonzero temperature T. In order
to circumvent renormalization as a function of the bare coupling we apply a
fixedscale approach. We study the influence of the magnetic field on the
critical temperature. At rather small pseudoscalar meson mass ($m_{\pi}
\approx 175 \mathrm{MeV} \approx T_c(B=0)$) we confirm a monotonic rise of the
quark condensate $<\bar{\psi} \psi>$ with increasing magnetic field strength,
i.e. magnetic catalysis, as long as one is staying within the confinement or
deconfinement phase. In the transition region we find indications for a
nonmonotonic behavior of $T_c(B)$ at low magnetic field strength ($qB<0.8
\mathrm{GeV}^2$) and a clear rise at stronger magnetic field. The conjectured
existence of a minimum value $T_c(B^{*}) < T_c(B=0)$ would leave a temperature
window for a decrease of $<\bar{\psi} \psi>$ with rising $B$ (inverse magnetic
catalysis) also in the present model.

Topological objects of $SU(3)$ gluodynamics are studied at the infrared scale
near the transition temperature with the help of zero and nearzero modes of
the overlap Dirac operator. We construct UV filtered topological charge
densities corresponding to three versions of the temporal boundary condition
applied to this operator, for which the zero mode is known to be located on
corresponding three constituent dyons (antidyons) in the reference case of an
analytical (anti)caloron solution. The clustering of the three topological
charge densities marks the positions of three types of dyons and antidyons
which can therefore be considered as present in equilibrium (Monte Carlo)
gluonic fields at the given resolution scale. We classify them either as
constituents of nondissociated (anti)calorons or as constituents of (anti)dyon
pairs or as isolated (anti)dyons. The pattern of the Polyakov loop describing
the centers and the interior of these clusters is observed after a limited
number of overimproved cooling steps and resembles the description known from
analytical caloron solutions.

In this talk we report on our study of twocolour lattice QCD with N_f=4
staggered fermion degrees of freedom with equal electric charge q in a
homogeneous magnetic field B at nonzero temperature T. We find indications for
a nonmonotonic behaviour of the critical temperature as a function of the
magnetic field strength and, as a consequence, for the occurence of `inverse
magnetic catalysis' within the transition region for magnetic fields in the
range 0 < qB < 0.7 GeV^2.

In this contribution we extend our unquenched computation of the Landau gauge
gluon and ghost propagators in lattice QCD at nonzero temperature. The study
was aimed at providing input for investigations employing continuum functional
methods. We show data which correspond to pion mass values between 300 and 500
MeV and are obtained for a lattice size 32**3 x 12. The longitudinal and
transversal components of the gluon propagator turn out to change smoothly
through the crossover region, while the ghost propagator exhibits only a very
weak temperature dependence. For a pion mass of around 400 MeV and the
intermediate temperature value of approx. 240 MeV we compare our results with
additional data obtained on a lattice with smaller Euclidean time extent N_t =
8, 10 and find a reasonable scaling behavior.

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.

We investigate the temperature dependence of the Landau gauge gluon and ghost
propagators in lattice QCD with two flavors of maximally twisted mass fermions.
For these propagators we provide and analyze data which corresponds to pion
mass values between 300 and 500 MeV. For the gluon propagator we find that both
the longitudinal and transversal component change smoothly in the crossover
region, while the ghost propagator exhibits only a very weak temperature
dependence. For momenta between 0.4 and 3.0 GeV we give a parametrization for
our lattice data. It may serve as input to studies which employ continuum
functional methods.

We continue the systematic computation of Landau gauge gluon and ghost
propagators of SU(2) gluodynamics using a sequence of increasing lattice sizes
L^4 up to L=112 with corresponding \betavalues chosen to keep the linear
physical size a(\beta)L ~ 9.6 fm fixed. To extremize the Landau gauge
functional we employ simulated annealing combined with subsequent
overrelaxation. Renormalizing the propagators at momentum \mu= 2.2 GeV we
observe quite strong lattice artifacts for the gluon propagator as well as for
the ghost dressing function within the momentum region q < 1.0 GeV. The
dependence on the lattice spacing for the gluon propagator at lowest accessible
physical momentum values does not yet allow a simple extrapolation to the
continuum limit. On the contrary, the running coupling derived from the bare
dressing functions seems less affected by lattice artifacts.