1. Introduction
We consider space–time fractional equations with general fractional operators in space and time. More precisely, we deal with a very general fractional space operator that covers a large class of non-local operators, such as fractional Laplacians.
This non-local space-operator can be related to time-changed processes where the time change is given by a subordinator (for definitions, examples and applications, see ref. [
1]) characterized by a symbol, which is a Bernstein function (see [
2]).
Additionally, the non-local time operator is very general, and it includes a huge class of convolution-type operators, such as the Caputo fractional derivatives. This non-local time-operator can be related to time-changed processes where the time change is given by an inverse to a subordinator characterized by a symbol, which is again a Bernstein function.
The literature on space–time fractional equations and their applications is very extensive. We mention here only some basic works [
3,
4,
5,
6,
7] and the references therein. Connections with Sturm–Liouville problems were investigated in [
8], whereas for the the higher-order counterpart, we refer to [
9,
10]. For the fractional Cauchy problem on manifolds, we reference the work in [
11]. Recently, the authors also obtained results on irregular domains in the case of randomly varying fractals [
12,
13].
The aim of this paper is to relate the asymptotic analysis of space–time fractional equations to the convergence of corresponding symbols (see Theorems 5 and 6). Our result extends Theorem 7 in [
9], where asymptotic properties for time-changed processes were investigated for pseudo-processes. The symbol of a subordinator may be approximated by the symbols of a continuous-time random walk (see [
9] (Theorem 5)). The base process is Markovian, but it is driven by a signed measure; that is, the governing equation of the base process is a higher-order equation.
We highlight that the results of the present paper provide a useful tool for studying the approximation of space–time fractional equations in several contexts since the theory of Dirichlet forms allows us to describe many structures in an appropriate functional environment.
For example, we can approximate space–time fractional equations related to relativistic
-stable processes, spherical symmetric
-stable processes, and gamma processes with suitable sequences of relativistic
-stable processes in
Moreover, we can use the results of the present paper for “denoising” variance gamma processes, that is, Brownian motions time changed using a gamma subordinator as a random time. Such a kind of “denoising” can be carried out by considering the asymptotic limits of the parameters characterizing the symbol of the subordinator (see
Section 7).
The paper is set out as follows. In
Section 2, we recall some basic facts about processes associated with Dirichlet forms. In
Section 3, we introduce symbols corresponding to Bernstein functions associated with subordinate processes. In
Section 4, we consider space fractional equations, and we recall asymptotic results via the convergence of symbols obtained in [
14]. In
Section 5, we introduce the time fractional equations associated with inverse processes. In
Section 6, we consider space–time fractional equations, and we prove asymptotic results via the convergence of symbols. Finally, in the last section, we provide some examples and applications.
2. Processes Associated with Dirichlet Forms
We now recall some basic facts about processes associated with Dirichlet forms (see [
15]).
We consider an m-symmetric right process
X on a Lusin space
Without loss of generality,
X can be assumed as an
m-symmetric Hunt process associated with a regular Dirichlet form
on a locally compact separable metric space
where
m is a Radon measure with full support on
E (by using quasi-homeomorphism, see [
16]). The
-infinitesimal generator
is a non-positive definite self-adjoint operator, and it has the following spectral representation:
with domain
Here,
is the spectral family of
it is a right continuous increasing sequence of orthogonal projections in
with
, and
the identity operator. The corresponding Dirichlet form
associated with
X is defined as follows:
for
, where
We highlight that the Dirichlet form
can be described by using spectral representation in the following way:
for
where
We now recall the definition of Mosco convergence (see [
17]). We consider a sequence of forms
with domain
and a form
with domain
The forms
,
can be defined in the whole of
by setting the following:
Definition 1. A sequence of forms M-converges to a form in if
For every converging weakly to u in For every , there exists converging strongly to u in , such that
The M-convergence of forms can be characterized in terms of convergence of the resolvent operators and semigroup operators.
Theorem 1. (see [
18]).
M-converges to in if and only if the sequence of the resolvent operators converges to the resolvent operator in the strong operator topology of . Theorem 2. (see [
18]).
M-converges to in if and only if for every the sequence of the semigroup operators converges to the semigroup operator associated with the strong operator topology of uniformly on every interval . 3. Symbols and Associated Subordinators
Here, we focus on subordination, which is a time change given by a subordinator. The corresponding semigroup is termed a subordinated semigroup.
We consider the following symbols corresponding to Bernstein functions
where
is a Lévy measure on
with
. We also recall that
and
is the so called tail of the Lévy measure (see [
2]). We highlight that the symbol
can be related to the Laplace exponent of a subordinator
H—that is, a one-dimensional almost surely increasing Lévy process—as follows:
(see [
1]).
Typical examples are the following:
, associated with stable subordinator;
with and associated with generalized stable subordinator;
with associated with inverse Gaussian subordinator;
with associated with gamma subordinator.
By using spectral representation, we have
For example, if H is a stable subordinator with symbol and X is a Brownian motion, then is the fractional Laplacian.
For the process
X with generator
and the independent subordinator
we define the time-changed process as follows:
for
.
The process
,
can be considered in order to solve the following equation:
as the probabilistic representation of the solutions to (
3) is given by
where
is the lifetime of
, which is the part process of
on
E.
4. Space Fractional Equations via Convergence of Symbols
We consider a subordinator H with Laplace exponent and subordinators with Laplace exponent
We suppose that the process
X is independent of
H and
, and we define the subordinate processes
and
We denote by the corresponding Dirichlet form associated with and by the Dirichlet form associated with
By spectral representation, we have
and
for
In a similar way,
and
for
Thus, the generator of
is
, where
with domain
and the generator of
is
, where
with domain
We recall the following results.
Lemma 1. (Lemma 4.2 of [
14]).
Assume thatThen, the Dirichlet form M-converges to
For an open subset D of we denote by the part process of on D and by the part process of on is the corresponding Dirichlet form associated with the part process , and is the corresponding Dirichlet form associated with part process
Theorem 3. (Theorem 4.3 of [
14]).
Assume thatThen, the Dirichlet form M-converges to
5. Inverse of Subordinators and Time Fractional Derivatives
We introduce the inverse process
L
and define (for
) the time-changed process
This process is strictly related to the following time fractional equation:
Here, the fractional time operator
is defined in the following way. For
and
we consider the set
of (piecewise) continuous function on
of exponential order
w, such that
. We define the operator
, such that
where
is the Laplace transform of
u.
By using
(the tail of the Lévy measure), we can also write the following:
We highlight that operator
was previously considered in [
19] (Remark 4.8) as the generalized Caputo derivative. In particular, we observe the following:
if the operator becomes the ordinary derivative;
if
, the operator
becomes the Caputo fractional derivative
with
if the operator becomes tempered fractional derivative;
if for , the operator becomes the telegraph fractional operator.
The probabilistic representation of the solution of time fractional Equation (
4) is given by the following:
where
is the lifetime of
, the part process of
on
E. In particular, the following theorem states the existence and uniqueness of a strong solution in
to (
4) (see [
12,
20,
21]).
Theorem 4. ([
12] (Theorem 5.2)).
The function (5) is the unique strong solution in to (4) in the sense that: - (1)
is such that and ;
- (2)
is such that ;
- (3)
, holds m-a.e in E;
- (4)
, as .
6. Time–Space Fractional Equations via Convergence of Symbols
As in
Section 4, we consider symbols
and
and their corresponding subordinators,
H and
We assume that the process
X is independent of
H and
and consider the subordinate processes
and
Moreover, we consider a symbol
and the corresponding inverse of its associated subordinator denoted again by
We examine the following time–space fractional equations
The probabilistic representation of the solution can be written in terms of the time-changed process
, that is,
Similarly, the probabilistic representation of the solution to
can be written in terms of the time-changed process
that is,
is the set of continuous functions from to , which are right continuous on with left limits on , where ∂ is the cemetery point. In the following Theorem 5, we prove the asymptotic results for space–time fractional equations via the convergence of symbols.
Proof. As
by Lemma 1, we observed that
M-converges to
By using the results of a recent paper [
12], we found the convergence of the time-changed processes from the M-convergence of the forms
More precisely, from the
M-convergence of the forms,
we found the strong convergence of the corresponding semigroups. Then, by using Theorem 17.25 in [
22] by means of which we know that the strong convergence of semigroups (Feller semigroups) is equivalent to the weak convergence of measures if
in distribution, we found that
in distribution in
.
By using results in [
23], we found that as
in distribution in
. □
As in
Section 4,
D is an open subset of
We use
to denote the part process of
on
D and
to denote the part process of
on
D.
Proof. By Theorem 3, following the same tools of previous proof, we obtained the result. □
Remark 1. We remark that a probabilistic interpretation in terms of the mean lifetime of the base and time-changed processes was given recently in [24]. 7. Examples and Applications
In this section, we present some examples to illustrate the main results of this paper. First, we consider the case where
X is a Brownian motion in
running twice as fast as the standard Brownian motion. We consider the following relativistic
-stable process in
with
where
is a sequence in
When
the time changed process
is a relativistic
-stable process in
If converges to some in as we find that tends to a relativistic -stable process in
If converges to some with as we find that tends to a spherical symmetric -stable process in .
If
converges to some
with
as
we find that
tends to
which is related to a gamma process.
Then, by using the results of the previous sections, we can approximate space–time fractional equations related to relativistic -stable processes, spherical symmetric -stable processes, and gamma processes with suitable sequences of relativistic -stable processes in
We highlight that other examples can be given in a similar way by replacing X with another kind of symmetric process, such as a spherically symmetric -stable process, symmetric Lévy process, and symmetric diffusions with infinitesimal generators of divergence form.
Another interesting example is the following. Consider the sequences of symbols
that are related to a gamma processes with parameters
and
The parameters
can be related to the plot of the observed data that may fit the path of a realization of the Laplace motion (also termed variance gamma processes, that is, the Brownian motion time-changed by a gamma subordinator).
If the parameters
characterizing the phenomenon satisfy
as
we find that
tends to
as
; that is we have a sort of “denoising”, and the underline (base) process appears.