Bled Workshops in Physics Vol. 17, No. 1 p.l A Proceedings of the Mini-Workshop Quarks, Hadrons, Matter Bled, Slovenia, July 3 - 10, 2016 Limits on hadron spectrum from bulk medium properties* Wojciech Broniowskia'b a Institute of Physics, Jan Kochanowski University 25-406 Kielce, Poland b The H. Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, 31-342 Cracow, Poland Abstract. We bring up the fact that the bulk thermal properties of the hadron gas, as measured on the lattice, preclude a very fast rising of the number of resonance states in the QCD spectrum, as assumed by the Hagedorn hypothesis, unless a substantial repulsion between hadronic resonances is present. If the Hagedorn growth continued above masses ~ 1.8 GeV, then the thermodynamic functions would noticeably depart from the measured lattice values at temperatures above 140 MeV, just below the transition temperature to quark-gluon plasma. In this talk we point out the sensitivity of thermal bulk medium properties (energy density, entropy, sound velocity...) to the spectrum of the hadron resonance gas. In particular, we explore the effects of the high-lying part of the spectrum, above ~ 1.8 GeV, where it is poorly known, on the thermal properties still below the cross-over transition to the quark-gluon plasma phase. Such investigations were carried out in the past by various authors, see [1-5] and references therein, where the reader may find more details and results. The presently established QCD spectrum reaches about 2 GeV, and it is a priori not clear what happens above. Does the growth continue, or is saturated? As is evident from Fig. 1, the Hagedorn hypothesis [7] works very well up to about 1.8 GeV [8]. In the following, we explore two models: 1) hadron resonance gas with the Breit-Wigner width, HRG(r), which takes into account all states listed in the Particle Data Group tables [6] with mass below 1.8 GeV, and 2) this model amended with the states above 1.8 GeV, modeled with the Hagedorn formula fitted to the spectrum at lower masses (see Fig. 1). In short, model 1) includes the up-to-now established resonances, and model 2) extends them according to the Hagedorn hypothesis. First, we recall the fact that the inclusion of widths of resonances [11], as listed in the Particle Data Group tables, affects the results noticeably and in fact improves them. This is shown in Fig. 2, where the hadron resonance gas calculation for the QCD trace anomaly, e — 3p divided by T4. Here e stands for the energy density, p for the pressure, and T for the temperature. In the calculation, * Supported by the Polish National Science Centre grants DEC-2015/19/B/ST2/00937 and DEC-2012/06/A/ST2/00390. Wojciech Broniowski o 1000 "öd _C2 $ 100 B 10 Hagedorn HRG(r) M [GeV] Fig. 1. Solid line, labeled HRG(F): Number of QCD states (mesons, baryons, and antibaryons combined) with mass below M. All stable particles and resonances from the Particle Data Group tables [6] are included and their Breit-Wigner width is taken into account. Dashed line: the fit with the Hagedorn formula for the density of states, p(m) = A exp(m/TH), with TH = 260 MeV. iL CL CO 2 I 2 1 • WB • HotQCD 300 T [MeV] Fig. 2. The QCD trace anomaly (divided by T4) plotted as a function of temperature T. The inclusion of width of resonances to the hadron resonance gas improves the agreement with the lattice data from the Wuppertal-Budapest (WB) [9] and Hot QCD [10] collaborations. 2 4 0.5 1.0 1.5 2.0 2.5 4 3 0 the hadrons are treated as components of an ideal gas of fermions and bosons. We note that the overall agreement with of the hadron resonance gas model HRG(F) with the lattice measurement is remarkable. The virial expansion of Kamerlingh Onnes yields p/T = p + B2(T)p2 + B3(T)p3 +____Correspondingly, for the partition function of a thermodynamic system including the 1 —» 1, 2 —» 2, etc., processes one has ln Z = ln Z(1) + ln Z(2) + . (1) Limits on hadron spectrum from bulk medium properties 3 The non-interacting term ln Z(1) = ^ ln Zfble = ^ fkV ' d3p (2n)3 ln 1 ± e-Ep/T ±1 (2) includes the sum over all stable particles, whereas the second-order virial term involves the sum over pairs of stable particles denoted as K, ln Z(2) = ^ fxV d5K(M) ndM dM d3 P (2n)3 ln 1 ± e-Ep/T ±1 (3) where SK(M) stands for the phase shift in the channel K. For narrow resonances the correction to the density of two-body states d5K(M)/(ndM) [12] can be accurately approximated with the Breit-Wigner form, which is a basis of the hadron resonance gas model. In Fig. 3 we show the result of extending the Hagedorn hypothesis above the present experimental limit on the QCD spectrum. We note that the inclusion of extra (non-interacting) states above M = 1.8 GeV has a quite dramatic effect on the trace anomaly 0£ = e — 3p, placing it way above the lattice data at T > 140 MeV (the model calculation is credible below T ~ 170 MeV, where a crossover to the quark-gluon plasma occurs). A similar conclusion is drawn for other thermodynamic quantities, such as the entropy (cf. Fig. 4) or the sound velocity (cf. Fig. 5). T [MeV] Fig. 3. Same as in Fig. 2 but with the lines denoting the hadron resonance gas model, HRG(F), up to M = 1.8 GeV, and this model amended with the Hagedorn spectrum above M = 1.8 GeV. Therefore, if the hadron resonances were non-interacting, there would be no room for extra states above 1.8 GeV in the QCD spectrum. This conclusion may be affected by repulsion between the states (e.g., the excluded volume corrections), which decreases the contribution to the partition function. The issue is 2 Wojciech Broniowski T [MeV] Fig. 4. Same as in Fig. 3 but for the entropy density divided by T3. T [MeV] Fig. 5. Same as in Fig. 3 but for the square of the sound velocity. discussed quantitatively in [3,4], where a reduction of contributions to the thermodynamic quantities is assessed. The excluded volume reduces the contribution of resonances, and this makes them possible to appear in the spectrum in an "innocuous" way. The effect is explicit in Eq. (3), as repulsion leads to a decrease of the phase shift with M, or a negative correction to the density of states d6K(M)/(ndM). An important example of such an explicit cancellation occurs in the case of the a meson, whose contribution to one-body observables is canceled by the isospin-2 channel [12]. The case of the trace anomaly is shown in Fig. 6. Note that the phase shift taken into account in this analysis automatically includes the short-distance repulsion in specific channel, hence there is no need to model it separately. The cancellation experienced by the a state may occur also for other states with higher mass. Limits on hadron spectrum from bulk medium properties 3 "L 0.2 CL CO A 0.1 0.0 0.3 0.05 0.10 0.15 T [GeV] Fig. 6. Contributions to the trace anomaly from the pions, p mesons, a meson, and the isospin-2 component of the pion-pion interaction. We note an almost perfect cancellation of the a and isospin-2 channels. In conclusion, the thermodynamic quantities offered by the modern lattice QCD calculations allow to place limits on the high-lying spectrum on the QCD resonances, but the interactions between the states, such as the short-range repulsion, must be properly taken into account, as the two effects: increasing the number of states and introducing repulsion works in the opposite way. References 1. F. Karsch, K. Redlich, and A. Tawfik, Eur. Phys. J. C29, 549 (2003), arXiv:hep-ph/0303108 [hep-ph] . 2. P. Huovinen and P. Petreczky, Nucl. Phys. A837, 26 (2010), arXiv:0912.2541 [hep-ph] . 3. A. Andronic, P. Braun-Munzinger, J. Stachel, and M. Winn, Phys. Lett. B718, 80 (2012), arXiv:1201.0693 [nucl-th] . 4. E. Ruiz Arriola, L. L. Salcedo, and E. Megias, 54th Cracow School of Theoretical Physics: QCD meets experiment: Zakopane, Poland, June 12-20, 2014, Acta Phys. Polon. B45, 2407 (2014), arXiv:1410.3869 [hep-ph] . 5. P. M. Lo, M. Marczenko, K. Redlich, and C. Sasaki, Phys. Rev. C92, 055206 (2015), arXiv:1507.06398 [nucl-th] . 6. K. A. Olive, Chin. Phys. C40,100001 (2016). 7. R. Hagedorn, Nuovo Cim. Suppl. 3,147 (1965). 8. W. Broniowski, W. Florkowski, and L. Ya. Glozman, Phys. Rev. D70, 117503 (2004), arXiv:hep-ph/0407290 [hep-ph] . 9. S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, and K. K. Szabo, Phys. Lett. B730, 99 (2014), arXiv:1309.5258 [hep-lat] . 10. A. Bazavov et al. (HotQCD), Phys. Rev. D90,094503 (2014), arXiv:1407.6387 [hep-lat]. 11. E. Ruiz Arriola, W. Broniowski, and P. Masjuan, Proceedings, Conference on Modern approaches to nonperturbative gauge theories and their applications (Light Cone 2012): Cracow, Poland, July 8-13,2012, Acta Phys. Polon. Supp. 6,95 (2013), arXiv:1210.7153 [hep-ph]. 12. W. Broniowski, F. Giacosa, and V. Begun, Phys. Rev. C92, 034905 (2015), arXiv:1506.01260 [nucl-th] .