Scaling properties of spectra in new exact solutions of rotating, multi-component fireball hydrodynamics

We describe fireballs that rehadronize from a perfect fluid of quark matter, characterized by the lattice QCD equation of state, to a chemically frozen, multi-component mixture, that contains various kinds of observable hadrons. For simplicity and clarity, we apply a non-relativistic approximation to describe the kinematics of this expansion. Unexpectedly, we identify a secondary explosion that may characterize fireball hydrodynamics at the QCD critical point. After rehadronization, the multi-component mixture of hadrons keeps on rotating and expanding together, similarly to a single component fluid. After kinetic freeze-out, the effective temperature $T_{i}$ of the single-particle spectra of hadron type $h_i$ is found to be a sum of the kinetic freeze-out temperature $T_f$ (that is independent of the hadron type $h_i$) and a term proportional to the mass $m_i$ of hadron type $h_i$. The coefficient of proportionality to $m_i$ is also found to be independent of the hadron type $h_i$ but be dependent on the radial flow and vorticity of collective dynamics.


Introduction
In this paper, the main theme of the analysis is the mass systematics of the single particle spectra as measured in relativistic heavy ion collisions, based on exact solutions of fireball hydrodynamics. Interestingly, the NA44 Collaboration observed signals of collective expansion in symmetric heavy-ion collisions in 200 AGeV S+S and 158 AGeV Pb+Pb reactions at CERN SPS: the transverse mass spectra of pions, kaons and protons was found to be nearly exponential in these reactions, with a slope parameter that increased nearly linearly with the mass of the observed hadrons [1]. This linear mass dependence of the slope parameters was interpreted as an evidence of collective transverse flow in these heavy ion induced central collisions.
The inverse slope parameters for positively and negatively charged pions and kaons, as well as for protons and anti-protons has subsequently also been measured in various heavy ion collisions, for example in √ s NN = 200 GeV Au+Au collisions by the PHENIX Collaboration in three centrality bins, as is shown very clearly on Figure 10 and in Table IV of ref. [2], reproduced and detailed on Fig. 1. These inverse slope parameters were also found to increase with increasing particle mass, in all centrality bins. This increase for central collisions was found to be more rapid for heavier particles, taken as another indication of collective expansion dynamics not only in central but also in peripheral heavy ion reactions at RHIC. Such a behavior was predicted -under certain conditions corresponding to a non-relativistic radial expansion dynamics -in refs. [3][4][5][6], for central collisions. A similar behaviour was derived for non-central heavy ion collisions as well, in refs. [7,8].
The early derivations were based on the blast-wave [3], or on the Buda-Lund model [4,5] that focused on a parameterization of the particle emission in phase-space around the time of kinetic freeze-out. It is remarkable that the transverse, radial flows in both the blast-wave and the Buda-Lund parameterization go back to precisely the same exact solution of non-relativistic fireball hydrodynamics, namely the Zimányi-Garpman-Bondorf solution of ref. [9]. A few years after the successfull hydrodynamical parameterizations of the freeze-out phase-space distribution were obtained, the time evolution was also accessed with the help of the first exact solution within the Buda-Lund family of exact solutions of fireball hydrodynamics. The first of such Buda-Lund hydro solutions was found for a spherically symmetric, radially expanding fireball with spatially homogeneous temperature profile and a Gaussian density profile in ref. [6], where the inverse slope parameter of the single particle spectra, T eff has also been shown to scale as an affine-linear function of the particle mass m: where m is the mass of a single kind of elementary particle with three kinetic degree of freedom, that consists the expanding medium. The average transverse flow is denoted by u T and the kinetic freeze-out temperature is denoted by T f . Although the above relation is typically quoted as a justification for the linear mass dependence of the effective temperature in a hot and dense, hadronic matter, actually the derivations had a much more limited scope, evaluating the mass dependence for one kind of hadrons only, while in the experiment, a mixture of various hadronic components is observed.
Such an experimental result is beautifully illustrated on Fig. 10. of ref. [2] and reproduced also on our Fig. 1, that compares measured slope parameters in relativistic heavy ion collisions at mid-rapidity with eq. (1), derived in a non-relativistic context. The success of this comparision suggests that the key scaling properties of the transverse mass spectra at mid-rapidity may perhaps be understood in the framework of non-relativistic kinematics. Inspired by this insight, we try to handle the equations of relativistic and non-relativistic hydrodynamics in a similar way in Section 2. For the sake of clarity and brevity, we present new, rotating ellipsoidal exact solutions of fireball hydrodynamics only in the non-relativistic limit, as detailed in Section 3.
Note, however, that all the earlier derivations dealt with a hydrodynamically expanding medium that had only a single component, i.e. only one kind of hadrons, with a given mass m. The linear rise in the data was compared with the results of the calculations by extrapolating the theoretical results as a continuous function of the mass m, that was considered as a smoothly varying parameter of the solution and the resulting effective temperature or slope parameter. However, the experimentally observed mass spectrum of the hadrons is essentially discrete, and not a smoothly varying function. In the data analysis, for example, the slope parameters were measured at the mass of pions, kaons and (anti)protons, as indicated on Fig. 1. One of the challenges considered in this work is to generalize the derivation from a single-component hadron gas to a multi-component hadron gas, that is a mixture of hadrons with different discrete values of their masses. Such a derivation seemed to be an almost hopeless theoretical challenge, as far as we know it was not even attempted before. This challenge however is positively solved in this manuscript. We derive a new class of solution for the equations of fireball hydrodynamics for a collectively expanding and rotating, strongly interacting perfect fluid that contains a multi-component hadronic matter. We also obtain the slope parameters for various kinds of hadrons that emerge from the collective expansion after kinetic freeze-out.
For the sake of clarity, we introduce the i index to distinguish the different type of hadrons thus the contribution of each kind of particle to the inverse slope can be separated. Note that in this way one can find what is really eye-catching in Fig. 1: the inverse slope T i and as well as the single particle spectra depends on the particle type only through the mass of the hadrons: the freeze-out temperature T f and the slope parameter u T is independent of the type i of the observed hadrons, furthermore the freeze-out temperature is a static characteristics of the medium, apparently independent of the centrality, while the coefficient of linearity is increasing with increasing centralities. The result is summarized as  Figure 10 and Table IV of ref. [2]. Dashed lines are drawn to guide the eye, corresponding to eq. (1) and also to earlier theoretical results for a single-component hadron gas in non-central heavy ion collisions as detailed in refs. [7,8]. Green symbols indicate our new result, based on our new exact solutions of fireball hydrodynamics for a collectively expanding and rotating, multi-component hadronic matter, as indicated by symbols with "MC theory" in the legends. These green symbols appear at discrete values of hadronic masses, corresponding to eq. (2), along the eye-guiding dashed lines.
In this eq. (2), the mass m i is not a smoothly varying parameter, in contrast to eq. (1), but takes on the discrete values, corresponding to the observable, PDG listed masses of those hadron species h i that constitute the mixture in the expanding, multi-component hadron gas.
As an extra difficulty, we also consider and are interested in describing the time evolution of fireball hydrodynamics using a realistic, lattice QCD based equation of state of ref. [10]. This excercise is a kind of academic part of our study, motivated by the fact that in lattice QCD simulations, the speed of sound and the ratio of pressure to energy density at vanishing baryochemical potential are known to depend significantly on the temperature [10]. We consider such an equation of state in this work, indicated also on Fig. 2, although in a future, even more realistic scenario relativistic kinematics will have to be considered as well. In this sense, our present study is also preparatory, yet difficult step, aiming at future studies using lattice QCD equation of state in the dynamics, together with relativistic kinematics and a multi-component hadronic matter in the final state.
By this token, the main research topic of this manuscript is given: can the lattice QCD Equation of State and the scaling behaviour, the mass systematics of the slope parameter of the single particle spectra be understood in a self-consistent, hydrodynamical picture?

Equations of relativistic as well as non-relativistic hydrodynamics
The basic equations of relativistic, perfect fluid hydrodynamics consists of a set of partial differential equations, that express the local consvervation of entropy, momentum and energy: where σ stands for the entropy density, T µν is the energy-momentum tensor of a perfect fluid in the co-moving frame: Here the metric tensor is denoted by g µν = diag(1, −1, −1, −1), and u µ = γ(1, v) is the four-velocity field, normalized as u µ u µ = 1. The coordinates are denoted as x µ = (t, r) = t, r x , r y , r z and the four-momentum is defined as Eq. (3) expresses local entropy conservation, while eq. (4) expresses local conservation of energy and momentum. This equation can be projected to a component parallel to u µ that yields the energy equation: while the component pseudo-orthogonal to the four-velocity field yields the relativistic Euler equation: Local conservation laws of densities cannot be used generally, because in the initial, hot temperature stage, at vanishing baryochemical potential, the only conserved charge is the net baryon numberwhich is zero in the considered case. As the medium cools down quark and gluon degrees of freedom get confined to hadrons. After hadro-chemical freeze-out, this system of equations is supplemented with the continuity equation of particle densities. One new conservation law is obtained for each of the hadrochemically frozen particle types: where n i is the particle density of the i th hadron, and j counts that how many kind of hadrons are frozen out hadrochemically.
In the high temperature phase, T > T chem , we have five equations, eqs. (3,6,7), for six independent variables, the three independent spatial components of the four-velocity field v and the energy, the pressure and the entropy density. The energy equation can be transformed to a temperature equation, given in Table 2, that introduces a new variable, the temperature, utilizing also the ε = κ(T)p equation of state, resulting in five equations for six independent quantities, but replacing the unknown energy density with the unknown temperature field and a known κ(T) function. This set of equations becomes closed with the fundamental equation of thermodynamics, and the assumption that the enthalpy is dominated by entropy density times the temperature at such a high temperatures, summarized in the left column of Table 1.
The state before the hadronization corresponds to the quark matter (QM) and after the transition the hadronic matter (HM) is formed. Different thermodynamical approaches may be used to describe the different phases. The general form of the fundamental equation of thermodynamics is: In heavy ion collisions at RHIC and LHC energies, the bariochemical potential (µ) of the strongly coupled quark-gluon plasma is close to zero. Consequently in this case we can use the approximation. As the medium cools down and a chemically frozen hadron gas is formed, the hadrochemical potential approaches the particle mass m i for hadron type i and the Tσ term becomes negligible compared to the energy of the rest-mass. In this way, in the chemically frozen hadronic matter one may approximate the enthalpy density as At lower temperatures, T ≤ T chem , we have additional relations, namely the continuity equations for each of the chemically frozen hadronic species. Suppose that there are j such chemically frozen hadrons (i = 1, ..., j), then we have j − 1 independent new equations, in addition to the entropy conservation. This way we have j + 4 equations for j + 5 unknowns, that we take as the j different densities, the temperature, the pressure, and the three independent spatial components of the four-velocity field. Similarly to the high temperature case, the energy equation can be transformed again to a temperature equation, given in the right column of Table 2. Indeed, this transformation utilizes also the ε = κ(T)p equation of state, similarly to the high temperature case, and results again in j + 4 equations for j + 5 unknowns. The κ(T) function is assumed to be known from lattice QCD. This set of equations is closed, similarly to the high temperature phase, with the fundamental equation of thermodynamics, and the assumption that the chemical potentials are given approximately by the masses of hadrons and that the sum over chemical potentials times the corresponding densities is the dominant term in this low temperature phase, as summarized in the right column of Table 1. Actually, such an approximation is justified in the non-relativistic kinematic region, as it leads naturally to the usual form of the non-relativistic Euler equation.
The medium behaves differently before and after the hadronization which is manifested in the different system of equations in the corresponding hydrodynamical description. The relativistic system of equations of the two phases are summarized in Table 2.
From the energy equation, we obtain a differential equation for the temperature (T) by utilizing the equation of state and the expressions for the pressure. Note that this procedure can be followed not only in the case of relativistic kinematics, as summarized in Table 2, but also exactly the same method can be Table 1. We assume, that at high temperature, the enthalpy density is dominated by the entropy term while below the chemical freeze-out temperature of hadrons, their rest mass dominates the same term. In these equations, j counts how many kind of different hadrons are hadrochemically frozen-out.
used to obtain the temperature equation in case of non-relativistic kinematics, and this way one obtains a striking similarity between the system of partial differential equations of fireball hydrodynamics both in the relativistic, and in the non-relativistic kinematic domain. Consequently in the |v| 2 1 approximation, the system of partial differential equations of relativistic hydrodynamics of Table 2 directly correspond to the system of partial differential equations of non-relativistic hydrodynamics, as summarized in Table 3.

Exact and analytic solutions of fireball hydrodynamics
In this section, we focus on the non-relativistic approximation, for the sake of clarity and simplicity, leaving the discussion of the relativistic kinematics to a follow-up publication. Although such a non-relativistic approximation limits the direct applicability of our results to a detailed comparision with particle production to low transverse momentum and to mid-rapidity, similar simple and exact non-relativistic fireball solutions provided already important insights to the mass systematics of the single-particle spectra at low p T and at mid-rapidity.
Linear mass dependencies of the slope parameters of the single particle spectra were obtained before in exact solutions of non-relativistic fireball hydrodynamics, where a spatially homogeneous temperature profile was matched with a Gaussian density profile, both for spherically [6] and for ellipsoidally symmetrically expanding fireballs [7,8], corresponding to central and non-central heavy ion collisions, respectively. However, these Gaussian solutions are readily generalized to an arbitrary but still spherically [11] or ellipsoidally symmetric temperature profiles [12], where maxima of the Table 3. The system of partial differential equations that describe non-relativistic fireball hydrodynamics of quark matter (QM) at high, T > T chem temperatures and a hadrochemically frozen, multi-component hadronic matter (HM) at lower, T ≤ T chem temperatures. Similarly to the relativistic case, we have obtained 5 equations for 5 unknowns in the high temperature, QM phase, while we have j + 4 equations for j + 4 unknowns in the chemically frozen, multi-component HM phase.
Tκ 1+κ n i T temperature profile function correspond to local minima in the matching density profile functions. Similarly, axially or ellipsoidally symmetric expanding fireballs are described by certain recently found exact solutions of perfect fluid hydrodynamics in the relativistic kinematic domain as well, that are characterized by a scaling variable and an arbitrary positive definite scaling function for the initial temperature profile, that is matched with a corresponding density profile [14,15]. We have summarized in Tables 1,2,3 the systems of partial differential equations of a fireball that is created in (non-)relativistic heavy-ion collisions. It is important to note that in the temperature equation the coefficient of the co-moving derivative of the temperature is the same in both Tables 2  and 3, and this similarity between the relativistic and non-relativistic kinematics may play a role in subsequent, future studies.
Given such a background, let us present in this section two new, rotating exact solutions of non-relativistic fireball hydrodynamics. Their validity can be straigthforwardly tested with the help of Table 3, so we do not detail their derivation here. These solutions listed in the following subsections are given in the laboratory frame K, assuming a collider type of experiment so that the center of the fireball is at rest in the laboratory frame K.
The exact and analytic solutions that are already known, can be grouped in a new manner as well, into three different kind of solutions, noted first in ref. [18].
1. The first kind of scaling solutions of fireball hydrodynamics are characterized by spatially homogeneous temperature profile and corresponding Gaussian (entropy) density profile. These solutions may also feature a realistic, temperature dependent energy density/pressure ratio, or a corresponding temperature dependent speed of sound, see for example refs. [6][7][8]13,18,20]. 2. The second kind of scaling solutions of fireball hydrodynamics are characterized by arbitrary, spatially inhomogeneous temperature profiles and corresponding, matching density profiles. The price for having the freedom of a realistic, arbitrary initial temperature profile is the need to have a temperature independent speed of sound, see refs. [11,12,[14][15][16][17][18]. 3. The existence of third kind of scaling solutions of fireball hydrodynamics was noted first in ref. [18], but elaborated only in ref. [21]. These third kind of solutions may gain further importance in the future, as they allow to use a lattice QCD equation of state in a realistic, spatially inhomogeneous, scaling solution of fireball hydrodynamics.

New solutions of the first kind:
Exact, parametric solutions of rotating and expanding fireballs with Gaussian density profiles The first kind of scaling solutions of fireball hydrodynamics correspond to homogeneous temperature profiles, with a corresponding Gaussian (entropy)density profile, and a temperature dependent κ = κ(T) = ε/p function.
In this subsection we generalize the results of ref. [20], to a quark matter rehadronizing to a multi-component hadron gas. Following ref. [20], we have written up the velocity field as a sum of two terms: a "Hubble-term" v H , and a "rotational term" v R . This form of the velocity field v and the scaling variable s satisfies the scale equation: The Hubble term has zero curl (and thus does not contribute to the vorticity of the flow), while the rotational term has zero divergence, hence it does not contribute to the overall expansion of the fluid: ∇v R = 0. Table 4. Equation of state independent part of new, exact, rotating solutions of fireball hydrodynamics, as given in the laboratory frame K, where a rotating ellipsoid with time dependent principal axis X, Y and Z is not only expanding but also spinning around the r y direction. The time dependent parameter ϑ stands for the angle between the principal axis X and the impact parameter direction, r x in the intertia frame K. The first of this kind of solutions was given in ref. [20], that we generalize in this work to quark matter and multi-component hadronic matter scenarios.

Relations valid both in Quark Matter (QM) and in Hadronic Matter (HM):
The terms in v that are proportional to the angular velocityθ contribute to the rotational flow velocity v R , which determines the vorticity of the solution as The vorticity vector is parallel with the axis of rotation, and the value of its only non-vanishing component in the laboratory frame, ω = (ω x , ω y , ω z ) , where the non-vanishing component is given as The equation of state specific parts of the solution are summarized in Table 5. In the same table the acceleration parameter D is related to m , the average mass of the particles in the multi-component hadronic matter: In such a multi-component hadronic matter, the average mass is assumed to be time independent, m = (∑ j i=1 n i,0 m i )/(∑ j i=1 n i,0 ).

Valid in QM only:
Valid in HM only: In this way we reduced the complicated set of partial differential equations to a set of ordinary differential equations. This result generalizes the triaxial, rotating and expanding fireball solution of ref. [20] from a single component hadron gas, characterized by mass m, to a perfect fluid of quark matter that hadronizes to a mixture of hadrons with average mass m . Table 5 indicates that this mixture expands and rotates together, because the scale parameters X, Y and Z as well as the parameters of the rotation ω are independent of the index i of hadron types h i . These solutions are all (rotating) Gaussian solutions, due to the exp(−s/2) type of terms in the entropy and chemically frozen density profiles given in Table 5, that correspond to a position independent, but time dependent temperature profiles.

Fireball explosion at the QCD Critical Point
The last line of Table 5 indicates a beautiful exact and analytic result: Although the lattice QCD EoS suggests that the function κ(T) is a smoothly varying function of the temperature as indicated on Fig. 2, the coefficient ofṪ/T in the last line of Table 5 actually changes. The difference of these coefficients is which vanishes at the temperature of T max ≈ 158 MeV, where the κ(T) function has a maximum and changes sign at this point: • If T chem > T max , then the coefficient ofṪ/T increases, and the rate of change of the temperature decreases at the same logarithmic derivative of the volumeV/V, characterizing the volume expansion at hadrochemical freeze-out. • If T chem = T max , the time evolution of the temperature changes smoothly at the hadrochemical freeze-out. • If T chem < T max , the coefficient ofṪ/T decreases, and the rate of change of the temperature increases at the same logarithmic derivative of the volumeV/V.
As we shall see below, not only the rate of change of the temperature with the increase of the volume is sensitive to quark confinement and subsequent hadrochemical freeze-out, but also the dimensionless measure of the rate of acceleration of the expansion changes at the chemical freeze-out, due to the changes in the expansion dynamics.
We assume that the hadronization and the hadrochemical freeze-out happens nearly simultaneously, as a consequence T chem ≈ T c in this case. For a spatially homogeneous temperature profile, and a simultaneous quark confinement and hadrochemical freeze-out, the boundary conditions are: Due to these special boundary conditions one can obtain an important relation for D that is a dimensionless measure of the rate of expansions of the scales in the equations of motion: The expression of the left side in eq. (24) characterizes the dimensionless acceleration parameter D Q in the quark phase, and the right side of the same relation relates D H to the acceleration of the hadronic matter. Given that D H > D Q , the acceleration of the scales jumps at t chem in every direction, so the directional Hubble constantsẊ/X,Ẏ/Y andŻ/Z have a break at these temperatures, and each of the scales X, Y and Z starts to expand faster as the temperature cools below the hadrochemical freeze-out temperature T chem . We propose to call this phenomenon as a second explosion in the Little Bangs or heavy ion collisions. The first explosion is the violent expansion that starts just after collision, due to large initial energy densities, temperatures and pressure gradients, while this second explosion starts just after the conversion to the hadronic phase and appears due to the role of hadrochemical freeze-out. The effect of the second explosion becomes smoothened, if the average mass in the fireball is increased at T ≈ T chem . We expect, but cannot detail in this paper that a smooth cross-over may further slow down the expansion dynamics, but the detailed discussion of such a phenomena goes beyond the scope of this manuscript. Our calculations seem to be laying the ground for similar calculations in the relativistic kinematic domain. Before looking for such new solutions of relativistic hydrodynamics, let us perform a consistency check and see if the single particle spectra from such solutions looks to be realistic or not.

Evaluation of the single-particle spectra
The single particle spectra has already been calculated from single-component (SC) hydrodynamical solutions [19,20]. In this section we follow the notations and conventions of these earlier works, without giving the full details, referring also to the Introduction for a motivation, and highlighting only the main, characteristic features of our results.
The slope parameters of the single particle spectra are apparently simple also in our multi-component (MC) generalization of rotating and expanding fireballs, when given in the K frame that co-rotates with the fireball: where k, l = {1, 2, 3} = {x, y, z} and (T ) −1 kl,i is the inverse of the effective temperature matrix or in other words the inverse slope matrix.
Similarly to the hydrodynamical solutions of a single-component hadron gas [19,20], this single particle spectra can be transferred to the laboratory frame and we obtain: We have collected the matrix elements for multi-component (MC) and single component (SC) hadronic matter (HM) in Table 6, where we denote these final state hadronic observables with the subscript f that denotes kinetic freeze-out. Note that for the non-rotating fireballs, the inverse temperature matrix is diagonal in the center of mass frame K, while for rotating fireballs, the inverse temperature matrix is diagonal only in the co-rotating frame K , so an additional transformation to K has to be made, similarly how it was performed in ref. [19] . Table 6. The comparison of the inverse slope parameters after kinetic freeze-out at temperature T f , for a Single Component (SC) and for a Multi-Component (MC) Hadronic Matter (HM) from exact and analytic solutions of fireball hydrodynamics, in the case of rotating and nonrotating ellipsoid as well. In case of non-rotating fireballs, the temperature matrix is diagonal in the center of mass frame K, T kl = diag(T x , T y , T z ). For rotating fireballs, the temperature matrix is diagonal only in the co-rotating frame K , where T kl = diag(T xx , T yy , T zz ). [8,19,20] MC HM

SC HM
We thus find that starting from a high temperature quark matter phase, followed by a nearly simultaneous quark confinement and hadrochemical freeze-out, we obtain apparently thermal single-particle spectra, where the inverse slope parameters have linear mass for each of the hadronic components h i and they depend on the particle types only through their masses m i , similarly to the experimental data summarized in the Introduction. After kinetic freeze-out, the effective temperature T i of the single-particle spectra of hadron type h i is found to be a sum of the kinetic freeze-out temperature T f (that is independent of the hadron type h i ) and a term proportional to the mass m i of hadron type h i . The coefficient of proportionality to m i is also found to be independent of the hadron type h i but be dependent on the radial flow (and vorticity) of collective dynamics.
As a conclusion, the presented new class of multi-component, exact solution of fireball hydrodynamics provides a clear-cut, but not-yet-relativistic explanation for the mass-scaling behaviour of the single particle spectra in heavy ion collisions. Scaling violations, corresponding to deviations from the linear mass scaling may be expected due to relativistic kinematics and temperature gradients, as discussed in ref. [5], or, due to possible non-equilibrium kinetic freeze-out effects [23].

Summary
We have presented new, exact, parametric solutions of fireball hydrodynamics for a quark matter hadronizing to a chemically frozen, multi-component hadronic matter. These solutions have spatially homogeneous initial temperature profiles and corresponding Gaussian (entropy) density profiles. We have identified a change in the time evolution of the temperature that is due to the quark confinement to hadrons and to the simultaneous hadrochemical freeze-out. In addition, we have found a new phenomenon at the hadrochemical freeze-out, that relates to the dynamics of the hadronizing fireball: a second explosion, that seems to be, rather counter-intuitively, the strongest if the energy density to pressure ratio is the largest (frequently referred to as the softest point of QCD). This second explosion is caused by the increase of the acceleration after hadrochemical freeze-out, which effect seems to appear due to the change of the thermodynamics from a quark matter where no conserved charges exist to a chemically frozen hadronic matter, where each of the chemically frozen hadronic types obeys a local continuity equation. We conjecture, that this second explosion is washed out if the transition from quark matter to hadronic matter corresponds to a continuous crossover, however, we could not detail this conjecture and elaborate the description of such a cross-over transition in this manuscript.
We have successfully generalized earlier exact and analytic hydrodynamical solutions from a single-component to a multi-component hadronic matter and we have shown that after the hadrochemical freeze-out of the multi-component hadronic matter, the slope parameters of the single particle spectra for various hadrons follow the experimentally observed affine-linear mass-scaling behaviour. Scaling violations, or deviations from such an affine-linear mass-scaling of the slope parameters may also be expected, due to relativistic kinematics and/or temperature gradients [5], or, further possible non-equilibrium effects [23].