Extreme states of nuclear matter - 1980

The theory of hot nuclear fireballs consisting of all possible finite-size hadronic constituents in chemical and thermal equilibrium is presented. As a complement of this hadronic gas phase characterized by maximal temperature and energy density, the quark bag description of the hadronic fireball is considered. Preliminary calculations of temperatures and mean transverse momenta of particles emitted in high multiplicity relativistic nuclear collisions together with some considerations on the observability of quark matter are offered.


Overview
I wish to describe, as derived from known traits of strong interactions, the likely thermodynamic properties of hadronic matter in two different phases: the hadronic gas consisting of strongly interacting but individual baryons and mesons, and the dissolved phase of a relatively weakly interacting quark-gluon plasma. The equations of state of the hadronic gas can be used to derive the particle temperatures and mean transverse momenta in relativistic heavy ion collisons, while those of the quarkgluon plasma are more difficult to observe experimentally. They may lead to recognizable effects for strange particle yields. Clearly, the ultimate aim is to understand the behavior of hadronic matter in the region of the phase transition from gas to plasma and to find characteristic features which will allow its experimental observation. More work is still needed to reach this goal. This report is an account of my long and fruitful collaboration with R. Hagedorn [1].
The theoretical techniques required for the description of the two phases are quite different: in the case of a The original address byline 1980 : Gesellschaft für Schwerionenforschung mbH, Darmstadt and Institut für Theoretische Physics der Universität Frankfurt/M; originally printed in GSI81-6 Orange Report, pp. 282-324, edited by R. Bock and R. Stock. b e-mail: rafelski@physics.arizona.edu hadronic gas, a strongly attractive interaction has to be accounted for, which leads to the formation of the numerous hadronic resonances -which are in fact bound states of several (anti) quarks. If this is really the case, then our intuition demands that at sufficiently high particle (baryon) density the individuality of such a bound state will be lost. In relativistic physics in particular, meson production at high temperatures might already lead to such a transition at moderate baryon density. As is currently believed, the quark-quark interaction is of moderate strength, allowing a perturbative treatment of the quark-gluon plasma as relativistic Fermi and Bose gases. As this is a very well studied technique to be found in several reviews [2][3][4][5][6][7][8], we shall present the relevant results for the relativistic Fermi gas and restrict the discussion to the interesting phenomenological consequences. Thus the theoretical part of this report will be devoted mainly to the strongly interacting phase of hadronic gas. We will also describe some experimental consequences for relativistic nuclear collisions such as particle temperatures, i.e., mean transverse momenta and entropy.
As we will deal with relativistic particles throughout this work, a suitable generalization of standard thermodynamics is necessary, and we follow the way described by Touschek [9]. Not only is it the most elegant, but it is also by simple physical arguments the only physical generalization of the concepts of thermodynamics to relativistic par-ticle kinematics. Our notation is such that = c = k = 1. The inverse temperature β and volume V are generalized to become four-vectors: where u μ , v μ , and w μ are the four-velocities of the total mass, the thermometer, and the volume, respectively. Usually, u μ = v μ = w μ . We will often work in the frame in which all velocities have a timelike component only. In that case we shall often drop the Lorentz index μ, as we shall do for the arguments V = V μ , β = β μ of different functions.
The attentive reader may already be wondering how the approach outlined here can be reconciled with the concept of quark confinement. We will now therefore explain why the occurrence of the high temperature phase of hadronic matter -the quark-gluon plasma-is still consistent with our incapability to liberate quarks in high energy collisions. It is thus important to realize that the currently accepted theory of hadronic structure and interactions, quantum chromodynamics [10], supplemented with its phenomenological extension, the MIT bag model [11], allows the formation of large space domains filled with (almost) free quarks. Such a state is expected to be unstable and to decay again into individual hadrons, following its free expansion. The mechanism of quark confinement requires that all quarks recombine to form hadrons again. Thus the quark-gluon plasma may be only a transitory form of hadronic matter formed under special conditions and therefore quite difficult to detect experimentally.
We will recall now the relevant postulates and results that characterize the current understanding of strong interactions in quantum chromodynamics (QCD). The most important postulate is that the proper vacuum state in QCD is not the (trivial) perturbative state that we (naively) imagine to exist everywhere and which is little changed when the interactions are turned on/off. In QCD, the true vacuum state is believed to a have a complicated structure which originates in the glue ("photon") sector of the theory. The perturbative vacuum is an excited state with an energy density B above the true vacuum. It is to be found inside hadrons where perturbative quanta of the theory, in particular quarks, can therefore exist. The occurrence of the true vacuum state is intimately connected to the glue-glue interaction. Unlike QED, these massless quanta of QCD, also carry a charge -color-that is responsible for the quark-quark interaction.
In the above discussion, the confinement of quarks is a natural feature of the hypothetical structure of the true vacuum. If it is, for example, a color superconductor, then an isolated charge cannot occur. Another way to look at this is to realize that a single colored object would, according to Gauss' theorem, have an electric field that can only end on other color charges. In the region penetrated by this field, the true vacuum is displaced, thus effectively raising the mass of a quasi-isolated quark by the amount BV field .
Another feature of the true vacuum is that it exercises a pressure on the surface of the region of the perturbative vacuum to which quarks are confined. Indeed, this is just the idea of the original MIT bag model [12]. The Fermi pressure of almost massless light quarks is in equilibrium with the vacuum pressure B. When many quarks are combined to form a giant quark bag, then their properties inside can be obtained using standard methods of many-body theory [2][3][4][5][6][7][8]. In particular, this also allows the inclusion of the effect of internal excitation through a finite temperature and through a change in the chemical composition.
A further effect that must be taken into consideration is the quark-quark interaction. We shall use here the first order contribution in the QCD running coupling constant α s (q 2 ) = g 2 /4π. However, as α s (q 2 ) increases when the average momentum exchanged between quarks decreases, this approach will have only limited validity at relatively low densities and/or temperatures. The collective screening effects in the plasma are of comparable order of magnitude and should reduce the importance of perturbative contributions as they seem to reduce the strength of the quark-quark interaction.
From this general description of the hadronic plasma, it is immediately apparent that, at a certain value of temperature and baryon number density, the plasma must disintegrate into individual hadrons. Clearly, to treat this process and the ensuing further nucleonisation by perturbative QCD methods is impossible. It is necessary to find a semi-phenomenological method for the treatment of the thermodynamic system consisting of a gas of quark bags. The hadronic gas phase is characterized by those reactions between individual hadrons that lead to the formation of new particles (quark bags) only. Thus one may view [13][14][15] the hadronic gas phase as being an assembly of many different hadronic resonances, their number in the interval (m 2 , m 2 +dm 2 ) being given by the mass spectrum τ (m 2 , b)dm 2 . Here the baryon number b is the only discrete quantum number to be considered at present. All bag-bag interaction is contained in the mutual transmutations from one state to another. Thus the gas phase has the characteristic of an infinite component ideal gas phase of extended objects. The quark bags having a finite size force us to formulate the theory of an extended, though otherwise ideal multicomponent gas.
It is a straightforward exercise, carried through in the beginning of the next section, to reduce the grand partition function Z to an expression in terms of the mass spectrum τ (m 2 , b). In principle, an experimental form of τ (m 2 , b) could then be used as an input. However, the more natural way is to introduce the statistical bootstrap model [13], which will provide us with a theoretical τ that is consistent with assumptions and approximations made in determining Z.
In the statistical bootstrap, the essential step consists in the realization that a composite state of many quark bags is in itself an "elementary" bag [1,16]. This leads directly to a nonlinear integral equation for τ . The ideas of the statistical bootstrap have found a very successful application in the description of hadronic reactions [17] over the past decade. The present work is an extension [1,15,18] and application [19] of this method to the case of a system containing any number of finite size hadronic clusters with their baryon numbers adding up to some fixed number. Among the most successful predictions of the statistical bootstrap, we record here the derivation of the limiting hadronic temperature and the exponential growth of the mass spectrum.
We see that the theoretical description of the two hadronic phases -the individual hadron gas and the quark-gluon plasma-is consistent with observations and with the present knowledge of elementary particles. What remains is the study of the possible phase transition between those phases as well as its observation. Unfortunately, we can argue that in the study of temperatures and mean transverse momenta of pions and nucleons produced in nuclear collisions, practically all information about the hot and dense phase of the collision is lost, as most of the emitted particles originate in the cooler and more dilute hadronic gas phase of matter. In order to obtain reliable information on quark matter, we must presumably perform more specific experiments. We will briefly point out that the presence of numerous s quarks in the quark plasma suggest, as a characteristic experiment, the observation Λ hyperons.
We close this report by showing that, in nuclear collisions, unlike pp reactions, we can use equilibrium thermodynamics in a large volume to compute the yield of strange and antistrange particles. The latter, e.g., Λ, might be significantly different from what one expects in pp collisions and give a hint about the properties of the quark-gluon phase.

Thermodynamics of the gas phase and the SBM
Given the grand partition function Z(β, V, λ) of a manybody system, all thermodynamic quantities can be determined by differentiation of ln Z with respect to its arguments. Here, λ is the fugacity introduced to conserve a discrete quantum number, here the baryon number. The conservation of strangeness can be carried through in a similar fashion leading then to a further argument λ s of Z. Whenever necessary, we will consider Z to be implicitly dependent on λ s . The grand partition function is a Laplace transform of the level density σ(p, V, b), where p μ is the four-momentum and b the baryon number of the many-body system enclosed in the volume V We recognize the usual relations for the thermodynamic expectation values of the baryon number and the energy-momentum four-vector which follow from the definition in eq. (2). The theoretical problem is to determine σ(p, V, b) in terms of known quantities. Let us suppose that the physical states of the hadronic gas phase can be considered as being built up from an arbitrary number of massive objects, henceforth called clusters, characterized by a mass spectrum τ (m 2 , b), where τ (m 2 , b)dm 2 is the number of different elementary objects (existing in nature) in the mass interval (m 2 , m 2 + dm 2 ) and having the baryon number b. As particle creation must be permitted, the number N of constituents is arbitrary, but constrained by four-momentum conservation and baryon conservation. Neglecting quantum statistics (it can be shown that, for T 40 MeV, Boltzmann statistics is sufficient), we have The sum over all allowed partitions of b into different b i is included and Δ is the volume available for the motion of the constituents, which differs from V if the different clusters carry their proper volume V ci The phase space volume used in eq. (4) is best explained by considering what happens for one particle of mass m 0 in the rest frame of Δ μ and β μ The density of states in eq. (4) implies that the creation and absorption of particles in kinetic and chemical equilibrium is limited only by four-momentum and baryon number conservation. These processes represent the strong hadronic interactions which are dominated by particle productions. τ (m 2 , b) contains all participating elementary particles and their resonances. Some remaining interaction is here neglected or, as we do not use the complete experimental τ , it may be considered as being taken care of by a suitable choice of τ . The short-range repulsive forces are taken into account by the introduction of the proper volume V of hadronic clusters.
One more remark concerning the available volume Δ is in order here. If V were considered to be given and an independent thermodynamic quantity, then in eq. (4), a further built-in restriction limits the sum over N to a certain N max , such that the available volume Δ in eq. (5) remains positive. However, this more conventional assumption of V as the independent variable would significantly obscure our mathematical formalism. It is important to realize that we are free to select the available volume Δ as the independent thermodynamic variable and to consider V as a thermodynamic expectation value to be computed from eq. (5) Here V μ c is the average sum of proper volumes of all hadronic clusters contained in the system considered. As already discussed, the standard quark bag leads to the proportionality between the cluster volume and hadron mass. Similar arguments within the bootstrap model [15], as for example discussed in the preceding lecture by R. Hagedorn [16], also lead to where 4B is the (at this point arbitrary) energy density of isolated hadrons in the quark bag model [11]. Since our hadrons are under pressure from neighbors in hadronic matter, we have in principle to take instead of 4B the energy density of a quark bag exposed to a pressure P (see eq. (54) below) ε bag = 4B + 3P.
Combining eqs. (7)-(9), we find, with ε(β, Δ, As we shall see, the pressure P in the hadronic matter never rises above 0.4B, see fig. 5a below, and arguments following eq. (60). Consequently, the inclusion of P above -the compression of free hadrons by the hadronic matter by about 10%-may be omitted for now from further discussion. However, we note that both ε and P will be computed as ln Z becomes available, whence eq. (9) is an implicit equation for Δ/ V . It is important to record that the expression in eq. (9) can approach zero only when the energy density of the hadronic gas approaches that of matter consisting of one big quark bag: ε → 4B, P → 0. Thus the density of states in eq. (4), together with the choice of Δ as a thermodynamic variable, is a consistent physical choice only up to this point. Beyond we assume that a description in terms of interacting quarks and gluons is the proper physical description. Bearing all these remarks in mind, we now consider the available volume Δ as a thermodynamic variable which by definition is positive. Inspecting eq. (4) again, we recognize that the level density of the extended objects in volume V can be interpreted for the time being as the level density of point particles in a fictitious volume Δ σ(p, V, b) = σ pt (p, Δ, b), (10) whence this is also true for the grand canonical partition function in eq. (2) Combining eqs. (2) and (4), we also find the important relation (12) This result can only be derived when the sum over N in eq. (4) extends to infinity, thus as long as Δ/ V in eq. (9) remains positive.
In order to continue with our description of hadronic matter, we must now determine a suitable mass spectrum τ to be inserted into eq. (4). For this we now introduce the statistical bootstrap model. The basic idea is rather old, but has undergone some development more recently making it clearer, more consistent, and perhaps more convincing. The details may be found in [15] and the references therein. Here a simplified naive presentation is given. We note, however, that our present interpretation is nontrivially different from that in [15].
The basic postulate of statistical bootstrap is that the mass spectrum τ (m 2 , b) containing all the "particles", i.e., elementary, bound states, and resonances (clusters), is generated by the same interactions which we see at work if we consider our thermodynamical system. Therefore, if we were to compress this system until it reaches its natural volume V c (m, b), then it would itself be almost a cluster appearing in the mass spectrum τ (m 2 , b). Since σ(p, Δ, b) and τ (p 2 , b) are both densities of states (with respect to the different parameters d 4 p and dm 2 ), we postulate that where = means "corresponds to" (in some way to be specified). As σ(p, Δ, b) is (see eq. (4)) the sum over N of Nfold convolutions of τ , the above "bootstrap postulate" will yield a highly nonlinear integral equation for τ . The bootstrap postulate (13) requires that τ should obey the equation resulting from replacing σ in eq. (4) by some expression containing τ linearly and by taking into account the volume condition expressed in eqs. (7) and (8).
We cannot simply put V = V c and Δ = 0, because now, when each cluster carries its own dynamically determined volume, Δ loses its original meaning and must be redefined more precisely. Therefore, in eq. (4), we tenta-tively replace Next we argue that the explicit factors m 2 and m 2 i arise from the dynamics and therefore must be absorbed into with where either H or m 0 may be taken as a new free parameter of the model, to be fixed later. (If m 0 is taken, then it should be of the order of the "elementary masses" appearing in the system, e.g., somwhere between m π and m N in a model using pions and nucleons as elementary input.) Finally, if clusters consist of clusters which consist of clusters, and so on, this should end at some "elementary" particles (where what we consider as elementary is fixed by convention). Inserting eq. (15) into eq. (4), the bootstrap equation (BE) then reads Clearly, the bootstrap equation (16) has not been derived.
We have made it more or less plausible and state it as a postulate. For more motivation, see [15]. In other words, the bootstrap equation means that the cluster with mass p 2 and baryon number b is either elementary (mass m b , spin isospin multiplicity g b ), or it is composed of any number N ≥ 2 of subclusters having the same internal composite structure described by this equation. The bar over m b indicates that one has to take the mass which the "elementary particle" will have effectively when present in a large cluster, e.g., in nuclear matter, m = m− E bind , and m N ≈ 925 MeV. That this must be so becomes obvious if one imagines eq. (16) solved by iteration (the iteration solution exists and is the physical solution). Then Hτ (p 2 , b) becomes in the end a complicated function of p 2 , b, all m b , and all g b . In other words, in the end a single cluster consists of the "elementary particles". As these are all bound into the cluster, their mass m should be the effective mass, not the free mass m. This way we may include a small correction for the long-range attractive meson exchange by choosing m N = m − 15 MeV.
Let us make a brief excursion to the bag model at this point. There the mass of a hadron is computed from the assumption of an isolated particle (= bag) with its size and mass being determined from the equilibrium between the vacuum pressure B and the internal Fermi pressure of the (valence) quarks. In a hadron gas, this is not true as a finite pressure is exerted on hadrons in matter. After a short calculation, we find the pressure dependence of the bag model hadronic mass We have already noted that the pressure never exceeds 0.4B in the hadronic gas phase, see fig. 5a below, and arguments following eq. (60). Hence we see that the increase in mass of constituents (quark bags) in the hadronic gas never exceeds 1.5% and is at most comparable with the 15 MeV binding in m. In general, P is about 0.1B and the pressure effect may be neglected.
Thus we can consider the "input" first term in eq. (16) as being fixed by pions, nucleons, and whenever necessary by the usual strange members of meson and baryon multiplets. Furthermore, we note that the bootstrap equation (16) makes use of practically all the same approximations as our description of the level density in eq. (4). Thus the solution of eq. (16) is particularly suitable for our use.
We solve the BE by the same double Laplace transformation which we used before eq. (2). We define Once the set of input particles is a known function, while Φ(β, λ) is unknown. Applying the double Laplace transformation to the BE, we obtain This implicit equation for Φ in terms of ϕ can be solved without regard for the actual β, λ dependence. Writing we can draw the curve ϕ(G) and then invert it graphically (see fig. 1) to obtain G(ϕ) = Φ(β, λ). G(ϕ) has a square root singularity at ϕ = ϕ 0 = ln(4/e) = 0.3863. Beyond this value, G(ϕ) becomes complex. Apart from this graphical solution, other forms of solution are known: The expansion in terms of (ϕ 0 −ϕ) n/2 has been used in our numerical work (12 terms yield a solution within computer accuracy) and the integral representation will be published elsewhere 2 . Henceforth, we consider Φ(β, λ) = G(ϕ) to be a known function of ϕ(β, λ). Consequently, τ (m 2 , b) is also in principle known. From the singularity at ϕ = ϕ 0 , it follows [1] that τ (m 2 , b) grows, for m m N b, exponentially ∼ m −3 exp(m/T 0 ). In some weaker form, this has been known for a long time [13,20,21].

The hot hadronic gas
The definition of Φ(β, λ) in eq. (18) in terms of the mass spectrum allows us to write a very simple expression for ln Z in the gas phase (passing now to the rest frame of the gas) We recall that eqs. (9) and (19) define (implicitly) the quantities Δ and Φ in terms of the physical variables V , β, and λ. Let us now introduce the energy density ε pt of the hypothetical pointlike particles as which will turn out to be quite helpful as it is independent of Δ. The proper energy density is while the pressure follows from Similarly, for the baryon number density, we find with (28) From eqs. (23)-(23), the crucial role played by the factor Δ/ V becomes apparent. We note that it is quite straightforward to insert eqs. (24) and (25) into eq. (9) and solve the resulting quadratic equation to obtain Δ/ V as an explicit function of ε pt and P pt . First we record the limit P B while the correct expression is (30) The last of the important thermodynamic quantities is the entropy S. By differentiating eq. (25), we find Considering Z as a function of the chemical potential, viz., with E being the total energy. From eqs. (31) and (33), we find the "first law" of thermodynamics to be Now quite generally Equations (25) and (33) now allow us to write The entropy density in terms of the already defined quantities is therefore We shall now take a brief look at the quantities P , ε, ν, Δ/ V . They can be written in terms of ∂Φ(β, λ)/∂β and its derivatives. We note that (see eq. (20)) and that ∂G/∂ϕ ∼ (ϕ 0 − ϕ) −1/2 near to ϕ = ϕ 0 = ln(4/e) (see fig. 1). Hence at ϕ = ϕ 0 , we find a singularity in the point particle quantities ε pt , ν pt , and P pt . This implies that all hadrons have coalesced into one large cluster. Indeed, from eqs. (24), (26), (27), and (29), we find We can easily verify that this is correct by establishing the average number of clusters present in the hadronic gas. This is done by introducing an artificial fugacity ξ N in eq. (4) in the sum over N , where N is the number of clusters. Denoting by Z(ξ) the associated grand canonical partition functions in eq. (22), we find which leads to the useful relation Thus as P V → 0, so must N , the number of clusters, for finite T . We record the astonishing fact that the hadron gas phase obeys an "ideal" gas equation, although of course N is not constant as for a real ideal gas but a function of the thermodynamic variables. The boundary given by thus defines a critical curve in the β, λ plane. Its position depends, of course, on the actually given form of ϕ(β, λ), i.e., on the set of "input" particles {m b , g b } assumed and the value of the constant H in eq. (15). In the case of three elementary pions π + , π 0 , and π − and four elementary nucleons (spin ⊗ isospin) and four antinucleons, we have from eq. (18) ϕ(β, λ) = 2πHT 3m π K 1 (m π /T ) and the condition (42), written in terms of T and μ = T ln λ, yields the curve shown in fig. 2, i.e., the "critical curve". For μ = 0, the curve ends at T = T 0 , where T 0 , the "limiting temperature of hadronic matter", is the same as that appearing in the mass spectrum [13,15,20,21] . The value of the constant H in eq. (15) has been chosen [19] to yield T 0 = 190 MeV. This apparently large value of T 0 seemed necessary to yield a maximal average decay temperature of the order of 145 MeV, as required by [22].
where the value of m 0 lies as expected between m π and m N ((m π m N ) 1/2 = 0.36 GeV). The critical curve limits the hadron gas phase. By approaching it, all hadrons dissolve into a giant cluster, which is not in our opinion a hadron solid [23]. We would prefer to identify it with a quark-gluon plasma. Indeed, as the energy density along the critical curve is constant (= 4B), the critical curve can be attained and, if the energy density becomes > 4B, we enter into a region which cannot be described without making assumptions about the inner structure and dynamics of the "elementary particles" {m b , g b } -here pions and nucleons-entering into the input function ϕ(β, λ). Considering pions and nucleons as quark-gluon bags leads naturally to this interpretation.

The quark-gluon phase
We now turn to the discussion of the region of the strongly interacting matter in which the energy density would be equal to or higher than 4B. As a basic postulate, we will assume that it consists of -relatively weakly-interacting quarks. To begin with, only u and d flavors will be considered as they can easily be copiously produced at T 50 MeV. Again the aim is to derive the grand partition function Z. This is a standard exercise. For the massless quark Fermi gas up to first order in the interaction [1][2][3][4][5][6][7][8]18], the result is valid in the limit m q < T ln λ q .
Here g = (2s + 1)(2I + 1)C = 12 counts the number of the components of the quark gas, and λ q is the fugacity related to the quark number. As each quark has baryon number 1/3, we find where as before λ allows for conservation of the baryon number. Consequently The glue contribution is We notice the two relevant differences with the photon gas: -The occurrence of the factor eight associated with the number of gluons.
-The glue-glue interaction as gluons carry color charge.
Finally, let us introduce the vacuum term, which accounts for the fact that the perturbative vacuum is an excited state of the "true" vacuum which has been renormalized to have a vanishing thermodynamic potential, Ω = −β −1 ln Z. Hence in the perturbative vacuum This leads to the required positive energy density B within the volume occupied by the colored quarks and gluons and to a negative pressure on the surface of this region. At this stage, this term is entirely phenomenological, as discussed above. The equations of state for the quark-gluon plasma are easily obtained by differentiating with respect to β, λ, and V . The baryon number density, energy, and pressure are, respectively Let us first note that, for T μ and P = 0, the baryon chemical potential tends to which assures us that interacting cold quark matter is an excited state of nuclear matter. We have assumed that, except for T , there is no relevant dimensional parameter, e.g., quark mass m q or the quantity Λ which enters into the running coupling constant α s (q 2 ). Therefore the relativistic relation between the energy density and pressure, viz., ε − B = 3(P + B), is preserved, which leads to a relation we have used occasionally before (see eq. (9)). From eq. (54), it follows that, when the pressure vanishes, the energy density is 4B, independent of the values of μ and T which fix the line P = 0. This behavior is consistent with the hadronic gas phase. This may be used as a reason to choose the parameters of both phases in such a way that the two lines P = 0 coincide. We will return to this point again below. For P > 0, we have ε > 4B. Recall that, in the hadronic gas, we had 0 < ε < 4B. Thus, above the critical curve of the μ, T plane, we have the quark-gluon plasma exposed to an external force. In order to obtain an idea of the form of the P = 0 critical curve in the μ, T plane for the quark-gluon plasma, we rewrite eq. (52) using eqs. (45) and (46) for P = 0 Here, the last term is the glue pressure contribution. (If the true vacuum structure is determined by the glueglue interaction, then this term could be modified significantly.) We find that the greatest lower bound on temperature T q at μ = 0 is about This result can be considered to be correct to within 20%. Its order of magnitude is as expected. Taking eq. (55) as it is, we find for α s = 1/2, T q = 0.88B 1/4 . Omitting the gluon contribution to the pressure, we find T q = 0.9B 1/4 . It is quite likely that, with the proper treatment of the glue field and the plasma corrections, and with larger B 1/4 ∼ 190 MeV, the desired value of T q = T 0 corresponding to the statistical bootstrap choice will follow. Furthermore, allowing some reasonable T , μ dependence of α s , we can then easily obtain an agreement between the critical curves. However, it is not necessary for the two critical curves to coincide, even though this would be preferable. As the quark plasma is the phase into which individual hadrons dissolve, it is sufficient if the quark plasma pressure vanishes within the boundary set for non-vanishing positive pressure of the hadronic gas. It is quite satisfactory for the theoretical development that this is the case. In fig. 3a, a qualitative picture of the two P = 0 lines is shown in the μ, T plane. Along the dotted straight line at constant temperature, we show in fig. 3b the pressure as a function of the volume (a P , V diagram). The volume is obtained by inverting the baryon density at constant fixed baryon number The behavior of P (V, T = const.) for the hadronic gas phase is as described before in the statistical bootstrap model. For large volumes, we see that P falls with rising V . However, when hadrons get close to each other so that they form larger and larger lumps, the pressure drops rapidly to zero. The hadronic gas becomes a state of few composite clusters (internally already consisting of the quark plasma). The second branch of the P (V, T = const.) line meets the first one at a certain volume V = V m . The phase transition occurs for T = const. in fig. 3b at a vapor pressure P v obtained from the conventional Maxwell construction: the shaded regions in fig. 3b are equal. Between the volumes V 1 and V 2 , matter coexists in the two phases with the relative fractions being determined by the magnitude of the actual volume. This leads to the occurrence of a third region, viz., the coexistence region of matter, in addition to the pure quark and hadron domains. For V < V 1 , corresponding to ν > ν 1 ∼ 1/V 1 , all matter has gone into the quark plasma phase.
The dotted line in fig. 3b encloses (qualitatively) the domain in which the coexistence between the two phases of hadronic matter seems possible. We further note that, at low temperatures T ≤ 50 MeV, the plasma and hadronic gas critical curves meet each other in fig. 3a. This is just the domain where, at present, our description of the hadronic gas fails, while the quark-gluon plasma also begins to suffer from infrared difficulties. Both approaches have a very limited validity in this domain.
The qualitative discussion presented above can be easily supplemented with quantitative results. But first we turn our attention to the modifications forced onto this simple picture by the experimental circumstances in high energy nuclear collisions. Fig. 4. a) The critical curve of hadron matter (bootstrap), together with some "cooling curves" in the T , μ plane. While the system cools down along these lines, it emits particles. When all particles have become free, it comes to rest on some point on these curves ("freeze out"). In the shaded region, our approach may be invalid. b) The critical curve of hadron matter (bootstrap), together with some "cooling curves" (same energy as in fig. 4a) in the variables T and ν/ν 0 = ratio of baryon number density to normal nuclear baryon number density. In the shaded region, our approach may be invalid.

Nuclear collisions and inclusive particle spectra
We assume that in relativistic collisions triggered to small impact parameters by high multiplicities and absence of projectile fragments [24], a hot central fireball of hadronic matter can be produced. We are aware of the whole problematic connected with such an idealization. A proper treatment should include collective motions and distribution of collective velocities, local temperatures, and so on [25][26][27][28], as explained in the lecture by R. Hagedorn [16]. Triggering for high multiplicities hopefully eliminates some of the complications. In nearly symmetric collisions (projectile and target nuclei are similar), we can argue that the numbers of participants in the center of mass of the fireball originating in the projectile or target are the same. Therefore, it is irrelevant how many nucleons do form the fireball -and the above symmetry argument leads, in a straightforward way, to a formula for the center of mass energy per participating nucleon where E k,lab /A is the projectile kinetic energy per nucleon in the laboratory frame. While the fireball changes its baryon density and chemical composition (π + p ↔ Δ, etc.) during its lifetime through a change in temperature and chemical potential, the conservation of energy and baryon number assures us that U in eq. (58) remains constant, assuming that the influence on U of pre-equilibrium emission of hadrons from the fireball is negligible. As U is the total energy per baryon available, we can, supposing that kinetic and chemical equilibrium have been reached, set it equal to the ratio of thermodynamic expectation values of the total energy and baryon number Thus we see that, through eq. (59), the experimental value of U in eq. (58) fixes a relation between allowable values of β, λ: the available excitation energy defines the temperature and the chemical composition of hadronic fireballs. In fig. 4a and b, these paths are shown for a choice of kinetic energies E k,lab /A in the μ, T plane and in the ν, T plane, respectively. In both cases, only the hadronic gas domain is shown. We wish to note several features of the curves shown in fig. 4a and b that will be relevant in later considerations: 1) Beginning at the critical curve, the chemical potential first drops rapidly when T decreases and then rises slowly as T decreases further ( fig. 4a). This corresponds to a monotonically falling baryon density with decreasing temperature (fig. 4b), but implies that, in the initial expansion phase of the fireball, the chemical composition changes more rapidly than the temperature.
2) The baryon density in fig. 4b is of the order of 1-1.5 of normal nuclear density. This is a consequence of the choice of B 1/4 = 145 MeV. Were B three times as large, i.e., B 1/4 = 190 MeV, which is so far not excluded, then the baryon densities in this figure would triple to 3-5ν 0 . Furthermore, we observe that, along the critical curve of the hadronic gas, the baryon density falls with rising temperature. This is easily understood as, at higher temperature, more volume is taken up by the numerous mesons. 3) Inspecting fig. 4b, we see that, at given U , the temperatures at the critical curve and those at about ν 0 /2 differ little (10%) for low U , but more significantly for large U . Thus, highly excited fireballs cool down more before dissociation ("freeze out"). As particles are emitted all the time while the fireball cools down along the lines of fig. 4a and b, they carry kinetic energies related to various different temperatures. The inclusive single particle momentum distribution will yield only averages along these cooling lines.
Another remark which does not follow from the curves shown is: 4) Below about 1.8 GeV, an important portion of the total energy is in the collective (hydrodynamical) motion of hadronic matter, hence the cooling curves at constant excitation energy do not properly describe the evolution of the fireball.
Calculations of this kind can also be carried out for the quark plasma. They are, at present, uncertain due to the unknown values of α s and B 1/4 . Fortunately, there is one particular property of the equation of state of the quarkgluon plasma that we can easily exploit. Combining eq. (54) with eq. (59), we obtain Thus, for a given U (the available energy per baryon in a heavy ion collision), eq. (60) describes the pressure-volume (∼ 1/ν) relation. By choosing to measure P in units of B and ν in units of normal nuclear density ν 0 = 0.14/ fm 3 , we find with γ := m N ν 0 4B = 0.56, for : B 1/4 = 145 MeV, ν 0 = 0.14/fm 3 .
Here, γ is the ratio of the energy density of normal nuclei (ε N = m N ν 0 ) and of quark matter or of a quark bag (ε q = 4B). In fig. 5a, this relation is shown for three projectile energies: E k,lab /A = 1.80 GeV, 3.965 GeV, and 5.914 GeV, corresponding to U = 1.314 GeV, 1.656 GeV, and 1.913 GeV, respectively. We observe that, even at the lowest energy shown, the quark pressure is zero near the baryon density corresponding to 1.3 normal nuclear density, given the current value of B. Before discussing this point further, we note that the hadronic gas branches of the curves in fig. 5a and b show a quite similar behavior to that shown at constant temperature in fig. 3b. Remarkably, the two branches meet each other at P = 0, since both have the same energy density ε = 4B and therefore V (P = 0) ∼ 1/ν = U/ε = U/4B. However, what we cannot see by inspecting fig. 5a and b is that there will be a discontinuity in the variables μ and T at this point, except if parameters are chosen so that the critical curves of the two phases coincide. Indeed, near to P = 0, the results shown in fig. 5a should be replaced by points obtained from the Maxwell construction. The pressure in a nuclear collision will never fall to zero. It will correspond to the momentary vapor pressure of the order of 0.2B as the phase change occurs.
A further aspect of the equations of state for the hadronic gas is also illustrated in fig. 5a. Had we ignored the finite size of hadrons (one of the van der Waals effects) in the hadron gas phase then, as shown by the dashdotted lines, the phase change could never occur because the point particle pressure would diverge where the quark pressure vanishes. In our opinion, one cannot say it often enough: inclusion of the finite hadronic size and of the finite temperature when considering the phase transition to quark plasma lowers the relevant baryon density (from 8-14ν 0 for cold point-nucleon matter) to 1-5ν 0 (depending on the choice of B) in 2-5 GeV/A nuclear collisions. The possible formation of quark-gluon plasma in nuclear collisions was first discussed quantitatively in ref. [3], see also ref. [29].
The physical picture underlying our discussion is an explosion of the fireball into vacuum with little energy being converted into collective motion, e.g., hydrodynamical flow, or being taken away by fast pre-hadronization particle emission. Thus the conserved internal excitation energy can only be shifted between thermal (kinetic) and chemical excitations of matter. "Cooling" thus really means that, during the explosion, the thermal energy is mostly convered into chemical energy, e.g., pions are produced.
While it is at present hard to judge the precise amount of expected deviation from the cooling curves shown in fig. 2, it is possible to show that they are entirely inconsistent with the notion of reversible adiabatic, i.e., entropy conserving, expansion. As the expansion proceeds along U = const. lines, we can compute the entropy per participating baryon using eqs. (36) and (37), and we find a significant growth of total entropy. As shown in fig. 5b, the entropy rises initially in the dense phase of the matter by as much as 50-100% due to the pion production and resonance decay. Amusingly enough, as the newly produced entropy is carried mostly by pions, one will find that the entropy carried by protons remains constant. With this remarkable behavior of the entropy, we are in a certain sense, victims of our elaborate theory. Had we used, e.g., an ideal gas of Fermi nucleons, then the expansion would seem to be entropy conserving, as pion production and other chemistry were forgotten. Our fireballs have no tendency to expand reversibly and adiabatically, as many reaction channels are open. A more complete discussion of the entropy puzzle can be found in [1].
Inspecting fig. 4a and b again, it seems that a possible test of the equations of state for the hadronic gas consists in measuring the temperature in the hot fireball zone, and doing this as a function of the nuclear collision energy. The plausible assumption made is that the fireball follows the "cooling" lines shown in fig. 4a and b until final dissociation into hadrons. This presupposes that the surface emission of hadrons during the expansion of the fireball does not significantly alter the available energy per baryon. This is more likely true for sufficiently large fireballs. For small ones, pion emission by the surface may influence the energy balance. As the fireball expands, the temperature falls and the chemical composition changes. The hadronic clusters dissociate and more and more hadrons are to be found in the "elementary" form of a nucleon or a pion. Their kinetic energies are reminiscent of the temperature found at each phase of the expansion.
To compute the experimentally observable final temperature [1,19], we shall argue that a time average must be performed along the cooling curves. Not knowing the reaction mechanisms too well, we assume that the temperature decreases approximately linearly with the time in the significant expansion phase. We further have to al-low that a fraction of particles emitted can be reabsorbed in the hadronic cluster. This is a geometric problem and, in a first approximation, the ratio of the available volume Δ to the external volume V ex is the probability that an emitted particle not be reabsorbed, i.e., that it can escape The relative emission rate is just the integrated momentum spectrum The chemical potential acts only for nucleons. In the case of pions, it has to be dropped from the above expression. For the mean temperature, we thus find where the subscript "c" on the integral indicates here a line integral along that particular cooling curve in fig. 4a and b which belongs to the energy per baryon fixed by the experimentalist.
In practice, the temperature is most reliably measured through the measurement of mean transverse momenta of the particles. It may be more practical therefore to calculate the average transverse momentum of the emitted particles. In principle, to obtain this result we have to perform a similar averaging to the one above. For the average transverse momentum at given T, μ, we find [14] p ⊥ (m, T, μ) p = The average over the cooling curve is then (66) We did verify numerically that the order of averages does not matter which shows that the mean transverse momentum is also the simplest (and safest) method of determining the average temperature (indeed better than fitting ad hoc exponential type functions to p ⊥ distributions).  . 4a). The mean temperatures are obtained by integrating along the cooling curves. Note that T N is always greater than Tπ.
In the presented calculations, we chose the bag constant B = (145 MeV) 4 , but we now believe that a larger B should be used. As a consequence of our choice and the measured pion temperature of T ex π = 140 MeV at highest ISR energies, we have to choose the constant H such that T 0 = 190 MeV (see eq. (43b)).
The average temperature, as a function of the range of integration over T , reaches different limiting values for different particles. The limiting value obtained thus is the observable "average temperature" of the debris of the interaction, while the initial temperature T cr at given E k,lab (full line in fig. 6) is difficult to observe. When integrating along the cooling line as in eq. (64), we can easily, at each point, determine the average hadronic cluster mass. The integration for protons is interrupted (protons are "frozen out") when the average cluster mass is about half the nucleon isobar mass. We have also considered baryon density dependent freeze-out, but such a procedure depends strongly on the unreliable value of B.
Our choice of the freeze-out condition was made in such a way that the nucleon temperature at E k,lab /A = 1.8 GeV is about 120 MeV. The model dependence of our freeze-out introduces an uncertainty of several MeV in the average temperature. In fig. 6, the pion and nucleon average temperatures are shown as a function of the heavy ion kinetic energy. Two effects contributed to the difference between the π and N temperatures: 1) The particular shape of the cooling curves ( fig. 4a).
The chemical potential drops rapidly from the critical curve, thereby damping relative baryon emission at lower T . Pions, which do not feel the baryon chemical potential, continue being created also at lower temperatures.
2) The freeze-out of baryons occurs earlier than the freeze-out of pions.
A third effect has been so far omitted -the emission of pions from two-body decay of long-lived resonances [1] would lead to an effective temperature which is lower in nuclear collisions. In fig. 7, we show the dependence of the average transverse momenta of pions and nucleons on the kinetic energy of the heavy ion projectiles.

Strangeness in heavy ion collisions
From the averaging process described here, we have learned that the temperatures and transverse momenta of particles originating in the hot fireballs are more reminiscent of the entire history of the fireball expansion than of the initial hot compressed state, perhaps present in the form of quark matter. We may generalize this result and then claim that most properties of inclusive spectra are reminiscent of the equations of state of the hadronic gas phase and that the memory of the initial dense state is lost during the expansion of the fireballs as the hadronic gas rescatters many times while it evolves into the final kinetic and chemical equilibrium state.
In order to observe properties of quark-gluon plasma, we must design a thermometer, an isolated degree of freedom weakly coupled to the hadronic matter. Nature has, in principle (but not in practice) provided several such thermometers: leptons and heavy flavors of quarks. We would like to point here to a particular phenomenon perhaps quite uniquely characteristic of quark matter. First we note that, at a given temperature, the quark-gluon plasma will contain an equal number of strange (s) quarks and antistrange (s) quarks, naturally assuming that the hadronic collision time is much too short to allow for light flavor weak interaction conversion to strangeness. Thus, assuming equilibrium in the quark plasma, we find the density of the strange quarks to be (two spins and three colors) neglecting for the time being the perturbative corrections and, of course, ignoring weak decays. As the mass m s of the strange quarks in the perturbative vacuum is believed to be of the order of 280-300 MeV, the assumption of equilibrium for m s /T ∼ 2 may indeed be correct. In eq. (68), we were able to use the Boltzmann distribution again, as the density of strangeness is relatively low. Similarly, there is a certain light antiquark density (q stands for either u or d ) where the quark chemical potential is μ q = μ/3, as given by eq. (46). This exponent suppresses the qq pair production.
What we intend to show is that there are many more s quarks than antiquarks of each light flavor. Indeed, The function x 2 K 2 (x) is, for example, tabulated in [30]. For x = m s /T between 1.5 and 2, it varies between 1.3 and 1. Thus, we almost always have more s than q quarks and, in many cases of interest, s/q ∼ 5. As μ → 0, there are about as many u and q quarks as there are s quarks. When the quark matter dissociates into hadrons, some of the numerous s may, instead of being bound in a qs kaon, enter into a q q s antibaryon and, in particular 3 , a Λ or Σ 0 . The probability for this process seems to be comparable to the similar one for the production of antinucleons by the antiquarks present in the plasma. What is particularly noteworthy about the s-carrying antibaryons is that they can conventionally only be produced in direct pair production reactions. Up to about E k,lab /A = 3.5 GeV, this process is very strongly suppressed by energy-momentum conservation because, for free pp collisions, the threshold is at about 7 GeV. We would thus like to argue that a study of the Λ and Σ 0 in nuclear collisions for 2 < E k,lab /A < 4 GeV could shed light on the early stages of the nuclear collisions in which quark matter may be formed.
Let us mention here another effect of importance in this context: the production rate of a pair of particles with a conserved quantum number like strangeness will usually be suppressed by the Boltzmann factor e −2m/T , rather than a factor e −m/T as is the case in thermomechanical equilibrium (see, for example, the addendum in [14]). As relativistic nuclear collisions are just on the borderline between those two limiting cases, it is important when considering the yield of strange particles to understand the transition between them. We will now show how one can describe these different cases in a unified statistical description [31].
As we have already implicitly discussed (see eq. (12)), the logarithm of the grand partition function Z is a sum 3 Σ 0 decays into Λ by emitting a photon and is always counted within the Λ abundance. over all different particle configurations, e.g., expressed with the help of the mass spectrum. Hence, we can now concentrate in particular on that part of ln Z which is exclusively associated with the strangeness.
As the temperatures of interest to us and which allow appreciable strangeness production are at the same time high enough to prevent the strange particles from being thermodynamically degenerate, we can restrict ourselves again to the discussion of Boltzmann statistics only.
The contribution to Z of a state with k strange particles is where the one-particle function Z 1 for a particle of mass m s is given in eq. (16). To include both particles and antiparticles as two thermodynamically independent phases in eq. (71), the sum over s in eq. (71) must include them both. As the quantum numbers of particles (p) and antiparticles (a) must always be present with exactly the same total number, not each term in eq. (71) can contribute. Only when n = k/2 = number of particles = number of antiparticles is exactly fulfilled do we have a physical state. Hence, We now introduce the fugacity factor f n to be able to count the number of strange pairs present. Allowing an arbitrary number of pairs to be produced, we obtain where I 0 is the modified Bessel function and We have to maintain the difference between the particles (p) and antiparticles (a), as in nuclear collisions the symmetry is broken by the presence of baryons and there is an associated need for a baryon fugacity (chemical potential μ) that controls the baryon number. We obtain for particles (+μ) and antiparticles (−μ), where W (x) = x 2 K 2 (x), x i = m i /T , and all kaons and hyperons are counted. In the quark phase, we have with T x s = m s ∼ 280 MeV. We note in passing that the baryon chemical potential cancels out in y of eq. (74) when eq. (76) is inserted in the quark phase (compare with eq. (68)). By differentiating ln Z s of eq. (73) with respect to f , we find the strangeness number present at given T and V For large y, that is, at given T for large volume V , we find n s = √ y ∼ e −m/T , as expected. For small y, we find n s = y ∼ e −2m/T . In fig. 8, we show the dependence of the quenching factor I 1 /I 0 = η as a function of the volume V measured in units of V h = 4π/3 fm 3 for a typical set of parameters: T = 150, μ = 550 MeV (hadronic gas phase).
The following observations follow from inspection of fig. 8: 1) The strangeness yield is a qualitative measure of the hadronic volume in thermodynamic equilibrium. 2) Total strangeness yield is not an indicator of the phase transition to quark plasma, as the enhancement ( η q /η = 1.25) in yield can be reinterpreted as being due to a change in hadronic volume. 3) We can expect that, in nuclear collisions, the active volume will be sufficiently large to allow the strangeness yield to correspond to that of "infinite" volume for reactions triggered on "central collisions". Hence, e.g., Λ production rate will significantly exceed that found in pp collisions.
Our conclusions about the significance of Λ as an indicator of the phase transition to quark plasma remain valid as the production of Λ in the hadronic gas phase will only be possible in the very first stages of the nuclear collisions, if sufficient center of mass energy is available.

Summary
Our aim has been to obtain a description of hadronic matter valid for high internal excitations. By postulating the kinetic and chemical equilibrium, we have been able to develop a thermodynamic description valid for high temperatures and different chemical compositions. In our work we have found two physically different domains: firstly, the hadronic gas phase, in which individual hadrons can exist as separate entities, but are sometimes combined into larger hadronic clusters, while in the second domain, individual hadrons dissolve into one large cluster consisting of hadronic constituents, viz., the quark-gluon plasma.
In order to obtain a theoretical description of both phases, we have used some "common" knowledge and plausible interpretations of currently available experimental observations. In particular, in the case of hadronic gas, we have completely abandoned a more conventional Lagrangian approach in favour of a semi-phenomenological statistical bootstrap model of hadronic matter that incorporates those properties of hadronic interaction that are, in our opinion, most important in nuclear collisions.
In particular, the attractive interactions are included through the rich, exponentially growing hadronic mass spectrum τ (m 2 , b), while the introduction of the finite volume of each hadron is responsible for an effective shortrange repulsion. Aside from these manifestations of strong interactions, we only satisfy the usual conservation laws of energy, momentum, and baryon number. We neglect quantum statistics since quantitative study has revealed that this is allowed above T ≈ 50 MeV. But we allow particle production, which introduces a quantum physical aspect into the otherwise "classical" theory of Boltzmann particles.
Our approach leads us to the equations of state of hadronic matter which reflect what we have included in our considerations. It is the quantitative nature of our work that allows a detailed comparison with experiment. This work has just begun and it is too early to say if the features of strong interactions that we have chosen to include in our considerations are the most relevant ones. It is important to observe that the currently predicted pion and nucleon mean transverse momenta and temperatures show the required substantial rise (see fig. 7) as required by the experimental results available at E k,lab /A = 2 GeV (BE-VALAC, see [24]) and at 1000 GeV (ISR, see [22]). Further comparisons involving, in particular, particle multiplicities and strangeness production are under consideration.
We also mention the internal theoretical consistency of our two-fold approach. With the proper interpretation, the statistical bootstrap leads us, in a straightforward fashion, to the postulate of a phase transition to the quark-gluon plasma. This second phase is treated by a quite different method. In addition to the standard Lagrangian quantum field theory of weakly interacting particles at finite temperature and density, we also introduce the phenomenological vacuum pressure and energy density B.
Perhaps the most interesting aspect of our work is the realization that the transition to quark matter will occur at much lower baryon density for highly excited hadronic matter than for matter in the ground state (T = 0). The precise baryon density of the phase transition depends somewhat on the bag constant, but we estimate it to be at about 2-4ν 0 at T = 150 MeV. The detailed study of the different aspects of this phase transition, as well as of possible characteristic signatures of quark matter, must still be carried out. We have given here only a very preliminary report on the status of our present understanding.
We believe that the occurrence of the quark plasma phase is observable and we have proposed therefore a measurement of theΛ/p relative yield between 2 and 10 GeV/N kinetic energies. In the quark plasma phase, we expect a significant enhancement ofΛ production which will most likely be visible in theΛ/p relative rate.
1980: Many fruitful discussions with the GSI/LBL Relativistic Heavy Ion group stimulated the ideas presented here. I would like to thank R. Bock and R. Stock for their hospitality at GSI during this workshop. As emphasized before, this work was performed in collaboration with R. Hagedorn. This work was in part supported by Deutsche Forschungsgemeinschaft. 2015: Also in part supported by the US Department of Energy, Office of Science, Office of Nuclear Physics under award number DE-FG02-04ER41318.

Open Access
This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.