Abstract
The present status of investigations on fluctuations and correlations seen in high energy multiparticle production processes made using the notion of nonextensivity is reviewed.
Fluctuations; Correlations; Nonextensivity
FLUCTUATIONS AND CORRELATIONS
Fluctuations, correlations and non-extensivity
Grzegorz Wilk
The Andrzej Soltan Institute for Nuclear Studies, Hoza 69, 00681 Warsaw, Poland
ABSTRACT
The present status of investigations on fluctuations and correlations seen in high energy multiparticle production processes made using the notion of nonextensivity is reviewed.
Keywords: Fluctuations; Correlations; Nonextensivity
I. INTRODUCTION
There are many examples that distributions of produced particles deviate from the expected exponential form towards its power-like generalization [1], called the q-exponential and know also as Tsallis distribution [2],
described by parameter q (in cases of interest to us X = pT or X = mTcoshy). Distributions of this type were proposed on phenomenological grounds long time ago as useful parametrization of data interpolating between results coming from soft and hard scattering [3]. On the other hand, when seen from the information-theory point of view or from the thermodynamical approach, such distributions arise from the use of non-extensive Tsallis entropy [4] (indexed by q) instead of the usual Shannon-Gibbs-Boltzmann one:
Notice that Sq is nonextensive in such sense that for two independent systems, A and B (in the usual sense, i.e., that the corresponding probabilities of their occurrence factorize, pij(A+B) = pi(A)pj(B)) one has that
The entropic index q is thus also a measure of nonextensivity in system under consideration and therefore is also called non-extensivity parameter.
What are the circumstances leading to q ¹ 1? The common consensus is that it happens whenever [4]: - there are long range correlations in the system (or system is small - notice that our Universe is small with respect to the gravitational interactions!); - there are memory effects of different kind; - the phase-space in which system operates is limited or has fractal structure and, finally, - there are intrinsic fluctuations in the system under consideration.
II. THE ROLE OF INTRINSIC FLUCTUATIONS
Our works, which I shall review now, were concerned with the last possibility mentioned above, namely that | 1-q | measures fluctuations in the hadronizing system under consideration. As shown in [5] in the case when these fluctuations can be described by gamma distribution,
where á ... ñ are the respective averages taken with respect to f. In general case one refers to the concept of the so called superstatistics discussed at length in [6]. This approach has been recently further generalized in [7].
The situation encountered can be visualized in the following way. Already in [8] it was suggested that the overwhelmingly success of thermodynamical models has its origin in the fact that out of the large number N of secondaries produced in a given event only one (very rarely two) is chosen for making distributions. Inevitably, the remain (N-1) act as a kind of heat bath which action can in most cases be described by single parameter T (called usually temperature because of associations with thermal models). In this case one has simple exponential distributions (i.e., q = 1). However, because such "heat bath" is in reality neither homogeneous nor infinite, it is natural to expect that it should be described by more parameters. The simplest extension is to assume that T fluctuates, replace T ® T0 = á T ñ and introduce another (second) parameter describing its fluctuations - this leads to q ¹ 1 in most natural way. Mathematically it can be realized by introducing to the well-known linear Langevin equation with additive noise, modelling the Brownian motion of some test particle (and leading to a Boltzmann distribution as stationary solution of the corresponding Fokker-Planck (FP) equation), a small multiplicative noise. This leads immediately to Tsallis distributions of relevant observable as exact stationary solution of the corresponding FP equation with (q-1) being directly connected with the strength of this new term (cf. [5, 7]). To close this point it must be mentioned that the role of correlations is so far not excluded but it seems that in the realm of high energy multiparticle production they are, at least at circumstances considered so far, not the dominant component [9, 10].
Experimental data of interest come usually in two categories: as distributions in longitudinal phase space (usually in rapidity y) averaged over transverse momenta, dN/dy, or distributions in transverse momenta, dN/dpT, averaged over y or taken for some window in y. We have therefore also two kinds of "temperatures", TL and TT, and their fluctuations are thus described by, respectively, qL and qT. Out of these two classes only exponential shape of dN/dpT spectra is usually regarded as indication of thermal-like character of such processes with a kind of local thermal equilibrium setting in and characterized by temperature T = TT, which can be deduced from the corresponding slopes of dN/dpT. Any deviation from them is then interpreted as signal for some dynamical (nonequilibrium) effects showing up (like, for example, the flow or decay of resonances, see [11] and references therein). Instead of trying to exclude them one can investigate the possibility that the observed nonexponential spectra result from some new form of equilibrium characteristic of nonextensive thermodynamics (incorporating, among others, also effects of the aforementioned factors). The plausibility of such approach is illustrated in Fig. 1 where examples of both types of distributions are shown [12]. The characteristic feature encountered here is that qL > qT and that TL >> TT (qL ranges in linear fashion from 1.05 for = 20 GeV to 1.33 for
= 1800 GeV whereas variation of qT is limited to between 1.095 for
= 200 GeV and 1.11 for
= 900 GeV; the respective changes of TL are from 1.76 GeV to 55.69 GeV whereas those of TT are much smaller varying from 0.134 GeV to 0.14 GeV). The immediate question arises: what is the meaning of qL, or, equivalently, what is the meaning of fluctuations of "partition temperature" TL? The answer is that qL = 1 + 1/k, where k is parameter defining (in addition to mean multiplicity á n ñ ) the Negative Binomial distributions of multiplicity observed in such collisions [16]. It is because TL ~ M/ á n ñ , where M is energy available for production of secondaries, therefore for M kept constant fluctuations of TL mean automatically fluctuations of á n ñ and those lead immediately to NBD [16]. Interestingly enough NBD arises also in a natural way when particles are produced with energies distributed according to Tsallis distribution (when they are distributed according to Boltzmann distributions one gets Poisson distribution instead, see [9]).
Let us proceed now to recent RHIC data on pT distributions, which consist our potential source of the parameter T, as mentioned above. As shown in [13], using approach based either on nonextensive statistics or on stochastic ideas one can successfully account for the whole range of the observed transverse momenta. This is because in both cases the resultant distributions are intrinsically non-exponential. In [14] the same data were analyzed using Hagedorn model with spectrum of resonances given by r (m):
As seen in Fig. 2, although r (m) introduces already some fluctuations to the system, experimental data can be fitted only with qT > 1 (of the order of qT ~ 1.00015). The new and potentially very interesting fact is that similarly good fit can be obtained with simple q-exponential (i.e., by putting r (m) = 1 in eq. (7)), however, in this case q is noticeably greater and equal to qT = 1.065. This effect is seen for all RHIC data and for all centralities. It strongly suggests that including some well identified sources of fluctuations to simple statistical model (represented here by function r (m) as given by eq. (8)) accounts for some fluctuations and lowers therefore the value of parameter qT. However, as seen in Fig. 2, data are very sensitive to q and, as has been shown in [14], to obtain good estimation of T0 = 1/ b0 one has to resort to nonextensive version of eq. (7) with qT > 1.
We would like to close this Section by mentioning that there is yet another source of fluctuation in the production process, which also leads to nonextensivity but not of Tsallis type, i.e., with gaussian distributions instead of eq. (1). Their source lies in the fact that usually only fraction K (called inelasticity) of the total energy of reaction is used for production of secondaries and this fraction can fluctuate. It turns out that such fluctuations are most naturally described in the frame of the so called Extended Gaussian Ensemble (EGE) by simple gaussian distribution, c (K) ~ exp[ -(K- á K ñ )2/(2s 2) ], (cf. [17] for details and references) and that they can reasonably fit distributions of c (K) for different energies obtained in [16], see Fig. 3. In fact, closer inspection of Fig. 3 shows that still better fit can be obtained when using the so called lorentzian distribution, which is nothing else as q-gaussian with q = 2. It would mean then that even here there are still some other fluctuations seen in data, which are undisclosed by the usual EGE approach.
III. SUMMARY
To summarize: we argue (since quite a time already) that whenever one finds in data a hint for eq. (1) the immediate suspicion should be that they hide some intrinsic fluctuations with strength given by q-1. What is the dynamical origin of these fluctuations is, however, out of the scope of the procedure discussed here as q accounts for all of them. The q > 1 results should then urge experimentalists to devise some special measurements devoted to search for fluctuations.
To allow the reader to make personal judgement in what concerns the possible role of correlations let us list recent attempts to interpret parameter q by some dynamical correlations caused either by incomplete phase space occupance [18] or by some specific changes introduced in generalized form of the Boltzmann equation (either by using corresponding collision rates nonlinear in the one-particle densities or by using nontrivial energy composition rules in the energy conservation constraint part, cf. [19] and references therein). It should also be added at this point that, as advocated in [20], one can also view q > 1 as a general (leading order) finite-size effect.
Finally, one should be aware of the fact that there is still an ongoing discussion on the meaning of the temperature in nonextensive systems. However, the small values of parameter qT deduced from data allow us to argue that, to first approximation, T0 can be regarded as the hadronizing temperature in such system. One must only remember that in general what we study here is not so much the state of equilibrium but rather some kind of stationary state (see [21] and references therein).
Let us close by saying that this subject is still an open issue for further research (like, for example, event-by-event analysis of data [22] or hydrodynamical models [23]).
Acknowledgements
The subject reviewed here has been investigated in collaboration with O.Utyuzh, Z.Wlodarczyk, F.S.Navarra, M.Biyajima, M.Kaneyama, T.Mizoguchi, N.Nakamija, N.Suzuki and T.Osada. Partial support (GW) of the Ministry of Science and Higher Education (grants Nr 621/E-78/SPB/CERN/P-03/DWM 52/2004-2006 and 1 P03B 022 30) is acknowledged.
[1] I. Bediaga, E. M. Curado, and J. M. de Miranda, Physica A 286, 156 (2000); W. M. Alberico, A. Lavagno, and P. Quarati, Eur. Phys. J. C 12, 499 (2000); G. Wilk and Z. Wlodarczyk, Nucl. Phys. B (Proc. Suppl.) A 75, 191 (1999); Physica A 305, 227 (2002); F. S. Navarra, O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Nuovo Cim. C 24, 725 (2001); T. Wibig and I. Kurp, Int. J. High Energy Phys. 0312, 039 (2003).
[2] C. Tsallis, Braz. J. Phys. 29, 1 (1999).
[3] C. Michael and L. Vanryckeghem, J. Phys. G 3, L151 (1977); C. Michael, Prog. Part. Nucl. Phys. 2, 1 (1979); G. Arnison et al. [UA1 Collaboration], Phys. Lett. B 118, 167 (1982); R. Hagedorn, Riv. Nuovo Cim. 6 (10), 1 (1983); P. Steinberg et al., Expression of interest for a comprehensive new detector at RHIC II, arXiv:nucl-ex/0503002.
[4] Cf., for example, C. Tsallis, Physica A 340, 1 (2004), A 344, 718 (2004) and A 365, 7 (2006) and references therein. See also special issue of Europhysics News, Nov-Dec. 2005 (EPSP) (http://www.europhysicsnews.com).
[5] G. Wilk and Z. Wlodarczyk, Phys. Rev. Lett. 84, 2770 (2000) and Chaos, Solitons and Fractals 13/3, 581 (2001).
[6] C. Beck and E. G. D. Cohen, Physica A 322, 267 (2003); F. Sattin, Eur. Phys. J. B 49, 219 (2006).
[7] T. S. Biro and A. Jakovác, Phys. Rev. Lett. 94, 132302 (2005).
[8] L. Van Hove, Z. Phys. C 21, 93 (1985) and C 27, 135 (1985).
[9] G. Wilk and Z. Wlodarczyk, Physica A 376, 279 (2007).
[10] K. E. Bassler, G. H. Gunarante, and J. L. McCauley, MarkoveProcesses, Hurst Exponents, and Nonlinear Diffusion Equations with application to finance, cond-mat/0602316.
[11] B. Müller and J. L. Nagle, Ann. Rev. Nuc. Part. Sci. 56, 93 (2006).
[12] F. S. Navarra, O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Physica A 340, 467 (2004) and A 344, 568 (2004).
[13] M. Biyajima, M. Kaneyama, T. Mizoguchi, and G. Wilk, Eur. Phys. J. C 40, 243 (2005).
[14] M. Biyajima, T. Mizoguchi, N. Nakajima, N. Suzuki, and G. Wilk, Eur. Phys. J. C 48, 597 (2006).
[15] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 91, 172302 (2003).
[16] F.S. Navarra, O.V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Phys. Rev. D 67, 114002 (2003).
[17] T. Osada, O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, Eur. Phys. J. B 50, 7 (2006).
[18] T. Kodama, H.-T. Elze, C.E. Augiar and T. Koide, Europ. Lett. 70, 439 (2005); V. Garcia-Morales and J. Pellicer, Physica A 361, 161 (2006).
[19] T. S. Biro and G. Kaniadakis, Eur. Phys. J. B 50, 3 (2006).
[20] A. S. Parvan, Phys. Lett. A 350, 331 (2006); Phys. Lett. A 360, 26 (2006).
[21] S. Abe, Physica A 368, 430 (2006).
[22] O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, J. Phys. G 26, L39 (2000).
[23] T. Osada and G. Wilk, in preparation.
Received on 23 October, 2006
- [1] I. Bediaga, E. M. Curado, and J. M. de Miranda, Physica A 286, 156 (2000);
- W. M. Alberico, A. Lavagno, and P. Quarati, Eur. Phys. J. C 12, 499 (2000);
- G. Wilk and Z. Wlodarczyk, Nucl. Phys. B (Proc. Suppl.) A 75, 191 (1999);
- Physica A 305, 227 (2002);
- F. S. Navarra, O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Nuovo Cim. C 24, 725 (2001);
- T. Wibig and I. Kurp, Int. J. High Energy Phys. 0312, 039 (2003).
- [2] C. Tsallis, Braz. J. Phys. 29, 1 (1999).
- [3] C. Michael and L. Vanryckeghem, J. Phys. G 3, L151 (1977);
- C. Michael, Prog. Part. Nucl. Phys. 2, 1 (1979);
- G. Arnison et al. [UA1 Collaboration], Phys. Lett. B 118, 167 (1982);
- R. Hagedorn, Riv. Nuovo Cim. 6 (10), 1 (1983);
- P. Steinberg et al., Expression of interest for a comprehensive new detector at RHIC II, arXiv:nucl-ex/0503002.
- [4] Cf., for example, C. Tsallis, Physica A 340, 1 (2004), A 344, 718 (2004) and A 365, 7 (2006) and references therein.
- See also special issue of Europhysics News, Nov-Dec. 2005 (EPSP) (http://www.europhysicsnews.com).
- [5] G. Wilk and Z. Wlodarczyk, Phys. Rev. Lett. 84, 2770 (2000) and Chaos,
- Solitons and Fractals 13/3, 581 (2001).
- [6] C. Beck and E. G. D. Cohen, Physica A 322, 267 (2003);
- F. Sattin, Eur. Phys. J. B 49, 219 (2006).
- [7] T. S. Biro and A. Jakovác, Phys. Rev. Lett. 94, 132302 (2005).
- [8] L. Van Hove, Z. Phys. C 21, 93 (1985) and C 27, 135 (1985).
- [9] G. Wilk and Z. Wlodarczyk, Physica A 376, 279 (2007).
- [10] K. E. Bassler, G. H. Gunarante, and J. L. McCauley, MarkoveProcesses, Hurst Exponents, and Nonlinear Diffusion Equations with application to finance, cond-mat/0602316.
- [11] B. Müller and J. L. Nagle, Ann. Rev. Nuc. Part. Sci. 56, 93 (2006).
- [12] F. S. Navarra, O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Physica A 340, 467 (2004) and A 344, 568 (2004).
- [13] M. Biyajima, M. Kaneyama, T. Mizoguchi, and G. Wilk, Eur. Phys. J. C 40, 243 (2005).
- [14] M. Biyajima, T. Mizoguchi, N. Nakajima, N. Suzuki, and G. Wilk, Eur. Phys. J. C 48, 597 (2006).
- [15] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 91, 172302 (2003).
- [16] F.S. Navarra, O.V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Phys. Rev. D 67, 114002 (2003).
- [17] T. Osada, O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, Eur. Phys. J. B 50, 7 (2006).
- [18] T. Kodama, H.-T. Elze, C.E. Augiar and T. Koide, Europ. Lett. 70, 439 (2005);
- V. Garcia-Morales and J. Pellicer, Physica A 361, 161 (2006).
- [19] T. S. Biro and G. Kaniadakis, Eur. Phys. J. B 50, 3 (2006).
- [20] A. S. Parvan, Phys. Lett. A 350, 331 (2006);
- Phys. Lett. A 360, 26 (2006).
- [21] S. Abe, Physica A 368, 430 (2006).
- [22] O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, J. Phys. G 26, L39 (2000).
Publication Dates
-
Publication in this collection
10 Aug 2007 -
Date of issue
June 2007
History
-
Received
23 Oct 2006