arxiv:math/ v1 [math.ca] 23 Jun 2002

Transcription

1 ON FRACTIONAL KINETIC EQUATIONS arxiv:math/ v1 [math.ca] 23 Jun 2002 R.K. SAXENA Department of Mathematics and Statistics, Jai Narain Vyas University Jodhpur , INDIA A.M. MATHAI Department of Mathematics and Statistics, McGill University, 805 Sherbooke Street West, Montreal, CANADA H3A 2K6 H.J. HAUBOLD Office for Outer Space Affairs, United Nations, P.O. Box 500, A-1400 Vienna, AUSTRIA Abstract. The subject of this paper is to derive the solution of generalized fractional kinetic equations. The results are obtained in a compact form containing the Mittag-Leffler function, which naturally occurs whenever one is dealing with fractional integral equations. The results derived in this paper provide an extension of a result given by Haubold and Mathai in a recent paper (Haubold and Mathai, 2000). 1 Introduction and Preliminaries In terms of Pochammer s symbol we can express the binomial series as (α) n = { 1,n=0 α(α+1)...(α+n 1),n N The Mittag-Leffler function is defined by (1 x) α (α) r x r =. (1) r! E α (x) := n=0 z n Γ(αn+1), (2) 1

2 This function was defined and studied by Mittag-Leffler(Mittag-Leffler, 1902, 1905). We note that this function is a direct generalization of an exponential function, since E 1 (z) := exp(z). It also includes the error functions and other related functions, for we have where The equation E 1/2 (±z 1/2 ) = e z [1+erf(±z 1/2 )] = e z erfc( z 1/2 ), (3) erf(z) := 2 π 1/2 z 0 e u2 du,erfc(z) := 1 erf(z),z C. (4) E α,β (z) := n=0 z n Γ(αn+β) gives a generalization of the Mittag-Leffler function. This generalization was studied by Wiman (1905), Agarwal (1953), Humbert (1953) and Humbert andagarwal(1953)andseveral others. Whenβ = 1, (5)reduces to(2). Both the functions defined by (2) and (5) are entire functions of order 1/α and type 1. A detailed account of these functions is available from the monograph of Erdélyi et al. (1955). The Laplace transform of E α,β (z) follows from the integral e pt t β 1 E α,β (λat α )dt = p β (1 ap α ) 1, (6) 0 where Re(p) > a 1/α,Re(β) > 0, which can be established by means of the Laplace integral e pt t ρ 1 dt = Γ(ρ)/p ρ, (7) 0 where Re(p) > 0, Re(ρ) > 0. The Riemann-Liouville operator of fractional integration is defined as adt ν f(t) = 1 t f(u)(t u) ν 1 du,ν > 0, (8) Γ(ν) a with a Dt 0 f(t) = f(t) (Oldham and Spanier, 1974; Miller and Ross, 1993; Srivastava and Saxena, 2001). By integrating the standard kinetic equation (5) d dt N i(t) = c i N i (t),(c i > 0), (9) 2

3 it is derived that (Haubold and Mathai, 2000) N i (t) N 0 = c i 0 D 1 t N i (t), (10) where 0 Dt 1 is the standard Riemann integral operator. Here we recall that, in the original paper of Haubold and Mathai (2000), the number density of species i,n i = N i (t), is a function of time and N i (t = 0) = N 0 is the number density of species i at time t = 0. By dropping the index i in (10), the solution of its generalized form is obtained (Haubold and Mathai, 2000) as N(t) N 0 = c ν 0D ν t N(t), (11) N(t) = N 0 k=0 ( 1) k (ct) νk Γ(νk +1), (12) By virtue of (2) we can rewrite (12) in terms of the Mittag-Leffler function in a compact form as N(t) = N 0 E ν ( c ν t ν ),ν > 0. (13) In this paper we investigate the solutions of three generalized forms of (11). The results are obtained in a compact form in terms of the generalized Mittag-Leffler function and given in the form of three theorems. 2 Generalized Fractional Kinetic Equations Theorem 1. If ν > 0,µ > 0, then the solution of the integral equation is given by N(t) N 0 t µ 1 = c ν 0D ν t N(t), (14) N(t) = N 0 Γ(µ)t µ 1 E ν,µ ( c ν t ν ), (15) where E ν,µ (.) is the generalized Mittag-Leffler function defined by (5). Proof. We know that (Erdélyi et al., 1954) the Laplace transform of the Riemann-Liouville fractional integral is given by L { 0D σ t f(t);p } = p σ F(p), (16) 3

4 where F(p) = u=0 e pu f(u)du. (17) Projecting the equation (14) to Laplace transform, it gives By virtue of the relation it is found that [ ] L 1 N 0 Γ(µ) p µ {1+(p/c) ν } N(p) = L{N(t);p} = N 0 Γ(µ) p µ {1+(p/c) ν }. (18) L 1{ p ρ} = tρ 1,Re(ρ) > 0, (19) Γ(ρ) { = N 0 Γ(µ)L 1 p µ (1) r [ (p/c) ν ] r } (r)! = N 0 Γ(µ) ( 1) r c νr L 1{ p µ νr} = N 0 Γ(µ) ( 1) r c rν t µ+rν 1 Γ(rν +µ) = N 0 Γ(µ)t µ 1 E ν,µ ( c ν t ν ). (20) The result (15) now readily follows by taking inverse Laplace transform of (18). For µ = 1, we obtain the result given by Haubold and Mathai (2000). Theorem 2. If ν > 0,c > 0,d > 0,µ > 0,Re(p) > d ν/α, c d then for the solution of the equation there holds the formula N(t) N 0 t µ 1 E ν,µ ( d ν t ν ) = c ν 0D ν t N(t), (21) N(t) = N 0 t µ ν 1 c ν d ν [E ν,µ ν( d ν t ν ) E ν,µ ν ( c ν t ν )]. (22) Proof. Projecting(21) to Laplace transform and using(5) and(16), we obtain p µ [1+(p/d) ν ] 1 N(p) = L{N(t);p} = N 0 [1+(p/c) ν ] = N 0p ν µ [ ( 1) r (1) r (d/p) rν ( 1) r (1) r (c/p) rν ] c ν d ν (r)! (r)! 4

5 Hence L 1 {N(p)} = N(t) [ N ] 0 = ( 1) r d rν L 1 (p (µ+rν ν) ) ( 1) r c rν L 1 (p (µ+rν ν) ) c ν d ν [ N 0 ( 1) r d rν t µ+rν ν 1 ( 1) r c rν t µ+rν ν 1 ] = c ν d ν Γ(rν +µ ν) Γ(rν +µ ν) = N 0t µ ν 1 c ν d ν [E ν,µ ν ( d ν t ν ) E ν,µ ν ( c ν t ν )] This completes the proof of (22). When µ = ν +1, theorem 2 reduces to Corollary 2.1. If ν > 0,c > 0,d > 0,c d,re(p) > d ν/α, then for the solution of N(t) N 0 t ν E ν,ν+1 ( d ν t ν ) = c ν 0Dt ν N(t), (23) the following result holds N(t) = N 0 c ν d ν [E ν( d ν t ν ) E ν ( c ν t ν )]. (24) On the other hand if d 0 in (22), we arrive at Corollary 2.2. If c > 0,ν > 0,µ > 0,Re(p) > d ν/α, then for the solution of the following result holds N(t) N 0t µ 1 Γ(µ) N(t) = N 0t µ ν 1 c ν Note: When µ = ν +1, then for the solution of there holds the formula = c ν 0D ν t N(t), (25) [ ] 1 Γ(µ ν) E ν,µ ν( c ν t ν ). (26) t ν N(t) N 0 Γ(ν +1) = cν 0Dt ν N(t), (27) N(t) = N 0 c ν [1 E ν( c ν t ν )], (28) 5

6 where c > 0,ν > 0. The case c = d is given by Theorem 3. If c > 0,ν > 0,µ > 0, then for the solution of the equation the following result holds N(t) N 0 t µ 1 E ν,µ ( c ν t ν ) = c ν 0D ν t, (29) N(t) = N 0 ν tµ 1 [E ν,µ 1 ( c ν t ν )+(1+ν µ)e ν,µ ( c ν t ν )]. (30) Proof. Proceeding in a similar manner, it is observed that N(t) = N 0 L 1 { p µ (1+c ν p ν ) 2} = N 0 = N 0 t µ 1 ( 1) r (r+1) Γ(rν +µ) (ct)rν = N 0 t µ 1 = N 0 t µ 1 ν ( 1) r[1 ν [ = R.H.S. of (30). 3 Conclusions ( 1) r (2) r c rν L 1 {p µ rν } (r)! {(rν +µ 1)+(1+ν µ)}(ct)rν Γ(rν +µ) ( c ν t ν ) r Γ(rν +µ 1) +(1+ν µ ( c ν t ν ) r ] Γ(rν +µ) The fractional kinetic equation (11) has been extended to generalized fractional equations (14), (21), (23), and (29). Their respective solutions are given in terms of the ordinary Mittag-Leffler function and their generalization, which can also be represented as FOX s H-functions. The ordinary and generalized Mittag-Leffler functions interpolate between a purely exponential law and power-like behavior of phenomena governed by ordinary kinetic equations and their fractional counterparts, respectively (Lang, 1999; Hilfer, 2000). A specific example for such behavior is the application of Tsallis statistics (Tsallis, 2002) to phenomena that may arise from fluctuations of temperature or energy dissipation rate (Lavagno and Quarati, 2002). The application of fractional kinetic equations to describe such phenomena has not been fully developed yet. 6

7 References Agarwal, R.P.: 1953, A propos d une note de M. Pierre Humbert, C.R. Acad. Sci. Paris 236, Erdélyi, A., Magnus, W., Oberhettinger, F., and Tricomi, F.G.: 1953, Higher Transcendental Functions, Vol.1, McGraw-Hill, New York-Toronto-London. Erdélyi, A., Magnus, W., Oberhettinger, F., and Tricomi, F.G.: 1954, Tables of Integral Transforms, Vol.1, McGraw-Hill, New York-Toronto-London. Erdélyi, A., Magnus, W., Oberhettinger, F., and Tricomi, F.G.: 1955, Higher Transcendental Functions, Vol.3, McGraw-Hill, New York-Toronto-London. Haubold, H.J. and Mathai, A.M.: 2000, The fractional kinetic equation and thermonuclear functions, Astrophysics and Space Science 327, Hilfer, R. (ed.): 2000, Applications of Fractional Calculus in Physics, World Scientific, Singapore. Humbert, P.: 1953, Quelques resultats relatifs a la fonction de Mittag-Leffler, C.R. Acad. Sci. Paris 236, Humbert, P. and Agarwal, R.P.: 1953, Sur la fonction de Mittag-Leffler et quelques-unes de ses generalisations, Bull. Sci. Math. (Ser.II) 77, Lang, K.R.: 1999, Astrophysical Formulae Vol. I (Radiation, Gas Processes and High Energy Astrophysics) and Vol. II (Space, Time, Matter and Cosmology), Springer-Verlag, Berlin-Heidelberg. Lavagno, A. and Quarati, P.: 2002, Classical and quantum non-extensive statistics effects in nuclear many-body problems, Chaos, Solitons and Fractals 13, Miller, K.S. and Ross, B.: 1993, An Introduciton to the Fractional Calculus and Fractional Differential Equations, John Wiley and Sons, New York. Mittag-Leffler, G.M.: 1903, Sur la nouvelle fonction E α (x), C.R. Acad. Sci. Paris, (Ser.II) 137, Mittag-Leffler, G.M.: 1905, Sur la representatin analytique d une branche uniforme d une fonciton monogene, Acta Math. 29, Oldham, K.B. and Spanier, J.: 1974, The Fractional Calculus: Theory and Applications of Differentation and Integration to Arbitrary Order, Academic Press, New York. 7

8 Srivastava, H.M. and Saxena, R.K.: 2001, Operators of fractional integration and their applications, Applied Mathematics and Computation 118, Tsallis, C.: 2002, Entropic nonextensivity: A possible measure of complexity, Chaos, Solitons and Fractals 13, Wiman, A.: 1905, Ueber den Fundamentalsatz in der Theorie der Funktionen E α (x), Acta Math. 29, Wiman, A.: 1905, Ueber die Nullstellen der Funktionen E α (x), Acta Math. 29,