arxiv:cond-mat/ v1 [cond-mat.stat-mech] 9 Jan 2004

Transcription

1 arxiv:cond-mat/446v [cond-mat.stat-mech] 9 Jan 24 Distributed-Order Fractional Kinetics I.M. Sokolov Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 5, D-2489 Berlin, Germany and A.V. Chechkin Institute for Theoretical Physics, National Science Center Kharkov Institute of Physics and Technology Akademicheskaya st., 68 Kharkov, Ukraine and J. Klafter School of Chemistry, Sackler Faculty of Exact Sciences, Tel Aviv University, Tel Aviv 69978, Israel Fractional diffusion equations are widely used to describe anomalous diffusion processes where the characteristic displacement scales as a power of time. For processes lacking such scaling the corresponding description may be given by distributed-order equations. In the present paper we consider different forms of distributed-order fractional kinetic equations and investigate the effects described by different classes of such equations. In particular, the equations describing accelerating and decelerating subdiffusion, as well as the those describing accelerating and decelerating superdiffusion are presented. PACS 5.4.+j; 2.5.-r. Introduction The diffusion equation proposed by Adolf Fick almost 5 years ago is a partial differential equation of parabolic type, with the first temporal Presented at the 6th Marian Smoluchowski Symposium on Statistical Physics: Fundamentals and Applications, September 6-, 23 ()

2 2 derivative on its l.h.s. and the second spatial derivative on its r.h.s.. The corresponding equation can be put down both for the particles concentration in a diffusing cloud and for a probability of a single particle s position: 2 p(x,t) = K t x2p(x,t) (.) The mean squared particle s displacement from its initial position given by the solution of this equation grows linearly in time: x 2 (t) = 2Kt. This scaling behavior follows immediately from the structure of Fick s equation, being second order in spatial coordinate and first order in time: Changing the spatial scale by a factor of 2 corresponds to changing the time scale by a factor of 4. In complex systems this kind of behavior is often violated, being replaced by an anomalous diffusion relationship, x 2 (t) t β, (.2) with β, or (in absence of the second moment) by related forms where the lower, not necessarily integer, moment of the distribution scales as a function of time. The continuous description of anomalous diffusion is given by generalizations of Fick s scheme based on fractional derivatives []. The fractional generalizations of the diffusion equation may have either a corresponding fractional derivative instead of the whole-number one, or an additional fractional derivative on the wrong side of the equation. In what follows we refer to these two possibilities as to normal and modified forms of a fractional diffusion equation. An example of the normal form is an equation for superdiffusion with the Riesz fractional spatial derivative instead of the second derivative on the r.h.s. [2]. The situation with an additional derivative on the wrong side is exemplified by the standard fractional diffusion equation for subdiffusion [3]. We note that the two forms, i.e. the form with a Caputo derivative on the l.h.s. for subdiffusive processes, and a form with an additional spatial derivative on the l.h.s. for a superdiffusive process, are equivalent to the commonly used ones. This equivalence will be discussed in detail in the Section 2 of the present work, introducing the corresponding fractional operators. Many physical processes, however, lack power-law scaling, Eq.(.2), over the whole time-domain. These processes can not be characterized by a single scaling exponent β. Examples of processes lacking scaling include several cases of decelerating subdiffusion (e.g. the Sinai diffusion [4]) and decelerating superdiffusion (as exemplified by truncated Levy-flights [5]). Such processes can be described by derivatives of distributed order [6, 7], introduced by Caputo [8]. The equivalence between the different forms of fractional diffusion equations is lost in this case: Different forms of the distributed-order

3 3 fractional equations describe different situations. Thus, the equations with the distributed-order derivatives on the proper side describe processes getting more anomalous in course of the time (accelerating superdiffusion and decelerating subdiffusion) [9], while the equations with the additional distributed-order on the wrong side describe the situations getting less anomalous (decelerating superdiffusion and accelerating subdiffusion). As an example a special model with two fractional derivatives of different orders is used throughout the work. In Section 2 and 3 we discuss the four forms of fractional diffusion equations for the processes showing scaling behavior. These forms correspond to possible permutations of fractional temporal/spatial derivative on the proper / wrong side of the equation. Three of the four forms discussed below are well-known. The fourth one, to our best knowledge, was not previously discussed. We first turn to the forms pertinent to temporal fractional equations, i.e. to subdiffusion, and then to the superdiffusive case. Sections 4 and 5 are devoted to the distributed-order generalizations of the corresponding equations. The results are summarized in Section Two forms of time fractional diffusion equations 2.. Riemann-Liouville form The time fractional diffusion equation (TFDE) in the Riemann-Liouville form (RL - form), which in our terminology corresponds to a modified form of the fractional diffusion equation, reads [2]: t p(x,t) = K β D β t x2p(x,t), (2.) p(x,t = ) = δ(x), < β. Here K β is a positive constant, [K β ] = cm 2 /sec β, and D µ t is the Riemann- Liouville fractional derivative on the right semi-axis, which, for a sufficiently well-behaved function φ(t) is defined as follows: t D µ t φ = d d φ(τ) dt J µ φ = dτ Γ( µ) dt (t τ) µ, µ <, (2.2) where J α φ(t) = t Γ(α) dτ(t τ)α φ(t), t >, α R + is the Riemann- Liouville fractional integral of the order α. In what follows, we omit the subscript, for brevity. Applying the Laplace transform, φ(s) L {φ(t)} = 2 dte st φ(t),

4 4 and Fourier-transform, g(k) Φ {g(x)} = dxe ikx g(x) in succession, and using the Laplace transform of the Riemann-Liouville derivative (2.2), L {D µ t φ(t)} = sµ φ(s) (2.3) we get from Eq.(2.) the form of the Laplace-transformed characteristic function f(k,t) of the distribution p(x,t): f(k,s) = s β s β + K β k2. (2.4) 2.2. Caputo form The TFDE in the Caputo form (C - form, corresponding to a normal form of the fractional diffusion equation) is written as follows: β 2 t β p(x,t) = K β x2p(x,t), (2.5) p(x,) = δ(x), where K β is the same constant as in Eq.(2.), and the time fractional derivative of order β, < β < is understood in the Caputo sense [], β t β φ Dβ t φ(t) = J β d dt φ = Γ( β) t dτ(t τ) β d φ(τ). (2.6) dτ Here, the sequence of temporal integration and differentiation is reversed with respect to a Riemann-Liouville operator. Recalling the Laplace transform of the Caputo derivative, { d µ } φ L dt µ = s µ φ(s) s µ φ(), (2.7) < µ <, and making the Fourier-Laplace transform of Eq.(2.5) we again arrive at Eq.(2.4). Thus, both forms of TFDE are equivalent. In the literature, the third form of TFDE is mentioned, this one also using the Riemann - Liouville derivative [2]: 2 D β t p(x,t) = K β x 2p(x,t) + t β p(x,), (2.8) Γ( β)

5 5 < β <. The equivalence of the RL- and this third forms can be easily shown by applying t dt... to both sides of Eq.(2.) and using the Riemann - Liouville derivative of a constant, D β t = t β /Γ( β). The equivalence of the normal (C-) and the third forms can also be shown easily, if one uses the relation between the Riemann-Liouville and Caputo derivatives which can be obtained straightforwardly: D β t β φ(t) = x β φ + t β φ() (2.9) Γ( β) with < β <. Then, starting, e.g., from the C - form of TFDE, Eq.(2.5), and using Eq.(2.9) we immediately arrive at Eq.(2.8). In this paper we are interested in the mean squared displacement (MSD) given by x 2 (t) ) } = L {( 2 f k 2. (2.) k= From Eq.(2.4) we get: x 2 (t) = L { 2K β s β } = 2K βt β Γ( + β). (2.) 3. Two forms of space fractional diffusion equations 3.. The normal form The normal form of space fractional diffusion equation reads as: p t = K α p α x α, (3.) p(x,t = ) = δ(x), < α 2, where K α is a positive constant, [K α ] = cm α /sec, and the Riesz fractional derivative α / x α (we adopt here the notation introduced in [2]) is defined for a sufficiently well-behaved function f(x) through the Liouville - Weil derivatives []: d α d x α f(x) = { 2cos(πα/2) [ D α + + D α ] α dxĥf(x) d, (3.2) α = where D α ± are the left - and right side Liouville - Weil derivatives, D α +φ = D α φ = Γ(2 α) dx 2 d 2 d 2 Γ(2 α) dx 2 x x φ(ξ)dξ (ξ x) α, φ(ξ)dξ (x ξ) α (3.3)

6 6 for < α < 2, α (for α = D± = ±d/dx), and H is the Hilbert transform operator, Hφ = φ(ξ)dξ π x ξ. In Fourier space the operators of fractional derivatives have a simple form: where Since Φ(D± α φ) = dxexp(ikx)d± α φ = ( ik)α φ(k), (3.4) ( ( ik) α = k α exp απi ) 2 sign k. Φ(Hφ) = isign kφ(k) (3.5) then, with the use of Eqs.(3.2) - (3.5) we get the expression, which is valid for the Fourier transform of the Riesz fractional derivative for all values of α: ( d α ) φ Φ d x α = k α φ(k). (3.6) Applying the Fourier-transform to Eq.(3.), and noting Eq.(3.6), we get the characteristic function for the PDF of Lévy flights, f(k,t) = exp ( D k α t). (3.7) 3.2. The modified form Let us turn to the fractional equation for superdiffusion with the additional spatial derivative on its l.h.s. 2 α p x 2 α t = K 2 α x2p, (3.8) where K α is the same as in Eq.(3.). Note the minus sign in Eq.(3.8). This sign gets clear when turning to a Fourier representation: applying Fourier transform and using Eq.(3.6), we arrive at the characteristic function (3.7). Thus, both forms of the space fractional diffusion equation, Eqs.(3.) and (3.8), are equivalent. Since the mean square displacement diverges for Lévy flights, their anomalous nature can be characterized by a typical displacement δx of the diffusing particle, δx x q /q, q < α, (3.9)

7 7 which differs, of course, from the MSD discussed in Eq.(2.). To get the q-th moment we use the following expression [3]: x q = 2 πq Γ( + q)sin π 2 dk ( Ref(k,t)) k q. (3.) Inserting Eq.(3.7) into Eq.(3.) and changing to a new variable, ξ = K α k α t, we obtain and, thus Here dk... = (K αt) /α α dξ( e ξ )ξ q/α = (K αt) q/α q ( Γ q ), α (3.) x q = C(q,α)(K α t) q/α, q < α. (3.2) C(q,α) = 2 ( πq ) ( πq sin Γ( + q)γ q ). (3.3) 2 α Note that for a = q = 2 Eqs.(3.2), (3.3) give x 2 = 2Kt. 4. Distributed-order time fractional diffusion equations 4.. Distributed-order time fractional diffusion equation in the RL form The fractional diffusion equation with a distributed-order Riemann- Liouville derivative reads: p t = dβw(β)k(β)d β i 2 p x2, (4.) p(x,) = δ(x), where K(β) = Kτ β, [K] = cm 2 /sec, [τ] = sec, w(β) is a dimensionless non-negative function, which should fulfill dβw(β) =. If we set w(β) = δ(β β ), < β <, then we arrive at time fractional diffusion equation in the RL form, see Eq.(2.), where K β = Kτ β. We now prove that the solution of Eq.(4.) is a PDF. The derivation here follows the method used in [4]. Its aim is to show that the random process whose PDF obeys Eq.(4.) is subordinated to the Wiener process. Making a Fourier-Laplace transform of Eq.(4.) and using Eq.(2.3) we get f(k,s) = ( si RL I RL + ), (4.2) k2 Kτ

8 8 where I RL (sτ) = We rewrite Eq.(4.2) as follows: f(k,s) = si RL where duexp [ u ( I RL + k2 Kτ )] = dβ(sτ) β w(β). (4.3) [ G RL (u,s) = si RL (sτ) exp u ] I RL (sτ) due uk2 Kτ GRL (u,s) (4.4) (4.5) is a Laplace transform of a function whose properties will be specified below. Now, with the help of Eqs.(4.4) and (4.5) the PDF p(x,t) can be written as dk ds p(x,t) = 2π e ikx 2πi est due uk2 Kτ GRL (u,s) = = Br du e x2 /4uKτ 4πuKτ G RL (u,t) (4.6) where Br denotes the Bromwich integration contour. In order to prove the positivity of p(x,t) we demonstrate that the function G RL (u,t) is the PDF providing the subordination transformation from time scale t to time scale u, that is, G RL (u,t) is positive and normalized with respect to u for any t. At first we demonstrate normalization, { } { } dug RL (u,t) = L du e u/i RL = L =. (4.7) si RL s Now, to prove the positivity of function G RL (u,t), according to the Bernstein theorem [5], it is enough to show that G RL (u,s), Eq.(4.5), is completely monotonic as a function of s on positive real axis, i.e., it is positive and the signs of its derivatives alternate. We do this for the special case when w(β) = B δ(β β ) + B 2 δ(β β 2 ) (4.8) with < β < β 2, B >, B 2 >, B + B 2 =. This choice allows us to show in a simple way the property of the diffusive behavior governed by distributed-order diffusion equations. This is why this case is repeatedly discussed in our article. Inserting Eq.(4.8) into Eq.(4.3) we have I RL (sτ) = b s β + b 2 s β 2 (4.9)

9 9 where b = B /τ β, b 2 = B 2 /τ β 2. From Eq.(4.5) it follows that G RL (u,s) is a product of two functions: GRL (u,s) = φ (s)φ 2 (u,s) with [ φ = si RL (sτ), φ 2 = exp u ]. (4.) I RL (sτ) We will prove that both functions φ and φ 2 are completely monotonic, and therefore, as a product, G RL (u,s) is completely monotonic too. Let us start from φ. It can be rewritten as φ = φ(h(s)), where φ(y) = /y is a completely monotonic function, and h(s) = b s β + b 2 s β 2 is a positive function with a completely monotonic derivative (this is evident by direct inspection). Therefore, φ (s) is completely monotonic according to Criterion 2 in [5], Chapter XIII, 4. Now turn to φ 2. Again, it can be rewritten as φ 2 = φ(ψ(s)), where φ(y) = exp( uy) is a completely monotonic function. According to the same Criterion, it is enough to show that the positive function ψ(s) = I RL (sτ) = b s β + b2 s β 2 = sβ2 b 2 + b b2 s β 2 β (4.) possesses completely monotonic derivative. By differentiating Eq.(4.) we get ψ (s) = sβ 2 b b 2 + b sβ2 b 2 (β 2 β )s β 2 β b2 s β ( ) 2 β b 2 = 2 + b b2 s β 2 β = β 2s β 2 b 2 + b b2 s β 2 β + b β b 2 β 2 s β 2 β + b b2 s β 2 β. (4.2) Denoting ξ = β 2 β, < ξ <, we note that this function is a product of three completely monotonic functions, and thus is itself completely monotone, too. Indeed, (i) the first function is s β 2 being a negative power of s; (ii) the second one is completely monotonic since it has the form g(h(s)), with g(y) = /( + y) being completely monotone, and h(ξ) = s ξ is a positive function with a completely monotone derivative; (iii) the last function has the same form with g(y) = ( + cy)/( + dy), where < c < d with c = b β /b 2 β 2, d = b /b 2. The function g(y) is completely monotonic since its n-th derivative (obtained using the Leibnitz rule) has the form g (n) (y) = ( )n (c d)n!d n ( + dy) n+, and the signs of the successive derivatives alternate. Thus, we have proved that ψ (s), Eq.(4.2), is a completely monotonic function. Therefore the

10 exponential φ 2, Eq.(4.), is a completely monotonic function, too, and the whole function, G RL (u,s), is completely monotone as a product of two completely monotonic functions, φ and φ 2. Therefore, the function G RL (u,t) is a PDF and, according to Eq.(4.6), the function p(x,t) is a PDF, too. Thus, we have proved that the solution of the distributed order time fractional diffusion equation in the Riemann-Liouville form is a PDF. We are interested in the MSD which is given by x 2 (t) ) } { } = L {( 2 f k 2 = KτL 2IRL (sτ). (4.3) s k= For the particular case (4.8) we get, by inserting Eq.(4.9) into Eq.(4.3) and making an inverse Laplace transform, x 2 (t) = 2κ Γ( + β ) ( ) t β + τ 2κ 2 Γ( + β 2 ) ( ) t β2 (4.4) τ where κ = B Kτ and κ 2 = B 2 Kτ. Since < β < β 2, at small times the first term in the r.h.s. of Eq.(4.4) prevails, whereas at large times the second one dominates. Thus, the overall behavior corresponds to accelerating subdiffusion Distributed-order time fractional diffusion equation in the C - form The distributed-order time fractional diffusion equation in the normal form can be written as dβτ β w(β) β p t β = K 2 p x2, (4.5) p(x, ) = δ(x), where τ is a positive constant representing some characteristic time of the problem (vide infra), [τ] = sec, K is the diffusion coefficient, [K] = cm 2 /sec, w(β) is a dimensionless non-negative function, and the time fractional derivative of order β, < β < is understood in the Caputo sense, Eq.(2.6). Note the difference between Eq.(4.5) and Eq.(4.). If we set w(β) = δ(β β ), < β <, we arrive at the normal form of a fractional diffusion equation, Eq.(2.5), with β = β and K β = Kτ β. Let us prove that the solution of Eq.(4.5) is a PDF. Applying the Laplace- and Fourier-transforms in succession, we get: f(k,s) = s I C (sτ) I C (sτ) + k 2 Kτ, (4.6)

11 where I C (sτ) = We rewrite Eq.(4.6) in the form analogous to Eq.(4.4), f(k,s) = I C s dβ(sτ) β w(β). (4.7) due u[ic+k2 Kτ] = due uk2 Kτ GC (u,s) (4.8) where G C (u,s) = I C(sτ) e ui C(sτ) (4.9) s is the Laplace transform of a function G C (u,t) whose properties will be specified below. Now, p(x,t) can be written in the form analogous to Eq.(4.6): dk ds p(x,t) = 2π e ikx 2πi est due uk2 Kτ GC (u,s) = = Br du e x2 /4uKτ 4πuKτ G C (u,t) (4.2) Similar to the RL case, the function G C (u,t) is the PDF providing the subordination transformation, from time scale t to time scale u. Indeed, at first we note that G C (u,t) is normalized with respect to u for any t. Using Eq.(4.9) we get { dug C (u,t) = L du I } { } RL s e u/i C = L =. (4.2) s To prove the positivity of G C (u,t) it is enough to show that its Laplace transform is completely monotone on the positive real axis [5]. The last statement is proved by noting that it is a product of two completely monotonic functions, I C /s and exp( ui C ). We again demonstrate this for the particular choice of w(β), see Eq.(4.8). For this choice we obtain from Eq.(4.7) I C (s) = b s β + b 2 s β 2, (4.22) where b = B τ β, b 2 = B 2 τ β 2. Then, I C /s is completely monotone as a sum of the negative powers of s, and I C itself is a positive function with a completely monotone derivative. Thus, exp( ui C ) is also completely monotone. Thus, we have proved that is completely monotone, and G RL (u,t) is a PDF, according to the Bernstein theorem. We are interested in the behavior of the MSD, i.e. in second moment of the PDF. Using Eq.(4.6), we have x 2 (t) = L {( 2 f k 2 ) k= } = 2KτL { }. (4.23) si C (sτ)

12 2 Inserting Eq.(4.22) into Eq.(4.23) one obtains: x 2 (t) { } = 2KτL s (b s β + b2 s β = 2Kτ L 2 ) b 2 { s β b b 2 + s β 2 β (4.24) Recalling the Laplace transform of the generalized Mittag-Leffler function E µ,ν (z), µ >, ν >, which can be conveniently written as [] we get from Eq.(4.24): L { t ν E µ,ν ( λt µ ) } = sµ ν s µ + λ, Res > λ /µ, (4.25) x 2 (t) = 2Kτ b 2 t β 2 E β2 β,β 2 + } ( b ) t β 2 β. (4.26) b 2 To obtain asymptotics at small t, we use an expansion, which is, in fact the definition of E µ,ν (z), see [6], Ch.XVIII, Eq.(9): z n E µ,ν (z) = Γ(µn + ν), (4.27) n= which yields in the main order for the MSD x 2 (t) = 2Kτ B 2 Γ(β 2 + ) ( t τ ) β2 t β 2. (4.28) For large t we use the following expansion valid on the real negative axis, see [6], Ch.XVIII, Eq.(2): E µ,ν (z) = which yields N n= z n ( Γ( µn + ν) + O z N), z (4.29) x 2 (t) = 2Kτ B Γ(β + ) ( ) t β t β. (4.3) τ Since β < β 2, we have retarded subdiffusion. It is worthwhile to note that the distributed-order equation (4.5) with the Caputo derivative always describes retarding, or slowing-down, sub-diffusive processes. Indeed, it is clearly seen from Eqs.(4.22), (4.23) that, for w(β) given by Eq.(4.8), for large s (short times) I(sτ) behaves as s β 2, so that, due to the Tauberian theorem, the MSD behaves as t β 2, while at long times (small s) I(sτ) s β and, respectively, x 2 t β, β < β 2. Thus, Eq.(4.5) cannot describe the accelerating sub-diffusive process. The same arguments applied to Eq.(4.) with the Riemann - Liouville derivative demonstrate that Eq.(4.) cannot describe retarded sub-diffusive processes..

13 3 5. Distributed-Order Space Fractional Diffusion Equation 5.. The normal form We write the distributed-order space fractional diffusion equation in the normal form as p t = 2 dαk(α) α p x α, p(x,) = δ(x) (5.) where K(α) is a (dimensional) function of the order of the derivative α, and the Riesz space fractional derivative α / x α is given by its Fourier transform, see Eq.(3.6). Setting K(α) = K α δ(α α ) we arrive at Eq.(3.). In the general case K(α) can be represented as K(α) = l α 2 Kw(α), (5.2) where l and K are dimensional positive constants, [l] = cm, [K] = cm 2 /sec and w is a dimensionless non-negative function of α. The equation for the characteristic function of Eq.(5.) has the solution f(k,t) = exp { Kt l 2 Note that the normalization condition, 2 is fulfilled. Let us consider the simple particular case, dαw(α)( k l) α }. (5.3) dxp(x,t) = f(k =,t) = (5.4) w(α) = A δ(α α ) + A 2 δ(α α 2 ), (5.5) where < α < α 2 2, A >, A 2 >. Inserting Eq.(5.5) into Eq.(5.3) we have f(k,t) = exp { a k α t a 2 k α 2 t} (5.6) where a = A K/l 2 α, a 2 = A 2 K/l 2 α 2. The characteristic function (5.6) is a product of two characteristic functions of Lévy stable PDFs with the Lévy indexes α, α 2, and the scale parameters a /α and a /α 2 2, respectively. Therefore, the inverse Fourier transformation of Eq.(5.6) gives the PDF which is a convolution of two stable PDFs, p(x,t) = a α a α 2 2 t α α 2 ( dx x x L α, (a t) α ) L α2, ( x (a 2 t) α 2 ), (5.7)

14 4 where L α, is the PDF of a symmetric Lévy stable law given by its characteristic function ˆL α, (k) = exp ( k α ). (5.8) The PDF given by Eq.(5.7) is, obviously, positive, as the convolution of two positive PDFs. The PDF will be also positive, if the function A(α) is represented as a sum of N delta-functions multiplied by positive constants, N is a positive integer. Moreover, if A(α) is a continuous positive function, then discretizing the integral in Eq.(5.) by a Riemann sum and passing to the limit we can also conclude on the positivity of the PDF in the general case. Let us consider the q-th moment, Eq.(3.), at small and large t. After inserting Eq.(5.6) into Eq.(3.) we get where x q = 2 π Γ( + q)sin ( πq 2 ) ( dk e tφ(k)) k q, (5.9) φ(k) = a k α + a 2 k α 2. (5.) In order to get the q-th moment at small t, we expand exp ( a k α t) in power series with subsequent integration over k. As the result we have the following expansion valid at q < α : x q = 2 πq (a 2t) q/α 2 sin + ( Γ q ( πq 2 α 2 ) ) Γ( + q)γ n= ( ) n+ α 2 n! ( qα2 ) (5.) nα a n a α 2 2 Γ ( nα q α 2 )t n ( α α 2 ) for t. We recall that the radius R of convergence of a series n= c nx n is determined from the equation /R = lim n c n /n. Then, by using the Stirling formula, one can see that in case of Eq.(5.) c n /n n α /α 2 for n, and therefore expansion (5.) is valid at all t. The leading term of a series (5.) gives for the characteristic displacement at small t, δx x q /q t /α 2. (5.2) In order to get the q-th moment at large t we integrate by parts the right hand side of Eq.(5.9): x q = 2 ( πq ) t π Γ( + q)sin J(t), (5.3) 2 q

15 5 where with J(t) = ζ(k) = φ (k) k q dkζ(k)e tφ(k) (5.4) = a α k α + a 2 α 2 k α2 k q+. Since at small k one has φ(k) a k α, ζ(k) a α k α q, then at large t we have exactly the case of Laplace asymptotic integral, see [7], Chapter III, Theorem 7.. For J(t) one immediately obtains J(t) = a q/α t q/α Γ( q/α ) and x q 2 ( πq (a t) q/α πq ) ) sin Γ( + q)γ ( qα, q < α. (5.5) 2 For the characteristic displacement at large t we get δx x q /q t /α. (5.6) Therefore, at small times the characteristic displacement grows as t /α 2, whereas at large times it grows as t /α. Since α < α 2, one encounters an accelerated superdiffusion The modified form Let us now consider the conjugated form of the equation, namely one with an additional distributed-order fractional operator in the l.h.s.: 2 dαw(α)l 2 α 2 α p x 2 α t = K 2 p dx2, p(x,) = δ(x), (5.7) where l, K and w have the same meaning as in Eq.(5.2). Setting w(α) = δ(α α ) we arrive at Eq.(3.8) where K α K α = Kl α 2. The equation for the characteristic function of the solution of Eq.(5.7) reads: [ ] Kt/l 2 f(k,t) = exp 2 dαw(α)( k (5.8) l)α (compare with Eq.(5.3)). Note that the normalization condition, Eq.(5.4), is fulfilled. As we did it throughout the present article, let us consider a particular case of w(α) = A δ(α α ) + A 2 δ(α α 2 ), (5.9)

16 6 < α < α 2 2. From Eq.(5.8) we get for the characteristic function, [ ] t f(k,t) = exp, (5.2) a k α + a 2 k α 2 where a,2 = A,2 l 2 α,2 /K. The proof of non-negativity of the PDF given by inverse Fourier transform of Eq.(5.2) follows along the same lines as for the accelerating subdiffusion case. The q-th moment is given by Eq.(5.9), where φ(k) = a k α + a2 k α 2. (5.2) We insert Eq.(5.2) into Eq.(5.9) and pass to a new variable ξ = tk α /α. For the integral over k we get ( ) t q/α dk... = (5.22) a α dξξ q/α exp ξ + a α2 2 ( α a ξ t ) α 2 α. a At small t we can neglect the term with t in the denominator in square brackets. Therefore, ( ) x q t q/α = C(q,α ), q < α, (5.23) α where C(q,α ) is given by Eq.(3.3), and δx x q /q t /α (5.24) for t. In order to get the q-th moment at large t, we again use the Laplace method. Turn to Eq.(5.3), where φ(k) is given by Eq.(5.2) and ζ(k) = φ (k) k q = k +q a α k α + a 2 α 2 k α2 (a k α + a2 k α 2 ) 2. (5.25) Since at small k we have φ(k) k α 2, ζ(k) (α 2 /a 2 ) k α2 q, then at large t we again have exactly the case of Laplace asymptotic integral [7]. For J(t) we have immediately J(t) = a q/α 2 t q/α2 Γ( q/α 2 ), and ( ) x q t q/α C(q,α 2 ), q < α, (5.26) α 2

17 7 where C(q,α) is given by Eq.(3.3), and δx x q /q t /α 2 (5.27) for t. Therefore, at small times the characteristic displacement grows as t /α, whereas at large times it grows as t /α 2. Since α < α 2, we encounter retarding superdiffusion. 6. Conclusions Distributed-order diffusion equations generalize the approach based on fractional diffusion equations to processes lacking temporal scaling. The typical forms of fractional diffusion equations can be classified with respect to the position of the fractional derivative in time/coordinate instead of or in addition to the first and second derivatives in the classical Fick s form. In the present paper we considered the corresponding forms of distributed-order fractional diffusion equations and elucidate the effects described by different classes of such equations. We show that equations with the distributed-order fractional operator replacing the corresponding wholenumber derivative describe processes getting more anomalous in course of the time, i.e. the accelerating superdiffusion or retarded subdiffusion. On the opposite, equations with additional fractional operators on the wrong side of the diffusion equation describe processes getting less anomalous, i.e. retarded superdiffusion and accelerating subdiffusion. 7. Acknowledgements The work was supported by an INTAS grant. IMS gratefully acknowledges the support by the Fonds der Chemischen Industrie. REFERENCES [] I.M. Sokolov, J. Klafter, A. Blumen, Physics Today 55, 48 (22). [2] R. Metzler and J. Klafter, Phys. Rep. 339, (2). [3] R. Metzler, E. Barkai, J. Klafter, Phys. Rev. Lett. 82, 3563 (999). [4] Ya. G. Sinai, Theor. Prob. Appl. 27, 256 (982). [5] R.N. Mantegna, H.E. Stanley, Phys. Rev. Lett (994). [6] A.V. Chechkin, J. Klafter and I.M. Sokolov, Europhys. Lett., (23). [7] I.M. Sokolov, A.V. Chechkin, J. Klafter, preprint cond-mat/ [8] M. Caputo, Fract. Calc. Appl. Anal. 4, 42 (2).

18 8 [9] A. V. Chechkin, R. Gorenflo, I. M. Sokolov, Phys. Rev. E66, 4629 (22). [] S.G.Samko, A.A. Kilbas, O.I. Marichev, Integraly i proizvodnye drobnogo poryadka i nekotorye ikh prilozheniya, Nauka i technika, Minsk, 987 (in Russian); Engl. Transl. Fractional Integrals and Derivatives, Theory and Applications, Gordon and Breach, Amsterdam, 993. [] R. Gorenflo and F. Mainardi, in A. Carpinteri, F. Mainardi (eds.) Fractals and Fractional Calculus in Continuum Mechanics, Springer Verlag, Wien and New York, 997, p.223. [2] A. I. Saichev, G. M. Zaslavsky, Chaos 7, 753 (997). [3] V.M. Zolotarev, Odnomernye ustoichivye raspredeleniya, Mir, Moskva, 983 (in Russian); Engl. Trans. One Dimensional Stable Distributions. Amer. Math. Soc., Providence, RI, 986. [4] I.M. Sokolov, Phys. Rev. E63, 56, (2). [5] W. Feller, An Introduction to Probability Theory and Its Applications, vol.ii, John Wiley and Sons, Inc., New York, 2-nd ed., 97. [6] A. Erdélyi, ed. Bateman Manuscript Project. Higher Transcendental Functions, vol.iii. McGraw-Hill, Inc., New York (955). [7] F.W.J. Olver, Asymptotics and Special Functions, Academic Press, New York, 974.