Study of Wavelet Methodology and Chaotic Behavior of Produced Particles in Different Phase Spaces of Relativistic Heavy Ion Collisions

We applied the wavelet methodology for our earlier published research work of the chaotic behavior so called multiplicity fluctuations of secondary charged particles produced during the nucleus-nucleus (A-A) collisions at an energy of the order of ≈ 409 GeV in a new fashion. We illustrated the wavelet coherency in a relation of chaotic behavior for above said data of secondary charged pions in different phase spaces of collisions such as: η-space, φ-space (in one dimension) and ηφ-space (in two dimensions) respectively. We have shown the changes in the wavelet coherence when there are different values of two parameters “q” and “p”. We discussed our new results for the comparison purpose and findings were in the good agreements.


Introduction
The main objective of relativistic heavy ion collisions lies in the investigating of properties of nuclear matter consisting of strongly interacting particles at extreme condition of temperature and density.The most important theory from strongly interacting matter i.e. quantum chromodynamics (QCD) [1] [2] suggests that at a critical temperature of approximately 200 -250 MeV, a phase transition from hadronic matter to a novel state of matter that is quark gluon plasma (QGP) occurred which consists of the deconfined quarks and gluons with in a small volume.The energy density for this new phase transition should be a few GeV/fm 3 (~2 GeV/fm 3 ) compared with 0.17 GeV/fm 3 of normal nuclear matter in its ground state [3] [4].Observations of such type of phase transition would confirm the predictions of QCD and it is also of prime importance in Astrophysics.According to Big Bang theory the initial state of matter consists of dense and hot plasma of unconfined quarks and gluons: after that cosmological expansion process takes place and the energy density decreases, the plasma cools down and transforms into colorless objects [1] [2] [3] [4] [5].This process is called harmonization.Therefore one can expect that the observation of quarks-gluon plasma (QGP) would reproduce in laboratory [5], the hypothetical state of matter that prevailed in early universe.According to phase diagram shown in Figure 1 of strongly interacting matter following points may be noted about the newly state of matter (QGP) [6].
(i) Low temperature and high nuclear energy density is believed to have existed in the early Universe and today it may exist in the heart of neuron star.One can achieve it directly.
(ii) High temperature and low nuclear energy density, which is applicable only by relativistic heavy ion collisions experiments in laboratory which is called Little Bang i.e.Relativistic Heavy Ion Collisions (RHIC) and Large Hadron Collider (LHC) and so many future experiments.
Already now there are plenty of studies that describe different experiments.In references [7] [8] showed that the high energy nucleus-nucleus (A-A) collisions provide a system in which the properties of hot, dense strongly interacting matter can be achieved.Also there have been numerous experimental results, which indicate that collisions of hadron-nucleus (h-A), nucleus-nucleus (A-A) interactions Figure 1.A schematic phase diagram of strongly interacting matter [6].cannot be completely understood in term of superposition of nucleon-nucleon scattering [9] [10] [11].There are many studies where such experimental results have been compared with the data generated with the computer code FRITIOF based on Lund Monte Carlo Model for high energy nucleus-nucleus collisions [12] [13].The measured phenomena, such as strangeness enhancement and J/ψ suppression show that extremely dense strongly interacting matter has been found [14].But still of such research is not enough in order to completely understand the nature of heavy ion collisions at ultra-relativistic energies.
One significant area of research is the study of the dynamics of the chaotic behavior in nuclear collisions at ultra-high energies.Here it is necessary to understand such things: the accuracy of the estimates of the chaotic behavior in relativistic heavy ion collisions in different phase spaces and the connection of such assessments with each other.These questions and is dedicated to present study.To solve this problem we use wavelets ideology [15] [16].This is due to the fact that of the chaoticity and/or errticity pattern in ultra-relativistic heavy ion collisions can have quite difficult structure, contain local features of the various form and time extent.
Finally, In this article, an attempt has been made to study the wavelet methodology over a chaotic system to understand the underlying relation of the chaotic behavior of relativistic shower particles produced in 28 Si-emulsion interactions at 14.6 AGeV.This will help to extract the valuable information's about the mechanism of multi-particle production in different phase spaces; η-space, φ-space (in one dimension) and ηφ-space (in two dimensions).The merits of wavelet analysis over chaotic behavior were foundto be in good agreements.

Brief Detail of Experimental Data
The track detectors (so called nuclear emulsion detector) have been extensively used in high energy physics over many decades.This nuclear emulsion detector also called 4π detector and used for recording of secondary charged particles and is accompanied by the emergence of observable traces (track producing particles) corresponding to the elementary particle trajectories.Therefore in the present experiment, a FUJI film detector has been used to pick the data samples.In this experiment a beam(projectile) of 28 Si nuclei at total energy ≈ 409 GeV hit to the heterogeneous mixture/various nuclei (fixed target) of nuclear emulsion at Alternating Gradient Synchro-phasotron (AGS) of Brookhaven National Laboratory (BNL), NewYork, USA.The other relevant details about the present experiments and target identifications may be seen in our earlier publications [17]- [26].
The projectile spectator fragments, singly and multiple charged, emitted inside fragmentation cone, usually we determine the number of alpha particles (n α ) and the number of fragments with Z > 2, separately [6].
The polar (θ) and azimuthal ( φ) emission angles of all tracks have been measured, and the pseudo-rapidity has been calculated by This study has been carried out for the experimental data along with the theoretical prediction of Ultra-relativistic Quantum Molecular Dynamics model (UrQMD) and Monte-Carlo (RanMC) simulation and the total number of events/data points were 10000.

Analytical Tool for Wavelet Methodology
We made an approach for the generalization of the constant wavelet transformation over time interval "t" which later on changed into an input time series such as; ( ) ( ) where, the usual meanings to symbols such as: s represents to normalization processes, "u" indicates to location parameter and "s" is a scale parameter.
The wavelet analysis transforms the original data to a classified structure by which we get a new set of wavelet coefficients.This is done using so-called multiresolution analysis method.The multiresolution analysis consists in splitting of an investigated number into two components-approximating and detailing, with their subsequent crushing for the purpose of change of level of expansion of a signal to the set level of expansion.The importance of wavelet analysis in the study of time series is to determine the fact that the method of wavelet analysis allows to discover the local features of the studied time series due to the decomposition of the input data [28] [29].And finally, this allows determining the presence of special characteristics of the analyzed data, as well as the point where these characteristics may arise [30].On the other hand, the wavelet analysis allows conducting a detailed analysis for the original data.
In practice, it is quite spread the use of discrete wavelet transformation (DWT) [15].It is connected with that application DWT becomes especially effective when the signal has high-frequency components of short duration and extensive low-frequency components.In the consent of time series study of the wavelet analysis that, it stimulates for the various methods of wavelet transformation and their application in different sectors such as: scaled analysis, cross wavelet transformation, wavelet coherence.
But the leading role of the wavelet transformation methods is to use for gene- ( ) where; Q is a smoothing operator.
In the present work, we used the Morlet wavelet that was a complex wavelet with good time-frequency localization, as a parent one [31] [32].The squared wavelet coherency coefficient was in the rangeof

Analysis for the Chaotic Behavior
To eliminate the effect of non-uniform density distribution, we used the proposed method of Bialas and Gazdzicki [33].In this method the original values of pseudo-rapidity (η) distribution is transformed into a new cumulative variable ( ( ) X ϖ ).And this was demonstrated by following mathematical relation: where the usual meanings of such above parameters were; ( ) ( ) is the single particle (η) distribution of the secondary charged pions and min ϖ and max ϖ are the two extreme points in the density distribution of ( ) ρ ϖ .The effect of this analysis the corresponding region of investigation the interval of pseudo-rapidity∆η, become changed in 0 to 1 order of a variable ( ) Further, for the study of event-to event multiplicity fluctuations in relativistic nuclear collisions, it is necessary to investigate the performance of the event factorial foment, e q F in small bins of variable ( ) X ϖ .The factorial moments (FMs) for the order "q" of an event was mathematically represented such as [34] [35]: where "M" is the partition number of variable ( ) X ϖ , n m is the number of pro- duced secondary charged particles into the m th bin and q = 2 -6 is the order of the moment.
Since the fluctuations of e q F from event-to-event, we get a distribution of e q F denoted by ( ) e q P F after a large number of events.The normalized factorial moments of a particular event can now be defined such as: where It is worth mentioning here that the scaled factorial moments (SFMs) defined by Bialas and Paschanski [34] [35] to study of intermittency in nuclear collisions is only an estimate of mean of the distribution ( ) e q P F .It should be realized that the averaging procedure, apart from its clear advantages, brings also a danger of losing some important information on spatial patterns from event-to-event.In particular, some interesting effects, if present only in a part of sample of events produced in high energy collisions, may be lost.A possible example of this kind is the quark gluon plasma (QGP), which is expected to be characterized by specific intermittency exponents [34] [35].It is therefore, essential to investigate the full shape of the distribution and the way it changes with the bin size.
It has been argued [34]- [39] that the event-to-event fluctuations can probe the dynamics of multiparticle production more deeply than the variables such as the multiplicity distribution and the average factorial moments.Therefore, we study the chaoticity or event-to-event fluctuations in the density of particle produced in 28 Si-emulsion collisions.

( ) ( ) ( )
, 1 where, the usual meaning of various symbols/and or parameters are such as: the "p" is any positive real number, ϖ -a normalized factorial moment of a single event that is associated with the component in 1-D phase space (η-space and φspace) ( [24] [25] [33] [36] and references therein] of ( ) X ϖ -variable ( ( ) X ϖ -a variable that characterizes the number of relativistic heavy ion collisions).

The Chaoticity in Relativistic Heavy Ion Collisions
First we have calculated the Chaoticity, so called chaotic behavior from the present experimental data.For this task, the event factorial moments e q F and the erraticity moments, , p q C were calculated for order of q = 2 -4 and parameter "p" = p = 0.5, 0.9, 1.2, 1.4 and 1.6 by using the Equation (7) and Equation ( 8), the values of "M" was varied from 1 to 40.And the outcomes of all above these calculations were represented in Figure 2, for different phase spaces; η-space, φ-space (in one dimension) and ηφ-space (in two dimensions) respectively.
Hence the dependence of lnC p,q (M) as function of lnM in different phase spaces for the produced secondary charged particle in the collisions of 28 Si (projectile) with fixed target of Emulsion nuclei at energy 14.6A GeV has been depicted in Figure 2. Some part of these results already published by the author M. Ayaz Ahmad et al., in reference [40].From the Figure 2, it can be seen that the expe-Figure 2. Variations of lnC p,q (M) as function of lnM in η-space (1D) (φ-space (1D) and ηφ-space (2D)) in the collisions of 28 Si+Em at energy 14.6A GeV.rimental data of C p,q lie well above those for the generated uncorrelated (Monte Carlo) data in the region of "M ≥ 13.It means that the contribution of the statistical fluctuations to C p,q values in present experimental data is small.Shaoshun and Zhaomin [41] and Jinghua Fu et al., [42]

Results from Wavelet Methodology
Thus, we have a few of time series.These time series should be compared using the ideology of wavelets.For this important task we will use the wavelet coherency.We consider the wavelet coherency for each triple of time series for the produced particles in the relativistic nuclear collisions for various phase spaces; η-space, φ-space (in one dimension) and ηφ-space (in two dimensions) respectively, where the factorial moments order was q = 2 to 4, and the parameter "p" was of the order of, 0.5, 0.9, 1.2, 1.4 and 1.6.
The findings of this analysis were depicted in Figures 3-17.One can check the results of wavelet coherence between selected time series in these Figures 3-17.
Each of the following figures indicated a separate group of time series that match each other.In this case the time scale (x-axis) is equal to the consistent change of values ( ) ln M (these changes represented a sequence number 1, 2, ).The correspondence between a sequence number and value of is shown in Table 1.
In Figures 3-17 (from left to right): the first picture (a) is wavelet coherence between η-space and φ-space, the second picture (b)-wavelet coherence between η-space and ηφ-space, the third picture (c)-wavelet coherence between φ-space and ηφ-space.
Further it has been find that on the vertical axis were the weighted features of the analyzed data series in the frequency space.Corresponding to each of the figures there was a significant scale and is obtainable as separate columns for reflections.The wavelet coherence illustrates the regions in the time-scale space where the central variables co-vary (but do not necessarily have high power).
The maximum reflection is addressed about the full coherence between the data which were analyzed.By definition, the regions inside the black lines plotted in warmer colors indicate a strong interdependence between the investigated time series.The colder color indicates relatively weak co-movement between these variables.The defined lines were a sign of localization for individual irregularities within studied time series according to importance of irregularities.The       In our case, these data suggest interrelations between different spaces (η-space, φ-space (in 1D) and ηφ-space (in 2D) respectively) and the impact of this relationship on the chaotic behavior in relativistic heavy ion collisions.Thus this implies that wavelet plots enable to appropriately identify both frequency bands and time intervals where data series move together.
The Figures 3-7 show the results of wavelet coherence for q = 2.We see that between φ-space and ηφ-space there is full consistency of the chaotic behavior in relativistic heavy ion collisions.Between η-space and φ-space and between η-space and ηφ-space wavelet coherence varies depending on the value  3-5) and then increases (Figure 6 and Figure 7).shows the results of wavelet coherence for q = 3.We see that Figures 3-7 coincide with Figures 8-12.But we can to observe and minor differences (Figure 6 and Figure 11).It is connected with an existing fluctuation of the original data (Figure 6).At the same time it confirms the validity of the results that were obtained.In general, can be seen that with increasing q (from q = 2 to q = 3) significant changes of the chaotic behavior in relativistic heavy ion collisions between different phase spaces does not occur.shows the results of wavelet coherence for q = 4.We see that between η-space and ηφ-space there is full consistency of the chaotic behavior in relativistic heavy ion collisions.Between η-space and φ-space and between φ-space and ηφ-space wavelet coherence varies depending on the value ( ) ln M .Other characteristics of the wavelet coherency have not changed.Thus, with increasing q (from q = 3 to q = 4) we can observe changes in the interaction phase spaces in the study of the chaotic behavior in relativistic heavy ion collisions.Therefore with the change in value of q there is a change influence of the chaotic behavior in relativistic heavy ion collisions, which is reflected in the change the values ( ) M of in the phase space.

Conclusions and Final Remark
The relativistic heavy ion collisions provide an experimental setting for the study of the exotic behavior of the matter.It has been suggested that the strongly interacting matter at high energy densities produced in these collisions may undergo a phase transition to quark gluon plasma (QGP).The produced particles in such nuclear collisions are believed to carry relevant information about the collision mechanism.
We have received different wavelet coherency for the chaotic behavior in relativistic heavy ion collisions between different phase spaces η-space, φ-space (in one dimension) and ηφ-space (in two dimensions) respectively).We have shown changes in the wavelet coherence when changing the values of q and p.This must be considered when analyzing the data.
The ideology of wavelet analysis and/or methodology was found an important tool for the wavelet coherency between the studied data series and it can be implemented very well for the study of relativistic heavy ion collisions.
ralized cross-reference analysis between different time series is wavelet coherence.The wavelet coherency simultaneously assess how the co-movement and causalities between two variables vary across different frequencies involved and change over time in a time-frequency window.Moreover, co-movement of the time series is recognizable among different time scales, which the standard approaches, failed to perform.It allows to calculate local correlation of two time series (x and y) in a region of time-frequency.Hence, one can use the following formalized model; e.g. the wavelet coherence as the squared absolute value of the smoothed cross wavelet spectra ( ) , xy W u s , should normalized by the product of the smoothed individual wavelet power spectra of each series [31] [32]: have studied the chaoticity in NA27 data on p-p collisions at 400 GeV/c.They pointed out that the observed chaoticity could be reproduced by the statistical fluctuations only.But the chaoticity observed in the present experimental work cannot be reproduced by the statistical fluctuations only.Thus the chaotic behavior observed in our data and has the dynamical origin.

Figure 17 .
Figure 17.Wavelet coherence for the experimental data for p = 1.6 and q = 4.
Thus, of the chaotic behavior in relativistic heavy ion collisions from the point of view of different spaces is changing.The smallest values the wavelet coherence is observed for the mean values ( ) ln M .With increasing values p, violation of wavelet coherence to first decreases (Figures and it was found that the values close to zero indicate weak correlation, while values close to one are evidences of strong correlation.Thus, wavelet coherency analysis enables interconnection between the studied time series and analyzes the frequency of such communications.The Monte Carlo methods are used here to properly find statistical level of significance of the wavelet coherence.

Table 1 .
The correspondence between a sequence number and value of