Abstract

Experimental results of the rapidity distributions of negatively charged pions produced in proton-proton (-) and beryllium-beryllium (Be-Be) collisions at different beam momentums, measured by the NA61/SHINE Collaboration at the super proton synchrotron (SPS), are described by a revised (three-source) Landau hydrodynamic model. The squared speed-of-sound parameter is then extracted from the width of rapidity distribution. There is a local minimum (knee point) which indicates a softest point in the equation of state (EoS) appearing at about  GeV/ (or 8.8 GeV) in excitation function (the dependence of on incident beam momentum (or center-of-mass energy)). This knee point should be related to the searching for the onset of quark deconfinement and the critical point of quark-gluon plasma (QGP) phase transition.

1. Introduction

Comparing with the relativistic heavy ion collider (RHIC) in USA [13] and the large hadron collider (LHC) in Switzerland [4, 5], the fixed target experiments performed at the super proton synchrotron (SPS) in Switzerland [6, 7] present relatively simple and clean collisions process. The multiplicity in collisions at the SPS is also low comparing with those at the RHIC and LHC. As one of the “first day” measurement quantities, the rapidity (pseudorapidity) distributions of charged particles are reported by experimental collaborations [17]. These distributions give us a chance to analyze longitudinal picture of particle productions. At the same time, based on the rapidity distributions, one can extract other information such as the penetrating (stopping) power of projectile (target) nucleus, energy and rapidity losses of projectile nucleus, and energy and particle densities of interacting region.

As a measurement of particle density and mean free path, the squared speed-of-sound parameter which characterizes partly the formation of matters in interacting region can be extracted from the width of Gaussian rapidity distribution described by the Landau hydrodynamic model [818]. Generally, the rapidity spectra obtained in collisions at the SPS, RHIC, and LHC are not simply Gaussian distributions. Instead, two-Gaussian distributions, one in the backward (target) region and the other in the forward (projectile) region, are needed to fit the experimental data [6, 19, 20]. However, the two-Gaussian distribution implies that there is no central source at midrapidity, which seems unbelievable. According to the three-source relativistic diffusion model [2124], a central source arising from interactions between low-momentum gluons in both target and projectile should be located at midrapidity, a target-like (projectile-like) source arising from interactions between valence quarks in the target (projectile) and low-momentum gluons in the projectile (target) are expected to appear in the backward (forward) rapidity region.

To extract the squared speed-of-sound parameter so that the formation of matters can be characterized by an alternative method, we need the Landau hydrodynamic model [810] and its improved or simplified version [1118] in the descriptions of rapidity distributions. In fact, the Landau hydrodynamic model is not good enough to fit the whole rapidity range. In our recent work [25, 26], we have revised the Landau hydrodynamic model [818] to a three-source situation. The rapidity distributions measured in experiments are then described by three Gaussian distributions. This picture is consistent with the three-source relativistic diffusion model [2124].

In this paper, we use the revised (three-source) Landau hydrodynamic model [25, 26] to describe the rapidity distributions of negatively charged pions produced in proton-proton (-) and beryllium-beryllium (Be-Be, exactly 7Be-9Be) collisions at different beam momentums measured by the NA61/SHINE Collaboration at the SPS. Combining with our previous works [25, 26], we observe the dependence of squared speed-of-sound on collision energy in a wide range.

The rest of this paper is structured as follows. The model is shortly described in Section 2. Results and discussion are given in Section 3. In Section 4, we summarize our main observations and conclusions.

2. The Model

The revised (three-source) Landau hydrodynamic model used in the present work can be found in our previous work [25, 26]. To give a short and clear description of the model, we introduce the main results of the model in the following.

The Landau hydrodynamic model [818] results approximately in a Gaussian rapidity distribution which does not exactly describe the experimental data [6, 19, 20]. We have revised the model to three sources: a central source which locates at midrapidity and covers the rapidity range as wide as possible, and a target (projectile) source which locates in the backward (forward) region and revises the rapidity distribution from the central source [25, 26]. The experimental rapidity distributions are then described by three Gaussian functions. And based on the three Gaussian functions, the experimental pseudorapidity distributions can be described by a method which makes a distinction between rapidity and pseudorapidity.

According to the Landau hydrodynamic model [818] and our revision [25, 26], the rapidity distribution, , of charged particles produced in a given source in high energy collisions can be described by a Gaussian function:where , , and denote the width, peak position, and normalization constant, respectively; , , and are for the central, target, and projectile sources, respectively. The experimental result is in fact a sum weighted by the three Gaussian functions.

The relation between and the squared speed-of-sound iswhere is the logarithmic Lorentz contraction factor which is independent of , denotes the center-of-mass energy per pair of nucleons in nucleus-nucleus collisions and is simplified to in - collisions, and denotes the rest mass of a proton. From (2), is expressed by using to beIn the above extractions of and from the rapidity distribution, we distinguish accurately the rapidity and pseudorapidity distributions [25, 26]. In the case of representing the pseudorapidity distribution in experiment, we also extract and from the hidden rapidity distribution. Our treatment ensures the method of extraction being concordant.

3. Results and Discussion

Figure 1 presents the rapidity distributions, , of negatively charged pions produced in - and 7Be-9Be collisions at the SPS. For - collisions, the incident momentums for Figures 1(a)1(e) are 20, 31, 40, 80, and 158 GeV/, respectively, which correspond to , 7.7, 8.8, 12.3, and 17.3 GeV, respectively. For 7Be-9Be collisions, the incident momentums for Figures 1(a)1(e) are , , , , and  GeV/, respectively, which correspond to , 7.7, 8.8, 11.9, and 16.8 GeV, respectively. The symbols represent the experimental data of the NA61/SHINE Collaboration [6, 7] and the curves are our results fitted by the revised (three-source) Landau hydrodynamic model. Different centrality classes (0–5%, 5–10%, 10–15%, and 15–20%) for 7Be-9Be collisions are presented by different symbols marked in the panels. The values of fit parameters, , (), , rapidity shift (), relative contribution of the target source (, the relative contribution of the projectile source), and normalization constant , are listed in Table 1 with the values of per degree of freedom (dof). The last two columns in Table 1 present the values of and [], which are calculated from (3) due to . Both and [] are in the units of , where is the speed of light in vacuum. One can see that the model describes approximately the experimental data in most cases.

To see clearly the relation between parameter and energy, Figures 2(a) and 2(b) give the dependences of main free parameter and extracted parameter on , respectively. The results corresponding to - collisions and 7Be-9Be collisions with different centrality classes are displayed by different symbols, which reflect a part of parameter values listed in Table 1. The line segments for - collisions and for 7Be-9Be collisions with centrality class 15–20% are given for guiding the eyes. One can see that increases with increase of , () does not show a change with , shows itself a local minimum (knee point) which indicates a softest point in the equation of state (EoS) at  GeV, and [] decreases with . The dependences of other parameters (, , , and ) on energy are not analyzed due to trivialness.

To study the knee point in detail, Figures 3(a) and 3(b) show and excitation functions (the dependences of and on ), respectively, where the results for - and Be-Be collisions at lower energies (≤17.3 GeV) are taken from Figure 2 (Table 1) and the results for -, proton-antiproton (-), copper-copper (Cu-Cu), gold-gold (Au-Au), and lead-lead (Pb-Pb) collisions at higher energies (≥19.6 GeV) are taken from our previous works [25, 26]. Different symbols represent different collisions, which are marked in the panels. The symbols with large size denote central collisions, and the symbols with small size denote noncentral collisions. The centrality classes for Cu-Cu collisions are from 0–3% to 50–55% [27]; for Au-Au collisions from 0–3% to 45–50% (maximum 40–45% at 19.6 GeV) [27]; and for Pb-Pb collisions from 0–5% to 20–30% [28]. The line segments for - collisions at SPS and higher energies are given for guiding the eyes, and the two dotted lines in Figure 3(b) show and 1/2, respectively. One can see that at higher energies increases with increase of and most are between 1/3 and 1/2 [a few are beyond 1/2 due to statistical fluctuations].

In particular, one can see some detailed laws in excitation function (the dependence of on ). (i) In the lower energy range, decreases with increase of , then the first knee point appears at  GeV, and afterwards increases with increase of . (ii) From 17.3 to 19.6 GeV, has a jump from ≤0.3 to 1/3–1/2. (iii) In the higher energy range, seems not related to , and a saturation appears. The first knee point is a local minimum which indicates a softest point in the EoS. As the predecessor of , increases always with increase of . We do not observe a softest point in the EoS in excitation function. These abundant phenomena should be related to the searching for the onset of deconfinement of the quarks and gluons in proton-proton collisions and the critical point of phase transition from hadronic matter to quark-gluon plasma (QGP) in nucleus-nucleus collisions.

The softest point (8.8 GeV) obtained in the present work is compatible with the previous works [29, 30] which used the Landau hydrodynamic model and the ultra-relativistic quantum molecular dynamics hybrid approach and indicated the softest point locating in the energy range from 4 to 9 GeV. Other works which study dependences of ratio of numbers of positive kaons and pions () [20, 31, 32], chemical freeze-out temperature () [31, 32], mean transverse mass minus rest mass () [31], and ratio of widths of experimental negative pion rapidity distribution and Landau hydrodynamic model prediction [ [32] on show knee point around 7-8 GeV. A wiggle in the excitation function of a specific reduced curvature of the net-proton rapidity distribution at midrapidity is expected in the energy range from 4 to 8 GeV [33, 34]. However, the searching for the onset of quark deconfinement and the critical point of QGP phase transition is a complex process [35]. The relation between the softest point and onset of deconfinement is still an open question. We need to analyze more other experimental data so that these results can be confirmed each other.

As for the jump of from ≤0.3 to 1/3–1/2 when increases from 17.3 to 19.6 GeV, we explain it as the result of changing the mean free path of produced particles, if the statistical fluctuation is excluded. At the SPS energies, the interacting system has a small density due to low collision energy and the produced particles have a large mean free path which results from a gas-like state and results in a small . The situation at the RHIC or LHC energies is opposite, where the interacting system has a large density due to high collision energy and the produced particles have a small mean free path which results from a liquid-like state and results in a large . Let denote the dimensionality of space. According to Refs. [36, 37], we have the relation of for massless particles. The particles at the gas-like state have a larger probability to appear in three-dimensional space, which results in the maximum being 1/3 which is the situation at the SPS energies. The particles at the liquid-like state have a larger probability to appear in two-dimensional space, which results in the maximum being 1/2 which is the situation at the RHIC and LHC energies.

To see the relations of and other parameters at the SPS energies which show the softest point in the EoS, Figures 4(a), 4(b), 4(c), and 4(d) present the dependences of on , , , and , respectively, where the closed circles, open circles, closed squares, open squares, and stars correspond to the collision energy being 6.3, 7.7, 8.8, 12.3 (11.9), and 17.3 (16.8) GeV, respectively, which are taken from Table 1. One can see that does not show an obvious dependence on and , though shows somehow a saturation. increases with increases of and due to the latter two increasing with increase of . The relation between and is the most important one among all the relations. Other relations such as the relations between and , , , as well as are less important. To present the relations between and , , , as well as is trivial due to being calculated from .

To study the most important relation between and in detail, Figure 5(a) displays a few examples of linear relations between and which reflect approximately the main area of parameter points in Figure 3(a). The solid lines corresponded to cases A, B, C, and D and can be expressed by , , , and , respectively. After conversion by (3), the four similar lines of show very different relations of given correspondingly in Figure 5(b). Cases A and B show obvious minimums, while cases C and D do not show any minimum. Decreasing the intercept in case A from 0.3 to 0.2 and 0.1, respectively, the corresponding results presented, respectively, by the dashed and dotted lines in Figure 5(a) have a small change, while the results presented, respectively, by the dashed and dotted curves in Figure 5(b) have a large change. Decreasing the slope in case B from 0.3 to 0.29 and 0.28, respectively, the corresponding results presented, respectively, by the dashed and dotted lines (curves) in Figures 5(a) and 5(b) have a small change. The results presented in Figure 5 show that the minimum in (3) appearing in some special conditions.

4. Conclusions

We summarize here our main observations and conclusions.(a)The rapidity distributions of negatively charged pions produced in - collisions and in 7Be-9Be collisions with different centrality classes at the SPS energies are analyzed by using the revised (three-source) Landau hydrodynamic model. The model results are in agreement with the experimental data measured by the NA61/SHINE Collaboration over an energy range from 6.3 to 17.3 GeV.(b)The values of squared speed-of-sound parameter are extracted from the widths of rapidity distributions. For the target (projectile) source, [] decreases obviously with increase of in the considered energy range, and () does not depend on . For the central source, shows a minimum (the softest point in the EoS) at 8.8 GeV when increases from 6.3 to 17.3 GeV, and increases obviously with increase of .(c)Combining with our previous works, one can see that is between 1/3 and 1/2 in most cases at higher energies (≥19.6 GeV). From 17.3 to 19.6 GeV, has a jump from ≤0.3 to 1/3–1/2. In the energy range from the minimum SPS energy to the maximum LHC energy, increases monotonously with increase of , though some fluctuations appear in the excitation function.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank Dr. Andrzej Wilczek at University of Silesia in Katowice, Poland, and Dr. Emil Aleksander Kaptur at CERN, Switzerland, for supplying us the numerical data of Be-Be collisions in Figure 1 which might be a little bit different from that published in ICPAQGP-2015 [7]. This work was supported by the National Natural Science Foundation of China under Grant no. 11575103 and the US DOE under contract DE-FG02-87ER40331.A008.