Mott-Hagedorn Resonance Gas and Lattice QCD Results

A compound effective model reproducing the equation of state of hadronic matter as obtained in recent lattice QCD simulations is presented.


Introduction
We are going to construct 1 a combined effective model reproducing the equation of state of hadronic matter as obtained in recent lattice QCD simulations [2,3]. The model should reproduce basic physical characteristics of processes encountered in the dense hadronic matter, from the hot QCD phase through the critical temperature region till the lower temperature hadron resonance gas phase. It has been shown in [4] that the equation of state derived from that time QCD lattice calculation [5] can be reproduced by a simple hadron gas resonance model below critical temperature T c . For higher temperatures the model is modified by introducing finite widths of heavy hadrons [6] with a heuristic ansatz for the spectral function which reflects medium modifications of hadrons. This Mott-Hagedorn type model is constructed to fit nicely the lattice data, also above T c where it does so because it leaves light hadrons below a mass threshold of m 0 = 1 GeV unaffected. The description of the lattice data at high temperatures is accidental because the effective number of those degrees of freedom approximately coincides with that of quarks and gluons.
In order to remove this unphysical aspect of the otherwise appealing model one has to extend the spectral broadening also to the light hadrons and thus describe their disappearance due to the Mott effect while simultaneously the quark and gluon degrees of freedom appear at high temperatures due to chiral symmetry restoration and deconfinement. In the present contribution we will report first results obtained by introducing a unified treatment of all hadronic resonances with a state-dependent width Γ i (T ) in accordance with the inverse hadronic collision time scale from a recent model for chemical freeze-out in a resonance gas [7]. The appearance of quark and gluon degrees of freedom is introduced by the Polyakov-loop improved Nambu-Jona-Lasinio (PNJL) model [8,9]. The model is further refined by adding perturbative corrections to O(α s ) for the high-momentum region above the three-momentum cutoff inherent in the PNJL model. One obtains eventually a good agreement with lattice QCD data, comparable with all important physical characteristics taken into account.

Extended Mott-Hagedorn resonance gas
We introduce the width Γ of a resonance in the statistical model through the spectral function where N m is the standard normalization factor. The model ansatz for the resonance width Γ is given by [6] where C Γ = 10 −4 , N m = 2.5, N T = 6.5 and T H = 165 MeV. For simplicity, we assume n S = 0 for the strangeness number density and n B = 0 for the baryon number density. Then the respective chemical potentials µ B = 0 and µ S = 0 always, so that the temperature is the only significant statistical parameter here.
The energy density of this model can then be cast in the form where m 0 = 1 GeV and the energy density per degree of freedom with a mass M is with the degeneracy g i . For mesons, δ i = −1 and for baryons δ i = 1. For our hadronic gas kept at fixed volume we have the relation where P = P (T ) and s = s(T ) are the pressure and entropy density.
Lattice data: Borsanyi et al. In the left panel of Fig. 1 the results for the pressure and energy density of the model at this stage are shown. Although providing us with an excellent fit of the lattice data, the high-temperature phase of this model is unphysical. Imposing that all mesons lighter than m 0 = 1 GeV are stable provides us with a SB limit at high temperatures which mimicks that due to quarks and gluons in the case for three flavors [10]. In reality, due to the chiral phase transition at T c , the quarks loose their mass and therefore the threshold of the continuum of quark-antiquark scattering states is lowered. The light meson masses, however, remain almost unaffected by the increase in the temperature of the system. Consequently, they merge the continuum and become unbound -their spectral function changes from a delta-function (on-shell bound states) to a Breit-Wigner-type (off-shell, resonant scattering states). This phenomenon is the hadronic analogue [11] of the Mott-Anderson transition for electrons in solid state physics (insulatorconductor transition). It has been first introduced for the hadronic-to-quarkmatter transition in [12]. Later, within the NJL model, a microscopic approach to the thermodynamics of the Mott dissociation of mesons in quark matter has been given in the form of a generalized Beth-Uhlenbeck equation of state [13], see also [14].
As a microscopic treatment of the Mott effect for all resonances is presently out of reach, we introduce an ansatz for a state-dependent hadron resonance width Γ i (T ) given by the inverse collision time scale recently suggested within an approach to the chemical freeze-out and chiral condensate in a resonance gas [7] Γ which is based on a binary collision approximation and relaxation time ansatz using for the in-medium hadron-hadron cross sections the geometrical Povh-Hüfner law [15]. In Eq. (6) the coefficient λ is a free parameter, n j (T ) is the partial density of the hadron j and the mean squared radii of hadrons r 2 i T obtain in the medium a temperature dependence which is governed by the (partial) restoration of chiral symmetry. For the pion this was quantitatively studied within the NJL model [16] and it was shown that close to the Mott transition the pion radius is well approximated by Here the Gell-Mann-Oakes-Renner relation has been used and the pion mass shall be assumed chirally protected and thus temperature independent.
For the nucleon, we shall assume the radius to consist of two components, a medium independent hard core radius r 0 and a pion cloud contribution r 2 N (T ) = r 2 0 + r 2 π (T ) , where from the vacuum values r π = 0.59 fm and r N = 0.74 fm one gets r 0 = 0.45 fm. A key point of our approach is that the temperature dependent hadronic radii shall diverge when hadron dissociation (Mott effect) sets in, driven basically by the restoration of chiral symmetry. As a consequence, in the vicinity of the chiral restoration temperature all meson radii shall behave like that of the pion and all baryon radii like that of the nucleon.
The resulting energy density and pressure behavior is shown in the right panel of Fig. 1. This part of the model we call Mott-Hagedorn-Resonance-Gas (MHRG). When all hadrons are gone at T ∼ 250 MeV, we are clearly missing degrees of freedom!

Quarks, gluons and hadron resonances
We improve the PNJL model over its standard versions [8,9] by adding perturbative corrections in O(α s ) for the high-momentum region above the threemomentum cutoff Λ. In the second step, the MHRG part is replaced by its final form, using the state-dependent spectral function for the description of the Mott dissociation of all hadron resonances above the chiral transition. The total pressure obtains the form P (T ) = P MHRG (T ) + P PNJL,MF (T ) + P 2 (T ) .
where P MHRG (T ) stands for the pressure of the MHRG model, accounting for the dissociation of hadrons in hot dense, matter.
The O(α s ) corrections can be split in quark and gluon contributions where P quark 2 stands for the quark contribution and P gluon 2 contains the ghost and gluon contributions. The total perturbative QCD correction to O(α s ) is where The corresponding contribution to the energy density is given in standard way by Eq. (5).
We will now include an effective description of the dissociation of hadrons due to the Mott effect into the hadron resonance gas model by including the state dependent hadron resonance width (6) into the definition of the HRG pressure From the pressure as a thermodynamic potential all relevant thermodynamical functions can be obtained. Combining the α s corrected meanfield PNJL model for the quark-gluon subsystem with the MHRG description of the hadronic resonances we obtain the results shown in the right panel of Fig. 1 where the resulting partial contributions in comparison with lattice QCD data from Ref. [2] are shown. We see that the lattice QCD thermodynamics is in full accordance with a hadron resonance gas up to a temperature of ∼ 170 MeV which corresponds to the pseudocritical temperature of the chiral phase transition. The lattice data saturate below the Stefan-Boltzmann limit of an ideal quark-gluon gas at high temperatures. The PNJL model, however, attains this limit by construction. The deviation is to good accuracy described by perturbative corrections to O(α s ) which vanish at low temperatures due to an infrared cutoff procedure. The transition region 170 ≤ T [MeV] ≤ 250 is described by the MHRG model, resulting in a decreasing HRG pressure which vanishes at T ∼ 250 MeV.
We have presented two stages of an effective model description of QCD thermodynamics at finite temperatures which properly accounts for the fact that in the QCD transition region it is dominated by a tower of hadronic resonances. To this end we have further developed a generalization of the Hagedorn resonance gas thermodynamics which includes the finite lifetime of hadronic resonances in a hot and dense medium by a model ansatz for a temperature-and mass dependent spectral function.