Temperature Dependent Hadronic Bag and QGP Phase Transition in Dual QCD

Based on the magnetic symmetry structure of non-Abelian gauge theories, a dual QCD gauge theory has been constructed which takes into account the local structure as well as the topological features of the color gauge group into its dynamics in a completely dual-symmetric way. Using such dual version of QCD in thermal domain following the partition function approach and the grand canonical ensemble formulation, the phase transition from hadron to QGP phase has been investigated within the framework of temperature dependent hadronic bag in the entire T − μ plane. The various thermodynamic properties like pressure, energy density, speed of sound and specific heat of the hadron/QGP phase have been evaluated and shown to lead an evidence for the first order phase transition. In the region around Tc < T < 4Tc , the specific heat and speed of sound are strongly influenced by the magnetically charged particles directly related to thermal monopoles evaporating from the magnetic condensate present at low temperature.


Introduction
Soon after the discovery of Quantum Chromodynamics (QCD) [1,2,3] and following the realization that QCD exhibits asymptotic freedom [4] it was recognized that normal hadronic matter undergoes a phase transition, where the individual hadrons dissolve into their constituents and produce a collective form of matter known as the Quark-Gluon Plasma (QGP) [5] under extreme conditions of high temperature and low chemical potential ( µ B ). The precise determination of the phase boundary between QGP and hadron gas (HG) at high T and bariochemical potential ( µ B ) has been a subject of intense research in recent years from experimental as well as theoretical point of view. In this connection, during past few years, the possibility of creating such high temperature QGP and studying QGP phase of matter by colliding heavy ions in the laboratory has been the main goal of experiments at Nuclotron-based ion collider facility (NICA) [6,7] at Dubna, the Relativistic Heavy-Ion Collider (RHIC) [8,9] at BNL and the Large Hadron Collider (LHC) [10,11,12] at CERN and needs some reliable theoretical explanation of various signals which depend on the pressure, entropy, transition temperature and the equation of state. Unfortunately, the theoretical description of the QGP from first principle is extremely difficult as the interaction between quarks and gluons as described by QCD, is strong, therefore, perturbative QCD is not applicable. Being in the non-perturbative domain, the lattice QCD [13,14,15,16,17] and the effective models remain the only tools to study the phase transition. The lattice-QCD simulation is, however, of limited practical use at finite baryon density due to the well known fermion sign problem. Therefore, the critical behavior of such matter with finite baryon chemical potential needs to be studied in the framework of QCD motivated phenomenological models such as MIT bag model [18] which provides a semi phenomenological description of an EOS that features a quark-hadron transition.
The present paper mainly deals with the investigation of thermodynamic properties of QGP using an infrared effective dual QCD formulation based on magnetic symmetry. In section "Dual QCD formalism with magnetic symmetry", the dual QCD formulation has been presented and analyzed to explain the confining behavior of QCD vacuum in its non-perturbative sector. In section "QGPphase transition dynamics and thermal hadronic bags", the QGP-phase transformation dynamics with a thermal version of hadronic bags has been studied in order to investigate the thermodynamics quantities like the pressure, the energy density, the specific heat and the square of speed of sound in QGP phase. Calculation of critical parameters involves the use of Gibbs's criteria of thermodynamic equilibrium. The results and their implications have been discussed in section "Evaluation of thermodynamic and transport properties of QGP".

Dual QCD formalism with magnetic symmetry
Based on the magnetic symmetry structure of non-Abelian gauge groups, a dual gauge formulation [19,20,21,22,23,24,25,26,27,28,29,30,31,32,33]. which produce the magnetic condensation dynamically for its vacuum and provides a gauge invariant investigation and topological ground of confinement, the resulting dual QCD has been a subject of utmost importance to analyze the Quark-gluon plasma phase of the hadronic matter. It is based on introducing the magnetic symmetry as an internal isometry H admitting some additional Killing vector fields which are internal such that H is a Cartan's subgroup of a gauge group G and commutes with it [20,29]. It, in turn, restricts the associated gauge potential which may be expressed in the following form, where W µ is the gauge potential of the gauge group G . It leads to the following form given as is the Abelian component and m may then be viewed to define the homotopy class of the mapping Π 2 (S 2 ) as, m : S 2 R → S 2 = SU(2)/U(1) , where, S 2 R is the two dimensional sphere of the three dimensional space and S 2 is the group coset space completely fixed by m . Using the dual magnetic potential B (d) µ for the magnetic part in order to remove the problem of singular behavior of the potential associated with monopoles, the dual QCD Lagrangian is constructed in the following form [29].
where V(ϕ) is the effective potential responsible for the dynamical breaking of magnetic symmetry and has its quadratic form reliable in the phase transition study of dual QCD vacuum as given by, Using the cylindrical symmetry the field equations associated with the Lagrangian (2) in the quenched approximation may be derived in the following form [29], Table 1.
The masses of vector and scalar glueballs as a functions of α s .
Imposing the asymptotic boundary conditions appropriate for the large-scale be- , where the function F(ρ) , in asymptotic limit may then be obtained in the following form, Utilizing the asymptotic solutions of the associated dual QCD fields [29,30], the energy per unit length of the resulting flux tube configuration may be given as, The associated vector and scalar glueball masses generated after the dynamical breaking of magnetic symmetry may be evaluated using relations, associated with the string tension and the Regge slope parameter. The numerical results of the glueball masses, characteristic length scales and Ginzburg-Landau (GL) parameter [20] are shown in table 1. The estimate for the GL parameter ( κ (d) QCD demonstrate different superconducting behavior (type-I, II) of the QCD vacuum intimately connected to the process of quark confinement in the low energy sector. It is therefore, strongly desired to construct a low energy effective theory of QCD which might explain the correct confining features of QCD on one hand and complete phase transition behavior on the other hand, in an effective way.

QGP-phase transition dynamics and thermal hadronic bags
Phase transitions are associated with the evolution of thermodynamic quantities such as pressure, energy density and entropy density, as well as a set of response functions, like, specific heat and speed of sound. In order to study the phase transition dynamics, let us begin with a relativistic quantum system where particles can be efficiently created and destroyed. Using the grand canonical ensemble formalism the corresponding partition function is given in the following form, where, β = 1/T and H and N is the Hamiltonian and the particle number operator respectively. The pressure P and the energy density ǫ for a spatial homogeneous system are given by the first derivatives of ln Z as, Furthermore, other observables of phenomenological interest derivable from the basic thermodynamic quantities P and ǫ are the normalized specific heat and square of speed of sound expressed in the following form, In order to calculate the various thermodynamic quantities let us start with the resulting expression for bosonic partition function written as [32], while that for fermions is [32], g b and g f are the degeneracy factors for bosons and fermions, respectively, and m b and m f denote the masses of the particles. Hence, using the above expression (12) with the degeneracy factor g f = 2 × 2 for nucleons, the pressure and energy density for the hadronic matter is expressed in the following form, Having a description of the hadronic phase at hand, it is straightforward to analyze the thermodynamic properties of the quark-gluon plasma and to study the dynamics of hadron to quark-gluon phase transition. Possibly, the simplest model to describe an approximate physics describing the matter where quark and gluons are the proper degrees of freedom of the system is the MIT bag model [18]. The non-trivial feature of the bag model EOS is that the vacuum pressure generates positive and large contribution to the bag energy density as given by, where, B represents the bag constant which gives the bag pressure. Hence, the grand canonical partition function for the QGP phase consisting of a perturbatively interacting gas of quarks and gluons can be written in the following form, where, in the r.h.s the first term is the contribution at free level while the second term represents the non-perturbative vacuum contribution in the form of a Tdependent bag constant B(T) which may be identified by the confining part of the energy expression recalculated after taking the multi-flux tube system as a periodic system on a S 2 -sphere [30,32] and is given by, where m

(T) B
B is the thermal vector glueball mass derived using path-integral formalism alongwith the mean-field approach [31] and is given by the following expression, where m (0) B is the non-thermal contribution to the vector glueball mass. With these considerations, resulting expressions for the energy density and pressure for QGP phase is then given as, The thermodynamical phase stability then requires the equality of pres-sure, chemical potential at the boundary of two phases i.e., P h = P p = P c , µ = 3µ q = µ c at the transition point T p = T h = T c . Since at the critical temperature the additional degrees of freedom carried by the quark-gluon plasma, are to be released which results in an increase in the thermodynamic quantities and is investigated in the next section.

Evaluation of thermodynamic and transport properties of QGP
The critical parameters of QGP phase transition within the framework of thermal bag may be evaluated using above mentioned phase equilibrium criteria due to Gibbs which may be used with the 3-d pressure as depicted in Figure 1(a) for different coupling values in infrared sector of QCD. The change of the quarkchemical potential from zero to non-zero values defines uniquely and precisely the QGP phase transition temperatures and their values are found to decrease with increasing values of chemical potential. At low values of µ and T the nuclear matter is composed of confined hadrons, but with increasing T or µ , the hadronic matter undergoes a phase transition towards a plasma of deconfined quarks and gluons. At vanishing chemical potential ( µ q = 0) the transition temperature of 0.187 GeV , 0.140 GeV , 0.116 GeV and 0.090 GeV for the case of couplings α s = 0.12; 0.22; 0.47 and 0.96 are obtained respectively.
In Figure 1(b), the variation of energy density for hadron and QGP phase has been depicted and it has been observed that the energy density increases abruptly at the phase transition temperature followed by a finite jump discontinuity indicating the rapid rise in the degree of freedom carried by quarks and gluons. The size of the discontinuity in the energy density, the so-called latent heat ( ∆ǫ ) is given by, Another important thermodynamic quantity characterizing the equation of state for a system undergoing phase transition is the specific heat which is a measure of energy uctuations in the system and these uctuations tend to rise sharply near a phase transition as shown in Figure 2. The thermal evolution of normalized specific heat shows a small upward cusp at T c and above the critical temperature specific heat rises with the temperature due to emergent chromomagnetic monopoles. In this scenario, magnetically charged particles are important component above the phase transition, possibly contributing to the physical properties of the strongly interacting QGP [34,35,36,37,38]. These component of the deconfined plasma are, in turn, directly related to thermal monopoles evaporating from the magnetic condensate present at low temperature. Furthermore, another observable of phenomenological interest is the square of speed of sound in the hot medium related to the speed of small perturbations produced in the QCD matter. The thermal evolution of the square of speed of sound with temperature and chemical potential has been depicted in Figure 3 and a first order phase transition is indicated as a result of a sudden drop of c 2 s around transition temperature alongwith a further rise with temperature approaching to the value c 2 s = 0.33 near T c . In the region around T c < T < 4T c , c 2 s decreases with temperature and in this scenario, the near T c QCD matter is considered as a semiquark-gluon-monopole plasma (sQGMP) [39] that contains not only electrically charged quasi-particles, quark and gluons, but also magnetically charged quasi particles, monopoles. In conclusion, this scenario emphasizes the change in chromodegrees of freedom, and recasts the QCD phase diagram into electrically and magnetically dominated regimes.

Conclusion
Based on the magnetic symmetry structure of non-Abelian gauge theories, a dual QCD gauge theory has been constructed which takes into account the local structure as well as the topological features of the color gauge group into its dynamics in a completely dual-symmetric way. Using such dual version of QCD in thermal domain following the partition function approach and the grand canonical ensemble formulation, the phase transition from hadron to QGP phase has been investigated within the framework of temperature dependent hadronic bag in the entire T − µ plane. The various thermodynamic properties like pressure, energy density, speed of sound and specific heat of the hadron/QGP phase have been evaluated and shown to lead an evidence for the first order phase transition. In the region around T c < T < 4T c , the specific heat and speed of sound are strongly influenced by the magnetically charged particles directly related to thermal monopoles evaporating from the magnetic condensate present at low temperature.