Journal of Modern Physics
Vol.09 No.05(2018), Article ID:83998,-2 pages

A Study of Multifractal Analysis in 16O-AgBr Collisions at 60A and 200A GeV

Nazeer Ahmad, Tufail Ahmad, Omveer Singh, Shakeel Ahmad

Department of Physics, Aligarh Muslim University, Aligarh, India

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: December 7, 2017; Accepted: April 22, 2018; Published: April 25, 2018


A multifractal analysis to study the multiparticle dynamics in 60A and 200A GeV, 16O-AgBr collisions has been performed in the pseudorapidity phase space. Multifractal moments, Gq, as the function of pseudorapidity bin size for different order of the moments, q, have been calculated. The power-law behaviour has been observed in the considered data sets. The variations of multifractal dimension, Dq, and the multifractal spectral function, f ( α q ) , with order of the moments, q, have been studied thoroughly. Dq is found to decrease with increasing order of moments, q, indicating thereby a self-similar behaviour in the multiparticle production for the considered collisions. We have also found a concave downward curve of multifractal spectral function with maximum at q = 0.


Multifractality, Pseudorapidity, Fluctuation, Generalized Dimension, Scaling

1. Introduction

The study of event-by-event fluctuations in the pseudorapidity windows of decreasing bin width in multiparticle production at high energy has divulged self-similar properties as speculated by Bialas and Peschanski [1] [2] which is known as intermittency. The concept of self-similarity is closely related to the fractal theory which is a natural consequence of the cascading mechanism prevailing in the multiparticle production process. The power law behaviour, scaled factorial moments, indeed implies the existence of some kind of fractal pattern [2] in the dynamics of particles produced in their final state of reaction.

In order to investigate the cascading mechanism [3] of multiparticle production in the framework of the multifractal technique, a formalism has been developed by R.C. Hwa [4] [5] for a systematic study of the fractal properties and it also provides an effective means of defining non-uniform rapidity distribution of produced particles in nuclear collisions. The aim of this technique is to verify the authenticity of basic scaling properties, which occur in multifractal theories, when applied to particle production for the dynamical mechanism responsible for the hadronisation process in nucleus-nucleus (AA) collisions [6] [7] [8] . The explanation of cascading as self-similar process in analogy with geometrical objects such as fractals has been allowed by scaling law [9] .

It is worth-mentioning that the issues of possible signals of Quark Gluon Plasma (QGP) formation in relativistic heavy-ion collisions are still on discussion stage. Fluctuations with different features are visualized to occur in different types of phase transitions. A study of multiplicity fluctuations in the final state of collisions might enable one to find out new observables linked to the characteristics of different types of phase transitions.

Some workers [10] [11] have proposed entropy based tools which are commonly used to describe the dynamics of complex systems. In the last few decades, it has been observed that Tsallis entropy and multifractal technique are found to be useful to characterize on large scale scaling behaviour. Furthermore, the Shannon, Reyni and Tsallis entropies [12] related formalisms are used as multifractal measures. However, the concept of divergence is being employed to quantify the departure between two distributions. It is one of the most used divergences related to the Tsallis entropy.

It is also interesting to mention that A. Deppman [13] describes a system which has a fractal structure in its thermodynamical functions called as thermofractals. The author has mentioned that its thermodynamics is more naturally described by the Tsallis statistics and a relation between fractal dimension and entropy index has been established. Furthermore, according to the Deppman such complex hadronic system exhibits thermofractal structure.

Although many attempts have been made to study the fractal properties, using μ p , p p ¯ and e + e data [14] [15] , using the method of multifractal moments, Gq, however AA collisions are less studied using this method. In the present paper, an attempt is made to investigate various interesting features of multifractality in 16O-AgBr collisions at 60A and 200A GeV. The Hwa’s proposal which is based on the general theory of multifractals [16] [17] allows us to estimate mass exponent, tq, generalized dimension, Dq and multifractal spectral function f ( α q ) for different order of the moments. These parameters serve as the measure of the multifractal structure in the considered data sets.

2. Theory

In order to examine the dependence of multifractal moments, Gq, on pseudorapidity, h, which is defined as; η = ln tan ( θ s / 2 ) , θ s is the space angle of a secondary particle with mean direction of primary, a given pseudorapidity range Δ η = η max η min is divided into M0 bins of width δ η = Δ η / M 0 . A multifractal moment is defined as;

G q = j = 1 M p j q (1)

where, p j = n j / n , such that n ( = n 1 + n 2 + + n M ) . M denotes the number of non-empty bins. q is a real number and may has both positive or negative values. Once Gq is calculated, its average over the entire sample is determined;

G q = 1 N e v 1 N e v G q (2)

where, Nev stands for the total number of events. If there is self-similarity in the production of particles, Gq moments can be written in the form of a power law;

G q = ( δ η ) τ q for δ η 0 (3)

where, tq stands for mass exponents. The linear dependence of l n G q on l n δ η over all the windows is related as;

τ q = lim δ η ( Δ ln G q / Δ ln δ η ) (4)

The multifractal spectrum f ( α q ) is related to the mass exponents tq and calculated from Legendre transform as:

α q = d τ q / d q , (5)

f ( α q ) = q α q τ q . (6)

The properties of f ( α q ) for multifractal behaviour are defined [4] [5] [16] [17] by;

d f ( α q ) d α q = q , d 2 f ( α q ) d α q 2 < 0 (7)

This downward concave form of curve of f ( α q ) has the following characteristics for multifractals;

1) f ( α q ) is downward concave

2) f ( α q = α 0 ) is maximal

3) f ( α q ) < f ( α 0 ) , for q 0

The behaviour of f ( α q ) shows that f ( α q ) = α q = 1 for all q as a special case if there are absolutely no fluctuations. The width of f ( α q ) is a measure of the size of the fluctuations and the value f ( α 0 ) < 1 is a measure of the number of empty bins.

The generalized dimensions may be defined as;

D q = τ q q 1 (8)

Different types of dimensions are given as under:

1) D 0 = f ( α 0 ) ; Capacity dimension: This shows how the data points of multifractal pattern fill the phase space domain

2) D 1 = f ( α 1 ) = α 1 ; Entropy dimension: This is a measure of order-disorder of the data points in the phase space domain under study. Larger the value higher the disorder.

3) D 2 = 2 α 2 f ( α 2 ) ; Correlation dimension: This quantifies the degree of clustering. Larger the value corresponds to higher level clustering.

3. Experimental Details

The data analyzed in the present paper were collected using two emulsion stacks exposed to Oxygen beams at 60A and 200A GeV at CERN, SPS (EMU01 Collaboration) [7] [18] [19] [20] . Along-the-track scanning method has been used to record the interactions in emulsion because of inherent high detection efficiency [7] [18] [21] [22] . In each event the polar and azimuthal angles of all the particles emitted in the interactions were measured with high magnification microscope. Tracks of the particles in the interactions were classified on the basis of their relative ionization, g = I / I 0 where I 0 and I are respectively the ionizations of singly charged particle and charged particle to be identified. The tracks having relative ionization g 10 are known as black tracks and their number in an interaction is denoted by N b . This ionization cut corresponds to relative velocity β < 0.3 . The grey tracks have relative ionization in the range 0.3 g 0.7 and N g denotes their number in an event. However, relativistic charged particles, called shower tracks have relative ionization, g < 1.4 and this ionization cut corresponds to the relative velocity, β > 0.7 .

Two data samples of 391 and 212 interactions of 16O with AgBr at 60A and 200A GeV respectively having N s 10 were used for analysis using standard emulsion criteria [23] . The pseudorapidity distribution for 60A and 200A GeV, 16O-AgBr interaction are plotted in Figure 1. In the figure, the red regions show the pseudorapidity ranges (1.6 - 3.6 ( η ± 1 ) for 60A and 2.1 - 4.1 ( η ± 1 ) for 200A GeV) which have been used in present analysis.

4. Results and Discussions

The variation of G q as a function of 1 / δ η for 16O-AgBr at 60A and 200A GeV are plotted in Figure 2 and Figure 3. The errors displayed in these figures represent the standard deviations from the mean values of Gq moments.

From the figures, it may be noted that the moments with negative q values

Figure 1. h distribution for 16O-AgBr interactions at 60A and 200A GeV. The red regions show the pseudorapidity ranges of analysis.

Figure 2. G q as a function of 1 / δ η for 6 q 0 .

Figure 3. G q as a function of 1 / δ η for 1 q 6 .

saturate as δ η decreases whereas for positive q values, this shows linearity over a wide range of δ η . This saturation could be due to decrease in number of particles as bin size is reduced. The variation of l n G q as a function of l n δ η also plotted in Figure 4. The data have been fitted using method of least squares for l n δ η 1. The linear rise of the multifractal moments with decreasing bin size of pseudorapidity shows a power-law behaviour in experimental data which is the indication of the self-similarity in the production mechanism of investigated reactions.

The mass exponents ( τ q ) have been calculated for the linear region of plots l n G q verses l n δ η and plotted as a function of q in Figure 5. It is observed from the figure that the values of τ q increases linearly with increasing order of moments and independent of the collision energy.

The generalized dimensions, Dq, have been calculated using Equation (8) and plotted as a function of q in Figure 6 for both data sets. At both the energies Dq decreases with increasing q, which shows the multifractal behaviour in multipartcle production. It may also be mentioned that for positive q values Dq increases with increasing beam energy for same projectile, whereas for negative q values it seems to be independent of the projectile beam energy.

The multifractal spectrum, f ( α q ) , can also give an idea about the presence of multifractal behaviour in multipartcle production. We have calculated the multifractal spectrum f ( α q ) and plotted in Figure 7. It is represented by a continuous concave downward curve with maximum at q = 0, f ( α ( 0 ) ) = D ( 0 ) = 1 and the dotted line represents a common tangent at an angle of 45˚ at α 1 = f ( α 1 ) .

Figure 4. ln G q as a function of l n δ η for 6 q + 6 .

Figure 5. Variations of tq with q for 16O-AgBr interactions at 60A and 200A GeV.

Figure 6. Dq versus q plots for 60A and 200A GeV, 16O-AgBr collisions.

These observations are in accordance with those reported earlier [24] [25] [26] [27] . However, all the spectra are wide enough to indicate the occurrence of multifractality in the multiparticle production process in 60A and 200A GeV 16O-AgBr collisions. It is also observed from Figure 7 that the width of f (α)

Figure 7. f ( α q ) versus α q plots for 60A and 200A GeV, 16O-AgBr collisions.

increases with increasing beam energy for the same projectile.

5. Conclusion

On the basis of the results discussed in the present paper, it is observed that, the moments having positive q values saturate with decrease in the bin size. Further, the variation of mass exponent, τ q , with the order of the moment, q, shows a power law behavior. As far as multifractality is concerned, a decreasing trend in the variation of multifractal dimension, D q , with q is observed. This indicates the presence of multifractality in the particle production process for the considered nuclear reactions. The observed behavior of multifractal spectral function f ( α q ) in pseudorapidity space manifests self similarity in the mechanism of multiparticle production


We are grateful to Professor Anju Bhasin, Jammu University, Jammu, INDIA for providing the experimental data used in the present analysis.

Cite this paper

Ahmad, N., Ahmad, T., Singh, O. and Ahmad, S. (2018) A Study of Multifractal Analysis in 16O-AgBr Collisions at 60A and 200A GeV. Journal of Modern Physics, 9, 1029-1036.


  1. 1. Bialas, A. and Peschanski, R.B. (1986) Nuclear Physics B, 273, 703.

  2. 2. Bialas, A. and Peschanski, R.B. (1988) Nuclear Physics B, 308, 857-867.

  3. 3. Azhinenko, I.V., et al. and EHS/NA22 Collaboration (1989) Physics Letters B, 222, 306-310.

  4. 4. Hwa, R.C. (1990) Physical Review D, 41, 1456.

  5. 5. Chiu, C.B. and Hwa, R.C. (1991) Physical Review D, 43, 100.

  6. 6. Holynski, R., et al. (1989) Physical Review Letters, 62, 733.

  7. 7. Adamovich, M.I., et al. (1990) Physical Review Letters, 65, 412.

  8. 8. Akesson, T., et al. and HELIOS Emulsion Collaboration (1990) Physics Letters B, 252, 303-310.

  9. 9. Florkowski, W. and Hwa, R.C. (1991) Physical Review D, 43, 1548.

  10. 10. Bershadsky, A. (1998) The European Physical Journal A, 2, 223-225.

  11. 11. Hanel, R. and Thurner, S. (2013) Entropy, 15, 5324-5337.

  12. 12. Angulo, J.M. and Esquivel, F.J. (2015) Entropy, 17, 5382-5401.

  13. 13. Deppman, A. (2016) Physical Review D, 93, Article ID: 054001.

  14. 14. Derado, I., Hwa, R.C., Jancso, G. and Schmitz, N. (1992) Physics Letters B, 283, 151-154.

  15. 15. Albajar, C., et al. (1992) Zeitschrift für Physik C, 56, 37-46.

  16. 16. Feder, J. (1988) Fractals. Plenum Press, New York and London.

  17. 17. Paladin, G. and Vulpiani, A. (1991) Physics Reports, 156, 147-225.

  18. 18. Adamovich, M.I., et al. (1996) Journal of Physics G, 22, 1469-1481.

  19. 19. Adamovich, M.I., et al. (1989) Physics Letters B, 227, 285-290.

  20. 20. Adamovich, M.I., et al. (1988) Physics Letters B, 201, 397-402.

  21. 21. Ahmad, S., Ahmad, A., Chandra, A., Zafar, M. and Irfan, M. (2013) Physica Scripta, 4, Article ID: 045201.

  22. 22. Ahmad, S., Khan, M.M., Khan, S., Khatun, A. and Irfan, M. (2014) International Journal of Modern Physics E, 23, Article ID: 1450065.

  23. 23. Hasan, R., Zafar, M., Ahmad, M.S. and Tufail, A. (1999) International Journal of Modern Physics A, 14, 3451-3465.

  24. 24. Zhang, D.H., et al. (2008) Radiation Measurments, 43, S229-S233.

  25. 25. Jain, P.L., Singh, G. and Mukhopadhyay, A. (1992) Physical Review C, 46, 721-726.

  26. 26. Ahmad, N., et al. (2014) Journal of Modern Physics, 5, 1288-1293.

  27. 27. Ahmad, N., et al. (2017) Multifractal Analysis of Multiplicity Fluctuations in 16O-AgBr Interactions at Relativistic Energies. Proceedings of the DAE Symposium on Nuclear Physics, Patiala, 20-24 December 2017, 62.