Next Article in Journal / Special Issue
Asymptotic Expansions of Fractional Derivatives andTheir Applications
Previous Article in Journal / Special Issue
Fractional Diffusion in Gaussian Noisy Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution of Generalized Space-Time Fractional Cable Equation

1
Department of Mathematics and Statistics, Jai Narain Vyas University, Jodhpur 342004, India
2
Department of Mathematics, University of Rijeka, Radmile Matejcic 2, 51000 Rijeka, Croatia
3
Faculty of Natural Sciences and Mathematics, Institute of Mathematics, Saints Cyril and Methodius University, 1000 Skopje, Macedonia
4
Radiation Safety Directorate, Partizanski odredi 143, P.O. Box 22, 1020 Skopje, Macedonia
5
Max Planck Institute for the Physics of Complex Systems, Noethnitzer Strasse 38, 01187 Dresden, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 7 March 2015 / Accepted: 2 April 2015 / Published: 9 April 2015
(This article belongs to the Special Issue Recent Advances in Fractional Calculus and Its Applications)

Abstract

:
In this paper, we consider generalized space-time fractional cable equation in presence of external source. By using the Fourier-Laplace transform we obtain the Green function in terms of infinite series in H-functions. The fractional moments of the fundamental solution are derived and their asymptotic behavior in the short and long time limit is analyzed. Some previously obtained results are compared with those presented in this paper. By using the Bernstein characterization theorem we find the conditions under which the even moments are non-negative.

Graphical Abstract

1. Introduction

The cable equation has been used to model electrotonic properties of neuronal dendrities [1]. In the work of Baer and Rinzel [2], for the first time, theoretical study of wave propagation mediated by dendritic spine is performed, where the dendritic tissue is modeled with the classical cable equation. The classical cable equation which models the membrane potential V = V(x,t) along the axial x-direction of a dendrite with diameter d, relative to the resting membrane potential Vrest, is given by
c m r m d V ( x * , t * ) d t * = d r m 4 r L d 2 V ( x * , t * ) d x * 2 V V r e s t r m + r m i e ( x * , t * )
where r m denotes the specific membrane resistance, r L is the longitudinal resistivity, c m denotes the membrane capacitance per unit area, and i e is the external injected current per unit area. The product τ = c m r m is the time constant for the dendrite.
The cable equation can be derived from the Nernst-Planck equation for electrodiffusive motion of ions [3]. There are generalizations of the cable Equation (1) where the non-local nature of the diffusive transport of living cells is taken into consideration. In the excellent works of Henry, Langlands and Wearne [4,5,6], time fractional cable models for spiny neuronal dendrities are introduced and investigated in detail. They derived the time fractional cable equations from fractional Nernst-Planck equations, and investigated the electrotonic effects of the trapping properties of spines [4]. In the work of Li and Deng [7], space-time fractional cable equation is introduced where the spatial non-local effects are taken into account by modifying the Ohm’s law to fractional one. Thus, Li and Deng considered the following space-time fractional cable equation in the infinite domain < x <
V ( x , t ) t = D t 1 α ( Δ μ V ( x , t ) ) λ 2 D t 1 β ( V ( x , t ) ) + f ( x , t )
where x and t are dimensionless parameters (see relation (10) from [7]), 0 < α 1 , 0 < β 1 , λ = d r m 4 r L is the space dendrite for the dendrite, i.e., the cable, f ( x , t ) = λ 2 D t 1 β ( V r e s t + r m i e ( x , t ) ) is the external source (external injected current). In Equation (2), D t 1 ν is the Riemann-Liouville (R-L) fractional derivative of order 1 ν [8,9]:
D t 1 ν f ( t ) = 1 Γ ( ν ) d d t 0 t f ( u ) ( t u ) 1 ν d u ,       0 < ν < 1
and Δ μ is the Riesz fractional operator defined by [9,10]
Δ μ = 1 2 cos ( μπ 2 ) ( μ x μ + μ ( x ) μ ) ,       1 < μ 2
where μ x μ and μ ( x ) μ are the left and right R-L fractional derivatives of order μ defined respectively by [8,9]
μ x μ f ( x ) = 1 Γ ( n μ ) n x n x f ( u ) ( x u ) μ n + 1 d u ,      n 1 < μ n ,       n N
μ ( x ) μ f ( x ) = 1 Γ ( n μ ) n x n x f ( u ) ( u x ) μ n + 1 d u ,      n 1 < μ n ,       n N
Here we note that the Riesz fractional operator (4) can be represented as ( Δ ) μ / 2 , and that it is a generalization of the Laplacian which is obtained for μ = 2 . It represents pseudodifferential operator with Fourier symbol | k | μ , i.e., F [ Δ μ f ( x ) ] ( k ) = | k | μ F [ f ( x ) ] ( k ) [10].
The boundary conditions for the cable Equation (2) are given by
lim x ± V ( x , t ) = 0
and the initial condition by
V ( x , 0 ) = g ( x ) .
Here we note that the fractional cable equation with nonlocal boundary conditions is recently considered in [11].
Equation (2) with λ = 0 and f ( x , t ) = 0 is the space-time fractional diffusion equation which describes the competition between the subdiffusion and Levy flights, and which can be obtained from the continuous time random walk (CTRW) theory for broad distribution of waiting times and long-tailed distribution of jump lengths [12]. Fractional diffusion equations in presence of external force were also of interest in many recent papers [12,13,14,15]. Analytical solutions of fractional differential equations and numerical methods of fractional cable equation have been considered by many authors [16,17,18,19,20,21,22,23,24], to name but a few.
In our work, we further generalize Equation (2) by introducing time fractional derivative of Caputo form of order 0 < γ 1 , defined by [8,9]
D t γ C f ( t ) = 1 Γ ( n λ ) 0 t f ( n ) ( u ) ( t u ) γ n + 1 d u ,      n 1 < γ n ,       n N
which has been used in different contexts.
The application of different forms fractional derivatives [8,9,10,25,26] have been found to be a useful tool for modeling systems with memory, due to the power-law memory kernel contained in the convolution integral. It is shown that the fractional derivatives appear in the diffusion equation describing anomalous diffusive processes where the mean square displacement has a power-law dependence on time [12], and that they can be used for modeling aquifer problems [13,27], non-exponential relaxation problems [28,29], etc. Furthermore, time fractional derivatives are equivalent to infinitesimal generators of generalized time fractional evolutions that arise in the transition from microscopic to macroscopic time scales [28,29]. The transition from first order time derivative to the fractional one arises in many physical problems as shown by Hilfer [25,30,31]. Contrary to the time fractional derivatives, the Riesz space fractional derivative has benn shown to represent suitable tool for modeling Levy flights for which the second moment diverges, and, thus, the fractional moments should be calculated in order to analyze the superdiffusive behavior of the particle [12,14].
Thus, the generalization of the cable equation by using time fractional derivatives and Riesz space fractional derivative take into consideration temporal memory effect and long range spatial correlations.
This paper is organized as follows. In Section 2, we consider generalized space-time fractional cable equation with a source term. The fundamental solution (Green function) is represented in terms of infinite series in H-functions. The fractional moments of the fundamental solution are derived in Section 3, and their asymptotic behaviors are analyzed. Summary is given in Section 4.

2. Solution of the problem

In this paper, we consider the following generalized space-time fractional cable equation
τ γ D t γ C V ( x , t ) = D t 1 α ( Δ μ V ( x , t ) ) λ 2 D t 1 β ( V ( x , t ) ) + f ( x , t )
where τ γ is a time parameter introduced for dimensional purposes. Without loss of generality we set τ γ = 1 . A few words about Equation (10) are in order. Let us apply R-L time fractional integral of order 0 < γ 1 defined by [8,9]
I t γ f ( t ) = 1 Γ ( γ ) 0 t f ( u ) ( t u ) 1 γ d u
to Equation (10). By using the relation I t γ D t γ C f ( t ) = f ( t ) f ( 0 + ) for 0 < γ 1 [8,9], and then applying first time derivative, we rewrite Equation (10) in form (2), i.e.,
d V ( x , t ) d t = D t 1 γ D t 1 α ( Δ μ V ( x , t ) ) λ 2 D t 1 γ D t 1 β ( V ( x , t ) ) + D t 1 γ f ( x , t )
where we use the relation d d t I t γ f ( t ) = D t 1 γ f ( t ) for 0 < γ 1 [8,9]. In order Equation (12) to be generalization of the classical cable equation we need 1 < γ + α 2 and 1 < γ + β 2 . From Equation (12) we see that we have fractional derivative of the source term, i.e., D t 1 γ f ( x , t ) = λ 2 D t 1 γ D t 1 β ( V r e s t + r m i e ( x , t ) ) . This additional parameter 0 < γ 1 may be used for better fitting experimental data.
Theorem 1. The solution of Equation (10) with boundary conditions (7) and initial condition (8) is given by
V ( x , t ) = G ( x ξ , t ) g ( ξ ) d ξ + 0 t G f ( x ξ , t τ ) f ( ξ , τ ) d τ d ξ ,
where
G ( x , t ) = 1 μ | x | π n = 0 ( λ 2 ) n t ( γ + β 1 ) n n ! H 2 , 3 2 , 1 [ | x | 2 t ( α + γ 1 ) / μ | ( 1 , 1 μ ) , ( 1 + ( γ + β 1 ) n , γ + α 1 μ ) ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ] ,
G f ( x , t ) = 1 μ | x | π n = 0 ( λ 2 ) n t ( γ + β 1 ) n + γ 1 n ! H 2 , 3 2 , 1 [ | x | 2 t ( α + γ 1 ) / μ | ( 1 , 1 μ ) , ( γ + ( γ + β 1 ) n , γ + α 1 μ ) ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ]
Re ( μ ) > 0 , and H p , q m , n [ z | ( a p , A p ) ( b q , B q ) ] is the Fox H-function [32,33,34].
Proof. By applying Laplace transform with respect to the time t to Equation (10) we obtain
s γ V ^ ( x , s ) s γ 1 g ( x ) = s 1 α Δ μ V ^ ( x , s ) λ 2 s 1 β V ^ ( x , s ) + f ^ ( x , s )
where we use the initial condition (7) and the following formula for the Laplace transform of the Caputo fractional derivative [8,9]
L [ D t γ V ( x , t ) ] ( x , s ) = s γ V ^ ( x , s ) V ( x , 0 + )
where L [ V ( x , t ) ] ( x , s ) = V ^ ( x , s ) , L [ f ( x , t ) ] ( x , s ) = f ^ ( x , s ) , and the Laplace transform of the R-L fractional derivative defined by [8,9]
L [ D t 1 ν V ( x , t ) ] ( x , s ) = s 1 ν V ^ ( x , s ) [ D t ν 1 V ( x , t ) ] t = 0
Here we suppose that [ D t ν 1 V ( x , t ) ] t = 0 = 0 , ( ν = { α , β } ) (please see [4]). The Fourier transform of (16) yields
s γ V ^ ˜ ( k , s ) s γ 1 g ˜ ( k ) = s 1 α | k | μ V ^ ˜ ( k , s ) λ 2 s 1 β V ^ ˜ ( k , s ) + f ^ ˜ ( k , s )
where F [ V ^ ( x , s ) ] ( k , s ) = V ^ ˜ ( k , s ) is the Fourier transform of V ^ ( x , s ) , g ˜ ( k ) = F [ g ( x ) ] ( k ) , f ^ ˜ ( k , s ) = F [ f ^ ( x , s ) ] ( k , s ) , and we use that F [ Δ μ f ( x ) ] ( k ) = | k | μ f ˜ ( k ) [10]. Thus, we obtain
V ^ ˜ ( k , s ) ( k , s ) = s γ 1 s γ + s 1 α | k | μ + λ 2 s 1 β g ˜ ( k ) + 1 s γ + s 1 α | k | μ + λ 2 s 1 β f ^ ˜ ( k , s ) .
Here
G ^ ˜ ( k , s ) = s γ 1 s γ + s 1 α | k | μ + λ 2 s 1 β = s γ + α 2 s γ + α 1 + | k | μ + λ 2 s α β
is the Green function (solution for g ( x ) = δ ( x ) i.e., g ˜ ( k ) = 1 , and f ( x , t ) = 0 ), and
G ^ ˜ f ( k , s ) = 1 s γ + s 1 α | k | μ + λ 2 s 1 β = s α 1 s γ + α 1 + | k | μ + λ 2 s α β .
Thus, the solution can be represented in a form given by (13). Expanding the right hand side of (21) and (22) in a power series in s (see the approach given in [8,9]), we arrive at
G ^ ˜ ( k , s ) = n = 0 ( λ 2 ) n s ( α β ) n + γ + α 2 ( s γ + α 1 + | k | μ ) n + 1 ,
G ^ ˜ f ( k , s ) = n = 0 ( λ 2 ) n s ( α β ) n + α 1 ( s γ + α 1 + | k | μ ) n + 1 .
Using the Laplace transform formula for the three parameter Mittag-Leffler (M-L) function E α , β δ ( z ) = n = 0 ( δ ) n Γ ( α n + β ) z n n ! ( ( δ ) n = Γ ( δ + n ) / Γ ( δ ) is the Pochhammer symbol) [35]:
L [ t β 1 E α , β δ ( a t α ) ] ( s ) = s α δ β ( s α + a ) δ ,       ( Re ( s ) > | a | 1 / α )
it yields
G ˜ ( k , t ) = n = 0 ( λ 2 ) n t ( γ + β 1 ) n E γ + α 1 , ( γ + β 1 ) n + 1 n + 1 ( | k | μ t γ + α 1 ) ,
G ˜ f ( k , t ) = n = 0 ( λ 2 ) n t ( γ + β 1 ) n + γ 1 E γ + α 1 , ( γ + β 1 ) n + γ n + 1 ( | k | μ t γ + α 1 ) .
The convergence of series in three parameter M-L functions of form (26) and (27) is shown in [36]. For various relations and applications of the M-L functions we refer to [37]. If we use the relation between the three parameter M-L function and the H-function [33]
E α , β δ ( z ) = 1 Γ ( δ ) H 1 , 2 1 , 1 [ z | ( 1 δ , 1 ) ( 0 , 1 ) , ( 1 β , α ) ]
we obtain
G ˜ ( k , t ) = n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n H 1 , 2 1 , 1 [ | k | μ t γ + α 1 | ( n , 1 ) ( 0 , 1 ) , ( ( γ + β 1 ) n , γ + α 1 ) ] ,
G ˜ f ( k , t ) = n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n + γ 1 H 1 , 2 1 , 1 [ | k | μ t γ + α 1 | ( n , 1 ) ( 0 , 1 ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 ) ] .
The inverse Fourier transform of above expressions
G ( x , t ) = 1 2 π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n e i k x H 1 , 2 1 , 1 [ | k | μ t γ + α 1 | ( n , 1 ) ( 0 , 1 ) , ( ( γ + β 1 ) n , γ + α 1 ) ] d k ,
G f ( x , t ) = 1 2 π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n + γ 1 e i k x H 1 , 2 1 , 1 [ | k | μ t γ + α 1 | ( n , 1 ) ( 0 , 1 ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 ) ] d k ,
can be found by using the cosine transform of the H-function with ρ = 1 [33]
0 k ρ 1 cos ( k x ) H p , q m , n [ a k μ | ( a p , A p ) ( b q , B q ) ] d k = 2 ρ 1 π x ρ H p + 2 , q m , n + 1 [ a ( 2 x ) μ | ( 2 ρ 2 , μ 2 ) , ( a p , A p ) , ( 1 ρ 2 , μ 2 ) ( b q , B q ) ]
where Re ( ρ ) + μ min 1 j m Re ( b j B j ) > 0 , Re ( ρ ) + μ max 1 j n Re ( a j 1 B j ) < 1 , Ω = j = 1 n A j j = n + 1 p A j + j = 1 m B j j = m + 1 q B j , Ω > 0 , | arg a | < πΩ 2 . Therefore, we find
G ( x , t ) = 1 π | x | n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n H 3 , 2 1 , 2 [ ( 2 | x | ) μ t γ + α 1 | ( 1 2 , μ 2 ) , ( n , 1 ) , ( 0 , μ 2 ) ( 0 , 1 ) , ( ( γ + β 1 ) n , γ + α 1 ) ] ,
G f ( x , t ) = 1 π | x | n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n + γ 1 H 3 , 2 1 , 2 [ ( 2 | x | ) μ t γ + α 1 | ( 1 2 , μ 2 ) , ( n , 1 ) , ( 0 , μ 2 ) ( 0 , 1 ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 ) ] .
If we use the following property of the H-function [33]
H p , q m , n [ z | ( a p , A p ) ( b q , B q ) ] = 1 σ H p , q m , n [ z 1 / σ | ( a p , A p / σ ) ( b q , B q / σ ) ] ,       σ > 0
then (32) and (33) yield
G ( x , t ) = 1 μ π | x | n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n H 3 , 2 1 , 2 [ 2 t ( γ + α 1 ) / μ | x | | ( 1 2 , 1 2 ) , ( n , 1 μ ) , ( 0 , 1 2 ) ( 0 , 1 μ ) , ( ( γ + β 1 ) n , γ + α 1 μ ) ] ,
G f ( x , t ) = 1 μ π | x | n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) n + γ 1 H 3 , 2 1 , 2 [ 2 t ( γ + α 1 ) / μ | x | | ( 1 2 , 1 2 ) , ( n , 1 μ ) , ( 0 , 1 2 ) ( 0 , 1 μ ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 μ ) ] .
If we further use the following transformation formula for the H-function [33]
H p , q m , n [ z | ( a p , A p ) ( b q , B q ) ] = H q , p n , m [ 1 z | ( 1 b q , B q ) ( 1 a p , A p ) ] ,
relations (35) and (36) transform into the form (14) and (15). Thus, we prove the theorem. Here we note that instead of solutions (14) and (15) one can uses the solutions (35) and (36).
Next, we consider several special cases of the fundamental solution of Equation (10).
Corollary 1. The Green function (14) in case where μ = 2 (second space derivative), becomes
G ( x , t ) = 1 2 | x | π n = 0 ( λ 2 ) n t ( γ + β 1 ) n n ! H 1 , 2 2 , 0 [ | x | 2 t ( α + γ 1 ) / 2 | ( 1 + ( γ + β 1 ) n , γ + α 1 2 ) ( 1 2 , 1 2 ) , ( 1 + n , 1 2 ) ] .
Corollary 2. The Green function (38) in case where γ = 1 is given by
G ( x , t ) = 1 2 | x | π n = 0 ( λ 2 ) n t β n n ! H 1 , 2 2 , 0 [ | x | 2 t α / 2 | ( 1 + β n , α 2 ) ( 1 2 , 1 2 ) , ( 1 + n , 1 2 ) ]
which, by using the properties of H-function [16], can be represented in a form given in [4], i.e.,
G ( x , t ) = 1 4 π t α n = 0 ( λ 2 ) n t β n n ! H 1 , 2 2 , 0 [ x 2 4 t α | ( 1 α 2 + β n , α ) ( 0 , 1 ) , ( 1 2 + n , 1 ) ] .
Corollary 3. The solution (40) in case of α = β = 1 reduces to the solution of the classical cable equation
G ( x , t ) = 1 4 π t n = 0 ( λ 2 ) n t n n ! H 0 , 1 1 , 0 [ x 2 4 t α | ( 0 , 1 ) ] = e x 2 4 t α 4 π t n = 0 ( λ 2 ) n t n n ! = e x 2 4 t λ 2 t 4 π t
as it should be. For λ = 0 , solution (41) turns to the Gaussian form for classical diffusion equation.
Corollary 4. The Green function (14) in case where γ = 1 and λ = 0 (space-time fractional diffusion equation), becomes
G ( x , t ) = 1 μ | x | π H 2 , 3 2 , 1 [ | x | 2 t α / μ | ( 1 , 1 μ ) , ( 1 , α μ ) ( 1 2 , 1 2 ) , ( 1 , 1 μ ) , ( 1 , 1 2 ) ] = 1 μ | x | H 3 , 3 2 , 1 [ | x | t α / μ | ( 1 , 1 μ ) , ( 1 , α μ ) , ( 1 , 1 2 ) ( 1 , 1 ) , ( 1 , 1 μ ) , ( 1 , 1 2 ) ]
which corresponds to the fundamental solution of space-time fractional diffusion equation (see for example [14,38]).
Example 1. The fundamental solution (Green function) of Equation (10), without external source term, with boundary conditions (7) and initial condition (8) is given by (14), and can be obtained if in solution (13) we set g ( x ) = δ ( x ) and f ( x , t ) = 0 .
If we set γ = 1 , we obtain the following result for the Green function given by Li and Deng [7]
G ( x , t ) = 1 μ | x | π n = 0 ( λ 2 ) n t β n n ! H 2 , 3 2 , 1 [ | x | 2 t α / μ | ( 1 , 1 μ ) , ( 1 + β n , α μ ) ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ] .
Example 2. The solution (13) in case of g ( x ) = 0 , and external source, given by f ( x , t ) = λ 2 ( V r e s t t β 1 Γ ( β ) + r m δ ( x ) t β 2 Γ ( β 1 ) ) , which corresponds to instantaneous input i e ( x , t ) = δ ( x ) δ ( t ) [7], is given by
V ( x , t ) = V r e s t λ 2 μ π Γ ( β ) n = 0 ( λ 2 ) n n ! 1 | x ξ | 0 t τ β 1 ( t τ ) ( γ + β 1 ) n + γ 1 × H 3 , 2 1 , 2 [ 2 ( t τ ) ( γ + α 1 ) / μ | x ξ | | ( 1 2 , 1 2 ) , ( n , 1 μ ) , ( 0 , 1 2 ) ( 0 , 1 μ ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 μ ) ] d τ d ξ + r m λ 2 μ π Γ ( β 1 ) n = 0 ( λ 2 ) n n ! 1 | x | 0 t τ β 2 ( t τ ) ( γ + β 1 ) n + γ 1 × H 3 , 2 1 , 2 [ 2 ( t τ ) ( γ + α 1 ) / μ | x | | ( 1 2 , 1 2 ) , ( n , 1 μ ) , ( 0 , 1 2 ) ( 0 , 1 μ ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 μ ) ] d τ .
Here we note that we use solutions (35) and (36) instead of (14) and (15) in order to apply the Euler transform formula [33]
0 t τ ρ 1 ( t τ ) σ 1 H p , q m , n [ b τ ν | ( a p , A p ) ( b q , B q ) ] d τ = Γ ( σ ) t ρ + σ 1 H p + 1 , q + 1 m , n + 1 [ b t ν | ( 1 ρ , ν ) , ( a p , A p ) ( b q , B q ) , ( 1 ρ σ , ν ) ] .
Therefore, for (44) we find
V ( x , t ) = V r e s t λ 2 μ π n = 0 ( λ 2 ) n n ! 1 | ξ | t ( γ + β 1 ) ( n + 1 ) × H 4 , 3 1 , 3 [ 2 ( t τ ) ( γ + α 1 ) / μ | ξ | | ( 1 ( γ + β 1 ) n γ , γ + α 1 μ ) , ( 1 2 , 1 2 ) , ( n , 1 μ ) , ( 0 , 1 2 ) ( 0 , 1 μ ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 μ ) , ( 1 ( γ + β 1 ) ( n + 1 ) , γ + α 1 μ ) ] d ξ + r m λ 2 μ | x | π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) ( n + 1 ) 1 × H 4 , 3 1 , 3 [ 2 t ( γ + α 1 ) / μ | x | | ( 1 ( γ + β 1 ) n γ , γ + α 1 μ ) , ( 1 2 , 1 2 ) , ( n , 1 μ ) , ( 0 , 1 2 ) ( 0 , 1 μ ) , ( 1 γ ( γ + β 1 ) n , γ + α 1 μ ) , ( 1 ( γ + β 1 ) ( n + 1 ) , γ + α 1 μ ) ] = V r e s t λ 2 μ π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) ( n + 1 ) 1 | ξ | × H 3 , 4 3 , 1 [ | ξ | 2 t ( γ + α 1 ) / μ | ( 1 , 1 μ ) , ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( ( γ + β 1 ) ( n + 1 ) , γ + α 1 μ ) ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ] d ξ + r m λ 2 μ | x | π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) ( n + 1 ) 1 × H 3 , 4 3 , 1 [ | x | 2 t ( γ + α 1 ) / μ | ( 1 , 1 μ ) , ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( ( γ + β 1 ) ( n + 1 ) , γ + α 1 μ ) ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ] ,
where we use relation (37). By using the Mellin transform of the Fox H-function
0 x ξ 1 H p , q m , n [ a x | ( a p , A p ) ( b q , B q ) ] d x = a ξ θ ( ξ ) ,
where θ ( s ) = j = 1 m Γ ( b j B j s ) j = 1 n Γ ( 1 a j + A j s ) j = m + 1 q Γ ( 1 b j + B j s ) j = n + 1 p Γ ( a j A j s ) , 0 n p , 1 m q , a i , b j C , A i , B j R + , i = 1 , 2 , , p , j = 1 , 2 , , q , finally we obtain
V ( x , t ) = 2 V r e s t λ 2 μ t γ + β 1 H 1 , 2 1 , 1 [ λ 2 t γ + β 1 | ( 0 , 1 ) ( 0 , 1 ) , ( 1 ( γ + β 1 ) , γ + β 1 ) ] + r m λ 2 μ | x | π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) ( n + 1 ) 1 × H 3 , 4 3 , 1 [ | x | 2 t ( γ + α 1 ) / μ | ( 1 , 1 μ ) , ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( ( γ + β 1 ) ( n + 1 ) , γ + α 1 μ ) ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ] .
Example 3. The solution (2.4) in case of g ( x ) = 0 , and external source given by f ( x , t ) = λ 2 ( V r e s t t β 1 Γ ( β ) + r m δ ( x ) t β 1 Γ ( β ) ) , which corresponds to instantaneous input i e ( x , t ) = δ ( x ) h ( t ) [7], where h ( t ) is the Heaviside step function which is equal to one for t 0 and zero otherwise, can be calculated in a same way as in Example 2. Thus, we find
V ( x , t ) = 2 V r e s t λ 2 μ t γ + β 1 H 1 , 2 1 , 1 [ λ 2 t γ + β 1 | ( 0 , 1 ) ( 0 , 1 ) , ( 1 ( γ + β 1 ) , γ + β 1 ) ] + r m λ 2 μ | x | π n = 0 ( λ 2 ) n n ! t ( γ + β 1 ) ( n + 1 ) × H 3 , 4 3 , 1 [ | x | 2 t ( γ + α 1 ) / μ | ( 1 , 1 μ ) , ( ( γ + β 1 ) n + γ , γ + α 1 μ ) , ( ( γ + β 1 ) ( n + 1 ) + 1 , γ + α 1 μ ) ( ( γ + β 1 ) n + γ 1 , γ + α 1 μ ) , ( 1 2 , 1 2 ) , ( 1 + n , 1 μ ) , ( 1 , 1 2 ) ] .

3. Fractional Moments

Next we calculate the fractional moments [12] of the Green function (14)
| x ( t ) | δ = 2 0 x δ G ( x , t ) d x ,       δ > 0
from where we analyze the anomalous diffusive behavior. For Levy flights the second moment does not exist, thus the fractional moment | x ( t ) | δ is introduced and then instead of second moment one finds | x ( t ) | δ 2 / δ .
Using the Mellin transform [33] for 0 < δ < μ 2 , from (47), we obtain
| x ( t ) | δ = 2 δ + 1 μ π n = 0 ( λ 2 ) n t ( γ + β 1 ) n + ( γ + α 1 ) δ μ n ! Γ ( δ μ ) Γ ( 1 2 + δ 2 ) Γ ( 1 + n + δ μ ) Γ ( 1 + ( γ + β 1 ) n + ( γ + α 1 ) δ μ ) Γ ( δ 2 ) .
From (51) we see that the second moment ( δ = 2 ) exists only if μ = 2 . We note that the even moments also exist only in case where μ = 2 .
The following special cases are worth to be mentioned. For δ 0 , recalling the asymptotic formula 1 Γ ( z ) ~ z when z < < 1 we obtain that the Green function is not normalized to one, i.e.,. it is time dependent quantity
lim δ 0 | x ( t ) | δ = E γ + β 1 ( λ 2 t γ + β 1 )
which means that the Green function is not a probability distribution function. From (52) we conclude that lim δ 0 | x | δ does not depend on parameter α . By using the asymptotic behavior of the M-L function [39,40,41] for ω > 0 , σ > 0
E σ ( ω t σ ) ~ t σ ω Γ ( 1 σ ) ,       t
it yields
lim δ 0 | x ( t ) | δ ~ t ( γ + β 1 ) λ 2 Γ ( 1 ( γ + β 1 ) ) ,       t .
Since γ + β 1 > 0 , G ( x , t ) d x shows power law decay in time.
Such non-conservation of the norm is important in certain cases, as described by the Hilfer idea of fractional generators of the dynamics [42]. The non-conservation of norm appears in the decaying of the charge density in semiconductors with exponential distribution of traps, as well as power law time decay of the ion-recombination isothermal luminescence in condensed media [43,44]. Here we note that non-conservation of norm has been observed in the analysis of fractional diffusion equation with composite time fractional derivative [13,14].
Graphical representation of | x ( t ) | 0 is given in Figure 1. From the figure one can see that for classical cable equation γ = β = 1 (solid line), | x ( t ) | 0 has an exponential decay to zero as it is expected from relation (52). For γ = β = 3 / 4 (dashed line), and γ = 3 / 4 , β = 1 / 2 (dot-dashed line) it shows slower power law decay to zero of form (55). Here we note that since 0 < γ + β 1 < 1 , | x ( t ) | 0 (52) is a completely monotone function and is always positive (see for example [45,46]).
Figure 1. Graphical representation of | x ( t ) | 0 (52) for λ = 1 .
Figure 1. Graphical representation of | x ( t ) | 0 (52) for λ = 1 .
Mathematics 03 00153 g001
For δ 2 and μ 2 , temporal behavior of the second moment x 2 ( t ) gives
lim δ 2 , μ 2 | x ( t ) | δ = x 2 ( t ) = 2 t γ + α 1 n = 0 ( 2 ) n Γ ( ( γ + β 1 ) n + γ + α ) ( λ 2 ) n t ( γ + β 1 ) n n ! = 2 t γ + α 1 E γ + β 1 , γ + α 2 ( λ 2 t γ + β 1 ) ,
which in the short time limit behaves as
x 2 ( t ) ~ 2 t γ + α 1 Γ ( γ + α )
and in the long time limit as
x 2 ( t ) ~ 2 t 1 + α 2 β γ λ 4 Γ ( 2 + α 2 β γ )
where we apply the asymptotic behavior of the three parameter M-L function [39,47,48].
E α , β δ ( ω t α ) ~ t αδ ω δ Γ ( β αδ ) ,       t .
The second moment (55) gives the rate of spreading of the Green function for μ = 2 .
Graphical representation of the second moment (55) is presented in Figure 2. Solid line corresponds to the case where γ = α = β = 1 , i.e., for the classical cable equation. Power law decay of the second moment can be observed for γ = 3 / 4 , α = β = 7 / 8 (dashed line), which is in accordance to the behavior (57). The dot-dashed line gives the power law behavior of the second moment for γ = 3 / 4 , α = 7 / 8 , β = 1 / 2 . The dotted line gives the behavior of the second moment for γ = 7 / 8 , α = 1 / 4 , β = 1 . We see that it has negative sign, which seems to have no physical meaning. In order to avoid such situations we should set 2 + α 2 β γ 0 , which can be obtained from the analysis of the Laplace transform of the second moment (55) (see relation (65), and [45,46]). However, such negativity of the second moment is interpreted in a way that for these values of the parameters the current switches directions [5]. This negativity of the second moment means that the Green function is not strictly positive for some values of parameters, i.e., it may have negative values as well.
Figure 2. Graphical representation of second moment (57) for λ = 1 .
Figure 2. Graphical representation of second moment (57) for λ = 1 .
Mathematics 03 00153 g002
For δ 4 and μ → 2, we calculate the fourth moment x 4 ( t ) which gives
lim δ 4 , μ 2 | x ( t ) | δ = x 4 ( t ) = 24 t 2 ( γ + α 1 ) n = 0 ( 3 ) n Γ ( ( γ + β 1 ) n + 2 ( γ + α 1 ) + 1 ) ( λ 2 ) n t ( γ + β 1 ) n n ! = 24 t 2 ( γ + α 1 ) E γ + β 1 , 2 ( γ + α 1 ) + 1 3 ( λ 2 t γ + β 1 ) .
The fourth moment reflects the convergence of the tails of the Green function. Its calculation is important in order to derive the non-Gaussian parameter, which is defined through the ration between the fourth moment and second moment [49,50]. The short time limit of the fourth moment yields
x 4 ( t ) ~ 24 t 2 ( γ + α 1 ) Γ ( 2 ( γ + α 1 ) + 1 )
and the long time limit
x 4 ( t ) ~ 24 t 1 + 2 α 3 β γ λ 6 Γ ( 2 + 2 α 3 β γ ) .
Graphical representation of the fourth moment (59) is given in Figure 3. Same parameters as for Figure 2 are used, i.e., γ = α = β = 1 (solid line), γ = 3 / 4 , α = β = 7 / 8 (dashed line), γ = 3 / 4 , α = 7 / 8 , β = 1 / 2 (dot-dashed line). The negativity of the forth moment (dotted line, which corresponds to γ = 7 / 8 , α = 1 / 4 , β = 1 ) appears since the inequality 3 ( γ + β 1 ) 2 ( γ + α 1 ) + 1 is not satisfied (see relation (65) from Remark 1).
Figure 3. Graphical representation of the fourth moment (59) for λ = 1 .
Figure 3. Graphical representation of the fourth moment (59) for λ = 1 .
Mathematics 03 00153 g003
Furthermore, we can calculate the even moments x 2 m ( t ) , m N , for μ = 2 , of the fundamental solution. They are given by three parameter M-L function in the following form
x 2 m ( t ) = 2 2 m π n = 0 ( λ 2 ) n t ( γ + β 1 ) n + ( γ + α 1 ) m n ! × Γ ( 1 2 + m ) Γ ( 1 + m + n ) Γ ( ( γ + β 1 ) n + ( γ + α 1 ) m + 1 ) . = ( 2 m ) ! t ( γ + α 1 ) m E γ + β 1 , m ( γ + α 1 ) + 1 m + 1 ( λ 2 t γ + β 1 ) .
From general result (62) directly follow relations (52), (55) and (59) for m = 0 , m = 1 and m = 2 , respectively.

Positivity of the Even Moments

Here we note that one can find the constraints of parameters under which the even moments (62) are non-negative. For this reason, we consider the following function
k ( t ) = t ( γ + α 1 ) m E γ + β 1 , m ( γ + α 1 ) + 1 m + 1 ( λ 2 t γ + β 1 ) .
In order (63) to be non-negative, from the Bernstein characterization theorem (see for example [45,46,51]), its Laplace transform
f ( s ) = L [ K ( t ) ] ( s ) = s ( γ + β 1 ) ( m + 1 ) ( γ + α 1 ) m 1 ( s γ + β 1 + λ 2 μ ) m + 1
should be completely monotone function (given function f : ( 0 , ) R is a completely monotone if ( 1 ) n f ( n ) ( s ) 0 for all n N { 0 } and s > 0 , see Definition 1.3 on page 2 from [45]). Since 0 < γ + β 1 1 then, from Definition 3.1 on page 15 from [45], s γ + β 1 is a Bernstein function ( f : ( 0 , ) R is a Bernstein function if f ( s ) 0 for all s > 0 , and ( 1 ) n 1 f ( n ) ( s ) 0 for all n N and s > 0 ), and thus s γ + β 1 + λ 2 is a Bernstein function too, as a sum of two Bernstein functions (see Corollary 3.7 on page 20 from [45]). From here we conclude that 1 s γ + β 1 + λ 2 is a completely monotone function (Exercise 9.9 from [52], see also [53]), and 1 ( s γ + β 1 + λ 2 ) m + 1 is a completely monotone as a product of ( m + 1 ) completely monotone functions (see Corollary 6 on page 5 from [45]). Therefore, s ( γ + β 1 ) ( m + 1 ) ( γ + α 1 ) m 1 should be completely monotone. Function of form sν, s > 0, is completely monotone for any ν ≤ 0 [45]. From here, we conclude that the following inequality
2 + α m β ( m + 1 ) γ 0
should be satisfied in order even moments (62) to be non-negative. The case with γ = 1 yields 1 + α m β ( m + 1 ) 0 . For the classical cable equation α = β = γ = 1 all the even moments are non-negative since (65) is satisfied for any m N . This can also be concluded from the Green function (41), which is non-negative function.

4. Conclusions

In this paper, a generalization of the space-time fractional cable equation is discussed. The new mathematical model takes into account the temporal memory effects and spatial non-locality. Further we derive the Green function of the generalized space-time fractional cable equation by the application of Laplace and Fourier transforms. The Green function of the proposed equation is obtained in an exact form in terms of infinite series in H-functions. We analyze the anomalous behavior of the derived model by calculating the fractional moments of the Green function. It is shown that the Green function is not normalized, i.e., it is a time dependent quantity which shows power law decay. Such non-conservation of norm has been observed in different fractional models. We analyze the even moments and we show that they may have negative sign, which means that the Green function is not always positive, and that the current switches directions. By using of Bernstein characterization theorem we show under which conditions of parameters the even moments are non-negative.

Acknowledgments

The author Zivorad Tomovski is supported under the European Commission and the Croatian Ministry of Science, Education and Sports Co-Financing Agreement No. 291823. In particular, ZT acknowledges project financing from the Maria Curie FP7-PEOPLE-2011-COFUND program NEWFELPRO Grant Agreement No. 37 – Anomalous diffusion.

Author Contributions

All authors have contributed equally.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Segev, I.; London, M. Untangling dendrites with quantitative models. Science 2000, 290, 744–750. [Google Scholar] [CrossRef] [PubMed]
  2. Baer, S.M.; Rinzel, J. Propagation of dendritic spikes mediated by excitable spines: A continuum theory. J. Neurophysiol. 1991, 65, 874–890. [Google Scholar] [PubMed]
  3. Qian, N.; Sejnowski, T.J. An electro-diffusion model for computing membrane potentials and ionic concentrations in branching dendrites, spines and axons. Biol. Cybern. 1989, 62, 1–15. [Google Scholar] [CrossRef]
  4. Henry, B.I.; Langlands, T.A.M.; Wearne, S.L. Fractional cable models for spiny neuronal dendrites. Phys. Rev. Lett. 2008, 100, 128103. [Google Scholar] [CrossRef] [PubMed]
  5. Henry, B.I.; Wearne, S.L. Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions. J. Math. Biol. 2009, 59, 761–808. [Google Scholar]
  6. Langlands, T.A.M.; Henry, B.I.; Wearne, S.L. Fractional cable equation models for anomalous electrodiffusion in nerve cells: Finite domain solutions. SIAM J. Appl. Math. 2011, 71, 1168–1203. [Google Scholar] [CrossRef]
  7. Li, C.; Deng, W. Analytical solutions, moments, and their asymptotic behaviors for the time-space fractional cable equation. Commun. Theor. Phys. 2014, 62, 54–60. [Google Scholar] [CrossRef]
  8. Podlubny, I. Fractional Differential Equations; Academin Press: San Diego, CA, USA, 1999. [Google Scholar]
  9. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, the Netherlands, 2006. [Google Scholar]
  10. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives, Theory and Applications; Gordon and Breach: London, UK, 1993. [Google Scholar]
  11. Bazhlekova, E.G.; Dimovski, I.H. Exact solution for the fractional cable equation with nonlocal boundary conditions. Cent. Eur. J. Phys. 2013, 11, 1304–1313. [Google Scholar] [CrossRef]
  12. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  13. Sandev, T.; Metzler, R.; Tomovski, Z. Fractional diffusion equation with a generalized Riemann-Liouville time fractional derivative. J. Phys. A: Math. Theor. 2011, 44, 255203. [Google Scholar] [CrossRef]
  14. Tomovski, Z.; Sandev, T.; Metzler, R.; Dubbeldam, J. Generalized space-time fractional diffusion equation with composite fractional time derivative. Phys. A 2012, 391, 2527–2542. [Google Scholar] [CrossRef]
  15. Lenzi, E.K.; Rossato, R.; Lenzi, M.K.; da Silva, L.R.; Goncalves, G. Fractional diffusion equation and external forces: solutions in a confined region. Z. Naturforschung Sect. A 2010, 65, 423–430. [Google Scholar]
  16. Huang, F.; Liu, F. The space-time fractional diffusion equation with Caputo derivatives. J. Appl. Math. Comput. 2005, 19, 179–190. [Google Scholar] [CrossRef]
  17. Huang, F.; Liu, F. The time fractional diffusion equation and advection-dispersion equation. ANZIAM J. 2005, 46, 317–330. [Google Scholar] [CrossRef]
  18. Huang, F.; Liu, F. The fundamental solution of the space-time fractional advection-dispersion equation. J. Appl. Math. Comput. 2005, 18, 339–350. [Google Scholar] [CrossRef]
  19. Liu, F.; Anh, V.; Turner, I. Numerical solution of the space fractional Fokker-Planck equation. J. Comput. Appl. Math. 2004, 166, 209–219. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, F.; Zhuang, P.; Anh, V.; Turner, I.; Burrage, K. Stability and Convergence of the difference Methods for the space-time fractional advection-diffusion equation. Appl. Math. Comput. 2007, 191, 12–20. [Google Scholar] [CrossRef]
  21. Liu, F.; Yang, Q.; Turner, I. Two new implicit numerical methods for the fractional cable equation. J. Comput. Nonlinear Dyn. 2010, 6, 01109. [Google Scholar] [CrossRef]
  22. Chen, C.; Liu, F.; Burrage, K. Numerical analysis for a variable-order nonlinear cable equation. J. Comput. Appl. Math. 2011, 236, 209–224. [Google Scholar] [CrossRef] [Green Version]
  23. Zheng, M.; Liu, F.; Turner, I.; Anh, V. A novel high order space-time spectral method for the time-fractional Fokker-Planck equation. SIAM J. Sci. Comput. 2015, 37, A701–A724. [Google Scholar] [CrossRef]
  24. Liu, F.; Chen, S.; Turner, I.; Burrage, K.; Anh, V. Numerical simulation for two-dimensional Riesz space fractional diffusion equations with a nonlinear reaction term. Cent. Eur. J. Phys. 2013, 11, 1221–1232. [Google Scholar] [CrossRef]
  25. Hilfer, R. Application of Fractional Calculus in Physics; World Scientific: Singapore, Singapore, 2000. [Google Scholar]
  26. Atangana, A.; Secer, A. A note on fractional order derivative and table of fractional derivative of some special function. Abstr. Appl. Anal. 2013, 2013, 279681. [Google Scholar] [CrossRef]
  27. Atangana, A. Drawdown in prolate spheroidal-spherical coordinates obtained via Green’s function and perturbation methods. Comun. Nonlin. Sci. Numer. Simul. 2014, 29, 1259–1269. [Google Scholar] [CrossRef]
  28. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. Chem. Phys. 2002, 284, 399–408. [Google Scholar] [CrossRef]
  29. Hilfer, R. On fractional relaxation. Fractals 2003, 11, 251–257. [Google Scholar] [CrossRef]
  30. Hilfer, R. Exact solutions for a class of fractal time random walks. Fractals 1995, 3, 211–216. [Google Scholar] [CrossRef]
  31. Hilfer, R. An extension of the dynamical foundations for the statistical equilibrium concept. Phys. A 1995, 221, 89–96. [Google Scholar] [CrossRef]
  32. Mathai, A.M.; Saxena, R.K. Distribution of a product and the structural set up of densities. Ann. Math. Statist. 1969, 40, 1439–1448. [Google Scholar] [CrossRef]
  33. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function: Theory and Applications; Springer: New York, NY, USA, 2010. [Google Scholar]
  34. Srivastava, H.M.; Gupta, K.C.; Goyal, S.P. The H-Functions of One and Two Variables with Applications; South Asian Publishers: New Delhi, Madras, India, 1982. [Google Scholar]
  35. Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  36. Sandev, T.; Tomovski, Z.; Dubbeldam, J.L.A. Generalized Langevin equation with a three parameter Mittag-Leffler noise. Phy. A 2011, 390, 3627–3636. [Google Scholar] [CrossRef]
  37. Srivastava, H.M.; Tomovski, Z. Fractional calculus with an integral operator containing a generalized Mittag-Leffler function in the kernel. Appl. Math. Comput. 2009, 211, 198–210. [Google Scholar] [CrossRef]
  38. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Solutions of fractional reaction-diffusion equations in terms of the H-function. Bull. Astr. Soc. India 2007, 35, 681–689. [Google Scholar]
  39. Saxena, R.K.; Mathai, A.M.; Haubold, H.J. Unified fractional kinetic equation and a fractional diffusion equation. Astrophys. Space Sci. Trans. 2004, 209, 299–310. [Google Scholar] [CrossRef]
  40. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelesticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  41. Seybold, H.; Hilfer, R. Numerical algorithm for calculating the generalized Mittag-Leffler function. SIAM J. Numer. Anal. 2008, 47, 69–88. [Google Scholar] [CrossRef]
  42. Hilfer, R. Classification theory for anequilibrium phase transitions. Phys. Rev. E 1993, 48, 2466. [Google Scholar] [CrossRef]
  43. Bisquert, J. Fractional diffusion in the multiple-trapping regime and revision of the equivalence with the continuous-time random walk. Phys. Rev. Let. 2003, 91, 010602. [Google Scholar] [CrossRef]
  44. Orenstein, J.; Kastner, M. Photocurrent transient spectroscopy: measurement of the density of localized states in a –As2Se3. Phys. Rev. Lett. 1981, 46, 1421. [Google Scholar] [CrossRef]
  45. Schilling, R.; Song, R.; Vondracek, Z. Bernstein Functions; De Gruyter: Berlin, Germany, 2010. [Google Scholar]
  46. Tomovski, Z.; Pogany, T.; Srivastava, H.M. Laplace type integral expressions for a certain three-parameter family of generalized Mittag-Leffler functions with applications involving complete monotonicity. J. Frankl. Inst. 2014, 351, 5437–5454. [Google Scholar] [CrossRef]
  47. Sandev, T.; Tomovski, Z. Langevin equation for a free particle driven by power law type noises. Phys. Lett. A 2014, 378, 1–9. [Google Scholar] [CrossRef]
  48. Sandev, T.; Metzler, R.; Tomovski, Z. Correlation functions for the fractional generalized Langevin equation in the presence of internal and external noise. J. Math. Phys. 2014, 55, 023301. [Google Scholar]
  49. Spanner, M.; Hoefling, F.; Schroeder-Turk, G.; Mecke, K.; Franosch, T. Anomalous transport of a tracer on percolation clusters. J. Phys.: Condens. Matter 2011, 23, 234120. [Google Scholar]
  50. Metzler, R.; Jeon, J.-H.; Cherstvy, A.G.; Barkai, E. Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and aging at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 2014, 16, 24128–24164. [Google Scholar] [CrossRef] [PubMed]
  51. Garra, R.; Gorenflo, R.; Polito, F.; Tomovski, Z. Hilfer-Prabhakar derivatives and some applications. Appl. Math. Comput. 2014, 242, 576–589. [Google Scholar] [CrossRef]
  52. Berg, C.; Forst, G. Potential Theory on Locally Compact Abelian Groups; Springer: Berlin, Germany, 1975. [Google Scholar]
  53. Fujita, Y. A generalization of the results of Pillai. Ann. Inst. Statist. Math. 1993, 45, 361–365. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Saxena, R.K.; Tomovski, Z.; Sandev, T. Analytical Solution of Generalized Space-Time Fractional Cable Equation. Mathematics 2015, 3, 153-170. https://0-doi-org.brum.beds.ac.uk/10.3390/math3020153

AMA Style

Saxena RK, Tomovski Z, Sandev T. Analytical Solution of Generalized Space-Time Fractional Cable Equation. Mathematics. 2015; 3(2):153-170. https://0-doi-org.brum.beds.ac.uk/10.3390/math3020153

Chicago/Turabian Style

Saxena, Ram K., Zivorad Tomovski, and Trifce Sandev. 2015. "Analytical Solution of Generalized Space-Time Fractional Cable Equation" Mathematics 3, no. 2: 153-170. https://0-doi-org.brum.beds.ac.uk/10.3390/math3020153

Article Metrics

Back to TopTop