KnE Energy | The 3rd International Conference on Particle Physics and Astrophysics (ICPPA) | pages: 226–233

1. Introduction

The NA61/SHINE experiment [1] is a multi-purpose fixed target experiment at the Super Proton Synchrotron (SPS) of the European Organization for Nuclear Research (CERN). The strong interactions programme of NA61/SHINE consists of studies of the onset of deconfinement (OD) [2] in nucleus-nucleus collisions and search for the critical point (CP) [3] of strongly interacting matter.

NA61/SHINE performs measurements of hadron production in collisions of protons and various nuclei (p+p, p+Pb, Be+Be, Ar+Sc, Xe+La, Pb+Pb) in a range of beam momenta (13A - 150/158A GeV/c). It is expected that there will be a non-monotonic dependence of fluctuations of a number of observables on energy and system size in this scan due to the phase transition of strongly interacting matter and the possible existence of the CP [4]. Some hints of such behaviour were already observed by the NA49 experiment [5].

2. PT-N fluctuation measures

The volume of the system created in ultrarelativistic ion collisions fluctuates from event to event. In order to suppress contribution from this "trivial" fluctuations strongly intensive observables are used [6,7]:


where stands for averaging over all events and ω[X] is the scaled variance of any quantity X defined as (X2-X2)/X . In case of joint fluctuations of the total transverse momentum PT and multiplicity N one defines: A=PT=i=1NpTi , B=N , CΔ=CΣ=Nω[pT] , where ω[pT] is the scaled variance of the inclusive pT spectrum.

The advantage of the strongly intensive quantities is that in the models of independent sources [8] and in the models of independent particle production the contribution from volume fluctuations is eliminated, allowing to probe the genuine CP signals. Moreover, in the models of independent particle production Δ[PT,N]=Σ[PT,N]=1 . These equalities also hold for the ideal Boltzmann gas both in the grand canonical and in the canonical ensemble formulations. Another trivial property of these quantities is that they vanish in case of absence of fluctuations of A and B .

In recent measurements by the NA61/SHINE collaboration [9,10] of the strongly intensive quantities Δ[PT,N] and Σ[PT,N] no anomaly attributable to the CP was observed neither in p+p nor in forward energy selected Be+Be and Ar+Sc collisions. This observation led to the extension of the analysis to studies of the pseudorapidity dependence of strongly intensive observables. By changing the width of the chosen pseudorapidity interval one can probe different values of baryochemical potential as evident from the different shapes of proton and antiproton rapidity spectra [11]. Therefore, the new analysis allows to probe additional regions of the phase diagram of strongly interacting matter.

In this paper new results were obtained for Be+Be collisions at 150A GeV/c selected for the smallest 8 % of forward energies. This selection of the most central events by the forward energy was done using information from the NA61/SHINE forward calorimeter PSD. This detector consists of 44 square modules. A specific set of modules is selected in order to maximize the signal from projectile nucleon spectators and to suppress background from produced particles. Analysis was performed for the 120k selected events containing 2000k quality selected tracks in total.

Fluctuations were studied in 9 pseudorapidity intervals defined in the laboratory reference frame using all charged particles produced in strong and electromagnetic processes from the primary vertex. The forward edge of the intervals was fixed at 5.2 units of pseudorapidity, with the backward edge changing from 3 up to 4.6 units. The choice of lower bounds was motivated by the small azimuthal angle acceptance at smaller pseudorapidities. The upper bound was introduced in order to suppress possible quasi-diffractive or electromagnetic effects which become important at larger pseudorapidities. Transverse momenta of all charged particles were restricted to 0<pT<1.5 GeV/c. Moreover, the NA61/SHINE acceptance map [12] was applied.

Figure 1 shows preliminary results on the dependence on the width of the pseudorapidity interval for Δ[PT,N] and Σ[PT,N] and comparisons to the EPOS1.99 model [13] predictions. Both quantities change monotonously for the data in contrast to the EPOS1.99 results for Δ[PT,N] which show a minimum for intermediate width and lie significantly above the measurements. EPOS1.99 also fails to describe Δ[PT,N] for the full NA61/SHINE pseudorapidity acceptance in inelastic p+p interactions [9]. Another observation is that for smaller windows these quantities approach unity (independent particle production limit) as the number of particles in the interval gets small. Moreover, the inequalities Δ[PT,N]<1 and Σ[PT,N]1 , previously observed for all studied systems at all collision energies [9,10], also hold when modifying the width of the pseudorapidity window of the measurement. In general, no traces of the possible critical point of strongly interacting matter are visible.

Figure 1

Dependence of Δ[PT,N] (a) and Σ[PT,N] (b) on the width of the pseudorapidity window. Red squares are preliminary NA61/SHINE measurements for 0-8 % Be+Be collisions at 150A GeV/c. Blue band represents the EPOS1.99 model predictions.

fig-1 fig-2
(a) a, (b) a,

Figure 2

Dependence of Σ[NF,NB] (a), Σ[NF,PTB] (b) and Σ[PTF,PTB] (c) on the separation between the two pseudorapidity windows. Red squares are preliminary NA61/SHINE measurements for 0-8 % Be+Be collisions at 150A GeV/c. Blue band represents the EPOS1.99 model predictions.


(a) a,


(b) a,


(c) a,

The results were not corrected for detector inefficiencies and trigger biases as simulations have shown that their effect estimated using the GEANT3 package does not exceed 5 % . In contrast the trigger bias modifies results significantly for p+p collisions [9].

The statistical uncertainties were determined using the sub-sample method. The systematic uncertainties were estimated to be smaller than 10 % for Δ[PT,N] and Σ[PT,N] . Further analysis of the systematic uncertainties will follow.

3. N-N , PT-N and PT-PT fluctuations in two rapidity intervals

By analogy with the PT-N case described in the previous section one can introduce a number of other strongly intensive observables. Joint fluctuations of two quantities in two separated pseudorapidity intervals are of special interest because they are closely connected with studies of forward-backward correlations that have a long history of measurements [14,15]. Forward-backward correlations are usually quantified by the correlation coefficient which is not strongly intensive and is sensitive to the employed centrality selection procedures in nucleus-nucleus collisions [16]. Use of strongly intensive quantities can suppress these trivial fluctuations and allows to study the intrinsic properties of particle emitting sources. Therefore, for studies of joint fluctuations of multiplicities ( NF , NB ) and the sum of transverse momenta ( PTF , PTB ) in two windows of rapidity the following quantities are defined [17]:






(5) Here, averaging over all tracks in the forward (backward) window is denoted as:


(6) Moreover, the scaled variance of inclusive pT spectra ω[pT]F,B is introduced for both forward and backward windows.

Typically, forward-backward correlations are studied not for the sum of transverse momenta but for the event mean transverse momentum [18], but in order to construct a strongly intensive quantity from two observables A and B both of these quantities have to be extensive [6].

The analysis was performed for 7 pairs of pseudorapidity intervals, with forward window fixed at (4.7,5.2) and the backward window moving from (3,3.5) to (4.2,4.7). Figure 2 shows preliminary results for the dependence of quantities (3)-(5) on the separation between the two windows. It is peculiar that all three quantities show qualitatively the same behaviour which is quite well reproduced by the EPOS1.99 model. Moreover, analogous behaviour was predicted by the quark-gluon string model [19] for p+p interactions at LHC energies. The PHOBOS collaboration studied joint fluctuations of multiplicities in two windows quantified by the variance of C=NF-NBNF+NB [20] which has similar properties like Σ[NF,NB] , although it is not strongly intensive. The variance of C showed a similar trend with separation between rapidity windows in Au+Au collisions at sNN=200 GeV.

Similarly to the results in the previous section, no corrections for detector inefficiencies and trigger biases were applied as simulations have shown that their effect does not exceed 5 % .

The statistical uncertainties were estimated using the sub-sample method, and the systematic uncertainties were estimated to be smaller than 10 % .

4. Conclusions

Dependence of the quantities Δ[PT,N] and Σ[PT,N] on the width of the pseudorapidity interval were measured by the NA61/SHINE experiment in forward energy selected Be+Be collisions at 150A GeV/c. No indications of the critical point of strongly interacting matter were observed. However, a significant qualitative discrepancy between the experimental results on Δ[PT,N] and the EPOS1.99 model predictions is clearly visible.

NA61/SHINE preliminary results for Σ[NF,NB] , Σ[NF,PTB] and Σ[PTF,PTB] on the separation of the forward and backward pseudorapidity intervals were also presented. It is the first measurement of this type at SPS energies. Experimental results are in good agreement with the EPOS1.99 model predictions. The observed trend to increase with separation is also similar to the predictions of the quark-gluon string model for LHC energies.


This work was supported by the Russian Science Foundation under grant 17-72-20045.



Abgrall N et al. [NA61 Coll.] 2014 JINST 9 06005


Gazdzicki M and Gorenstein M 1999, Acta Phys. Polon. B 30 2705


Fodor Z and Katz S 2004 J. High Energy Phys. JHEP0404(2004)050


M. Stephanov K. Rajagopal E. Shuryak Event-by-event fluctuations in heavy ion collisions and the QCD critical pointPhysical Review D: Particles, Fields, Gravitation and Cosmology1999601110.1103/physrevd.60.114028114028


K. Eskola H. Paukkunena C. Salgado EPS09—global NLO analysis of nuclear PDFs and their uncertaintiesNuclear Physics A20098301–4599c602c10.1016/j.nuclphysa.2009.10.054


M. I. Gorenstein M. Gaździcki Strongly intensive quantitiesPhysical Review C: Nuclear Physics201184110.1103/PhysRevC.84.014904


M. Gazdzicki M. I. Gorenstein M. Mackowiak-Pawlowska Normalization of strongly intensive quantitiesPhysical Review C nuclear physics20138822-s2.0-8488422954210.1103/PhysRevC.88.024907024907


A. Biallas M. Bleszynski W. Czyz Multiplicity distributions in nucleus-nucleus collisions at high energiesNuclear Physics B197611146147610.1016/0550-3213(76)90329-1


A. Aduszkiewicz Acta Physica Polonica B201243463510.5506/APhysPolB.43.635


E. Andronov Transverse momentum and multiplicity fluctuations in Ar+Sc collisions at the CERN SPS from NA61/shineActa Physica Polonica B20171034494532-s2.0-8503754224310.5506/APhysPolBSupp.10.449


Aduszkiewicz A et al [NA61 Coll.] 2017 Eur. Phys. J. C 77 no. 10 p 671


Czopowicz T Acceptance map used in paper


Pierog T and Werner K 2009 Nucl. Phys. Proc. Suppl. 196 102-105


Derrick M et al. 1986 Phys. Rev. D 34 3304


Adam J et al. [ALICE Coll.] 2015 J. High Energy Phys. JHEP1505(2015)097


Abelev B et al. [STAR Coll.] 2009 Phys. Rev. Lett. 103 172301


Andronov E 2015 Theor. Math. Phys. 185 no.1 p 1383


M. A. Braun R. S. Kolevatov C. Pajares V. V. Vechernin Correlations between multiplicities and average transverse momentum in the percolating color strings approachThe European Physical Journal C20043245355462-s2.0-244249307610.1140/epjc/s2003-01443-6


Vechernin V and Andronov E 2017 Strongly intensive observable for multiplicities in forward and backward windows in string model (XII Workshop on Particle Correlations and Femtoscopy, Amsterdam, 12-16 June 2017)


Back B et al. [PHOBOS Coll.] 2006 Phys. Rev. C 74 011901



  • Downloads 15
  • Views 67



ISSN: 2413-5453