arxiv:8.0889v [nucl-th] 2 Nov 208 Strongly intensive fluctuations and correlations in ultrarelativistic nuclear collisions in the model with string fusion Vladimir Kovalenko, Saint Petersburg State University Introduction Abstract. The several types of strongly intensive correlation variables are studied in nuclear collisions at LHC energy. These quantities are expected not to depend on centrality class width. They were calculated in dipole-based partonstring Monte Carl model with string fusion. The centrality dependence of mean transverse momentum correlation coefficient and strongly intensive quantity Σ between multiplicity and P T were obtained. Dynamical charge fluctuation ν dyn was also calculated and compared with experimental data. It is shown that string fusion improves agreement with the experiment. The exploration of the QCD phase diagram and a search for new phenomena is one of the main subjects of heavy ion physics. The experimental investigation of the phase diagram of a strongly interacting matter is related to the study of the fluctuations and correlations of observables in the nuclear collisions at high energy[]. Around the critical point, anomalies, such as enhanced fluctuations are expected to appear [2]. Physical program of NA6/SHINE at SPS [3, 4] and MPD at NICA [5] is focused on such investigations. Correlations and fluctuations are also studied at LHC energy [6 8] because they are sensitive to the very first stages of a collision [9]. So-called long-range correlations [0] between observables in separated rapidity intervals play a special role. Causality implies that they are developed in the initial state. In addition, a large enough rapidity gap can suppress the background contribution of resonance decays. In the present work, the string-parton Monte Carlo model [ 3] is applied to the relativistic nuclear collisions at high energy for calculation of long-range correlations and fluctuations. The model includes the string fusion [4 7] effects and the finite rapidity length of strings, and also charge conservation in string fragmentation. We calculated several types of strongly intensive event-by-event variables and long-range correlations for the experiment. 2 Model description The present model describes nucleon-nucleon collisions on the partonic level. For the nuclear density of Pb and Xe the Woods-Saxon distribution is taken. Each nucleon contain a valence quark-diquark pair and a certain number of sea quark-antiquark pairs, distributed around the center of nucleon according to two-dimensional Gauss distribution. The number of sea pairs e-mail: v.kovalenko@spbu.ru
is distributed according to Poisson law. An elementary interaction is realized in the model of colour dipoles [8, 9]. Then it is assumed, that if there is a collision between two dipoles, two quark-gluon strings are stretched between the ends of the dipoles, and the process of string fragmentation gives observable particles. The particle production of a string is assumed to go uniformly between the string rapidity ends y min and y max [2, 3], with the mean number of charged particles per rapidity µ 0, and independently in each rapidity interval, with Poisson distribution. Each parton can interact with another one only once, forming a pair of quarkgluon strings, hence, producing particles. The transverse position of a string is assigned to the arithmetic mean of the transverse coordinates of the partons at the ends of the string. Due to the finite transverse size of the strings they overlap, that in the framework of string fusion model [4 7] gives a source with higher tension. Mean multiplicity of charged particles and mean p t originated from the cell where k strings are overlapping are the following: µ k = µ 0 k, pt k = p 0 4 k. Here µ 0 and p 0 are mean charged multiplicity from one single string per rapidity unit and mean transverse momentum from one single string. For the overlapping of the strings of a finite rapidity length, the rapidity space is divided into regions where the number of overlapped strings is an integer, and they are processed separately. All parameters of the model are constrained from the data on inelastic cross section and multiplicity in pp, p-pb and Pb-Pb collisions [3]. The charge differentiation of the particles in this model is performed in the following way. If we neglect the production of double charge at string fragmentation, each emitted particle can change the charge going along rapidity only by 0 or (emitting of neutral or charged pion). After throwing away the neutral particles, we see that the charge of particles emitted from one string (or a cluster of fused strings) must alternate between to - along the rapidity axis. This procedure is implemented for each string cluster. 3 Definition of observables Firstly, we define, correlation coefficient of long range correlations between observables in two separated rapidity windows as follows: b = BF B F F 2 F 2, () where each of (B, F) could be either N ch the number of charged particles in the rapidity window or p t the mean event transverse momentum of charged particles in the given window: p t = N ch p ti. (2) N ch Correspondingly, three types of correlation coefficients can be studied n-n, pt-n, pt-pt correlations. Unfortunately, n-n and pt-n correlation coefficients are shown to be very dependent on the centrality class width and details of the centrality determination method [20] (note pt-pt has no this property). To overcome such issue, so-called, strongly intensive variables can be used [2 24]. In this paper we use only Σ[F, B] defined as: i= Σ(F, B) F ω B + B ω F 2 ( FB F B ) F + B, (3)
hb2 i hb2 i hf 2 i hf 2 i, ωb = scaled (or relative) variance. hfi hbi Variables F and B here must be intensive. For them we used total multiplicity in a rapidity window (N), and total transverse momentum (PT ). For the matter of study the charge fluctuations, we used dynamical fluctuation measure, νdyn, defined as *!2 +! N+ N. (4) νdyn = + hn+ i hn i hn+ i hn i It can be shown that the variable Nch νdyn = hn+ + N iνdyn is also strongly intensive. where ωf 4 Results Figure shows the centrality dependence of the mean transverse momentum correlation coefficient. We calculated it for two colliding systems Pb-Pb and Xe-Xe at the LHC energy to demonstrate also the system size dependence of pt-pt correlation coefficient. In case of Pb-Pb collisions, the model predicts [20] non-monotonic behaviour with centrality, which was also confirmed by experimental data [25], but it is hard to reproduce in the models [26]. The position of the peak can be attributed to the point of the saturation of the color field. In the case of the smaller system, the saturation region seems not to be reached because even in most central Xe+Xe collisions the required density is not achieved. Similar behavior also was obtained in p-pb collisions [27]. Figure : Mean transverse momentum correlation coefficient as a function of centrality at LHC energy, calculated in the Monte Carlo model with string fusion for Pb-Pb and Xe-Xe nuclei. In figure 2 the forward-backward strongly intensive measure Σ is shown for sufficiently large separated rapidity intervals (-3.8, -3), (3, 3.8). We see that above the very peripheral collisions the value of Σ decreases with centrality. Such behavior has been recently observed in the experimental data; however other models cannot reproduce such effect [28, 29]. Figure 3 shows the centrality dependence of the two other types of the strongly intensive quantities between observables in separated rapidity intervals, namely, Σ[PT, N] and Σ[PT, PT ]. Σ[PT, N] is rather flat, what reminds one the dependence of the average transverse momentum with multiplicity [30]. In pp and p-pb collisions the latter is less trivial so it worth to check the multiplicity dependence of a strongly intensive pt-n measure there.
,2,5,,05 0,95 0,9 0 0 20 30 40 50 60 70 80 90 00 Centrality, % Figure 2: Strongly intensive quantity Σ NN as a function of centrality, calculated in the Monte Carlo model with string fusion for Pb-Pb collisions at 2.76 TeV. Two centrality class binning (2% and 5%) are shown by different markers. 5% 2% 7 0,6 6 0,5 5 0,4 4 3 0,3 2 5% 0,2 5% 2% 0, 2% 0 0 0 20 30 40 50 60 70 80 90 00 Centrality, % 0 0 0 20 30 40 50 60 70 80 90 00 Centrality, % Figure 3: Strongly intensive quantities Σ PT N (left) and Σ PT P T (right) as a function of centrality, calculated in the Monte Carlo model with string fusion for Pb-Pb collisions at 2.76 TeV. As in fig. 2, two centrality class binning (2% and 5%) are shown by different markers. The strongly intensive transverse momentum variable Σ[P T, P T ] (right plot of figure 3) shows particular non-monotonic behavior. This reminds us the behavior of b pt p t. However, the position of the minimum here is at a different point than the maximum of b pt p t (here it lies in the range of less central collision). To summarize, all three variants of the forward-backward strongly intensive measure indeed doesn t depend on centrality; they provide additional non-trivial information which can not be extracted only from correlation coefficient and thus should be widely experimentally studied, at LHC energy and also at lower such as MPD at NICA. The figure 4 shows the charge fluctuations in terms of ν dyn (left plot) and N ch ν dyn (right plot) with comparison to the experimental data [6]. The results show that taking into account the string fusion improves the agreement with the experimental data. Without fusion, the model behaves as hadron resonance gas, wheres an adding of a string collectivity modifies N ch ν dyn and makes it centrality-dependent.
Figure 4: Dynamical charge fluctuation ν dyn as a function of centrality (number of participant nucleons) for Pb-Pb collisions at 2.76 TeV, calculated in the Monte Carlo model with string fusion and compared with experimental data [6]. The solid line corresponds to the model with string fusion, dashed without fusion. The right panel shows ν dyn multiplied by multiplicity. Areas that correspond to hadron resonance gas (HRG) and quark-gluon plasma (QGP) are hatched. 5 Conclusions Various strongly intense quantities are studied in AA collisions at the LHC energy, depending on centrality compared (where possible) with experimental data of ALICE. All analyzed strongly intensive variables and correlations between intensive observables show robust behavior against the volume fluctuations and the details of centrality determination. It is shown that the string fusion model correctly qualitatively describes the Σ NN dependence on the centrality. Charge fluctuations in Pb-Pb collisions at an energy of 2.76 TeV have been investigated. The correct description of the experimental ALICE data for ν dyn was obtained depending on the centrality of the collisions (string fusion effects improves the agreement with the experiment). String fusion model explains the non-monotonic behavior of pt-pt forward-backward correlations in Pb-Pb collisions and predicts that at smaller system this should not be observed. The studied observables provide a reference to the properties of the initial state of AA collisions and are useful to study both at LHC and NICA energies. 6 Acknowledgments The research was supported by the grant of the Russian Foundation for Basic Research (project 8-32-0055 mol_a). References [] L. Jiang, P. Li, H. Song, Phys. Rev. C 94, 02498 (206), arxiv: 52.0664 [nucl-th] [2] X. Luo, N. Xu, Nucl. Sci. Tech. 28, 2 (207), arxiv: 70.0205 [nucl-ex] [3] M. Gazdzicki, Eur. Phys. J. CST 55, 37 (2008), arxiv: 080.499 [nucl-ex] [4] M. Gazdzicki (NA49, NA6/SHINE Collaborations), J. Phys. G 38, 24024 (20), arxiv: 07.2345 [nucl-ex]
[5] V. Toneev, PoS CPOD07, 057 (2007), arxiv: 0709.459 [nucl-ex] [6] B. Abelev et al. (ALICE Collaboration), Phys. Rev. Lett. 0, 5230 (203), arxiv: 207.6068 [nucl-ex] [7] G. Aad et al. (ATLAS Collaboration), JHEP 207, 09 (202), arxiv: 203.300 [nucl-ex] [8] B. Abelev et al. (STAR Collaboration), Phys. Rev. Lett. 03, 7230 (2009), arxiv: 0905.0237 [nucl-ex] [9] P. Tribedy, J. Phys. Conf. Ser. 832, 0203 (207) [0] I. Altsybeev, J. Phys. Conf. Ser. 798, 02056 (207), arxiv: 6.0090 [nucl-ex] [] V. Kovalenko, V. Vechernin, PoS ISHEPP202, 077 (202), arxiv: 22.2590 [2] V. Kovalenko, Phys. Atom. Nucl. 76, 89 (203), arxiv: 2.6209 [hep-ph] [3] V. Kovalenko, PoS QFTHEP203, 052 (203) [4] M. Braun, C. Pajares, Phys. Lett. B 287, 54 (992) [5] M. Braun, C. Pajares, Nucl. Phys. B 390, 542 (993) [6] N. Amelin, M. Braun, C. Pajares, Phys. Lett. B 306, 32 (993) [7] M. Braun, C. Pajares, J. Ranft, Int. J. Mod. Phys. A 4, 2689 (999), arxiv: hep-ph/9707363 [8] C. Flensburg, G. Gustafson, L. Lonnblad, Eur. Phys. J. C 60, 233 (2009), arxiv: 0807.0325 [hep-ph] [9] G. Gustafson, Acta Phys. Polon. B 40, 98 (2009), arxiv: 0905.2492 [hep-ph] [20] V. Kovalenko, V. Vechernin, EPJ Web Conf. 66, 0405 (204), arxiv: 308.668 [nucl-th] [2] M.I. Gorenstein, M. Gazdzicki, Phys. Rev. C 84, 04904 (20), arxiv: 0.4865 [nucl-th] [22] M. Gazdzicki, M. Gorenstein, M. Mackowiak-Pawlowska, Phys. Rev. C 88, 024907 (203), arxiv: 303.087 [nucl-th] [23] E.V. Andronov, Theor. Math. Phys. 85, 383 (205) [Teor. Mat. Fiz. 85, no., 28 (205)] [24] E. Andronov, V. Vechernin, arxiv: 808.09770 [hep-ph] (208) [25] I. Altsybeev, KnE Energ. Phys. 3, 304 (208), arxiv: 7.04844 [nucl-ex] [26] V. Kovalenko, V. Vechernin, J. Phys. Conf. Ser. 798, 02053 (207), arxiv: 6.07274 [nucl-th] [27] V. Kovalenko, AIP Conf. Proc. 606, 74 (204), arxiv: 402.028 [nucl-th] [28] I. A. Sputowska (For ALICE Collaboration), Forward-backward correlations and multiplicity fluctuations in Pb-Pb collisions at s NN = 2.76 TeV from ALICE at the LHC (Hot Quarks 208. https://indico.cern.ch/event/70305/contributions/309537/, 208) [29] I. Sputowska, Correlations in Particle Production in Nuclear Collisions at LHC Energies. Ph.D. thises. (Cracow, INP, 206) [30] B.B. Abelev et al. (ALICE Collaboration), Phys. Lett. B 727, 37 (203), arxiv: 307.094 [nucl-ex]