Threshold Collision Energy of the QCD Phase Diagram Tricritical Endpoint
Using the most advanced formulation of the hadron resonance gas model we analyze the two sets of irregularities found at chemical freeze-out of central nuclear-nuclear collisions at the center of mass energies 3.8-4.9 GeV and 7.6-9.2 GeV. In addition to previously reported irregularities at the collision energies 4.9 GeV and 9.2 GeV we found sharp peaks of baryonic charge density. Also we analyze the collision energy dependence of the modified Wroblewski factor and the strangeness suppression factor. Based on the thermostatic properties of the mixed phase of a 1-st order phase transition and the ones of the Hagedorn mass spectrum we explain, respectively, the reason of observed chemical equilibration of strangeness at the collision energy 4.9 GeV and above 8.7 GeV. It is argued that the both sets of irregularities possibly evidence for two phase transitions, namely, the 1-st order transition at lower energy range and the 2-nd order transition at higher one. In combination with a recent analysis of the light nuclei number fluctuations we conclude that the center of mass collision energy range 8.8-9.2 GeV may be in the nearest vicinity of the QCD tricritical endpoint. The properties of the phase existing between two phase transitions are revealed and discussed.
K. A. Bugaev111E-mail: firstname.lastname@example.org, R. Emaus, V. V. Sagun A. I. Ivanytskyi, L. V. Bravina,
D. B. Blaschke, E. G. Nikonov, A. V. Taranenko, E. E. Zabrodin,
G. M. Zinovjev
Bogolyubov Institute for Theoretical Physics, Metrologichna str. 14, Kiev 03680, Ukraine
Department of Physics, University of Oslo, PB 1048 Blindern, N-0316 Oslo, Norway
CENTRA, Instituto Superior Tcnico, Universidade de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal
Institute of Theoretical Physics, University of Wroclaw, pl. M. Borna 9, 50-204 Wroclaw, Poland
Bogoliubov Laboratory of Theoretical Physics, JINR Dubna, Joliot-Curie str. 6, 141980 Dubna, Russia
National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), Kashirskoe Shosse 31, 115409 Moscow, Russia
Laboratory for Information Technologies, JINR, Joliot-Curie str. 6, 141980 Dubna, Russia
Skobeltzyn Institute of Nuclear Physics, Moscow State University, 119899 Moscow, Russia
PACS: 25.75.q, 25.75.Nq
Theoretical and experimental searches for the (tri)critical endpoint ((3)CEP) of the quantum chromodynamics (QCD) phase diagram form one of the most important directions of physics of heavy ion collisions [1, 2]. Despite the claims that the QCD CEP parameters can be extracted from the RHIC data , the real situation is not so clear . In part, such a situation exists because the rigorous and reliable definition of finite volume analog of (3)CEP does not exist. In our opinion, therefore, the question is not whether the work  can justify the finite size scaling for the isothermal compressibility as it is assumed in  or for the isentropic compressibility as its is argued in , but whether there is sufficient experimentally established signals of the (3)CEP formation. Furthermore, we are sure that to safely locate the threshold collision energy of the (3)CEP one needs a simultaneous match of thermodynamical, statistical and hydrodynamical signals of its formation.
Recently the thermodynamical and hydrodynamical signals of the mixed quark-gluon-hadron phase formation were revealed [6, 7] from the analysis of different irregularities observed at chemical freeze-out (CFO) [8, 9, 10]. Moreover, such signals were observed at two ranges of the center-of-mass collision energies, namely GeV and GeV [6, 7, 10]. The thermodynamic signals include two sharp peaks of the trace anomaly (here , and denote, respectively, the energy density of the system, its pressure and temperature) observed at GeV and GeV, while the hydrodynamic signals consist of the 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 were found at the two ranges of the center-of-mass collision energies discussed above [6, 7]. Thermodynamical and hydrodynamical signals in the low collision energy region can be easily interpreted since they can be related with each other via the generalized shock adiabat model [11, 12] which shows that there is one-to-one correspondence between the peak of at CFO and the peak of at the shock adiabat [6, 7]. Moreover, at CFO the low energy signals are accompanied by the strong jumps of the pressure and the effective number of degrees of freedom [6, 7, 8]. In contrast to these findings, the high energy signals are less pronounced and, hence, their interpretation is not straightforward [6, 7, 10] and, therefore, one needs additional analysis of the collision energy range GeV.
In a recent work on the analysis of critical fluctuations of light nuclei  at CFO it was argued that at the collision energy GeV the QCD matter has strongest fluctuations and, hence, during its evolution this matter passes through the CEP. This work puts forward a strong argument in favor of the hypothesis that at the vicinity of the point GeV, i.e. almost at the location of the second peak of trace anomaly , there occurs another phase transformation. Moreover, the work  concludes that this is a vicinity of CEP, i.e. it has to be a 2-nd order phase transition. Thus, in addition to the thermodynamic and hydrodynamic signals of a possible phase transformation at the collision energy range GeV [6, 7, 10] the work , first of all, gives an additional and independent evidence in favor of such a transformation and, moreover, it provides us with the necessary statistical (fluctuation) signal. However, in combination with the previous findings on the signals of 1-st order phase transition at the collision energies GeV [6, 7, 10] one should conclude that these are the signals of two phase transitions and, therefore, at the vicinity of collision energies GeV there may exist not a CEP, but a 3CEP. It should be stressed that this is highly nontrivial conclusion, since neither the lattice QCD nor the QCD inspired field-theoretical models can presently tell us for sure whether QCD has the CEP or 3CEP.
Furthermore, if we observe the signals of two phase transitions there are two basic questions related to this fact: (I) can one determine which of these transitions is deconfinement and which one is the chiral symmetry restoration? and (II) what kind of phase can be probed at the collision energies GeV? The present work is an attempt to formulate the proper answers to these two major questions.
The work is organized as follows. In the next section we present the details of hadron resonance gas model and the fitting procedure of hadronic multiplicities. Section 3 is devoted to an analysis of new irregularities at CFO. In Section 4 we discuss the properties of explicit thermostat and apply them to explain the strangeness equilibration observed at the collision energy GeV. Also in this section we estimate the number of degrees of freedom of a phase possibly formed at GeV. The equation of state of (almost) massless particles with the relativistic treatment of hard-core repulsion is analyzed in Section 5. Our conclusions are summarized n Section 6.
Analysis of hadron multiplicities
First we perform the cross-check of the previously obtained results on the irregularities at CFO [6, 7, 10]. For this purpose we employ the newest version of the hadron resonance gas model (NHRGM) [14, 15, 16] which allows us to safely go beyond the usual Van der Waals approximation for the multicomponent case, i.e. for several different hard-core radii of hadrons. Also below we reveal the new irregularities and give their explanation based on the hypothesis of existence of two phase transitions.
The new key element of the NHRGM is an inclusion of the surface tension which is generated by the hard-core repulsion. Therefore, this equation of state is called the induced surface tension one. For the system with multicomponent hard-core repulsion in the grand canonical ensemble this equation of state is given by the system of two coupled equations
where the sums are running over all particles (and antiparticles) with the chemical potentials , the hard-core radii , the proper volumes and the proper surfaces . Here denotes the partial pressure of the point-like particles of sort with the degeneracy and the mass which in case of the Boltzmann statistics is
Here is the energy of particle with the 3-momentum and is the effective chemical potential.
In the grand canonical ensemble the one component Van der Waals equation of state (EoS) [17, 8, 18] can be obtained from the pressure of ideal gas, classical or quantum, by the modification of the chemical potential as . In the one component NHRGM the free energy associated with the excluded volume of particle is split up into a sum of the volume part , which is proportional to the pressure , and the surface part , which is proportional to the induced surface tension coefficient (here is the proper surface of particle having the hard-core radius ). Such a splitting allows one not only to account for the second virial coefficient, but also to reproduce the third and the fourth ones of the gas of hard-spheres of the radius with a good accuracy [14, 15]. The system (Analysis of hadron multiplicities) and (Analysis of hadron multiplicities) is a generalization of such an approach to a multicomponent case. The parameter switches between the excluded-volume and the proper-volume regimes. Numerical analysis  shows that for the system (Analysis of hadron multiplicities) and (Analysis of hadron multiplicities) correctly reproduces the one component  and multicomponent  versions of the well-known Carnahan-Starling EoS up to the packing fractions (here is the particle number density of hadron species ). Note that the Carnahan-Starling EoS and its multicomponent version  are well-known in the theory of simple liquids since they are able to reproduce the EoS of hard spheres up to its transition to a liquid state [22, 23].
The particle density of hadrons of sort ,
is expressed in terms of the auxiliary coefficients
At CFO the ratio of total hadronic multiplicities of sorts and is given by
This expression accounts for the hadrons which appear in strong decays of resonances with the branching ratios . The total chemical potential of hadron of the sort is and it depends on the baryonic chemical potential , the strange chemical potential , the isospin third projection chemical potential and the corresponding charges of this hadron. More details on the NHRGM can be found in .
Similarly to , here we fit the high quality experimental multiplicity ratios measured at AGS for GeV [24, 25, 26, 27, 28, 29, 30, 31, 32], the NA49 data measured at SPS energies GeV [33, 34, 35, 36, 37, 38, 39, 40, 41] and the STAR data measured at RHIC energies GeV . These experimental data allow us to construct 111 independent ratios measured at 14 collision energies [8, 18]. In order to verify the stability of the results on the new signals of phase transition between the hadron and quark gluon matters found in  with the help of the HRGM based on the Van der Waals approximation and to improve the description of (anti)-hyperons and (anti)protons we added 10 experimental data points of the ratio (on its exclusive role see a discussion in ). Although the ratio can be expressed via other ratios, i.e. it is dependent, but at AGS and SPS energies this ratio has rather small error bars compared to other ratios and this is important for the fit stability . Hence, in present work we analyze 121 ratios.
In addition, compared to works [10, 14, 15] we used a new set of global fitting parameters to minimize the mean square deviation , where the experimental ratios are denoted as , the theoretical ones as and the summation is carried out over all data points with the weights defined by the experimental error . The minimization of showed that compared to previous results [7, 10, 15] the values of local fitting parameters , , and are practically the same. We found the best description of the hadronic multiplicity ratios for the hard-core radius of pions fm, (anti)-hyperons fm, (anti)protons fm, other baryons fm and other mesons fm. This set of radii, i.e. the new radii afterwards, corresponds to a global minimum of which provides practically the same quality of fit as the set of global fitting parameters reported in [10, 15].
The set of old radii fitted in  provided and it included the hard-core radii of pions =0.15 fm, kaons =0.395 fm, -hyperons =0.085 fm, baryons =0.365 fm and mesons =0.42 fm. Comparing these two sets of hard-core radii and their fit quality one concludes that there are slight changes only which lead to little improvements. This can be seen from the ratios shown in Figs. 1 and 2. Therefore, our main conclusion is that the results previously obtained within the multicomponent Van der Waals EoS are almost unchanged and, hence, there is no need to revise the signals of phase transitions suggested earlier in [6, 7, 10].
Concerning the small values of the hard-core radii of pions and (anti) hyperons we would like to mention that these radii are an effective ones. It is possible that small value of the pion hard-core radius is a reflection of necessity to account for its relativistic nature at temperatures above 120 MeV [43, 44]. On the other hand the small hard-core radius of (anti) hyperons may reflect some subtleties of their interaction with hadronic medium. Hence, we hope that their values provide an important information for microscopic models of hadronic interaction which will have to explain our results.
New irregularities at CFO
Further results concerning the trace anomaly and baryon density are shown in Fig. 3. From the both panels of this figure one can see that each peak of the trace anomaly is accompanied by a strong peak of the total baryonic density . Moreover, the trace anomaly peak at GeV has a small amplitude, while its counterpart in the baryonic density has a large amplitude accompanied by small error bars making this peak quite pronounced.
where in the numerator and denote, respectively, the number of strange valence quarks and antiquarks in the hadron of sort , whereas and in the denominator denote, respectively, the number of and quarks in it (apparently, and play the same role for antiquarks). Note that the denominator in Eq. (New irregularities at CFO) differs from the traditional Wroblewski factor  because it accounts for a whole number of and valence quarks and antiquarks, and not only the ones which are paired to their valence antiquarks (quarks).
As one can see from Fig. 4 the factor demonstrates a jump right in the collision energy region which is associated with the mixed phase formation, i.e. for GeV, and in addition it shows a change of slope at GeV. Note that a similar behavior was found in  for the ratio and we confirm that in the present analysis it behaves similarly. The found collision energy dependence of further supports the hypothesis of Ref.  that such a behavior of and, hence, of is an indicator of two phase transformations. Indeed, the observed jump of these ratios is located in the collision energy range of the mixed phase formation (i.e. it is a 1-st order phase transition), then a change of their slope at GeV can be naturally associated with a weak 1-st order or a 2-nd order phase transition.
The collision energy behavior of the strangeness suppression factor  shown in the right panel of Fig. 4 demonstrates more complicated dependence. Here one can see two peaks for GeV which evidence about the strangeness enhancement (over-saturation) and the strangeness equilibrium dale at GeV. At and above GeV one observes a chemical equilibration of strange charge, since . One might think that compared to the peak at GeV the smaller one seen at GeV is not statistically significant. However, a thorough analysis of two alternative approaches to describe the chemical non-equilibrium of strangeness performed in  shows that the peak in the vicinity of GeV exists practically unmodified (see the upper panel of Fig. 2 in ). At the same time the regime of reaching the strangeness chemical equilibrium, i.e. , at GeV is also observed in .
Based on this discussion we can further refine a hypothesis on the cause of chemical equilibrium of strange charge in central nuclear collisions at the collision energy GeV formulated in . In our opinion, the most natural explanation of this phenomenon is that at GeV there appeared the quark-gluon bags with the exponential mass spectrum proposed by R. Hagedorn . As it was predicted in  and shown numerically in [50, 51, 52], the exponential mass spectrum acts as a perfect thermostat and a perfect particle reservoir. In other words, all particles which appear from such bags at their hadronization will be born in a state of full thermal and chemical equilibrium [48, 50, 51, 52].
Furthermore, an existence of two peaks of strongly enhances the hypothesis about existence of two phase transformations. Indeed, since the right slope of the higher peak contains the region of the 1-st order phase transition at GeV, then it is naturally to assume that the right slope of the lower peak also contains the phase transition region. Moreover, from Fig. 4 one can see that the end of the 1-st order phase transition is shifted on 1.1 GeV to a higher collision energy than the peak existing at GeV. Then in accordance with our hypothesis the end of the other phase transition should be also shifted on about 1.1 GeV to a higher collision energy than the peak at GeV, i.e. one would expect it at GeV. Thus, again we independently arrive to the same conclusion as from the analysis of the trace anomaly and the baryonic density peaks that the vicinity of GeV corresponds to a phase transformation. Since this conclusion is based on the assumption that the 1-st order phase transition can be reached at GeV, then it is logically to expect that QCD has a 3CEP and it is located close to the collision energy range GeV. Thus, quite independently to Ref.  we came to a similar conclusion about the endpoint, but, in contrast to , we consider it as the 3CEP.
Note that the collision energy GeV as the onset of deconfined phase appears not only in Ref. . Thus, in Ref.  it was shown that the horn in ratio can be naturally explained, if one assumes that the partonic phase, or the quark gluon plasma (QGP) production threshold is GeV. Note that we came to a similar conclusion above on the basis of the Hagedorn thermostat hypothesis. But now we are facing a question, if the QGP is formed in heavy ion collisions at GeV, then what kind of phase transition occurs at lower collision energies GeV? From the discussion above it is evident that such a phase cannot consist from the bags of QGP which are filled with the quarks and gluons, otherwise one would find .
Strangeness equilibrium dale and mixed phase as an explicit thermostat
Before discussing the possible interpretations of the above results, we would like to reinforce our arguments about the mixed phase formation at GeV. Using the thermostatic properties of the mixed phase of the 1-st order phase transition  below we would like to explain the appearance of the strangeness equilibrium dale, i.e. , at GeV. In  the examples of the explicit thermostat and explicit particle reservoir were discussed and it was argued that under the constant pressure condition the mixed phase of the 1-st order phase transition, i.e. two pure phases being in a thermal and chemical equilibrium with each other, represent both thermostat and particle reservoir as long as it has enough energy resources to keep a constant temperature. In other words, under the constant pressure condition an explicit thermostat keeps a constant temperature despite transmitting out (in) some amount of heat. The latter only changes the volume fractions of two phases: the phase with higher heat capacity (for definiteness, a liquid) condenses a certain amount of gas under external cooling or it partly evaporates into a gas under external heating. Apparently, for finite systems the amount of imparted heat is also finite and depends on the masses of both phases and their heat capacities. Similarly, one can add or remove some amount of each phase, but under the constant pressure condition the system will continue to keep a constant temperature and a full chemical equilibrium, i.e. the equal values of chemical potentials for both pure phases. For our discussion it is important that up to some maximal value any amount of the removed phase, namely the gas of hadrons, for definiteness, will be, by definition, in a full chemical and thermal equilibrium with the mixed phase. Again the maximal amount of removed phase depends on the masses of both phases and their energy densities. Note that to justify such a picture we have to assume that the mixed phase has sufficiently large volume, so the finite size effects are not significant.
Above we discussed a static picture, but the question is whether the condition of constant pressure can be achieved in central nuclear collisions at GeV. To extend this picture to the case of central collisions of heavy ions we refer to the compressional shock model of these collisions [54, 55, 56, 57, 58, 59] and its generalization to the mixed phase with the anomalous thermodynamic properties [11, 12, 6, 7]. This model of the collision process, which neglects the nuclear transparency, can be reasonably well justified at intermediate collision energies per nucleon A GeV 20 A GeV [55, 56, 57, 58, 59]. Moreover, the direct comparison of the results of the compressional shock model with those of the three-fluid model  made in Ref.  demonstrates once more that at laboratory energies up to 30 A GeV this model can be used for quantitative estimates, while at higher energies (up to about 50 A GeV in laboratory frame) it provides a qualitative description only.
The compressional shock model of central nuclear collisions allows one to determine the initial conditions for the subsequent hydrodynamic evolution. The latter are given by the shock adiabat , each point of which corresponds to a certain energy of collision (more details can be found in [59, 11]). The shock adiabat OABC (see Fig. 5) which was found in [6, 7] shows that in the mixed phase (the segment AAB in Fig. 5) of the 1-st order phase transition the pressure at the central longitudinal rapidity region practically does not depend on the collision energy. Note, however, that the shock transitions into the states belonging to the segment ABC for the shock adiabat OABC of Fig. 5 are mechanically unstable [62, 63] due to anomalous thermodynamic properties of the mixed phase (the segment AB) and in this case more complicated flow patters should appear. The latter for the relativistic shock waves were determined in [11, 12] and were calculated in [6, 7]. In particular, the shock adiabat region of mixed phase AB should be replaced by the the Poisson adiabat  AD (see the right panel of Fig. 5; here and below the variable denotes the generalized specific volume ). It is remarkable that from the point A to the point D on the Poisson adiabat the pressure changes from about 440 MeVfm to about 450 MeVfm, i.e. the relative change of pressure is about or 2.27%. The flow configuration which corresponds to the segment AD consists of the compressional shock wave and compressional simple wave moving outwards (see Fig. 6). The shock wave describes the transition from the cold nuclear matter state (the point O in Fig. 5) to the state , i.e. the tangent point A to the shock adiabat, while the compressional simple wave describes the transformation from the state to any state on the considered segment of the Poisson adiabat with and (see Refs. [11, 12, 6, 7] for more details). In other words, by increasing the collision energy, at the central rapidity region of the collision one creates the states with higher energy density, but, practically, with the constant pressure. Also due to the fact that the energy density of created dense phase (the phase of massless particles or (PMP) hereafter) is about 10 times higher than the one of the hadron gas [6, 7], the major part of energy of a matter formed after the compressional shocks between the states and disappear is concentrated right at the central rapidity region of collision. Therefore, just the mixed phase states with practically the same pressure and the same value of entropy per baryon in each 3-dimensional point of the formed matter should define the pattern of its hydrodynamic expansion. Apparently, the hadronic matter which will be born from the mixed phase will be in a full thermal and chemical equilibrium with it. Its further hydrodynamic evolution is out of the scope of the present work, but according to the contemporary paradigm the hadronic matter, including the strange hadrons, born off the mixed phase will remain in the full equilibrium till the moment of CFO.
If the collision energy increases above the point D on the Poisson adiabat, then the matter formed at the central rapidity region of collision will correspond to the segment and the part of mixed phase will gradually decrease vanishing completely at the point C of the shock adiabat. Simultaneously, the role of mixed phase as the explicit thermostat and particle reservoir will be gradually diminished. Counting the discussion above, we can explain the dale seen in the right panel of Fig. 4 by the mixed phase formation at the collision energies GeV and its steady disappearance at GeV together with the corresponding appearance of chemical nonequilibrium of strangeness. Furthermore, such a hypothesis allows one to simultaneously understand the strong increase of the modified Wroblewski factor within the collision energy interval GeV and the nontrivial collision energy dependence of the factor at these energies (compare the both panels of Fig. 4).
Although it is expected that in QCD there may exist two phase transitions, i.e. deconfinement and chiral symmetry restoration, it is hard to unambiguously identify the nature of phase transition at GeV, since presently there is no reliable information from lattice QCD at high baryonic densities. In Refs. [6, 7, 10] it was concluded that at these collision energies there may exists the QGP which is probed by the lattice QCD approach. According to our present hypothesis in heavy ion collisions the two phase transitions may be observed. From the discussion above it is clear that at the collision energies GeV there may exist transition to the QGP, i.e. the phase of matter studied by the lattice QCD, and, hence, at the collision energies GeV the transition to QGP is impossible. Thus, the latter collision energy range may correspond to another phase which we call the phase of massless particles (PMP). Although at present the amount of experimental data at the collision energy range GeV is rather limited, we can make an educated guess about its properties. Our guess is based on the EoS of the PMP which was determined in Refs. [6, 7] from fitting the entropy per baryon along the shock adiabat [11, 12]. This EoS is similar to the MIT-Bag model 
but the constants , , , and are rather different from what is predicted by the perturbative QCD for massless gluons and (anti)quarks. Despite this fact, an additional condition that the pseudocritical temperature value at zero baryonic density is about 150 MeV was used in the fitting procedure of Ref. . Note that the value MeV is in agreement with the lattice QCD data .
In Ref.  the difference of coefficients , and from the ones of perturbative QCD , and was interpreted as the and dependence of non-perturbative pressure at high baryonic densities. Then the pressure (Strangeness equilibrium dale and mixed phase as an explicit thermostat) can be rewritten as . This was a way to avoid the contradiction with perturbative QCD. However, now, in accordance with the hypothesis of two phase transitions, one can consider Eq. (Strangeness equilibrium dale and mixed phase as an explicit thermostat) differently, namely as a source of information about the PMP which is formed at the collision energies 4.9 GeV GeV.
Recall that the EoS for massless bosons and fermions (with the chemical potential ) has the form 
where the number of the massless degrees of freedom is given in terms of degeneracy factors of massless bosons and massless (anti)fermions . Comparing the coefficients of powers of temperature in Eqs. (Strangeness equilibrium dale and mixed phase as an explicit thermostat) and (Strangeness equilibrium dale and mixed phase as an explicit thermostat) one can estimate the number of the massless degrees of freedom of this phase from the coefficient . Then one obtains the number of massless boson and fermion degrees of freedom , where we used the Planck constant in order to get a dimensionless value of . Note that is essentially larger than the number of massless degrees of freedom of QCD (the total degeneracy of gluons, quarks and antiquarks), but it is comparable to the total number of spin-isospin configurations of known hadronic states . Assuming that fermions have the baryonic charge , one can use the coefficient to estimate the number of massless fermionic and antifermionic degrees of freedom , then one finds that it is still far larger than the expected number of quarks and antiquarks. Also one can use the coefficient to get the number of massless fermionic and antifermionic degrees of freedom which is inconsistent with found from the coefficient . Clearly, it would be a great surprise, if the coefficients and of Eq. (Strangeness equilibrium dale and mixed phase as an explicit thermostat) would completely agree with the corresponding coefficients in Eq. (Strangeness equilibrium dale and mixed phase as an explicit thermostat). Hence, we conclude that interaction between the constituents of PMP may differ from the simple bag model parameterization. Nevertheless, our hypothesis is that the PMP consist of non-strange hadronic states with the mass which is much smaller than , i.e. . In this case the PMP can be considered as the hadronic state with the restored unitary symmetry. Our hypothesis on the PMP existence is strongly supported by the results of the microscopic model Parton-Hadron-String Dynamics (PHSD) [68, 69] which show that at the vicinity of the collision energy GeV the effects of baryonic and partonic mass reduction, i.e. the effects of chiral symmetry restoration in the hadronic phase, are very strong.
The question about existence of almost massless strange hadrons in the PMP cannot be answered at the moment with high confidence, i.e. now it is unclear whether the unitary symmetry in PMP is restored completely or in the non-strange sector only. However, looking at the peaks at GeV and GeV in the right panel of Fig. 4 one may guess that the same mechanism of the strangeness enhancement may be responsible for them, while this mechanism disappears at GeV due to the formation of quark-gluon bags with the Hagedorn mass spectrum. Therefore, if in the hadronic phase, i.e. at GeV the strange hadrons are massive, then we believe it is reasonable to assume that in the PMP they are massive too. Although the PHSD microscopic model predicts an essentially weaker reduction of the constituent strange quark compared to the mass of constituent u- and d-quarks , we believe that the questions of whether and how the masses of strange hadrons in the PMP are reduced compared to the vacuum values, will have to be answered experimentally on the accelerators of new generation, i.e. NICA JINR and FAIR GSI.
It is interesting that a long time ago an existence of the hadronic phase with abundant number of baryon-antibaryon pairs was discussed in Ref. . In this model the non-linear interaction of scalar meson provides an essential reduction of baryonic masses and generates an additional solution to the usual hadronic matter which in Ref.  is called as the abnormal matter. Similar additional solutions for the self-interacting mesonic fields are also discussed in [71, 72, 73]. It is intriguing, however, that some thermodynamic parameters of the normal/abnormal matter are close to our findings. For example, on the isotherm MeV shown in Fig. 3 of  the maximal pressure of the usual hadronic matter (normal state according to Ref. ) is about 420 MeVfm. Just compare it with the pressure of point A on the boundary of hadronic and mixed phases in Fig. 5. On the other hand, from Fig. 3 of  one can see that along the isotherm MeV the abnormal matter can acquire higher values of pressure which are comparable to the ones shown in Fig. 5 of this work for the compressional and generalized shock adiabats.
Of course, the existence of PMP with nearly massless non-strange hadrons is not the only possibility. Another possibility could be a presence of tetraquark condensate, which may generate two 1-st order chiral phase transitions . One more possibility could be an existence of the quarkyonic matter , which has an appropriate collision energy range. Hopefully, the nature and properties of the PMP will be studied at the accelerators NICA JINR and FAIR GSI. However, despite of its origin, if at high temperatures one takes the EoS (Strangeness equilibrium dale and mixed phase as an explicit thermostat) literally, then one apparently faces a contradiction with the lattice QCD thermodynamics. Indeed, if one fixes and increases , then despite of the large value of constant pressure in Eq. (Strangeness equilibrium dale and mixed phase as an explicit thermostat) at high temperatures the pressure of PMP will exceed the pressure of massless (anti)quarks and gluons. The first possibility to avoid such a contradiction is to claim that the PMP pressure (Strangeness equilibrium dale and mixed phase as an explicit thermostat) can be used just in the vicinity of the phase transition, since it does not correctly account for interaction of constituents at higher energy density. Basically we agree with this statement, but, in addition, we consider a possibility to get rid of the discussed contradiction with lattice QCD on a physical ground.
Phase of massless particles with relativistic hard-core repulsion
If, according to our hypothesis, the PMP consists of hadrons, then one cannot ignore their hard-core repulsion. Nowadays there are many articles discussing various applications of the hard-core repulsion in high energy nuclear physics [8, 9, 14, 15, 16, 17, 18, 76, 77, 78, 43, 44], however, the rigorous approach to relativistic treatment of Lorentz contraction of rigid spheres was developed in Refs. [43, 44]. The ultra-relativistic limit of Lorentz contracted rigid spheres analyzed in  shows that in contrast to the non-relativistic case the EoS with relativistic Van der Waals hard-core repulsion is causal at high pressures, i.e. its speed of sound does not exceed the speed of light.
Our suggestion to avoid a contradiction with the lattice QCD thermodynamics is to account for the Lorentz contraction of rigid spheres in the PMP pressure (Strangeness equilibrium dale and mixed phase as an explicit thermostat) in the spirit of Ref. . For simplicity we consider the Boltzmann limit of a one component gas, i.e. for equal masses of bosons and (anti)fermions and the same hard-core radius . Note that such a simplification is in line with the main property of PMP that particle masses are negligible compared to system temperature and simultaneously it provides us with a qualitatively correct solution of this complicated task. The pressure of the Lorentz contracted rigid spheres having the temperature , the chemical potential and the relativistic excluded volume per particle is given by a solution of the system 
where is the mass of hadron, is its degeneracy factor, denotes the thermal density of particles per single degree of freedom, is relativistic energy of particle having the 3-momentum and is its baryonic chemical potential. In the ultra-relativistic limit the relativistic excluded volume per particle can be cast as 
where the non-relativistic proper volume of particle with the hard-core radius is . Eq. (Phase of massless particles with relativistic hard-core repulsion) is valid for ; to use it for one has to make a replacement in (Phase of massless particles with relativistic hard-core repulsion). The coordinate system is chosen in such a way that the angle between the 3-vectors of particles’ momenta and coincides with the usual spherical angle of spherical coordinates (for more details see Appendix A in ). For definiteness, the OZ-axis of the momentum space coordinates of the second particle is chosen to coincide with the 3-vector of the momentum of the first particle.
It is evident that the volume strongly depends not only on particle momenta, but on the angle between them. Due to this property at high pressure the contribution of configurations with larger volume are suppressed and, hence, the main contribution to the double integral in Eq. (Phase of massless particles with relativistic hard-core repulsion) is given by the configurations with the minimal value of relativistic excluded volume. It is remarkable that with an accuracy of about 7% the ultra-relativistic expression (Phase of massless particles with relativistic hard-core repulsion) recovers the excluded volume of two non-relativistic rigid spheres  and, therefore, it can be safely used at low particle densities too.
Using the pressure of Lorentz contracted rigid spheres (Phase of massless particles with relativistic hard-core repulsion), for , one can generalize EoS (Strangeness equilibrium dale and mixed phase as an explicit thermostat) to high energy densities as
The first term on the right hand side of Eq. ( ‣ Phase of massless particles with relativistic hard-core repulsion) corresponds to the Boltzmann gas of chargeless particles, while the second and third ones correspond to the Boltzmann gas of particles with charges +1 and -1, respectively. Expanding the exponents which are staying on the right hand side of ( ‣ Phase of massless particles with relativistic hard-core repulsion) up to terms , one can readily check that for and the right hand side of ( ‣ Phase of massless particles with relativistic hard-core repulsion) recovers the one of Eq. (Strangeness equilibrium dale and mixed phase as an explicit thermostat). The term on the right hand side of ( ‣ Phase of massless particles with relativistic hard-core repulsion) “compensates” an extra contribution coming from the expansion of fermionic exponents . It is easy to verify that for massless particles and a zero value of hard-core radius the and dependent terms of pressures (Strangeness equilibrium dale and mixed phase as an explicit thermostat) and ( ‣ Phase of massless particles with relativistic hard-core repulsion) differ from each other by less than 10% for and any .
Apparently, Eq. ( ‣ Phase of massless particles with relativistic hard-core repulsion) is not a unique way to extrapolate the PMP pressure (Strangeness equilibrium dale and mixed phase as an explicit thermostat) to high energy densities, but, besides the fact that the pressure ( ‣ Phase of massless particles with relativistic hard-core repulsion) is one of the simplest possibilities, one can also check that for MeV and it has only a single positive solution for any value of . Therefore, we believe that Eq. ( ‣ Phase of massless particles with relativistic hard-core repulsion) is well suited for demonstrating how the Lorentz contracted relativistic excluded volume can modify the PMP pressure at high energy densities. A sufficient condition for an existence of a single positive solution of ( ‣ Phase of massless particles with relativistic hard-core repulsion) can be found as follows. For the fixed values of and with the restriction one can solve Eq. ( ‣ Phase of massless particles with relativistic hard-core repulsion) graphically. Noting that the right hand side of ( ‣ Phase of massless particles with relativistic hard-core repulsion) is a monotonously decreasing function of the variable due to excluded volume effect, one can conclude that such a function can have a single intersection with the straight line , i.e. with the left hand side of ( ‣ Phase of massless particles with relativistic hard-core repulsion), if at the point the inequality
Assuming that the condition (Phase of massless particles with relativistic hard-core repulsion) is obeyed, we consider the asymptotic behavior of the pressure ( ‣ Phase of massless particles with relativistic hard-core repulsion) for two cases, namely high limit and high limit. Apparently, both of them correspond to a high pressure limit and, therefore, the bag pressure and the last term on the right hand side of Eq. ( ‣ Phase of massless particles with relativistic hard-core repulsion) can be neglected. Applying the results of Ref.  to Eqs. (Phase of massless particles with relativistic hard-core repulsion), one can easily show that in the limit the dimensionless variable . Using this result for the pressure ( ‣ Phase of massless particles with relativistic hard-core repulsion), in the limit and finite values of one can find
where is a slowly decreasing dimensionless function of hadronic mass . The effective number of degrees of freedom now is . Moreover, the temperature dependence of pressure ( ‣ Phase of massless particles with relativistic hard-core repulsion) is essentially weaker than the one for massless quarks and gluons (see Eq. (Strangeness equilibrium dale and mixed phase as an explicit thermostat)).
In the limit of fixed and one finds for Eqs. (Phase of massless particles with relativistic hard-core repulsion) and, similarly, for ( ‣ Phase of massless particles with relativistic hard-core repulsion) that 
the pressure of relativistic hard-spheres does not depend on the number of degrees of freedom at all and its dependence is again essentially weaker than dependence of massless (anti)quarks. Therefore, these two examples show that the Lorentz contraction of relativistic hard-spheres, in principle, is able to resolve the problem of high energy density extrapolation of PMP EoS (Strangeness equilibrium dale and mixed phase as an explicit thermostat) in such a way that it does not contradict to QCD phenomenology.
In this work we perform a thorough inspection of the irregularities of thermodynamic quantities observed at CFO and discuss the possible signals of two QCD phase transitions. The most remarkable irregularities include two sets of correlated quasi-plateaus found in [6, 7] which are located at the collision energy ranges GeV and GeV, and two peaks of trace anomaly observed at the maximal energy of each set of quasi-plateaus. Using the most advanced version of the hadron resonance gas model which allows one to safely go beyond the standard Van der Waals approximation, here we found two strong peaks of the baryonic charge density located exactly at the collision energies of the trace anomaly peaks, i.e. at GeV and GeV. The usage of the NHRGM which is not restricted by the Van der Waals approximation was necessary, since at the mentioned peaks the baryonic charge density is high and corresponds to values of the normal nuclear density fm. Moreover, during the fitting even higher values of baryonic charge density were analyzed.
In addition in this work we closely studied the collision energy dependence of the ratio of total number of strange quarks and antiquarks to the number of non-strange quarks and antiquarks (a modified Wroblewski factor ) and the strangeness suppression factor . With the help of NHRGM we find a sizable jump of the factor, when the collision energy increases from GeV to GeV, while at the collision energies above GeV we observe a saturation of the factor. Similarly to previous findings [8, 9, 15] we observe rather complicated energy dependence of factor: at GeV and GeV there exist two maxima, while at GeV and at GeV the strangeness seems to reach a chemical equilibrium.
Since the low energy set of quasi-plateaus along with the trace anomaly peak at GeV are unambiguously explained by the anomalous thermodynamic properties of mixed phase of the 1-st order phase transition formed at GeV [6, 7], we conclude that other irregularities observed at GeV, i.e. the peak of baryonic charge density, the factor jump and the strangeness equilibration can be the other signals of this phase transition. In this work we give a straightforward explanation to the strangeness equilibration at GeV and at GeV. Based on the unique thermostatic properties of Hagedorn mass spectrum we follow Ref.  and suggest that at GeV the formation of quark-gluon bags with the exponential mass spectrum forces the hadrons born at the moment of hadronization to be in a full thermal and chemical equilibrium with the bags. Similarly, in this work we notice that under the condition of the constant external pressure any mixed phase of the 1-st order phase transition is an explicit thermostat and explicit particle reservoir and, hence, anything born out of it should be in a full thermal and chemical equilibrium with the mixed phase. Since the mixed phase EoS employed in Refs. [6, 7], indeed, leads to an approximate collision energy independence of initial pressure of a system formed in the collision process, we came to a natural conclusion that the formation of mixed phase at GeV is responsible for the chemical equilibration of strangeness at this energy. Thus, the explicit thermostatic properties of the mixed phase cause an appearance of the strangeness equilibrium dale at GeV. On the other hand, due to an absence of an alternative explanation of such a dale, its appearance, in turn, is an independent and strong argument in favor of the mixed phase formation at this collision energy.
As it is argued here and in Refs. [6, 7] the above mentioned irregularities observed at the collision energy range GeV are the signals of the 1-st order phase transition, then, according to our hypothesis, a set of similar irregularities found at the collision energies GeV should be associated with another phase transition. Such a hypothesis is supported by the detailed analysis of ratio of Ref.  and by the one of light nuclei fluctuations of Ref. . In addition to thermodynamic and hydrodynamic signals of phase transition discussed here the work  provides us with the fluctuation signal of the 2-nd order phase transition at the vicinity of GeV. Therefore, we argue that the collision energy range GeV may be the nearest vicinity of 3CEP.
The EoS (Strangeness equilibrium dale and mixed phase as an explicit thermostat) used in [6, 7] to model the properties of high density matter formed at the initial stage of collision for the energies GeV provides us with a unique opportunity to reveal the total number of massless bosonic and fermionic degrees of freedom. The found number is about 1770, which is too large for elementary degrees of freedom of QCD, but it has the same order of magnitude as the total number of experimentally known hadronic spin-isospin states. Although, at the moment there is insufficient data to discuss a possible origin of this phase, our hypothesis that such a phase corresponds to a gas of interacting, but almost massless hadrons, is strongly supported by the results of PHSD microscopic model which accounts for the effects of hadronic and partonic mass reduction in a very similar energy range [68, 69]. Therefore, we explicitly show that, if one accounts for the excluded volume of these hadronic degrees of freedom in a relativistic fashion, then at high energy densities the EoS (Strangeness equilibrium dale and mixed phase as an explicit thermostat) can be modified in such a way that it does not contradict to the QCD thermodynamics.
The above results and hypotheses are rather speculative, but they reflect the present state of art of heavy ion collision phenomenology, where at the moment there is insufficient experimental and theoretical information right at the most promising regions of the QCD phase diagram. Therefore, we hope that this work will stimulate both theoreticians and experimentalists to study the discussed problems in more details in order to prove or disprove our results.
Acknowledgments. The authors appreciate the valuable comments of E. E. Kolomeitsev, S. N. Nedelko and A. I. Titov. K.A.B., A.I.I., V.V.S. and G.M.Z. acknowledge a partial support from the program “Nuclear matter under extreme conditions” launched by the Section of Nuclear Physics of the National Academy of Sciences of Ukraine. The work of K.A.B. and L.V.B. was performed in the framework of COST Action CA15213 “Theory of hot matter and relativistic heavy-ion collisions” (THOR). K.A.B. is thankful to the COST Action CA15213 for a partial support. V.V.S. acknowledges a partial support by grants from “Fundação para a Ciência e Tecnologia”. The work of D.B. was supported in part by the Polish National Centre (NCN) under contract number UMO-2011/02/A/ST2/00306. The work of A.V.T. was partially supported by the Ministry of Science and Education of the Russian Federation, grant No 3.3380.2017/4.6, and by National Research Nuclear University “MEPhI” in the framework of the Russian Academic Excellence Project (contract No 02.a03.21.0005, 27.08.2013).
- 1. Stephanov M., Rajagopal K. and Shuryak E. Event-by-event fluctuations in heavy ion collisions and the QCD critical point // Phys. Rev. D —- 1999. —- V. 60, no. 114028. —- P. 1–32 and references therein.
- 2. NICA White Paper // —- 2013. —- V. 9.01 and references therein.
- 3. Lacey R. Indications for a Critical End Point in the Phase Diagram for Hot and Dense Nuclear Matter // Phys. Rev. Lett. —- 2015. —- V. 114, no. 14. —- P. 142301 and references therein.
- 4. Antoniou N.G., Davis N. and Diakonos F.K., What is the gamma gamma resonance at 750 GeV? —- 2016. —- arXiv:1607.01326v1 [nucl-th] see a discussion therein.
- 5. Ladrem M. and Ait-El-Djoudi A., Finite-size effects and scaling for the thermal QCD deconfinementphase transition within the exact color-singlet partition function // Eur. Phys. J. C —- 2005. —- V. 44, no. 2. —- P. 257â265.
- 6. Bugaev K.A., Ivanytskyi A.I., Oliinychenko D.R., Sagun V.V., Mishustin I.N., Rischke D.H., Satarov L.M. and Zinovjev G.M. Thermodynamically Anomalous Regions As A Mixed Phase Signal // Phys. Part. Nucl. Lett. —- 2015. —- V. 12, no. 2 —- P. 238â245, arXiv:1405.3575.
- 7. Bugaev K.A., Ivanytskyi A.I., Oliinychenko D.R., Sagun V.V., Mishustin I.N., Rischke D.H., Satarov L.M. and Zinovjev G.M. Thermodynamically anomalous regions and possible new signals of mixed-phase formation // Eur. Phys. J. A —- 2016, —- V. 52, no. 6, —- P. 175, arXiv:1412.0718.
- 8. Bugaev K.A., Oliinychenko D.R., Cleymans J., Ivanytskyi A.I., Mishustin I.N., Nikonov E.G. and Sagun V.V. Chemical Freeze-out of Strange Particles and Possible Root of Strangeness Suppression // Europhys. Lett. —- 2013, —- V. 104,—- P. 22002 and references therein.
- 9. Sagun V.V. -Anomaly in the Hadronic Chemical Freeze-out // Ukr. J. Phys. —- 2014, —- V. 59, no. 8,—- P. 755-763.
- 10. Bugaev K.A., Sagun V.V., Ivanytskyi A.I., Oliinychenko D.R., Ilgenfritz E.-M., Nikonov E.G., Taranenko A.V. and Zinovjev G.M. New Signals of Quark-Gluon-Hadron Mixed Phase Formation // Eur. Phys. J. A —- 2016, —- V. 52, —- P. 227, arXiv:1510.03099.
- 11. Bugaev K.A., Gorenstein M.I., Kämpfer B.,Zhdanov V.I. Generalized shock adiabatics and relativistic nuclear collisions // Phys. Rev. D —- 1989, —- V. 40, —- P. 2903.
- 12. Bugaev K.A., Gorenstein M.I., Rischke D.H. Pion multiplicity in heavy-ion collisions: possible signature of the deconfinement transition // Phys. Lett. B —- 1991, —- V. 255, —- P. 18.
- 13. Kai-Jia Sun, Lie-Wen Chen, Che Ming Ko and Zhangbu Xu Probing QCD critical fluctuations from light nuclei production in relativistic heavy-ion collisions —- 2017. —- arXiv:1702.07620v1 [nucl-th].
- 14. Bugaev K.A., Sagun V.V., Ivanytskyi A.I., Yakimenko I.P., Nikonov E.G., Taranenko A.V., Zinovjev G.M. Hadron Resonance Gas Model for An Arbitrarily Large Number of Different Hard-Core Radii —- 2016. —- arXiv:1611.07349v2 [nucl-th].
- 15. Sagun V.V., Bugaev K.A., Ivanytskyi A.I., Yakimenko I.P., Nikonov E.G., Taranenko A.V., Greiner C., Blaschke D.B., Zinovjev G.M. Hadron Resonance Gas Model with Induced Surface Tension —- 2017. —- arXiv:1703.00049 [nucl-th].
- 16. Bugaev K.A., Ivanytskyi A.I., Sagun V.V., Nikonov E.G., Zinovjev G.M. Equation of State of Quantum Gases Beyond the Van der Waals Approximation —- 2017. —- arXiv:1704.06846 [nucl-th].
- 17. Andronic A., Braun-Munzinger P. and Stachel J. Hadron production in central nucleusânucleus collisions at chemical freeze-out // Nucl. Phys. A —- 2006, —- V. 722, —- P. 167.
- 18. Bugaev K.A., Oliinychenko D.R., Sorin A.S., Zinovjev G.M. Simple solution to the Strangeness Horn description puzzle // Eur. Phys. J. A —- 2013, —- V. 49, —- P. 30.
- 19. Rischke D.H., Gorenstein M.I., Stöcker H. and Greiner W. Excluded volume effect for the nuclear matter equation of state // Z. Phys. C —- 1991, —- V. 51, —- P. 485.
- 20. Carnahan N.F., Starlin K.E. Equation of State for Nonattracting Rigid Spheres // J. Chem. Phys. —- 1969, —- V. 51, —- P. 635.
- 21. Mansoori G.A., Carnahan N.F., Starling K.E., Leland T. W. Jr. Equation of State for Nonattracting Rigid Spheres // J. Chem. Phys. —- 1971, —- V. 54, —- P. 1523.
- 22. Hansen J.P., McDonaldI.R. Theory of simple liquids. —- Academic, London, 2006.
- 23. Salacuse J.J., Stell G. Polydisperse systems: Statistical thermodynamics, with applications to several models including hard and permeable spheres // J. Chem. Phys. —- 1982, —- V. 77, —- P. 3714.
- 24. Klay J.L. et al. (E895 Collaboration) Charged pion production in 2A to 8AGeV central Au+Au Collisions // Phys. Rev. C —- 2003, —- V. 68, —- P. 054905.
- 25. Ahle L. et al. (E866 Collaboration) Excitation function of and production in Au+Au reactions at 2â10 AGeV // Phys. Lett. B —- 2000, —- V. 476, —- P. 1.
- 26. Ahle L. et al. (E866 Collaboration) An excitation function of and production in Au+Au reactions at the AGS // Phys. Lett. B —- 2000, —- V. 490, —- P. 53.
- 27. Pinkenburg C. et al. (E895 Collaboration) Production and collective behavior of strange particles in Au+Au collisions at 2â8 A GeV // Nucl. Phys. A —- 2002, —- V. 698, —- P. 495.
- 28. Barrette J. et al. (E877 Collaboration) Proton and pion production in Au+Au collisions at 10.8 AGeV/c // Phys. Rev. C —- 2000, —- V. 62, —- P. 024901.
- 29. Ahle L. et al. (E802 Collaboration) Proton and deuteron production in Au+Au reactions at 11.6 AGeV/c // Phys. Rev. C —- 1999, —- V. 60, —- P. 064901.
- 30. Ahle L. et al. (E802 Collaboration) Kaon production in Au+Au collisions at 11.6 AGeV/c // Phys. Rev. C —- 1998, —- V. 58, —- P. 3523.
- 31. Albergo S. et al. Spectra in 11.6 AGeV/c Au-Au Collisions // Phys. Rev. Lett. —- 2002, —- V. 88, —- P. 062301.
- 32. Back B.B. et al. (PHOBOS Collaboration) Ratios of Charged Antiparticles-to-Particles near Mid-Rapidity in Au+Au Collisions at = 130 GeVs // Phys. Rev. Lett. —- 2001, —- V. 87, —- P. 102301.
- 33. Bearden I.G. et al. (NA44 Collaboration) Particle production in central Pb+Pb collisions at 158 AGeV/c // Phys. Rev. C —- 2002, —- V. 66, —- P. 044907.
- 34. Afanasiev S.V. et al. (NA49 Collaboration) Energy dependence of pion and kaon production in central Pb+Pb collisions // Phys. Rev. C —- 2002, —- V. 66, —- P. 054902.
- 35. Afanasiev S.V. et al. (NA49 Collaboration) Energy and centrality dependence of deuteron and proton production in Pb+Pb collisions at relativistic energies // Phys. Rev. C —- 2004, —- V. 69, —- P. 024902.
- 36. Alt C. et al. (NA49 Collaboration) Energy and centrality dependence of anti-p and p production and the anti-Lambda/anti-p ratio in Pb+Pb collisions between 20/A-GeV and 158/A-Gev // Phys. Rev. C —- 2006, —- V. 73, —- P. 044910.
- 37. Anticic T. et al. (NA49 Collaboration) and Production in Central Pb-Pb Collisions at 40, 80, and 158 AGeV // Phys. Rev. Lett. —- 2004, —- V. 93, —- P. 022302.
- 38. Afanasiev S.V. et al. (NA49 Collaboration) and production in central Pb+Pb collisions at 158 GeV/c per nucleon // Phys. Rev. B —- 2002, —- V. 538, —- P. 275.
- 39. Alt C. et al. (NA49 Collaboration) and Production in Central Pb+Pb collisions at 40 and 158 AGeV // Phys. Rev. Lett. —- 2005, —- V. 94, —- P. 192301.
- 40. Antinori F. et al. (NA57 Collaboration) Energy dependence of hyperon production in nucleusânucleus collisions at SPS // Phys. Lett. B —- 2004, —- V. 595, —- P. 68.
- 41. Antinori F. et al. (NA57 Collaboration) Rapidity distributions around mid-rapidity of strange particles in PbâPb collisions at 158 A GeV/c // J. Phys. G —- 2005, —- V. 31, —- P. 1345.
- 42. Adams J. et al. (STAR Collaboration) Identified Particle Distributions in and Collisions at GeV // Phys. Rev. Lett. —- 2004, —- V. 92, —- P. 112301.
- 43. Bugaev K.A., Gorenstein M.I., Stöcker H., Greiner W. Van der Waals Excluded Volume Model for Lorentz Contracted Rigid Spheres // Phys. Lett. B —- 2000, —- V. 485, —- P. 121.
- 44. Bugaev K.A. Van-der-Waals Gas EOS for the Lorentz Contracted Spheres // Nucl. Phys. A —- 2008, —- V. 807, —- P. 251.
- 45. Wroblewski A.K. On the Strange Quark Suppression Factor in High Energy Collisions // Acta Phys. Pol. B —- 1985, —- V. 16, —- P. 379.
- 46. Rafelski J. Strange anti-baryons from quark-gluon plasma // Phys. Lett. B —- 1991, —- V. 262, —- P. 333.
- 47. Bugaev K.A., Oliinychenko D.R., Ivanytskyi A.I., Cleymans J., Mironchuk E.S., Nikonov E.G., Taranenko A.V., Zinovjev G.M. Separate Chemical Freeze-Outs of Strange and Non-Strange Hadrons and Problem of Residual Chemical Non-Equilibrium of Strangeness in Relativistic Heavy Ion Collisions // Ukr. J. Phys. —- 2016, —- V. 61, —- P. 659.
- 48. Hagedorn R. Statistical thermodynamics of strong interactions at high energies // Nuovo Cim. Suppl. —- 1965, —- V. 3, —- P. 147.
- 49. Moretto L.G., Bugaev K.A., Elliott J.B., Phair L. Hagedorn Thermostat // Europhys. Lett. —- 2006, —- V. 76, —- P. 402.
- 50. Beitel M., Gallmeister K., Greiner C. Thermalization of Hadrons via Hagedorn States // Phys. Rev. C —- 2014, —- V. 90, —- P. 045203, arXiv:1402.1458.
- 51. Beitel M., Gallmeister K., Greiner C. Equilibration of hadrons in HICs via Hagedorn States // J. Phys. Conf. Ser. —- 2016, —- V. 668, —- P. 012057.
- 52. Beitel M., Greiner C., Stoecker H. Fast dynamical evolution of a hadron resonance gas via Hagedorn states // Phys. Rev. C —- 2016, —- V. 94, —- P. 021902.
- 53. Nayak J.K., Banik S., Alam J. The ”horn” in the kaon-to-pion ratio // Phys. Rev. C —- 2010, —- V. 82, —- P. 024914.
- 54. Galitskij V.M., Mishustin I.N. The influence of a phase transition on the shock wave dynamics in nuclear matter // Phys. Lett. B —- 1978, —- V. 72, —- P. 285.
- 55. Stöcker H., Graebner G., Maruhn J.A., Greiner W. Hot, dense hadronic and quark matter in relativistic nuclear collisions // Phys. Lett. B —- 1980, —- V. 95, —- P. 192.
- 56. Kämpfer B. Phase transitions in dense nuclear matter and consequences for neutron stars // J. Phys. G —- 1983, —- V. 9, —- P. 1487.
- 57. Stöcker H., Greiner W. High energy heavy ion collisionsâprobing the equation of state of highly excited hardronic matter // Phys. Rep. —- 1986, —- V. 137, —- P. 277.
- 58. Barz H.W., Csernai L.P., Kämpfer B., Lukacs B. Stability of detonation fronts leading to quark-gluon plasma // Phys. Rev. D —- 1985, —- V. 32, —- P. 115.
- 59. Merdeev A.V., Satarov L. M., Mishustin I. N. Hydrodynamic modeling of the deconfinement phase transition in heavy-ion collisions in the NICAâFAIR energy domain // Phys. Rev. C —- 2011, —- V. 84, —- P. 014907.
- 60. Ivanov Yu. B., Russkikh V. N., Toneev V. D. Relativistic heavy-ion collisions within three-fluid hydrodynamics: Hadronic scenario // Phys. Rev. C —- 2006, —- V. 73, —- P. 044904.
- 61. Landau L.D., Lifshitz E.M. Fluid Mechanics. —- Pergamon, New York, 1979.
- 62. Zel’dovich Y.B., Raiser Y.P. Physics of shock waves and high temperature hydrodynamic phenomena. —- Academic, New York, 1967.
- 63. Rozhdestvensky B.L., Yanenko N.N. Systems of Quasi-Linear Equations. —- Nauka, Moscow, 1978.
- 64. Chodos A. et al. New extended model of hadrons // Phys. Rev. D —- 1974, —- V. 9, —- P. 3471.
- 65. Aoki Y. et al. The QCD transition temperature: Results with physical masses in the continuum limit // Phys. Lett. B —- 2006, —- V. 643, —- P. 46.
- 66. Chin S.A. Transition to hot quark matter in relativistic heavy-ion collision // Phys. Lett. B —- 1978, —- V. 78, —- P. 552.
- 67. Particle Data Group // Phys. Lett. B —- 2004, —- V. 591, —- P. 1.
- 68. Cassing W., Palmese A., Moreau P., Bratkovskaya E. L. Chiral symmetry restoration versus deconfinement in heavy-ion collisions at high baryon density // Phys. Rev. C —- 2016 —- V. 93, —- P. 014902.
- 69. Palmese A. et. al.Chiral symmetry restoration in heavy-ion collisions at intermediate energies // Phys. Rev. C —- 2016 —- V. 94, —- P. 044912.
- 70. Glendenning N.K. Hot metastable state of abnormal matter in relativistic nuclear field theory // Nucl. Phys. A —- 1987 —- V. 469, —- P. 600.
- 71. Garpman S.I.A., Glendenning N.K., Karant Y.J. Thermodynamic behavior of non-strange baryonic matter // Nucl. Phys. A —- 1979 —- V. 322, —- P. 382.
- 72. Theis J. et al. Phase transition of the nucleon-antinucleon plasma in a relativistic mean-field theory // Phys. Rev. D —- 1983 —- V. 28, —- P. 2286.
- 73. Nakai T., Takagi S. Influence of the to Density Isomer and a Phase Transition of the Baryon-Antibaryon Plasma at Finite Temperatures // Prog. Theor. Phys. —- 1984 —- V. 71, —- P. 1118.
- 74. Pisarski R.D., Skokov V.V. How tetraquarks can generate a second chiral phase transition // Phys. Rev. D —- 2016 —- V. 94, —- P. 054008.
- 75. Andronic A. et al. Hadron production in ultra-relativistic nuclear collisions: Quarkyonic matter and a triple point in the phase diagram of QCD // Nucl. Phys. A —- 2010 —- V. 837, —- P. 65.
- 76. Lattimer J.M., Swesty F.D. A Generalized equation of state for hot, dense matter // Nucl. Phys. A —- 1991 —- V. 535, —- P. 331.
- 77. Shen L.H., Toki H., Oyamatsu K., Sumiyoshi K. Relativistic equation of state of nuclear matter for supernova and neutron star // Nucl. Phys. A —- 1998 —- V. 637, —- P. 435.
- 78. Typel S. Variations on the excluded-volume mechanism // Eur. Phys. J. A —- 2016 —- V. 52, —- P. 16.