Thermal model description of p–Pb collisions at sNN=5.02\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}} = 5.02$$\end{document} TeV

The ALICE data on light flavor hadron production obtained in p–Pb collisions at sNN=5.02\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}} = 5.02$$\end{document} TeV are studied in the thermal model using the canonical approach with exact strangeness conservation. The chemical freeze-out temperature is independent of centrality except for the lowest multiplicity bin, with values close to 160 MeV but consistent with those obtained in Pb–Pb collisions at sNN=2.76\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}} = 2.76$$\end{document} TeV. The value of the strangeness non-equilibrium factor γs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma _s$$\end{document} is slowly increasing with multiplicity from 0.9 to 0.96, i.e. it is always very close to full chemical equilibrium.


Introduction
Over the last three decades, the hadron resonance gas (thermal model for short) in its grand-canonical (GC) and canonical formulations has been very successful in describing the abundances of light flavored hadron obtained in heavy-ion collisions [1][2][3] up to the highest beam energies.
The thermal model has a long history going back to papers by Heisenberg [4], Fermi [5] and Hagedorn [6]. The model uses a minimal number of parameters, mainly, the chemical freeze-out temperature T ch , the baryon chemical potential μ B and the volume V . The physical description emerging from this is a very consistent one, the freeze-out temperature rises rapidly at low beam energies to reach a plateau at around 160 MeV while the baryon-chemical potential decreases steadily from a very high value to become compatible with zero at the highest energies. The implication is that all hadrons are produced at one temperature and within one volume. It is therefore of interest to test the region of applicability of the In memory of Helmut Oeschler. a e-mail: natasha.sharma@cern.ch b e-mail: jean.cleymans@uct.ac.za c e-mail: kumar.lokesh@cern.ch thermal model and this is done in the present paper for p-Pb collisions at the LHC.
In the present paper we investigate in detail the particle abundances in p-Pb collisions [7][8][9][10] in the hope that this contributes to the understanding and the status of the model description.
At LHC energies, the particle abundances in central heavyion collisions have been well described [3,11,12] over nine orders of magnitude with only two parameters T ch and V , with μ B being restricted to zero because of the particleantiparticle symmetry at the LHC.
Our results show that there is no clear dependence on the centrality, as measured by the values of d N ch /dη, for the chemical freeze-out temperature. This confirms results obtained earlier on the dependence of thermal parameters on the size of the system obtained [13].
In Sect. 2 we briefly review the main features of the model. Section 3 presents our results for p-Pb collisions. Section 4 contains a discussion of results and Sect. 5 presents our conclusions.

The model
The identifying feature of the thermal model is that all the resonances listed in [14] are assumed to be in thermal and chemical equilibrium. This assumption drastically reduces the number of free parameters as this stage is determined by just a few thermodynamic variables namely, the chemical freeze-out temperature T ch , the various chemical potentials μ determined by the conserved quantum numbers and by the volume V of the system. It has been shown that this description is also the correct one [15][16][17] for a scaling expansion as first discussed by Bjorken [18]. After integration over p T these authors have shown that: where N 0 i is the particle yield as calculated in a fireball at rest. Hence, in the Bjorken model with longitudinal scaling and radial expansion the effects of hydrodynamic flow cancel out in ratios.
In general, if the number of particles carrying quantum numbers related to a conservation law is small, then the grand-canonical description no longer holds. In such a case the conservation of quantum numbers has to be implemented exactly in the canonical ensemble [19,20]. Here, we refer only to strangeness conservation and consider charge and baryon number conservation to be fulfilled on the average in the grand canonical ensemble because the number of charged particles and baryons is much larger than that of strange particles [21].
In the Grand-Canonical Ensemble (GCE), the volume V , temperature T ch and the chemical potentials μ determine the partition function Z (T, V, μ). In the hadronic fireball of noninteracting hadrons, ln Z is the sum of the contributions of all i-particle species where μ = (μ B , μ S , μ Q ) are the chemical potentials related to the conservation of baryon number, strangeness and electric charge, respectively. Z 1 i is the one-particle partition function.
The partition function contains all information needed to obtain the number density n i of particle species i. Introducing the particle's specific chemical potential μ i , one gets Any resonance that decays into species i contributes to the yields eventually measured. Therefore, the contributions from all heavier hadrons that decay to hadron i are included. In the GCE the particle yields are determined by the volume of the fireball, its temperature and the chemical potentials.
If the number of particles is small, conservation laws have to be implemented exactly. We refer here only to the Strangeness Canonical Ensemble (SCE) in which strangeness conservation is considered and charge and baryon number are conserved on the average. The density of strange particle i carrying strangeness s can be obtained from [21], where Z C S=0 is the canonical partition function where Z 1 i is the one-particle partition function calculated for μ S = 0 in the Boltzmann approximation. The arguments of the Bessel functions I s (x) and the parameters a i are introduced as, where S s is the sum of all Z 1 k (μ S = 0) for particle species k carrying strangeness s.
In the limit where x n < 1 (for n = 1, 2 and 3) the density of strange particles carrying strangeness s is well approximated by [21] From these equations it is clear that in the canonical ensemble the strange particle density depends explicitly on the volume through the arguments of the Bessel functions. It has been suggested that this volume might be different from the overall volume V [22,23]. In our analysis we will not entertain this possibility and work instead with an overall parameter γ s to describe the deviation of strange particles from chemical equilibrium [24]. This amounts to replacing each particle density by where |s| is the sum of the number of strange and anti-strange quarks. Note that the φ meson picks up a factor γ 2 s since it contains a strange and an anti-strange quark. This makes the φ meson behave in a similar way to the baryon which contains two strange quarks, a behavior which is supported by the data. The changes for the other strange mesons and baryons are as follows:

Scaling of yields with strange quark content
The experimental results were taken from [7][8][9], these are obtained for central rapidity in an interval y = 1. One of the basic features of the SC approach is that the particle ratios exhibit an approximate scaling with the difference  Fig. 1. The left pane shows the ratios of various strange particle to pion yields i as a function of multiplicity. Following references [7][8][9] we use K /π as an abbreviation for the ratio (K + + K − )/(π + + π − ) and similarly for all other ratios, except that φ/π is an abbreviation for φ/(π + + π − )/2. Scaling these ratios with the power of the corresponding number of strange quarks one observes a common trend (right pane of Fig. 1). As can be seen from Fig. 1 scaling with the number of strange quarks is very well supported by the experimental data. This scaling behavior with strange quark content is motivated by the SCE description including the factor γ s . In particular, the yield of φ mesons closely follows the one of e.g.
baryons. This argues against the use of a different correlation radius for the strangeness suppression as this would not affect the φ meson since it has strangeness zero.
We have used THERMUS [25] to perform these fits in the strangeness canonical ensemble (SCE) with μ B and μ Q fixed to 0 and using only one radius R for the system. The fits to the bins with the most central collisions are shown in Fig. 2. In the following figures π refers to (π + + π − )/2, K refers to K + + K − )/2 and similarly p refers to ( p +p)/2, refers to ( +¯ )/2, refers to ( − +¯ )/2 and finally refers to ( +¯ )/2. All fits were performed including the φ meson.
The upper part of the figures compares the experimental results, shown as a black circle, to the fits obtained from the thermal model as described in the previous section. As the There are no very outstanding deviations from the thermal model fits. In particular, the φ yield is reproduced reasonably well.
In Fig. 2 we present our results for the two most central multiplicity bins, i.e. 0-5 and 5-10%.
The pattern observed for the bins with the highest multiplicity reproduces itself for the next two bins i.e. 10-20 and 20-40% as shown in Fig. 3. The deviations from the thermal model are very similar to the previous ones. In particular, the yield of φ mesons is again reasonably well reproduced. This confirms the use of the strangeness non-equilibrium factor γ s . The use of a different correlation radius for the strangeness suppression would not affect the φ meson since it has net zero strangeness.
All yields are compatible with the thermal model description within two standard deviations.The next multiplicity bins are shown in Fig. 4. The left-hand side is consistent with the results shown in Figs. 2 and 3. The right-hand pane is a peripheral one corresponding to 60-80% centrality and shows the discrepancies with the thermal model fits. Especially there are clear deviations for the and baryons. The φ is reasonably well described.
The largest deviations appear for the most peripheral collisions, 80-100%, and are shown in Fig. 5. Especially the values for data/fit are very high for the baryon at 2.45 and for the it is 4.6. To summarize the results obtained above we show the chemical freeze-out temperature T ch and the corresponding radius R as a function of charged particle multiplicity in Fig. 6. The values of T ch are remarkably independent of multiplicity except for the most peripheral bin. The radius determined from the yields fixes the normalization of the yields and shows a steady increase with the centrality of the collision. The more particles there are in the final state, the larger the radius and, correspondingly, the volume of the system. It should be noted that the volume increases linearly with the multiplicity. This indicates that the density of the fireball remains constant.
The values of the strangeness non-equilibrium factor γ s are shown on Fig. 7. It can be seen that γ s is very slowly Fig. 7 Values of the strangeness non-equilibrium factor γ s as a function of the charged particle multiplicity d N ch /dη |η|<0.5 increasing with multiplicity, starting from a value of about 0.9 for the most peripheral collision reaching about 0.96 for the bin with the highest multiplicity. This shows that the system is very close to chemical equilibrium, the largest deviation being about 10%. Again, being so close to chemical equilibrium is a remarkable property of the fireball produced.
Apart from the most peripheral collisions, the chemical freeze-out temperature shows a remarkable consistency for all multiplicity classes. This shows that the final state hadrons are consistent with being produced in a single fireball having a temperature slightly about 160 MeV. The best fits are achieved for the most central collisions, the worst fits are for the two most peripheral bins. This is quantitatively reflected in the values of the χ 2 shown in Table 1.
The information presented above is summarized in Table 1I below. In addition the values of χ 2 are also listed. As can be seen by far the worst fits are obtained for the most peripheral collisions. Fig. 6 The temperature (left pane) and the radius (right pane) at chemical freeze-out as a function of charged particle multiplicity d N ch /dη |η|<0.5

Conclusions
The fireball being produced in p-Pb collisions has simple properties, it is close to full chemical equilibrium, the temperature is independent of the multiplicity in the final state and the radius increases smoothly which is to be expected for a system that increases with the multiplicity in the final state.
It is remarkable that the system is always close to full chemical equilibrium with γ s never below 0.9 even for the most peripheral bins. The chemical freeze-out temperature is also remarkably constant, except (again) for the most peripheral bin. The radius decreases when going to more peripheral collisions, a feature which is to be expected.
Fits to pp collisions have recently been made [26] as a function of the multiplicity of charged particles in the event. It was shown that for the highest multiplicities the thermal model works very well and that a grand canonical description becomes valid.
For such a simple model based on only three parameters it is not surprising that a perfect agreement is not always reached, this can be clearly seen in the underestimation of the pion yield and the overestimate of the proton one. This indicates that the thermal model as presented here is, in our opinion, a very good first order description of the particle composition but that it is open to further improvements. Some questions about this approach have been raised recently, e.g.: -what is the effect of an incomplete list of hadronic resonances [27][28][29][30]? -are there separate freeze-out temperature for strange particles [31,32]? -is the theoretical description complete [33]?
The first question has to be answered by new experimental results on hadronic resonances. The latter two cases naturally lead to the introduction of new extra parameters which reduce the simplicity of the model.
In conclusion, in this paper we have performed an analysis of the particle composition in p-Pb collisions at 5.02 TeV using the thermal model. This model provides a good description of the hadronic yields. The largest deviations occurring for the most peripheral collisions. The chemical freeze-out temperature is independent of centrality T ch = 162±3 MeV, except for the lowest multiplicity bin, this value is consistent with values obtained in Pb-Pb collisions. The value of the strangeness non-equilibrium factor γ s is slowly increasing with multiplicity from 0.9 to 0.96, i.e. it is always very close to full chemical equilibrium.