Thermodynamic Investigation of the QCD Phase Diagram with 2+1 Quark Flavors †

: This work deals with the deconﬁnement phase transition from a hadronic gas (HG) phase consisting of massive pions, to a quark–gluon plasma (QGP) phase consisting of gluons, massless up and down quarks and massive strange quarks, in addition to their antiquarks. Based on the Bag and coexistence models, we study the variations of pressure characterizing both HG and QGP phases. For the latter, we calculate the partition function of the color-singlet QGP within the projection method using a density of states containing the volume term only. We investigate the phase diagram of the strongly interacting matter, in the µ – T plane, in several cases: in the HG phase, we consider massless pions then we account for their masses, and in the QGP phase, ﬁrst we take it consisting of two massless u and d quarks, then we consider additional massive strange quarks.


Introduction
Quantum chromodynamics (QCD) is the basic theory that describes the strong interactions between quarks and gluons in the hadrons [1]. QCD predicts the asymptotic freedom of short-range quarks (distances < 1 fm) inside the hadrons, in addition to the confinement property, which means that under ordinary conditions, quarks and gluons are confined to the hadrons. However, at high temperatures and/or density, hadronic matter undergoes a phase transition: quarks and gluons become deconfined, i.e., free [2], in a new phase called the quark-gluon plasma (QGP).The existence of this phase of matter was proposed in the mid-seventies, just ten years after the birth of the quark model of hadrons [3,4].
To illustrate the phase transition from the hadronic gas (HG) phase to the QGP phase, we will investigate the phase diagram of the strongly interacting matter, in the µ-T plane, based on the Gibbs criterion setting the mechanical equilibrium between the two phases at the transition, with massless and massive particles in the two phases.

Phase Diagram in the µ-T Plane
According to the Gibbs criterion, the HG and the QGP phases are in equilibrium when their pressures P, temperature T and chemical potential µ are equal, which reads: Thus, we have to calculate the pressure of both phases, using the thermodynamic relation: from the individual partition functions Z HG and Z QGP of the two phases. For the HG phase, we consider an ideal gas of pions (π + , π − and π o ), which are the lightest bosons, and the obtained partition function reads: where k is the momentum, d π the pion degeneracy factor, m π the pion mass and β = 1 T . The QGP being considered asan ideal gas of gluons, quarks and their antiquarks, and assuming non-interaction between its constituents, we can write Z QGP as: where: , accounts for the real vacuum pressure B exerted on the perturbative vacuum, B being the bag constant. The calculation of Z quark and Z gluon gives: Setting the condition P HG = P QGP , with B 1/4 = 200 MeV, we obtain the phase diagram in the µ − T plane shown in Figure 1, characterized by a critical line separating the HG and the QGP phases and giving at each point the transition parameters µ c , T c , in all the studied cases of massless and massive particles. When T < T c , the HG is favored, and when T > T c , the QGP is favored. There are two extreme cases: the first at µ = 0, obtainable in URHIC, and the second at T = 0, which could be reached in the core of certain neutron stars. It can also be noted by comparing the red dashed curve (with massless pions in the HG phase) and the solid blue curve (with massive pions), in Figure 1b, that there is a slight difference in T c at µ = 0, meaning that the pion mass has to be considered for accurate investigation of the transition temperature, while the small mass of u and d quarks can simply be neglected. Moreover, by examining the dashed-dotted magenta line, obtained by considering massive s quarks additionally to massless u and d quarks in the QGP phase, a clear mismatch appears compared to the two other curves, showing that the number of flavors considerably affects the transition parameters, and particularly the T c value at µ = 0.

Phase Transition to a Color-Singlet QGP
To implement the color-singletness constraint into the quantum statistical description of the QGP, we use the group theoretical projection formulated by Turko and Redlich [5]. The projected partition function of a QGP on the color-singlet SU(3) representa-

Phase Transition to a Color-Singlet QGP
To implement the color-singletness constraint into the quantum statistical description of the QGP, we use the group theoretical projection formulated by Turko and Redlich [5]. The projected partition function of a QGP on the color-singlet SU(3) representation, in a volume V, at temperature T and chemical potential µ reads: g 1 , g 2 and g 3 result for massless u and d quarks, massive s quarks and massless gluons, respectively, In the following, we study the variations of the pressures P HG and P QGP with temperatures T, at µ = 0, for B 1/4 =200 MeV, in different cases of massless and massive particles in both HG and QGP phases. Figure 2 shows the curve of the HG pressure, for massless pions (red dotted line) and massive pions (blue dotted line) as well as that of the QGP for massless u and d quarks only (dashed lines), then when adding the massive s quarks (solid lines), for a small and a large volume in both of the last cases.  In Figure 2, we studied the effect of the particle masses for small and large volume of the QGP system on the deconfinement phase transition temperature at µ = 0, when considering the color-singletness requirement for the QGP phase, by the extraction of the value of from the intersection point of the curves and . The obtained values of in the different studied cases are shown in Table 1, and we can easily see that as the volume increases, decreases until the volume reaches 10,000 fm 3 , where begins to stabilize. We can also notice that (massless particles) > (massive particles).  In Figure 2, we studied the effect of the particle masses for small and large volume of the QGP system on the deconfinement phase transition temperature T c at µ = 0, when considering the color-singletness requirement for the QGP phase, by the extraction of the value of T c from the intersection point of the curves P HG and P QGP . The obtained values of T c in the different studied cases are shown in Table 1, and we can easily see that as the volume increases, T c decreases until the volume reaches 10,000 fm 3 , where T c begins to stabilize. We can also notice that T c (massless particles) > T c (massive particles).