arxiv: v1 [cond-mat.stat-mech] 10 Nov 2009

Transcription

1 APS/3-QED Generalization of symmetric -stable Lévy distributions for > Sabir Umarov Department of Mathematics, Tufts University, Medford, MA 055, USA arxiv:09.009v [cond-mat.stat-mech] 0 ov 009 Constantino Tsallis Centro Brasileiro de Pesuisas Fisicas ational Institute of Science Technology for Complex Systems Rua Dr. Xavier Sigaud 50, Rio de Janeiro, RJ, Brazil Santa Fe Institute, 399 Hyde Park Road, Santa Fe, M 8750, USA Murray Gell-Mann Santa Fe Institute, 399 Hyde Park Road, Santa Fe, M 8750, USA Stanly Steinberg Department of Mathematics Statistics, University of ew Mexico, Albuuerue, M 873, USA The -stable distributions introduced by Lévy play an important role in probabilistic theoretical studies their various applications, e.g., in statistical physics, life sciences, economics. In the present paper we study seuences of long-range dependent rom variables whose distributions have asymptotic power law decay, which are called (, )-stable distributions. These seuences are generalizations of i.i.d. -stable distributions, have not been previously studied. Longrange dependent (, )-stable distributions might arise in the description of anomalous processes in nonextensive statistical mechanics, cell biology, finance. The parameter controls dependence. If = then they are classical i.i.d. with -stable Lévy distributions. In the present paper we establish basic properties of (, )-stable distributions, generalize the result of Umarov, Tsallis Steinberg (008), where the particular case =, [, 3), was considered, to the whole range of stability nonextensivity parameters (0,] [, 3), respectively. We also discuss possible further extensions of the results that we obtain, formulate some conjectures. PACS numbers: I. ITRODUCTIO The central limit theorem (CLT) -stable distributions have rich applications in various fields including the Boltzmann-Gibbs (BG) statistical mechanics. The nonextensive statistical mechanics [,, 3, 4, 5, 6] characterized by the nonextensivity index (which recovers the BG theory in the case = ) studies, in particular, strongly correlated rom states, mathematical models of which can be represented by specific long-range dependent rom variables. The -central limit theorem consistent with nonextensive statistical mechanics was established in paper [7]. The main objective of [7] was to study the scaling limits (attractors) of sums of -independent rom variables with a finite ( )- variance. The mapping F : G [] G z() [], () Electronic address: Sabir.Umarov@tufts.edu Electronic address: tsallis@cbpf.br Electronic address: mgm@santafe.edu Electronic address: stanly@math.unm.edu where F is the -Fourier transform (see Section ) z(s) = ( + s)/(3 s), G [] is the set of -Gaussians up to a constant factor (see, e.g. [, 3]), was essentially used in the description of attractors. The number in the notation will soon become transparent in the context of the current paper. In the present work we study a -analog of the - stable Lévy distributions. In this sense, the present paper is a conceptual continuation of [7]. The classic theory of -stable distributions was originated by Paul Lévy developed by Lévy, Gnedenko, Feller others; for details history see, for instance, [8-4] references therein. Distributions with asymptotic powerlaw decay (-stables, particularly -Gaussians, in the first place) found a huge number of applications in various practical studies (see, e.g., [4-3], just to mention a few), confirming the freuent nature of these distributions. As it will become clear later on, (, )-stable distributions unify both of them. Indeed, (, )-stable distributions correspond to the -stable ones, the (, )-stable distributions correspond to the -Gaussian ones. All (, )-stable distributions, except Gaussians ((, )-distributions), exhibit asymptotic power-laws. In practice the researcher is often interested in identification of a correct attractor of correlated states, which plays a major role in the adeuate modeling of physical phe-

2 nomenon itself. This motivates the study of seuences of (, )-stable distributions their attractors, as focused in the present paper. For simplicity we will consider only symmetric (, )- stable distributions in the one-dimensional case (see [8] for the multivariate -CLT). We denote the class of rom variables with (, )-stable distributions by L []. A rom variable X L [] has a symmetric density f(x) with asymptotics f C x + +( ), x, where <, 0 < <, C is a positive constant. Hereafter g(x) h(x), x a, means that g(x) lim x a h(x) =. Linear combinations properly scaling limits of seuences of -independent rom variables with (, )-stable distributions are again rom variables with (, )-stable distributions, justifying that L [] form a class of stable distributions. To this end, we note that L [] shares the same asymptotic behavior with the set L sym (γ) of symmetric Lévy distributions centered at 0, where γ = γ(, ) = ( ) + ( ). However, there is an essential difference between (, )- stability of -independent rom variables the classic stability of -stable distributions. amely, - independence exhibits a special long-range correlation between rom variables (see the exact definition in Section 3). In practice this notion reflects physical states (arising e.g. in nonextensive statistical mechanics), which are strongly correlated. The term global correlation instead of strong correlation is also used widely in physics literature; see, e.g. [3, 4]. Examples of such systems include earthuakes [4], cold atoms in optical dissipative lattices [5], dusty plasma [6]. A decomposition of nonextensive processes with strong correlation into independent states can not adeuately reflect their evolution. Likewise, (, )-stable distributions can not be captured by the existing theory of -stable distributions, which is heavily based on the concept of independence (or weak dependence). This distinction ends up with essential implication: attractors of (, )-stable distributions are different from the attractors of -stable distributions, unless =. If = then correlation disappears, that is - independence becomes usual probabilistic independence, γ(, ) =, implying L [] = L sym []. Following the method established in [7], we will apply F -transform for the study of seuences of -independent (, )-stable distributions. Parameter controls correlation. We will classify (, )-stable distributions depending on parameters < (or euivalently Q < 3, Q = ) 0 <. We establish the mapping F : G L[] G [], () where G [] is the set of functions {be β ξ, b > 0, β > 0}, L = 3 + Q +, Q =, i.e., L = + + ( ). The particular case = Q = recovers L = 3+ +, already known in the literature []. Denote Q = {(Q, ) : Q < 3, = }, Q = {(Q, ) : Q < 3, 0 < < } Q = Q Q. ote that the case (Q, ) Q for -independent rom variables with a finite Q-variance was studied in [7]. For (Q, ) Q the Q-variance is infinite. We will focus our analysis namely on the latter case. ote that the case =, in the framework of this classification like the classic -stable distributions, becomes peculiar. In the scope of second classification we study the attractors of scaled sums, exp the results of paper [7] to the region Q generalizing the mapping () to F ζ() : G [] G z()[], <, 0 <, (3) where ζ (s) = ( ) z (s) = + +. ote that, if =, then ζ () = z () = ( + )/(3 ), thus recovering the mapping (), conseuently, the result of [7]. These two classifications of (, )-stable distributions based on mappings () (3) respectively, can be unified to the scheme L [] F G [] F ζ() G ζ() [] (4) F G L[], which gives the full picture of interrelations of (, )- stable distributions with parameters [, ) (0, ) (see details in Section VIII). II. PRELIMIARIES AD AUXILIARY RESULTS A. Basic operations of -algebra In this section we briefly recall the basic operations of -algebra. Indeed, the analysis we will conduct is entirely based on the -structure of nonextensive statistical mechanics (for more details see [3, 4, 5] references therein). To this end, we recall the well known fact that the classical Boltzmann-Gibbs entropy S BG = i p i lnp i satisfies the additivity property. amely, if A B are two independent subsystems, then S BG (A + B) = S BG (A) + S BG (B). However, the - generalization of the classic entropy introduced in []

3 3 given by S = P i p i with R S = S BG, does not possess this property if. Instead, it satisfies the pseudo-additivity (or -additivity) [,, 4] S (A + B) = S (A) + S (B) + ( )S (A)S (B). Inherited from the right h side of this euality, the -sum of two given real numbers, x y, is defined as x y = x + y + ( )xy. The -sum is commutative, associative, recovers the usual summing operation if = (i.e. x y = x + y), preserves 0 as the neutral element (i.e. x 0 = x). By inversion, x y +( )y. we can define the -subtraction as x y = The -product for x, y is defined by the binary relation x y = [x + y ] +. Here the symbol [x] + means that [x] + = x if x 0, [x] + = 0 if x < 0. This operation also commutative, associative, recovers the usual product when =, preserves as the unity. The -product is defined if x + y. Again by inversion, it can be defined the -division: x y = (x y + ). B. -generalization of the exponential cyclic functions ow let us recall the main properties of two functions, -exponential -logarithm, which will be essentially used in this paper. Let e x ln x denote respectively the functions e x = [ + ( )x] +, ln x = x, (x > 0). The entropy S then can be conveniently rewritten in the form S = i p i ln p i. For the -exponential the relations e x y = e x e y e x+y = e x e y hold true. These relations can be written euivalently as follows: ln (x y) = ln x+ln y, ln (xy) = (ln x) (ln y). The -exponential -logarithm have the asymptotics e x = + x + x + o(x ), x 0, (5) ln ( + x) = x x + o(x ), x 0, (6) respectively. The -product -exponential can be extended to complex numbers z = x+iy (see [7, 7, 8]). In addition, for the function e z can be analytically extended to the complex plain except the point z 0 = /( ) defined as the principal value along the cut (, z 0 ). If <, then, for real y, e iy e iy K (+y ) ( ), y, with K = ( ) /( ). Similarly, if >, then 0 < e iy e iy 0 if y. For complex z it is not hard to verify the power series representation e z = + z + z n=0 A n () (n + )! zn, z <, (7) where A n () = n k=0 a k(), a k () = k( ). Let I = ( /, / ). Then it follows from (7) that for arbitrary real number x I the euation e ix = { x n=0 + i{x x n=0 ( ) n A n () x n } (n + )! ( ) n A n+ () x n+ } (n + 3)! holds. Define for x I the functions -cos -sin by formulas cos (x) = x sin (x) = x x n=0 n=0 ( ) n A n () x n, (8) (n + )! ( ) n A n+ () x n+. (9) (n + 3)! In fact, cos (x) sin (x) can be defined for all real x by using appropriate power series expansions. Properties of -sin, -cos, corresponding -hyperbolic functions, were studied in [9]. Here we note that the -analogs of Euler s formulas read cos (x) = eix e ix = cos (x) + i sin (x), + e ix, sin (x) = eix e ix It follows from the definitions of cos (x) sin (x), from the euality (e x ) = e x (+)/ (see Lemma. in [7]), that cos (x) = e ( )x sin (x). (0) Denote Ψ (x) = cos x. Then euation (0) implies Ψ (x) = (e ( )x ) sin (x). () The following two properties of Ψ will be used later on. Proposition II. Let. Then. Ψ (x) 0;. Ψ (x) = x + o(x 3 ), x 0. Proof. Assume. Since e ( )x, then () immediately implies that Ψ (x) 0. Further, sin (x) can be written in the form (see [9]) sin (x) = ρ (x) sin[ϕ (x)], where ρ (x) = (e ( )x ) / i.

4 4 ϕ (x) = arctan( )x. A simple calculation yields Ψ (x). Using the asymptotic relation (5), we get e ( )x = ( )x + o(x 3 ), x 0. () In turn, it follows from (9) that sin (x) = x + o(x 3 ), x 0. (3) ow (), () (3) imply the second part of the statement. Representation (7) shows the behaviour of - exponential near the origin. It is not hard to verify that in the case > for x > ( ) the representation holds. e x = [( )x] { ( ) x + ( ) 4 x n=0 ( ) n A n () (n + )!( ) n ( x )n } C. -Fourier transform for symmetric densities The -Fourier transform for was introduced in [7] used as a basic tool in establishing the -analog of the stard central limit theorem. Formally the - Fourier transform for a given function f(x) is defined by F [f](ξ) = e ixξ f(x)dx. (4) For discrete functions f k, k = 0, ±,..., this definition takes the form F [f](ξ) = e ikξ f(k). (5) k= In the future we use the same notation in both cases. We also call (4) or (5) the -characteristic function of a given rom variable X with an associated density f(x), using the notations F [X] or F [f] euivalently. It should be noted that, if in the formal definition (4) f is compactly supported then integration has to be taken over this support, although, in contrast with the usual analysis, the function e ixξ f(x) under the integral does not vanish outside the support of f. This is an effect of the -product. The -Fourier transform for nonnegative f(x) can be written in the form F [f](ξ) = f(x)e ixξ(f(x)) dx. (6) We note that, if the -Fourier transform of f(x) defined by (4) exists, then it coincides with (6). The -Fourier transform determined by the formula (6) has an advantage to compare to the formal definition: it does not use the -product, which is, as we noticed above, restrictive in use (for ). Proposition II. Let f(x) be an even function. Then its -Fourier transform can be written in the form F [f](ξ) = f(x)cos (xξ[f(x)] )dx. (7) Proof. otice that, because of the symmetry of f, e ixξ f(x)dx = Taking this into account, we have F [f](ξ) = ( e ixξ ow due to (6) we obtain F [f](ξ) = which coincides with (7). Further, denote e ixξ f(x) + e ixξ f(x) eixξ[f(x)] f(x)dx. f(x) ) dx. + e ixξ[f(x)] dx, H, = {f L : f(x) C x + +( ), x }. For a given f H, the constant C = C f is defined uniuely by f. It is readily seen that φ(, ) = + +( ) > for all (0, ) [, ). Moreover, φ(, )( ) < 3 for all (0, ) [, ), which implies σ (f) =. otice also φ(, ) = + (, ), where = (, ) = ( ) + ( ). On the other h, the density g of any -stable Lévy distributions has the asymptotic behaviour g(x) C/ x +, x. Hence, for a fixed [, ) H, is asymptotically euivalent to the set of densities of - stable Lévy distributions. The following proposition plays a key role in our further analysis. Proposition II.3 Let f(x), x R, be a symmetric probability density function. Further, let either (i) the ( )-variance σ (f) <, (associated with =, < ), or (ii) f(x) H,, where (, ) Q. Then, for the -Fourier transform of f(x), the following asymptotic relation holds true: where µ, = F [f](ξ) = µ, ξ + o( ξ ), ξ 0, (8) σ ν, if = ; (+( ))C f 0 with ν (f) = [f(x)] dx. Ψ (y) y + dy, if (, ) Q. (9)

5 5 Proof. First, assume that =. By Proposition II. F [f](ξ) = = (e ixξ ) f(x)dx f(x)cos (xξ[f(x)] )dx. (0) Making use of the asymptotic expansion (5) we can rewrite the right h side of (0) in the form F [f](ξ) = f(x){ + ixξ[f(x)] /x ξ [f(x)] ( ) }dx + o(ξ 3 ) = (/)ξ σ ν + o(ξ 3 ), ξ 0, () from which the first part of Proposition follows. ow, we assume (, ) Q. We apply Proposition II. to obtain F [f](ξ) = = + 0 f(x)[cos (xξ[f(x)] ) ]dx f(x)ψ ( xξ[f(x)] )dx f(x)ψ ( xξ[f(x)] )dx, where is a sufficiently large finite number. In the first integral we use the asymptotic relation Ψ( x ) = x + o(x 3 ), x 0, which follows from Proposition II., get 0 f(x)ψ ( xξ[f(x)] )dx = ξ x f (x)dx + o(ξ 3 ), ξ 0, () 0 that is a uantity of order o( ξ δ ), ξ 0, for any δ <. In the second integral taking into account the hypothesis of the proposition with respect to f(x), we have = C f f(x)ψ ( xξ[f(x)] )dx (3) x + +( ) We use the substitution x +( ) = in the last integral, obtain where (+)( ) Ψ ( x +( ) ξ C f y C f ξ )dx. (4) f(x)ψ ( xξ[f(x)] )dx = µ, ξ + o( ξ ), ξ 0, (5) µ, = ( + ( ))C f 0 Ψ (y) dy. y+ Hence, the obtained asymptotic relations () (we take δ (, )) (5) complete the proof. For stable distributions µ, must be positive. We have seen (Proposition II.) that if, then Ψ (x) 0 (not being identically zero), which yields µ, > 0. ote also that the condition for f(x) to be symmetric was not reuired in [7], if σ (f) <. III. WEAK COVERGECE OF CORRELATED RADOM VARIABLES Let us start this section by introducing the notion of -independence. We will also introduce two types of convergence, namely, -convergence weak -convergence establish their euivalence. By definition, two rom variables X Y are said to be (,, )-independent if F [X + Y ](ξ) = F [X](ξ) F [Y ](ξ). (6) In terms of densities, euation (6) can be rewritten as follows. Let f X f Y be densities of X Y respectively, let f X+Y be the density of X + Y. Then e ixξ f X+Y (x)dx = F [f X ](ξ) F [f Y ](ξ). (7) If all three parameters, coincide, i.e. = =, then we call simply -independent. For = the condition (6) turns into the well known relation F[f X f Y ] = F[f X ] F[f Y ] between the convolution (noted ) of two densities the multiplication of their (classical) characteristic functions, holds for independent X Y. If, then (,, )-independence describes a specific class of correlations. Let X be a seuence of identically distributed rom variables. Denote Y = X +...+X. By definition, X is said to be (,, )-independent (or (,, )-i.i.d.) if the relations F [Y ](ξ) = F [X ](ξ)... F [X ](ξ) (8) hold for all =, 3,... For = = = the condition (8) turns into the condition for the seuence X to be usual i.i.d. If = = then we call the seuence X simply a - i.i.d. Consider example of an (,, )-i.i.d. seuence of rom variables, where (, 3) = (3 )/(+ ). Assume X is the seuence of identically distributed rom variables with the associated Gaussian density β G (β, x) = C e βx, where C is the normalizing constant (see, e.g. [7]). Further, assume the sums X X, =, 3,...,

6 6 are distributed according to the density G (, x), where = β. Then the seuence X satisfies (8) for all =, 3,..., with =, thus being (,, )- independent identically distributed seuence of rom variables. For the sake of simplicity in this paper we will consider only -i.i.d. rom variables. By definition, a seuence of rom variables X is said to be -convergent to a rom variable X if lim F [X ](ξ) = F [X ](ξ) locally uniformly in ξ. Evidently, this definition is euivalent to the weak convergence (denoted by ) of rom variables if =. For denote by W the set of continuous functions φ satisfying the condition φ(x) C( + x ), x R. A seuence of rom variables X with the density f is called weakly -convergent to a rom variable X with the density f if R f (x)dm R f(x)dm d for arbitrary measure m defined as dm (x) = φ (x)dx, where φ W. We denote the -convergence by the symbol. Proposition III. Let >. Then X X 0 yields X X0. The proof of this statement immediately follows from the obvious fact that W is a subset of the set of bounded continuous functions. Recall that a seuence of probability measures µ is called tight if, for an arbitrary ǫ > 0, there is a compact K ǫ an integer ǫ such that µ (R d \ K ǫ ) < ǫ for all ǫ. Proposition III. Let < <. Assume a seuence of rom variables X, defined on a probability space with a probability measure P, associated densities f, is -convergent to a rom variable X with an associated density f. Then the seuence of associated probability measures µ = P(X ) is tight. Proof. Assume that < < X is a - convergent seuence of rom variables with associated densities f associated probability measures µ. We have R [f ](ξ))dξ = R R( F R = R ( R R R R R ( e ixξf It is not hard to verify that ( )dξ ) R f e ixξf dx)dξ dµ (x). (9) R e ixξt dξ = sin (Rx( )t). (30) R R Rx( )t It follows from (9) (30) that R R R ( F [f ](ξ))dξ ( sin (x( )Rf Rx( )f ) ) dµ (x). (3) Since < < by assumption, > as well. It is known [9, 30, 3] that for any > the properties sin (x) (sin (x))/x, x 0 hold. Moreover, (sin (x))/x, x R. Suppose, lim x x f = L,. Divide the set { 0 } into two subsets A = { j 0 : L j > } B = { k 0 : L k }. If A, since sin, there is a number a > 0 such that R R R ( F [f ](ξ))dξ ( ) x a R x ( )f dµ (x) Cµ ( x a), C > 0 A, for R small enough. ow taking into account the - convergence of X to X, if necessary, taking R smaller, for any ǫ > 0, we obtain µ ( x a) CR R R ( F [f 0 ](ξ))dξ < ǫ, A. If B, then there exist constants b > 0, δ > 0, such that f (x) L + δ + δ, x b, B. x x Hence, we have µ ( x > b) = x >b ( + δ) f (x)dx x >b x dx, B. Since, /( ) >, for any ǫ > 0 we can select a number b ǫ b such that µ ( x > b ǫ ) < ǫ, B. As far as A B = { 0 } the proof of the statement is complete. Further, we introduce the function D (t) = D (t; a) = te iat = t( + i( )at ), (3) defined on [0, ], where < < a is a fixed real number. Obviously, D (t) is continuous on [0, ] differentiable in the interval (0, ). In accordance with the classical Lagrange average theorem for any t, t, 0

7 7 t < t there exists a number t, t < t < t such that D (t ) D (t ) = D (t )(t t ), (33) where D means the derivative of D (t) with respect to t. Consider the following Cauchy problem for the Bernoulli euation y t ia( ) y = y, y(0) = 0, (34) t It is not hard to verify that y(t) = D (t) is a solution to problem (34). Proposition III.3 For D (t) the estimate D (t; a) C( + a ), t (0, ], a R, (35) holds, where constant C does not depend on t. Proof. It follows from (3) (34) that y (t) t y + ia( )y = e iat + ia( )t (e iat ) = + ia( )t C( + a ), t (0, ]. ow we are in a position to formulate the following two theorems on the relationship between -convergence weak -convergence. Theorem III.4 Let < < a seuence of rom vectors X be weakly -convergent to a rom vector X. Then X is -convergent to X. Proof. Assume X, with associated densities f, is weakly -convergent to a X, with an associated density f. The difference F [f ](ξ) F [f](ξ) can be written in the form F [f ](ξ) F [f](ξ) = (D (f (x)) D (f(x))) dx, (36) R d where D (t) = D (t; a) is defined in (3) with a = xξ. It follows from (33) (35) that F [f ](ξ) F [f](ξ) C ( + x ) (f (x) f(x)) dx, R d which yields F [f ](ξ) F [f ](ξ) for all ξ R d. Theorem III.5 Let < < a seuence of rom vectors X with the associated densities f is - convergent to a rom vector X with the associated density f F [f](ξ) is continuous at ξ = 0. Then X weakly -converges to X. Proof. Suppose that f converges to f in the sense of -convergence. It follows from Proposition III. that the corresponding seuence of induced probability measures µ = P(X ) is tight. This yields relatively weak compactness of µ. Theorem III.4 implies that each weakly convergent subseuence {µ j } of µ converges to µ = P(X ). Hence, µ µ, or the same, X X. ow applying Proposition III. we complete the proof. IV. SYMMETRIC (, )-STABLE DISTRIBUTIOS AD THEIR PROPERTIES In this section we introduce the symmetric (, )-stable distributions classify them on the base of mapping (). In this classification takes any value in [, ), however we distinguish the cases = 0 < <. Definition IV. A rom variable X is said to have a (, )-stable distribution if its -Fourier transform is represented in the form e β ξ, with β > 0. We denote the set of rom variables with (, )-stable distributions by L []. Denote G [] = {b e β ξ, b > 0, β > 0}. In other words X L [], if F [X] G [] with b =. ote that if =, then G [] represents the set of -Gaussians L [] - the set of rom variables whose densities are -Gaussians, where = (3 )/( + ). Further, from the asymptotic relation (5) we have e β ξ = β ξ + o( ξ ). This Proposition II.3 imply that the associated density of any (, )-stable distribution belongs to H,. Proposition IV. Let -independent rom variables X j L [], j =,.., m. Then for constants a,..., a m, Proof. Let m a j X j L []. j= F [X j ](ξ) = e βj ξ, j =,..., m. Using the properties e x e y = e x+y F [ax](ξ) = F [X](a ξ), it follows from the definition of the - independence that m F [ a j X j ] = e β ξ, β = j= m β j a ( ) > 0. j= Proposition IV. justifies the stability of distributions in L []. Recall that if = then -independent rom variables are independent in the usual sense. Thus, if =, 0 < <, then L [] L sym [], where L sym [] is the set of -stable Lévy distributions.

8 8 Moreover, the appropriately scaling limit of seuences of -independent rom variables with (, )-stable distributions has again a (, )-stable distribution. To this end consider the sum Z = s (, ) (X X ), =,,... where s (, ) is a scaling parameter specified below. First we prove a general result. Theorem IV.3 Assume (, ) Q. Let X be symmetric -independent rom variables all having the same probability density function f(x) H,. Then Z, with s (, ) = (µ, ) ( ), is -convergent to a (, )-stable distribution, as. Proof. Assume (Q, ) Q. Let f be the density associated with X. First we evaluate F [X ] = F [f(x)]. Using Proposition II.3 we have F [f](ξ) = µ, ξ + o( ξ ), ξ 0. (37) Denote Y j = X j, j =,,... Then Z = Y Y. Further, it is readily seen that for a given rom variable X real a > 0, the euality F [ax](ξ) = F [X](a ξ) holds. It follows from this relation that ξ F [Y j ] = F [f]((µ,) ), j =,,... Moreover, it follows from the -independence of X, X,..., the as- / sociativity of the -product that F [Z ](ξ) ξ = F [f]( (µ, ) )... }{{} factors ξ F [f](( ). (38) (µ, ) Further, making use of the expansion (6) for the - logarithm, euation (38) implies ln F [Z ](ξ) = ln F [f]((µ, ) ξ) = ln ( ξ + o( ξ )) = ξ + o(),, (39) locally uniformly by ξ. Hence, locally uniformly by ξ, lim F [Z ] = e ξ G []. (40) Thus, Z is -convergent to a rom variable with (, )-stable distribution, as. Since the density of X L [] is in H [] it follows immediately the following Corollary from Theorem IV.3. Corollary IV.4 Assume (, ) Q. Let X be a seuence of symmetric -independent (, )-stable rom variables. Then Z, with the same s (, ) in Theorem IV.3, -weakly converges to a (, )-stable distribution. ote that = is not included to Q in Theorem IV.3. The case =, in accordance with the first part of Proposition II.3, coincides with Theorem of [7]. Recall that in this case L [] consists of rom variables whose densities are in G [], where = 3 +. Theorem IV.3 also allows to establish a connection between the classic Lévy distributions -Gaussians. L Indeed, for a X L [], its density function f has asymptotics f C f /x (+)/(+( )), x. It is not hard to verify that there exists a L -Gaussian, which is asymptotically euivalent to f. Let us now find L. Any L -Gaussian behaves asymptotically C / x η = C / x /(L ), C j = const, j =,, i.e. η = /( L ). Hence, we obtain the relation + + ( ) = (4). Solving this euation with respect to L, we have L L = 3 + Q, Q =, (4) + linking three parameters:, the parameter of the - stable Lévy distributions,, the parameter of correlation, L, the parameter of attractors in terms of L - Gaussians. Euation (4) identifies all (Q, )-stable distributions with the same index of attractor G L, proving the following proposition. Proposition IV.5 Let Q < 3 (Q = ), 0 < <, 3 + Q + = L, (43) Then the density of X L [] is asymptotically euivalent to L -Gaussian. In the particular case Q =, we recover the known connection between the classical Lévy distributions ( = Q = ) corresponding L -Gaussians. In fact, putting Q = in euation (4), we obtain L = 3 +, 0 < <. (44) + When increases between 0 (i.e. 0 < < ), L decreases between 3 5/3 (i.e. 5/3 < L < 3). It is useful to find the relationship between η = L which corresponds to the asymptotic behaviour of, the attractor depending on (, Q). Using formula (4), we obtain η = ( + ) + (Q ). (45)

9 9 Proposition IV.6 Let X L Q [], Q < 3, 0 < <. Then the associated density function f X has asymptotics f X (x) x η, x, where η = η(q, ) is defined in (45). If Q = (classic Lévy distributions), then (45) implies the well-known fact η = +. Analogous relationships can be obtained for other values of Q. We call, for convenience, a (Q, )-stable distribution a Q-Cauchy distribution, if =. We obtain the classic Cauchy-Poisson distribution if Q =. For Q-Cauchy distributions (43) (45) imply respectively. L (Q) = 3 + Q η = 4 Q +, (46) V. SCALIG LIMITS OF SUMS OF (, )-STABLE DISTRIBUTIOS In this section we generalize the -central limit theorem established in [7] for -Gaussians, that is in the case of =, to symmetrical (, )-stables with any (0, ]. Let < <, f G [], 0 <. It follows from the definition of the -exponential that f C f x, Cf > 0, as x. Analogously, if g G [], then g C g x, Cg > 0, as x. Comparing orders of asymptotics we can easily verify that for a fixed (0, ] for any (, ) there exists a one-to-one mapping M, : G [] G [], + ( ) =, such that the image of a density f G [] is again density. Analogously, there is a one-to-one mapping K, : G [] G [], with the same, such that it maps f(x) = e β x, an element of G [] with the coefficient b = onto the element g(x) = e β x with the same coefficient b =. We notice that if =, then = both operators coincide with the identity operator. Let F be an operator defined as F = K z( ), F M,, where z( ) = + 3. It is readily seen that in the particular case = it coincides with the -Fourier transform, F = F. We call F a generalized -Fourier transform. Proposition V. Assume 0 < let the numbers, be connected through the relationships = ( ) Then the mapping holds. = + ( ) + ( ). (47) F : G [] G (48) Proof. We use the scheme G [] F G [] M, K z( ), (49) G [] F G z( )[] for the proof. Let a density f G [], i.e. asymptotically f(x) C f x /( ), x with some C f > 0. Its image M, [f](x), a -Gaussian G (β; x), in order to be asymptotically euivalent to f, necessarily Hence, G (β; x) C x = + ( ) C f, x. x = + ( ). Further, it follows from Corollary.0 of [7], that where F : G [] G [], = + + ( ) = 3 ( ). ow taking into account the asymptotic euality (the right vertical line in (49)) we obtain G (β ; x) C x = ( ) ( ) C 3 x, x, = + ( ) ( ). Thus, the mapping (48) holds with in euation (47). Let us now introduce two functions that are important for our further analysis: z (s) = s (s ) (s ) where 0 <, s < +, ζ (s) = + (s ) = + = + (s ) (s ), (50) (s ), 0 <. (5) It can be easily verified that ζ (s) = s if =. The inverse, z (t), t (, ), of the the first function reads z t + (t ) (t ) (t) = = + + (t ) + (t ). (5)

10 0 The function z(s) possess the properties: z ( z ) = (s) s z ( s ) =. If we denote z (s), = z (), = z (), then z (, ) = z ( ) =,. (53) Proposition V. implies that for 0 < < min{, + } the following mappings hold: (i) F : G [] G z()[], (ii) F : G z()[] G [], where F is the inverse to F. It should be noted that as Hilhorst [3] noticed - Fourier transform in general is not one-to-one in the space of densities. In paper [33] the invertibility of F in the set of -Gaussians is established. Since mappings M, K, are one-to-one, relationship (49) yields invertibility of F in G z()[] validity of property (ii). Further, we introduce the seuence,n = z,n () = z(z,n ()), n =,,..., with a given = z 0 (), < +. We can extend the seuence,n for negative integers n =,,... as well, setting, n = z, n () = z (z, n()), n =,,.... It is not hard to verify that,n = + ( ) n( ) = n( ) n( ), (54) for all integer n satisfying < n [ ]. The restriction n [/( )] implies the necessary condition,n >, since -Fourier transform is defined for. ote that,n is a function of (, n/), that,n for all n = 0, ±, ±,..., if =, that lim n ± z,n () = for all. Euation (54) can be rewritten as follows:,n n =, n = 0, ±, ±,... (55) We note that the latter coincides with euation (3) of [34], once we identify with the uantity z therein defined, which was obtained through a uite different approach (related to the renormalization of the index emerging from summing a specific expression over one degree of freedom). We also note an interesting property of,n. If we have a -Gaussian in the variable x / ( ), i.e., a - exponential in the variable x, its successive derivatives integrations with respect to x precisely correspond to,n -exponentials in the same variable x. Further, we introduce the seuence,n = ζ(,n), which can be written in the form,n = + ( ) n( ) for n = 0, ±,..., or, euivalently, + (n )( ) =, (56) n( ),n + n =, n = 0, ±,... (57) It follows from Proposition V. definitions of seuences,n,n that F,n : G,n [] G,n+, < n [ ]. (58) Proposition V. For all n = 0, ±, ±,... the following relations hold. Proof. We notice that,n+,n +,n+ =, (59) = On the other h, by (56),n =,n, (60) ( ) (n )( ). ( ) (n )( ) =,n, which implies (59) immediately. The relation (60) can be checked easily. The property,n =,n shows that the seuences (54) (56) coincide if =. Hence, the mapping (58) takes the form F,n : G,n () G,n+ (), recovering Lemma.6 of [7]. Moreover, in this case the duality (59) holds for the seuence,n as well. If < then the values of,n are distinct from the values of,n. The difference is given by,n,n = ( )( ) + n( ), vanishing for =,, or for =,. In the latter case,n =,n. Further, we define for n = 0, ±,..., k =,,..., n + k [ ] +, the operators F k n(f) = F,n+k... F,n [f] = F,n+k [...F,n+ [F,n [f]]...], Fn k (f) = F,n k... F,n [f] = F,n k [...F,n [F,n [f]]...]. In addition, we assume that F k [f] = f, if k = 0 for any appropriate. Summarizing the above mentioned relationships, we obtain the following assertions. Proposition V.3 The following mappings hold:. F,n : G,n [] G,n+ [], < n [ ];

11 . F k n : G,n [] G,k+n [], k =,,..., n = 0, ±,..., < n + k [ ] + ; 3. lim k F k ng [] = G[], n = 0, ±,..., where G() is the set of densities of classic symmetric -stable Lévy distributions. Theorem V.4 Assume 0 < a seuence,n, < n [/( )], is given as in (54) with 0 = [, min{, + }). Let X be a symmetric,k -independent (for some < k [/( )] (0, ])) rom variables all having the same probability density function f(x) satisfying the conditions of Proposition II.3. Then the seuence Z = X X (µ,k,) (,k ) is,k -convergent to a (,k, )-stable distribution, as. Proof. The case = coincides with Theorem of [7]. For k = 0, the first part of Theorem (convergence) is proved in Section 5 of the present paper. The same method can be applied for k. For the readers convenience we proceed the proof of the first part also in the general case, namely for arbitrary k. Suppose that 0 < <. We evaluate F,k (Z ). Denote Y j = X j /s (,k ), j =,,..., where s (,k ) = ( (µ,k,),k). Then Z = Y Y. Again using the relationship F [ax](ξ) = F [X](a ξ), we obtain F,k (Y ) = F,k [f]( ). Further, it ξ follows (µ,k,) from,k -independence of X, X,... the associativity property of the -product that F,k [Z ](ξ) = ξ,k F,k [f]( ), (6) (µ,k,) the right h side of which exhibits the,k -product of ξ identical factors F,k [f]( ). Hence, making (µ,k,) use of the properties of the -logarithm, from (6) we obtain ξ ln,k F,k [Z ](ξ) = ln,k F,k [f]( ) (µ,k,) = ln,k ( ξ + o( ξ )) = ξ + o(),, (6) locally uniformly by ξ. Conseuently, locally uniformly by ξ, lim F,k (Z ) = e ξ,k G,k (). (63) Thus, Z is,k -convergent., To show the second part of Theorem we use Proposition V.3. In accordance with this lemma there exists a density f(x) G,k [], such that F,k [f] = e ξ,k. Hence, Z is,k -convergent to a (,k, )-stable distribution, as. VI. In paper [7] the formula SCALIG RATE AALYSIS ( 3 k ) β k = k 4 k C. (64) k k was obtained for the -Gaussian parameter β of the attractor. It follows from this formula that the scaling rate in the case = is δ = k = k+, (65) where k is the -index of the attractor. Moreover, if we insert the evolution parameter t, then the translation of a -Gaussian to a density in G [] changes t to t /. Hence, applying these two facts to the general case, 0 <, taking into account that the attractor index in our case is,k, we obtain the formula for the scaling rate δ = (,k (66) ). In accordance with Proposition V.,,k = /,k+. Conseuently, δ =,k+ = (k )( ) (k + )( ). (67) Finally, in terms of Q = the formula (67) takes the form δ = (k )(Q ) (k + )(Q ). (68) In [7] it was noticed that the scaling rate in the nonlinear Fokker-Planck euation can be derived from the model corresponding to the case k =. Taking this fact into account we can conjecture that the scaling rate in the fractional generalization of the nonlinear Fokker-Planck euation is δ = + Q, which can be derived from (68) setting k =. In the case = we get the known result δ = /(3 Q) obtained in [35].

12 VII. O ADDITIVE AD MULTIPLICATIVE DUALITIES In the nonextensive statistical mechanical literature, there are two transformations that appear uite freuently in various contexts. They are sometimes referred to as dualities. The multiplicative duality is defined through µ() = /, (69) the additive duality is defined through ν() =. (70) They satisfy µ = ν =, where represents the identity, i.e., () =,. We also verify that (µν) m (νµ) m = (νµ) m (µν) m = (m = 0,,,...). Consistently, we define (µν) m (νµ) m, (νµ) m (µν) m. Also, for m = 0, ±, ±,...,, (µν) m () = ν(µν) m () = m (m ) m + m m + (m + ) m + m = + m( ) + m( ), (7) = + m( ) + m( ) (µν) m µ() = m + + m m( ) = m + (m + ) m( ). We can easily verify, from euations (54) (7), that the seuences,n (n = 0, ±, ±4,...),n (n = 0, ±, ±,...) coincide with the seuence (µν) m () (m = 0, ±, ±,,...). VIII. CLASSIFICATIO OF (, )-STABLE DISTRIBUTIOS AD SOME COJECTURES The -CLT formulated in [7] states that the appropriately scaling limit of sums of k -independent rom variables with a finite ( k )-variance is a k -Gaussian, which is a k -Fourier preimage of a k-gaussian. Here k k are seuences defined as k = k( ), k = 0, ±,..., k( ) k = k, k = 0, ±,... Schematically -CLT in [7] can be represented as {f : σ k (f) < } F k G k [] F k G k []. (7), We have also noticed that -CLT can be described by the triplet (P att, P cor, P scl ), where P att, P cor P scl represent parameters of the attractor, the correlation the scaling rate, respectively. We found that (see details in [7]) for -CLT this triplet (P att, P cor, P scl ) ( k, k, k+ ). (73) Schematically Theorem 3 of the current paper can be represented as L [] F F G [] G L[], 0 < <, (74) where L [] is the set of (, )-stable distributions, G L[] is the set of L -Gaussians with index L defined as L = L 3 + ( ) () =. + Recall that the case = was peculiar we agree to refer to the scheme (7) in this case. Theorem 4 generalizes the -CLT (which corresponds to = ) to the whole range 0 <. Schematically this theorem can be represented as L,k [] F F,k,k G,k [] G,k [], 0 <, (75) generalizing the scheme (7). The seuences,k,k in this case read,k =,k + k( ), k = 0, ±,..., + k( ) ( ) =, k = 0, ±,... + k( ) ote that the triplet (P att, P cor, P scl ) mentioned above, in this case, takes the form (P att, P cor, P scl ) (,k,,k, (/),k+), (76) recovering the triplet (73) in the case =. In connection with the above discussion about triplets, we note that the existence of a -triplet, namely ( sen, rel, stat ), related respectively to sensitivity to the initial conditions, relaxation, stationary state was conjectured in [36]. Later it was observed in the solar wind at the distant heliosphere [37, 38]. The triplet in (73) obtained theoretically might be useful hint for its understing. Finally, unifying the schemes (74) (75) we obtain the general picture for the description of (, )-stable distributions: L,k [] F F,k,k G,k [] G,k [] (77) F G L,k [],

13 3 The rectangle { < Q < 3; 0 < < }, at the right of the classic Lévy line, is covered by non-intersecting curves C L {(Q, ) : 3 + Q + = L }, 5/3 < L < 3. This family of curves describes all (Q, )-stable distributions based on the mapping (74) with -Fourier transform. The constant L is the index of the L -Gaussian attractor corresponding to the points (Q, ) on the curve C L. For example, the green curve corresponding to L = describes all Q-Cauchy distributions, recovering the classic Cauchy-Poisson distribution if = (the green box in the figure). Every point (Q, ) lying on the brown curve corresponds to L =.5. The second classification of (Q, )-stable distributions presented in the current paper, based on the mapping (75) with -Fourier transform leads to a covering of Q by curves distinct from C L. amely, in this case we have the following family of straight lines L {(Q, ) : 4 Q + = 3 }, < 3, (78) where FIG. : (Q, )-regions. L,k = L (,k ) = 3 + (,k ) + In Fig. the dependence of L on parameters (Q, ) Q in the case k = 0 is represented. If Q = = (the blue box in the figure), then the rom variables are independent in the usual sense have finite variance. The stard CLT applies, the attractors are classic Gaussians. If Q belongs to the interval (, 3) = (the blue straight line on the top), the rom variables are not independent. If the rom variables have a finite Q- variance, then -CLT [7] applies, the attractors belong to the family of -Gaussians. ote that runs in [, 5/3). Thus, in this case, attractors ( -Gaussians) have finite classic variance (i.e., -variance) in addition to finite -variance. If Q = 0 < < (the vertical green line in the figure), we have the classic Lévy distributions, rom variables are independent, have infinite variance. Their scaling limits-attractors belong to the family of -stable Lévy distributions. It follows from (44) that in terms of -Gaussians classic symmetric -stable distributions correspond to 5/3<<3 G [8]. If 0 < <, Q belong to the interval (, 3) we observe the rich variety of possibilities of (, )-stable distributions. In this case rom variables are not independent, have infinite variance infinite Q-variance.. which are obtained from (56) replacing n = = Q. For instance, every (Q, ) on the line F-I (the blue diagonal of the rectangle in the figure) identifies -Gaussians with = 5/3. This line is the frontier of points (Q, ) with finite infinite classic variances. amely, all (Q, ) above the line F-I identify attractors with finite variance, points on this line below identify attractors with infinite classic variance. Two bottom lines in Fig. reflect the sets of corresponding to lines { Q < 3; = )} (the top boundary of the rectangle in the figure) { Q < 3; = 0.6)} (the brown horizontal line in the figure). Some conjectures. Both classifications of (Q, )- stable distributions are restricted to the region Q = { Q < 3, 0 < }. This limitation is caused by the tool used for these representations, namely, Q-Fourier transform is defined for Q. However, at least two facts, the positivity of µ, in Proposition II.3 for > max{0, /} (or, the same, Q > max{, /}) continuous extensions of curves in the family C L, strongly indicate to following conjectures, regarding the region {Q < } on the left to the vertical green line (the classic Lévy line) in Fig.. In this region we see three frontier lines, F-II, F-III F-IV. Conjecture VIII. The line F-II splits the regions where the rom variables have finite infinite Q- variances. More precisely, the rom variables corresponding to (Q, ) on above the line F-II have a finite Q-variance,, conseuently, -CLT [7] applies. Moreover, as seen in the figure, the L -attractors corresponding to the points on the line F-II are the classic Gaussians, because L = for these (Q, ). It follows from this fact, that L -Gaussians corresponding to points

14 4 above F-II have compact support (the blue region in the figure), L -Gaussians corresponding to points on this line below have infinite support. Conjecture VIII. The line F-III splits the points (Q, ) whose L -attractors have finite or infinite classic variances. More precisely, the points (Q, ) above this line identify attractors (in terms of L -Gaussians) with finite classic variance, the points on this line below identify attractors with infinite classic variance. Conjecture VIII.3 The frontier line F-IV with the euation Q+ = 0 joining the points (, 0) (, ) is related to attractors in terms of -Gaussians. It follows from (78) that for (Q, ) lying on the line F- IV, the index =. Thus the horizontal lines corresponding to < can be continued only up to the line 4 Q+ F-IV with (, 3 ) (see the dashed horizontal brown line in the figure). If 0, the Q-interval becomes narrower, but -interval becomes larger tending to (, 3). Results confirming or refuting any of these conjectures would be an essential contribution to deeper understing of the nature of (Q, )-stable distributions, nonextensive statistical mechanics, in particular. Finally we note that Fig. corresponds to the case k = 0 in the description (77). The cases k 0 can be treated in the same way. Acknowledgement. We acknowledge thoughtful remarks by R. Hersh, E.P. Borges S.M.D. Queiros. Financial support by the Fullbright Foundation, SI International IH grant P0 GM (USA Agencies), CP Faperj (Brazilian Agencies) are acknowledged as well. [] Tsallis, C. (988) Possible generalization of Boltzmann- Gibbs statistics. J. Stat. Phys. 5, [] Prato, D., Tsallis, C. (999) onextensive foundation of Levy distributions. Phys. Rev. E 60, [3] Gell-Mann, M., Tsallis, C. (004) onextensive Entropy - Interdisciplinary Applications Oxford University Press, ew York. [4] Boon J.P., Tsallis, C. (005) onextensive Statistical Mechanics: ew Trends, ew Perspectives, Europhysics ews 36, (6). European Physical Society. [5] Tsallis, C. (009) Entropy. Encyclopedia of Complexity Systems Science, Springer, Berlin. [6] Tsallis, C. (009) Introduction to onextensive Statistical Mechanics - Approaching a Complex World, Springer, ew York. [7] Umarov, S., Tsallis, C., Steinberg, S. (008) On a - central limit theorem consistent with nonextensive statistical mechanics. Milan J. Math. 76, [8] Umarov, S., Tsallis, C. (007) On multivariate generalizations of the -central limit theorem consistent with nonextensive statistical mechanics. In Complexity, Metastability onextensivity, eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda C. Tsallis, American Institute of Physics Conference Proceedings 965, 34 4, ew York. [9] Gnedenko, B.V., Kolmogorov A.. (954) Limit Distributions for Sums of Independent Rom Variables. Addison- Wesley, Reading. [0] Feller, W. (966) An Introduction to Probability Theory its Applications II, John Wiley Sons, Inc, ew York London Sydney. [] Samorodnitsky, G., Tau M.S. (994) Stable non- Gaussian Rom Processes. Chapman Hall, ew York. [] Uchaykin, V.V., Zolotarev V.M. (999) Chance Stability. Stable Distributions their Applications, VSP, Utrecht. [3] Meerschaert, M.M., Scheffler, H.-P. (00) Limit Distributions for Sums of Independent Rom Vectors. Heavy Tails in Theory Practice. John Wiley Sons, Inc. [4] Zaslavsky, G (00) Chaos, fractional kinetics, anomalous transport. Physics Reports, 37, [5] Metzler, R., Klafter, J. (000) The rom walk s guide to anomalous diffusion: a fractional dynamics approach. Physics Reports 339,, 77. [6] Melbrot, B. (997) Fractals Scaling in Finance. Springer. [7] McCulloch, J. H. (996) Financial applications of stable distributions. In G. S. Maddala C. R. Rao (Eds.), Hbook of Statistics Volume 4. ew York: orth-holl. [8] Mittnik, S., Rachev, S. (993). Modeling asset returns with alternative stable distributions. Economics Reviews, [9] Montroll, E. W., Shlesinger M.F. (98). On /f noise other distribu- tions with long tails. Proc. at. Acad. Sci. USA 79, [0] Montroll, E. W., Bendler, J.T. (984). On Lévy (or stable) distributions the Williams-Watt model of dielectric relaxation. Journal of Statistical Physics 34, 9 6. [] Pierce, R. D. (996) RCS characterization using the alphastable distribution. In Proceedings of the ational Radar Conference, IEEE Press, [] Beck, C., Schloegel, F. (993) Thermodynamics of Chaotic Systems: An Introduction. Cambridge University Press, Cambridge. [3] Schmitt, F.G., Seuront, L. (00) Multifractal Rom Walk in Copepod Behavior. Physica A 30, [4] Caruso, F., Pluchino, A., Latora, V., Vinciguerra, S., Rapisarda, A. (007). Analysis of self-organized criticality in the Olami-Feder-Christensen model in real earthuakes. Phys. Rev. E 75, [5] Douglas, P., Bergamini, S., Renzoni, F. (006) Tunable Tsallis distributions in dissipative optical lattices. Phys. Rev. Lett. 96, 060. [6] Liu, B., Goree, J. (008) Superdiffusion non-gaussian statistics in a driven-dissipative D dusty plasma. Phys. Rev. Lett. 00, [7] Tsallis, C., Queiros, S.M.D. (007) onextensive statistical mechanics central limit theorems I - Convolution of independent rom variables -product. In Complexity, Metastability onextensivity, eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda C. Tsallis. American Institute of Physics Conference Proceedings 965, 8 0, ew York. [8] Queiros, S.M.D., Tsallis, C. (007) onextensive statistical mechanics central limit theorems II - Convolution of - independent rom variables. In Complexity, Metastability onextensivity, eds. S. Abe, H.J. Herrmann, P. Quarati, A. Rapisarda C. Tsallis. American Institute of Physics

15 5 Conference Proceedings 965, 33, ew York. [9] Borges, E.P. (998) A -generalization of circular hyperbolic functions. J. Phys. A: Math. Gen. 3, [30] ivanen, L., Le Mehaute, A., Wang, Q.A. (003) Generalized algebra within a nonextensive statistics, Rep. Math. Phys. 5, [3] Borges, E.P. (004) A possible deformed algebra calculus inspired in nonextensive thermostatistics.physica A 340, [3] Hilhorst H.J. (009) Central limit theorems for correlated variables: some critical remarks. Brazilian Journal of Physics. 39, A, [33] Umarov, S., Tsallis, C. (008) On a representation of inverse F-transform. Physics Letters A. 37, 9, [34] Mendes, R.S., Tsallis, C. (00) Renormalization group approach to nonextensive statistical mechanics, Phys. Lett. A 85, [35] Tsallis, C., Bukman, D.J. (996) Anomalous diffusion in the presence of external forces: exact time-dependent solutions their thermostatistical basis, Phys. Rev. E 54, R97. [36] Tsallis, C. (004) Dynamical scenario for nonextensive statistical mechanics. ews Expectations in Thermostatistics, eds. G. Kaniadakis M. Lissia, Physica A 340, 0. [37] Burlaga, L.F., Vinas, A.F. (005) Triangle for the entropic index of non-extensive statistical mechanics observed by Voyager in the distant heliosphere. Physica A 356, [38] Burlaga, L.F., ess,.f., Acuna, M.H. (007) Magnetic fields in the heliosheath distant heliosphere: Voyager observations during Astrophys. J. 668,