## Centrality Dependence of Multiplicity Fluctuations in Ion-Ion Collisions from the Beam Energy Scan at FAIR

Anuj Chandra, Bushra Ali and Shakeel Ahmad 1

Department of Physics, Aligarh Muslim University Aligarh INDIA

## Abstract

Multiplicity distributions and event-by-event multiplicity fluctuations in AuAu collisions at energies in future heavy-ion experiment at the Facility for Anti-proton and Ion Research ( FAIR ) are investigated. Events corresponding to FAIR energies are simulated in the frame work of Ultra Relativistic Quantum Molecular Dynamics ( URQMD ) model. It is observed that the mean and the width of multiplicity distributions monotonically increase with beam energy. The trend of variations of dispersion with mean number of participating nucleons for the centrality-bin width of 5% are in accord with the Central Limit Theorem. The multiplicity distributions in various centrality bins as well as for full event samples are observed to obey Koba, Nielsen and Olesen ( KNO ) scaling. The trends of variations of scaled variance with beam energy are also found to support the KNO scaling predictions for larger collision centrality. The findings also reveal that the statistical fluctuations in 5% centrality-bin width appear to be under control.

PACS numbers:

25.75.q, 25.75.Gz

KEY-WORDS: Beam Energy Scan, Event-by-event fluctuations, Relativistic heavy-ion collisions.

1 email: Shakeel.Ahmad@cern.ch Authors declare that there is no conflict of interest

## 1. Introduction

Any physical quantity measured in an experiment is subject to fluctuations. These fluctuations depend on the property of the system and are expected to provide important information about the nature of the system under study [1,2]. As regards relativistic heavy-ion (AA) collisions, the system so created is a dense and hot fireball consisting of partonic and (or) hadronic matter [1,2]. To investigate the existence of partonic matter in the early life of fireball is one of the main goals of AA collisions. Study of fluctuations in AA collisions would help check the idea that fluctuations of a thermal system are directly related to various susceptibilities and could be an indicator for the possible phase transitions [1,2,3]. Fluctuations in experimental observables, such as charged particle multiplicity, particle ratios, mean transverse momentum and other global observables are related to the thermodynamic properties of the system, like, entropy, specific heat, chemical potential, etc. [4,5,6,7]. Event-by-event (ebe) fluctuations of these quantities are regarded as an important mean to understand the particle production dynamics which, in turn, would lead to understand the nature of phase transition and the critical fluctuations at the QCD phase boundary. A non-monotonic behavior of the fluctuations as a function of collision centrality and energy of the colliding beam may signal the onset of confinement and may be used to probe the critical point in the QCD phase diagram [7]. The multiplicity of charged particles produced in heavy-ion collisions is the simplest and day-one observable which provides a mean to investigate the dynamics of highly excited multi-hadron system. Studies involving multiplicity distributions (MDs) of the relativistic charged particles produced would allow finding the deviations from a simple superposition of multiple independent nucleon-nucleon ( nn ) collisions. Such studies, if carried out in limited rapidity space are envisaged to provide useful information on dynamical fluctuations [8,9,10,11]. It has been stressed that moments of MDs in full and limited rapidity bins would lead to make some interesting remarks about the production mechanisms involved. Dependence of MDs and their moments on collision centrality is also expected to lead to some interesting conclusions because of the fact that in narrow centrality windows the geometrical fluctuations may be treated as under control, whereas, such windows, if correspond to most central collisions, may be of additional importance because of the extreme conditions of temperature and excitation energy [7]. An attempt is, therefore, made to study the multiplicity fluctuations in the narrow centrality windows in AuAu collisions for the Beam Energy Scan ( BES ) at FAIR energies (for E lab = 10, 20, 30 and 40A GeV) in the frame work of URQMD model, using the code, urqmd-v3.4 [12, 13]. The number of events simulated at these energies are 2.3, 2.3, 2.1 and 2.2M (M = 10 ) respectively. 6 The analysis is carried out in the pseudorapidity and p T intervals with

-1 0 . &lt; η &lt; 1 0 and 0 2 . . &lt; p T &lt; . 5 0 GeV/c respectively.

## 2. The URQMD Model

Multiparticle production in relativistic collisions have been described earlier within the hydrodynamic approach [14]. At a later stage the Regge theory [15] and multiperipheral models were developed [15,16]. Although the difficulties attributed to the statistical models wereover come in these models yet the inconvenience of this approach is the large number of free parameters which are to be fixed by comparison with the experiments. Subsequently various quark-parton models motivated by QCD were introduced and as a consequence a large variety of models for hadronic and heavy-ion collisions were proposed. These models may be classified into macroscopic (statistical and thermodynamic) models [17] and microscopic (string, transport, cascade, etc.) models, like URQMD , VENUS , RQMD , etc. The microscopic models describe the individual hadron-hadron collisions.

URQMD model is based on the co-variant propagation of constituent quarks and di-quarks but

|   Centrality(%) | < N part >      | < N ch >      | σ            | ω               |
|-----------------|-----------------|---------------|--------------|-----------------|
|               5 | 348.00 ± 0.0020 | 231.03 ± 0.07 | 25.07 ± 0.05 | 2.7198 ± 0.0009 |
|              10 | 289.90 ± 0.0020 | 182.37 ± 0.06 | 21.63 ± 0.05 | 2.5650 ± 0.0009 |
|              15 | 238.45 ± 0.0020 | 144.36 ± 0.05 | 19.34 ± 0.05 | 2.5898 ± 0.0010 |
|              20 | 195.27 ± 0.0020 | 113.62 ± 0.04 | 17.26 ± 0.05 | 2.6214 ± 0.0012 |
|              25 | 159.21 ± 0.0020 | 89.70 ± 0.04  | 15.19 ± 0.05 | 2.5727 ± 0.0013 |
|              30 | 127.17 ± 0.0020 | 70.10 ± 0.04  | 13.85 ± 0.05 | 2.7352 ± 0.0016 |
|              35 | 100.08 ± 0.0020 | 53.49 ± 0.03  | 12.18 ± 0.05 | 2.7734 ± 0.0019 |
|              40 | 77.97 ± 0.0010  | 39.72 ± 0.03  | 10.76 ± 0.05 | 2.9171 ± 0.0022 |
|              45 | 58.89 ± 0.0010  | 28.99 ± 0.02  | 9.22 ± 0.05  | 2.9297 ± 0.0027 |
|              50 | 44.44 ± 0.0010  | 20.62 ± 0.02  | 7.84 ± 0.05  | 2.9803 ± 0.0032 |
|              55 | 31.99 ± 0.0010  | 13.86 ± 0.02  | 6.46 ± 0.05  | 3.0073 ± 0.0039 |
|              60 | 22.27 ± 0.0010  | 9.13 ± 0.01   | 5.17 ± 0.05  | 2.9287 ± 0.0049 |
|              65 | 14.83 ± 0.0010  | 5.95 ± 0.01   | 4.14 ± 0.05  | 2.8798 ± 0.0059 |
|              70 | 10.15 ± 0.0010  | 3.67 ± 0.01   | 3.21 ± 0.05  | 2.8056 ± 0.0071 |
|              75 | 7.06 ± 0.0020   | 2.17 ± 0.01   | 2.42 ± 0.05  | 2.6991 ± 0.0089 |
|              80 | 5.53 ± 0.0050   | 1.20 ± 0.00   | 1.77 ± 0.05  | 2.6062 ± 0.0111 |

TABLE I: Values of &lt; N part &gt; &lt; N , ch &gt; , dispersion ( σ ) and scaled variance ( ω ) in various centrality bins at E lab = 10A GeV/c

has been accompanied by baryonic and mesonic degrees of freedom. At low energies, √ s NN &lt; 5

TABLE II: Values of the same variables, as in Table I, but for E lab = 20A GeV/c

|   Centrality(%) | < N part >      | < N ch >      | σ            | ω               |
|-----------------|-----------------|---------------|--------------|-----------------|
|               5 | 348.00 ± 0.0020 | 288.81 ± 0.07 | 27.91 ± 0.05 | 2.6980 ± 0.0008 |
|              10 | 289.90 ± 0.0020 | 227.70 ± 0.06 | 24.58 ± 0.05 | 2.6535 ± 0.0008 |
|              15 | 238.45 ± 0.0020 | 179.76 ± 0.06 | 22.37 ± 0.04 | 2.7839 ± 0.0010 |
|              20 | 195.27 ± 0.0020 | 141.24 ± 0.05 | 20.27 ± 0.04 | 2.9100 ± 0.0012 |
|              25 | 159.21 ± 0.0020 | 111.38 ± 0.05 | 18.10 ± 0.03 | 2.9414 ± 0.0014 |
|              30 | 127.17 ± 0.0020 | 86.96 ± 0.04  | 16.52 ± 0.03 | 3.1384 ± 0.0017 |
|              35 | 100.08 ± 0.0020 | 66.34 ± 0.04  | 14.62 ± 0.03 | 3.2211 ± 0.0020 |
|              40 | 77.97 ± 0.0010  | 49.44 ± 0.03  | 12.90 ± 0.02 | 3.3637 ± 0.0024 |
|              45 | 58.89 ± 0.0010  | 36.03 ± 0.03  | 11.05 ± 0.02 | 3.3876 ± 0.0030 |
|              50 | 44.44 ± 0.0010  | 25.75 ± 0.02  | 9.49 ± 0.02  | 3.4943 ± 0.0035 |
|              55 | 31.99 ± 0.0010  | 17.36 ± 0.02  | 7.80 ± 0.01  | 3.5043 ± 0.0042 |
|              60 | 22.27 ± 0.0010  | 11.53 ± 0.02  | 6.32 ± 0.01  | 3.4638 ± 0.0055 |
|              65 | 14.83 ± 0.0010  | 7.55 ± 0.01   | 5.11 ± 0.01  | 3.4627 ± 0.0066 |
|              70 | 10.15 ± 0.0010  | 4.71 ± 0.01   | 4.00 ± 0.01  | 3.4008 ± 0.0080 |
|              75 | 7.06 ± 0.0020   | 2.79 ± 0.01   | 3.06 ± 0.01  | 3.3583 ± 0.0100 |
|              80 | 5.53 ± 0.0050   | 2.25 ± 0.21   | 2.62 ± 0.15  | 3.0484 ± 0.3079 |

GeV, the collisions are described in terms of interactions between hadrons and their excited states [17], whereas at higher energies ( &gt; 5 GeV), the quark and gluon degrees of freedom are considered and the concept of color string excitation is introduced with their subsequent fragmentation into hadrons [13]. In a transport model, AA collisions are considered as the superposition of all possible binary nn collisions. Every nn collision corresponding to the impact parameter, b ≤ √ σ tot /π is considered, where σ tot represents the total cross section. The two colliding nuclei are described by Fermi gas model [17] and hence the initial momentum of each nucleon is taken at random between zero and Thomas-Fermi momentum. The interaction term includes more than 50 baryon and 45 meson species. The model can treat the intermediate fireball both in and out of a local thermal and chemical equilibria. The URQMD model, thus, provides an ideal framework to study heavy-ion collisions. Although, the phase transition from a hadronic to partonic phase are not explicitly included in the model, thus a clear suggestion about the location of critical point cannot be made. The study, however, might help in the interpretation of the experimental data since it will permit subtraction of simple dynamical and geometrical effects from the expected Quark Gluon Plasma ( QGP ) signals [18].

TABLE III: Values of the same variables, as in Table I, but for E lab = 30A GeV/c

|   Centrality(%) | < N part >      | < N ch >      | σ            | ω               |
|-----------------|-----------------|---------------|--------------|-----------------|
|               5 | 348.00 ± 0.0020 | 327.23 ± 0.09 | 31.63 ± 0.06 | 3.0566 ± 0.0009 |
|              10 | 289.90 ± 0.0020 | 257.60 ± 0.08 | 27.73 ± 0.05 | 2.9857 ± 0.0010 |
|              15 | 238.45 ± 0.0020 | 203.53 ± 0.07 | 25.22 ± 0.05 | 3.1245 ± 0.0011 |
|              20 | 195.27 ± 0.0020 | 159.94 ± 0.06 | 22.70 ± 0.04 | 3.2205 ± 0.0013 |
|              25 | 159.21 ± 0.0020 | 126.11 ± 0.06 | 20.26 ± 0.04 | 3.2551 ± 0.0016 |
|              30 | 127.17 ± 0.0020 | 98.62 ± 0.05  | 18.56 ± 0.03 | 3.4930 ± 0.0019 |
|              35 | 100.08 ± 0.0020 | 75.28 ± 0.04  | 16.47 ± 0.03 | 3.6032 ± 0.0023 |
|              40 | 77.97 ± 0.0010  | 56.15 ± 0.04  | 14.48 ± 0.03 | 3.7316 ± 0.0027 |
|              45 | 58.89 ± 0.0010  | 40.95 ± 0.03  | 12.40 ± 0.02 | 3.7530 ± 0.0033 |
|              50 | 44.44 ± 0.0010  | 29.36 ± 0.03  | 10.63 ± 0.02 | 3.8484 ± 0.0040 |
|              55 | 31.99 ± 0.0010  | 19.87 ± 0.02  | 8.83 ± 0.02  | 3.9211 ± 0.0048 |
|              60 | 22.27 ± 0.0010  | 13.20 ± 0.02  | 7.19 ± 0.01  | 3.9103 ± 0.0063 |
|              65 | 14.83 ± 0.0010  | 8.68 ± 0.02   | 5.81 ± 0.01  | 3.8935 ± 0.0076 |
|              70 | 10.15 ± 0.0010  | 5.44 ± 0.01   | 4.60 ± 0.01  | 3.8949 ± 0.0094 |
|              75 | 7.06 ± 0.0020   | 3.23 ± 0.01   | 3.51 ± 0.01  | 3.8103 ± 0.0116 |
|              80 | 5.53 ± 0.0050   | 3.01 ± 0.25   | 3.13 ± 0.18  | 3.2661 ± 0.2933 |

## 3. Results and discussion

The URQMD model gives the value of impact parameter, b on ebe basis which allows to determine the collision centrality and mean number of participating nucleons, 〈 N part 〉 using the Glauber model [7,19]. Values of number of participating nucleons, mean charged particle multiplicities and dispersion of MDs ( σ ) for various collision centralities at the four energies are estimated and listed in Tables I - IV. The centrality selection is made from the MDs of charged particles for the minimum bias events in the considered η and p t ranges. This is illustrated in FIG.1, where the multiplicity distribution of minimum bias events for E lab = 40A GeV is displayed. The shaded regions show 10% centrality cross-section bins. Variations of 〈 N ch 〉 and σ with 〈 N part 〉 for the centrality bin width = 2, 5 and 10% are presented in FIGs.2 and 3. The statistical errors associated with these parameters are too small to be noticed in the figure. It may be noted from the figure that 〈 N ch 〉 and σ increase smoothly with 〈 N part 〉 or collision centrality. The lines in FIG.2 are due to the best fits to the data obtained using the equation

$$< N _ { c h } > = a + b < N _ { p a r t } > + c < N _ { p a r t } > ^ { 2 }$$

TABLE IV: Values of the same variables, as in Table I, but for E lab = 40A GeV/c

|   Centrality(%) | < N part >      | < N ch >      | σ            | ω               |
|-----------------|-----------------|---------------|--------------|-----------------|
|               5 | 348.00 ± 0.0020 | 353.86 ± 0.10 | 35.23 ± 0.07 | 3.5065 ± 0.0010 |
|              10 | 289.90 ± 0.0020 | 278.50 ± 0.08 | 30.55 ± 0.06 | 3.3503 ± 0.0011 |
|              15 | 238.45 ± 0.0020 | 219.93 ± 0.07 | 27.86 ± 0.05 | 3.5289 ± 0.0013 |
|              20 | 195.27 ± 0.0020 | 172.92 ± 0.07 | 24.98 ± 0.05 | 3.6082 ± 0.0015 |
|              25 | 159.21 ± 0.0020 | 136.43 ± 0.06 | 22.09 ± 0.04 | 3.5766 ± 0.0017 |
|              30 | 127.17 ± 0.0020 | 106.50 ± 0.05 | 20.16 ± 0.04 | 3.8159 ± 0.0020 |
|              35 | 100.08 ± 0.0020 | 81.49 ± 0.05  | 17.75 ± 0.03 | 3.8673 ± 0.0024 |
|              40 | 77.97 ± 0.0010  | 60.80 ± 0.04  | 15.78 ± 0.03 | 4.0966 ± 0.0029 |
|              45 | 58.89 ± 0.0010  | 44.31 ± 0.04  | 13.50 ± 0.03 | 4.1148 ± 0.0036 |
|              50 | 44.44 ± 0.0010  | 31.76 ± 0.03  | 11.55 ± 0.02 | 4.1974 ± 0.0042 |
|              55 | 31.99 ± 0.0010  | 21.58 ± 0.02  | 9.61 ± 0.02  | 4.2835 ± 0.0051 |
|              60 | 22.27 ± 0.0010  | 14.31 ± 0.02  | 7.78 ± 0.01  | 4.2327 ± 0.0066 |
|              65 | 14.83 ± 0.0010  | 9.43 ± 0.02   | 6.28 ± 0.01  | 4.1880 ± 0.0079 |
|              70 | 10.15 ± 0.0010  | 5.90 ± 0.01   | 4.95 ± 0.01  | 4.1634 ± 0.0097 |
|              75 | 7.06 ± 0.0020   | 3.50 ± 0.01   | 3.77 ± 0.01  | 4.0707 ± 0.0120 |
|              80 | 5.53 ± 0.0050   | 2.04 ± 0.03   | 2.97 ± 0.02  | 4.3207 ± 0.0620 |

whereas, in FIG.3 the lines are due to the least square fits to the data of the form

$$\sigma = p + q \sqrt { < N _ { p a r t } > }$$

The values of coefficients, occurring in Eqs.1 and 2 are listed in Tables V and VI respectively. As described in ref.7, the centrality dependence of the moments may be understood by the Central Limit Theorem (CLT), according to which, 〈 N ch 〉 ∝ N part and σ ∝ √ N part . However, in the present study the mean multiplicity is observed to grow with 〈 N part 〉 , as given by Eq.1., i.e. a slight deviation from linearity is exhibited by the data irrespective of the fact that how large or small the centrality bins are chosen. The variations of σ with 〈 N part 〉 , shown in FIG.3, is seen to be nicely fitted by Eq.2 for 5% centrality bin width, while for the centrality bin widths of 2% and 10% the data are seen to be fitted only for centrality &gt; 20%, as indicated by the lines in this figure; the lines are drawn for the range of centrality for which the fits of the data have been performed. Similar deviations from CLT have also been observed in AuAu collisions at RHIC and lower energies [7]. In order to extract dynamical fluctuations arising from physical processes, fluctuations in mean number of participating nucleons are to be minimized. To achieve the same, centrality bins

TABLE V: Values of parameters, a, b and c, occurring in Eq.1 at different energies.

| E lab   | Fit Par.   | Centrality 10%   | Centrality 5%   | Centrality 2%   |
|---------|------------|------------------|-----------------|-----------------|
| 10A GeV | a × 10 - 1 | -17.425 ± 0.032  | -19.211 ± 0.071 | -13.718 ± 0.004 |
| 10A GeV | b × 10 - 2 | 49.595 ± 0.018   | 49.581 ± 0.020  | 43.791 ± 0.015  |
| 10A GeV | c × 10 - 4 | 4.895 ± 0.008    | 4.452 ± 0.007   | 5.350 ± 0.005   |
| 20A GeV | a × 10 - 1 | -15.317 ± 0.051  | -18.215 ± 0.063 | -17.031 ± 0.045 |
| 20A GeV | b × 10 - 2 | 60.103 ± 0.024   | 60.250 ± 0.021  | 53.923 ± 0.017  |
| 20A GeV | c × 10 - 4 | 6.668 ± 0.010    | 6.670 ± 0.007   | 6.900 ± 0.006   |
| 30A GeV | a × 10 - 1 | -16.633 ± 0.061  | -18.670 ± 0.066 | -17.852 ± 0.050 |
| 30A GeV | b × 10 - 2 | 68.297 ± 0.029   | 68.483 ± 0.024  | 60.889 ± 0.019  |
| 30A GeV | c × 10 - 4 | 7.458 ± 0.012    | 7.428 ± 0.008   | 7.872 ± 0.007   |
| 40A GeV | a × 10 - 1 | -18.499 ± 0.063  | -19.733 ± 0.065 | -18.947 ± 0.049 |
| 40A GeV | b × 10 - 2 | 74.059 ± 0.030   | 73.900 ± 0.026  | 65.885 ± 0.020  |
| 40A GeV | c × 10 - 4 | 7.986 ± 0.013    | 8.080 ± 0.010   | 8.468 ± 0.007   |

TABLE VI: Values of parameters, p and q, occurring in Eq.2 at different energies.

| E lab   | Fit Par.   | Centrality 10%   | Centrality 5%   | Centrality 2%   |
|---------|------------|------------------|-----------------|-----------------|
| 10A GeV | p × 10 - 1 | -21.171 ± 0.061  | -10.181 ± 0.087 | 8.087 ± 0.055   |
| 10A GeV | q × 10 - 1 | 16.436 ± 0.017   | 13.244 ± 0.001  | 11.704 ± 0.011  |
| 20A GeV | p × 10 - 1 | -19.790 ± 0.082  | -8.694 ± 0.008  | -6.706 ± 0.067  |
| 20A GeV | q × 10 - 1 | -19.050 ± 0.019  | 15.308 ± 0.001  | 13.634 ± 0.013  |
| 30A GeV | p × 10 - 1 | -21.280 ± 0.097  | -8.795 ± 0.008  | -6.936 ± 0.073  |
| 30A GeV | q × 10 - 1 | 21.334 ± 0.023   | 17.131 ± 0.002  | 15.295 ± 0.015  |
| 40A GeV | p × 10 - 1 | -22.988 ± 0.101  | -11.275 ± 0.008 | -8.658 ± 0.075  |
| 40A GeV | q × 10 - 1 | 23.292 ± 0.026   | 18.901 ± 0.002  | 16.814 ± 0.016  |

considered should be kept narrow because the fluctuations in the particle multiplicities are directly related to the fluctuations in the mean number of participating nucleons. The inherent fluctuations may be reduced by choosing narrow centrality bins; the inherent fluctuations are the fluctuations which arise due to the difference in the geometry even within the selected centrality bin. A very narrow centrality bin, if considered, would, therefore, minimize this effect but may cause additional fluctuations due to statistics. Centrality resolution of the detectors also demands that the chosen centrality bins should not be too narrow. Thus, our observations from FIG.3, tend to suggest that fluctuations effects dominate if the centrality bin width is somewhat larger or quite small.

FIG. 1: An example of centrality selection from the multiplicity distribution of minimum bias simulated events at E lab = 40A GeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-14358a6a-c085-4c66-b335-4f3a51659615.jpg)

Multiplicity distributions of relativistic charged particles for minimum bias events for | η | &lt; 1 0 . and p T = 0 2 . -5 0 . GeV/c are displayed in FIG.4. It may be noted from the figure that MDs at the four beam energies considered, acquire nearly similar shapes and it is expected that the maximum values of N ch become higher with increasing energies. Similar trends in MDs have also been reported by S. Ghosh et al [17] at the same beam energies predicted by URQMD model. MDs of relativistic charged particles for various centrality groups at the four beam energies have also been examined. Distributions for E lab = 40A GeV are presented in FIG.5 along-with the distribution of full sample of events (minimum bias). It is evidently clear from the figure that MD of minimum bias sample is a convolution of MDs with different centrality classes.

Yet another way to examine and predict the MDs, is to plot MDs in terms of KNO scaling variable Z (= N ch / N 〈 ch 〉 ). It has been observed that MDs in hadron-hadron collisions exhibit a universal behavior in a wide range of incident energies if plotted as 〈 N ch 〉 P N ( ch ) against the variable Z [20,21,22,23,24,25]. It was shown that MDs corresponding to pp collisions in the energy range ∼ (50 - 303) GeV are nicely reproduced by the functional form given by Slattery [22]. MDs in pp collisions, for non single diffractive events at ISR energies have also been observed to exhibit KNO scaling [26]. Since the width of MDs for a given centrality gives the extent of fluctuations, the

FIG. 2: Variations of mean multiplicity with mean number of participating nucleons. The lines are due to fits obtained using Eq.1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-3771b818-ff57-4cca-8a9f-6e208010c767.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-d26383f7-5f18-49c2-9005-83d4778a6f61.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-2402d641-061f-4873-9422-9d420c748dd3.jpg)

origin of the fluctuations are, thus, inherent in the width of MDs. To understand this behavior, MDs should be plotted for different centrality bins in terms of KNO scaling variable. MDs for 10, 30 and 50% centrality are plotted in terms of KNO scaling variable in FIG.6. For clarity sake, each next distribution is shifted up on y-scale by a factor of 10. It is observed that the distributions become wider with increasing collision centrality, but exhibits a perfect scaling behavior. MDs, plotted in terms of KNO variable for full event sample in FIG.7, are also noticed to show a perfect KNO scaling.

The scaled variance, ω of the MDs defined as,

σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-24af7b54-4389-4f2b-9e6c-c1b12be867f3.jpg)

σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-d27479e6-13bb-421d-822e-0a85a1c7ac8c.jpg)

σ

FIG. 3: Variations of dispersion with mean number of participating nucleons. The lines are due to fits obtained using Eq.2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-08a11bbc-125f-4f54-abc4-b48df81fb20a.jpg)

$$\omega = \frac { \sigma ^ { 2 } } { \langle N _ { c h } \rangle } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

here ω is regarded as a quantitative measure of the particle number fluctuations [7,18,27,28,29]. The scaled variance, ω is an intensive quantity which does not depend on the volume of the system within the grand canonical ensemble (GCE) of statistical mechanics or on the number of sources within models of independent source, like wounded nucleon model. The value of scaled variance will be zero in the absence of fluctuations in MDs and unity for Poisson MDs. Since the volume of the system created in AA collisions fluctuates from event to event, and ω would depend on volume fluctuations, it becomes important to reduce the fluctuation effects in fluctuation studies

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-ea5b0fc4-533e-47a4-ad56-c7f4b6da8b29.jpg)

ch

FIG. 4: Multiplicity distributions of charged particles for AuAu collisions at 10A, 20A 30A and 40A GeV in the range p T = 0 2 . -5 0 GeV/c and . | η | = 1 0. .

[28]. As mentioned earlier, one way to reduce the fluctuation effects is to reduce the number of participating nucleons by selecting the narrow centrality bins. However, the choices of centrality should be such that it does not introduce additional fluctuations due to finite multiplicity and detector resolutions. Once the statistical fluctuation part is under control, the fluctuation effects present will be mostly of dynamical origin which may contain interesting physics associated with the collisions, like hydrodynamic expansion, hadronization at freeze-out, etc.

Variation of scaled variance with c.m. energy for different centrality bins are plotted in FIG.8. It may be noted from the figure that ω increases with beam energy as well as in centrality bin widths. It may also be noted that increase of ω with c.m. energy becomes linear for the centrality classes 35% and above. If the data obey the KNO scaling [21], it is predicted that ω should increase linearly with mean charge multiplicity [29]. It may also be noticed in FIG.8 that increase of ω with beam energy is somewhat weaker for the central collisions. Similar trends of variations of ω with energy have also been reported in pp collisions by NA61 collaboration [29].

Centrality dependence of scaled variance at the four incident energies are exhibited in FIG.9. It is observed that for 10% centrality bins ω increases with centrality bin widths, whereas for 5%

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-33fcb958-5614-4d15-8094-dc48274c2c01.jpg)

ch

FIG. 5: Multiplicity distributions of relativistic charged particles for various centrality classes at 40A GeV in the range p T = 0 2 . -5 0 GeV/c and . | η | = 1 0. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-5282a5d2-51da-4513-80a0-2636ca556bcb.jpg)

ch

FIG. 6: Scaled Multiplicity distributions of relativistic charged particles for the centrality bins 0-10%, 3040% and 50-60%. Distributions corresponding to 30-40% and 50-60% are shifted up on the y-scale by factors 10 and 100 for clarity sake only.

FIG. 7: Scaled Multiplicity distributions of relativistic charged particles for the minimum bias events at the four beam energies considered.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-79f190ad-2ccd-4c99-8c2b-474686ecadb6.jpg)

ω

FIG. 8: Dependence of scaled variance, ω on beam energy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-5f6cb9be-98ba-41cb-8192-706ab64ea8b7.jpg)

and 2% this parameter slowly decreases with increasing centrality and thereafter tends to acquire nearly constant values. This observation, thus, supports that statistical fluctuations arising due to fluctuations in N part becomes visible if the centrality bin width is ∼ 10% or more and hence

ω

FIG. 9: Dependence of scaled variance, ω on &lt; N part &gt; in different centrality-bin widths.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00233v1-figure-8929c57b-a20b-40c8-a19e-57cffd5c2c05.jpg)

considering a bin as wide as 5%, would help arrive at some meaningful conclusions on dynamical fluctuations, if present.

## 4. Conclusions

MDs and ebe multiplicity fluctuations in AuAu collisions from the beam energy scan in future heavy-ion experiment at the Facility for Antiproton and Ion Research ( FAIR ) are examined in the frame work of Ultra-Relativistic Quantum Molecular Dynamics model, URQMD . The mean values of MDs are observed to shift towards the higher multiplicity and the width of the distributions are found to become wider from central to peripheral collisions. The MDs are also observed to obey KNO scaling in various centrality windows as well as for full event (minimum bias) samples. Centrality-bin width dependence of the 2 nd moments and scaled variance gives the idea of bin width

effect and centrality window-width selection, where the statistical fluctuations may be treated as under control.

## References

- 1 S.A. Voloshin, V. Koch and H. G. Ritter, 'Event-by-event fluctuations in collective quantities', Phys. Rev. C60 (1999) 024091.
- 2 Shakeel Ahmad et al, 'A study of event-by-event fluctuations in relativistic heavy-ion collisions', Int. J. Mod. Phys. E23 (2014) 1450065.
- 3 J.L. Albacete et al, 'The initial state of heavy ion collisions', Int. J. Mod. Phys. A28 (2013) 1340010.
- 4 M.A. Stephanov et al, 'Event-by-event fluctuations in heavy ion collisions and the QCD critical point', Phys. Rev. D60 (1999) 114028.
- 5 M.A. Stephanov et al, 'Signatures of the tri-critical point in QCD', Phys. Rev. Lett. 81 (1998) 4816.
- 6 M.A. Stephanov, 'QCD phase diagram and the critical point', Int. J. Mod. Phys. A20 (2005) 4387.
- 7 M. Mukherjee et al, 'Fluctuations in charged particle multiplicities in relativistic heavy-ion collisions', J Phys. G: Nucl. Part. Phys. 43 (2016) 085102, arXiv:1603.02083v3[nucl-ex].
- 8 J. B¨chler et al (NA35 Collaboration), 'Fluctuations of multiplicities in rapidity windows a in sulphur-sulphur collisions at 200A GeV', Z. Phys. C56 (1992) 347.
- 9 A. Bialas and R. Peschanski, 'Moments of rapidity distributions as a measure of short-range fluctuations in high-energy collisions', Nucl. Phys. B273 (1986) 703.
- 10 R.C. Hwa, 'A proposed analysis of multiplicity fluctuations in high-energy heavy-ion collisions', Phys. Lett. B201 (1988) 165.
- 11 R.C. Hwa, 'Enhanced multiplicity fluctuation as a possible signature of quark matter', Prog. Part. Nucl. Phys. 41 (1988) 277.
- 12 S.A. Bass et al, 'Microscopic models for ultrarelativistic heavy ion collisions', Z. Phys. C38 (1998) 255; arXiv:9803035v2 [nucl-th].
- 13 M. Bleicher et al, 'Relativistic hadron-hadron collisions in the ultra-relativistic quantum molecular dynamics model', J. Phys . G25 (1999) 1859; arXiv:9909407v1 [hep-ph].

- 14 S.Z. Belen kji and L.D. Landau, 'Hydrodynamic theory of multiple production of particles', ′ Nuovo Cim. Supp. 3 (1956) 15.
- 15 T. Regge, 'Introduction to complex orbital momenta', Nuovo Cim. 14 (1959) 951.
- 16 H. St¨cker and W. Greiner, 'High energy heavy ion collisions-probing the equation of state o of highly excited hardronic matter', Phys.Rep. 137 (1986) 277.
- 17 S. Ghosh et al, 'Net-charge fluctuation in Au+Au collisions at energies available at the Facility for Antiproton and Ion Research using the UrQMD model' Phys. Rev. C96 (2017) 024912.
- 18 V.P. Konchakovski et al, 'Multiplicity fluctuations in nucleus-nucleus collisions: Dependence on energy and atomic number', Phys. Rev. C76 (2008) 024906, arXiv:0712.2044v2 [nucl-th]
- 19 M.L. Miller et al, 'Glauber modeling in high-energy nuclear collisions', Annu. Rev. Nucl. Post. Sci. 57 (2007) 205.
- 20 A. Shakeel et al, 'Scaling of multiplicity distribution of charged shower particles in protonnucleus interactions at 400 GeV', Phys. Scr. 29 (1984) 435.
- 21 Z. Koba, H. B. Neclson and P. Olesen, 'Scaling of multiplicity distributions in high energy hadron collisions', Nucl. Phys. B40 (1972) 317.
- 22 P. Slattery, 'Evidence for the onset of semi-inclusive scaling in proton-proton collisions in the 50 - 300 GeV/c momentum range', Phys. Rev. Lett. 29 (1972) 1624.
- 23 J.W. Martin et al, 'Scaling of multiplicity and angular distributions in p-emulsion interactions at 30, 67 and 200 GeV', Nuovo Cimento A25 (1975) 447.
- 24 P. Olesen, 'The 200 GeV multiplicity distribution and scaling', Phy. Lett. B41 (1972) 602.
- 25 A.J. Buras et al, 'Multiplicity scaling at low energies, a generalized Wroblewski-formula and the leading particle effect', Phys. Lett. B47 (1973) 251.
- 26 A. Breakstone et al, 'Charged multiplicity distribution in pp interactions at CERN ISR energies', Phys. Rev. D30 (1984) 528.
- 27 M. Mukherjee, 'Multiplicity distributions and fluctuations in proton-proton and heavy-ion collisions', Euro. Phys J. Web. Conf. 112 (2016) 04004.

- 28 Andrey Seryakov, NA61/SHINE Coll., 'Rapid change of multiplicity fluctuations in system size dependence at SPS energies', KnE Energy and Phys., ICPPA 2017, The 3rd Int. Conf. on Part. Phys. and Astrophy. , (2018) 170; arXiv:1712.03014v1 [hep-ex].
- 29 Maja. Mackowiak-Pawlowska and A. Wilczek, NA61 Coll., 'Multiplicity fluctuations of identified hadrons in p+p interactions at SPS energies', J. of Phys Conf. Series 509 (2014) 012044; arXiv:1402.0707v1 [hep-ex].
Anuj Chandra Bushra Ali Shakeel Ahmad
Centrality Dependence of Multiplicity Fluctuations in Ion-Ion Collisions from the Beam Energy Scan at FAIR
Multiplicity distributions and event-by-event multiplicity fluctuations in
AuAu collisions at energies in future heavy-ion experiment at the Facility for
Anti-proton and Ion Research (FAIR) are investigated. Events corresponding to
FAIR energies are simulated in the frame work of Ultra Relativistic Quantum
Molecular Dynamics (URQMD) model. It is observed that the mean and the width of
multiplicity distributions monotonically increase with beam energy. The trend
of variations of dispersion with mean number of participating nucleons for the
centrality-bin width of 5\% are in accord with the Central Limit Theorem. The
multiplicity distributions in various centrality bins as well as for full event
samples are observed to obey Koba, Nielsen and Olesen (KNO) scaling. The trends
of variations of scaled variance with beam energy are also found to support the
KNO scaling predictions for larger collision centrality. The findings also
reveal that the statistical fluctuations in 5\% centrality-bin width appear to
be under control.