## Deconfinement and freezeout boundaries in equilibrium thermal models

Abdel Nasser Tawfik ∗

Nile University, Egyptian Center for Theoretical Physics (ECTP), Juhayna Square of 26th-July-Corridor, 12588 Giza, Egypt and Goethe University, Institute for Theoretical Physics (ITP), Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany

Muhammad Maher, A. H. El-Kateb, and Sara Abdelaziz Helwan University, Faculty of Science, Physics Department, 11795 Ain Helwan, Egypt (Dated: January 30, 2020)

In different approaches, the temperature-baryon density plane of QCD matter is studied for deconfinement and chemical freezeout boundaries. Results from various heavy-ion experiments are compared with the recent lattice simulations, the effective QCD-like Polyakov linear-sigma model, and the equilibrium thermal models. Along the entire freezeout boundary, there is an excellent agreement between the thermal model calculations and the experiments. Also, the thermal model calculations agree well with the estimations deduced from the Polyakov linear-sigma model (PLSM). At low baryonic density or high energies, both deconfinement and chemical freezeout boundaries are likely coincident and therefore the agreement with the lattice simulations becomes excellent as well, while at large baryonic density, the two boundaries become distinguishable forming a phase where hadrons and quark-gluon plasma likely coexist.

PACS numbers: 12.40.Ee, 12.40.Yx, 05.70.Ce

Keywords: Statistical models, Hadron mass models and calculations, Thermodynamic functions and equa- tions of state

## I. INTRODUCTION

Strongly interacting matter under extreme conditions is characterized by different phases and different types of the phase transitions [1]. The hadronic phase, where stable baryons build up a great part of the Universe and the entire everyday life, is a well known phase. At high temperatures and/or densities, other phases appear. For instance, at temperatures of a few MeV, chiral symmetry restoration and deconfinement transition take place, where quarks and gluons are conjectured to move almost freely within colored phase known as the quark-gluon plasma (QGP) [2]. At low temperatures but large densities, the hadronic (baryonic) matter forming compact interstellar objects such as neutron stars is indubitably observed in a conventional way and very recently gravitational waves from neutron star mergers have been detected, as well [3]. At higher densities, extreme interstellar objects such as quark stars are also speculated [4]. In lattice quantum chromodynamics (QCD), different orders of chiral and deconfinement transitions have be characterized, especially at low baryon densities.

The program of heavy-ion collision experiments dates back to early 1980's. Past (AGS, SIS, SPS), current (RHIC, LHC), and future facilies (FAIR, NICA) help in answering essential questions about the thermodynamics of the strongly interacting matter and in mapping out the temperature-baryon density plane [2]. The unambiguous evidence on the formation of QGP is an example of a great imperical achievement [5, 6]. The colliding nuclei are conjectured to form a fireball that cools down by rapid expansion and finally hadronizes into individual uncorrelated hadrons. The present script focuses on the temperature-baryon density plane, concretely near the hadron-QGP boundaries, in framework of equilibrium thermal model [7]. To this end, we put forward a basic assumption that both directions, hadron-QGP and QGP-hadron phase are quantum-mechanically allowed [8]. In other words, the picture drawn so far seems in fundamental conflict with the time arrow. The concept of arrow of time prevents the reverse direction, especially if the change in the degrees of freedom or entropies

∗ Electronic address: atawfik@nu.edu.eg

aren't following the causality principle; second law of thermodynamics. The statistical thermal approaches work well near to both deconfinement and chemical freezeout boundaries [2, 9]. This could be understood in the light of the themal nature of an arbitrary small part of the highly entangled fireball states. Following the Eigenstate Thermalization Hypothesis [8, 10], the corresponding probability distribution of the projection of these states remains thermal. We follow the line that the thermal models reproduce well the particle yields and the thermodynamic properties of the hadronic matter including the chiral and freezeout temperatures. We compare our calculations with reliable lattice QCD simulations, an effective QCD-like approach, and available experimental results.

The present script is organized as follows. In section II approaches for deconfinement and freezeout boundaries in equilibrium thermal models are introduced. The results are discussed in section III. Section IV is devoted to the conclusions and outlook.

## II. EQUILIBRIUM THERMAL MODELS

It was conjectured that the formation of the hadron resonances follows the bootstrap picture, i.e. the hadron resonances or the fireballs are composed of further resonances or fireballs, which in turn are consistent of lighter resonances or smaller fireballs and so on [11, 12]. The thermodynamic quantities of such a system can be deduced from the partition function Z T,µ,V ( ) of an ideal gas. In a grand canonical ensemble, this reads [2, 13-17]

$$Z ( T, V, \mu ) = \text{Tr} \left [ \exp \left ( \frac { \mu N - H } { T } \right ) \right ],$$

where H is Hamiltonian combining all relevant degrees of freedom and N is the number of constituents of the statistical ensemble. Eq. (1) can be expressed as a sum over all hadron resonances taken from recent particle data group (PDG) [18] with masses up to 2 5 GeV [19], .

$$\ln Z ( T, V, \mu ) = \sum _ { i } \ln Z _ { i } ( T, V, \mu ) = V \frac { g _ { i } } { 2 \pi ^ { 2 } } \int _ { 0 } ^ { \infty } \pm p ^ { 2 } d \text{p} \ln \left [ 1 \pm \lambda _ { i } \exp \left ( \frac { - \varepsilon _ { i } ( p ) } { T } \right ) \right ],$$

where the pressure reads T∂ ln Z T,V,µ /∂V ( ) , ± stands for fermions and bosons, respectively. ε i = ( p 2 + m 2 i ) 1 / 2 is the dispersion relation and λ i is the fugacity factor of i -th particle [2],

$$\lambda _ { i } ( T, \mu ) = \exp \left ( \frac { B _ { i } \mu _ { \mathfrak { b } } + S _ { i } \mu _ { \mathfrak { s } } + Q _ { i } \mu _ { \mathfrak { q } } } { T } \right ),$$

where B µ i ( b ), S i ( µ S ), and Q µ i ( Q ) are baryon strangeness, and electric charge quantum numbers (their corresponding chemical potentials) of the i -th hadron, respectively. From phenomenological point of view, the baryon chemical potential µ b can be related to the nucleon-nucleon center-of-mass energy √ s NN [20]

$$\mu _ { \text{b} } \, = \, \frac { a } { 1 + b \sqrt { s _ { \text{BM} } } },$$

where a = 1 245 . ± 0 049 GeV and . b = 0 244 . ± 0 028 GeV . -1 . The number and energy density, respectively, can be derived as

$$n _ { i } ( T, \mu ) \, & = \, \sum _ { i } \frac { \partial \ln Z _ { i } ( T, V, \mu ) } { \partial \, \mu _ { i } } = \sum _ { i } \frac { g _ { i } } { 2 \pi ^ { 2 } } \int _ { 0 } ^ { \infty } p ^ { 2 } d p \frac { 1 } { \exp \left [ \frac { \mu _ { i } - \varepsilon _ { i } ( p ) } { T } \right ] \pm 1 }, \\ \nu \quad & \quad \varepsilon \, \quad \varepsilon \, \partial \ln Z _ { i } ( T, V, \mu ) \quad \varepsilon \, \quad q _ { i } \quad f ^ { \infty } \quad \varepsilon \,. \quad - \varepsilon _ { i } ( p ) \pm \mu _ { i } \quad \nu \quad \nu \quad \nu.$$

$$\rho _ { i } ( T, \mu ) \ = \ \sum _ { i } \frac { \partial \ln Z _ { i } ( T, V, \mu ) } { \partial ( 1 / T ) } = \sum _ { i } \frac { g _ { i } } { 2 \pi ^ { 2 } } \int _ { 0 } ^ { \infty } p ^ { 2 } d p \frac { \binom { \iota } { - \varepsilon _ { i } ( p ) \pm \mu _ { i } } } { \exp \left [ \frac { \mu _ { i } - \varepsilon _ { i } ( p ) } { T } \right ] \pm 1 }. \\ \cdot \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \.. \.$$

Likewise, the entropy and other thermodynamic quantities can be derived straightforwardly.

Both temperature T and the chemical potential µ = B µ i b + S µ i S + · · · are related to each other and to √ s NN [2]. As an overall thermal equilibrium is assumed, µ S is taken as a dependent variable to be estimated due to the strangeness conservation, i.e. at given T and µ b , the value assigned to µ S is the one assuring 〈 n S 〉-〈 n ¯ S 〉 = 0. Only then, µ S is combined with T and µ b in determining the thermodynamic quantities, such as the particle

number, energy, entropy, etc. The chemical potentials related to other quantum charges, such as the electric charge and the third-component isospin, etc. can also be determined as functions of T , µ b , and µ S and each of them must fulfill the corresponding laws of conversation.

This research intends to distinguish between deconfinement and freezeout boundaries in equilibrium thermal models. The latter is characterized by T χ and µ b , which are conditioned to one of the universal freezeout conditions [21], such as constant entropy density normalized to T 3 χ [22, 23], constant higher-order moments of the particle multiplicity [24, 25], constant trace anomaly [26] or an analogy of the Hawking-Unruh radiation [27]. The experimental estimation for T χ and µ b , as shown in Fig. 1 proceeds through statistical fits for various particle ratios calculated in statistical thermal models. The former, the deconfinement transition, is conditioned to line-of-constant-physics, such as constant energy density, ρ [28]. The inclusion of the strange quarks seems to affect the critical temperatures, as these come up with extra hadron resonances and their thermodynamic contributions, where the mass of strange quarks is of the order of the critical temperature.

## III. RESULTS

Fig. 1: The freezout and deconfinement parameters T χ and µ b as deduced from different experimental results (symbols with error bars) are depicted and confronted with the thermal model calculations for the chemical freezeout (dot-dashed and dashed curves) and the deconfinement parameters (solid and long-dashed curves) with and without strange quarks.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00426v2-figure-f5122405-6f15-4944-b51f-d525e1c210bb.jpg)

Figure 1 depicts the freezout and deconfinement parameters T χ and µ b as determined from the measurements gained from different experiments (symbols with error bars); HADES (Ar+KCal) [29], NA61/SHINE (Au+Au, p+p) [30-36], NA49/SHINE (Pb+Pb) [37-39], and ALICE [40-47] as well as analysis at RHIC/SPS/AGS [31] and RHIC/SPS/AGS/SIS energies [48-54] are combined with each others and compared with the thermal model calculations. The latter take into account both freezeout (dot-dashed and dashed curves) and deconfinement boundaries (solid and long-dashed curves) with and without strange quarks.

With the experimental results, we mean the parameters obtained when measured particle yields and/or ratios are fitted to calculations based on statisical thermal models, in which the parameters T χ and µ b are taken as independent variables, where the baryon-chemical potential µ b , for instance, can directly be fixed at a given center-of-mass energy √ s NN , Eq. (4). For the freezeout parameters T χ , µ b , µ S , etc. the thermodynamic quantities which fulfill one of the freezeout conditions reviewed in refs. [16, 21], such as constant entropy density normalized to temperature cubed. Conditions for deconfinement phase transitions have been discussed in refs. [16, 28], line-of-constant-physics, such as constant energy density with varying µ b , µ S , √ s NN , etc. Details of

the various approaches (curves) become in order now. With freezeout q l and freezeout q l + q s we mean QCD phase boundaries as determined under freezeout conditions, where in the HRG model is configured to have only hadrons whose constituents are light quarks ( q l ) and to have only hadrons with light and strange quarks q l + q s . Under these conditions, it is likely to characterize the chemical freezeout T -µ b -plane. Furthermore, with deconfine q l and deconfined q l + q s , we mean the line-of-constant-physics, which is defined by constant energy density, for instance. Such a line is mapped out in the HRG model in which only hadrons with light quarks ( q l ) and only hadrons with light and strange quarks q l + q s are included.

It is obvious that both sets of freezeout parameters seem identical, for instance, at low µ b or high √ s NN , where the slight difference could be tolerated. At large µ b or low √ s NN , the difference between the temperatures of freezeout and deconfinement becomes larger. Such a difference would be understood based on the assumption that the chemical freezeout takes place very late after the phase of hadronization. The latter is QCD confinement transition. Its order as simulated in recent lattice QCD is a likely crossover, i.e. there a wide range of temperatures within which QGP hadronizes or hadrons go through QGP. The time span becomes longer with the increase in µ b or the decrease in √ s NN . The conjecture of the existence of a mixed phase is probably another possibility. In this phase, both types of degrees of freedom, hadrons and QGP, live together until the system goes through deconfinement to colored QGP or finally entirely freezes out to uncorrelated colorless hadrons.

The co-existence of different QCD phases was discussed in litrature, for instance [55, 56]. The mixed QCD phases can be formed in macroscopic, mesoscopic, and microscopic mixture. As shown in Fig. 1, these mixed phases start being produced at √ s NN ranging between ∼ 5 and ∼ 12 GeV, i.e. µ b /similarequal 320 to /similarequal 560 MeV.

For the freezeout parameters, it is apparent that the agreement between the thermal model calculations and the experimental results is very convincing. This covers the entire µ b -range and can - among other evidence -be interpreted based on the fact that the freezeout stage is the latest along the temporal evolution of the high-energy collision, where the number of produced particles is entirely fixed. The time elapsed from the stage of the chemical freezeout to the time of detection is likely shorter than the time from any other QCD processes, such as hadronization, chiral symmetry breaking, etc. and therefore it is apparently the most accurate one.

In the present calculations, full quantum statistics [13-17] and hadron resonances with masses up to 2 5 GeV . [19] are taken into account. The strangeness degrees of freedom play an important role, especially at low µ b or high √ s NN .

For the sake of completeness, we have also checked the same calculations but with the inclusion of a large number of possible missing states [57, 58]. We found that the thermodynamic quantities, especially the ones to which the present script is limited, show sensitivity to these missing states [59]. They are entering our calculations in the same manner as done for the PDG hadrons and resonances.

The missing states are resonances predicted, theoretically, but not yet confirmed, experimentally. Their quantum numbers and physical characteristics are theoretically well known [60]. Basically, they are conjectured to greatly contribute to the fluctuations and the correlations, i.e. higher derivatives of the partition function, estimated in recent lattice QCD simulations [60]. These are the occasions where their contributions becomes unavoidable [57]. Another reason for adding the missing states is that they come up with additional degrees of freedom and considerable decay channels even to the hadrons and resonances which are subject of this present study.

For T χ and µ b , a comprehensive comparison between the thermal model calculations (curves) and the results deduced from the lattice QCD simulations (bands) [61, 62] and the Polyakov linear-sigma model (symbols with error bars) [63] is presented in Fig. 2. Within their statistical and systematic certainities, there is an excellent agreement between the lattice QCD simulations (bands) [61, 62] and the Polyakov linear-sigma model (symbols) [63]. The reason why the lattice QCD simulations are limited to µ /T q χ ≤ 1 is the so-called sign problem and the difficulties which arise because of the importance of sampling becomes no longer possible. There are various attempts to anticipate this limitation; continuation from imaginary chemical potential, reweighting methods, applying complex Langevin dynamics, and Taylor expansions in the quark chemical potential µ q [64].

We also find an excellent agreement between the thermal model calculations for the chemical freezeout parameters and the predictions deduced from the Polyakov linear-sigma model, especially at µ b &gt; ∼ 300 MeV. At lower µ b , the thermal model calculations seem to slightly overestimate T χ .

This observed agreement would be taken as an evidence supporting the conclusion that the first-principle calculations likely result in T χ -µ b plane similar to that of the Polyakov linear-sigma model, especially at large µ b , where the first-principle calculations are no longer applicable.

It is in order now to highlight a few details of the linear-sigma model, which is much simpler than QCD, but based on QCD symmetries as well [63, 65-70]. Originally, it was intended to describe the pion-nucleon interactions and the chiral degrees of freedom. A spinless scalar field σ a and triplet pseudoscalar fields π a are introduced in theory of quantized fields to the linear-sigma model, which is a low-energy effective model, in

which the generators T a = λ / a 2 with Gell-Mann matrices λ a and the real classical field forming an O (4) vector, /vector Φ = T a ( /vector σ , i/vector a π a ) are included. The chiral symmetry is explicitly broken by 3 × 3 matrix field H = T h a a , where h a are the external fields. Accordingly, under SU (2) L × SU (2) R chiral transformation, such as Φ → L + Φ R , the σ a sigma fields acquire finite vacuum expectation values, which in turn break SU (2) L × SU (2) R down to SU (2) L + R . These transformations produce massive sigma particle and nearly massless Goldstone bosons, the pions. Therefore, the constituent quarks gain masses, as well; m q = gf π , where g is coupling and f π is the pion decay constant. Also, the fermions can be introduced either as nucleons or as quarks. The σ fields under chiral transformations exhibit the same behaviour as that of the quark condensates and thus, σ can be taken as order parameters for the QCD chiral phase transition.

With the incorporation of the Polyakov-loop potential so that the Lagrangian of the PLSM reads L = L ¯ ψψ + L m -U ( φ, φ 0 , T ), where the first term stands for Lagrangian density of fermions with N c color degreesof-freedom, the second term gives the contributions of the mesonic fields, and finally the third term represents the Polyakov-loops potential incorporating the gluonic degrees-of-freedom and the dynamics of the quark-gluon interactions, i.e. deconfinement is also incorporated in this chiral model.

The questions which arise now are why PLSM reproduces well the non-perturbative lattice QCD simulations and why the PLSM agrees well with the thermal model calculations, especially for the freezeout boundary? The first question can be directly answered. PLSM incorporates both chiral and deconfinement QCD symmetries. On the other hand, it seems that both types of transitions are nearly coincident, especially at vanishing or small baryon chemical potentials. Within this region, both calculations are in excellent agreement with each others. At high temperatures, both chiral symmetry restoration and deconfinement transition produce almost free quarks and gluons, e.g. QGP. The reliability of the chiral effective model, PLSM, seems crucial, especially where lattice field theory is unavailable or the experimental results are not accessible yet.

The second question about the reasons why PLSM agrees well with the freezeout parameters deduced from the thermal model calculations can be answered as follows. First, at µ b &gt; ∼ 300 MeV, where lattice field theory likely suffers from the sign problem, it seems that both chiral and deconfinement boundaries become more and more distinguishable. It might be obvious that the critical temperature of the chiral phase transition would be smaller than that of the deconfinement transition, which in turn differs from the freezeout temperatures. Within these two limits, which should be subject of further studies, a temperature region is created, in which a phase of mixed hadron-QGP likely takes place. Last but not least, the T χ -µ b plane of the Polyakov linearsigma model [63] was determined under the condition of constant entropy density normalized to T 3 [21-23], i.e. likely manifesting the freezeout boundary. A future phenomenological study should be conducted in order to find out whether the condition of line-of-constant-physics gives results in agreement with the T χ -µ b plane for deconfinement.

Now it is in order to summarize some details about the two lattice QCD calculations (yellow and grey bands). [61, 62]. In ref. [61], the crossover boundary separating the hadron gas from the quark gluon plasma phases at small baryon chemical potentials was calculated using a four times stout smeared staggered fermion action with dynamical up, down, strange and charm quarks (2 + 1 + 1). This means that the two light quarks are degenerate. The masses of light quarks and that of the strange quark mass are tuned such that the physical pion and Kaon mass over the pion decay constant are reproduced for every lattice spacing. For the gauge action, the tree-level Symanzik improvement was used. In order to overcome the sign problem, imaginary values of the chemical potential have been considered, which are then translated through analytic continuation to the real values of the baryon-chemical potentials. The continuum extrapolation is based on lattices with 10, 12, and 16 temporal dimensions. The curvature of x = T c ( µ b ) /T c ( µ b = 0) was estimated according to anzaetze 1 + ax , 1 + ax + bx 2 , (1 + ax / ) (1 + bx )) and (1 + ax + bx ) -1 , where a and b are fit parameters. The corresponding critical temperature at µ b = 0 was estimated as 157 MeV, where vanishing strange density was assured.

In ref. [62], the critical temperatures of chiral crossovers at vanishing and finite values of baryon ( ), b strangeness ( s ), electric charge ( Q ), and isospin ( I ) chemical potentials obtained in the continuum limit from lattice QCD calculations carried out for 2 + 1 highly improved staggered quarks (HISQ) and the tree-level improved Symanzik gauge action with two degenerate up and down dynamical quarks and a dynamical strange quark, with physical quark masses corresponding to physical pion and Kaon masses are presented. The temporal extents are varied from N τ = 6, 8, 12, and 16, going towards progressively finer lattice spacing. The critical temperatures have been parameterized as T c ( µ x ) = T c ( µ x = 0) [ 1 -κ x 2 ( µ /T x c ( µ x )) 2 -κ x µ /T 4 ( x c ( µ x = 0)) 4 ] , where κ x 2 and κ x 4 are determined from Taylor expansions of chiral observables in µ x . The corresponding critical temperature at µ b = 0 was estimated as 156 5 . ± 1 5 MeV. At the chemical freeze-out of relativistic heavy-ion . collisions and under thermal conditions, such as, µ s ( T, µ b ) and µ Q ( T, µ b ) determined from strangeness neutrality and isospin imbalance, κ b 2 = 0 012 . ± 0 004 and . κ b 4 = 0 000 . ± 0 0004. .

Fig. 2: The same as in Fig. 1 but here a comparison between lattice QCD simulations [61, 62] and Polyakov linear-sigma model [63, 65-70] with the thermal model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00426v2-figure-14f489a9-4c86-46db-81d1-b1115b5735da.jpg)

## IV. CONCLUSIONS AND OUTLOOK

Among the various phases which take place in the strongly interacting matter under extreme conditions, we focused on the deconfinement and chemical freezeout boundaries. The authors compared results on T χ and µ b deduced from various heavy-ion experiments with recent lattice simulations, effective QCD-like Polyakov linear-sigma model, and equilibrium thermal model. Along the entire freezeout boundary, we conclude that an excellent agreement between the thermal model calculations and the experiments is found. Also, the estimations deduced from the Polyakov linear-sigma model excellently agree with the thermal model calculations. It should be noted that at low baryon density or high energies, both deconfinement and chemical freezeout boundaries are likely coincident. Accordingly, we can also conclude that the lattice calculations for the deconfinement transition agree well with the Polyakov linear-sigma model, where in both approaches QCD symmetries are included. At large baryon density or low energies, the two boundaries become distinguishable and probably form a phase in which hadrons and quark-gluon plasma likely coexist.

Based on the fact that Polyakov linear-sigma model agrees well with the lattice QCD simulations, at least within µ b -range of reliable simulations, a future phenomenological study should be conducted on Polyakov linear-sigma model to find out whether the condition of line-of-constant-physics gives results in agreement with the T χ -µ b plane for deconfinement. Furthermore, it intends to characterize the phase of mixed hadron-QGP and its possible predictions at the future facilities FAIR and NICA as well as its astrophysical consequences.

- [1] T. Banks and A. Ukawa, Nucl. Phys. B225 , 145 (1983).
- [2] A. N. Tawfik, Int. J. Mod. Phys. A29 , 1430021 (2014).
- [3] B. P. Abbott et al., Phys. Rev. Lett. 119 , 161101 (2017).
- [4] D. D. Ivanenko and D. F. Kurdgelaidze, Astrophysics 1 , 251 (1965).
- [5] M. Gyulassy and L. McLerran, Nucl. Phys. A750 , 30 (2005).
- [6] W. Busza, K. Rajagopal, and W. van der Schee, Ann. Rev. Nucl. Part. Sci. 68 , 339 (2018).
- [7] C. Greiner, P. Koch-Steinheimer, F. M. Liu, I. A. Shovkovy, and H. Stoecker, J. Phys. G31 , S725 (2005).
- [8] B. Mueller and A. Schaefer, arXiv:1712.03567 [nucl-th] (2017).
- [9] V. Vovchenko, M. I. Gorenstein, C. Greiner, and H. Stoecker, Phys. Rev. C99 , 045204 (2019).

- [10] L. D'Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, Adv. Phys. 65 , 239 (2016).
- [11] G. Fast, R. Hagedorn, and L. W. Jones, Nuovo Cim. 27 , 856 (1963).
- [12] G. Fast and R. Hagedorn, Nuovo Cim. 27 , 208 (1963).
- [13] F. Karsch, K. Redlich, and A. Tawfik, Eur. Phys. J. C29 , 549 (2003).
- [14] F. Karsch, K. Redlich, and A. Tawfik, Phys. Lett. B571 , 67 (2003).
- [15] K. Redlich, F. Karsch, and A. Tawfik, J. Phys. G30 , S1271 (2004).
- [16] A. Tawfik, Phys. Rev. D71 , 054502 (2005).
- [17] A. Tawfik and D. Toublan, Phys. Lett. B623 , 48 (2005).
- [18] M. Tanabashi et al., Phys. Rev. D98 , 030001 (2018).
- [19] J. Beringer et al., Phys. Rev. D86 , 010001 (2012).
- [20] A. N. Tawfik and E. Abbas, Phys. Part. Nucl. Lett. 12 , 521 (2015).
- [21] A. Tawfik, M. Y. El-Bakry, D. M. Habashy, M. T. Mohamed, and E. Abbas, Int. J. Mod. Phys. E25 , 1650018 (2016).
- [22] A. Tawfik, Nucl. Phys. A764 , 387 (2006).
- [23] A. Tawfik, Europhys. Lett. 75 , 420 (2006).
- [24] A. Tawfik, Nucl. Phys. A922 , 225 (2014).
- [25] A. Tawfik, Adv. High Energy Phys. 2013 , 574871 (2013).
- [26] A. Tawfik, Phys. Rev. C88 , 035203 (2013).
- [27] A. N. Tawfik, H. Yassin, and E. R. A. Elyazeed, Phys. Rev. D92 , 085002 (2015).
- [28] A. Tawfik, J. Phys. G31 , S1105 (2005).
- [29] G. Agakishiev et al., Physical Review C 84 , 014902 (2011).
- [30] V. V. Begun, V. Vovchenko, and M. I. Gorenstein, J. Phys. Conf. Ser. 779 , 012080 (2017).
- [31] M. Lorenz, Nucl. Phys. A931 , 785 (2014).
- [32] N. Abgrall et al., Eur. Phys. J. C74 , 2794 (2014).
- [33] G. Agakishiev et al., Phys. Rev. C90 , 015202 (2014).
- [34] G. Agakishiev et al., Phys. Rev. C92 , 024903 (2015).
- [35] T. Anticic et al., Phys. Rev. C80 , 034906 (2009).
- [36] T. Anticic et al., Physical Review C 79 , 044904 (2009).
- [37] V. Vovchenko, M. Gorenstein, L. Satarov, and H. St¨cker, Chemical freeze-out conditions in hadron resonance gas, o in New Horizons in Fundamental Physics , pages 127-137, Springer, 2017.
- [38] C. Alt et al., Phys. Rev. Lett. 94 , 192301 (2005).
- [39] I. Kraus, J. Phys. G31 , S147 (2005).
- [40] M. Floris, Nucl. Phys. A931 , 103 (2014).
- [41] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Nature 561 , 321 (2018).
- [42] B. Abelev et al., Physical Review Letters 111 , 222301 (2013).
- [43] D. Alexandre et al., Nuclear Physics A 931 , 1093 (2014).
- [44] B. Abelev et al., Physical Review C 91 , 024609 (2015).
- [45] A. Collaboration et al., arXiv preprint arXiv:1506.08453 (2015).
- [46] A. Collaboration et al., arXiv preprint arXiv:1710.07531 (2017).
- [47] B. Abelev et al., Physical Review Letters 109 , 252301 (2012).
- [48] J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, Physical Review C 73 , 034905 (2006).
- [49] C. Adler et al., Physical review letters 89 , 202301 (2002).
- [50] C. Adler et al., Physical Review Letters 87 , 112303 (2001).
- [51] J. Adams et al., Physics Letters B 637 , 161 (2006).
- [52] J. Velkovska, Nucl. Phys. A698 , 507 (2002).
- [53] J. Adams et al., Nucl. Phys. A757 , 102 (2005).
- [54] H. Oeschler, Lect. Notes Phys 516 , 1 (1999).
- [55] V. I. Yukalova and E. P. Yukalov, PoS Baldin-ISHEPP-XXI , 046 (2012).
- [56] K. A. Bugaev et al., Eur. Phys. J. A52 , 227 (2016).
- [57] P. Man Lo, M. Marczenko, K. Redlich, and C. Sasaki, Eur. Phys. J. A52 , 235 (2016).
- [58] J. Noronha-Hostler, Implications of Missing Resonances in Heavy Ions Collisions, in Workshop on Excited Hyperons in QCD Thermodynamics at Freeze-Out (YSTAR2016) Mini-Proceedings , pages 118-127, 2016.
- [59] S. Capstick and N. Isgur, Phys. Rev. D34 , 2809 (1986), [AIP Conf. Proc.132,267(1985)].
- [60] A. Bazavov et al., Phys. Rev. Lett. 113 , 072001 (2014).
- [61] R. Bellwied et al., Phys. Lett. B751 , 559 (2015).
- [62] A. Bazavov et al., Phys. Lett. B795 , 15 (2019).
- [63] A. Tawfik, N. Magdy, and A. Diab, Phys. Rev. C89 , 055210 (2014).
- [64] O. Philipsen, Lattice QCD at non-zero temperature and baryon density, in Modern perspectives in lattice QCD: Quantum field theory and high performance computing. Proceedings, International School, 93rd Session, Les Houches, France, August 3-28, 2009 , pages 273-330, 2010.
- [65] A. N. Tawfik and A. M. Diab, Phys. Rev. C91 , 015204 (2015).
- [66] A. M. Abdel Aal Diab and A. N. Tawfik, EPJ Web Conf. 177 , 09005 (2018).

- [67] A. Nasser Tawfik and A. Magied Diab, Int. J. Mod. Phys. A30 , 1550059 (2015).
- [68] A. N. Tawfik, A. M. Diab, and M. T. Hussein, J. Phys. G45 , 055008 (2018).
- [69] A. N. Tawfik, A. M. Diab, and M. T. Hussein, Int. J. Mod. Phys. A31 , 1650175 (2016).
- [70] A. N. Tawfik, A. M. Diab, and M. T. Hussein, Chin. Phys. C43 , 034103 (2019).
Abdel Nasser Tawfik Muhammad Maher A. H. El-Kateb Sara Abdelaziz
Deconfinement and freezeout boundaries in equilibrium thermal models
In different approaches, the temperature-baryon density plane of QCD matter
is studied for deconfinement and chemical freezeout boundaries. Results from
various heavy-ion experiments are compared with the recent lattice simulations,
the effective QCD-like Polyakov linear-sigma model, and the equilibrium thermal
models. Along the entire freezeout boundary, there is an excellent agreement
between the thermal model calculations and the experiments. Also, the thermal
model calculations agree well with the estimations deduced from the Polyakov
linear-sigma model (PLSM). At low baryonic density or high energies, both
deconfinement and chemical freezeout boundaries are likely coincident and
therefore the agreement with the lattice simulations becomes excellent as well,
while at large baryonic density, the two boundaries become distinguishable
forming a phase where hadrons and quark-gluon plasma likely coexist.