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

1. Introduction

Nowadays numerous experimental and theoretical investigations of high energy nucleus-nucleus collisions show that the quark-gluon plasma could exist in nature. Moreover, the results of CERN LHC [1] and BNL RHIC [2] experiments and the observation of the transition between hadronic matter and quark-gluon plasma at CERN SPS energies [3], [4], [5] revealed that the key question in nuclear and particle physics now is to determine the structure of the phase transition between the hadron gas and the deconfined matter. It is expected that in the phase diagram of strongly interacting matter the hadron gas and quark-gluon plasma regions are separated by a first order phase transition line at high baryo-chemical potentials and moderate temperatures. A crossover between both phases is assumed for high temperatures and low baryo-chemical potentials. The first order phase transition line then ends in a critical point. But the exact location of the critical end-point in the phase diagram is unknown. Moreover, some lattice QCD calculations suggest that there might be no critical point at all and only a crossover separates the two phases.

The strong interactions programme of NA61/SHINE, a fixed-target experiment at the CERN SPS [6], includes the study of the properties of deconfinement and the search for the critical behaviour of strongly interacting matter. The main strategy of the NA61/SHINE collaboration in the search for the critical point is to perform a comprehensive two-dimensional scan of the phase diagram of strongly interacting matter by changing the collision energy and the system size [7]. If the critical point exists it is expected that at some values of these parameters a region of increased fluctuations should be observed. At the top of this “shark fin” hill the value of the critical fluctuations is expected to be a maximum. However, the critical signal could be shadowed since results on fluctuations are sensitive to conservation laws, resonance decays, volume fluctuations in the system of the colliding nuclei, quantum statistical effects and the limited acceptance of the experiment. Hence one should try to reduce the contribution from trivial fluctuations. This led to the idea to use intensive and strongly intensive quantities as probes of the critical behaviour [8], [9].

2. Quantities of interest

In order to make proper comparison of the results from different colliding systems, one should choose so-called intensive variables which are independent of the system size. Since in the vicinity of a critical point central second moments of distributions of extensive event quantity A are believed to diverge [10], the scaled variance ω[A] , an intensive quantity, was chosen for the analysis. The normalization results in ω[A] = 0 in the absence fluctuations of A and ω[A] = 1 in the case of a Poisson distribution of A [11]. One should keep in mind that ω[A] is still sensitive to fluctuations of the volume.

Due to the imperfect centrality determination in A+A collisions, one should expect event-by-event volume fluctuations. Consequently, to eliminate the influence of usually poorly known distributions of the system volume, it was suggested to use strongly intensive quantities which are independent both of the volume and fluctuations of the volume within the statistical model of the ideal Boltzmann gas in the grand canonical ensemble formulation [9], [11]. Two families of strongly intensive variables Δ[A,B] and Σ[A,B] were suggested which are functions of two extensive event quantities A and B .

The normalization of these variables can be chosen such that [11]:

  • Δ[A,B]=Σ[A,B]=0 in the absence of fluctuations of A and B

  • Δ[A,B]=Σ[A,B]=1 in the Independent Particle Model

3. Analysis details

The main goal of this work is to extend the investigation of the phase diagram by measuring the pseudorapidity dependence of fluctuations. Analysis of proton-proton collisions is the baseline for future investigations of nucleus-nucleus collisions.

This paper presents results referring to all charged hadrons with pT<1.5 GeV/c produced in the analysis acceptance of the NA61/SHINE experiment [12] in inelastic proton-proton collisions at beam momenta of 20, 31, 40, 80 and 158 GeV/c. Only data-based corrections for off-target interactions were performed. Simulation-based corrections using the event generator EPOS1.99 [13] and the NA61/SHINE detector simulation chain are in progress for other biases.

Figure 1

Sketch of the uncorrected pseudorapidity distribution of charged hadrons with chosen windows (pseudorapidity intervals) a) for one window analysis, b) for two windows analysis





Figure 2

Dependence on the width of the rapidity window (see Fig.1a) of a) scaled variance ω[N] , b) strongly intensive quantity Δ[PT,N] in the one window analysis at SPS energies for all charged hadrons produced in inelastic p+p collisions in the NA61/SHINE acceptance.





Figure 3

Dependence on the width of the pseudorapidity window (see Fig.1a) of strongly intensive quantities a) Δ[PT,N] , b) Σ[PT,N] in the one window analysis for data and EPOS1.99 at beam momentum of 158 GeV/c for all charged hadrons produced in inelastic p+p collisions in the NA61/SHINE acceptance.





Figure 4

Dependence of the strongly intensive fluctuation measure Σ[NF,NB] for all charged hadrons produced in inelastic p+p collisions in the NA61/SHINE acceptance at beam momenta of a) 20 GeV/c, b) 158 GeV/c on the distance between the Forward and the Backward pseudorapidity window in the two windows analysis (see Fig.1b)





The analysis has two parts. In the one window analysis the intensive quantity ω[N] and the strongly intensive quantities Δ[PT,N] and Σ[PT,N] are studied as functions of the width of the chosen pseudorapidity interval. This corresponds to changing the rapidity averaged baryo-chemical potential at the freeze-out stage [14]. The two windows analysis performs a search for short- and long-range correlations via study of the dependence of the strongly intensive quantity Σ[NF,NB] on the distance between two separated pseudorapidity windows in which the Forward and Backward multiplicities are evaluated.

The one window analysis

Here A was taken as the event multiplicity of charged hadrons N and B as the scalar sum PT of their transverse momenta. With denoting the average value over all events and ¯ the inclusive average value (over all particles and all events), one can define the following intensive quantities [8]:


and strongly intensive quantities [8]:


with normalization [11]:


The two windows analysis

Taking the extensive event quantities A as the multiplicity NF in the Forward pseudorapidity window and B as the multiplicity NB in the Backward pseudorapidity window, the above formulae will transform [15] into:



Definitions of pseudorapidity intervals

All results are presented as functions of Δη/Δηmax in the lab system according to the Fig.1. In the following the lower edges of the pseudorapidity intervals are moving from ybeamlab/2 to ybeamlab corresponding to the range from ycms=0 to ybeamcms . ηlab is restricted to this range to exclude the influence of the bad acceptance coverage at small ηlab and to reduce contributions from elastic processes at ηlab>ybeamlab .

4. Results

The presented preliminary results refer to all charged hadrons with pT<1.5 GeV/c produced in the acceptance of the NA61/SHINE experiment [12] in inelastic p+p collisions.

The studied fluctuation measures significantly depend on the width and the location of the pseudorapidity intervals. Results for ω[N] and Δ[PT,N] depend on the collision energy (see Fig.2). The values of ω[N] rise with increasing beam energy and width of the rapidity interval. The values of Δ[PT,N] are less than one with deviation increasing with beam energy and width of the rapidity interval. On the other hand the values of Σ[PT,N] are similar at all beam momenta and show the same rising tendency with increasing width of the rapidity interval (see conference slides). The azimuthal acceptance changes only slightly for different beam momenta. A significant discrepancy between data and EPOS1.99 calculations is observed for Δ[PT,N] at all collision energies (see example for 158 GeV/c presented in Fig. 3, left). It is more pronounced for the larger width of the pseudorapidity interval. However, EPOS1.99 well describes ω[N] and Σ[PT,N] (see Fig. 3, right and conference slides). In the two windows analysis there is an increase of the value of Σ[NF,NB] with the distance between forward and backward pseudorapidity intervals which is more pronounced for the higher energy (see Fig. 4). EPOS1.99 predictions are in a good agreement with the data.


This work was supported by the Saint-Petersburg State University research grants, 11.41.1385.2017.



Schukraft J et al [for ALICE Collaboration] 2011 J. Phys. G 38 124003


Adams J et al [for STAR Collaboration] 2005 Nucl. Phys. A 757 102


Alt C et al [NA49 Collaboration] 2008 Phys. Rev. C 77 024903


M. Gazdzicki M. Gorenstein P. Seyboth Onset of deconfinement in nucleus-nucleus collisions: review for pedestrians and expertsActa Physica Polonica B201142230735110.5506/aphyspolb.42.3072-s2.0-79952027646


Gazdzicki M, Gorenstein M I and Seyboth P 2014 Int. J. Mod. Phys. E 23 1430008


Abgrall N et al [NA61/SHINE Collaboration] 2014 JINST 9 P06005


M. Gazdzicki P. Seyboth Search for critical behaviour of strongly interacting matter at the cern super proton synchrotronActa Physica Polonica B2016475120112362-s2.0-8500847820810.5506/aphyspolb.47.120110.5506/APhysPolB.47.1201


Aduszkiewicz A et al [NA61/SHINE Collaboration] 2016 Eur. Phys. J. C 76 635


Gorenstein M and Gazdzicki M 2011 Phys. Rev. C 84 014904


V. Vovchenko D. V. Anchishkin M. I. Gorenstein Particle number fluctuations for the van der Waals equation of stateJournal of Physics A: Mathematical and General20154830305001, 1610.1088/1751-8113/48/30/305001MR3367963


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


Czopowicz T Acceptance map used in paper


K. Werner F.-M. Liu T. Pierog Parton ladder splitting and the rapidity dependence of transverse momentum spectra in deuteron-gold collisions at the BNL Relativistic Heavy Ion ColliderPhysical Review C: Nuclear Physics20067441110.1103/PhysRevC.74.044902044902


F. Becattini J. Manninen M. Gaździcki Energy and system size dependence of chemical freeze-out in relativistic nuclear collisionsPhysical Review C: Nuclear Physics200673410.1103/PhysRevC.73.044905


Andronov E V 2015 Theor. Math. Phys. 185, no. 1 1383 1390



  • Downloads 15
  • Views 96



ISSN: 2413-5453