Abstract

We have studied the multifractality of pion emission process in 16O-AgBr interactions at 2.1 AGeV  and  60 AGeV, 12C-AgBr  and  24Mg-AgBr interactions at 4.5 AGeV, and 32S-AgBr interactions at 200 AGeV using Multifractal Detrended Fluctuation Analysis (MFDFA) method which is capable of extracting the actual multifractal property filtering out the average trend of fluctuation. The analysis reveals that the pseudorapidity distribution of the shower particles is multifractal in nature for all the interactions; that is, pion production mechanism has inbuilt multiscale self-similarity property. We have employed MFDFA method for randomly generated events for 32S-AgBr interactions at 200 AGeV. Comparison of expt. results with those obtained from randomly generated data set reveals that the source of multifractality in our data is the presence of long range correlation. Comparing the results obtained from different interactions, it may be concluded that strength of multifractality decreases with projectile mass for the same projectile energy and for a particular projectile it increases with energy. The values of ordinary Hurst exponent suggest that there is long range correlation present in our data for all the interactions.

1. Introduction

The study of correlation and fractality is an active area of research in many fields including heavy ion collisions [15]. Natural systems, which have irregular pattern at different scales, exhibit fractal nature. Fractals are generally classified into two categories: (i) monofractal and (ii) multifractals. For monofractals scaling properties of the system are identical throughout the system; on the other hand, “multifractals” are more complicated self-similar structure that consist of a number of weighted fractals with different noninteger dimensions. 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 [6]. Investigation of multifractality is of great importance as its origin may be associated with the presence of long range correlation in the system. Correlation study has the potential to provide information about the characteristics of system evolution. Moreover, multifractal analysis is effective in understanding the underlying dynamics of any complex system such as pionisation in high energy nucleus-nucleus interactions. To get both qualitative and quantitative idea concerning the multiparticle production mechanism [7] multifractal analysis is expected to be very fruitful. Such a behavior has been observed for vast majority of high energy multiparticle production experiments [810].

The investigation of fractal dimension in hadronic multiparticle production was carried out probably for the first time by Carruthers and Duong-Van [11]. But there was no formalism developed for a systematic fractal study. A systematic approach for the fractal study was suggested by Hwa ( moment) [12]. But those moments are found to be influenced by statistical fluctuations especially for the low multiplicity events. Later in order to avoid large statistical fluctuations and exclude low multiplicity events Hwa and Pan [13] proposed a modified moment method introducing a step function, which acts as a filter to the low multiplicity events [1416]. Also these modified moments suffer from the demerit that they are defined only for positive orders and hence they are unable to explore the whole multifractal spectrum. Afterwards, a number of techniques [1720] were developed for the fractal study of multiparticle data. The techniques developed by Hwa [12] and Takagi [17] are the most popular and have been used in many cases [2, 2123] to analyze multipion production process. However, none of these methods can disentangle the dynamical “signal” from the “background.” Trending behaviors usually give rise to spurious multifractal effects for the analyzed series. Therefore, it is essential to study the intrinsic fluctuations characterizing the dynamical process after filtering out the average trending behavior.

Sophisticated methods have been invented to characterize the actual fluctuations extracted from the average behavior and the fractal nature of nonstationary time series. These include detrended fluctuation analysis (DFA) and its variance [24, 25], the wavelet transform [26, 27] based multiresolution analysis [28, 29], and Multifractal Detrended Fluctuation Analysis (MFDFA) [30]. DFA technique [24] was developed in order to determine minutely the presence of any long range correlation [24, 31] in a nonstationary series. However, despite a multitude of real-data analyses, a proper detection of the multifractality in the experimental data still presents much difficulty and is not always reliable [32]. MFDFA technique [30] is actually a generalization of standard DFA technique for the characterization of multifractal nature of a series. One main reason to employ MFDFA method is to avoid fallacious detection of correlations leading to multifractality which are artifacts evolving due to the nonstationarity of the signal. Thus, MFDFA is a powerful technique which has been applied successfully to characterize fluctuation in a variety of fields like finance [3336], medicine [37, 38], natural science [39, 40], solid state physics [41, 42], and so forth.

The spectrum of references of application of MFDFA technique is not a complete one. Recently, MFDFA method has been applied to analyze the pseudorapidity and azimuthal angle distribution of the pions produced in Au + Au interactions at 200 GeV/nucleon by Zhang et al. [43]. They studied a sample of only 10 events. Another group, Wang and Yang [44], studied the same interactions for UrQMD generated data using the same method. The DFA and MFDFA methods are also used by Mali et al. [45] to characterize the particle density fluctuation for 28Si + Ag/Br interactions at 14.5 GeV/nucleon and 32S + Ag/Br interactions at 200 GeV/nucleon. These analyses [4345] suggest that the MFDFA approach is a reasonably good technique for the multifractal analysis of multiparticle production process in high energy nucleus-nucleus (A-A) interactions and hence should be applied for understanding the dynamics of the process.

In this paper we have studied the pseudorapidity distribution of the pions produced in 16O-AgBr interactions at 2.1 AGeV and 60 AGeV, 12C-AgBr  and  24Mg-AgBr interactions at 4.5 AGeV, and 32S-AgBr interactions at 200 AGeV in the framework of MFDFA. The analysis is expected to reveal a comparative view of the genuine multifractal parameters for interactions initiated by projectiles of various energies and masses.

2. Experimental Section

The present analysis is performed on the interactions of 16O beam at 2.1 AGeV  and  60 AGeV, 12C beam  and  24Mg both at 4.5 AGeV, and 32S beam at 200 AGeV with AgBr being the target present in nuclear emulsion.

ILFORD G5 nuclear photographic emulsion plates were irradiated horizontally with a beam of 16O nuclei at energy of 2.1 AGeV obtained from BEVALAC Berkley [46]. The data for 12C-AgBr [47] and 24Mg-AgBr [48] interactions were obtained by exposing NIKFI BR2 emulsion plates to the beams of 12C and 24Mg nucleus, of energy 4.5 AGeV at JINR, Dubna, Russia. The data of 16O-AgBr interactions at 60 AGeV and 32S-AgBr interactions at 200 AGeV were obtained by exposing the stacks of ILFORD G5 emulsion plates to a beam of 16O nucleus of energy 60 AGeV and 32S nucleus of energy 200 AGeV, respectively, using the super proton synchrotron (SPS) at CERN [49]. A Leitz Ortholux microscope with a 10x objective and 25x ocular lens provided with a Brower travelling stage was used to scan the plates of 16O-AgBr interactions at 2.1 AGeV. A Leitz Metalloplan microscope with a 10x objective and 10x ocular lens provided with a semiautomatic scanning stage has been used to scan the other four interaction plates. Two observers scanned each plate independently so that the biases in detection, counting, and measurement could be minimized and consequently the scanning efficiency could be increased. After finding a primary interaction induced by the incoming projectile, the number of secondary tracks in an event belonging to each category was counted using oil immersion objectives. Measurements were carried out with the help of an oil immersion objective of 100x magnification. The measuring system fitted with both the microscopes has 1 μm resolution along the -axis and -axis and 0.5 μm resolution along the -axis.

Events were chosen according to the criteria given below.

(a) The incident beam track would have to lie within 3° from main beam direction. (b) Events that occurred within 20 μm from the top and bottom surfaces of the pellicle were rejected. (c) All the primary beam tracks were followed in the backward direction to ensure that the events chosen did not include interactions from the secondary tracks of the other interactions.

According to nuclear emulsion terminology [50], particles emitted after interactions can be classified as the shower, grey, and black particles. Shower particles are mostly (about more than 90%) due to pions with a small admixture of -mesons and hyperons having ionization and velocity greater than 0.7c, where is the minimum ionization of a singly charged particle produced in the emulsion medium. Grey particles are mainly fast target recoil protons with energies up to 400 MeV. They have ionization . They have velocities lying between 0.3c and 0.7c. Black particles consist of both singly and multiply charged fragments of the target nucleus with ionization and velocity less than 0.3c.

Along with the above stated three kinds of particle, there could also be a few projectile fragments. These projectile fragments are the spectator parts of the incident projectile nuclei that do not directly participate in an interaction.

In the experiments with the nuclear emulsion track detectors, interactions may be with three different types of targets, for example, hydrogen (H), light nuclei (CNO), and heavy nuclei (AgBr), present in the emulsion medium. Events, with , occur because of the collision between hydrogen and the projectile beam. Events with are due to collisions of projectile with light nuclei and events with are due to collisions with heavy nuclei. Here , the number of heavy tracks, is the total number of black and grey tracks. In our study, only events having a number of heavy tracks greater than 8 have been selected to exclude the H and CNO events.

According to the selection procedure mentioned above, we have chosen 730 events of 16O-AgBr interactions at 2.1 AGeV [46], 800 events of 12C-AgBr [51] and 24Mg-AgBr [52] interactions at 4.5 AGeV, 250 events of 16O-AgBr interactions at 60 AGeV [53], and 140 events of 32S-AgBr interactions at 200 AGeV [54]. The present analysis has been performed on the pion tracks only. The emission angle was measured for each pion track with respect to the beam direction by taking readings of the coordinates of the interaction point , coordinates of a point at some distance away from the interaction point on each secondary track, and coordinates of a point on the incident beam. From “” the pseudorapidity variable , which may be treated as a convenient substitute of the rapidity variable of a particle when the rest mass of the particle can be neglected in comparison to its energy or momentum, was calculated for each pion track.

3. MFDFA Method

The MFDFA technique was developed by Kantelhardt et al. [30] as a generalization of standard DFA method to analyze nonstationary time series. If is the signal of length , MFDFA consists of the following steps [30] among which the first three steps include the ordinary DFA technique.

Step 1. Calculation of the signal profile, which is the cumulative sum of the signal to be analyzed, is according to

Step 2. Consider division of the profile into numbers of nonoverlapping segments of length . Since is not an integer multiple of , a small part of the signal will be left at the end. In order to include that part, the same process is repeated starting from the other end. Thus, one obtains altogether number of segments.

Step 3. Consider determination of the local trend associated with each of such segments by a least square polynomial fit of the series in a particular segment, as well as calculation of the variance of the series. The variance of the series relative to the local trend in a given segment of length can be calculated aswhere is the fitting polynomial of order for the th segment.

Step 4. th order fluctuation function is then given by [32]

Step 5. To investigate the scaling behavior, one has to calculate for various values. It should be remembered that the sample size of the smallest segment (or scale) should be much larger compared to the polynomial order in order to prevent an overfitted trend.

3.1. Scaling Behavior and Multifractality

Presence of any long range correlation in the system will appear in the multifractal behavior of the system and then will show a power law dependence on scale likeThe scaling exponent is known as generalized Hurst exponent. For it reduces to the ordinary Hurst exponent [55]. indicates the presence of long range correlation in the system. An exponent corresponds to positive correlation in the system and corresponds to the presence of anticorrelation in the system. For a monofractal series is independent of and the fluctuation function is similar for all subsignals. On the other hand, if large and small fluctuations scale differently then there will be a notable dependence of on . Particularly for multifractal series decreases with order ; for +ve values describes the scaling behavior of the segments with large fluctuations and for −ve ’s it describes the scaling behavior of the segments with small fluctuations.

The multifractal exponent can be derived from [30, 55] using (6). ConsiderThe singularity spectrum can be obtained from through a Legendre transform [20, 55, 56]:The width of the multifractal spectrum is the difference between the maximum and minimum value of ; that is,It has been proposed by Ashkenazy et al. [57] and Shimizu et al. [58] that the width of a multifractal spectrum is a measure of the degree of multifractality. The broader the spectrum, the richer the multifractality [59].

Though originally MFDFA technique was developed for nonstationary time series analysis, it can also be used equally for nonuniform distribution like rapidity distribution of high energy multiparticle production process [4345]. In principle, multifractality of a natural system may originate either (a) due to broad probability density function (non-Gaussian distribution) of the concerned parameter or (b) due to the simultaneous presence of dissimilar characteristics of long range correlations for large as well as small fluctuations or sometimes (c) due to both reasons (a) and (b) [60]. MFDFA technique allows one to identify the source of multifractality when it is applied to the randomly shuffled distribution of the same. For shuffled distribution though the probability distribution function remains unchanged, all possible correlations are wiped out from the distribution. Thus, if multifractality is originated from type (b), then randomly generated distribution will show an arbitrary behavior with . On the other hand, if the source of multifractality is type (a), then for the generated distribution generalized Hurst exponent will be identical to that of the original one (). Furthermore, for reason (c) the generated series will show a weaker multifractality than that of the experimental one.

4. Results and Discussions

In the current analysis we have focused on 16O-AgBr interactions at 2.1 AGeV  and  60 AGeV, 12C-AgBr  and  24Mg-AgBr interactions at 4.5 AGeV, and 32S-AgBr interactions at 200 AGeV. For each interaction we have considered only those events which have multiplicities greater than average multiplicity of the corresponding interaction. Pseudorapidity distribution of the pions produced in each such event was obtained for 0.1 rapidity interval. MFDFA technique as discussed has been utilized to analyze the pseudorapidity distribution of the shower particles corresponding to each selected event. Here, the term scale corresponds to the width of each segment in the pseudorapidity space. The smallest scale is so chosen that at least 10 data points are within it; that is, it covers one pseudorapidity unit in the pseudorapidity space. The fluctuation function has been calculated for different “” values with “” varying from −5 to +5 with integer values only using (3) for each selected event. Event average of fluctuation function () is obtained by averaging over all selected events. The plot of with for , 0, and −5 is shown in Figure 1 for all interactions. For any particular order the event average fluctuation function manifests power law dependence on and consequently all the plots are of straight line nature with different slope. The point symbols in Figure 1 represent the values of for different whereas the solid lines represent the regression fit of the same. for different “” values shows a convergent trend as the scale increases. The slope for linear fit of versus graph provides the generalized Hurst exponent . values vary for different ’s.

The variation of with has been depicted in Figure 2. For all the interactions is dependent and it decreases with , signifying multifractal nature of the rapidity distribution as manifested by the different scaling behaviors for large and small fluctuations. values for seem to be higher than those for . This typical feature of is consistent with the observation by Manimaran et al. [29] for multifractal time series. The ordinary Hurst exponent , that is, for , is given in Table 1 for all the interactions. , as mentioned earlier, gives an idea about the correlation present in the system. For all the considered interactions indicates that long range correlation exists in our data for all the interactions.

We have also calculated for different values according to (5). The plot of against is shown in Figure 3. The nonlinear variation of with also confirms the fact that the pion density fluctuation in the rapidity space is multifractal in nature; that is, pion production mechanism has an inbuilt multiscale self-similarity property.

In order to shed light on the possible origin of the demonstrated multifractal behavior we have randomly shuffled the rapidity distribution of each selected event for all interactions and have employed MFDFA method to the shuffled data using the same approach as before. Results were found to be identical for all the interactions. Here we present the result of comparison of experimental and shuffled data only for 32S-AgBr interactions at 200 AGeV.

Figure 4 shows the plot of with along with versus . values remain around 0.5 (average 0.56) with a minor dependence on with . suggests that multifractality in our data is mainly sourced from the presence of long range correlation in the data. versus plot (Figure 5) for the shuffled and experimental data further strengthens this fact. shows an almost linear behavior instead of nonlinear one as for the experimental data. Similar behaviors were observed for the other interactions also.

To have a quantitative idea about the multifractality we have computed the multifractal spectrum corresponding to each interaction. The multifractal spectrums corresponding to different interactions have been displayed in Figure 6. Figure 7 shows the spectrum for 32S-AgBr interactions for both expt. and shuffled distributions. For shuffled data spectrum has negligible width () with respect to experimental data () and also the spectrum shifted to the left of the experimental spectrum as expected when multifractality is mainly sourced [60] from type (b) cause as discussed earlier.

To compare the strengths of multifractality, we have calculated the width of the multifractal spectrum for each of the interactions. These values have been presented in Table 2. The spectrum width is found to vary for different interactions suggesting that the strength of multifractality depends upon the mass and energy of the projectile beam. The spectrum width values reveal the following notable features of the heavy ion interaction process:(i)Comparison of 16O-AgBr interactions at 2.1 AGeV  and  at 60 AGeV indicates that, for a particular projectile, strength of multifractality increases with the energy of the projectile.(ii)From 24Mg-AgBr interactions  and  12C-AgBr interactions at the same energy 4.5 AGeV, it is evident that the strength of multifractality decreases with projectile mass.(iii)For projectiles with different mass and energies, multifractality decreases if there are simultaneous decrease in energy and increase in mass of the projectile beam (e.g., the pair of interactions (i) 16O-AgBr at 60 AGeV  and  24Mg-AgBr at 4.5 AGeV and (ii) 12C-AgBr at 4.5 AGeV  and  16O-AgBr at 2.1 AGeV). This observation agrees very well with the previous two observations.(iv)If projectile mass and energy simultaneously change in the same manner, that is, either both increase or both decrease, the strength of multifractality is predominantly affected by the mass of the projectile if energy is of the same order. This fact is revealed if we consider the pair of interactions (i) 16O-AgBr at 2.1 AGeV  and  24Mg-AgBr at 4.5 AGeV and (ii) 16O-AgBr at 60 AGeV  and  32S-AgBr at 200 AGeV. On the other hand, if mass difference is little, strength of multifractality is mainly determined by the energy difference of the projectile. This is demonstrated by the interactions 12C-AgBr at 4.5 AGeV  and  16O-AgBr at 60 AGeV pair.(v)Finally, if we go from relativistic energy regime to the ultrarelativistic one, influence of energy is much more prominent compared to mass. This is noticed comparing 16O-AgBr at 2.1 AGeV, 24Mg-AgBr at 4.5 AGeV, and 12C-AgBr at 4.5 AGeV interactions with the 32S-AgBr interactions at 200 AGeV.

5. Conclusions

We have presented a systematic study on pion density fluctuation in pseudorapidity spectra in the framework of sophisticated MFDFA technique using various heavy ion projectiles covering a wide range of energy starting from a few GeV to a few hundred GeV. The main observations of the analysis may be summarized as follows:(i)Both and spectra speak in favor of multifractal pion density fluctuation in pseudorapidity space for all the considered A-A interactions.(ii)Multifractal parameters and and the multifractal spectrum for the shuffled data clearly demonstrate that the observed multifractality of the experimental data does not originate from the trivial broad probability density distribution but occurs due to dynamics.(iii)The study suggests multifractal pion production dynamics in heavy ion interactions.(iv)The strength of multifractality is influenced both by mass and by energy of the projectile beam.(v)For the same projectile beam multifractality increases with projectile energy. On the other hand, it decreases with mass of the projectile beams having the same energy.(vi)When both mass and energy change, effect of mass is dominant if their energy does not differ much. But multifractality is much more influenced by energy of the projectile beam if mass difference is small.(vii)From the value of ordinary Hurst exponent we can conclude that long range correlation is present among the pions for all of the considered heavy ion interactions irrespective of projectile mass and energy.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors are grateful to Professor P. L. Jain, State University of Buffalo, Buffalo, USA, Professor P. S. Young of Mississippi State University, Mississippi, and Professor K. D. Tolstov, JINR, Dubna, U.S.S.R., for providing them with exposed and developed emulsion plates used for this analysis. One of the authors, Miss Gopa Bhoumik, acknowledges CSIR for a Junior Research Fellowship (Sanction no. 09/096(0802)/2013-EMR-I) during this work.