# Asymptotic solutions of a nonlinear diffusive equation in the framework of -generalized statistical mechanics

###### Abstract

The asymptotic behavior of a nonlinear diffusive equation obtained in the framework of the -generalized statistical mechanics is studied. The analysis based on the classical Lie symmetry shows that the -Gaussian function is not a scale invariant solution of the generalized diffusive equation. Notwithstanding, several numerical simulations, with different initial conditions, show that the solutions asymptotically approach to the -Gaussian function. Simple argument based on a time-dependent transformation performed on the related -generalized Fokker-Planck equation, supports this conclusion.

###### pacs:

05.20.DdKinetic theory and 05.20.-yClassical statistical mechanics and 05.90.+mOther topics in statistical physics## 1 Introduction

In the framework of non-equilibrium thermostatistics, irreversible processes can be often described
by means of Fokker-Planck equations (FPE’s) Frank () whose time evolutions
are characterized by monotonically non-increasing Lyapunov functionals.
In a previous work WS07 (), we
derived a non-linear FPE in the picture of a generalized statistical mechanics based on
the -entropy k-entropy1 (); k-entropy2 (), and discussed its relation with
the associate Lyapunov functional or Bregman type divergence Bregman ().
Based on the monotonic behavior of the Lyapunov functional, one can show that
any initial localized state, which evolves according to
the -generalized FPE, converges to a stationary solution which
minimizes the Lyapunov functional.
In addition, in the linear drift case, the corresponding
stationary solution is nothing but the -Gaussian, which is
a generalization of Gaussian by replacing the standard exponential
with its -generalized version WS06 ().

As well known, the diffusive equations play a fundamental rôle in
the description of several physical phenomena. In particular, the linear diffusive equation arising in the Boltzmann-Gibbs (BG) statistical mechanics has been widely studied jointly with its Gaussian self-similar solution Risken ().
Similarly, the -generalized diffusive equation arises, in a natural way,
in the framework of the Tsallis statistical physics.
This is equivalent to the porous medium equation (PME),
which is a transport equation of gases or fluids in a porous medium.
The main properties of PME are now well known PME ()
and its solutions have been widely investigated in literature Carrillo (); Witelski (), in particular
its self-similar solution is given by the so-called Barenblatt
solution Barenblatt ().

In this work we study a different nonlinear diffusive equation obtained in the picture of the
-generalized statistical mechanics.
As known, a useful method to obtain the solutions of
a partial differential equation is
based on the study of its Lie symmetries and the related group invariant solutions Olver (); Hydon (). We accomplished such analysis for the -diffusive equation considering only the classical Lie symmetries whose generators are functions of the independent and the dependent variables. We will leave the study of the generalized Lie symmetries, whose generators are functions also of the derivatives of the depending variables, to a future work.

In particular, we show that the only
group-invariant solutions for the -diffusive equation are the kink-like solutions (a non-normalizable family of self-similar solutions) and the traveling wave solutions (physically irrelevant because they are divergent).
Next, we explore numerically the evolutions of a localized initial
state, whose spreading is governed by the -diffusive equation.
It is shown that, independently from the initial states, the numerical solutions
approach to a shape which is asymptotically well fitted by the -Gaussian distribution, although this one is not a group invariant solution of the equation under investigation.

The paper is organized as follows. In the next section 2 we briefly review the -generalized thermostatistics and its associated non-linear FPE. In section 3, after presenting an alternative derivation of the -diffusive equation, we classify its classical Lie symmetries and their related group invariant solutions. Section 4 deals with the numerical analysis and the final section 5 is devoted to summary.

## 2 -generalized thermostatistics

We briefly review the generalized thermostatistics based on the -entropy k-entropy1 (); k-entropy2 () given by

(1) |

where is the -logarithm defined us

(2) |

The -entropy is an extension of the BG entropy
by means of a real parameter .

The inverse function of , namely , given by

(3) |

is called -exponential.

The -logarithm and -exponential functions are
the building blocks of the statistical mechanics based on the -entropy. In particular, the -exponential is a monotonic increasing and convex function, fulfilling the relationship . Moreover,
, and , as the ordinary exponential function does.
For , the -exponential
is well approximated by the standard exponential, whereas for ,
it asymptotically approaches to a power-law: .
In the limit, both and reduce to
the standard logarithm and exponential functions, respectively.
Accordingly, the -entropy reduces to the BG entropy.

Maximizing under the constraints of the linear kinetic energy average and the normalization

(4) |

leads to the -Gaussian function

(5) |

Here is the constant fixing the normalization of the distribution and depends on the Lagrange multiplier which controls the dispersion of the distribution. Finally, the parameters and are -dependent constants given by

(6) |

respectively, and are related each to the other through the relation

(7) |

The -entropy and the corresponding statistical mechanics preserve
several properties of the BG theory k-entropy1 ().
In particular, the Legendre structure of the thermo-statistics theory has been investigated in SW06 () through the introduction of several thermodynamic potentials.

For instance, it is found that the generalized partition function,
as a function of , given by

(8) |

satisfies the Legendre relation

(9) |

and that the -entropy and the generalized partition function are related through the relationship

(10) |

where

(11) |

The function in equation (8) is defined by

(12) |

where the quantity , defined as

(13) |

fulfills the following proprieties and . In the limit, it reduces to the unity: .

Both functions and play an important rôle in the development of a theory based on the -entropy.

Starting from the definition of the partition function , we can introduce a -generalization of the free-energy, according to

(14) |

In this way, it was shown that the -free-energy satisfies the Legendre transformation structures SW06 () summarized by the following relationships

(15) |

### 2.1 -generalized Fokker-Planck equation

In a previous work WS07 (), we have studied the nonlinear FPE associated with the -entropy. It can be introduced as a continuity equation for the density field

(16) |

with the current density

(17) |

and

(18) |

As shown in WS07 () the quantity is a -generalized Lyapunov functional,
which is monotonically non-increasing in time.

The -generalized FPE reads

(19) |

where is a constant diffusion coefficient.

We remark that the current density is the sum of a linear drift current , which describe the standard Uhlenbeck-Ornstein process and a nonlinear diffusive current, given by , which reduces to the standard Fick’s current in the limit. Consequently, in the same limit, equation (19) reduces to the standard linear FPE

(20) |

Since , the functional is minimized for the stationary solutions of equation (19) according to

(21) |

It is easy to verify that the stationary solution of the -generalized FPE is equal to the optimal maximal entropy distribution , i.e.

(22) | ||||

(23) |

with the position .

Without the external potential the drift therm disappear and equation (19) reduces to a purely -diffusive equation whose explicit form is given by

(24) |

Equation (24) is the main subject of our investigation.

## 3 -generalized diffusive equation

### 3.1 A simple physical derivation

First, let us describe a physical derivation of the -diffusive equation (24). In this respect, we recall that the linear diffusive equation

(25) |

can be obtained from the following physical relations:

(26) |

where denotes the average velocity
of the fluid and denotes the pressure of the fluid, which is related to the density of the fluid through a suitable equation of state , with .

The continuity equation assures the conservation of the zero-th moment

(27) |

which corresponds to the total mass of the fluid.

The PME PME () can be derived from the same physical relations by posing, in the equation of state, :

(28) |

Here, is a positive real parameter characterizing the degree
of porosity of the medium and
characterizing also the power-law dependence of polytropic fluid.
In other words, is a polytropic index.

The corresponding PME is given by

(29) |

with .

Interestingly, the same settings can be also applied to the -diffusive equation (24), merely replacing, in the equation of state, the logarithm with its -generalization: .

In this way, the -diffusive equation (24) can be derived according to

(30) |

which describe a special type of barotropic fluid barotropic (). “A special type” means as follows. From the definition of -logarithm (2), the last relation can be written as

(31) |

i.e., a pressure of this type of barotropic fluid depends on the difference between one term with a polytropic index and the other term with .

### 3.2 Classical Lie symmetry analysis

In the following we classify the Lie symmetries of the -diffusive equation.

For the sake of comparison, we first consider the PME (29). In this case, the symmetry generators are given by (please refer to Bluman () for more details)

(32) | |||

(33) | |||

(34) | |||

(35) |

This means that if is a solution of PME, then

(36) | |||

(37) | |||

(38) | |||

(39) |

are also solutions.

Generators (32) and (33) describe translations invariance in - and -coordinates, respectively and their group invariant solutions (GIS) are merely the constant functions.

Generators (34) and (35) describe scaling invariance of independent and dependent variables, respectively. In particular, the GIS for the symmetry has a kink-like shape and in the limit, it reduces to the Erf function.

Further, GIS can be obtained by means of linear combination of
the vectors . For instance, from the generator we obtain the traveling wave solutions which are physically irrelevant because they diverge. Differently, the most general self-similar generator is provided by vector . Although the corresponding GIS can not be given in terms of elementary functions, explicit solutions can be obtained for special values of the constant . In particular, for we derive
the Barenblatt scale invariant solution Barenblatt ()

(40) |

where is the normalization constant. Note that this solution is nothing but a -Gaussian with gq ().

Next, we study the symmetries for the more complicated -diffusive equation. It is worth noting that the PME is characterized by the single-diffusive term , which leads to the GIS (39),
while the -diffusive equation contains two different powers
and .
This double nonlinearity in the diffusive term destroys the scaling invariance for the dependent variable .

Thus, the only Lie symmetry generators are

(41) | |||

(42) | |||

(43) |

so that, if is a solution of the -generalized diffusive equation then

(44) | ||||

(45) | ||||

(46) |

are also solutions.

Again, the GIS corresponding to the traveling wave generator are physically irrelevant because merely constants or divergent.
Differently, the only self-similar generator is and the corresponding GIS has a kink-like shape as shown in figure 1.
These solutions, although not normalizable, can describe the evolution of an order parameter belonging to the system. However, as compared to the case of PME, there is no scaling generator
like equation (35) and consequently, the -Gaussian is not an exact solution for
the -diffusive equation.

## 4 Asymptotical analysis

### 4.1 Numerical simulations

To further study the long time behavior of the -diffusion,
we performed numerical simulations against various initial probability
distribution functions with different shapes.

We used a numerical method developed in GT06 ()
in order to study the approach to the equilibrium of a given Cauchy problem for the one-dimensional generalized diffusive equation

(47) |

where is an initial localized state and
.
For computational simplicity, the diffusion coefficient is re-scaled in the independent variables according to the symmetry (43), which still holds whatever the function is.

In general, a solution of any diffusive equation spreads
out as time increases. Then, simple conventional
methods of discretizations in the -coordinate, which are often employed in the numerical simulations, leads to much computational costs, especially for a long time numerical simulation.

In order to overcome the difficulty due especially to the expanding support of ,
Gosse and Toscani GT06 () developed an explicit numerical scheme applicable to one-dimensional
differential equations of the type given in (47).

We here explain the basic points of this numerical scheme.
First, we introduce the cumulative distribution function

(48) |

associated with the probability distribution and its inverse relation

(49) |

By using equation (47) and integrating by part we obtain

(50) |

Next, since is non-decreasing in the variable, one can define its pseudo-inverse function for a given value as

(51) |

Then, we have

(52) |

Note that does not depend on time . In other words, at each time there exists only one point of whose cumulative distribution function equals to . Therefore, differentiating both sides of equation (52) w.r.t. we obtain

(53) |

In the case of we derive the time-evolution equation for as

(54) |

In short, we compute a time-evolution of the pseudo-inverse function
for a given initial condition .
Since the computational domain is now restricted we
can go around the above expanding support issue for any large time.

For the -diffusion equation, we set

(55) |

The top plot in figure 2 shows a typical time-evolution of with an initial state concentrated around .

The bottom plot in the same figure 2 shows the time-evolution
for an initial state with double peaks. Note that the double peaks
immediately disappear as time evolves and
become a single peaked shape, which seems to be
asymptotically approaching to the -Gaussian.

We have run several simulations with different initial shapes
to confirm this asymptotic behavior.

A further indication on this asymptotic behavior can be obtained by studying the time evolution
of the function

(56) |

which is the logarithmic ratio between the numerical solution of equation (24) and the trial fitting function

(57) |

This is a -Gaussian where and are the best fitted
parameters for the numerical solution at each time .

In figure 3, we plotted the time evolution of for
an initial with a triangle shape. It is clear from this picture that as , i.e. the function (56) gradually decreases to zero as time evolves, which gives a strong evidence that the numerical solution is asymptotically approaching to the -Gaussian function.

### 4.2 Time dependent transformation

In order to further confirm these asymptotic behaviors let us consider a time-dependent transformation for the -generalized FPE (19), defined by

(58) |

where

(59) |

and

(60) |

Preliminarily, we observe that the transformation (58) maps
a -Gaussian into a -Gaussian.

By transforming equation (19) according to equation (58) we obtain

(61) |

that is the -diffusive equation (24) with the source term .

In the large time limit, the source term becomes negligible and equation (61) can be well approximated by equation (24). On the other hand, we recall that the solution of -generalized FPE (19) approaches asymptotically to the stationary state which is
the -Gaussian (23). Therefore, this solution is transformed, by means of equation (58), into another one that well approximates the solution of the -diffusive equation. In other words, this solution asymptotically approaches to a -Gaussian function.

## 5 Summary

We have studied the asymptotic behaviors of the -diffusive equation,
which is naturally related with the -generalized statistical physics
and the -FPE (19).
The analysis based on the standard Lie symmetry showed that the -Gaussian function is not a scale invariant solution
of the -diffusion equation. Notwithstanding, the -Gaussian
plays a relevant rôle in the study of the time evolution of localized perturbations.
In fact, by performing several numerical simulations, with different
initial shapes, we have found strong evidence that these solutions,
for large time, are well approximated by the -Gaussian. Arguments based on a time-dependent transformation performed on the -generalized FPE also support this result.

Further confirmations of this behavior could be obtained
by studying whether a -Gaussian is invariant under asymptotic
symmetries Gaeta () of the -diffusive equation.

Acknowledgment

This research was partially supported by the Ministry of Education,
Science, Sports and Culture (MEXT), Japan, Grant-in-Aid for
Scientific Research (C), 19540391, 2008.

## References

- (1) T.D. Frank, Nonlinear Fokker-Planck Equations, (Springer-Verlag Berlin Heidelberg 2005).
- (2) T. Wada and A.M. Scarfone, AIP Conf. Proc. 965, 177 (2007).
- (3) G. Kaniadakis, Physica A 296, 405 (2001); Phys. Rev. E 66, 056125 (2002).
- (4) G. Kaniadakis and A.M. Scarfone, Physica A 305, 69 (2002).
- (5) L.M. Bregman, USSR Comp. Math. Phys. 7, 200 (1967).
- (6) T. Wada and H. Suyari, Phys. Lett. A 348, 89 (2006).
- (7) H. Risken, The Fokker Planck equation, (Springer Series in Synergetics, Springer-Verlag, New York, 1989).
- (8) J.L. Vazquez, The Porous Medium Equation: Mathematical Theory (Oxford University Press 2006).
- (9) J.A. Carrillo and G. Toscani, Indiana Univ. Math. J. 49, 113 (2000).
- (10) T.S. Witelski and A.J. Bernoff, Stud. Appl. Math. 100, 153 (1998).
- (11) G.I. Barenblatt, Prikl. Mat. Mekh., 16, 679 (1952).
- (12) P.J. Olver, Application to Lie groups to differential equations, (Graduate Text in Mathematics, Springer-Verlags, New York 1993).
- (13) P.E. Hydon, Symmetry methods for differential equations, (Cambridge Text in Applied Mathematics, Cambridgge University Press 2000).
- (14) A.M. Scarfone and T. Wada, Prog. Theor. Phys. Suppl. 162, 45 (2006).
- (15) A barotropic fluid is one in which the pressure depends only on the density. A common barotropic fluid used in astrophysics is a polytropic fluid, for which the equation of state is written as with a constant and the polytropic index .
- (16) G.W. Bluman and S. Kumei, Symmetries and differential equations, (Applied Mathematical Sciences, Springer-Verlags, New York 1989).
- (17) C. Tsallis and D. Bukman, Phys. Rev. E 54, R2197 (1996).
- (18) L. Gosse, G. Toscani, SIAM Journal on Numerical Analysis, 43, 2590 (2006).
- (19) G. Gaeta, Phys. Rev. A 72, 033419 (2005).