Abstract
First results on \(\hbox {K}/\pi \), \(\hbox {p}/\pi \) and K/p fluctuations are obtained with the ALICE detector at the CERN LHC as a function of centrality in \(\text{ Pb--Pb }\) collisions at \(\sqrt{s_\mathrm{{NN}}} =2.76\hbox { TeV}\). The observable \(\nu _{\mathrm{dyn}}\), which is defined in terms of the moments of particle multiplicity distributions, is used to quantify the magnitude of dynamical fluctuations of relative particle yields and also provides insight into the correlation between particle pairs. This study is based on a novel experimental technique, called the Identity Method, which allows one to measure the moments of multiplicity distributions in case of incomplete particle identification. The results for \(\hbox {p}/\pi \) show a change of sign in \(\nu _{\mathrm{dyn}}\) from positive to negative towards more peripheral collisions. For central collisions, the results follow the smooth trend of the data at lower energies and \(\nu _{\mathrm{dyn}}\) exhibits a change in sign for \(\hbox {p}/\pi \) and K/p.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
The theory of strong interactions, Quantum Chromodynamics (QCD), predicts that at sufficiently high energy density nuclear matter transforms into a deconfined state of quarks and gluons known as Quark–Gluon Plasma (QGP) [1, 2]. One of the possible signatures of a transition between the hadronic and partonic phases is the enhancement of fluctuations of the number of produced particles in the hadronic final state of relativistic heavy-ion collisions [3,4,5]. Event-by-event fluctuations and correlations may show critical behaviour near the phase boundary, including the crossover region where there is no thermal singularity, in a strict sense, associated with the transition from a QGP phase to a hadron-gas phase. A correlation analysis of event-by-event abundances of pions, kaons and protons produced in \(\text{ Pb--Pb }\) collisions at LHC energies may provide a connection to fluctuations of globally conserved quantities such as electric charge, strangeness and baryon number, and therefore shed light on the phase structure of strongly interacting matter [6].
In view of the predicted criticality signals at crossover for vanishing net-baryon densities [7], event-by-event fluctuations of relative particle yields are studied using the fluctuation measure \(\nu _\mathrm{dyn}[A,B]\) [8] defined in terms of moments of particle multiplicity distributions as
where \(N_{A}\) and \(N_{B}\) are the multiplicities of particles A and B measured event-by-event in a given kinematic range. The \(\nu _\mathrm{dyn}[A,B]\)Footnote 1 fluctuation measure contrasts the relative strength of fluctuations of species A and B to the relative strength of correlations between these two species. It vanishes when the particles A and B are produced in a statistically independent way [8, 9].
This study at LHC energies is of particular importance for establishing the energy and system size dependence of \(\nu _{\mathrm{dyn}}\) in order to understand the trend observed at lower collision energies from the RHIC Beam Energy Scan (BES) results reported by the STAR collaboration [10]. Furthermore, the advantage of this fluctuation measurement is its robustness against non-dynamical contributions such as those stemming from participant nucleon fluctuations and finite particle detection efficiencies [8, 11]. Measurements of the \(\nu _{\mathrm{dyn}}\) observable for net-charge fluctuations were already published by ALICE [12]. Moreover, for identified particles, it was measured at the Super Proton Synchrotron (SPS) [13] and at the Relativistic Heavy-Ion Collider (RHIC) [10] in \(\text{ Pb--Pb }\) and Au–Au collisions, respectively. The ALICE detector at the LHC is ideally suited to extend these measurements to higher collision energies. In particular, the excellent charged-particle tracking and particle identification (PID) capabilities in the central barrel of the detector allow for a precise and differential event-by-event analysis at midrapidity and low transverse momentum (\(p_{\mathrm{T}}\)).
The paper is organized as follows. In Sect. 2, details about the ALICE detector setup and the dataset are given. Section 3 discusses the event and track selection criteria, particle identification procedure, and the analysis method. Estimates of statistical and systematic uncertainties are given in Sect. 4. Results on \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\), \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) and \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) in \(\text{ Pb--Pb }\) collisions at \(\sqrt{s_\mathrm{{NN}}} =2.76\hbox { TeV}\) are presented in Sect. 5, and finally Sect. 6 summarizes the measurements presented in this paper.
2 Experimental setup and dataset
ALICE is a general-purpose detector system designed, in particular, for the study of collisions of heavy ions at the LHC. The design, components, and performance of the ALICE detector have been reported elsewhere [14, 15]. The ALICE detector is comprised of several detector components organized into a central barrel detection system and forward/backward detectors. The main tracking and PID devices in the central barrel of the experiment are the Inner Tracking System (ITS) and the Time Projection Chamber (TPC), which are operated inside a large solenoidal magnet with \(B=0.5\) T. Two forward scintillator arrays V0-A and V0-C are located on either side of the interaction point and cover the pseudorapidity (\(\eta \)) intervals \(2.8<\eta <5.1\) and \(-3.7<\eta <-1.7\). The V0 detectors and the two neutron Zero Degree Calorimeters (ZDC), placed at ± 114 m from the interaction point, were used for triggering and event selection.
The ITS-TPC tracking system covers the midrapidity region and provides charged-particle tracking and momentum reconstruction down to \(p_{\mathrm{T}} =100\,\hbox {MeV}/c\). The ITS is employed to reconstruct the collision vertex with high precision and to reject charged particles produced in secondary vertices.
The analysis presented in this paper is based on about 13 million minimum-bias \(\text{ Pb--Pb }\) collisions at \(\sqrt{s_\mathrm{{NN}}} =2.76\hbox { TeV}\) collected in the year 2010. The minimum-bias trigger condition is defined by the coincidence of hits in both V0 detectors. In the offline event selection, V0 and ZDC timing information is used to reject beam-gas background and parasitic beam-beam interactions. The definition of the collision centrality is based on the charged-particle multiplicity measured in the V0 detectors [14], which can be related to collision geometry and the number of participating nucleons through a Monte-Carlo (MC) simulation based on a Glauber model [16].
3 Data analysis
3.1 Event and track selection
Charged particles reconstructed in the TPC with full azimuthal acceptance and in the pseudorapidity range of \(|\eta |< 0.8\) were used in this analysis. The momentum range was restricted to \(0.2<p<1.5\hbox { GeV}/c\) in order to minimize systematic uncertainties arising from the overlap of the \(\hbox {d}E/\hbox {d}x\) distributions. Furthermore, the following track selection criteria were applied to guarantee optimal \(\hbox {d}E/\hbox {d}x\) and momentum resolution, which are crucial for precise particle identification. Charged-particle tracks were accepted in this analysis when they have at least 80 out of a maximum of 159 reconstructed space points in the TPC, and the \(\chi ^{2}\) per space point from the track fit is less than 4. Daughter tracks from reconstructed secondary weak-decay kink topologies were rejected. Additional suppression of secondary particles was achieved by restricting the distance-of-closest-approach (DCA) of the extrapolated trajectory to the primary vertex position to less than 2 cm along the beam direction. In the transverse plane the restriction in the DCA depends on \(p_{\mathrm{T}}\) in order to take into account the \(p_{\mathrm{T}}\) dependence of the impact parameter resolution [17]. The remaining contamination after the DCA cuts is typically less than 10% for the momentum range covered in this work [18].
3.2 Identity method
The standard approach of finding the moments \(\langle N_{A} \rangle \), \(\langle N_{B} \rangle \), \(\langle N_{A}(N_{A}-1) \rangle \) and \(\langle N_{B}(N_{B}-1) \rangle \) is to count the number of particles \(N_{A}\) and \(N_{B}\) event-by-event and calculate averages over the dataset. However, this approach suffers from incomplete particle identification due to overlapping \(\hbox {d}E/\hbox {d}x\) distribution functions, which could be circumvented by either selecting suitable phase-space regions or by using additional detector information such as time-of-flight measurements. These procedures reduce the overall phase-space coverage and detection efficiencies. The present study is based on the Identity Method [19,20,21] which overcomes the misidentification problem.
The Identity Method was proposed in Ref. [19] as a solution to the misidentification problem for the analysis of events with two different particle species. In Ref. [20], the method was developed further to calculate the second moments of the multiplicity distributions of more than two particle species. Subsequently, in Ref. [21], it was generalized to the first and higher moments of the multiplicity distributions for an arbitrary number of particle species. The first experimental results using the Identity Method were published by the NA49 collaboration [13].
Instead of counting every detected particle event-by-event, the Identity Method follows a probabilistic approach using two basic experimental per-track and per-event observables, \(\omega \) and W, respectively. They are defined as
where \(x_{i}\) stands for the \(\hbox {d}E/\hbox {d}x\) of a given track i, \(\rho _{j}(x)\) is the \(\hbox {d}E/\hbox {d}x\) distribution of particle species j within a given phase-space bin and N(n) is the number of tracks in the nth event. The quantity \(\omega _{j}(x_{i})\) represents the probability that particle i is of type j. Thus, in case of perfect particle identification, one expects \(W_{j} = N_{j}\), while this does not hold in case of overlapping \(\hbox {d}E/\hbox {d}x\) distributions. Figure 1 shows the \(\omega \) and W distributions for pions, kaons and protons in the momentum interval of \(0.3<p<0.8\) GeV/c. The W distribution of protons shows a discrete structure because proton \(\hbox {d}E/\hbox {d}x\) distributions have the least overlap.
The moments of the W distributions can be constructed directly from experimental data. The Identity Method calculates the moments of the particle multiplicity distributions by unfolding the moments of the W distributions with the following matrix operation
where \(\langle \overrightarrow{W} \rangle \) and \(\langle \overrightarrow{N} \rangle \) are the vectors of the moments of W quantities and unknown true multiplicity distributions, respectively. The response matrix \(\text {A}\) is defined by the \(\omega \) quantities. A detailed description of the technique and a demonstration of its robustness can be found in Refs. [20, 21].
The \(\hbox {d}E/\hbox {d}x\) measurements used as the only input for the Identity Method are obtained from the TPC, which provides a momentum resolution of better than 2% and a single-particle detection efficiency of up to 80% for the kinematic range considered in this paper [14]. The Identity Method employs fits of inclusive \(\hbox {d}E/\hbox {d}x\) distributions for the calculation of the \(\omega \) probabilities entering Eq. 3. Since the overlap regions in the \(\hbox {d}E/\hbox {d}x\) distributions are also properly taken into account, a very good description of the inclusive \(\hbox {d}E/\hbox {d}x\) spectra, and therefore an excellent understanding of the TPC detector response, is required over the full momentum range covered in this analysis. To this end, the \(\hbox {d}E/\hbox {d}x\) distributions of pre-selected samples of pions, protons and electrons, identified by the reconstruction of \(\mathrm {K}_{S}^{0}\) and \(\varLambda \) decays and photon conversions, were fitted with a generalized Gaussian function of the form:
where A, \(\mu \), \(\sigma \), \(\alpha \) and \(\beta \) stand for the abundance, mean, width, skewness and kurtosis of the distribution, respectively. The detector response functions obtained in this way were used later to fit the inclusive \(\hbox {d}E/\hbox {d}x\) spectra. To cope with the dependencies of the \(\hbox {d}E/\hbox {d}x\) on the track angle and particle multiplicity, fits were performed over the entire pseudorapidity range of \(|\eta |<0.8\) in steps of 0.1 units for each centrality class. Moreover, the momentum intervals were chosen narrow enough to minimize the effect of the momentum dependence on \(\hbox {d}E/\hbox {d}x\), most particularly at low momenta where the magnitude of \(\hbox {d}E/\hbox {d}x\) varies rapidly with the momentum. An example of a \(\hbox {d}E/\hbox {d}x\) distribution in a given phase-space bin and the corresponding fits are shown in Fig. 2.
4 Statistical and systematic uncertainties
The statistical uncertainties were determined by the number of events in this analysis and the finite number of tracks in each event. The number of events also affects the uncertainty of the shape of the inclusive \(\hbox {d}E/\hbox {d}x\) spectra, which is determined by a fit. This uncertainty enters the calculation of \(\omega \) and W, and finally the computation of the moments of multiplicity distributions with the Identity Method. Since standard error propagation is impractical given the rather complicated numerical derivation of the final result, the subsample method was chosen to evaluate the statistical uncertainties. To this end, the data set was subdivided into \(n=\) 25 random subsamples i. The \(\nu _{\mathrm{dyn}}\) values were reconstructed for each subsample and the statistical uncertainty was obtained according to
where
The summary of all sources of systematic uncertainties is shown in Table 1 and in the next paragraphs the main contributors to the systematics are detailed.
The largest contribution to the total systematic uncertainty is from the fits of the measured particle \(\hbox {d}E/\hbox {d}x\) distributions. The quality of the fits was monitored by Kolmogorov–Smirnov (K–S) and \(\chi ^{2}\) tests. To study the influence of possible systematic shifts in the fit parameters on \(\nu _{\mathrm{dyn}}\), the fit parameters of each particle in the overlap regions were varied by about ± 0.5 %, which defines the boundaries where the K–S test fails at 90% confidence level. The observed maximum variations range from about 7% to 15% for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) and \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\), respectively.
Even though \(\nu _{\mathrm{dyn}}\) is known to be robust against detection-efficiency losses, it may show an explicit dependence if the detector response functions differ from Binomial or the efficiencies exhibit large variations with detector occupancy [8]. Therefore, one also has to investigate the uncertainty resulting from the detection efficiency losses. For that, the \(\nu _{\mathrm{dyn}}\) results reconstructed from a full Monte Carlo simulation of HIJING [22, 23] events employing a GEANT3 [24] implementation of the ALICE detector were compared to the analysis at the generator level, where in both generated and reconstructed levels perfect PID information was used. The resulting systematic uncertainty from the finite tracking efficiency is less than 6%.
The systematic uncertainties due to the track selection criteria were estimated by a variation of the selection ranges. The systematics from contamination of weak decays and other secondary particles were obtained by varying the DCA cuts. Other contributions to the total systematic uncertainty arise from the cuts applied on the maximum distance of the reconstructed vertex to the nominal interaction point along the beam axis, the number of required TPC space points per track and the \(\chi ^{2}\) per degree of freedom of the track fit. Moreover, the effect of the magnetic field polarity was investigated by separate analyses of data taken under two polarities. Neither of these contributions to the total systematic uncertainty exceeds 5%. The total systematic uncertainty was obtained by adding in quadrature the individual maximum systematic variations from these different contributions.
5 Results
5.1 Centrality dependence and comparison to models
In this section, the results are presented as a function of collision centrality and compared to calculations with the HIJING [22, 23] and AMPT [25] models. The unscaled values of \(\nu _{\mathrm{dyn}}\) for different combinations of particles in each centrality class, together with the final statistical and systematic uncertainties, are given in Table 2. Due to the intrinsic multiplicity dependence of \(\nu _{\mathrm{dyn}}\), discussed in Refs [26, 27], the values of \(\nu _{\mathrm{dyn}}\) were scaled further by the charged-particle multiplicity density at midrapidity, \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \). The fully corrected experimental \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \) values were taken from Ref. [18]. Figure 3 shows measured values of \(\nu _{\mathrm{dyn}}\) scaled by \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \) as a function of the collision centrality expressed in terms of \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \). The values for \(\nu _{\mathrm{dyn}}\) and \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \) for HIJING and AMPT were calculated by using corresponding particle multiplicities at the generator level within the same experimental acceptance. A flat behaviour is expected in this representation if a superposition of independent particle sources is assumed, as in the Wounded Nucleon Model (WNM) [28].
Measured values of \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\) and \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) are positive across the entire centrality range, while \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) is negative for the most peripheral collisions and changes sign at mid-central collisions. The centrality dependencies observed in \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) and \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) are similar in shape, being flat from central to mid-central collisions and systematically decreasing for the most peripheral ones. In contrast, \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\) is almost independent of centrality from most peripheral to mid-central collisions and rises as the centrality increases. A similar qualitative behaviour is also observed for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\) within the kinematic range of \(|\eta |<1\) and \(0.2<p<0.6\hbox { GeV}/c\) as measured in Au–Au collisions at \(\sqrt{s_\mathrm{{NN}}} =200\hbox { GeV}\) by the STAR collaboration. The difference in the absolute values is, to a large extent, due to the increase in \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \) by almost a factor of two between the two collision energies. The same argument holds true for the most central STAR data at \(62.4\hbox { GeV}\), although the centrality dependence is rather flat in this case [27]. The overall behaviour is defined by the interplay between correlation and fluctuation terms encoded in the definition of the \(\nu _{\mathrm{dyn}}\) observable. To disentangle these terms, one needs a dedicated study focusing on separate charge combinations, which also makes it possible to investigate contributions from resonance decays and global charge conservations.
An important characteristic of HIJING is that it treats nucleus-nucleus collisions as an independent superposition of nucleon-nucleon interactions. As such, it does not incorporate mechanisms for final-state interactions among the produced particles and therefore phenomena such as equilibrium and collectivity do not occur. The AMPT calculations are performed with three different settings including (1) string melting, (2) hadronic rescattering, and (3) string melting and hadronic rescattering. All three versions of the AMPT model presented here use hard minijet partons and soft strings from HIJING as initial conditions. Partonic evolution is described by Zhang’s parton cascade (ZPC) [29] which is followed by a hadronization process. In the last step, hadronic rescattering and the decay of resonances takes place. In the default AMPT model, after minijet partons stop interacting with other partons, they are combined with their parent strings to form excited strings, which are then converted to hadrons according to the Lund string fragmentation model [25]. However, in the string melting scenario, instead of employing the Lund string fragmentation mechanism, hadronization is modeled via a quark coalescence scheme by combining two nearest quarks into a meson and three nearest quarks (antiquarks) into a baryon (antibaryon). This ultimately reduces the correlation between produced hadrons.
HIJING produces positive values for the three particle pair combinations and does not exhibit any non-monotonic behaviour as a function of centrality, even though it implements exact global conservation laws. In contrast, hadronic rescattering produces additional resonances at the hadronization phase thereby introducing additional correlations between particles [25]. Consequently, the AMPT configuration with hadronic rescattering drives the \(\nu _{\mathrm{dyn}}\) results towards negative values as the collision centrality increases. In particular, for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\), contrary to the data, it predicts negative values. On the other hand, the AMPT version with string melting shows a weak centrality dependence for the three particle pair combinations. None of the models investigated in this work give a reasonable quantitative description of the measured data.
5.2 Energy dependence
Values of \(\nu _{\mathrm{dyn}}\) measured in this work for the most central Pb–Pb collisions were compared to NA49 and STAR data in Fig. 4. Measurements from NA49 and STAR show a smooth evolution of \(\nu _{\mathrm{dyn}}\) with collision energy and do not reveal any indications for critical behaviour in the range \(6.3< \sqrt{s_\mathrm{NN}} < 200\) GeV. The apparent differences between NA49 and STAR data for \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) and \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\) at \(\sqrt{s_\mathrm{NN}} <10\) GeV were traced back in Ref. [13] to the dependence of \(\nu _{\mathrm{dyn}}\) on the detector acceptance. Above this energy, both experiments report positive values for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\), and a weak dependence on the collision energy, whereas \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) is negative and approaches zero as the collision energy increases.
ALICE data are positive for the three particle pair combinations and follow the trend observed at lower energies, involving a sign change for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) and \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) as a function of energy. Such a change of sign has been predicted by transport models HSD [30] and UrQMD [31] in the RHIC energy regime [10]. Since neither HSD nor UrQMD explicitly include the quark and gluon degrees of freedom, this observation can be attributed to the particular realization of the string and resonance dynamics used in the models [30]. Additionally, HIJING and AMPT model calculations at LHC energies predict positive values except for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) in the AMPT configuration with hadronic rescattering and without string melting. To understand the difference between the STAR and ALICE results, the acceptance dependence of \(\nu _{\mathrm{dyn}}\) was also investigated with the ALICE data by varying the phase-space coverage. Opening the pseudorapidity window from \(|\eta |<0.8\) up to \(|\eta |<1\) yields a reduction in \(\nu _{\mathrm{dyn}}\) of 10-20%. However, this reduction is insufficient to explain the difference between the ALICE and STAR results, most particularly the sign change with increasing energy.
6 Summary
In summary, measurements of \(\nu _{\mathrm{dyn}}\) in \(\text{ Pb--Pb }\) collisions at \(\sqrt{s_\mathrm{{NN}}} =2.76\hbox { TeV}\) for three specific particle pair combinations using the Identity Method were presented. Values of \(\nu _{\mathrm{dyn}}\), scaled by the charged-particle density at midrapidity \(\hbox {d}N_{\mathrm{ch}}/\hbox {d}\eta \), exhibit finite variations with collision centrality. This is in contrast to predictions by HIJING, which, for all three pair combinations, show essentially constant as well as positive values. The results for \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\) and \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) are positive across the entire centrality range, while \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\) changes sign from positive to negative towards more peripheral collisions suggesting differences in the production dynamics of these pairs. The centrality dependence of \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {K}]\) shows a similar behaviour, increasing with centrality, as measured in Au–Au collisions at \(\sqrt{s_\mathrm{{NN}}} =200\hbox { GeV}\) by the STAR collaboration, while the data at \(\sqrt{s_\mathrm{{NN}}} =62.4\hbox { GeV}\) shows no centrality dependence. Comparisons with calculations from the AMPT model, using three distinct configurations, show that AMPT is unable to reproduce measured data in this work. Calculations with quark coalescence show only a very slight centrality dependence and no sign changes. On the other hand, AMPT values with hadronic rescattering and no quark coalescence decrease significantly with increasing collision centrality and exhibit a sign change towards central collisions in the case of \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\). The evolution of \(\nu _{\mathrm{dyn}}\) with collision energy shows that the particle production dynamics changes significantly from that observed at lower energies. Values of \(\nu _{\mathrm{dyn}}\) measured with all three pair combinations follow a smooth continuation of the data measured by STAR and exhibit a change in sign for \(\nu _{\mathrm{dyn}}[\mathrm {p},\mathrm {K}]\) and \(\nu _{\mathrm{dyn}}[\pi ,\mathrm {p}]\). The analysis of \(\nu _{\mathrm{dyn}}\) with enlarged acceptance shows that the magnitude of \(\nu _{\mathrm{dyn}}\) depends on the kinematical limits but the change appears too small to explain the difference with the STAR results. A more detailed analysis of fluctuations with charge and species specific pairs is required to fully characterize the particle production dynamics in heavy-ion collisions and understand, in particular, the origin of the sign changes reported in this work.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The numerical values of the data points will be uploaded to HEPData.]
Notes
In this study, \(\nu _\mathrm{dyn}[A,B]\) was taken to be \(\nu _\mathrm{dyn}[A+\overline{A},B+\overline{B}]\), where \(\overline{A}\) and \(\overline{B}\) are the anti-particles of A and B, respectively.
References
J.C. Collins, M.J. Perry, Superdense matter: neutrons or asymptotically free quarks? Phys. Rev. Lett. 34, 1353 (1975)
E.V. Shuryak, Quantum chromodynamics and the theory of superdense matter. Phys. Rep. 61, 71–158 (1980)
M.A. Stephanov, K. Rajagopal, E.V. Shuryak, Signatures of the tricritical point in QCD. Phys. Rev. Lett. 81, 4816–4819 (1998). arXiv:hep-ph/9806219 [hep-ph]
E.V. Shuryak, M.A. Stephanov, When can long range charge fluctuations serve as a QGP signal? Phys. Rev. C 63, 064903 (2001). arXiv:hep-ph/0010100 [hep-ph]
HotQCD Collaboration, A. Bazavov, et al., Fluctuations and Correlations of net baryon number, electric charge, and strangeness: a comparison of lattice QCD results with the hadron resonance gas model. Phys. Rev. D 86, 034509 (2012). arXiv:1203.0784 [hep-lat]
V. Koch, Hadronic fluctuations and correlations, in R. Stock (ed.), Relativistic Heavy Ion Physics (Springer, Heidelberg, 2010), pp. 626–652. http://materials.springer.com/lb/docs/sm_lbs_978-3-642-01539-7_20
A. Bazavov et al., The chiral and deconfinement aspects of the QCD transition. Phys. Rev. D 85, 054503 (2012). arXiv:1111.1710 [hep-lat]
C. Pruneau, S. Gavin, S. Voloshin, Methods for the study of particle production fluctuations. Phys. Rev. C 66, 044904 (2002). arXiv:nucl-ex/0204011 [nucl-ex]
P. Christiansen, E. Haslum, E. Stenlund, Number-ratio fluctuations in high-energy particle production. Phys. Rev. C 80, 034903 (2009). arXiv:0902.4788 [hep-ex]
STAR Collaboration, N.M. Abdelwahab et al., Energy Dependence of \(K/\pi \), \(p/\pi \), and \(K/p\) Fluctuations in Au+Au Collisions from \(\sqrt{s_{NN}}\) = 7.7 to 200 GeV. Phys. Rev. C 92(2), 021901 (2015). arXiv:1410.5375 [nucl-ex]
P. Braun-Munzinger, A. Rustamov, J. Stachel, Bridging the gap between event-by-event fluctuation measurements and theory predictions in relativistic nuclear collisions. Nucl. Phys. A 960, 114–130 (2017). arXiv:1612.00702 [nucl-th]
ALICE Collaboration, B. Abelev et al., Net-charge fluctuations in Pb–Pb collisions at \(\sqrt{s}_{NN} = 2.76\) TeV. Phys. Rev. Lett. 110(15), 152301 (2013). arXiv:1207.6068 [nucl-ex]
T. Anticic et al., Phase-space dependence of particle-ratio fluctuations in Pb + Pb collisions from 20 A to 158 A GeV beam energy. Phys. Rev. C 89(5), 054902 (2014). arXiv:1310.3428 [nucl-ex]
ALICE Collaboration, B.B. Abelev et al., Performance of the ALICE Experiment at the CERN LHC. Int. J. Mod. Phys. A 29, 1430044 (2014). arXiv:1402.4476 [nucl-ex]
ALICE Collaboration, K. Aamodt et al., The ALICE experiment at the CERN LHC. JINST 3, S08002 (2008)
C. Loizides, J. Nagle, P. Steinberg, Improved version of the PHOBOS Glauber Monte Carlo, SoftwareX (1–2), 13–18 (2015) arXiv:1408.2549 [nucl-ex]. https://doi.org/10.1016/j.softx.2015.05.001
ALICE Collaboration, K. Aamodt, et al., Suppression of charged particle production at large transverse momentum in central Pb–Pb collisions at \(\sqrt{s_{NN}} =\) 2.76 TeV. Phys. Lett. B 696, 30–39 (2011). arXiv:1012.1004 [nucl-ex]
ALICE Collaboration, B. Abelev, et al., Centrality dependence of \(\pi \), K, p production in Pb–Pb collisions at \(\sqrt{s_{NN}}\) = 2.76 TeV. Phys. Rev. C 88, 044910 (2013). arXiv:1303.0737 [hep-ex]
M. Gazdzicki, K. Grebieszkow, M. Mackowiak, S. Mrowczynski, Identity method to study chemical fluctuations in relativistic heavy-ion collisions. Phys. Rev. C 83, 054907 (2011). arXiv:1103.2887 [nucl-th]
M.I. Gorenstein, Identity method for particle number fluctuations and correlations. Phys. Rev. C 84, 024902 (2011). arXiv:1106.4473 [nucl-th]
A. Rustamov, M.I. Gorenstein, Identity method for moments of multiplicity distribution. Phys. Rev. C 86, 044906 (2012). arXiv:1204.6632 [nucl-th]
M. Gyulassy, X.-N. Wang, HIJING 1.0: a Monte Carlo program for parton and particle production in high-energy hadronic and nuclear collisions. Comput. Phys. Commun. 83, 307 (1994). arXiv:nucl-th/9502021 [nucl-th]
W.-T. Deng, X.-N. Wang, R. Xu, Hadron production in p+p, p+Pb, and Pb+Pb collisions with the HIJING 2.0 model at energies available at the CERN Large Hadron Collider. Phys. Rev. C83, 014915 (2011). arXiv:1008.1841 [hep-ph]
R. Brun, F. Bruyant, F. Carminati, S. Giani, M. Maire, A. McPherson, G. Patrick, L. Urban, GEANT detector description and simulation tool, CERN-W5013. Tech. rep., CERN (1994). http://cds.cern.ch/record/1082634
Z.-W. Lin, C.M. Ko, B.-A. Li, B. Zhang, S. Pal, A multi-phase transport model for relativistic heavy ion collisions. Phys. Rev. C 72, 064901 (2005). arXiv:nucl-th/0411110 [nucl-th]
V. Koch, T. Schuster, On the energy dependence of K/pi fluctuations in relativistic heavy ion collisions. Phys. Rev. C 81, 034910 (2010). arXiv:0911.1160 [nucl-th]
S.T.A.R. Collaboration, B.I. Abelev et al., K/pi fluctuations at relativistic energies. Phys. Rev. Lett. 103, 092301 (2009). arXiv:0901.1795 [nucl-ex]
A. Bialas, M. Bleszynski, W. Czyz, Multiplicity distributions in nucleus–nucleus collisions at high-energies. Nucl. Phys. B 111, 461–476 (1976)
B. Zhang, ZPC 1.0.1: a Parton cascade for ultrarelativistic heavy ion collisions. Comput. Phys. Commun. 109, 193–206 (1998). arXiv:nucl-th/9709009 [nucl-th]
M.I. Gorenstein, M. Hauer, V.P. Konchakovski, E.L. Bratkovskaya, Fluctuations of the K/pi ratio in nucleus–nucleus collisions: statistical and transport models. Phys. Rev. C 79, 024907 (2009). arXiv:0811.3089 [nucl-th]
M. Bleicher et al., Relativistic hadron hadron collisions in the ultrarelativistic quantum molecular dynamics model. J. Phys. G 25, 1859–1896 (1999). arXiv:hep-ph/9909407 [hep-ph]
Acknowledgements
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Universidade Federal do Rio Grande do Sul (UFRGS), Financiadora de Estudos e Projetos (Finep) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil; Ministry of Science & Technology of China (MSTC), National Natural Science Foundation of China (NSFC) and Ministry of Education of China (MOEC) , China; Ministry of Science, Education and Sport and Croatian Science Foundation, Croatia; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research | Natural Sciences, the Carlsberg Foundation and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Energie Atomique (CEA) and Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Research and Religions, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE), Department of Science and Technology, Government of India (DST), University Grants Commission, Government of India (UGC) and Council of Scientific and Industrial Research (CSIR), India; Indonesian Institute of Science, Indonesia; Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionale di Fisica Nucleare (INFN), Italy; Institute for Innovative Science and Technology , Nagasaki Institute of Applied Science (IIST), Japan Society for the Promotion of Science (JSPS) KAKENHI and Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education and National Science Centre, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics and Romanian National Agency for Science, Technology and Innovation, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education and Science of the Russian Federation and National Research Centre Kurchatov Institute, Russia; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba and Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Spain; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; National Science and Technology Development Agency (NSDTA), Suranaree University of Technology (SUT) and Office of the Higher Education Commission under NRU project of Thailand, Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America.
Author information
Authors and Affiliations
Consortia
Additional information
See Appendix A for the list of collaboration members.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Acharya, S., Adamová, D., Adolfsson, J. et al. Relative particle yield fluctuations in \(\text{ Pb--Pb }\) collisions at \(\sqrt{s_\mathrm{{NN}}} =2.76\hbox { TeV}\). Eur. Phys. J. C 79, 236 (2019). https://doi.org/10.1140/epjc/s10052-019-6711-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-019-6711-x