Abstract

From the integration of nonsymmetrical hyperboles, a one-parameter generalization of the logarithmic function is obtained. Inverting this function, one obtains the generalized exponential function. Motivated by the mathematical curiosity, we show that these generalized functions are suitable to generalize some probability density functions (pdfs). A very reliable rank distribution can be conveniently described by the generalized exponential function. Finally, we turn the attention to the generalization of one- and two-tail stretched exponential functions. We obtain, as particular cases, the generalized error function, the Zipf-Mandelbrot pdf, the generalized Gaussian and Laplace pdf. Their cumulative functions and moments were also obtained analytically.

1. Introduction

The convenience of generalizing the logarithmic function has attracted the attention of researchers since long ago [1] and particularly in the last years [27]. In Physics, several one-parameter generalizations of the logarithmic function have been proposed in different contexts such as nonextensive statistical mechanics [813], relativistic statistical mechanics [14, 15], and quantum group theory [16]. Also, more sophisticated such as two-parameter [17] and three-parameter [18] generalizations have been proposed, each one including previous situations as particular cases. Examples of the convenience of these generalizations have been seen in different fields, for instance, psychophysics [19], neuroeconomics [20, 21], econophysics [22, 23], complex networks [24, 25], population dynamics [26, 27], and so forth.

Here, our main objective is to show that the generalized stretched exponential function, written as probability density function (pdf), is suitable to generalize a wide range of one- and two-tail pdfs. This approach, which stresses the emergence of several probability distributions, is motivated by a mathematical curiosity. In Section 2, we show that from the integration of nonsymmetrical hyperboles, one obtains a one-parameter generalization of the logarithmic function, which we call -logarithm. This generalization coincides with the one obtained in the context of nonextensive thermostatistics [9, 10]. Inverting the -logarithm, one obtains the generalized exponential (-exponential) function. Some properties of these generalized functions are presented. In Section 3, the very reliable rank distribution obtained by Naumis and Cocho [28] is conveniently described by the -exponential function, which permits us to detect the effect of finite sample size in the description. In Section 4, we first show that the Zipf-Mandelbrot function, which is a fingerprint of complex systems, can be conveniently written in terms of the -exponential. Raising the -exponential argument to a given power, one obtains a function that generalizes the stretched exponential function. Its generating differential equation is then presented. In Section 5, we consider the pdfs for continuous variables. First we consider the one-tail stretched exponential generalization and obtain analytically its cumulative function and moments. One obtains the generalized error function as a particular case and the Zipf-Mandelbrot pdf as another. Next, we consider the two-tail generalized stretched exponential pdf and obtain analytically its cumulative function and moments. One has the generalized Gaussian and the generalized Laplace pdf as particular cases. The characteristic function is analytically calculated. Our final remarks are drawn in Section 6.

2. The -Generalized Functions

From the integration of nonsymmetrical hyperboles, we obtain a one-parameter generalization of the logarithmic function, which coincides with the one obtained in the context of nonextensive thermostatistics [9, 10]. Inverting this function, one obtains the generalized exponential function. Some properties of these generalized functions are presented.

2.1. -Generalized Logarithm Function

In the one-parameter generalization we address here, the -logarithm function is defined as the value of the area underneath , in the interval [29] which is defined for all real values of , but only for positive . This is exactly the same function obtained from the nonextensive statistical mechanics context [8, 9], but here, only simple geometrical arguments have been used in the derivation.

The usual natural logarithm () is retrieved for and a linear function for . Scaling and deformation of the variable are given by , so that for , one has and for the particular case , . Notice that using instead of , as used in [8, 9], makes simpler to handle the scaling and deformation operations.

2.2. -Generalized Exponential Function

The -exponential function is defined as the -value, in such a way that the area underneath , in the interval , is . This is the inverse of the -logarithm function and it is given by Notice that is real only if .

This is a nonnegative function (, for all and ) and , independently of the value. For , one retrieves the usual exponential function and for , a linear function. Notice that taking the surface underneath to be unitary, one generalizes the Euler’s number: . An interesting property is that , meaning that the -exponential argument scaling corresponds to a power of a different -exponential function. For , one has that .

The derivative of the -exponential function with respect to is , so that it is the solution of the following nonlinear first-order differential equation: , which is a particular case of Bernoulli’s differential equation: , with and . Notice that has the dimension of over the dimension of . This means that it sets up a scale (inversion dimension of ) to the problem only if . One important application of the above differential equation concerns reaction kinetics [30].

3. Beta-Like Distribution

Let us turn our attention to discrete random variables, the rank distribution in particular. We show that the rank distribution obtained by Naumis and Cocho [28] is conveniently described by the -exponential function. In this way, we are able to quantify the finite size effects. This rank distribution is very reliable since it has an underlying microscopic model and has been validated by a wide range of experimental data [28].

To simultaneously fit the beginning, body and tail of experimental rank distributions of complex systems, Naumis and Cocho [28] consider independent subsystems with a large number of internal states. The rank of a system property dependent on the internal states of the subsystems decays as a two-free-parameter beta-like function: , where , with being the maximal value of , is the normalization factor, and the two free parameters are and . As noticed by the authors, if , then , where is the beta function [31] and , is the gamma function [31].

Finite size effects are described by factor , for with . In this way, we see that the rank distribution of is in fact a generalization of the standard technique of multiplying the power-law () by an exponential cutoff: .

If one writes then meaning that the -exponential function can properly take into account finite size effects. The exponential cutoff is retrieved where and grow, but must grow faster than , to have and , with .

4. Special Functions and Processes

In what follows, we first show that the Zipf-Mandelbrot function, which is a fingerprint of complex systems, can be written in terms of the -exponential. Next, raising the -exponential argument to a given power, one obtains a function that generalizes the stretched exponential function. Finally, we obtain the process (differential equation) of which the generalized stretched exponential is the solution.

4.1. Zipf-Mandelbrot Function

The envelope of a typical rank distribution of complex systems can be well described by the Zipf-Mandelbrot function [32], which may be written in terms of -exponential function: with , , and .

We remark that (4.1) has another interesting application in time-dependent luminescence spectroscopy. In this case the relaxation processes are known as Becquerel decay function [33].

4.2. Generalized Stretched Exponential Function

A deformation in the argument of (4.1) leads to the generalized stretched exponential function

The usual stretched exponential function, also known as the Kohlrausch function [3436], is obtained from (4.2) in the limit . Although the (usual) stretched exponential function has been used to describe relaxation processes in time-dependent luminescence spectroscopy [36], a generalization of the form seems to be more convenient to fit experimetal data [37]. In this case we stress the -exponential figures as the argument of the usual exponential function.

The stretched exponential function can also be obtained as the solution of the following differential equation: , which can be written in terms of the logarithmic function (relative variation) as . If one replaces the logarithmic function in the differential equation obtained above by the -logarithm, one obtains or equivalently which is a particular case of Bernoulli’s equation: , with and and which solution is precisely the generalized stretched exponential function of (4.2).

5. Probability Functions

Considering the factor and using the nonnegativeness of (4.2), we write the generalization of the stretched exponential pdf and study its properties. We consider one- and two-tail distributions and obtain some known pdfs (generalized Gaussian) and new ones (generalized error function and generalized Laplace pdf) as particular cases.

5.1. One-Tail PDF

If the considered independent variable is constrained to nonnegative (or eventually nonpositive) values, then one uses the one-tail pdf. We consider the generalization of the stretched exponential pdf and analytically obtain its cumulative function and moments. From this pdf, one obtains the generalized error function as a particular case and the Zipf-Mandelbrot pdf as another.

The normalization factor of (4.2) is

The integral of (5.1) does not diverge only if , , , and one has the generalized stretched exponential pdf:

For our purposes, it is more convenient to write : which is depicted in Figures 1 and 2.

As , , since , one retrieves the stretched exponential function or .

The cumulative function of (5.3) is where

is the hypergeometric function [31] and

is the Pochhammer symbol.

The moments of are where one sees that they are finite only if .

If , the mean value and variance are finite and, respectively, given by Notice that the ratio depends only on and , but not on .

Particular values of lead to a -generalization of the error function and to the Zipf-Mandelbrot pdf.

5.1.1. Generalized Error Function

To generalize the error function, consider and (or ) in (5.4) and one has As , and one retrieves the standard error function: .

An alternative way to obtain (5.9) can be found in [38].

5.1.2. Zipf-Mandelbrot PDF

For in (5.2), one obtains the Zipf-Mandelbrot’s pdf where the mean value and the variance are which is finite for and, from (5.4), one obtains its cumulative function: where the upper-tail distribution is simply given by which is more suitable for fitting the model to real data than the pdf (5.10) itself.

5.2. Two-Tail PDF

If the domain of the considered independent variable is not bounded, it is interesting to consider its absolute value in (5.3) and one has a symmetric pdf about the line . Notice that in this case, the normalization factor must be halved since the domain has been doubled in a symmetrical way, the generalized stretched exponential pdf is then which is a convenient function to use in wavelets [39]. Its cumulative function is where is given by (5.4).

On one hand, due to its symmetry around , the odd moments of this pdf vanish , with On the other hand, the even moments are finite only if :

In the following, we retrieve the generalized Gaussian as a particular case of (5.14). Also, as a new result, we propose to consider another particular case of (5.14) to generalize the Laplacian pdf. The characteristic function of both particular cases is analytically calculated.

5.2.1. Generalized Gaussian

An interesting particular case of the two-tail generalized stretched exponential function (5.14) is when , which leads to the -Gaussian [40] Due to symmetry, all odd moments vanish. Even moments are given by (5.16) and the variance is finite only for

Using that for , , when in (5.17), and one has a Guassian, with variance : , as a particular case. Another particular case is when and one retrieves the Lorentzian (Cauchy pdf) .

The characteristic function of (5.17) has an analytical closed form where is the -modified Bessel function [31]: For the Gaussian one has also a Gaussian but for the Lorentzian, one has the Laplace function .

5.2.2. Generalized Laplace PDF

For and , (5.14) leads to the Laplace pdf . For arbitrary , one has the generalized Laplace pdf and its cumulative function is where is given by (5.10).

The odd moments of (5.22) vanish and the even ones are finite if Its characteristic function has an analytical closed form where and is the hypergeometric function [31], with being the Pochhammer symbol.

6. Conclusion

We have shown that the -generalization of the exponential is suitable to generalize the stretched exponential function. The -generalized stretched exponential function has the generalized error function, the generalized Laplace pdf and the already known generalized Gaussian as special cases. Further, we have used the -exponential to write the very reliable rank distribution obtained by Naumis and Cocho. Since these distributions are the solution of differential equations that describe the complex systems, the -generalization brings many different systems to be described by the same underlying process.

Acknowledgments

A. S. Martinez acknowledges the Brazilian agencies CNPq (303990/2007-4 and 476862/2007-8) for support. R. S. González also acknowledges CNPq (140420/2007-0) for support.