Abstract

Entropy production in multiparticle systems is investigated by analyzing the experimental data on ion-ion collisions at AGS and SPS energies and comparing the findings with those reported earlier for hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions. It is observed that the entropy produced in limited and full phase space, when normalized to maximum rapidity, exhibits a kind of scaling which is nicely supported by Monte Carlo model HIJING. Using Rényi’s order information entropy, multifractal characteristics of particle production are examined in terms of generalized dimensions, . Nearly the same values of multifractal specific heat, , observed in hadronic and ion-ion collisions over a wide range of incident energies suggest that the quantity might be used as a universal characteristic of multiparticle production in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions. The analysis is extended to the study of spectrum of scaling indices. The findings reveal that Rényi’s order information entropy could be another way to investigate the fluctuations in multiplicity distributions in terms of spectral function , which has been argued to be a convenient function for comparison sake not only among different experiments but also between the data and theoretical models.

1. Introduction

The main objective of investigating the collisions of heavy nuclei at relativistic energies is to understand the nature of phase transition from normal hadronic matter to quark-gluon plasma (QGP) [1]. Dedicated experiments have been carried out at AGS and RHIC at BNL and SPS and LHC at CERN to search for the QGP phase and explore the QCD phase diagram. Heavy-ion collision experiments provide a unique opportunity to test the predictions of QCD and at the same time to understand the soft processes involved, which dominate even at LHC energies, as well as the hard scattering or the small cross-sectional physics [24]. Relativistic nucleus-nucleus (AA) collisions serve as a mean to study the novel regime of QCD where parton densities are high enough while the strong coupling constant between the partons is small and decreases further with decreasing interpartonic distances [2]. The parton densities in the early stage of the collision can be related to the density of the produced charged hadrons in the final state. Dominance of hard process (jets and minijets) in the hadron production increases with increasing collision energy and hence provides a unique opportunity to study the interplay between effects. In this scenario, the perturbative QCD (pQCD) lends a good basis for high energy dynamics and has achieved significant success in describing the hard processes involved in high energy collisions [2, 5]. Hadronic-jet production in annihilations [6, 7] and large jet production in hadron-hadron (hh) collisions [8, 9] are the examples of these hard processes. However, in soft processes like hadron production with sufficiently low in hadronic and heavy-ion collisions, the interactions become so strong that the pQCD is not applicable anymore [2]. Due to incapability of pQCD in this regime, phenomenological models based on experimental inputs have been proven to be an alternative tool to understand the dynamics of particle production in AA collisions.

Relativistic charged particle multiplicity is the simplest observable and by studying its distribution for a given data sample, information on soft QCD processes as well as on hard scattering can be extracted [3]. Furthermore, the charged hadron multiplicity of an event is taken as a direct measure of its inelasticity and can describe important features of particle production [3]. Investigations involving multiparticle production in hh, hA, and AA collisions at relativistic energies have been carried out by numerous groups during the last four decades [1015]. However, a complete understanding of particle production mechanism still remains elusive. Multiplicity distributions (MD) of relativistic charged particles produced in hh collisions have been observed to deviate from a Poisson distribution and are expected to provide information regarding the underlying production mechanism [13, 16]. Asymptotic scaling of MD in hh collisions referred to as the KNO scaling [17], predicted in 1972, was regarded as a useful phenomenological framework for comparing the MD at different energies ranging from  GeV to ISR energies [18]. It was, however, pointed out [19, 20] that KNO scaling is not strictly followed for inelastic hh collisions. This scaling law was observed to breakdown when collision energy reached SPS range [13, 16, 18, 21]. After the observations of KNO scaling violation in collisions at = 540 GeV, it was remarked that the observed scaling up to ISR energies was approximate and accidental [21]. A new empirical regularity, in place of KNO scaling, was then proposed [21] to predict the multiplicity distributions at different energies. It was shown that MD at different energies in full and limited rapidity () windows may be nicely reproduced by negative binomial distributions (NBD) [18, 21, 22]. These observations lead to revival of interest in investigations involving MD and new scaling laws. Simak et al. [23], by introducing a new variable, the information entropy, showed that MD of charged particles produced in full and limited phase space in hh collisions exhibits a new type of scaling law in the energy range, to 900 GeV.

Sinyukov and Akkelin [24] have proposed a method to estimate entropy of thermal pions in AA collisions and have studied the average phase space densities and entropy of such pions against their multiplicities and beam energies. Their findings apparently suggest the presence of deconfinement and chiral phase transition in relativistic AA collisions. Moreover at RHIC energies, entropy per unit rapidity at freeze-out has been extracted with minimal model dependence from the available measurements of particle yields, spectra, and source sizes, determined from two-particle interferometry [25]. The extracted entropy per unit rapidity was observed to be consistent with the lattice gauge theory for thermalized QGP with an energy density calculated from the transverse energy production at RHIC energies.

Analyses of the experimental data on , , and collisions over a wide energy range (up to = 900 GeV) carried out by several workers [23, 26, 27] indicate that entropy increases with beam energy while the entropy per unit rapidity appears to be an energy independent quantity. These results indicate the presence of entropy scaling up to a few TeV energy. Presence of a similar scaling behaviour in collisions at LHC energies has also been reported by Mizoguchi and Biyajima [28] and Das et al. [11, 13]. Analyses of AA collision data at AGS and SPS energies carried out by other workers [10, 15, 29, 30] too suggest that entropy produced in limited pseudorapidity () windows when normalized to the maximum rapidity is, essentially, independent of projectile and target mass as well as the beam energy, indicating the presence of entropy scaling.

The occurrence of unusual large particle density fluctuations in narrow phase space bins, observed in cosmic ray JACEE events [31] and in accelerator experiments [32, 33], have generated considerable interest in the study of nonlinear phenomena in hadronic and heavy-ion collisions. Such fluctuations may be taken as an indication of a phase transition from ordinary hadronic matter to QGP, predicted by QCD to occur in relativistic AA collisions. These fluctuations are also envisaged to arise either due to minijets produced at very high energies or (and) because of some other collective phenomena [34]. Such rare fluctuations of dynamical origin are to be identified and extracted from the statistical ones. Various methods for the identification of these fluctuations have been proposed which estimate the total deviation of the measured rapidity distribution from an idealized smooth distribution [34]. The analysis is based on the estimation of corresponding statistical probability through Monte Carlo (MC) simulations. The method of scaled factorial moments (SFMs), proposed by Bialas and Peschanski [35], has been proven to be well suited for not only to search for the dynamical fluctuations in narrow phase space bins but also to investigate the pattern of fluctuations which could lead to a physical interpretation of their origin [35]. The concept of SFMs was put forward in analogy with the phenomenon referred to as the “intermittency” in hydrodynamics of turbulent fluid flow [34]. This phenomenon is characterized by the presence of fluctuations in a small part of the available phase space. The method of SFMs was first applied to the pseudorapidity distribution of a single high multiplicity JACEE event [31] and the findings indicated the presence of intermittent pattern (large particle density fluctuations in narrow pseudorapidity bins) [35]. SFMs analysis, since its introduction, has been widely used to search for the nonlinear phenomena in hadronic and ion-ion collisions. It has been observed [3639] that the presence of such a nonlinear phenomena might be rare but not impossible [40]. Such fluctuations, if not arising due to the statistical reasons, are envisaged to occur due to the dynamical correlations among the produced particles. These correlations might arise due to the phase transition from QGP to normal hadronic matter. As pointed out by Bialas and Peschanski [35], if the intermittency exists, the SFMs of multiplicity distributions should exhibit a power law dependence of the rapidity-bin width, as . Such a study would help understand the chaotic behaviour of rapidity distribution on event-by-event (ebe) basis instead of investigating the average phenomena [41].

The observations of intermittency in annihilation [42], hh [43], hA [34], and AA collisions [41] attracted a great deal of attention towards the investigations involving the power law behaviour of SFMs of the form where the exponent increases with order of the moment [41].

Although the intermittency analysis in terms of SFMs has been successfully applied to the hadronic and heavy-ion collisions, yet the dynamical explanation of its origin in some cases is not clear [44]. The concept of self-similarity is closely related to the fractal theory [45, 46] which is a natural consequence of the cascading mechanism prevailing in the multiparticle production. A formalism for treating the fractal dimensions and its generalization has been developed and applied to the study of turbulent fluids and other transitions to chaos [41, 47]. However, in the case of multiparticle production, the dynamics is not well known and one has to discover the proper dynamics that gives rise to the chaotic structure. The fractal dimensions, in multiparticle production, were first studied by Carruthers and Duong-Van [48]. Later on Dermin [49] suggested the study of correlation dimensions, while Lipa and Buschbeck [50] considered other generalized dimensions. However, in none of these investigations a formalism for systematic studies of fractal properties has been presented. Considering inelastic collisions as purely geometrical objects with noninteger dimensions, a formalism for treating the fractal dimensions and its generalization was developed and successfully applied to the study of intermittent behaviour and other transitions to chaos [44, 47, 5153]. In order to reveal the multifractal nature of multiparticle production, Hwa [41], Chiu and Hwa [54], and Florkowski and Hwa [55] introduced moments, which can be estimated from the particle multiplicities in limited rapidity windows that might lead to inferences on the nature of the dimensions, , that are the generalizations of the fractal dimensions for multifractal sets. Using the moments properties of scaling indices spectrum may be investigated which, in turn, would indicate as to how it can provide effective means to describe a highly nonuniform rapidity distribution.

As the studies involving AA collisions at relativistic energies are concerned, the main goal is to study the properties of strongly interacting matter under extreme conditions of nuclear density and temperature, where formation of quark-gluon plasma (QGP) is predicted to take place [10, 15, 56, 57]. Fluctuations in physical observables in AA collisions are regarded as one of the important signals, for QGP formation because of the idea that, in many body systems, phase transition may result in significant changes in the quantum fluctuations of an observable from its average behaviour [15, 29, 56]. For example, when a system undergoes a phase-transition, heat capacity changes abruptly, whereas the energy density remains the smooth function of temperature [10, 15, 25, 58, 59]. Entropy is regarded yet another important characteristics of a system with many degrees of freedom [15, 6063]. Processes in which particles are produced may be considered as the so-called dynamical systems [6064] in which entropy is generally produced. Systematic measurements of local entropy produced in AA collisions may provide direct information about the internal degrees of freedom of the QGP medium and its evolution [10, 15, 27]. It has been pointed out [65] that in high energy collisions particle production occurs on the maximum stochasticity; that is, they should follow the maximum entropy principle. This type of stochasticity may also be quantified in terms of information entropy which may be regarded as a natural and more general parameter to measure the chaoticity in branching processes [66].

Bialas et al. [60, 61, 67] proposed that Rényi entropies may also be used as a tool for studying the dynamical systems and are closely related to the thermodynamic entropy of the system, the Shannon entropy. Bialas et al. [68] have also suggested that Rényi’s entropies may serve as a useful tool for examining the correlation among the particles produced in high energy collisions. Furthermore, the generalization of Rényi’s order- information entropy contains information on the multiplicity moments that can be used to investigate the multifractal characteristics of particle production [14, 69, 70]. It should be mentioned that this method of multifractal studies is not related to the phase space bin-width or the detector resolution but to the collision energy [69]. It is, therefore, considered worthwhile to carry out a well focused study of entropy production and subsequent scaling in AA collisions by analyzing the several sets of experimental data on AA collisions at AGS and SPS energies. Rényi’s information entropies are also estimated to investigate the multifractal characteristics of multiparticle production. Moreover, the analysis is further extended to study the fractal dimensions, multifractal spectrum, and their dependence on the energy and mass of the projectile.

2. Details of Data

Six sets of events produced in collisions of 16O, 28Si, 32S, and 197Au beams with AgBr group of nuclei in emulsion at AGS and SPS energies, available in the laboratory, are used in the present study. Details of these samples are presented in Table 1. These events are taken from the emulsion experiments performed by EMU01 Collaboration [7174]. The other relevant details of the data, like selection criteria of events, classification of tracks, extraction of AgBr group of interactions, method of measuring the emission angle, of relativistic charged particles, and so on, may be found elsewhere [10, 15, 71, 7577].

From the measured values of the emission angle , the pseudorapidity variable, , was calculated using the relation, . It is worthwhile to mention that the conventional emulsion technique has two main advantages over the other detectors: (i) its solid angle coverage and (ii) emulsion data being free from biases due to full phase space acceptance. In the case of other detectors, only a fraction of charged particles is recorded due to the limited acceptance cone. This not only reduces the multiplicity but also may distort some of the event characteristics like particle density fluctuations [15, 75, 78]. In order to compare the findings of the present work with the predictions of Monte Carlo model HIJING [79, 80], event samples corresponding to the experimental ones are simulated using the code HIJING-1.35; the number of events in each simulated sample is equal to that in the corresponding real event sample. The events are simulated by taking into account the percentage of interactions which occur in the collision of projectile with various targets in emulsion which constitute the AgBr group of nuclei [71, 74, 75]. The value of impact parameter for each data sample was so set that the mean multiplicity of the relativistic charged particles becomes nearly equal to those obtained for the real data sets.

3. Formalism

For entropy of the charged particle multiplicity distribution, Shannon’s information entropy is calculated using [10, 23] and its generalization, Rényi’s order information entropy, is estimated as [14, 26, 70] where for   , and is the probability of production of charged particles. The generalized dimensions of order may be estimated as [14, 69, 70] denotes the maximum rapidity in the centre-of-mass frame, represents the center-of-mass energy, is the pion rest mass, denotes the average number of participating nucleons, and is maximum multiplicity of relativistic charged particles produced for a given pair of colliding nuclei at a given energy.

For the particle production process to follow self-similar behaviour, the moments of order , defined as [81] should follow the power law of the formwhere is the pseudorapidity bin-width. The parameter is related to the dimension , for all , bywhere , , and are usually referred to as the fractal dimensions, information dimension, and correlation dimension, respectively [41]. It has been pointed out by Hwa [41] that values of moments obtained from various experiments can not be compared as they depend on the number of events in the data samples and on , that is, on detector resolution. It is the dependence of that is the aim of studying ; in particular data analysis should aim to extract the generalized dimensions , where The meaning of function becomes rather more obvious after performing the Legendre transformation from independent variables and to the variables and : is fractal dimension of a subset composed from bins whose occupancy probability lies in the interval ( to ). Thus by estimating the moments, the continuous scaling function can be constructed [44].

The thermodynamical interpretation of these relationships means that can be interpreted with an inverse temperature whereas the spectrum and play the role of entropy and energy (per unit volume), respectively [8184].

4. Results and Discussion

Probability of production of charged particles in a pseudorapidity window of fixed width is calculated by selecting a window of width . This window is so chosen that its midposition coincides with the center of symmetry of distribution, . Thus, all the relativistic charged particles having their values lying in the range are counted to evaluate . The window size is then increased in steps of 0.5 until the region, , is covered. Values of entropy, , for different -windows are calculated by using (2), while the value of maximum rapidity is estimated from (4). Variations of entropy normalized to maximum rapidity, , with -window width, also normalized to maximum rapidity, , for the experimental and HIJING data sets are plotted in Figure 1. It is observed in the figure that first increases up to and thereafter acquires nearly a constant value. It is interesting to note that the data points for various samples of events overlap to form a single curve. This indicates the presence of entropy scaling in AA collisions at AGS and SPS energies. Results from HIJING simulated events also support the presence of entropy scaling.

Similar entropy scaling in AA collisions has been reported by us [10, 15] and also by the other workers [29, 70] for central and minimum bias events. In our earlier work, attempt was made to look for whether the observed entropy scaling is of dynamical nature. For this purpose the correlation free Monte Carlo events samples (Mixed events) corresponding to each of the real data samples were generated and analyzed. These findings revealed that the entropy scaling observed was the distinct feature of the data and was of dynamical origin [10, 15].

According to (4) the quantity, , is equal to the maximum charged particle multiplicity at a given center of mass energy. It would, therefore, be convenient to examine the mean multiplicity versus entropy in limited windows as well as in full range; the entropy in the entire range, , is calculated using (2). Variations of with for the experimental and HIJING events are exhibited in Figure 2: denotes the total energy of the beam nucleus. It may be noticed in the figure that increases linearly with for both real and simulated event samples. HIJING, however, predicts somewhat smaller values of as compared to the corresponding values estimated from the real data. These findings are, thus, in agreement with those reported for 16O-nucleus collision at 60A and 200A GeV/c and S-nucleus interaction at 200 GeV/c [70]. Shown in Figure 3 are the variations of with for the six data sets considered. It is evident from the figure that the data points corresponding to various types of collisions almost overlap to form a single curve. It may also be noted that as . These observations too support the presence of entropy scaling in AA collisions at AGS and SPS energies.

It has been reported [14, 81, 85, 86] that the constant specific heat (CSH) approximation widely used in standard thermodynamics is applicable to the multifractal data too. Nearly constant value of multifractal specific heat has been observed [8790] from the analysis of the experimental data on hadronic and heavy-ion collisions at different energies. These investigations have been carried out by following the methods proposed by Hwa [41] and Takagi [88], where the multiplicities of particles in narrow -bins, , are involved. This introduces a limit to analysis because would depend on the detector resolution (see (6)). However, since the quantity scales with improving relative resolution [70] like , the Rényi’s order- information entropy may also be used to understand the multifractal nature of particle production and estimation of multifractal specific heat. The advantage of this method is that it does not depend on the phase space bin width and hence on detector resolution; rather it accounts for the fractal resolution, which is related to the collision energy [14]. From the definitions of Rényi’s information entropy, , and generalized dimensions, (see (3) and (4)), it is evident that, for a given , = . The highest entropy is achieved for the greatest “chaos” of a uniformly distributed probability function [71] and (4) gives .

Variations of with for various event samples are shown in Figure 4. The values of have been obtained from the spectrum (Figure 6) and discussed in the coming part of the text. It is observed that monotonously decreases with increasing order and the trend of decrease for the real and HIJING events is nearly the same except that HIJING predicts somewhat smaller values of as compared to that for experimental data. The spectrum for order , which for multifractals is decreasing function of , can be related to the scaling behaviour of point correlation integrals [14, 47]. Thus, the trend of variations of against observed in the present study indicates the multifractal nature of multiplicity distributions in full phase space in ion-ion collisions at the energies considered. Similar variations of with have also been reported earlier [70] in collisions of protons (800 GeV) and heavy ions (4He at 11A, Si at 14.5A,O at 60A and 200A, and S at 200A GeV/c) collisions. Although the presence of multifractality predicts the decrease of with , yet no further useful information about the -dependence of spectrum can be extracted which may lead to making remarks on the scaling properties of -correlation integrals [14]. It has been suggested [81, 85, 90] that in constant heat approximation, , dependence on acquires the following simple form:where is the information dimension, , while denotes the multifractal specific heat. The linear trend of variation with , given by (10), is expected to be observed for multifractals. On the basis of classical analogy with specific heat of gases and solids, the value of is predicted [89] to remain independent of temperature in a wide range of . In order to test the validity of (10), values are plotted against in Figure 5. It is evident from the figure that increases linearly with . The lines in the figure are due to the best fits to the data obtained using (10). The values of the parameters “” and “,” thus, obtained are listed in Table 2. It is interesting to note from the table that the values of multifractal specific heat, , for all the data sets, are nearly the same, ~0.2 there by indicating its independence of the beam energy and projectile mass. It may also be noticed from the table that the experimental values of are quite close to those predicted by HIJING. It is worthwhile to mention that the values of multifractal specific heat obtained in the present study are close to those obtained by us [87] by analysing some of these data sets using Takagi’s approach [88]. Incidentally similar values of have been reported by Du et al. [90] for 10.6A GeV/c 197Au-nucleus collisions. In -nucleus interactions too, the value of multifractal specific heat has been observed to be ~0.25 in the energy range, 200–800 GeV [70, 85, 90]. However, in the case of collisions the values of this parameter have been found to be ~0.08 in the wide energy range 25–1800 GeV [14].

These findings, thus, indicate that the constant-specific heat approximation is applicable to the multiparticle production in relativistic hadronic and ion-ion collisions. Moreover, nearly the same values of multifractal specific heat, , observed in the present study as well as other workers using the data on hh, hA, and AA collisions involving various projectiles/targets at widely different energies do indicate that the parameter may be regarded as a universal characteristics of hadronic and heavy-ion collisions.

In order to construct the multifractal spectrum, values of , , and are evaluated using (7) and (9) for to 6. spectra for various data sets are displayed in Figure 6. It is evidently clear from the figure that these spectra are represented by continuous curves, thus, characterising a qualitative manifestation of the multiplicity fluctuations. Spectra for all the data sets (real and HIJING ) are noticed to follow the general characteristics of occurrence of peaks at and a common target at an angle of at . The spectrum is noticed to have a peak at and is concave downwards everywhere. As expected the values of decrease with increasing as decreases with increasing . The region corresponds to positive and the curves in this region have positive slope, whereas the region would correspond to negative and the slope of the curves in this region is negative.

The values of for would give the fractal dimensions (), the information dimension (), and the correlation dimension (). The values of information dimension are shown in Figure 4 to distinguish these points in the figure; the data points are encircled. The widths of the spectrum () for various data sets (real and HIJING ) are also calculated and presented in Table 3. The width of the spectrum has been suggested as a measure of the degree of multifractality [9193]. The broader the spectrum, the higher the multifractality [93, 94].

It may be observed from Figure 6 that the spectrum becomes wider with increasing beam energy (16O-AgBr data at 14.5A, 60A, and 200A GeV/c) as well as with increasing projectile mass (16O-AgBr and 28Si-AgBr data at 14.5A GeV/c and 16O-AgBr and 32S-AgBr data at 200A GeV/c). Similar dependence of spectrum has also been observed in 12C-, 24Mg-, 16O-, and 32S-AgBr collisions at 2.1A, 4.5A, 60A, and 200A GeV/c [93].

These observations, thus, tend to suggest that multifractality is more pronounced if the beam energy or mass of the colliding nuclei increases. However, it has been pointed out that at higher energies influence of energy is more prominent as compared to the mass [93]. Hwa [41] has also pointed out that since increases with energy for all values, and will increase with energy too and hence the entire spectrum () would be a broader once. This, in turn, implies that the rapidity distribution will become more jagged and irregular with prominent sharp peaks and deep valleys. As the analysis method suggested by Hwa [41] is based on the multiplicity fluctuations in limited rapidity bins, the remark made in [41] that the highly chaotic behaviour of MD in narrow rapidity bins, , can be observed by a smooth function is quite interesting. It has also been remarked that averaging over all events would be devoid of any information about fluctuation and intermittency and also that if the experimental/detector resolution is not good enough to capture the sharp peaks and deep valleys, the spectrum would appear narrower than it should otherwise be [41]. In contrast to the standard intermittency and multifractal analysis which are based on the multiplicities in phase space bins, the present analysis, as mentioned earlier, does not depend on the experimental resolution. It rather considers fractal resolution related to the total energy available [14]. Assuming that, during collisions with many particles produced, energy dissipates into discrete sites. The site labeled by is occupied by the probability . Since most of the sites are unoccupied, the overlay of many inelastic events can be visualized as a fractal with overall extent . The sufficient condition to produce self-similar structure is that would exhibit some type of scaled invariant behaviour. At sufficiently high energies, which correspond to high enough resolution , acquires a power law dependence on the resolution and hence the quantity will scale with as [14]The independence of generalized dimensionon energy (resolution ) would refer the presence of multifractality in MD of relativistic charged particles produced in high energy collisions [14, 95].

5. Conclusions

On the basis of the findings of the present work the following conclusions may be reached:(1)The entropy normalized to maximum rapidity exhibits a saturation beyond indicating thereby the presence of large amount of entropy around midrapidity region.(2)Overlapping data points for various sets of events in vs plots exhibit entropy scaling in AA collisions at AGS and SPS energies.(3)The scaling observed with the experimental data is nicely supported by HIJING model.(4)The decreasing values of with increasing order number may be taken as a signal of multifractal nature of multiplicity distributions of relativistic charged particle produced in hadronic and heavy-ion collisions.(5)Besides the information dimension, , there is yet another parameter, , which may be used as a universal parameter for multiparticle production in high energy hadronic and heavy-ion collisions.(6)Rényi’s order -information entropy may also be used to construct the multifractal spectrum, which in turn would help study the multiplicity fluctuations and scaling in high energy hadronic and heavy-ion collisions.

Conflicts of Interest

The authors declare that there are no conflicts of interest.