Some Aspects of Extended Kinetic Equation

Motivated by the pathway model of Mathai introduced in 2005 [Linear Algebra and Its Applications, 396, 317–328] we extend the standard kinetic equations. Connection of the extended kinetic equation with fractional calculus operator is established. The solution of the general form of the fractional kinetic equation is obtained through Laplace transform. The results for the standard kinetic equation are obtained as the limiting case.


Introduction
The chemical evolution of a star like sun could be effectively explained by kinetic equations.The kinetic equations explain the rate of change of chemical composition of a star in terms of the thermonuclear reaction rates for destruction and production of the species involved.An arbitrary reaction is characterized by the rate of change dN dt of a time dependent quantity N (t) between the destruction rate d and production rate p.Here the destruction or production at time t depends not only on N (t) but also on the past history N (τ ), τ < t of the variable N .This may be formally represented by following [1,2] where N t denotes the function defined by N t (t * ) = N (t − t * ), t * > 0. It should be noted that d and p are functionals and Equation ( 1) represents a functional-differential equation.If we consider a simplified form of Equation (1) we could consider the decay rate of a radio-active substance which is given by a homogeneous differential equation where N is the number density of the radio-active substance and λ is the decay constant.The solution of this differential equation with initial condition N = N 0 at t = 0 is If we consider a more general form of the differential Equation ( 2) for the decay rate of a radio-active substance as we have the solution of the form where a is a constant.One may get the solution in Equation (3) from Equation ( 5) as α → 1.These types of problems arise in many experimental situations where one needs to switch from one family of functions to another family.In 2005, Mathai [3,4] introduced the pathway model by which one can switch among three different families of functions, say, type-1 beta families, type-2 beta families and gamma families.We get three different functional forms by varying the pathway parameter α.The pathway model in the real scalar case is defined as where a > 0, δ > 0, γ > 0, η > 0. c 1 , c 2 and c 3 are the normalizing constants when we consider the functions as statistical densities.The three different functional forms are respectively generalized type-1 beta, generalized type-2 beta and generalized gamma forms.By writing 1 − α = −(α − 1), the generalized type-2 beta form can be obtained from generalized type-1 beta form.Both generalized type-1 beta form and generalized type-2 beta form reduce to generalized gamma form as α → 1.
Due to this switching property, the pathway model has been widely used in many areas.In this paper, we use the pathway model to extend kinetic equations.The present paper is organized as follows: In the next section we discuss the extended kinetic equation and its solution with a brief description of the extended reaction rate probability integral.Connection of the extended kinetic equation with fractional calculus is examined in Section 3. In Section 4 we try to solve fractional kinetic equations and their various generalizations.Concluding remarks are given in Section 5.

Extended Kinetic Equations
The following discussion is based on [1,2].If we consider a production and destruction of nuclei in the proton-proton chain reaction, we can describe it by the equation where N i is the number density of the species i over time.Here the summation is taken over all reactions, productions or destructions of the species i.The number density N i of the species i can be expressed by the relation A i where ρ is the mass density, X i is the mass abundance, N A is the Avogadro number and A i is the mass of species i in mass units.The mean life time τ j (i) of species i for destruction by species j is given by the relation [2] where λ j (i) is the decay rate of i for interaction with j. σv ij denotes the reaction probability for an interaction involving species i and j defined as where µ is the reduced mass of the particles given by µ is the kinetic energy of the particles in the center of mass system.Consider the cross section σ(E) for low-energy non-resonant reactions given by where Z i and Z j are the atomic numbers of the nuclei i and j, e is the quantum of electric charge, h is the Planck's quantum of action, B the nuclear barrier height, S(E) is the cross-section factor which is a slowly varying function of energy over a limited energy range and which can be characterized depending on the nuclear reaction.The density function of the relative velocities of the nuclei for a non-degenerate and non-relativistic gas is assumed to be Maxwell-Boltzmann given as By substituting Equation (11) and Equation (10) in Equation ( 9) the reaction probability σv ij is obtained as If a deviation from the thermodynamic equilibrium with regard to their velocities is considered then it results in a deviation from the Maxwell-Boltzmann velocity.In this context, we consider a more general density function than Maxwell-Boltzmann density function by using the pathway model defined in Equation (6).The pathway energy density function has the form (πkT ) Replacing the Maxwell-Boltzmann density Function (11) by the pathway energy density Equation ( 13), we get the extended thermonuclear reaction probability integral in the form we get where Following [5], by taking the Mellin transform of Equation ( 16) and simplifying, we get where (s) > 0, (ν + 1 + s 2 ) > 0, By taking the inverse Mellin transform we get, where L is a suitable contour which separates the poles of Γ(s) and Putting s = 2s and using Legendre's duplication formula [6] we get where G 3,1 1,3 (.) is the G-function originally introduced by C.S. Meijer in 1936, see [5,7,8].The G 3,1 1,3 (.) used in Equation ( 21) converges for all (α−1)x 2
The following discussion is adapted from [1].The solution of the differential Equation (4) with initial condition N i (t) = 1 when t = 0 is When c i in Equation ( 32) is a constant, the total number of reactions in the time interval 0 ≤ t ≤ t 0 is obtained as 1 is the probability that the lifetime of species i is ≤ t 0 when t follows a distribution with density function or we have Equation ( 34) is the Tsallis statistics for α > 1, see [13,14] which can also be seen as a particular case of the pathway model Equation ( 6) for α > 1.If c i in Equation ( 32) is a function of time, say c i (t), then it should be replaced by c i (t)dt.When where The density in Equation ( 34) is the lifetime density of the destruction of the species i, with the expected mean value where E(•) is the expected value of (•).The mean value of the lifetime density function given in Equation ( 37) is Now as α → 1 we get the expected mean lifetime of the lifetime density function where E ν (−c ν t ν ) is the Mittag-Leffler function, introduced by M. G. Mittag-Leffler in 1903 [18] as If we consider a generalization of the fractional kinetic equation considered by Mathai, Haubold and Saxena [16,17] in the form where

Conclusions
The linear and non-linear kinetic equations establish a connection between the Boltzmann-Gibbs statistical mechanics and Tsallis non-extensive statistical mechanics.The pathway parameter α plays a key role in switching between these two cases.The theory of extended reaction rates and its closed form solutions can be seen in Haubold and Kumar [10] and Kumar and Haubold [12].Further, the fractional diffusion equation and its solution help us to understand the connection with the classical Laplace transform.In 2013, the author has solved the fractional kinetic equations discussed here by P α -transform [20].Various fractional differential equations and their solution by various transforms are studied by many authors, see [21,22].It should be noted that the Mittag-Leffler function arises in the solution of a fractional diffusion equation whereas the exponential function arises in the solution of its classical counterpart.A possible connection of the extended kinetic equation to fractional calculus can be established through the procedure adopted here.