Contents lists available at ScienceDirect

Physics Letters B

www.elsevier.com/locate/physletb

Hadronization conditions in relativistic nuclear collisions and the QCD pseudo-critical line

Francesco Becattini3*, Jan Steinheimerb, Reinhard Stockc, Marcus Bleicherb

CrossMark

a Università di Firenze and ¡NFN Sezione di Firenze, Firenze, Italy b Frankfurt Institute for Advanced Studies (FIAS), Frankfurt, Germany c Institut für Kernphysik, Goethe-Universität and FIAS, Frankfurt, Germany

A R T I C L E I N F 0

Article history:

Received 13 June 2016

Received in revised form 15 November 2016

Accepted 17 November 2016

Available online 22 November 2016

Editor: J.-P. Blaizot

A B S T R A C T

We compare the reconstructed hadronization conditions in relativistic nuclear collisions in the nucleon-nucleon centre-of-mass energy range 4.7-2760 GeV in terms of temperature and baryon-chemical potential with lattice QCD calculations, by using hadronic multiplicities. We obtain hadronization temperatures and baryon chemical potentials with a fit to measured multiplicities by correcting for the effect of post-hadronization rescattering. The post-hadronization modification factors are calculated by means of a coupled hydrodynamical-transport model simulation under the same conditions of approximate isothermal and isochemical decoupling as assumed in the statistical hadronization model fits to the data. The fit quality is considerably better than without rescattering corrections, as already found in previous work. The curvature of the obtained "true" hadronization pseudo-critical line k is found to be 0.0048 ± 0.0026, in agreement with lattice QCD estimates; the pseudo-critical temperature at vanishing ib is found to be 164.3 ± 1.8 MeV.

© 2016 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license

(http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

1. Introduction

It is the goal of Quantum Chromo-Dynamics (QCD) thermodynamics to study the phase diagram of strongly interacting matter. Its most prominent feature, the transition line between hadrons and partons, in the plane spanned by the baryon-chemical potential i^b and the temperature T, is located in the non-perturbative sector of QCD. Here, the theory can be solved on the lattice and has recently led to calculations of the curvature of the parton-hadron boundary line [1-8]. This line can also be studied experimentally, in relativistic collisions of heavy nuclei, where apparently local thermodynamical equilibrium is achieved at a temperature well above the (pseudo-)critical QCD temperature Tc . Expansion and cooling then take the system down to the phase boundary where hadronization occurs. We have lately demonstrated [9-11] that post-hadronization inelastic rescattering, chiefly baryon-antibaryons annihilation, is an important feature of the process, which drives the system slightly out of equilibrium from the primordial hadronization equilibrium, implying an actual distinction between hadronization and chemical freeze-out. This

* Corresponding author.

E-mail address: becattini@fi.infn.it (F. Becattini).

rescattering stage is taken into account in state-of-the art simulations of the QGP expansion [12-14], where the local equilibrium particle distribution (through the so-called Cooper-Frye formula) at some critical values of T and ib is used to generate hadrons and resonances which subsequently undergo collisions and decay. By calculating the modification of the multiplicities brought about by the rescattering stage - the so-called afterburning - it is possible to reconstruct the hadronization point by means of a fit to the multiplicities in the framework of the Statistical Hadroniza-tion Model (SHM). Strictly speaking, this method allows to pin down the latest chemical equilibrium point [11] henceforth denoted as LCEP - i.e. the point when the primordial chemical equilibrium starts being distorted by the afterburning. As equilibrium is an intrinsic feature of hadronization [15-17] - as shown by the analysis of elementary collisions - most likely LCEP coincides with hadronization itself, as the maintaining of full chemical equilibrium in a rapidly expanding hadronic system, for the time needed to produce a measurable temperature shift, is highly unlikely.

As the primordial system temperature (baryon-chemical potential) shifts upward (downward) with increasing collision energy, an ascending sequence of experimental energies can, thus, map a sequence of LCEPs or hadronization points along the QCD transition line. This was the main point of ref. [10] where we showed that,

http://dx.doi.org/10.1016/j.physletb.2016.11.033

0370-2693/© 2016 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

indeed, the reconstructed LCEPs seem to follow the extrapolated lattice QCD pseudo-critical line in the (j^B , T) plane as determined in, e.g., ref. [3]. The agreement between lattice QCD calculations and the reconstructed hadronization points in relativistic heavy ion collisions seem to imply that, in the examined energy range (Vsnn > 7 GeV), the pseudo-critical line has indeed been crossed, and, thus, those energies lie above the so-called "onset of decon-finement" [18].

This conclusion is less straightforward than it may seem at a glance because, as has been mentioned, hadron formation is evidently a universal statistical process [15-17] in all kinds of collisions at the same hadronization temperature, with a difference in the strangeness sector, whose phase space appears to be only partially filled in elementary collisions [19-23].1 Indeed, even if the strangeness phase space was fully saturated in elementary collisions, if hadron production process in a nuclear collision was fully consistent with a picture of subsequent and independent elementary hadronic reactions, strange particle production would be strongly suppressed by the exact strangeness conservation over the typical small volumes of an elementary collision (canonical suppression) and subsequent hadronic inelastic collisions would not be able to raise multi-strange particle abundance to the measured one, as it is shown by transport calculations [24-26]. Hence, the observation of a strange particle production in agreement with the prediction of the SHM for a coherent, large volume, and the agreement between lattice QCD extrapolations of the pseudo-critical line and the reconstructed hadronization point is a strong evidence that the pseudo-critical line has been overcome. However, this must cease to happen at some sufficiently low centre-of-mass energy, implying the failure of at least one of the above conditions. Estimates remain uncertain at present, pointing to a region between 4 and 8 GeV.

In this paper, we reexamine the hadronization conditions in rel-ativistic heavy ion collisions over the energy range from low SPS (7.6 GeV) to LHC (2.76 TeV) by using hadronic multiplicities. We also include the highest AGS energy point at = 4-5 GeV to probe the aforementioned deconfinement conditions further down in energy. For this purpose, we take advantage of an improved initialization of the afterburning process by enforcing a particle generation stage - or hydrodynamical decoupling - in UrQMD [27-30] at fixed values of energy density corresponding to mean temperatures and chemical potentials equals to those determined in ref. [10] and show how this leads to a further remarkable improvement of the fit quality compared to the plain statistical model fits [10,11,31]. Finally, we compare the resulting curvature of the LCEP-hadronization curve in the (/¿B, T) plane with the predictions of lattice QCD, reporting a good agreement.

2. Afterburning and modification factors

As has been mentioned, we studied the effect of post-hadroni-zation rescattering on hadron multiplicities and on the associated SHM fits in previous publications [9-11] employing a hybrid model [34,35] with a hydrodynamic expansion of the QCD plasma terminated at a predefined point where local equilibrium particle generation is assumed (hadronization), followed by a hadronic rescat-tering stage modelled by UrQMD [29] (afterburning). For the fluid dynamical simulation we employed an equation of state which follows from a so-called "combined hadron-quark model". It is based

Table 1

Energy densities used to implement hydrodynamic decoupling or Cooper-Frye parti-clization, at the different collision energies. Also quoted are the corresponding mean temperatures and baryon-chemical potentials.

VSnn (GeV) Energy density (MeV/fm3) Tcf (MeV) I^bcf (MeV)

4.75 508 135 563

7.6 435 155 426

8.7 435 161 376

17.3 435 163 250

2760 362 165 0

on a chiral hadronic model which provides a satisfactory description of nuclear matter properties. The quark phase is introduced as a PNJL type model. The transition from the hadronic to the quark phase occurs at about Tc & 165 MeV for fiB = 0, as a smooth crossover as shown in [32].

For each hadronic species a so-called modification factor is extracted which is defined as the ratio between the final multiplicity after the actual chemical (and kinetic) freeze-out (which is now species-dependent) and its value without afterburning (at hadronization):

f • = —

J J im

1 It is worth pointing out here that the strangeness undersaturation is still ob-

served in nuclear collisions at high energy but it can be accounted for by residual

nucleon-nucleon collisions nearby the outer edge of the nuclear overlapping region, see [15] and references therein.

The modification factors are then used as additional multiplicative factors to the theoretical equilibrium multiplicity yields in the SHM fit, ready to be compared to the data. Note that in the calculation of the modification factors, all weak decays are turned off, but all strong and EM decays are turned on; this limits the data analysis to measurements of multiplicities corrected for the weak decay feed-down.

In our previous studies, the hydrodynamic decoupling procedure was inspired by the so-called "inside-outside cascade" mechanism: the transition from the fluid dynamical phase to the hadronic transport part occurs in successive transverse slices, of thickness 0.2 fm, whenever all fluid cells of that slice fall below a critical energy density, that is six times the nuclear ground state density e & 850 MeV/fm3. In fact, in the present investigation, we have implemented an approximate isothermal termination of the hydrodynamical stage at some pre-established temperature TCF (the subscript CF stands for Cooper-Frye). This is certainly in much better accordance with the underlying picture of a statistical hadronization as well as with the previously discussed concept of LCEP, which is determined at a fixed value of the proper temperature. For the decoupling, the hypersurface is defined by a fixed energy density - at LHC energy - of approximately 0.360 GeV/fm3 which corresponds to a mean hadronization temperature close to 165 MeV at zero baryon-chemical potential (for the lower energies, see Table 1). The cell-to-cell temperature and chemical potential fluctuations corresponding to such a hydrodynamic decoupling procedure, at some given collision energy, are small. For instance, the dispersion of the temperature at LHC energy is ~1.5 MeV, the dispersion of the chemical potential at the SPS energy is of the order of 10 MeV; these values are comparable or smaller than the parameter fit errors (see Table 5).

The UrQMD model employs the hypersurface finder outlined in ref. [30], which is used in the Cooper-Frye prescription and sampled to produce hadrons in accordance with global conservation of charge strangeness baryon number and the total energy.

It should be pointed out that the calculated modification factors do depend on the chosen temperature TCF ending the hydro-dynamical expansion [33]. Ideally, this should coincide with the actual TLCEP at each energy, which is a priori unknown, except for a reasonable lower bound set by the chemical freeze-out temperature as determined in the traditional, plain, SHM fit. One may then

Table 2

Afterburning modification factors determined by means of UrQMD (see text for explanation). The input decoupling temperatures and chemical potentials are reported in Table 5.

Particle VSNN = 4-75 GeV VSNN = 7-6 GeV ySNN = 8-7 GeV VSNN = 17-3 GeV ySNN = 2760 GeV

n+ 1.01 0.998 1.01 1.03 1.05

n— 0.980 0.980 0.997 1.03 1.05

K+ 0.947 0.927 0.918 0.928 0.918

K— 0.890 0.901 0.868 0.891 0.919

p 0.975 0.978 0.979 0.956 0.758

p 0.341 0.369 0.340 0.533 0.754

A 0.981 0.935 0.932 0.927 0.816

A 0.432 0.478 0.422 0.601 0.821

3 — 1.01 0.979 0.977 0.970 0.886

3+ 0.573 0.575 0.531 0.698 0.876

a 0.888 0.882 0.808 0.873 0.789

a 0.596 0.573 0.566 0.706 0.778

0 - - - - 0.75

wonder to what extent T LCEP, which is the outcome of the subsequent SHM fit - corrected for afterburning - is affected by the chosen TCF. In general, larger TCF involve larger deviations of the modification factors from unity, so one could expect that the fit is influenced to such an extent that the corrected SHM fit tends to reproduce the initially chosen value TCF, making the whole method non-predictive. However, this would be the case only if the final particle multiplicities after freeze-out were independent of TCF. It was checked that in hybrid simulations this does not happen and final particle yields do depend on the chosen decoupling condition. In fact it appears that the change in the extracted TLCEP tends to be smaller than in the input TCF (i.e. independent of the exact values of the modification), and TLCEP shows a trend toward a definite value. For instance, for Pb-Pb collisions at VSNN = 8-7 GeV (see next section for details), for an input TCF = 144 MeV we obtained TLCEP = 155 MeV whereas for TCF = 161 MeV we obtained TLCEP = 163 MeV. In summary, the method converges.

The optimal situation, as has been mentioned, is TCF = T LCEP, which could be achieved by an iterative procedure; however, it would be computationally expensive and not worth the effort when - in view of the above observation - the difference between TCF and TLCEP is only few MeV's. Altogether, the small differences between TCF and TLCEP in our analysis (see Table 5) make us confident that the fitted thermal parameters are fully significant, with only a marginal dependence on the difference (TCF — TLCEP).

3. Data analysis

In this section we present the results of our analysis including 5 centre-of-mass energy points: VSNN = 2-76 TeV at the LHC, ^SNN = 17-3, 8-7, 7-6 GeV at the SPS and VSNN = 4-75 GeV at the AGS.

The modification factors for the strongly stable hadrons have been calculated according to the method described in the previous section and are quoted in Table 2. The decoupling conditions in terms of temperature, that is TCF, and baryon-chemical potential are those determined as LCEP's in our previous papers for the most central collisions at the LHC [11] and for SPS points [10] (see Table 5 further below). For the AGS point we did not have any clue about the TLCEP value, so we iterated the procedure until the fitted TLCEP was reasonably close to the TCF. The modification factor for the 0 is more difficult to extract than for other, strongly stable, particles as it is not the 0 itself which is absorbed, but its decay products which rescatter, making 0 reconstruction hardly feasible. Assuming that any rescattering of a decay product will lead to a loss of a 0, the lower bound of the survival probability has been estimated at 2.76 TeV to be about 0.75 [41-43]. At all lower energies, for afterburning is expected to be less important for this

0 100 200 300 400 500 600 700 Hb (MeV)

Fig. 1. (Color online.) x2/dof of the SHM multiplicity fits at the 5 different energies. Blue dots: plain SHM fit without afterburning corrections. Green dots: SHM fit with afterburning correction with the isochronous decoupling method (points from ref. [10]). Red dots: SHM fits with afterburning corrections with the new, approximately isothermal decoupling.

meson, we have used an educated guess of 0.875 which is the mean value between 0.75 and 1, varying between these bounds in order to check the stability of the best fit solutions.

The particle set used in the analysis is the intersection between the available measured multiplicities and the set of particles for which a modification factor was calculated, see Table 3. All data refer to central collisions of Au+Au (at the AGS) and Pb+Pb (at SPS and LHC). As has been mentioned, we have confined ourselves to data sets where weak feed-down was subtracted in order to make a proper comparison between corrected (with modification factors) and non-corrected fits.

The SHM, the formulas for primary and final multiplicities, the fitting procedure with and without modification factors have been described in detail elsewhere [9]. Herein, we simply summarize the obtained results in Table 5. The corrected fits are of remarkable better quality with respect to the plain SHM fits, as it is shown in Fig. 1, confirming previous findings. One exception stands out, the SPS point at 8.7 GeV, which is the point where the ratio K+ /n+ attains its maximum observed value [51]. Indeed, the measured ratio K+/n + overshoots the statistical model prediction [52] by more than 2a (see Table 4), a discrepancy which is not cured by the afterburning correction.

Notably, the x2/dof is of full statistical significance once the modification factors are introduced in the two highest energy points. Moreover, there is a further improvement of the fit quality with respect to the previous, non-isothermal, fits [10,11].

Table 3

Particle multiplicities measured at various centre-of-mass energies employed in our analysis. The data are 4n multiplicities except at Vsnn = 2760 GeV where they are midrapidity densities.

Particle VsNN = 4.75 GeV VsNN = 7.6 GeV v^SNN = 8.7 GeV •v/sNN = 17.3 GeV ySNN = 2760 GeV

n+ 133.7 ± 9.9 [36] 241 ± 12 [44] 293 ± 15.3 [44] 619 ± 35.4 [44] 733 ± 54 [45]

n - - 274 ± 14 [44] 322 ± 16.3 [44] 639 ± 35.4 [44] 732 ± 52 [45]

K+ 23.7 ± 2.86 [37] 52.9 ± 3.6 [44] 59.1 ± 3.55 [44] 103 ± 7.1 [44] 109 ± 9 [45]

K- 3.76 ± 0.47 [37] 16 ± 0.45 [44] 19.2 ± 1.12 [44] 51.9 ± 3.55 [44] 109 ± 9 [45]

p 1.23 ± 0.13a [38] - - - 34 ± 3 [45]

p - 0.26 ± 0.04b - 4.23 ± 0.35c 33 ± 3 [45]

A 18.1 ± 1.9 [39] 36.9 ± 3.3 [44] 43.1 ± 4.32 [44] 48.5 ± 8.6 [44] 26.1 ± 2.8 [46]

A 0.017 ± 0.005 [40] 0.39 ± 0.045 [44] 0.68 ± 0.076 [44] 3.32 ± 0.34 [44] -

a- - 2.42 ± 0.345 [44] 2.96 ± 0.41[44] 4.40 ± 0.64 [44] 3.57 ± 0.27c [47]

a+ - 0.120 ± 0.036 [44] 0.13 ± 0.022 [44] 0.71 ± 0.1 [44] 3.47 ± 0.26c [47]

Q - - 0.14 ± 0.05b [44] 0.59 ± 0.11 [44] 1.26 ± 0.22d,e [47]

Q - - - 0.260 ± 0.067 [44] -

<P - 1.84 ± 0.36 [44] 2.55 ± 0.25 [44] 8.46 ± 0.50 [44] 13.8 ± 1.77 [48]

Bf 363 ± 10 [37] 349 ± 5.1 [44] 349 ± 5.1 [44] 362 ± 8 [44] -

a This is the ratio p/n+.

b Our extrapolation [9] based on measurements in ref. [49].

c Our extrapolation [9] of spectra measured in ref. [50]. The NA49 data compilation [44] quotes 4.25 ± 0.28 by M. Utvic.

d Q + Q.

e Interpolation to 0-5% centrality quoted in ref. [11]. f Number of participants = net baryon number of the fireball.

Table 4

Measured vs fitted K+ multiplicities at VsNN = 7.6 and 8.7 GeV, in the so-called horn region, with corresponding deviations in units of a within round brackets.

v^nn (GeV) Measured Plain fit Modified fit

7.6 52.9 ± 3.6 47.1 (-1.6) 45.3 (-2.1)

8.7 59.1 ± 3.6 51.4 (-2.2) 49.6 (-2.7)

Table 5

Results of the SHM fits with and without afterburning corrections. In the last col-

umn we quote the corresponding decoupling temperatures TCF and chemical poten-

tials 1B cf employed to calculate the modification factors.

Parameters Without afterburner With afterburner CF in URQMD

Au-Au VSNN : = 4.75 GeV

T (MeV) 122.1 ± 4.0 130.5 ± 12.3 135

1b (MeV) 563 ± 15 588 ± 32 563

Ys 0.638 ± 0.074 0.71 ± 0.12 1.0

X2/dof 4.5/3 4.9/3

Pb-Pb ysNN = 7.6 GeV

T (MeV) 139.6 ± 3.7 157.7 ± 4.3 155

1b (MeV) 437 ± 20 424 ± 11 426

Ys 0.922 ± 0.075 0.871 ± 0.059 1.0

X2/dof 22.6/7 12.8/7

Pb-Pb ysNN = 8.7 GeV

T (MeV) 148.2 ± 3.8 163.3 ± 5.0 161

1b (MeV) 385 ± 11 371 ± 12 376

Ys 0.783 ± 0.062 0.773 ± 0.055 1.0

X2/dof 17.6/7 20.2/7

Pb-Pb ysNN = = 17.3 GeV

T (MeV) 150.4 ± 3.9 162.3 ± 2.7 163

1b (MeV) 265 ± 10 244 ± 6 250

Ys 0.914 ± 0.052 0.885 ± 0.029 1.0

X2/dof 26.9/9 9.1/9

Pb-Pb vsNN = = 2760 GeV

T (MeV) 155.0 ± 3.7 163.8 ± 3.3 165

1b (MeV) 0 (fixed) 0 (fixed) 0

Ys 1.07 ± 0.05 1.02 ± 0.04 1.0

X2/dof 15.2/8 4.7/8

The quoted errors in Table 5 are the fit errors. There are additional small systematic uncertainties on the fit parameters related to the errors on the modification factors. These errors stem from the uncertainties on the cross-sections used in UrQMD and from

finite Monte Carlo statistics. The former are difficult to estimate, whereas the latter are simpler; in our runs they are of the order of few percent for all particles, thus they do not imply any significant variation of the best fit parameter values. The only largely unknown modification factor is, as has been mentioned, the 0's, for which we could obtain an estimate of 0.75 at the top energy point = 2760 GeV. As the rescattering of a neutral meson

is expected to diminish in a lower multiplicity environment, one can reasonably set a lower bound of 0.75 at all lower energies. To estimate the effect of the uncertainty, we have varied the 0 modification factor to 0.75 and 1 in turn at each energy point. The resulting variation of the fit parameters is within 1 MeV for the temperature and few MeV's for the baryon chemical potentials so that the relative systematic error is always less than 1%, thus below the fit error.

4. Curvature of the pseudo-critical line

Finally, we have used the LCEP points in Table 5 to fit the curvature of the pseudo-critical line in the (jB, T) plane. The curvature parameter k is defined by the equation:

Tc (1b ) = Tc (0)

1b Tc (0)

which is the same formula used in lattice calculations. As the LCEP points in (jB, T) plane have errors on both coordinates, we have minimized the x2:

X2 = (Zi - Z0i)TC-\Zi - Z0i)

Zi = (jBi, Ti) Z0i = (j°Bi, Tc(jBi))

In the above equation, jBi and Ti are the output of the SHM fit of the i-th energy point, while Ci is their corresponding covariance matrix. The j0Bi's are free parameters which represent the "true" values of the chemical potential in the fitted curve, Tc(j0Bi) being the corresponding "true" temperatures. Therefore, the free parameters in this fit are the chemical potentials j0Bi, whose position is

Table 6

Best fit parameters and x2 values for the fit to the reconstructed LCEPs and chemical freeze-out points in the (/¿B, T) plane.

Fit method Tc(0) (MeV) K X x2

4 points 164.3 ± 1.8 0.0048 ± 0.0026 - 0.47/2

4 points, freezeout 157.4 ± 6.2 0.013 ± 0.0072 - 1.0/2

5 points 167.7 ± 4.0 0.0111 ± 0.0055 - 2.8/3

5 points, freezeout 162.1 ± 4.4 0.020 ± 0.004 - 2./3

Quartic 5 points 164.4 ± 2.7 0 ± 0.0091 0.0109 ± 0.00047 0.97/2

Table 7

Comparison between the curvature k in lattice QCD calculations and our estimate.

0 100 200 300 400 500 600 700

Mb (MeV)

Fig. 2. (Color online.) Reconstructed LCEPs (red squared dots) vs plain chemical freeze-out fitted points (blue round dots) in the (/¿B, T) plane. The solid lines are the 4 points quadratic fits quoted in Table 6; the dashed line is the 5 point quadratic fit including the lowest energy AGS point.

strongly constrained by the "measured" values ¡iBi, the value of the pseudo-critical temperature Tc (0) and k .

It should be pointed out that the equation (2) is a quadratic approximation of the actual pseudo-critical line, hence deviations are expected at large values of the chemical potentials. Therefore, we have first excluded the lowest energy point and made a fit to the four highest energy points at our disposal. We have also compared with the freezeout points, for which many systematic studies have been done in the past [53]. The fitted values of Tc(0) and k are reported in Table 6 while the fitted curves are shown in Fig. 2. It can be seen that the fit quality is excellent for the LCEP points and satisfactory for the plain freeze-out points. The systematic error on the curvature due to the uncertainty in the 0 meson modification factors has been obtained repeating the fit with the varied (j^B, T) points and turned out to be 0-0006.

The lowest energy point falls below the fit curve in both cases. There are three possible explanations for this:

1. the (mundane) effect of having excluded it from the fit;

2. the quadratic approximation in (2) falls short at such large chemical potential values;

3. the lowest energy point did not reach the pseudo-critical transition line, and so the onset of deconfinement can be located between 4.5 and 7.6 GeV.

The latter hypothesis is indeed the most intriguing, but its very consideration requires the ruling out the first two. If the AGS point is included, the fit quality is not as good as the 4 point fit, still it is within statistical significance, as it can be seen in Table 6. Yet, there is some tension between the fitted curve and the two extremal points (LHC and AGS) which both undershoot the curve by 4 and 15 MeV respectively, as it can be seen in Fig. 2. On the other hand, including a quartic term X(j^B/Tc(0))4 improves the fit but the limited number of points and the limited range does not

Reference K

This work - 4 points 0.0048 ± 0.0026

This work - 5 points 0.0111 ± 0.0055

[1] 0.0066 ± 0.0005

[2] (0.0033-0.0123)

[4] - Max 0.020 ± 0.002

[4] - Min 0.0066 ± 0.00020

[5] 0.0149 ± 0.0021

[7] 0.0135 ± 0.0020

[8] 0.020 ± 0.004

allow to pin down both the quadratic and the quartic term at the same time; indeed, the fit has multiple solutions and the best fit is awkwardly found for k ~ 0 (see Table 6).

Finally, returning to Fig. 2, which is our main result showing the estimated QCD pseudo-critical curve in the (jB, T) plane, it is appropriate to compare it with recent lattice QCD calculations (see Table 7). Because of the pseudo-critical nature of the transition, both Tc (0) and k parameters depend on the observable used to define it [4]. It could be therefore expected that these parameters will somewhat differ from those extracted with the fluctuation of conserved charges [54], which can be directly calculated in lattice QCD but are definitely less robust observables in relativistic heavy ion collisions with respect to mean multiplicities [55].

In our comparison, we have quoted all recent literature on the subject. We find that our main value of 0.0048 is in slightly better agreement with lower estimates [1,2,4]. We also note that our value is compatible with a recent estimate based on a comparison between lattice QCD and data [56].

5. Conclusions

In summary, we have determined the hadronization conditions (strictly speaking the latest chemical equilibrium points) in rela-tivistic heavy ion collisions, by using an improved calculation of the post-hadronization rescattering correction. The quality of the statistical model fit considerably improves with respect to traditional fits without afterburning corrections as well as with respect to our previous calculations. The pseudo-critical temperature at zero jB, determined with hadronic multiplicities, turns out to be Tc = 164 MeV, which is significantly higher than lattice QCD calculations based on different observables. We find good agreement between the extracted curvature of the hadronization curve and the corresponding QCD lattice calculations for the pseudo-critical line, with a preference for the lower estimates.

At this stage, it is not possible to make a definite statement about the crossing of the pseudo-critical line at the lowest energy point at V^nn = 4.7 GeV. This is due to lack of appropriate data and this issue will be tackled by future experiment in that energy range (NA61 at SPS and the facilities NICA and FAIR). Our observations might have interesting implications for the location of the critical point [57]; we note that a recent analysis [58] locate it at T ~ 165 MeV and jB ~ 95 MeV which just sits on our hadroniza-tion curve.

Acknowledgements

We are grateful to C. Bonati and V. Koch for useful suggestions. This work was supported by GSI and BMBF and by the University of Florence grant Fisica dei plasmi relativistici: teoria e appli-cazioni moderne. The computational resources were provided by the LOEWE Frankfurt Center for Scientific Computing (LOEWE-CSC).

References

[1] O. Kaczmarek, et al., Phys. Rev. D 83 (2011) 014504.

[2] P. Hegde, et al., Bielefeld-BNL-CCNU Collaboration, arXiv:1511.03378 [hep-lat].

[3] G. Endrodi, Z. Fodor, S.D. Katz, K.K. Szabo, J. High Energy Phys. 1104 (2011) 001.

[4] S. Borsanyi, G. Endrodi, Z. Fodor, S.D. Katz, S. Krieg, C. Ratti, K.K. Szabo, J. High Energy Phys. 1208 (2012) 053.

[5] R. Bellwied, S. Borsanyi, Z. Fodor, J. Günther, S.D. Katz, C. Ratti, K.K. Szabo, Phys. Lett. B 751 (2015) 559.

[6] C. Bonati, M. D'Elia, M. Mariti, M. Mesiti, F. Negro, F. Sanfilippo, Phys. Rev. D 90 (11) (2014) 114025.

[7] C. Bonati, M. D'Elia, M. Mariti, M. Mesiti, F. Negro, F. Sanfilippo, Phys. Rev. D 92 (5) (2015) 054503.

[8] P. Cea, L. Cosmai, A. Papa, Phys. Rev. D 93 (2016) 014507.

[9] F. Becattini, M. Bleicher, T. Kollegger, M. Mitrovski, T. Schuster, R. Stock, Phys. Rev. C 85 (2012) 044921.

[10] F. Becattini, M. Bleicher, T. Kollegger, T. Schuster, J. Steinheimer, R. Stock, Phys. Rev. Lett. 111 (2013) 082302.

[11] F. Becattini, E. Grossi, M. Bleicher, J. Steinheimer, R. Stock, Phys. Rev. C 90 (5) (2014) 054907.

[12] X. Zhu, F. Meng, H. Song, Y.X. Liu, Phys. Rev. C 91 (3) (2015) 034904.

[13] J. Auvinen, H. Petersen, Phys. Rev. C 88 (6) (2013) 064908.

[14] I. Karpenko, M. Bleicher, P. Huovinen, H. Petersen, arXiv:1601.00800.

[15] F. Becattini, An introduction to the statistical hadronization model, arXiv: 0901.3643;

F. Becattini, R. Fries, Landolt-Bornstein 23 (2010) 208, arXiv:0907.1031.

[16] R. Stock, Phys. Lett. B 456 (1999) 277.

[17] H. Satz, Eur. Phys. J. ST 155 (2008) 167.

[18] M. Gazdzicki, M. Gorenstein, P. Seyboth, Acta Phys. Pol. B 42 (2011) 307.

[19] F. Becattini, Z. Phys. C 69 (3) (1996) 485;

F. Becattini, A. Giovannini, S. Lupia, Z. Phys. C 72 (1996) 491.

[20] F. Becattini, P. Castorina, J. Manninen, H. Satz, Eur. Phys. J. C 56 (2008) 493.

[21] L. Ferroni, F. Becattini, Eur. Phys. J. C 71 (2011) 1824.

[22] F. Becattini, G. Passaleva, Eur. Phys. J. C 23 (2002) 551.

[23] P. Castorina, H. Satz, Adv. High Energy Phys. 2014 (2014) 376982.

[24] H. Petersen, M. Mitrovski, T. Schuster, M. Bleicher, Phys. Rev. C 80 (2009) 054910.

[25] M. Bleicher, F.M. Liu, A. Keranen, J. Aichelin, S.A. Bass, F. Becattini, K. Redlich, K. Werner, Phys. Rev. Lett. 88 (2002) 202501.

[26] H.J. Drescher, J. Aichelin, K. Werner, Phys. Rev. D 65 (2002) 057501.

[27] S.A. Bass, et al., Prog. Part. Nucl. Phys. 41 (1998) 255.

[28] M. Bleicher, et al., J. Phys. G 25 (1999) 1859.

[29] H. Petersen, J. Steinheimer, G. Burau, M. Bleicher, H. Stöcker, Phys. Rev. C 78 (2008) 044901.

[30] P. Huovinen, H. Petersen, Eur. Phys. J. A 48 (2012) 171.

[31] J. Stachel, A. Andronic, P. Braun-Munzinger, K. Redlich, J. Phys. Conf. Ser. 509 (2014) 012019.

[32] J. Steinheimer, S. Schramm, H. Stöcker, Phys. Rev. C 84 (2011) 045208.

[33] J. Steinheimer, J. Aichelin, M. Bleicher, Phys. Rev. Lett. 110 (2013) 042501.

[34] S.A. Bass, A. Dumitru, Phys. Rev. C 61 (2000) 064909.

[35] D. Teaney, J. Lauret, E.V. Shuryak, Nucl. Phys. A 698 (2002) 479.

[36] L. Ahle, et al., E-802 Collaboration, Phys. Rev. C 59 (1991) 2173.

[37] L. Ahle, et al., E-802 Collaboration, Phys. Rev. C 60 (1999) 044904.

[38] L. Ahle, et al., E802 Collaboration, Phys. Rev. C 60 (1999) 064901.

[39] S. Ahmad, et al., Phys. Lett. B 382 (1996) 35.

[40] S. Albergo, et al., Phys. Rev. Lett. 88 (2002) 062301.

[41] J. Steinheimer, J. Aichelin, M. Bleicher, EPJ Web Conf. 36 (2012) 00002.

[42] J. Steinheimer, M. Bleicher, EPJ Web Conf. 97 (2015) 00026.

[43] A.G. Knospe, C. Markert, K. Werner, J. Steinheimer, M. Bleicher, Phys. Rev. C 93 (1) (2016) 014911.

[44] NA49 collaboration compilation of numerical results https://edms.cern.ch/ ui/file/1075059/4/na49_compil_20130801.pdf and references therein.

[45] B. Abelev, et al., ALICE Collaboration, Phys. Rev. C 88 (2013) 044910.

[46] B.B. Abelev, et al., ALICE Collaboration, Phys. Rev. Lett. 111 (2013) 222301.

[47] B.B. Abelev, et al., ALICE Collaboration, Phys. Lett. B 728 (2014) 216;

B.B. Abelev, et al., Phys. Lett. B 734 (2014) 409, Erratum.

[48] B.B. Abelev, et al., ALICE Collaboration, Phys. Rev. C 91 (2015) 024609.

[49] C. Alt, et al., NA49 Collaboration, Phys. Rev. C 73 (2006) 044910.

[50] T. Anticic, et al., NA49 Collaboration, Phys. Rev. C 83 (2011) 014901.

[51] C. Alt, et al., NA49 Collaboration, Phys. Rev. C 77 (2008) 024903.

[52] F. Becattini, M. Gazdzicki, A. Keranen, J. Manninen, R. Stock, Phys. Rev. C 69 (2004) 024905.

[53] J. Cleymans, H. Oeschler, K. Redlich, S. Wheaton, Phys. Rev. C 73 (2006) 034905.

[54] F. Karsch, Cent. Eur. J. Phys. 10 (2012) 1234;

A. Bazavov, et al., Phys. Rev. Lett. 109 (2012) 192302;

P. Alba, W. Alberico, R. Bellwied, M. Bluhm, V. Mantovani Sarti, M. Nahrgang,

C. Ratti, Phys. Lett. B 738 (2014) 305;

S. Borsanyi, Z. Fodor, S.D. Katz, S. Krieg, C. Ratti, K.K. Szabo, Phys. Rev. Lett. 113 (2014) 052301.

[55] P. Braun-Munzinger, V. Koch, T. Schäfer, J. Stachel, Phys. Rep. 621 (2016) 76.

[56] A. Bazavov, et al., Phys. Rev. D 93 (2016) 014512.

[57] M.A. Stephanov, Prog. Theor. Phys. Suppl. 153 (2004) 139, Int. J. Mod. Phys. A 20 (2005) 4387.

[58] R.A. Lacey, Phys. Rev. Lett. 114 (14) (2015) 142301.