Properties of Chemical and Kinetic FreezeOuts in HighEnergy Nuclear Collisions
View this Special IssueResearch Article  Open Access
Multifractal Analysis of Charged Particle Multiplicity Distribution in the Framework of Renyi Entropy
Abstract
A study of multifractality and multifractal specific heat has been carried out for the produced shower particles in nuclear emulsion detector for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c in the framework of Renyi entropy. Experimental results have been compared with the prediction of UltraRelativistic Quantum Molecular Dynamics (UrQMD) model. Our analysis reveals the presence of multifractality in the multiparticle production process in high energy nucleusnucleus interactions. Degree of multifractality is found to be higher for the experimental data and it increases with the increase of projectile mass. The investigation of quarkhadron phase transition in the multiparticle production in ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5 AGeV/c in the framework of GinzburgLandau theory from the concept of multifractality has also been presented. Evidence of constant multifractal specific heat has been obtained for both experimental and UrQMD simulated data.
1. Introduction
The study of nonstatistical fluctuations and correlations in relativistic and ultrarelativistic nucleusnucleus collisions has become a subject of major interest among the particle physicists. Bialas and Peschanski [1, 2] proposed a new phenomenon called intermittency to study the nonstatistical fluctuations in terms of the scaled factorial moment. In high energy physics intermittency is defined as the power law behavior of scaled factorial moment with the size of the considered phase space [1, 2]. This method has its own advantage that it can extract nonstatistical fluctuations after extricating the normal statistical noise [1, 2]. The study of nonstatistical fluctuations by the method of scaled factorial moment leads to the presence of selfsimilar fractal structure in the multiparticle production of high energy nucleusnucleus collisions [3]. The selfsimilarity observed in the power law dependence of scaled factorial moments reveals a connection between intermittency and fractality. The observed fractal structure is a consequence of selfsimilar cascade mechanism in multiparticle production process. The close connection between intermittency and fractality prompted the scientists to study fractal nature of multiparticle production in high energy nucleusnucleus interactions. To get both qualitative and quantitative idea concerning the multiparticle production mechanism fractality study in heavyion collisions is expected to be very resourceful.
2. FractalsMultifractals and Monofractals
The term “fractal” was coined by Mandelbrot [4] from the Latin word fractus which means broken or fractured. Mandelbrot introduced the new geometry called the fractal geometry to look into the world of apparent irregularities. A fractal pattern is one that scales infinitely to reproduce itself such that the traditional geometry does not define it. In other words a fractal is generally a rough or fragmented geometrical shape that can be split into parts, each of which is at least approximately reduced size copy of the whole [4]. Sierpinski triangle, Koch snow flack, Peano curve, Mandelbrot set, and Lorentz attractor are the wellknown mathematical structures that exhibit fractal geometry [4]. Fractals also describe many real world objects such as clouds, mountains, coastlines, etc.; those do not corresponds to simple geometrical shape [4]. Fractals can be classified into two categories: multifractals and monofractals [4]. Multifractals are complicated selfsimilar objects that consist of differently weighted fractals with different noninteger dimensions. The fundamental characteristic of multifractality is that the scaling properties may be different in different regions of the systems [4]. Monofractals are those whose scaling properties are the same in different regions of the system [4]. As the scaling properties are dissimilar in different parts of the system, multifractal systems require at least more than one scaling exponent to describe the scaling behavior of the system [5]. A distinguishing feature of the processes that have multifractal characteristics is that various associated probability distributions display power law properties [6, 7]. Multifractal theory is essentially rooted in probability theory though draws on complex ideas from each of physics, mathematics, probability theory, and statistics [6]. Apart from multiparticle production in high energy physics, multifractal analysis has proved to be a valuable method of capturing the underlying scaling structure present in many types of systems including diffusion limited aggregation [8–10], fluid flow through random porous media [11], atomic spectra of rareearth elements [12], clustercluster aggregation [13], and turbulent flow [14]. In physiology, multifractal structures have been found in heart rate variability [15] and brain dynamics [16]. Multifractal analysis has been helpful in distinguishing between healthy and pathological patients [17]. Multifractal measures have also been found in manmade phenomena such as the Internet [18], art [19], and the stock market [20–22]. Multifractals have also been used in a wide range of application areas like the description of dynamical systems, rainfall modelling, spatial distribution of earthquakes and insect populations, financial time series modelling, and Internet traffic modelling [6, 7].
It should be mentioned here that the most important property of fractals is their dimensions [4, 6, 7]. Fractal dimension is used to describe the size of the fractal sets [4, 6, 7]. For example, the dimension of an irregular coastline may be greater than one but less than two, indicating that it is not like a simple line and has space filling characteristics in the plane. Likewise, the surface area of a snowflake may be greater than two but less than three, indicating that its surface is more complex than regular geometrical shapes and is partially volume filling [4, 6, 7]. Fractal dimension can be calculated by taking the limit of the quotient of the log change in object size and the log change in measurement scale, as the measurement scale approaches to zero [4]. The differences came in what is meant by the measurement scale and how to get an average number out of many different parts of the geometrical objects [4, 6, 7]. Fractal dimension quantifies the static geometry of an object [4].
Generalized fractal dimension is a wellknown parameter which reflects the nature of fractal structure [4]. From the dependence of generalized fractal dimension on the order q a distinguishable characterization of fractality is possible [4]. Decrease of generalized fractal dimension with the order of moment signals the presence of multifractality. On the other hand if remains constant with the increase of order q monofractality occurs. It has been pointed out in [23] that if QuarkGluonPlasma (QGP) state is created in hadronic collisions, a phase transition to hadronic matter will take place. The hadronic system will show monofractality in contrast to an order dependence of generalized fractal dimension , when a cascade process occurs [23].
Hwa [24] was the first to provide the idea of using multifractal moments , to study the multifractality and selfsimilarity in multiparticle production. According to the method enunciated by Hwa [24] if the particle production process exhibits selfsimilar behavior, the moments show a remarkable power law dependence on phase space bin size. However, if the multiplicity is low, the moments are dominated by statistical fluctuations. In order to suppress the statistical contribution, a modified form of moments in terms of the step function was suggested by Hwa and Pan [25]. Takagi [26] also proposed a new method called Takagi moment method ( moment) for studying the fractal structure of multiparticle production. Both the and the techniques have been applied extensively to analyze several high energy nucleusnucleus collision data [27–35]. Very recently some sophisticated methods have also been applied to study the fractal nature of multiparticle production process [36–45].
3. Entropy and Fractality
In high energy nucleusnucleus collisions, entropy measurement of produced shower particles may provide important information in studying the multiparticle production mechanism [46]. In high energy collisions, entropy is an important parameter and it is regarded as the most significant characteristic of a system having many degrees of freedom [46]. As entropy reflects how the effective degrees of freedom change from hadronic matter at low temperature to the quarkgluon plasma state at high temperature, it is regarded as a useful probe to study the nature of deconfinement phase transition [46]. Entropy plays a key role in the evolution of the high temperature quarkgluon plasma in ultrarelativistic nucleusnucleus interactions in Relativistic HeavyIon Collider (RHIC) experiments and in Large Hadron Collider (LHC) experiments [46]. In nucleusnucleus collisions, entropy measurement can be used to not only search for the formation of QuarkGluonPlasma (QGP) state but it may also serve as an additional tool to investigate the correlations and eventbyevent fluctuations [47]. Different workers have investigated the evolution of entropy in high energy nucleusnucleus collisions at different times [48–64]. The entropy of produced particles can easily be calculated from the multiplicity distribution of the data. If is the probability of producing n particles in a high energy interaction the entropy S is defined by the relation [65]. Mathematically this entropy is called the Shannon entropy [66].
Apart from studying the wellknown Shannon entropy, people are interested to explore the hidden physics of Renyi entropy [66–69] as well, mainly motivated by the inspiration of A. Bialas and W. Czyż [51, 64, 70]. According to C.W Ma and Y.G. Ma [71] the difference between qth order Renyi entropy and order Renyi entropy is found just to be a q dependent constant but which is very sensitive to the form of probability distribution. Renyi entropy can play a potential role to investigate the fractal characteristics of multiparticle production process [72, 73]. The advantage of this method to study the fractal properties of multiparticle production process is that it is not related to the width and resolution of the phase space interval [72, 73]. This method can be applied to events having higher as well as lower multiplicity. This method never suffers from the drawback of lower statistics.
In terms of the probability of multiplicity distribution , the q th order Renyi entropy can be defined as [51, 64, 70]If = , then (1) can be written asGeneralized fractal dimension can be evaluated from the concept of Renyi entropy according to the relationwhere is the central rapidity value in the centre of mass frame and is given by Here is the centre of mass energy of the concerned collision process, is the rest mass of pions, and denotes the average number of participating nucleons.
The generalized fractal dimension is related to the anomalous fractal dimension by a simple mathematical relation [74]Goal of our present study is to carry out an investigation of multifractality and multifractal specific heat in shower particle multiplicity distribution from the concept of Renyi entropy measurements in ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c. We have compared our experimental results with the prediction of UltraRelativistic Quantum Molecular Dynamics (UrQMD) model. Importance of this study is that so far very few attempts have been made to explore the presence of multifractality in multiparticle production process in the framework of higherorder Renyi entropy in high energy nucleusnucleus interactions.
4. Experimental Details
In order to collect the data used for the present analysis, stacks of NIKFIBR2 emulsion pellicles of dimensions 20cm10 cm600 were irradiated by the ^{16}O, ^{28}Si, and ^{32}S beam at 4.5 AGeV/c obtained from the Synchrophasotron at the Joint Institute of Nuclear Research (JINR), Dubna, Russia [75–78]. According to Powell [79], in nuclear emulsion detector particles emitted and produced from an interaction are classified into four categories, namely, the shower particles, the grey particles, the black particles, and the projectile fragments. Details of scanning and measurement procedure of our study along with the characteristics of these emitted and produced particles in nuclear emulsion can be found from our earlier publications [75–78].
One of the problem encountered in interpreting results from nuclear emulsion is the nonhomogeneous composition of emulsion which contain both light (H,C,N, and O) and heavy target nuclei (Ag, Br). In emulsion experiments it is very difficult to identify the exact target nucleus [75]. Based on the number of heavy tracks () total number of inelastic interactions can be divided into three broad target groups H, CNO, and AgBr in nuclear emulsion [75]. Detailed method of target identification has been described in our earlier publication [75]. For the present analysis we have not considered the events which are found to occur due to collisions of the projectile beam with H and CNO target present in nuclear emulsion. Our analysis has been carried out for the interactions of three different projectile ^{16}O, ^{28}Si, and ^{32}S at 4.5 AGeV/c with the AgBr target only. Applying the criteria of selecting AgBr events ( >8) we have chosen 1057 events of ^{16}OAgBr interactions, 514 events of ^{28}SiAgBr, and 434 events of ^{32}SAgBr interactions at 4.5 AGeV/c [75]. Our analysis has been performed on the shower tracks only. We have calculated the average multiplicity of shower tracks in each interaction and presented the values in Table 1 [75–78].

5. Analysis and Results
In a very recent paper [80] we have investigated the Renyi entropy of second order of shower particles using ^{16}O, ^{28}Si, and ^{32}S projectiles on interaction with AgBr and CNO target present in nuclear emulsion at an incident momentum of 4.5 AGeV/c. In this paper we have extended our analysis of Renyi entropy to the study of fractality in multiparticle production of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c.
In order to study the fractal nature of multiparticle production from the concept of Renyi entropy we have calculated the Renyi entropy values of order q=25 from relations (1) and (2) for all the three interactions. The calculated values of Renyi entropy of different orders of the produced shower particles in ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c have been presented in Table 2. It may be mentioned here that the values of the secondorder Renyi entropy have been taken from our recent publication [80]. From the table it can be noted that for all the interactions Renyi entropy values are found to decrease as the order number increases. The variation of Renyi entropy with order q has been presented in Figure 1 for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. Error bars drawn to every experimental point are statistical errors only. Using (3) and (4) we have calculated the values of generalized fractal dimension for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. Calculated values of generalized fractal dimension for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions have been presented in Table 3. From the table it may be noted that the values of generalized fractal dimension decrease with the increase of order number suggesting the presence of multifractality in multipion production mechanism. Presence of multifractality during the production of shower particles indicates the occurrence of cascade mechanism in particle production process. From Table 3 it may also be noted that the values of the generalized fractal dimension remain almost the same within the experimental error for ^{16}OAgBr and ^{28}SiAgBr interactions. But for ^{32}SAgBr interactions the values are higher in comparison to the other two interactions. The variation of against the order q has been shown in Figure 2 for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions.


Now it will be interesting to see how the analysis will look like if one uses Shannon entropy instead of Renyi entropy. However, Shannon entropy cannot be calculated for different orders and hence we can only calculate the values of information fractal dimension for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. We have calculated the values of Shannon entropies derived from the concept of GibbsBoltzmann theories of entropy and tabulated the values in Table 4. The values of information fractal dimension for the three interactions have also been calculated and presented in the same table. Comparing Tables 2 and 4 it may be noticed that the values of Shannon entropy for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions are little higher than those of Renyi entropies.

From the values of the generalized fractal dimension calculated from Renyi entropy values we have evaluated the values of anomalous fractal dimension and hence the ratio of has been calculated. In Table 5 the calculated values of have been presented for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. It is worthwhile to point out that the spiky structure of density distribution of shower particles can also be investigated with the help of a set of bunching parameters [81]. The higherorder bunching parameters can be expressed in terms of lowerorder parameters resulting in a linear expression for the anomalous fractal dimension [82].so thatA nonzero value of the implies the multifractal behavior [82]. We have applied this theory to our study in order to quantify the fractal nature of shower particle production. We have plotted the variation of against the order number in Figure 3 for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. The calculated values of the slope parameter have been presented in Table 6 for our data. The value of signifies the degree of multifractality. From the table it may be seen that for all the interactions the value is greater than zero. This reconfirms the multifractal nature of multiparticle production mechanism. Moreover, value characterizing the degree of multifractality is found to depend on the mass number of the projectile beam. Degree of multifractality is found to increase with the increase of projectile mass as evident from Table 6.


R.C Hwa suggested that [83] from the concept of multifractality a qualitative and quantitative investigation of quarkhadron phase transition in high energy nucleusnucleus collisions is possible. In analogy with the photo count problem at the onset of lasing in nonlinear optics, the coherent state description in high energy nucleusnucleus interactions can be used in the frame work of GinzburgLandau theory [83, 84]. A quantity in terms of the ratio of higherorder anomalous fractal dimension to the secondorder anomalous fractal dimension can be defined by the following relation [83, 84]:According to GinzburgLandau model [84] is related to (q1) by the relationRelation (8) and relation (9) are found to be valid for all systems which can be described by the GinzburgLandau (GL) theory and also are independent of the underlying dimension of the parameters of the model [84]. If the value of the scaling exponent is equal to or close enough (within the experimental error) to 1.304 then a quarkhadron phase transition is expected for the experimental data [84]. If the measured value of is different from the critical value 1.304 considering the experimental errors then the possibility of the quarkhadron phase transition has to be ruled out [84]. is a universal quantity valid for all systems describable by the GinzburgLandau (GL) theory. It is independent of the underlying dimension or the parameters of the model. The critical exponent is an important parameter to investigate the possibility of quarkhadron phase transition, since neither the transition temperature nor the other important parameters are known there [3, 84]. If the signature of quarkhadron phase transition depends on the details of the heavyion collisions, e.g., nuclear sizes, collision energy, transverse energy, etc., even after the system has passed the thresholds for the creation of quarkgluon plasma, such a signature is likely to be sensitive to the theoretical model used [3, 84]. But the critical exponent ν is independent of such details. The value of the critical exponent depends only on the validity of the GinzburgLandau (GL) description of the phase transition for the problem concerned. Here lies the importance of this critical exponent.
We have calculated the values of using relation (8) and presented the values in Table 5 for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. In order to search for the GinzburgLandau secondorder phase transition we have studied the variation of with (q1) in Figures 4, 5, and 6 in case of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5 AGeV/c for the experimental data. The variations of with (q1) have been fitted with the function in order to extract the critical exponent . In Table 7 we have shown the calculated values of the critical exponent for^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions for experimental events. Table 7 reflects that the calculated values of critical exponents obtained from our analysis increase with the increase of projectile mass. The experimental values of the critical exponent for ^{16}OAgBr and ^{28}SiAgBr interactions are found to be lower than the critical value 1.304 while for ^{32}SAgBr interaction the critical exponent is higher than the critical value signifying the absence of quarkhadron phase transition in our data.

Interpretation of multifractality from the thermodynamical point of view allows us to study the fractal properties of stochastic processes with the help of standard concept of thermodynamics. In thermodynamics the constantspecificheat approximation is widely applicable in many important cases; for example, the specific heat of gases and solids is constant, independent of temperature over a larger or smaller temperature interval [85]. This approximation is also applicable to multifractal data of multiparticle production processes proposed by Bershadskii [86]. Bershadskii pointed out that [87] regions of the temperature where the constantspecificheat approximation is applicable are usually far away from the phase transition regimes. For the considered interactions the phasetransitionlike phenomena can occur in the vicinity of q=0, where q is the inverse of temperature. Such situation is expected at the onset of chaos of the dynamical attractors [87]. Barshadskii argued that [86] multiparticle production in high energy nucleusnucleus collisions is related to phasetransitionlike phenomena [3]. He introduced a multifractal Bernoulli distribution which appears in a natural way at the morphological phase transition from monofractality to multifractality. Multifractal Bernoulli distribution plays an important role in multiparticle production at higher energies. It has been pointed out that [11] multifractal specific heat can be derived from the relation = + if the monofractal to multifractal transition is governed by the Bernoulli distribution. The slope of the linear best fit curve showing the variation of against has been designated as multifractal specific heat. The gap at the multifractal specific heat at the multifractality to monofractality transition allows us to consider this transition as a thermodynamic phase transition [88, 89].
We have studied the variation of against for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions in Figure 7. The experimental points have been fitted with the best linear behavior and the slope of the best linear behavior reflects the values of multifractal specific heat. The linear behavior in Figure 7 indicates approximately good agreement between the experimental data and the multifractal Bernoulli representation. The calculated values of the multifractal specific heat for all the three interactions have been presented in Table 8. From Table 8 it may be noted that within the experimental error multifractal specific heat remains almost constant for the three interactions.

The experimental results have been compared with those obtained by analyzing events generated by the UltraRelativistic Quantum Molecular Dynamics (UrQMD) model. UrQMD is a hadronic transport model and this model can be used in the entire available range of energies to simulate nucleusnucleus interactions. For more details about this model, readers are requested to consult [81, 90]. In our earlier papers [78, 80] we have utilized the UrQMD model to simulate ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5 AGeV/c. As described in our previous papers [78, 80] we have generated a large sample of events using the UrQMD code (UrQMD 3.3p1) for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c. We have also calculated the average multiplicities of the shower tracks for all the three interactions in case of the UrQMD data sample [78]. Average multiplicities of the shower tracks in case of UrQMD data sample have been presented in Table 1 along with the average multiplicity values of shower particles in the case of the experimental events. Table 1 shows that the average multiplicities of the shower tracks for the UrQMD events are comparable with those of the experimental values for all the interactions [78]. We have calculated the values of Renyi entropy of order q=25 for all the three interactions using the UrQMD simulated data. The calculated values of Renyi entropy have been presented in Table 2 along with the experimental values. For the UrQMD simulated data also secondorder Renyi entropy values have been taken from our recent publication [80]. From the table it may be seen that the experimentally calculated values of Renyi entropy are higher than those of their UrQMD counterparts signifies the presence of more disorderness for the experimental data. The variation of with order q for the UrQMD simulated data has been presented in Figure 1 along with the experimental plots for each interaction. For the UrQMD simulated events we have calculated the values of the generalized fractal dimension from the values of Renyi entropy using relations (3) and (4). The calculated values of generalized fractal dimension for the UrQMD simulated events of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c have been presented in Table 3 along with the corresponding experimental values. From the table it is seen that for the UrQMD simulated events also the generalized fractal dimension decreases with order q signifying the presence of multifractality for the simulated events. But the values of for the simulated events are lower than the corresponding experimental counterparts for all the interactions. The variation of with order q has been shown in Figure 2 for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions along with the experimental plots. We have calculated the values of Shannon entropies and information dimension for the UrQMD data set of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. The calculated values of Shannon entropy and information dimension have been presented in Table 4. From Table 4 it may be noticed that the Shannon entropy and information dimension for the UrQMD simulated events in case of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions are lower than the corresponding experimental values.
As in the case of experimental data in case of UrQMD simulated data also we have calculated the values of and presented the values in Table 5. From the table it may be noticed that the experimentally obtained values of are little higher than the corresponding UrQMD simulated values. To quantify the multifractality in case of UrQMD simulated data we have studied the variation of with order q in Figure 3 for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. From the slope of the plot of against q we have calculated the value of which signifies the degree of multifractality. The extracted value has been presented in Table 6 for the UrQMD simulated events. From the table it may be seen that the value characterizing the degree of multifractality calculated for the UrQMD simulated events is significantly lower than the corresponding experimental value for all the three interactions. This observation signifies the presence of stronger multifractality for the experimental data. We have also studied the GinzburgLandau phase transition with UrQMD simulated events of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions. The variations of with (q1) for the UrQMD simulated data have been presented in Figures 4, 5, and 6 in case of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions, respectively, at 4.5 AGeV/c along with the experimental plot. The variations of with (q1) have been fitted with the function in order to extract the critical exponent . In Table 7 we have shown the calculated values of the critical exponents for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions for the experimentally simulated events. From Table 7 it may be seen that the critical exponent values for the simulated events are lower than those of the experimental events and there is no evidence of quarkhadron phase transition for the simulated data also. The values of the critical exponent calculated from the multifractal analysis in case of UrQMD simulated data remain almost constant with respect to the mass number of the projectile beam.
We have studied the variation of against in Figure 7 for the UrQMD data sample of ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions in order to calculate the multifractal specific heat. From the slope of the best linear behavior of the plotted points the multifractal specific heat of the shower particles for the UrQMD data sample has been evaluated and presented in Table 8 for all the three interactions. From the table it may be seen that the values of the multifractal specific heat for the experimental and simulated data agree reasonably well with each other. This observation signifies the constancy of multifractal specific heat for the UrQMD simulated events also.
6. Conclusions and Outlook
In this paper we have presented an analysis of multifractality and multifractal specific heat in the frame work of Renyi entropy analysis for the produced shower particles in nuclear emulsion detector for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions at 4.5AGeV/c. Experimental results have been compared with the prediction of UltraRelativistic Quantum Molecular Dynamics (UrQMD) model. Qualitative information about the multifractal dynamics of particle production process has been extracted and reported in the present analysis. The significant conclusions of this analysis are as follows:(1)Renyi entropy values of all the interactions decrease with the order number q for both experimental and UrQMD simulated data. Renyi entropy values for UrQMD data are less than the corresponding experimental values.(2)Generalized fractal dimension calculated from Renyi entropy for both experimental and UrQMD simulated data decreases with the increase of order q suggesting the presence of multifractality in multiparticle production process.(3)The values of Shannon entropy for ^{16}OAgBr, ^{28}SiAgBr, and ^{32}SAgBr interactions are little higher than those of Renyi entropies.(4)Degree of multifractality is found to be higher for the experimental data in comparison to the simulated data and it increases with the increase of projectile mass for the experimental data.(5)The experimental values of the critical exponents for ^{16}OAgBr and ^{28}SiAgBr interactions are lower than the critical value 1.304 required for a quarkhadron phase transition to occur while for ^{32}SAgBr interaction the experimentally obtained values of critical exponent are higher than the critical value 1.304 signifying the absence of quarkhadron phase transition. Absence of quarkhadron phase transition is prominent for the simulated events also.(6)The calculated values of critical exponents obtained from our analysis increase with the increase of projectile mass for the experimental data. UrQMD predicted values of the critical exponent remain almost constant with the increase of projectile mass.(7)Multifractal specific heat for the simulated data agrees well with the experimental data. Constancy of multifractal specific heat is reflected from our analysis.
It is true that there are many papers available in the literature where presence of multifractality has been tested experimentally in multiparticle production in high energy nucleusnucleus interactions by different methods. But the method adopted in this paper to study multifractality seems to be simple and interesting and in this regard our study deserves attention. The observed multifractal behavior of the produced shower particles may be viewed as an experimental fact.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
There are no conflicts of interest in publishing the paper
Acknowledgments
The authors are grateful to Professor Pavel Zarubin of the Joint Institute of Nuclear Research (JINR), Dubna, Russia, for providing them the required emulsion data. Dr. Bhattacharyya acknowledges Professor Dipak Ghosh, Department of Physics, Jadavpur University, and Professor Argha Deb, Department of Physics, Jadavpur University, for their inspiration in the preparation of this manuscript.
References
 A. Bialas and R. Peschanski, “Moments of rapidity distributions as a measure of shortrange fluctuations in highenergy collisions,” Nuclear Physics B, vol. 273, no. 34, pp. 703–718, 1986. View at: Publisher Site  Google Scholar
 A. Bialas and R. Peschanski, “Intermittency in multiparticle production at high energy,” Nuclear Physics B, vol. 308, no. 4, pp. 857–867, 1988. View at: Publisher Site  Google Scholar
 E. A. De Wolf, I. M. Dremin, and W. Kittel, “Scaling laws for density correlations and fluctuations in multiparticle dynamics,” Physics Reports, vol. 270, no. 12, pp. 1–141, 1996. View at: Google Scholar
 B. B. Mandelbrot, The Fractal Geometry of Nature, Henry Holt and Company, San Francisco, Calif, USA, 1982. View at: MathSciNet
 Z. Chen, P. Ivanov, k. Hu, and H. E. Stanley, “Effect of nonstationarities on detrended fluctuation analysis,” Physical Review E, vol. 65, Article ID 041107, 2002. View at: Google Scholar
 D. S. Harte, Multifractals: Theory and Applications, Chapman & Hall/crc, 2001.
 K. J. Falconer, The Geometry of Fractal Sets, vol. 85 of Cambridge Tracts in Mathematics, Cambridge University Press, Cambridge, UK, 1985. View at: MathSciNet
 T. C. Halsey, P. Meakin, and I. Procaccia, “Scaling structure of the surface layer of diffusionlimited aggregates,” Physical Review Letters, vol. 56, no. 8, pp. 854–857, 1986. View at: Publisher Site  Google Scholar
 W. G. Hanan and D. M. Heffernan, “Multifractal analysis of the branch structure of diffusionlimited aggregates,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 85, no. 2, 2012. View at: Publisher Site  Google Scholar
 W. G. Hanan and D. M. Heffernan, “Global structure and finitesize effects in the f(α) of diffusionlimited aggregates,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 77, no. 1, 2008. View at: Publisher Site  Google Scholar
 H. E. Stanley and P. Meakin, “Multifractal phenomena in physics and chemistry,” Nature, vol. 335, no. 6189, pp. 405–409, 1988. View at: Publisher Site  Google Scholar
 A. Cummings, G. O'Sullivan, W. G. Hanan, and D. M. Heffernan, “Multifractal analysis of selected rareearth elements,” Journal of Physics B: Atomic, Molecular and Optical Physics, vol. 34, 2001. View at: Publisher Site  Google Scholar
 W. G. Hanan, D. M. Heffernan, and J. C. Earnshaw, “Leftsided Multifractality of the Harmonic Measure on 2D Cluster  Cluster Aggregates,” Chaos, Solitons & Fractals, vol. 9, no. 6, pp. 875–880, 1998. View at: Publisher Site  Google Scholar
 B. B. Mandelbrot, “Intermittent turbulence in selfsimilar cascades: divergence of high moments and dimension of the carrier,” Journal of Fluid Mechanics, vol. 62, pp. 331–358, 1974. View at: Google Scholar
 P. C. Ivanov, L. A. Nunes Amaral, A. L. Goldberger et al., “Multifractality in human heartbeat dynamics,” Nature, vol. 399, no. 6735, pp. 461–465, 1999. View at: Publisher Site  Google Scholar
 Y. Zheng, J. Gao, J. C. Sanchez, J. C. Principe, and M. S. Okun, “Multiplicative multifractal modeling and discrimination of human neuronal activity,” Physics Letters A, vol. 344, no. 24, pp. 253–264, 2005. View at: Google Scholar
 E. A. Ihlen, “Introduction to multifractal detrended fluctuation analysis in matlab,” Front. Physiol, vol. 3, p. 141, 2012. View at: Google Scholar
 A. Feldmann, A. C. Gilbert, and W. Willinger, “Data networks as cascades: investigating the multifractal nature of Internet WAN traffic,” in Proceedings of the the ACM SIGCOMM '98 conference on Applications, technologies, architectures, and protocols for computer communication, pp. 42–55, Vancouver, British Columbia, Canada, August 1998. View at: Publisher Site  Google Scholar
 N. Sala, “Fractal Geometry in the Arts: An Overview Across the Different Cultures,” in Thinking in Patterns: Fractals and Related Phenomena in Nature, M. M. Novak, Ed., World Scientific, Singapore, 2004. View at: Google Scholar
 B. B. Mandelbrot and R. L. Hudson, The (Mis)Behavior of Markets: A Fractal View of Risk, Ruin, and Reward, Basic Books, New York, NY, USA, 2004. View at: MathSciNet
 Ł. Czarnecki and D. Grech, “Multifractal dynamics of stock markets,” Acta Physica Polonica A, vol. 117, no. 4, pp. 623–629, 2010. View at: Publisher Site  Google Scholar
 A. Turiel and C. J. PérezVicente, “Multifractal geometry in stock market time series,” Physica A: Statistical Mechanics and its Applications, vol. 322, no. 14, pp. 629–649, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 A. Bialas and R. C. Hwa, “Intermittency parameters as a possible signal for quarkgluon plasma formation,” Physics Letters B, vol. 253, no. 34, pp. 436–438, 1991. View at: Publisher Site  Google Scholar
 R. C. Hwa, “Fractal measures in multiparticle production,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 41, no. 5, pp. 1456–1462, 1990. View at: Publisher Site  Google Scholar
 R. C. Hwa and J. Pan, “Fractal behavior of multiplicity fluctuations in highenergy collisions,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 45, no. 5, pp. 1476–1483, 1992. View at: Publisher Site  Google Scholar
 F. Takagi, “Multifractal structure of multiplicity distributions in particle collisions at high energies,” Physical Review Letters, vol. 72, no. 1, pp. 32–35, 1994. View at: Publisher Site  Google Scholar
 R. Saha, A. Deb, and D. Ghosh, “Study of phase transition and its dependence on target excitation using fractal properties of pionization process,” Canadian Journal of Physics, vol. 95, no. 8, pp. 699–710, 2017. View at: Publisher Site  Google Scholar
 N. Parashar, “Nuclear effects and intermittency in two dimensions,” Journal of Physics G: Nuclear and Particle Physics, vol. 22, p. 59, 1996. View at: Google Scholar
 M. K. Ghosh, A. Mukhopadhyay, and G. Singh, “Intermittency and multiplicity moments of charged particles produced in ^{32}S–Ag/Br interaction at 200 A GeV/c,” Journal of Physics G: Nuclear and Particle Physics, vol. 32, p. 2293, 2006. View at: Google Scholar
 S. Ahmad and M. A. Ahmed, “A comparative study of multifractal moments in relativistic heavyion collisions,” Journal of Physics G: Nuclear and Particle Physics, vol. 32, p. 1279, 2006. View at: Google Scholar
 S. Bhattacharyya, A. T. Neagu, and M. Haiduc, “Multifractality analysis of crack images from indirect thermal drying of thinfilm dewatered sludge,” Physica A: Statistical Mechanics and its Applications, vol. 390, no. 14, pp. 2678–2685, 2011. View at: Publisher Site  Google Scholar
 S. Ahmad, A. R. Khan, M. Zafar, and M. Irfan, “On multifractality and multifractal specific heat in ionion collisions,” Chaos, Solitons & Fractals, vol. 42, no. 1, pp. 538–547, 2009. View at: Publisher Site  Google Scholar
 D. Ghosh, A. Deb, and M. B. Lahiri, “Pion production and multifractality in ultrarelativistic nuclear collision – evidence of constant specific heat,” Canadian Journal of Physics, vol. 83, pp. 1169–1175, 2005. View at: Google Scholar
 A. Kamal, S. Ahmad, and N. Ahmad, “Evidence of deconfined phase transition in multiparticle production in relativistic nuclear collisions,” International Journal of Modern Physics A, vol. 23, Article ID 1450016, 2014. View at: Google Scholar
 S. Bhattacharyya, M. Haiduc, A. T. Neagu, and E. Firu Fractals, “A study of multifractality and phase transition in heavyion collisions — experimental data versus model simulation,” Fractals, vol. 26, Article ID 1850015, 2018. View at: Google Scholar
 S. Dutta, D. Ghosh, and S. Chatterjee, “Multifractal detrended fluctuation analysis of pseudorapidity and azimuthal distribution of pions emitted in high energy nuclear collisions,” International Journal of Modern Physics A, vol. 29, Article ID 1450084, 2014. View at: Google Scholar
 Y. X. Zhang, W. Y. Qian, and C. B. Yang, “Multifractal structure of pseudorapidity and azimuthal distributions of the shower particles in Au + Au collisions at 200 A geV,” International Journal of Modern Physics A, vol. 23, no. 18, pp. 2809–2816, 2008. View at: Publisher Site  Google Scholar
 S. Bhaduri and D. Ghosh, “Fractal study of pion void probability distribution in ultrarelativistic nuclear collision and its target dependence,” Modern Physics Letters A, vol. 31, no. 27, Article ID 1650152, 2016. View at: Publisher Site  Google Scholar
 S. Bhaduri and D. Ghosh, “A new approach of chaos and complex network method to study fluctuation and phase transition in nuclear collision at high energy,” The European Physical Journal A, vol. 53, p. 135, 2017. View at: Publisher Site  Google Scholar
 G. Bhoumik, “Comparative Multifractal Detrended Fluctuation Analysis of Heavy Ion Interactions at a Few GeV to a Few Hundred GeV,” Advances in High Energy Physics, vol. 2016, Article ID 7287803, 9 pages, 2016. View at: Google Scholar
 P. Mali, A. Mukhopadhyay, and G. Singh, “Multifractal detrended moving average analysis of particle density functions in relativistic nuclear collisions,” Physica A, vol. 450, pp. 323–332, 2016. View at: Publisher Site  Google Scholar
 P. Mali, “Multifractal analysis of multiparticle emission data in the framework of visibility graph and sandbox algorithm,” Physica A: Statistical Mechanics and its Applications, vol. 493, pp. 253–266, 2018. View at: Google Scholar
 P. Mali, A. Mukhopadhyay, and G. Singh, “Detrended analysis of shower track distribution in nucleusnucleus interactions at CERN SPS energ,” Chaos Solitons & Fractals, vol. 94, p. 86, 2017. View at: Google Scholar
 P. Mali, A. Mukhopadhyay, S. K. Manna, P. K. Haldar, and G. Singh, “Multifractal analysis of charged particle distributions using horizontal visibility graph and sandbox algorithm,” Modern Physics Letters A, vol. 32, Article ID 1750024, 2017. View at: Google Scholar
 P. Mali, S. Sarkar, S. Ghosh, A. Mukhopadhyay, and G. Singh, “Multifractal detrended fluctuation analysis of particle density fluctuations in highenergy nuclear collisions,” Physica A: Statistical Mechanics and its Applications, vol. 424, pp. 25–33, 2015. View at: Publisher Site  Google Scholar
 K. Yamazaki, T. Matsui, and g. Baym, “Entropy in the quark–hadron transition,” Nuclear Physics A, vol. 933, pp. 245–255, 2015. View at: Google Scholar
 S. Ahmad, A. Ahmad, A. Chandra, M. Zafar, and M. Irfan, “Entropy analysis in relativistic heavyion collisions,” Advances in High Energy Physics, vol. 2013, Article ID 836071, 10 pages, 2013. View at: Publisher Site  Google Scholar
 S. Pal and S. Pratt, “Entropy production at RHIC,” Physics Letters, vol. 578, no. 3, pp. 310–317, 2004. View at: Google Scholar
 A. Bialas, W. Czyz, and J. Wosiek, “Studying thermodynamics in heavy ion collisions,” Acta Phys. Pol .B, vol. 30, pp. 107–118, 1999. View at: Google Scholar
 K. Fialkowski and R. Wit, “Entropy in hadronhadron collisions,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 62, no. 11, 2000. View at: Publisher Site  Google Scholar
 A. Bialas and W. Czyz, “Event by event analysis and entropy of multiparticle systems,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 61, no. 7, 2000. View at: Publisher Site  Google Scholar
 M. Y. Sinyukov and S. V. Akkelin Acta, “Pion entropy and phasespace densities in A + A collisions,” Acta Physica Hungarica Series A, vol. 22, no. 12, pp. 171–178, 2005. View at: Google Scholar
 D. Ghosh, A. Deb, and S. Bhattacharyya, “Entropy scaling in highenergy nucleus–nucleus interactions,” Journal of Physics G: Nuclear and Particle Physics, vol. 22, pp. 263–1114, 2005. View at: Publisher Site  Google Scholar
 M. Reiter, A. Dumitru, J. Brachmann, J. A. Maruhn, H. Stöcker, and W. Greiner, “Entropy production in collisions of relativistic heavy ions — a signal for quarkgluon plasma phase transition?” Nuclear Physics A, vol. 643, pp. 99–112, 1998. View at: Publisher Site  Google Scholar
 D. Ghosh, A. Deb, P. K. Haldar, S. R. Sahoo, and D. Maity, “Validity of the negative binomial multiplicity distribution in case of ultrarelativistic nucleusnucleus interaction in different azimuthal bins,” Europhysics Letters, vol. 65, p. 311, 2004. View at: Google Scholar
 Y. K. Vermani and R. K. Puri, “Entropy and light cluster production in heavyion collisions at intermediate energies,” Nuclear Physics A, vol. 847, no. 34, pp. 243–252, 2010. View at: Google Scholar
 M. R. Atayan, “Entropy Analysis in ^{+}p and K^{+}p Collisions at √s = 22 GeV,” Acta physica polonica B, vol. 36, p. 2969, 2005. View at: Google Scholar
 K. Eskola, H. Paukkunena, and C. Salgado, “From 0 to 5000 in 2 × 10^{−2}4 seconds: Entropy production in relativistic heavyion collisions,” Nuclear Physics A, vol. 830, no. 14, pp. 519c–522c, 2009. View at: Google Scholar
 Y. G. Ma, “Application of Information Theory in Nuclear Liquid Gas Phase Transition,” Physical Review Letters, vol. 83, no. 18, pp. 3617–3620, 1999. View at: Publisher Site  Google Scholar
 D. Ghosh, A. Deb, S. Bhattacharyya, and U. Datta, “Multiplicity dependence of entropy in different rapidity bins in highenergy nucleus–nucleus interactions,” Physica Scripta, vol. 85, Article ID 065205, 2012. View at: Google Scholar
 D. Ghosh, A. Deb, and S. Bhattacharyya, “Entropy scaling in highenergy nucleus–nucleus interactions,” Journal of Physics G: Nuclear and Particle Physics, vol. 38, Article ID 065105, 2011. View at: Google Scholar
 S. Ahmad, A. Chandra, M. Zafar, and M. Irfan, “Entropy scaling in ionion collisions at AGS and SPS energies,” International Journal of Modern Physics E, vol. 22, Article ID 1350088, 2013. View at: Google Scholar
 M. Mohisin Khan, N. Ahmad, and M. Irfan, “Evidence of Entropy Scaling in Relativistic Nucleus–Nucleus Collisions,” International Journal of Modern Physics E, vol. 21, Article ID 1250068, 2012. View at: Google Scholar
 A. Bialas, W. Czyż, and A. Białas, “Renyi Entropies in Multiparticle Production,” Acta Physica Polonica B, vol. 31, p. 2803, 2000. View at: Google Scholar
 V. Simak, “Entropy in Multiparticle Production and Ultimate Multiplicity Scaling,” Physics Letters B, vol. 206, pp. 159–162, 1988. View at: Google Scholar
 C. E. Shannon, The Mathematical Theory of Communication, The University of Illinois Press, Urbana, Ill, USA, 1948. View at: MathSciNet
 A. Rényi, “On Measurements of Entropy and Information,” in Proceedings of the 4th Berkeley Symp on Math. Stat. Prob, pp. 547–561, University of California Press. View at: Google Scholar  MathSciNet
 C. Beck and F. Schloegl, Thermodynamics of chaotic systems, vol. 4 of Cambridge Nonlinear Science Series, Cambridge University Press, Cambridge, UK, 1993. View at: Publisher Site  MathSciNet
 I. M. Dremin and J. W. Gary, “Hadron multiplicities,” Physics Reports, vol. 349, pp. 301–393, 2001. View at: Google Scholar
 A. Bialas and W. Czyz, “Measurement of entropy of a multiparticle system: a dolist,” High Energy Physics  Phenomenology, vol. 31, p. 687, 2000. View at: Google Scholar
 C. W. Ma and Y. G. Ma, “Shannon information entropy in heavyion collisions,” Progress in Particle and Nuclear Physics, vol. 99, pp. 120–158, 2018, rogress in Particle and Nuclear Physics. View at: Google Scholar
 M. Pachr, V. Simak, M. Sumbera, and I. Zborovsky, “Entropy, dimensions and other multifractal characteristics of multiplicity distributions,” Modern Physics Letters A, vol. 07, no. 25, pp. 2333–2339, 1992. View at: Publisher Site  Google Scholar
 A. Mukhopadhyay, P. L. Jain, and G. Singh, “Entropy and fractal characteristics of multiparticle production at relativistic heavy ion interactions,” Physical Review C, vol. 47, p. 410, 1993. View at: Publisher Site  Google Scholar
 W. Ochs and J. Wosiek, “Multiparton correlations and fractal structure in QCD jets,” Physics Letters B, vol. 305, no. 12, pp. 144–150, 1993. View at: Publisher Site  Google Scholar
 S. Bhattacharyya, M. Haiduc, A. T. Neagu, and E. Firu, “Target dependence of clan model parameters at Dubna energychaotic pion production,” Journal of Physics G: Nuclear and Particle Physics, vol. 40, Article ID 025105, 2013. View at: Google Scholar
 S. Bhattacharyya, M. Haiduc, A. T. Neagu, and E. Firu, “Eventbyevent pseudorapidity fluctuations in high energy nucleus–nucleus interactions,” Physics Letters B, vol. 726, no. 13, pp. 194–205, 2013. View at: Publisher Site  Google Scholar
 S. Bhattacharyya, M. Haiduc, A. T. Neagu, and E. Firu, “Forward–backward multiplicity correlation in highenergy nucleus–nucleus interactions at a few AGeV/c,” Journal of Physics G: Nuclear and Particle Physics, vol. 41, Article ID 075106, 2014. View at: Google Scholar
 S. Bhattacharyya, M. Haiduc, A. T. Neagu, and E. Firu Euro, “Multiplicity distribution of shower particles in nucleusnucleus collisions at 4.14.5 A GeV/c,” The European Physical Journal Plus, vol. 132, p. 229, 2017. View at: Google Scholar
 C. F. Powell, P. H. Fowler, and D. H. Perkins, The study of elementary particles by photographic method, Oxford, 1959.
 S. Bhattacharyya, M. Haiduc, A. T. Neagu, and E. Firu, “An investigation of Renyi entropy in highenergy nucleus–nucleus collisions,” Canadian Journal of Physics, vol. 95, p. 715, 2017. View at: Publisher Site  Google Scholar
 M. Bleicherdag, E. Zabrodin, C. Spieles, and S. A. Bass, “Relativistic hadronhadron collisions in the ultrarelativistic quantum molecular dynamics model,” Journal of Physics G: Nuclear and Particle Physics, vol. 25, no. 9, p. 1859, 1999. View at: Google Scholar
 S. V. Chekanov and V. I. Kuvshinov Acta, “Bunching Parameter and Intermittency in HighEnergy Collisions,” Acta Physica Polonica B, vol. 25, p. 1189, 1994. View at: Google Scholar
 R. C. Hwa, “Scaling exponent of multiplicity fluctuation in phase transition,” Physical Review D, vol. 47, p. 2773, 1993. View at: Google Scholar
 R. C. Hwa and M. T. Nazirov, “Intermittency in secondorder phase transitions,” Physical Review Letters, vol. 69, no. 5, pp. 741–744, 1992. View at: Publisher Site  Google Scholar
 L. D. Landau and E. M. Lifshitz, Statistical Physics, Pergamon Press, 1980.
 A. Bershadskii, “Anomalous distributions in heavy ion collisions at high energies,” Physical Review C: Nuclear Physics, vol. 59, no. 1, pp. 364–368, 1999. View at: Publisher Site  Google Scholar
 A. Bershadskii, “Multifractal statistics of multiparticle production at high energies,” Journal of Physics A, vol. 2, p. 225, 1998. View at: Google Scholar
 A. Bershadskii, “Multifractal critical thermodynamics in nucleusnucleus collisions,” Journal of Physics A, vol. 4, p. 283, 1999. View at: Google Scholar
 R. Peschanski, “Intermittency in Particle Collisions,” International Journal of Modern Physics A, vol. 6, no. 18, p. 325681, 1991. View at: Publisher Site  Google Scholar
 S. A. Bass, M. Belkacem, M. Bleicher, and M. Brandstetter, “Microscopic Models for Ultrarelativistic Heavy Ion Collisions,” Nuclear Theory, vol. 41, p. 225, 1998. View at: Google Scholar
Copyright
Copyright © 2018 Swarnapratim Bhattacharyya et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.