Next Article in Journal
A New Perspective on Parametric Metric Spaces
Next Article in Special Issue
On the Matrix Mittag–Leffler Function: Theoretical Properties and Numerical Computation
Previous Article in Journal
The Truncation Regularization Method for Identifying the Initial Value on Non-Homogeneous Time-Fractional Diffusion-Wave Equations
Previous Article in Special Issue
Modeling Heavy Metal Sorption and Interaction in a Multispecies Biofilm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Note on the Generalized Relativistic Diffusion Equation

Dipartimento di Scienze Statistiche, “Sapienza” Università di Roma, P. le A. Moro 5, 00185 Roma, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2019, 7(11), 1009; https://doi.org/10.3390/math7111009
Submission received: 13 September 2019 / Revised: 21 October 2019 / Accepted: 22 October 2019 / Published: 24 October 2019
(This article belongs to the Special Issue Advanced Mathematical Methods: Theory and Applications)

Abstract

:
We study here a generalization of the time-fractional relativistic diffusion equation based on the application of Caputo fractional derivatives of a function with respect to another function. We find the Fourier transform of the fundamental solution and discuss the probabilistic meaning of the results obtained in relation to the time-scaled fractional relativistic stable process. We briefly consider also the application of fractional derivatives of a function with respect to another function in order to generalize fractional Riesz-Bessel equations, suggesting their stochastic meaning.

1. Introduction

In recent papers, relativistic diffusion equations have been investigated both from the physical [1] and mathematical [2,3] points of view. It is well known that these kinds of space-fractional equations are strictly related to relativistic stable processes (see, e.g., [2]).
In [2], a fractional relativistic stable process was considered, connected with a time-fractional relativistic diffusion equation involving a derivative in the sense of Caputo. In [3], the Fourier transform of the fundamental solution and the probabilistic interpretation of the solution for the Cauchy problem of a time-fractional relativistic equation was discussed. On the other hand, the fractional derivative of a function with respect to another function is a useful mathematical tool that is recently gaining more interest in relation to models involving time-varying coefficients, see, e.g., [4,5,6] and the references therein.
The aim of this short paper is to consider the generalization of the time-fractional relativistic diffusion equation by means of a Caputo fractional derivative of a function with respect to another function.
We are able to find the Fourier transform of the fundamental solution for this new class of generalized relativistic diffusion equations. Moreover, we discuss its probabilistic meaning, showing the role of the time-scaling involved in the application of time-fractional operators. A time-scaled fractional tempered stable process is considered in connection to this class of equations.
The time-scaled fractional relativistic stable process is also analyzed by evaluating the covariance and correlation coefficient and the connections with fractional-type equations.
We finally apply fractional derivatives with respect to another function in order to generalize the fractional Bessel-Riesz motion recently considered in [7].

2. Preliminaries on Fractional Relativistic Stable Processes and Fractional Operators

2.1. Fractional Relativistic Diffusion and Relativistic Stable Processes

The tempered stable subordinator (TSS) T μ , α ( t ) , t > 0 , with stability index α ( 0 , 1 ) and tempering parameter μ > 0 , is a Lévy process with Laplace transform
E ( e s T μ , α ( t ) ) = e t ( s + μ ) α μ α
and transition density
f μ , α ( x , t ) : = e μ x + μ α t h α ( x , t ) ,
where h α ( x , t ) = P r { H α ( t ) d x } and H α ( t ) is the α -stable subordinator. Moreover, let us denote with L α ( t ) : = inf { s : H α ( s ) > t } , t 0 , α ( 0 , 1 ) the inverse of the α -stable subordinator. For more details about TSS and applications we refer, for example, to [8,9].
The relativistic α -stable process is defined as a Brownian motion (hereafter denoted by B ( t ) , t 0 ) subordinated via an independent TSS, i.e.,
X μ , α ( t ) : = B ( T μ , α ( t ) ) , t > 0 .
Following the notation used, for example, in [3], in this paper we generalize the relativistic diffusion equation
u t = H α , m u ,
where the spatial differential operator
H α , m : = m m 2 α Δ α 2 , α ( 0 , 2 ) , m > 0 ,
is a relativistic diffusion operator, i.e., the infinitesimal generator of the relativistic stable process. In [3] the time-fractional generalization has been considered by replacing the first-order time-derivative with a Caputo fractional derivative. The related stochastic process has interesting properties, as discussed in [2,3].
We here briefly recall the notion of fractional tempered stable (TS) process introduced in [2].
Definition 1.
Let L β ( t ) , t 0 be the inverse of the stable subordinator, then the fractional TS process is defined as
T μ , α β ( t ) : = T μ , α ( L β ( t ) ) , t 0 , μ 0 , α , β ( 0 , 1 )
where L β is independent of the tempered stable subordinator (TSS) T μ , α .
The density of the fractional TS process T μ , α β ( t ) satisfies the fractional equation (see [2] Theorem 6)
D t β f = μ α μ + x α f , α ( 0 , 1 ) , μ 0 ,
under initial-boundary conditions
f ( x , 0 ) = δ ( x ) , f ( 0 , t ) = 0 .
We here denote by D t β the Caputo fractional derivative of order β ( 0 , 1 ) , i.e.,
D t β f ( t ) = 1 Γ ( 1 β ) 0 t ( t τ ) β f τ d τ .
Definition 2.
The fractional relativistic stable process is defined as the time-changed Brownian motion
X μ , α β ( t ) : = B ( T μ , α β ( t ) ) , t 0 .
The density of the fractional relativistic stable process X μ , α β ( t ) coincides with the solution of the time-fractional relativistic diffusion equation ([2], Theorem 16)
D t β g = μ α ( μ Δ ) α g ,
under initial and boundary conditions
g ( x , 0 ) = δ ( x ) , lim | x | g ( x ) = 0 , lim | x | x g ( x ) = 0 .
We will here present some generalizations of Equations (7) and (11) obtained by using fractional derivatives of a function with respect to another function.

2.2. Fractional Derivatives of a Function with Respect to Another Function

Fractional derivatives of a function with respect to another function have been considered in the classical monograph by Samko et al. [10] and Kilbas et al. [11] (Section 2.5). Recently they were reconsidered by Almeida in [4] where the Caputo-type regularization of the existing definition and some interesting properties are provided. We recall the basic definitions and properties for the reader’s convenience.
Let ν > 0 , f ( t ) C 1 ( [ a , t ] ) an increasing function such that f ( t ) 0 in [ a , t ] , the fractional integral of a function g ( t ) with respect to another function f ( t ) is given by
I a + ν , f g ( t ) : = 1 Γ ( ν ) a t f ( τ ) ( f ( t ) f ( τ ) ) ν 1 g ( τ ) d τ .
Observe that for f ( t ) = t β we recover the definition of Erdélyi-Kober fractional integral recently applied, for example, in connection with the Generalized Grey Brownian Motion [12].
The corresponding Caputo-type evolution operator (according to our notation) is given by
O t ν , f g ( t ) : = I a + n ν , f 1 f ( t ) d d t n g ( t ) ,
where n = [ ν ] + 1 . Here we have used the symbol O t ν , f ( · ) in order to underline the generic integro-differential nature of the time-evolution operator, beyond the fractional calculus theory.
A relevant property of Equation (14) is that, if g ( t ) = ( f ( t ) f ( a ) ) β 1 with β > 1 , then (see Lemma 1 of [4])
O t ν , f g ( t ) = Γ ( β ) Γ ( β ν ) ( f ( t ) f ( a ) ) β ν 1 .
As a consequence, the composite Mittag-Leffler function (see [6,13] for Mittag-Leffler functions)
g ( t ) = E ν ( λ ( f ( t ) f ( a ) ) ν )
is an eigenfunction of the operator O t ν , f . In the next section we will consider a = 0 .
Hereafter we denote by f ^ ( k ) the Fourier transform of a given function f ( x ) .

3. On the Generalized Relativistic Diffusion Equation

Theorem 1.
Let f ( t ) C 2 [ 0 , T ] , such that f ( 0 ) = 0 and f 0   t ( 0 , T ] , then the Fourier transform of the fundamental solution of the generalized time-fractional relativistic diffusion equation
O t ν , f u ( x , t ) = m ( m 2 / α Δ ) α / 2 u ( x , t ) , α ( 0 , 2 ) , ν ( 0 , 1 ) ,
is given by
u ^ ( k , t ) = E ν ( θ ( | k | ) ( f ( t ) ) ν ) ,
where
θ ( | k | ) = ( m α 2 + | k | 2 ) α / 2 m .
Moreover, in the case ν = 1 , the fundamental solution of the equation
1 f ( t ) t u ( x , t ) = m ( m 2 / α Δ ) α / 2 u ( x , t ) , α ( 0 , 2 )
is given by
u ^ ( k , t ) = e θ ( | k | ) f ( t ) .
Proof. 
By taking the spatial Fourier transform of Equation (17), we have that
O t ν , f u ^ ( k , t ) = { ( m α 2 + | k | 2 ) α / 2 m } u ^ ( k , t ) = θ ( | k | ) u ^ ( k , t ) ,
under the initial condition u ^ ( k , 0 ) = 1 . Then, by using the fact that
O t ν , f E ν ( θ ( | k | ) ( f ( t ) ) ν ) = θ ( | k | ) E ν ( θ ( | k | ) ( f ( t ) ) ν ) ,
we obtain the claimed result.
The case ν = 1 can be checked by simple calculations. □
Remark 1.
Observe that, for m = 0 , we obtain the Fourier transform of the fundamental solution of the generalized space-time-fractional diffusion equation
O t ν , f u ( x , t ) = ( Δ ) α / 2 u ( x , t ) , α ( 0 , 2 ) , ν ( 0 , 1 ) ,
while, for α = 2 , m = 0 , we have the Fourier transform of the generalized time-fractional heat equation
O t ν , f u ( x , t ) = Δ u ( x , t ) , ν ( 0 , 1 ) .
In the latter, some interesting cases have been considered in the literature. In particular
  • if f ( t ) = t , we recover the time-fractional diffusion equation, which is widely studied in the literature, (we refer, for example, to [14,15]).
  • if f ( t ) = ( 1 e β t ) / β , we have the fractional Dodson diffusion equation studied in [6].
  • if f ( t ) = ln t , we have the time-fractional diffusion involving a regularized Hadamard fractional derivative.
Remark 2.
In Theorem 1, we take for simplicity the condition that a = 0 in the definition of the fractional operator O t ν , f and we consider the renstriction to functions f such that f ( 0 ) = 0 . In the more general case a > 0 , f ( a ) 0 , we have that the Fourier transform of the fundamental solution for Equation (17) is given by
u ^ ( k , t ) = E ν ( θ ( | k | ) ( f ( t ) f ( a ) ) ν ) .
Much care should be given, for example, to the case of logarithmic functions, where "a" must be greater than zero.
Remark 3.
Equation (20) with α = 1 can be euristically interpreted as a relativistic Schrödinger-type equation with time-dependent mass.

4. Time-Scaled Fractional Tempered Stable Process

Let us consider the following equation
O t β , f h = μ α ( μ + x ) α h ,
under the initial-boundary conditions
h ( x , 0 ) = δ ( x ) , h ( 0 , t ) = 0 .
Let us denote by f ˜ the Laplace transform of a function f, i.e., f ˜ ( η ) : = 0 + e η x f ( x ) d x . We now evaluate the Laplace transform, with respect to the space argument, of Equation (27), which reads
O t β , f h ˜ ( η , t ) = ψ α , μ ( η ) h ˜ ( η , t ) , η > 0 ,
with initial condition h ˜ ( η , 0 ) = 1 and where ψ α , μ ( η ) : = ( η + μ ) α μ α (see [2] for details). Thus, we can write its solution as
h ˜ ( η , t ) = E β ( ψ α , μ ( η ) [ f ( t ) f ( 0 ) ] β ) .
By comparing Equation (29) with Equation (20) in [2], we can derive the following relationship holding for the corresponding process T μ , α β , f ( t ) , for any t > 0
E e η T μ , α β , f ( t ) = E e η T μ , α β ( f ( t ) ) = E β ( ψ α , μ ( η ) [ f ( t ) f ( 0 ) ] β )
and, by the unicity of the Laplace transform, we can write the following equality of the one-dimensional distributions
T μ , α β , f ( t ) = d T μ , α β ( f ( t ) ) ,
which holds for any t > 0 . Thus the process T μ , α β , f can be called time-scaled fractional tempered stable process.
Therefore, we have the following:
Theorem 2.
The density of the time-scaled fractional TS process T μ , α β ( f ( t ) ) satisfies the fractional Equation (27) under the initial-boundary Conditions (28).
From the finiteness of the moments of the standard tempered stable process and by taking into account Equation (6), we can draw the conclusion that the moment generating function of T μ , α β , f exists, for any f and α , β ( 0 , 1 ) , and is equal to Equation (30), for κ = η . By differentiating
E e κ T μ , α β , f ( t ) = E β ( ψ α , μ ( κ ) [ f ( t ) f ( 0 ) ] β ) ,
we then obtain the first and second moments of T μ , α β , f as follows
E T μ , α β , f ( t ) = α μ α 1 [ f ( t ) f ( 0 ) ] β Γ ( β + 1 )
and
E T μ , α β , f ( t ) 2 = α ( 1 α ) μ α 2 [ f ( t ) f ( 0 ) ] β Γ ( β + 1 ) + 2 α 2 μ 2 α 2 [ f ( t ) f ( 0 ) ] 2 β Γ ( 2 β + 1 ) ,
which, for β = 1 and f ( t ) = t , for any t, coincide with those of the standard tempered stable process (see [8]). Since, in this case, the process is not Lévy, we also loose the stationarity of increments.
We study the tails’ behavior of the distribution of the time-scaled fractional TS process T μ , α β , f ( t ) , i.e., we obtain an estimate for P ( T μ , α β , f ( t ) > x ) , for x + , by means of the asymptotic result holding for the usual TS subordinator T μ , α ( t ) : It is well-known that, for x + ,
P ( T μ , α ( t ) > x ) e μ x + μ α t x α Γ ( 1 + α ) sin ( π α ) t α π
(see [16]). Let f ( 0 ) = 0 , for simplicity, and let l β ( x , t ) : = P { L β ( t ) d x } be the transition density of the inverse β -stable subordinator, then we can write, by considering Equations (6) and (29) together, that
P T μ , α β , f ( t ) > x = 0 + P ( T μ , α ( z ) > x ) l β ( z , f ( t ) ) d z = [ by ( 33 ) ] e μ x x α Γ ( 1 + α ) sin ( π α ) α π 0 + e μ α z z l β ( z , f ( t ) ) d z = e μ x x α Γ ( 1 + α ) sin ( π α ) α π k 0 + e k z l β ( z , f ( t ) ) d z k = μ α = e μ x x α Γ ( 1 + α ) sin ( π α ) α π k E β ( k f ( t ) β ) k = μ α e μ x x α Γ ( 1 + α ) sin ( π α ) f ( t ) β α π β E β , β ( μ α f ( t ) β ) .
Taking the derivative out of the integral, in the previous lines, is allowed by the standard arguments for the moment generating function of the inverse stable subordinator
0 + e k z l β ( z , f ( t ) ) d z = E β ( k f ( t ) β )
(see also [17]).
Going back to the problem considered in the previous section, we are now able to provide a stochastic interpretation for the fundamental solution of the generalized fractional relativistic diffusion equation.
Theorem 3.
The fundamental solution of the generalized time-fractional relativistic diffusion Equation (17) coincides with the density of the time-scaled fractional relativistic stable process X m 2 / α , α 2 ν ( f ( t ) ) = B ( T m 2 / α , α 2 ν ( f ( t ) ) .
We omit the complete proof, since it can be easily derived from [2], Corollary 13. Indeed, we here underline that, by using a time-fractional derivative with respect to a function, we obtain a deterministic time-change; therefore the related process is the time-scaled counterpart of the original process.
We can evaluate, by conditioning the expected value and variance of the time-scaled fractional relativistic stable process X m , α / 2 ν . Indeed, we can write
E X m 2 / α , α / 2 ν ( t ) = E B ( T μ , α ν , f ( t ) ) = 0
and
v a r X m 2 / α , α / 2 ν ( t ) = E E B ( T μ , α ν , f ( t ) ) T μ , α ν , f ( t ) 2 = E T μ , α ν , f ( t ) = α μ α 1 [ f ( t ) f ( 0 ) ] ν Γ ( ν + 1 ) .
This confirms that the mean-square displacement of the process behaves as f ( t ) ν and thus crucially depends on the choice of f .
We also evaluate the covariance and correlation coefficient of the process B ( T μ , α ν , f ( t ) ) , at least in the special case where f is a sublinear or linear function. Let moreover f ( 0 ) = 0 , for simplicity. We prove, under these assumptions, that the process displays a long-range dependence (LRD). We need the following definition for LRD of non-stationary processes (see [9,18]):
Definition 3.
Let s > 0 and t > s . If, for a process Z ( t ) , t 0 , the following asymptotic behavior of the correlation function holds
c o r r ( Z ( s ) , Z ( t ) ) c ( s ) t d , t + ,
where c ( s ) is a constant (for fixed s) and d ( 0 , 1 ) , then Z ( t ) is said to have the LRD property.
We restrict our analysis to the class of increasing, positive functions f , such that lim t + f ( t ) t ρ = c o n s t , or equivalently f ( t ) t ρ , for ρ ( 0 , 1 ] . It corresponds to considering sublinear functions, for which | f ( t s ) | f ( t ) f ( s ) , for any t s , (recalling that the last difference is always positive for an increasing function). The equality holds in the special case of linear functions. We start by recalling that the Brownian motion has stationary increments and the same holds for the tempered stable and the inverse stable subordinators (see, for example, [16,19], respectively). Moreover, it is easy to check that the stationarity of increments is preserved under composition of processes. As a consequence, taking into account Equations (6) and (10), we can write, for s t , that
c o v ( B ( T μ , α ν , f ( s ) ) , B ( T μ , α ν , f ( t ) ) ) = 1 2 E B ( T μ , α ν , f ( s ) ) 2 + E B ( T μ , α ν , f ( t ) ) 2 E B ( T μ , α ν , f ( t ) ) B ( T μ , α ν , f ( s ) ) 2 = 1 2 E B ( T μ , α ν , f ( s ) ) 2 + 1 2 E B ( T μ , α ν , f ( t ) ) 2 1 2 E B ( T μ , α ν , f ( t s ) ) 2 = 1 2 α μ α 1 f ( t ) ν Γ ( ν + 1 ) + 1 2 α μ α 1 f ( s ) ν Γ ( ν + 1 ) α μ α 1 f ( t s ) ν 2 Γ ( ν + 1 ) 1 2 α μ α 1 f ( t ) ν Γ ( ν + 1 ) j = 1 ν j f ( t ) j ( f ( s ) ) j + f ( s ) ν f ( t ) ν = 1 2 α μ α 1 f ( t ) ν Γ ( ν + 1 ) ν f ( s ) f ( t ) l = 0 ν 1 l f ( t ) l ( f ( s ) ) l l + 1 + f ( s ) ν f ( t ) ν = 1 2 α μ α 1 f ( t ) ν Γ ( ν + 1 ) ν f ( s ) f ( t ) + f ( s ) ν f ( t ) ν + O ( f ( t ) 2 ) 1 2 α ν μ α 1 f ( s ) Γ ( ν + 1 ) f ( t ) ν 1 .
Thus we easily check that the correlation coefficient is bounded below by the following asymptotic expression f ( s ) 1 ν / 2 f ( t ) ν / 2 1 f ( s ) 1 ν / 2 t ρ ν / 2 ρ (for any fixed s and t + ). This means that B ( T μ , α ν , f ( t ) ) exhibits LRD behavior with d = ρ ( 1 ν / 2 ) ( ρ / 2 , 1 ) , at least inside the class of sublinear functions (asymptotically behaving as fractional powers of t).

5. An Application to Generalized Fractional Bessel-Riesz Motion

In the recent paper [7], the authors give a stochastic representation for the solution of the following fractional Cauchy problem
D t β p = λ ( Δ ) α / 2 ( I Δ ) γ / 2 p , , γ 0 , α ( 0 , 2 ] , β ( 0 , 1 ]
p ( x , 0 ) = δ ( x ) ,
where ( Δ ) α / 2 and ( I Δ ) γ / 2 are the inverses of the Riesz and the Bessel potential respectively. The Fourier transform of the Cauchy Equation (34), has been obtained in the form
G ^ ( k , t ) = E β λ t β k α ( 1 + k 2 ) γ / 2 .
In [7], it was proved that the fundamental solution for Equation (34) coincides with the density of the stochastic process
X B R ( t ) = W ( L α , γ ( L β ( t ) ) ) , t 0
where W is the n-dimensional Brownian motion, L β ( t ) is the inverse stable subordinator (according to our notation) and L α , γ is the Lévy subordinator with Laplace exponent
Φ ( s ) = λ s α / 2 ( 1 + s ) γ / 2 , s > 0 .
All these three processes are jointly independent.
Let us consider the generalized fractional counterpart of the Cauchy Equation (34) involving fractional derivatives with respect to another function, i.e.,
O t ν , f u ( x , t ) = λ ( Δ ) α / 2 ( I Δ ) γ / 2 u ( x , t ) .
Assuming that f ( t ) is a suitable smooth function such that f ( 0 ) = 0 , we have that the fundamental solution of Equation (39) given by
u ( x , t ) = E β λ ( f ( t ) ) β k α ( 1 + k 2 ) γ / 2
and therefore it coincides with the density of the time-scaled stochastic processes depending on f
X B R ( t ) = W ( L α , γ ( L β ( f ( t ) ) ) ) , t 0

6. Conclusions and Open Problems

We have analyzed here some applications of the time-fractional derivative with respect to a function, in the context of the so-called fractional relativistic diffusion equation. It is well-known that the most interesting case from the physical point of view corresponds to α = 1 , where the operator H 1 , m appearing in Equation (5) represents the free energy of the relativistic Schrödinger equation. A physical discussion about the utility and meaning of the time-fractional generalizations of the relativistic diffusion is still missing. In view of the wide exhisting literature on the time-fractional Schrödinger equation, it is worth investigating this topic for future studies, also for the stochastic interpretation discussed here.
From the mathematical point of view, another interesting topic is the connection between higher order equations and time-fractional equations (see, e.g., [20]). In this case the relation is not trivial and it probably works only for specific values of the real order of the derivative.
Finally, a probabilistic interpretation of the fractional integrals with respect to a function can be developed, following the discussion presented in [21].

Author Contributions

The authors contributed equally to this work.

Acknowledgments

The work of R.G. was carried out in the framework of the activities of the GNFM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baeumer, B.; Meerschaert, M.M.; Naber, M. Stochastic models for relativistic diffusion. Phys. Rev. E 2010, 82, 011132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Beghin, L. On fractional tempered stable processes and their governing differential equations. J. Comput. Phys. 2015, 293, 29–39. [Google Scholar] [CrossRef]
  3. Shieh, N.R. On time-fractional relativistic diffusion equations. J. Pseudo-Differ. Oper. Appl. 2012, 3, 229–237. [Google Scholar] [CrossRef]
  4. Almeida, R. A Caputo fractional derivative of a function with respect to another function. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 460–481. [Google Scholar] [CrossRef]
  5. Colombaro, I.; Garra, R.; Giusti, A.; Mainardi, F. Scott-Blair models with time-varying viscosity. Appl. Math. Lett. 2018, 86, 57–63. [Google Scholar] [CrossRef] [Green Version]
  6. Garra, R.; Giusti, A.; Mainardi, F. The fractional Dodson diffusion equation: A new approach. Ricerche di Matematica 2018, 1–11. [Google Scholar] [CrossRef]
  7. Anh, V.V.; Leonenko, N.N.; Sikorskii, A. Stochastic representation of fractional Bessel-Riesz motion. Chaos Solitons Fractals 2017, 102, 135–139. [Google Scholar] [CrossRef] [Green Version]
  8. Gajda, J.; Kumar, A.; Wylomańska, A. Stable Lévy process delayed by tempered stable subordinator. Stat. Probab. Lett. 2019, 145, 284–292. [Google Scholar] [CrossRef]
  9. Kumar, A.; Gajda, J.; Wylomanska, A.; Poloczanski, R. Fractional Brownian motion delayed by tempered and inverse tempered stable subordinators. Methodol. Comput. Appl. Probab. 2019, 21, 185–202. [Google Scholar] [CrossRef]
  10. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science Publishers: Yverdon-les-Bains, Switzerland, 1993. [Google Scholar]
  11. Kilbas, A.A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science Limited: Amsterdam, The Netherlands, 2006. [Google Scholar]
  12. Pagnini, G. Erdélyi-Kober fractional diffusion. Fract. Calc. Appl. Anal. 2012, 15, 117–127. [Google Scholar] [CrossRef]
  13. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014. [Google Scholar]
  14. Mainardi, F. The fundamental solutions for the fractional diffusion-wave equation. Appl. Math. Lett. 1996, 9, 23–28. [Google Scholar] [CrossRef] [Green Version]
  15. Orsingher, E.; Beghin, L. Fractional diffusion equations and processes with randomly varying time. Ann. Probab. 2009, 37, 206–249. [Google Scholar] [CrossRef]
  16. Gupta, N.; Kumar, A.; Leonenko, N. Mixtures of tempered stable subordinators. arXiv 2019, arXiv:1905.00192v1. [Google Scholar]
  17. Beghin, L.; Macci, C.; Martinucci, B. Random time-changes and asymptotic results for a class of continuous-time Markov chains on integers with alternating rates. arXiv 2018, arXiv:1802.06434v1. [Google Scholar]
  18. Ayache, A.; Cohen, S.; Lévy Véhel, J. The covariance structure of multifractional Brownian motion, with application to long range dependence. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing—ICASSP 2000, Istanbul, Turkey, 5–9 June 2000. [Google Scholar]
  19. Kaj, I.; Martin-Lof, A. Scaling Limit Results for the Sum of Many Inverse Levy Subordinators. arXiv 2012, arXiv:1203.6831. [Google Scholar]
  20. D’Ovidio, M. On the fractional counterpart of the higher-order equations. Stat. Probab. Lett. 2011, 81, 1929–1939. [Google Scholar] [CrossRef] [Green Version]
  21. Tarasov, V.E.; Tarasova, S.S. Probabilistic Interpretation of Kober Fractional Integral of Non-Integer Order. Prog. Fract. Differ. Appl. 2019, 5, 1–5. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Beghin, L.; Garra, R. A Note on the Generalized Relativistic Diffusion Equation. Mathematics 2019, 7, 1009. https://doi.org/10.3390/math7111009

AMA Style

Beghin L, Garra R. A Note on the Generalized Relativistic Diffusion Equation. Mathematics. 2019; 7(11):1009. https://doi.org/10.3390/math7111009

Chicago/Turabian Style

Beghin, Luisa, and Roberto Garra. 2019. "A Note on the Generalized Relativistic Diffusion Equation" Mathematics 7, no. 11: 1009. https://doi.org/10.3390/math7111009

APA Style

Beghin, L., & Garra, R. (2019). A Note on the Generalized Relativistic Diffusion Equation. Mathematics, 7(11), 1009. https://doi.org/10.3390/math7111009

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop