• ### Dyons near the transition temperature in $SU(3)$ lattice gluodynamics(1711.04476)

Nov. 15, 2017 hep-lat
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.
• ### Observation of deconfinement in a cold dense quark medium(1711.01869)

Nov. 6, 2017 hep-th, hep-ph, nucl-th, hep-lat
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 quark-gluon plasma is in essence a weakly interacting gas of quarks and gluons without a magnetic screening mass in the system, sharply different from a quark-gluon plasma at large temperature.
• ### Transversal and longitudinal gluon spectral functions from twisted mass lattice QCD with $N_f=2+1+1$ flavors(1710.04386)

Oct. 12, 2017 hep-lat
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 off-shell kinetics of the quark-gluon 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 net-baryon densities. There is evidence that a Quark-Gluon-Plasma (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 net-baryon densities and lower temperatures, it is expected that the QCD phase diagram exhibits a rich structure, such as a first-order 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 high-energy heavy-ion 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 net-baryon densities, because it is designed to run at unprecedented interaction rates. High-rate operation is the key prerequisite for high-precision measurements of multi-differential 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 baryon-chemical potentials (mu_B > 500 MeV), effects of chiral symmetry, and the equation-of-state 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 high-density QCD matter.
• ### Physical properties of Polyakov loop geometrical clusters in SU(2) gluodynamics(1606.04710)

Feb. 14, 2017 hep-lat
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.
• ### Geometrical clusterization of Polyakov loops in SU(2) lattice gluodynamics(1701.00536)

Dec. 31, 2016 hep-ph, hep-lat
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 liquid-gas phase transition. This finding may evidence for the fact that the SU(2) gluodynamics and this exactly solvable model of nuclear liquid-gas phase transition are in the same universality class.
• ### Study of the phase diagram of dense two-color QCD within lattice simulation(1605.04090)

Dec. 16, 2016 hep-ph, nucl-th, hep-lat
In this paper we carry out a low-temperature scan of the phase diagram of dense two-color 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 Bose-Einstein condensation of scalar diquarks. In this phase the system is still in confinement in conjunction with non-zero 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$.
• ### Dyons and Roberge - Weiss transition in lattice QCD(1611.07789)

Nov. 23, 2016 hep-lat
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.
• ### Dyons near the transition temperature in lattice QCD(1512.03217)

Dec. 22, 2015 hep-lat
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 time-like Abelian monopoles and specific KvBLL pattern of Polyakov line correlate with topological clusters.
• ### Study of QCD Phase Diagram with Non-Zero Chiral Chemical Potential(1512.05873)

Dec. 18, 2015 hep-lat
In this paper we report on lattice simulations of SU(3)-QCD with non-zero 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.
• ### New Signals of Quark-Gluon-Hadron Mixed Phase Formation(1510.03099)

Oct. 11, 2015 nucl-th
Here we present several remarkable irregularities at chemical freeze-out 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 freeze-out at the center of mass energy 4.9 GeV and two sets of highly correlated quasi-plateaus 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.8-4.9 GeV and 7.6-10 GeV. The low energy set of quasi-plateaus was predicted a long time ago. On the basis of the generalized shock-adiabat model we demonstrate that the low energy correlated quasi-plateaus give evidence for the anomalous thermodynamic properties inside the quark-gluon-hadron mixed phase. It is also shown that the trace anomaly sharp peak at chemical freeze-out corresponds to the trace anomaly peak at the boundary between the mixed phase and quark gluon plasma. We argue that the high energy correlated quasi-plateaus 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.
• ### Towards the quark--gluon plasma Equation of State with dynamical strange and charm quarks(1510.02262)

Oct. 8, 2015 hep-lat
We present an ongoing project aimed at determining the thermodynamic Equation of State (EoS) of quark--gluon 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 fixed--lattice 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.
• ### Topological susceptibility from $N_f=2+1+1$ lattice QCD at nonzero temperature(1510.02265)

Oct. 8, 2015 hep-lat
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 over--improved 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.
• ### Landau gauge ghost propagator and running coupling in SU(2) lattice gauge theory(1302.5943)

Oct. 8, 2015 hep-lat
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.
• ### Two-Color QCD with Non-zero Chiral Chemical Potential(1503.06670)

March 23, 2015 hep-lat
The phase diagram of two-color QCD with non-zero 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.
• ### The equation of state of quark-gluon matter from lattice QCD with two flavors of twisted-mass Wilson fermions(1412.6748)

Dec. 21, 2014 hep-lat
We report on lattice QCD results for the thermodynamic equation of state of quark-gluon matter obtained with Nf=2 degenerate quark flavors. For the fermion field discretization we are using the Wilson-twisted 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 time-like 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.
• ### Two-Color QCD with Chiral Chemical Potential(1411.5174)

Nov. 19, 2014 hep-lat
The phase diagram of two-color 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.
• ### Dyon structures in the deconfinement phase of lattice gluodynamics: topological clusters, holonomies and Abelian monopoles(1410.4632)

Oct. 17, 2014 hep-lat
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.
• ### Magnetic catalysis (and inverse catalysis) at finite temperature in two-color lattice QCD(1310.7876)

Feb. 8, 2014 nucl-th, hep-lat
Two-color 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 non-zero temperature T. In order to circumvent renormalization as a function of the bare coupling we apply a fixed-scale approach. We study the influence of the magnetic field on the critical temperature. At rather small pseudo-scalar 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 non-monotonic 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.
• ### Topology near the transition temperature in lattice gluodynamics analyzed by low lying modes of the overlap Dirac operator(1309.7850)

Feb. 5, 2014 hep-th, hep-lat
Topological objects of $SU(3)$ gluodynamics are studied at the infrared scale near the transition temperature with the help of zero and near-zero 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.
• ### Two-colour QCD at non-zero temperature in the presence of a strong magnetic field(1402.0112)

Feb. 1, 2014 hep-lat
In this talk we report on our study of two-colour lattice QCD with N_f=4 staggered fermion degrees of freedom with equal electric charge q in a homogeneous magnetic field B at non-zero temperature T. We find indications for a non-monotonic 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.
• ### Landau gauge gluon and ghost propagators from two-flavor lattice QCD at T > 0(1312.1914)

Dec. 6, 2013 hep-lat
In this contribution we extend our unquenched computation of the Landau gauge gluon and ghost propagators in lattice QCD at non-zero 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.
• ### Topology across the finite temperature transition studied by overimproved cooling in gluodynamics and QCD(1304.0935)

April 3, 2013 hep-th, hep-lat
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.
• ### Landau gauge gluon and ghost propagators from lattice QCD with Nf=2 twisted mass fermions at finite temperature(1212.1102)

March 15, 2013 hep-ph, hep-lat
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.
• ### On the continuum limit of Landau gauge gluon and ghost propagators in SU(2) lattice gauge gluodynamics(1303.3423)

March 14, 2013 hep-lat
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 \beta-values 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.