Multiplicity dependence of charged particle, ϕ meson, and multi-strange particle productions in p+p collisions at s = 200 GeV from PYTHIA simulation

Multiplicity dependence of charged particle, ϕ meson, and multi-strange particle productions in p+p collisions at <math id="M1"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV from PYTHIA simulation 1 first-author 1 1 corresp Corresponding author, ephy@ustc.edu.cn Corresponding author, ephy@ustc.edu.cn ephy@ustc.edu.cn 1 1 1 1 1 University of Science and Technology of China, Hefei 230026, China University of Science and Technology of China, Hefei 230026, China We report the multiplicity dependence of charged particle production for the π ± , K ± , p , <math id="M2"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> , and ϕ mesons at | y | < 1.0 in p+p collisions at <math id="M3"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV from a PYTHIA simulation. The impact of multiple parton interactions and gluon contributions is studied and found to be a possible source of the splitting of the particle yields as a function of p T with respect to the multiplicity. No obvious particle species dependence of the splitting is observed. The multiplicity dependence of the ratios K - / π - , K + / π + , <math id="M4"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> , p/π + , and <math id="M5"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> at mid-rapidity in p+p collisions is found to follow a tendency similar to that in Au+Au collisions at <math id="M6"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV at the Relativistic Heavy Ion Collider, indicating similar underlying initial production mechanisms despite the differences in the initial colliding systems. particle production heavy-ion collisions small system multiple parton interactions gluon contributions 1 Introduction 1 Searching for a novel form of nuclear matter with deconfined quarks and gluons created in ultrarelativistic heavy-ion collisions is the main goal of high-energy nuclear physics. The properties of this strongly coupled form of matter, the so-called quark–gluon plasma (QGP) [ 1 1 - 3 3 1-3 ], emerge by experimental comparison of elementary particle collisions. Measurements of particle production in proton–proton (p+p) collisions are critical to provide a baseline for understanding the interactions in the QGP created in heavy-ion collisions. In recent decades, there have been many measurements of the multiplicity or centrality dependence of particle production in heavy-ion collisions [ 4 4 - 7 7 4-7 ]. In particular, the production mechanism of hadrons containing strangeness is believed to be a signature of QGP formation in heavy-ion collisions [ 8 8 - 10 10 8-10 ]. To support this, it is worth studying the multiplicity dependence of strangeness production in elementary particle collisions without any medium effect. Recently, strangeness enhancement in high-multiplicity p+p collisions was observed in the ALICE experiment [ 11 11 ], which is similar to that observed in heavy-ion collisions [ 12 12 ], where a hot dense medium is created. This contradicts our knowledge that elementary particle collisions create cold, tiny systems; thus, it has attracted considerable interest. However, the p+p colliding system is always treated as a fundamental particle collision, and the particle production with respect to multiplicity is rarely studied, especially at a colliding energy of a few hundred gigaelectron volts at the Relativistic Heavy Ion Collider (RHIC). Thus, it is also of interest to investigate the multiplicity dependence of particle production in p+p collisions at RHIC energies to see if there is any similarity to that in Au+Au collisions. In this paper, we report the multiplicity dependence of particle production for π ± , K ± , p , and <math id="M7"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> in p+p collisions at <math id="M8"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV on the basis of a PYTHIA simulation. We study the multiplicity dependence of particle yields as a function of the transverse momentum ( p T ) for π ± , K ± , p , <math id="M9"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> , and ϕ and the ratios K - / π - , K + / π + , <math id="M10"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> , and <math id="M11"><mrow><mi>p</mi><mo>/</mo><msup><mi>π</mi><mo>+</mo></msup></mrow></math> . There are three important sources of the splitting of the p T distribution with multiplicity: jet fragmentation, multiple parton interactions (MPIs), and gluon contributions. Because the jet fragmentation effect has been discussed in Ref. [ 13 13 ], we focus on the effects of MPIs and gluon contributions. The particle ratios K - / π - , K + / π + , <math id="M12"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> , and p/π + with respect to multiplicity are presented and compared to those in Au+Au collisions at <math id="M13"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>N</mtext><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV measured by RHIC-STAR [ 4 4 , 5 5 ]. The <math id="M14"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> ratio in p+p collisions at <math id="M15"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV at RHIC energies is also presented and compared with experimental results in Au+Au collisions. The paper is organized as follows: Section 2 presents the simulation process and detailed PYTHIA settings. Simulation results for the charged particle and ϕ meson p T spectra, average 〈 p T 〉, and particle yield ratios, and related discussions, are presented in Sect. 3. Finally, Section 4 gives a summary. 2 Simulation processes 1 The PYTHIA program is widely used for event generation in high-energy physics to simulate multiparticle production in collisions between elementary particles [ 14 14 ]. In this work, PYTHIA version 6.416 is used. The p+p collision events at <math id="M16"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV are generated and categorized in six multiplicity bins, as shown in Fig. 1 Fig. 1 , to study the multiplicity dependence of π ± , K ± , p , <math id="M17"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> , ϕ , and multi-strange particle production. The initial configurations in PYTHIA are as follows: MSEL(1) (a switch to select between full user control and some preprogrammed alternatives) includes ISUB = 11 q i q j → q i q j , 12 <math id="M19"><mrow><msub><mi>q</mi><mi>i</mi></msub><msub><mrow><mover accent="true"><mi>q</mi><mo stretchy="true">¯</mo></mover></mrow><mi>i</mi></msub><mo>→</mo><msub><mi>q</mi><mi>k</mi></msub><msub><mrow><mover accent="true"><mi>q</mi><mo stretchy="true">¯</mo></mover></mrow><mi>k</mi></msub></mrow></math> , 13 <math id="M20"><mrow><msub><mi>q</mi><mi>i</mi></msub><msub><mrow><mover accent="true"><mi>q</mi><mo stretchy="true">¯</mo></mover></mrow><mi>i</mi></msub><mo>→</mo><mi>g</mi><mi>g</mi></mrow></math> , 28 q i g → q i g , 53 <math id="M21"><mrow><mi>g</mi><mi>g</mi><mo>→</mo><msub><mi>q</mi><mi>k</mi></msub><msub><mrow><mover accent="true"><mi>q</mi><mo stretchy="true">¯</mo></mover></mrow><mi>k</mi></msub></mrow></math> , 68 gg → gg , and 96 semihard QCD 2 → 2. If all processes are considered, an event is called a minimum bias (minibias) event. The gluon contribution is studied by changing the ISUB 28, 53, and 68 process setups. MSTP(81, 1): Turn multiple interactions on. MSTP(61, 1), MSTP(71, 1): Initial- and final-state QCD radiation is added to the above processes. MSTP(51, 7): CTEQ5L parton distribution function MSTP(33, 1): A common K factor is used, as stored in PARP(31). PARP(31, 1.5): ( D = 1.5) Common K factor for multiplying the differential cross section for hard parton–parton processes 3 Results and discussion 1 3.1 Ratios of particle production yields 2 The p+p collision events generated by PYTHIA produce charged particles ( π ± , K ± , p , <math id="M22"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> , e ± , and μ ± ) within | η | < 0.5 at <math id="M23"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV. The number of events is called the charged multiplicity, N ch . The multiplicity distribution is shown in Fig. 1 Fig. 1 . Six multiplicity bins are chosen so that there is an equal number of events in each bin. The events without any multiplicity selection are called minibias events. The ratio of the particle production yield normalized by 〈 N ch 〉 (the mean of N ch ) in each multiplicity bin over that in minibias events is used to study the event activity with respect to different multiplicities, which can be defined as <math display="block" id="M24"><mrow><msub><mi>R</mi><mrow><mtext>pp</mtext></mrow></msub><mo>=</mo><mfrac><mrow><mtext>d</mtext><mi>N</mi><mo>/</mo><mtext>d</mtext><msub><mi>p</mi><mtext>T</mtext></msub><mo stretchy="false">(</mo><mtext>mult</mtext><mn>,</mn><msub><mi>p</mi><mtext>T</mtext></msub><mo stretchy="false">)</mo><mo>/</mo><mo>〈</mo><msub><mi>N</mi><mrow><mtext>ch</mtext></mrow></msub><mo stretchy="false">(</mo><mtext>mult</mtext><mo stretchy="false">)</mo><mo>〉</mo></mrow><mrow><mtext>d</mtext><mi>N</mi><mo>/</mo><mtext>d</mtext><msub><mi>p</mi><mi>T</mi></msub><mo stretchy="false">(</mo><mtext>minibias</mtext><mn>,</mn><msub><mi>p</mi><mtext>T</mtext></msub><mo stretchy="false">)</mo><mo>/</mo><mo>〈</mo><msub><mi>N</mi><mrow><mtext>ch</mtext></mrow></msub><mo stretchy="false">(</mo><mtext>minibias</mtext><mo stretchy="false">)</mo><mo>〉</mo></mrow></mfrac><mn>.</mn></mrow></math> The multiplicity dependence of the transverse momentum spectra of the π - meson with and without MPIs and gluon contributions in p+p collisions at <math id="M25"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV within | y | < 1.0 is shown in Fig. 2 Fig. 2 . 〈 N ch 〉 for different multiplicity bins and different configurations is shown in Table 1 Table 1 . We observe a clear hardening of the p T spectra of the π - meson from low to high multiplicity, which is consistent with the results in Ref. [ 15 15 ], although the collision energies are different. Further, the p T spectra of kaons, protons, and ϕ mesons [ 16 16 ] versus the multiplicity are similar to those of the pion despite having different integrated yields, which are not shown here. 1 1 Mechanism 1 1 〈 N ch1 〉 1 1 〈 N ch2 〉 1 1 〈 N ch3 〉 1 1 〈 N ch4 〉 1 1 〈 N ch5 〉 1 1 〈 N ch6 〉 1 1 〈 N ch (minibias)〉 1 1 minibias w MPIs 1 1 1.50 1 1 3.41 1 1 5.38 1 1 7.65 1 1 10.80 1 1 14.99 1 1 3.14 1 1 minibias wo MPIs 1 1 1.48 1 1 3.38 1 1 5.31 1 1 7.45 1 1 10.52 1 1 14.84 1 1 2.55 1 1 gluon off w MPIs 1 1 1.51 1 1 3.40 1 1 5.36 1 1 7.62 1 1 10.77 1 1 14.95 1 1 2.99 1 1 gluon off wo MPIs 1 1 1.46 1 1 3.34 1 1 5.27 1 1 7.35 1 1 10.42 1 1 15.04 1 1 2.28 The R pp distributions of π ± , K ± , p , <math id="M27"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> , and ϕ mesons in p+p collisions at <math id="M28"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV are extracted using formula (1). The pion R pp distributions as a function of p T in different multiplicity bins for different initial production mechanisms are shown in Fig. 3 Fig. 3 . We observe clear R pp splitting with different production mechanisms, and the splitting becomes more obvious as p T increases. The MPIs suppress the splitting in the higher p T range. The contribution of gluon processes to the change in the relative momentum shape with respect to the multiplicity is small. We find similar conclusions for the R pp distributions of kaons, protons, and ϕ mesons, which are not shown here. The integrated R pp distributions for different particle species at p T > 1.3 GeV/c are shown in Fig. 4 Fig. 4 . The increase in the integrated R pp could be due to jet fragmentation, as described in Ref. [ 13 13 ]. As we can see, MPIs are the dominant source of R pp splitting suppression. The reason may be that the particle momenta become softer after multiple scatterings as energy is transferred to surrounding partons. Thus, the MPIs are the main source competing with jet fragmentation. We also find that gluon contributions have little impact on the R pp splitting. Furthermore, qualitatively, no obvious particle–antiparticle dependence of the R pp splitting is observed. More quantitative studies of the production difference between particle and antiparticle species are presented in Sect. 3.3. 3.2 Average transverse momenta 〈 p T 〉 2 In this section, the multiplicity dependence of 〈 p T 〉 for π - , K - , and <math id="M32"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> in p+p collisions is presented and compared with those in d+Au and Au+Au collisions at <math id="M33"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>N</mtext><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV from the STAR experiment [ 4 4 ]. Here the charged particle rapidity density (d N ch /d y ) is used instead of 〈 N ch 〉 in each multiplicity bin for an apple-to-apple comparison with the data. It is extracted by summing up the rapidity density of pions, kaons, protons, and antiprotons within | y | < 0.1. Fig. 5 Fig. 5 shows the 〈 p T 〉 distributions of π - , K - , and <math id="M34"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> as a function of d N ch /d y for different initial production mechanisms in p+p collisions at <math id="M35"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV. The 〈 p T 〉 distributions of pions, kaons, and protons in p+p collisions increase significantly with the multiplicity. The tendency is the same as that in Au+Au collisions but with a slightly larger slope. In p+p collisions, there are only hadronic processes, whereas the situation is much more complicated in Au+Au collisions; in addition to jet fragmentation and MPIs, radial flow and energy-loss mechanisms can modify the particle momentum, thus changing the 〈 p T 〉 distribution. However, the similar increasing tendency as a function of the charged particle density in both p+p and Au+Au collisions suggests that the baseline contributions of jet fragmentation and MPIs hold in different collision systems. 3.3 Particle ratios 2 In this section, the multiplicity dependence of the particle ratios π - , K - / K + , <math id="M38"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><mi>p</mi></mrow></math> , K + / π + , K - / π - , p / π + , and <math id="M39"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> in p+p collisions are presented and compared with those in Au+Au collisions. The antiparticle-to-particle ratios ( π - / π + , K - / K + , <math id="M40"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><mi>p</mi></mrow></math> ) as a function of the charged particle multiplicity within | y | < 0.1 for different initial production mechanisms in p+p collisions at <math id="M41"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV are shown in Fig. 6 Fig. 6 . In order to further understand the antiparticle and particle production mechanisms, the antiparticle-to-particle ratios are studied in different p T regions as a function of d N ch /d y . We observe that the π - π + ratio is independent of the multiplicity, and the K - / K + and <math id="M42"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><mi>p</mi></mrow></math> ratios depend slightly on the multiplicity in the entire p T region. However, in the high- p T region ( p T > 2 GeV/ c ), the dependence of the multiplicity in the low-multiplicity region is obvious for the π - / π + ratio and is stronger for the K - / K + and <math id="M43"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><mi>p</mi></mrow></math> ratios. The dependence becomes smaller when we switch off the gluon contributions. The ratios K + / π + , K - / π - , p / π + , and <math id="M45"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> for different multiplicities with different initial production mechanisms in p+p collisions at <math id="M46"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV are shown in Fig. 7 Fig. 7 , along with a comparison to those in Au+Au collisions at <math id="M47"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>N</mtext><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV [ 4 4 , 5 5 ]. We also compare them with the published STAR results in minimum-bias p+p collisions at the same collision energy. The ratios are consistent with the published STAR results within the uncertainties in p+p collisions [ 4 4 , 5 5 ]. The similar tendency of the particle ratios in p+p and Au+Au collisions within | y | < 0.1 indicates that similar production mechanisms hold in different collision systems. Likewise, we also study the ratios K + / π + , K - / π - , p / π + , and <math id="M50"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> with respect to the multiplicity in different p T regions. The results are shown in Figs. 8 Figs. 8 and 9 9 . Weak multiplicity dependence is observed for the integrated particle yield ratios in the entire p T region, where the low- p T d N /d y dominates. However, in the high- p T region, we see a clear difference between the particle ( K + , p )-to- π + and antiparticle ( K - , <math id="M51"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> )-to- π - ratios in the low-multiplicity region. This difference could be due to gluon contributions because it becomes smaller after the gluon contributions are switched off. This may suggest that the high- p T gluon jet compositions are different in terms of the particle-to- π + and antiparticle-to- π - ratios. Fig. 10 Fig. 10 shows the ratio <math id="M55"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>S</mtext><mn>0</mn></msubsup></mrow></math> as a function of p T at high multiplicity in p+p collisions at <math id="M56"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV, in comparison with those in central and peripheral Au+Au collisions with the same collision energy [ 6 6 ]. The ratio at high multiplicity in p+p collisions is consistent with the data for peripheral Au+Au collisions. However, the enhancements of the <math id="M57"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> ratio measured in central Au+Au collisions at the RHIC and p+p collisions at the Large Hadron Collider cannot be reproduced in the PYTHIA simulation, which is supposed that no hot medium is created. This provides additional support for coalescence hadronization [ 17 17 ] in the nuclear medium created in Au+Au collisions at the RHIC or even p+p collisions in TeV at reference to in the ALICE Nature article about enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions [ 11 11 ]. Furthermore, the distribution of the <math id="M60"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>S</mtext><mn>0</mn></msubsup></mrow></math> ratio as a function of d N ch /d y is studied and compared with data for Au+Au collisions. The d N ch /d y used for the Au+Au data is extracted from the published STAR results by polynomial extrapolation to full the number of participating nucleons ( N part ) coverage [ 4 4 ]. The distributions of the <math id="M61"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>S</mtext><mn>0</mn></msubsup></mrow></math> ratio as a function of d N ch /d y in p+p and Au+Au collisions are shown in Fig. 11 Fig. 11 . The ratio decreases as a function of d N ch /d y in p+p collisions, whereas it remains constant in Au+Au collisions at the same collision energy. We find that the <math id="M62"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> ratio in high-multiplicity events in 200 GeV p+p collisions in the PYTHIA simulation is comparable to that in most peripheral Au+Au collisions at the RHIC. The multiplicity dependence of the <math id="M63"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> ratio in p+p collisions shows a smooth connection to that in Au+Au collisions at <math id="M64"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>N</mtext><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV, which may indicate that some similar production mechanisms transfer p+p collisions to a hotter and denser system. 4 Summary 1 The multiplicity dependence of the π ± , K ± , p , <math id="M67"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover></mrow></math> , and ϕ meson production yields at | y | < 1.0 in proton–proton collisions at <math id="M68"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV from a PYTHIA simulation were presented. The R pp distributions showed obvious splitting as a function of multiplicity for all particle species, which is caused by jet fragmentation, MPIs, and gluon contributions. MPIs obviously suppress the splitting, whereas the gluon contributions play a less important role. The distributions of the K + / π + , K - / π - , p / π + , and <math id="M69"><mrow><mover accent="true"><mi>p</mi><mo stretchy="true">¯</mo></mover><mo>/</mo><msup><mi>π</mi><mo>−</mo></msup></mrow></math> ratios with respect to the multiplicity in p+p collisions were compared with those in Au+Au collisions within | y | < 0.1 and were found to show a similar tendency. This suggests that there are similar underlying initial production mechanisms in p+p and Au+Au collisions. We also presented the multiplicity dependence of the <math id="M70"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> ratio at | y | < 0.5 in p+p collisions at <math id="M71"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 200 GeV on the basis of a PYTHIA simulation. A smooth connection of the <math id="M72"><mrow><mi>Λ</mi><mo>/</mo><msubsup><mi>K</mi><mtext>s</mtext><mn>0</mn></msubsup></mrow></math> ratio in p+p collisions to that in Au+Au collisions was observed, which may provide a hint that a smaller hot medium may be created in such fundamental particle collisions with a sufficiently high initial energy density or that some similar production mechanisms make p+p collisions a hotter and denser system. References S. Edward , Strongly coupled quark-gluon plasma in heavy ion collisions . Rev. Mod. Phys . 89 , 035001 ( 2017 ). doi: 10.1103/RevModPhys.89.035001 http://doi.org/10.1103/RevModPhys.89.035001 Strongly coupled quark-gluon plasma in heavy ion collisions J. Adams et al . (STAR Collaboration) , Experimental and Theoretical Challenges in the Search for the Quark Gluon Plasma: The STAR Collaboration’s Critical Assessment of the Evidence from RHIC Collisions . Nucl. Phys. A 757 , 102 - 183 ( 2005 ). doi: 10.1016/j.nuclphysa.2005.03.085 http://doi.org/10.1016/j.nuclphysa.2005.03.085 Experimental and Theoretical Challenges in the Search for the Quark Gluon Plasma: The STAR Collaboration’s Critical Assessment of the Evidence from RHIC Collisions B. Muller , J. Schukraft , B. Wyslouch , First Results from Pb+Pb Collisions at the LHC . Ann. Rev. Nucl. Part. Sci . 62 , 361 - 386 ( 2012 ). doi: 10.1146/annurev-nucl-102711-094910 http://doi.org/10.1146/annurev-nucl-102711-094910 First Results from Pb+Pb Collisions at the LHC B.I. Abelev et al . (STAR Collaboration) , Systematic measurements of identified particle spectra in pp, d + Au, and Au + Au collisions at the STAR detector . Phys. Rev. C 79 , 034909 ( 2009 ). doi: 10.1103/PhysRevC.79.034909 http://doi.org/10.1103/PhysRevC.79.034909 Systematic measurements of identified particle spectra in pp, d + Au, and Au + Au collisions at the STAR detector J. Adams et al . (STAR Collaboration) , Identified Particle Distributions in pp and Au + Au Collisions at <math id="M73"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 92 , 112301 ( 2004 ). doi: 10.1103/PhysRevLett.92.112301 http://doi.org/10.1103/PhysRevLett.92.112301 Identified Particle Distributions in pp and Au + Au Collisions at sNN = 200 GeV G. Agakishiev et al . (STAR Collaboration) , Strangeness Enhancement in Cu-Cu and Au-Au Collisions at <math id="M74"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 108 , 072301 ( 2012 ). doi: 10.1103/PhysRevLett.108.072301 http://doi.org/10.1103/PhysRevLett.108.072301 Strangeness Enhancement in Cu-Cu and Au-Au Collisions at sNN = 200 GeV Abelev, B . et al . (ALICE Collaboration) , Corrigendum to ?Multi-strange baryon production at mid-rapidity in Pb?Pb collisions at <math id="M75"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 2.76 TeV . Phys. Lett. B 728 , 216 - 227 ( 2014 ); Corrigendum to ?Multi-strange baryon production at mid-rapidity in Pb?Pb collisions at sNN = 2.76 TeV Abelev, B . et al . (ALICE Collaboration) , Corrigendum to ?Multi-strange baryon production at mid-rapidity in Pb?Pb collisions at <math id="M76"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 2.76 TeV . Phys. Lett. B 734 , 409410 ( 2014 ). doi: 10.1016/j.physletb.2014.05.052 http://doi.org/10.1016/j.physletb.2014.05.052 Corrigendum to ?Multi-strange baryon production at mid-rapidity in Pb?Pb collisions at sNN = 2.76 TeV J. Rafelski , B. Muller , Strangeness Production in the Quark-Gluon Plasma . Phys. Rev. Lett . 48 , 1066 ( 1982 ). Strangeness Production in the Quark-Gluon Plasma J. Rafelski , B. Muller , Strangeness Production in the Quark-Gluon Plasma . Phys. Rev. Lett . 56 , 2334 ( 1986 ). doi: 10.1103/PhysRevLett.56.2334 http://doi.org/10.1103/PhysRevLett.56.2334 Strangeness Production in the Quark-Gluon Plasma J. Rafelski , in Elementary Hadronic Processes and Heavy Ion Interactions, ed. by J. Tran Thanh Van. Proceedings, 17th Rencontres De Moriond, Les Arcs, France, March 14-26, 1982. Vol. 2. Formation and Observables of the Quark-Gluon Plasma . Phys. Rept . 88 , 331 - 347 ( 1982 ). in Elementary Hadronic Processes and Heavy Ion Interactions, ed. by J. Tran Thanh Van. Proceedings, 17th Rencontres De Moriond, Les Arcs, France, March 14-26, 1982. Vol. 2. Formation and Observables of the Quark-Gluon Plasma P. Koch , B. Muller , J. Rafelski , Strangeness in relativistic heavy ion collisions . Phys. Rept . 142 , 167262 ( 1986 ). doi: 10.1016/0370-1573(86)90096-7 http://doi.org/10.1016/0370-1573(86)90096-7 Strangeness in relativistic heavy ion collisions J. Adam et al . (ALICE Collaboration) , Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions , Nature Physics 13 , 535 - 539 ( 2017 ). doi: 10.1038/nphys4111 http://doi.org/10.1038/nphys4111 Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions B. Abelev et al . (ALICE Collaboration) , <math id="M77"><mrow><msubsup><mi>K</mi><mi>S</mi><mn>0</mn></msubsup></mrow></math> and Λ Production in Pb-Pb Collisions at <math id="M78"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mi>N</mi><mi>N</mi></mrow></msub></mrow></msqrt></mrow></math> = 2.76 TeV . Phys. Rev. Lett . 111 , 222301 ( 2013 ). doi: 10.1103/PhysRevLett.111.222301 http://doi.org/10.1103/PhysRevLett.111.222301 KS0 and Λ Production in Pb-Pb Collisions at sNN = 2.76 TeV X. N. Wang , R. C. Hwa , Effect of jet production on the multiplicity dependence of average transverse momentum . Phys. Rev. D 39 , 187 ( 1989 ). doi: 10.1103/PhysRevD.39.187 http://doi.org/10.1103/PhysRevD.39.187 Effect of jet production on the multiplicity dependence of average transverse momentum T. Sjstrand , S. Mrenna , P. Skands , PYTHIA 6.4 Physics and Manual . JHEP 0605 026 ( 2006 ). doi: 10.1088/1126-6708/2006/05/026 http://doi.org/10.1088/1126-6708/2006/05/026 PYTHIA 6.4 Physics and Manual R. D. de Souza , Identified particle production in pp collisions at <math id="M79"><mrow><msqrt><mi>s</mi></msqrt></mrow></math> = 7 and 13 TeV measured with ALICE . J. Phys. Conf. Ser . 779 , 012071 ( 2017 ). doi: 10.1088/1742-6596/779/1/012071 http://doi.org/10.1088/1742-6596/779/1/012071 Identified particle production in pp collisions at s = 7 and 13 TeV measured with ALICE Y.J. Ye , J.H. Chen , Y.G Ma et al ., Φ-meson production at forward/backward rapidity in high-energy nuclear collisions from a multiphase transport model . Phys. Rev. C 93 , 044904 ( 2016 ). doi: 10.1103/PhysRevC.93.044904 http://doi.org/10.1103/PhysRevC.93.044904 Φ-meson production at forward/backward rapidity in high-energy nuclear collisions from a multiphase transport model V. Greco , C.M. Ko , Hadronization via Coalescence . Acta Phys.Hung. A 24 , 235 - 240 ( 2005 ). doi: 10.1556/APH.24.2005.1-4.32 http://doi.org/10.1556/APH.24.2005.1-4.32 Hadronization via Coalescence 10.1007/s41365-018-0469-y This work was supported by the Major State Basic Research Development Program in China (No. 2014CB845400), the National Natural Science Foundation of China (No. 11375184), the Youth Innovation Promotion Association fund of CAS (No. CX2030040079), the Ministry of Science and Technology (MoST) of China (No. 2016YFE0104800), and the Science and Technological Fund of Anhui Province for Outstanding Youth (No. 1808085J02). 2018-05-14 2018-05-29 2018-06-02 2018-08-02 2018-09-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0469-y http://dx.doi.org/10.1007/s41365-018-0469-y . 2018 Nuclear Science and Techniques 09 29