One of the most important directions of physics of heavy ion collisions is related to a location of the (tri)critical endpoint ((3)CEP) of the quantum chromodynamics (QCD) phase diagram. At present the lattice QCD cannot tell us whether at high baryonic charge densities the chiral symmetry restoration (CSR) phase transition (PT) and the deconfinement one of color degrees of freedom are two different PTs or a single one. On the other hand, a few advances approaches support an idea that at finite values of baryonic chemical potential a phase with a partial CSR occurs before the deconfinement PT [1,2,3]. However, one of the hardest problem of QCD phenomenology is to determine the collision energy threshold of these PTs using the existing experimental data obtained in the central nucleus-nucleus collisions.
|No and Type||Signal||C.-m. energy (GeV) Status||C.-m. energy (GeV) Status|
|1. Hydrodynamic||Highly correlated quasi-plateaus in entropy/baryon, thermal pion number/baryon and total pionnumber/baryon. Suggested in [11,12].||Seen at 3.8-4.9 GeV [4,5]. Explained by the shock adiabat model [4,5]. Require an explanation.||Seen at 7.6-9.2 GeV [4,5].|
|2. Thermodynamic||Minimum of the chemical freeze-out volume .||In the one component HRGM it is seen at 4.3-4.9 GeV . In the multicomponent HRGM it is seen at 4.9 GeV . Explained by the shock adiabat model [4,5].||Not seen.|
|3. Hydrodynamic||Minimum of the generalized specific volume at chemical freeze-out.||Seen at 4.9 GeV . Explained by the shock adiabat model [4,5].||Seen at 9.2 GeV . Require an explanation|
|4. Thermodynamic||Peak of the trace anomaly .||Strong peak is seen at 4.9 GeV . Is generated by the peak on the shock adiabat at high density end of the mixed phase .||Small peak is seen at 9.2 GeV . Require an explanation|
|5. Thermodynamic||Peak of the baryonic density .||Strong peak is seen at 4.9 GeV . Is explained by .||Strong peak is seen at 9.2 GeV . Require an explanation|
|6. Thermodynamic||Apparent chemical equilibrium of strange charge.||is seen at 4.9 GeV . Explained by thermostatic properties of mixed phase at .||is seen at 8.8 GeV [10,13]. Explained by thermostatic properties of QG bags with Hagedorn mass spectrum .|
|7. Fluctuational (statistical mechanics)||Enhancement of fluctuations||N/A||Seen at 8.8 GeV . Can be explained by CEP  or 3CEP formation .|
|8. Microscopic||Strangeness Horn ( ratio)||N/A||Seen at 7.6 GeV. Can be explained by the onset of deconfinement at /above  8.7 GeV.|
Fortunately, during last few years an essential progress in resolving such a problem was achieved [4,5][7,8]. In particular, two sets of remarkable hydrodynamic and thermodynamic signals of two PTs at the center of mass collision energies GeV and GeV were found in  and the hypothesis of their possible observation at these energies of collision was first formulated in [4,5]. In the works [7,8] a very good description of the large massive of experimental data on nuclear collisions was first achieved with the Parton-Hadron-Sring-Dynamics (PHSD) model by assuming an existence of CSR PT at about GeV in a hadronic phase and a deconfinement one at GeV. In 2017 the group of scientists analyzed the fluctuations of light nuclei and came to a conclusion that the vicinity of collision energy GeV is a nearest vicinity of the critical endpoint of the QCD phase diagram . The arguments of such a statistical signal of the endpoint were essentially enhanced in  and the conclusion of the tricritical endpoint existence in QCD at or slightly above GeV was first formulated in . The summary of found signals is given in Table 1.
Although the PHSD model provides us with some hints about the properties of the phase existing at the collision energy range GeV, the question is whether one can independently get the properties of this matter. In this work we briefly show how one can get them from the equation of state (EoS) which is obtained in [4,5] from fitting the data.
2. Hadron Resonance Gas Model with Hard-Core Repulsion
The possible PTs signals 1-6 presented in Table 1 were obtained with the help of the multicomponent Hadron Resonance Gas Model (HRGM) [4,5][18,19], which, in contrast to the HRGM with one or two hard core radii of hadrons  has the following hard-core radii of pions =0.15 fm, kaons =0.395 fm, -hyperons =0.085 fm, other baryons =0.365 fm and other mesons =0.42 fm. Thus, having only 2 or 3 additional global fitting parameters compared to the usual HRGM , one can get extremely good description of the hadronic multiplicity ratios measured at AGS, SPS and RHIC energies with a high quality [6,19], including traditionally the most problematic ones for the usual HRGM , i.e. , and ratios.
A high quality fit of hadronic multiplicity ratios achieved by the multicomponent HRGM gives us a high confidence that the EoS of hadronic matter is now fixed with high accuracy in the wide range of chemical freeze-out (CFO) temperature and baryonic chemical potential . This conclusion was thoroughly verified recently with the newest version [18,19] of the multicomponent HRGM which allows one to go beyond the Van der Waals approximation traditionally used in HRGM. In the simplest case of a single hard-core radius of hadrons the HRGM pressure in the grand canonical ensemble is
(1) where the sum is running over all particles (and antiparticles) with the chemical potentials , is the excluded volume of hadrons and is their proper volume. Here denotes the partial pressure of the point-like hadrons of sort with the degeneracy and the mass , while is the energy of particle with the 3-momentum and is the effective chemical potential. The parameter defines the Fermi ( ), the Bose ( ) or the Boltzmann ( ) statistics. Then the thermal density of particles of sort is defined as
(2) where denotes the particle number density of point-like hadrons of sort .
3. Necessity of Multicomponent HRGM
There are two main reasons of why the HRGM of Eqs. (1) and (2) with the Van der Waals is used to determine the CFO parameters. The main reason is that for the hard-core repulsion the energy per particle is same as in the ideal gas and, therefore, there is no need to transform the potential energy of the system into the kinetic energy of particles. Of course, one could add the attractive term to the pressure (1), but in this case one would face a hard mathematical problem to convert the interacting gas into free streaming hadrons [20,21,22] measured in the experiments.
The second reason is that only the hard-core repulsion provides the consistence with the lattice QCD results. In other words, if one takes into account all hadrons as the point-like particles with , then it is well-known that at high and their pressure will dramatically exceed the pressure of quarks and gluons. However, in order to provide a high quality fit of the data such a simplified model should be modified in two respects. First of all one should remember that the quantum second virial coefficients of particle of sort interacting with the particle of sort is 
(3) where the term is the virial coefficient due to quantum statistics of hadron of sort which is expressed in terms of the densities of auxiliary Boltzmann hadrons of the same sort , and the term is due to attractive interaction. This equation shows that the gas pressure (1) with the hard-core repulsion, indeed, accounts for the quantum properties of hadrons, if . It also shows that, if one introduces the different hard-core radii of hadrons, then one can even account for the attraction between them at the level of the second virial coefficient which is sufficient for the low particle densities at CFO. Of course, for all hadrons the second virial coefficients (3) are temperature dependent, but fortunately, at high temperatures such a dependence is not strong  and, hence, to a leading order one can restrict the treatment by the constant hard-core radii.
In addition, the HRGM pressure corresponds to the hadron resonances of vanishing width. This is, of course, a rough approximation because at CFO the density is sufficiently low that the inelastic reactions between hadrons can be neglected and, hence, the hadrons and their resonances should get their vacuum masses and vacuum widths before going into detector. The other reason to introduce the widths is the practical one. Thus, using the Briet-Wigner parameterization of resonance width of all hadronic resonances one can describe the hadron rations essentially better than with the Gaussian one or with the vanishing width . Therefore, it seems that the most efficient way to account for the residual attractive interaction between hadrons at CFO and to achieve a high quality of the hadron multiplicity description is to generalize the one component HRGM (1), (2) to the multicomponent case, i.e. to account for different hard-core radii and then to determine these radii from the fit of experimental data. This is exactly what was done in Refs. [4,5][18,19] during last five years. It is evident that the hard-core radii determined in this way are effective ones by construction.
4. EoS of Hadronic Matter with CSR
Using the multicomponent HRGM in Refs. [4,5] it was possible from fitting the entropy per baryon along the shock adiabat [11,12] to determine the EoS of the phase existing at the collision energy range GeV. This EoS is similar to the MIT-Bag model
(4) but the coefficients , , , and are rather different from what is predicted by the perturbative QCD for massless gluons and (anti)quarks. In Ref.  the EoS (4) was suggested to find out the number of bosonic and fermionic degrees of freedom of this phase. Recalling that first three terms of the EoS (4) correspond to the gas of massless particles and noting that the coefficient is small and its value is comparable to its own error, we could determine the numbers of total , bosonic and fermionic degrees of freedom as
(5) Since the numbers and are much larger than the corresponding number of degrees of freedom in perturbative QCD, but at the same time they are close to the total number of spin-isospin degeneracies of all known hadrons, in Ref.  we, independently of the works [7,8], concluded that the EoS (4) corresponds to the gas of massless hadrons with strong attraction given by the vacuum pressure .
Here we present a summary of possible signals of CSR and deconfinement PTs which may be, respectively, probed at the collision energies at GeV and above GeV. Also these signals may evidence for an existence of the tricritical endpoint of QCD phase diagram at the collision energy around or slightly above GeV. The EoS of the hadronic matter with CSR is discussed and the number of bosonic and fermionic degrees of freedom is found.
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ç ao para a Ciência e Tecnologia".