Next Article in Journal / Special Issue
Green’s Functions for a Fractional Boundary Value Problem with Three Terms
Previous Article in Journal
Telling the Wave Function: An Electrical Analogy
Previous Article in Special Issue
Detailed Error Analysis for a Fractional Adams Method on Caputo–Hadamard Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Caputo Fractional Evolution Equations in Discrete Sequences Spaces

by
Alejandro Mahillo
and
Pedro J. Miana
*,†
Calle Pedro Cerbuna 12, Departamento de Matemáticas, IUMA, Universidad de Zaragoza, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 30 August 2022 / Revised: 20 September 2022 / Accepted: 26 September 2022 / Published: 11 October 2022
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inclusions)

Abstract

:
In this paper, we treat some fractional differential equations on the sequence Lebesgue spaces p ( N 0 ) with p 1 . The Caputo fractional calculus extends the usual derivation. The operator, associated to the Cauchy problem, is defined by a convolution with a sequence of compact support and belongs to the Banach algebra 1 ( Z ) . We treat in detail some of these compact support sequences. We use techniques from Banach algebras and a Functional Analysis to explicity check the solution of the problem.

1. Introduction

The main objective of this paper is to study the following semidiscrete Cauchy differential equation
t u ( n , t ) = B u ( n , t ) + g ( n , t ) , n N 0 , t > 0 , u ( n , 0 ) = φ ( n ) , n N 0 ,
where B is a convolution operator in the discrete variable, i.e.,
B u ( n , t ) = j 0 b ( n j ) u ( j , t ) ,
and the sequence b belongs to the Banach algebra 1 ( Z ) . A first example is the one-dimensional discrete Laplacian, Δ d , which is defined by b = δ 1 2 δ 0 + δ 1 , where δ i denotes the discrete Dirac measure given by the Kronecker delta, i.e., δ i ( n ) = 1 if i = n and 0 in other case. Equation (1) is usually called the lattice diffusion equation or the semidiscrete heat equation.
These classes of equations have received a wide interest in the mathematical literature in the last years. They appear in diverse areas of knowledge. For example, in probability theory, the function u ( n , t ) of (1) with B = Δ d , expresses the probability that a continuous-time symmetric random walk arrives at point n at time t ; ([1], [Section 4]). In chemical physics, (1) describes the flow of a liquid in an infinite row of tanks where two neighbors are always connected [2], [Section 3]. Another amazing application takes place in transport theory. Equation (1) expresses the dynamics of an infinite chain of cars, each of them being coupled to its two neighbors. The function u ( n ; t ) is the displacement of car n at time t from its equilibrium point ([3], [Example 1]). Quite recently, Slavik [4] studied the asymptotic behavior of solutions of (1) when B = Δ d , showing that a bounded solution approaches the average of the initial values if the average exists. In the case that b = δ 1 δ 0 in (2), we obtain the forward difference operator B = Δ and then the Equation (2) describes the semidiscrete transport system, treated recently by Abadias et al. in [5].
Other interesting references, such as [6,7], present fundamental solutions for (1) and the second-order semidiscrete equation
t t u ( n , t ) = B u ( n , t ) + g ( n , t ) , n N 0 , t > 0 , u ( n , 0 ) = φ ( n ) , u t ( n , 0 ) = ϕ ( n ) , n N 0 ,
when B = ( Δ d ) α is the fractional power of discrete fractional Laplacian. In the particular case of [7], the authors apply operator theory techniques and some of the properties of the Bessel functions to obtain a theory of uni-parametric operators ( C 0 -semigroups and cosine operators) generated by Δ d and ( Δ d ) α on the Lebesgue space p ( Z ) . Moreover, note that the fractional forward difference operator B = ( Δ ) α has been treated in [5] where the maximum and comparison principles in the context of Fourier Analysis are shown.
However, there is no attempt (to the best of our knowledge) to present explicitly fundamental solutions of the general Equation (1) on the sequence space p ( N 0 ) for p 1 instead of on the sequence space p ( Z ) for p 1 ([8]).
The main technique in this paper is that we apply our knowledge from Banach algebras and Functional Analysis to fractional differential systems. This useful approach that we follow in this paper, allows us to obtain a completely new point of view. We prove results by introducing this new method and describing both the qualitative and quantitative behavior of the fundamental solutions of (1) in a unified way.
More generally, and to present simultaneously our studies of the subdiffusive and superdiffusive cases connected to Equations (1) and (3), in this article, we deal with the representation of the fundamental solutions for the following semidiscrete system:
D t β u ( n , t ) = B u ( n , t ) + g ( n , t ) , n N 0 , t > 0 , u ( n , 0 ) = φ ( n ) , n N 0 ,
in case 0 < β 1 and
D t β u ( n , t ) = B u ( n , t ) + g ( n , t ) , n N 0 , t > 0 , u ( n , 0 ) = φ ( n ) , u t ( n , 0 ) = ϕ ( n ) , n N 0 ,
in case 1 < β 2 . In both cases, B is the convolution operator B f ( n ) : = ( b f ) ( n ) defined for f p ( N 0 ) , p [ 1 , ] , b 1 ( Z ) and β ( 0 , 2 ] is a real number. The symbol D t β denotes the Caputo fractional derivative of order β > 0 .
The paper is organized as follows. In the first section, we introduce the main results about the Banach algebras and, in particular, about the spaces p ( N 0 ) and 1 ( Z ) . In the second section, we consider some particular finite difference operators in p ( N 0 ) , mainly
D f ( n ) = f ( n + 1 ) f ( n 1 ) = ( ( δ 1 δ 1 ) f ) ( n ) , Δ 2 f ( n ) = f ( n ) 2 f ( n + 1 ) + f ( n + 2 ) = ( ( δ 0 2 δ 1 + δ 2 ) f ) ( n ) , 2 f ( n ) = f ( n ) 2 f ( n 1 ) + f ( n 2 ) = ( ( δ 0 2 δ 1 + δ 2 ) f ) ( n ) ,
for f p ( N 0 ) . Finally, we present Theorem 6 where we include the representation of the fundamental solutions for semidiscrete Caputo fractional differential equations.
This paper contains part of the results included in the Master Thesis of the first author, entitled “Semigrupos y operadores coseno generados por operadores de diferencias finitas en espacios de sucesiones de Lebesgue”, Universidad de La Rioja, (2021).
Notation The usual set numbers N , N 0 = N { 0 } , Z , R and C are used. We write as T = { e i θ : θ [ π , π ) } the unit circumference (or also called torus) and D = { z C : | z | 1 } . The Dirac measures δ 0 and δ n are δ n ( j ) = 0 if n j and δ n ( n ) = 1 for n , j Z . We denote by χ I the indicator function on set I (i.e., χ I ( n ) = 1 if n I and χ I ( n ) = 0 if n I ).
Furthermore, I n and J n are the Bessel functions and H n the usual Hermite polynomials; Γ is the Gamma function, Φ β the Wright function and E α , β the Mittag-Leffler function.
Given X Banach space, X is the dual of Banach space and B ( X ) the set of linear bounded operators on X; given A B ( X ) , we write A B ( X ) is the adjoint of the operator A.

2. A Banach Algebra Framework

Given 1 p , we recall that the Banach spaces ( p ( N 0 ) , p ) are formed by infinite sequences f = ( f ( n ) ) n N 0 C such that
| | f | | p : = n = 0 | f ( n ) | p 1 p < , 1 p < ; | | f | | : = sup n 0 | f ( n ) | < .
We remind that the natural embeddings 1 ( N 0 ) p ( N 0 ) ( N 0 ) , for 1 p , and that the dual of p ( N 0 ) is identified with p ( N 0 ) where 1 p + 1 p = 1 for 1 < p < and p = 1 if p = .
In the case that f 1 ( N 0 ) and g p ( N 0 ), we define
( f g ) ( n ) : = j = 0 n f ( n j ) g ( j ) , n Z .
From Young’s Inequality, it follows that f g p ( N 0 ) , and f g p f 1 g p . We denote by T g B ( p ( N 0 ) ) defined by T g ( f ) = f g for g p ( N 0 ) . The element g is called the symbol of the convolution operator T g .
Note that ( 1 ( N 0 ) , ) is a commutative Banach algebra with identity that we denote by δ 0 : = χ { 0 } . We observe that δ 1 δ 1 = δ 2 and, in general, δ n δ m = δ n + m for n , m N 0 . As usual, we write f 2 = f f and f n = f n 1 f for n 2 .
The Gelfand transform associated to ( 1 ( N 0 ) , ) , is the Z -transform, Z : 1 ( N 0 ) C ( D ) , (or Taylor series), where
Z ( f ) ( z ) : = n 0 f ( n ) z n , z D .
We recall that the resolvent set of f, denoted as
ρ 1 ( N 0 ) ( f ) : = { λ C : ( λ δ 0 f ) 1 1 ( N 0 ) }
and the spectrum of f, σ 1 ( N 0 ) ( f ) = C \ ρ 1 ( N 0 ) ( f ) .
In what follows, we apply the general theory of Functional Analysis and commutative Banach Algebra as framework. In the following theorem, we collect some results that will be of our interest, see [9].
Theorem 1.
The following properties hold:
 (i) 
The spectrum Spec ( 1 ( N 0 ) ) is compact and, consequently, homeomorphic to the unit complex circle, D : = { z C : | z | 1 } .
 (ii) 
σ 1 ( N 0 ) ( f ) { z C ; | z | < f 1 } and
( λ δ 0 f ) 1 = n 0 λ n 1 f n , f 1 < | λ | .
 (iii) 
The algebra 1 ( N 0 ) is a semi-simple regular Banach algebra and the Z -transform is injective.
 (iv) 
Z ( f g ) = Z ( f ) Z ( g ) and
σ 1 ( N 0 ) ( f ) = Z ( f ) ( D ) , f 1 ( N 0 ) .
 (v) 
Given b 1 ( N 0 ) and the linear convolution operator T b ( f ) : = b f for f p ( N 0 ) for p 1 . Then,
σ B ( p ( N 0 ) ) ( T b ) σ 1 ( N 0 ) ( f ) = Z ( f ) ( D ) .
We recall that the Banach algebra ( 1 ( Z ) , 1 ) is formed by bi-infinite sequences f = ( f ( n ) ) n Z C such that
f 1 : = n = | f ( n ) | < .
Given f , g 1 ( Z ) , the product in the algebra 1 ( Z ) is the usual convolution product given by
( f g ) ( n ) : = j = f ( n j ) g ( j ) , n Z .
Note that ( 1 ( Z ) , ) is also a commutative Banach algebra with identity and δ n δ m = δ n + m for n , m Z . The Gelfand transform associated to ( 1 ( Z ) , * ) , is the Fourier series (or discrete Fourier transform), F : 1 ( Z ) C ( T ) , where
f ^ ( θ ) : = F ( f ) ( e i θ ) : = n Z f ( n ) e i n θ , θ T .
The spectrum Spec ( 1 ( Z ) ) is homeomorphic to the unit torus, T : = { z C : | z | = 1 } and σ 1 ( Z ) ( b ) = F ( b ) ( T ) ([8], (Theorem 2.1)).
Definition 1.
Given α , β > 0 , we define the vector-valued Mittag-Leffler function, E α , β : 1 ( Z ) 1 ( Z ) , by
E α , β ( a ) : = j = 0 a j Γ ( α j + β ) , a 1 ( Z ) .
Note that
E 1 , 1 ( a ) = j = 0 a j j ! = e a ; E 2 , 1 ( a ) = j = 0 a j ( 2 j ) ! .
The set exp ( 1 ( Z ) ) : = { e a ; a 1 ( Z ) } is the connected component of δ 0 in the set of regular elements in 1 ( Z ) ([9], Theorem 6.4.1).
Now, we remind the usual terminology in semigroup theory: an element a 1 ( Z ) is called the generator of the entire group given by the exponential function ( e z a ) z C . The cosine function is expressed by its generator in terms of Mittag-Leffer function C o s ( z , a ) : = E 2 , 1 ( z 2 a ) , see [10], (Sections 3.1 and 3.14). Moreover, the Laplace transform of an entire group or a cosine function is connected with the resolvent of its generator as follows:
( λ a ) 1 = 0 e λ s e a s d s , λ > a 1 , λ ( λ 2 a ) 1 = 0 e λ s C o s ( s , a ) d s , λ > a 1 ,
see, for example, ([10], p. 213).
Example 1.
For α , β > 0 , we have that
E α , β ( z δ 0 ) = E α , β ( z ) δ 0 ; E α , β ( z δ 1 ) = j = 0 z j δ j Γ ( α j + β ) .
In particular, e z δ 1 = j = 0 z j δ j j ! and C o s ( z , δ 1 ) = j = 0 z 2 j δ j ( 2 j ) ! are generated by δ 1 .
In the next proposition, we present some technical properties of these Mittag-Leffler functions in the Banach algebra 1 ( Z ) . As usual, we consider vector-valued integration (in the sense of Bochner) in the Banach space 1 ( Z ) , see, for example, ([11], Section 1.2).
Proposition 1.
For α , β > 0 and a 1 ( Z ) , we have that
 (i) 
E α , β ( a ) 1 E α , β ( a 1 ) .
 (ii)
F ( E α , β ( a ) ) = E α , β ( F ( a ) ) ; in particular, F ( e a z ) = e z F ( a ) and F ( C o s ( z , a ) ) = C o s ( F ( z ) , a ) for z C .
 (iii) 
σ 1 ( Z ) ( E α , β ( a ) ) = E α , β ( σ 1 ( Z ) ( a ) ) .
 (iv) 
The following Laplace transform formula holds
0 e λ t t α k + β 1 E α , β ( k ) ( t α a ) d t = k ! λ α β ( λ α a ) 1 ( k + 1 ) , ( λ ) > a 1 1 / α ,
for k N { 0 } .
Theorem 2.
Given b 1 ( Z ) and p 1 , we define the operator T b : p ( N 0 ) p ( N 0 ) by
T b ( f ) ( n ) : = b f ( n ) = j 0 b ( n j ) f ( j ) , n N 0 .
 (i) 
T b B ( p ( N 0 ) ) and b p T b B ( p ( N 0 ) ) b 1 .
 (ii) 
σ B ( p ( N 0 ) ) ( T b ) σ 1 ( Z ) ( b ) = F ( b ) ( T ) .
Proof. (i) It is clear that T b B ( p ( N 0 ) ) and T b b 1 . Now, take f = δ n , and T b ( δ n ) ( j ) = b ( j n ) for n , j 0 . Then,
T b ( δ n ) p p = j = 0 | b ( j n ) | p = l = n | b ( l ) | p
for n 0 . We conclude that b p T b B ( p ( N 0 ) ) .
(ii) Now, take λ ρ 1 ( Z ) ( b ) . Then, ( λ b ) 1 1 ( Z ) B ( p ( N 0 ) ) and λ ρ B ( p ( N 0 ) ) ( T b ) . We conclude that σ B ( p ( N 0 ) ) ( T b ) σ 1 ( Z ) ( b ) . □

3. Some Finite Difference Operators in 1 ( Z )

Sequences of compact support, i.e., elements in the set
c c ( Z ) : = { a 1 ( Z ) : m Z + : a ( n ) = 0 , | n | > m ) }
are an important case of finite difference operators. In such a case, the discrete Fourier Transform of a c c ( Z ) is a trigonometric polynomial
F ( a ) ( e i θ ) = j = m m a ( j ) e i j θ .
It is interesting to observe that if j = m m a ( j ) = 0 , then 0 σ 1 ( Z ) ( a ) . This follows immediately from σ 1 ( Z ) ( b ) = F ( b ) ( T ) .
Definition 2.
For f p ( N 0 ) , with 1 p , we define the following operators
1.
  Δ f ( n ) = f ( n ) f ( n + 1 ) = ( ( δ 0 δ 1 ) f ) ( n ) .
2.
  f ( n ) = f ( n ) f ( n 1 ) = ( ( δ 0 δ 1 ) f ) ( n ) .
3.
  D f ( n ) = f ( n + 1 ) f ( n 1 ) = ( ( δ 1 δ 1 ) f ) ( n ) .
4.
  Δ d f ( n ) = f ( n + 1 ) 2 f ( n ) + f ( n 1 ) = ( ( δ 1 2 δ 0 + δ 1 ) f ) ( n ) .
5.
  Δ 2 f ( n ) = f ( n ) 2 f ( n + 1 ) + f ( n + 2 ) = ( ( δ 0 2 δ 1 + δ 2 ) f ) ( n ) .
6.
  2 f ( n ) = f ( n ) 2 f ( n 1 ) + f ( n 2 ) = ( ( δ 0 2 δ 1 + δ 2 ) f ) ( n ) .
7.
  Δ d d f ( n ) = f ( n + 2 ) 2 f ( n ) + f ( n 2 ) = ( ( δ 2 2 δ 0 + δ 2 ) f ) ( n ) .
for n N 0 and f ( n ) = 0 for n < 0 .
The above operators are often used in the context of a numerical analysis. The operators Δ and ∇ are connected with the Euler scheme of approximation. The discrete Laplacian Δ d is the second-order central difference approximation for the second-order derivative. The double Laplacian, the operator Δ d d , is introduced in Bateman’s seminar paper ([12], Page 506), to treat the equations of Born and Karman on crystal lattices in vibration. Other operators Δ 2 , D and 2 are also considered in [12].
To consider the action of these operators in B ( p ( N 0 ) ) , we study these operators as elements in the Banach algebra 1 ( Z ) as Theorem 2 shows. Operators Δ , ∇, Δ d and Δ d d have been studied in detail in ([8], Theorem 3.2, 3.3, 3.4 and 3.5). In the following subsections, we treat D , Δ 2 and 2 .
In Table 1, we collect some basic results of the finite difference operators given in Definition 2. In Figure 1, we also plot the spectrum σ 1 ( Z ) ( a ) for these finite difference operators.

3.1. The Operator D

This operator is a finite difference operator of order 1 defined by D f ( n ) = f ( n + 1 ) f ( n 1 ) . We present some of these properties in the next theorem.
Theorem 3.
The operator D f = a f with a = δ 1 δ 1 , has the following properties.
 (i) 
The norm of a is equal to 2, a 1 = 2 .
 (ii) 
The discrete Fourier transform of a is given by F ( a ) ( z ) = z 1 z with z T .
 (iii) 
The spectrum of a is σ 1 ( Z ) ( D ) = [ 2 i , 2 i ] .
 (iv) 
The group generated by a is e z a ( n ) = J n ( 2 z ) , z C , n Z and e t a 1 e 2 t , with t R and t > 0 .
 (v) 
For λ C \ [ 2 i , 2 i ] ,
( λ δ 0 + a ) 1 = 2 n ( λ 2 + 4 λ ) n λ 2 + 4 , n Z .
Proof. Items (i) and (ii) are straightforward. To show (iii), we have that
σ 1 ( Z ) ( D ) = F ( a ) ( T ) = { z C : z = w 1 w , | w | = 1 } = { z C : z = 2 i ( w ) , | w | = 1 } = [ 2 i , 2 i ] .
We define for a while e z 1 ( Z ) with z C by e z ( n ) = J n ( 2 z ) . Note that
e z 1 = n Z | J n ( 2 z ) | n Z I n ( 2 | z | ) = e 2 | z | .
where we have applied the generating function of modified Bessel function of the first kind I n with n Z , see, for example, [13], [Appendix].
We apply the discrete Fourier transform to obtain
F ( e z ) ( w ) = n Z J n ( 2 z ) w n = e z 1 w w = e z F ( a ) ( w ) = F ( e z a ) ( w ) ,
with w T . As the discrete Fourier transform is injective, we conclude that e z ( n ) = e z a ( n ) and the item (iv) is proved.
Finally, to show item (v), we have that
( λ δ 0 + a ) 1 ( n ) = 0 e λ t e t a ( n ) d t = 0 e λ t J n ( 2 t ) d t = 2 n ( λ 2 + 4 λ ) n λ 2 + 4 ,
where we have applied ([14], Formula 6.623) for ( λ ) > 0 and n Z . By analytic prolongation, we extend the equality for λ C \ [ 2 i , 2 i ] . □

3.2. The Operator Δ 2

The operator Δ 2 f ( n ) = f ( n ) 2 f ( n + 1 ) + f ( n + 2 ) is a finite difference operator of order 2. Note that Δ 2 = ( Δ ) 2 . In the next theorem, we present some properties of Δ 2 .
Theorem 4.
The operator Δ 2 f = a f with a = δ 0 2 δ 1 + δ 2 verifies the following properties.
 (i) 
The norm of a is equal to 4, a 1 = 4 .
 (ii) 
The discrete Fourier transform F ( a ) ( z ) = ( z 1 ) 2 with z T .
 (iii) 
The spectrum of a in 1 ( Z ) is equal to σ 1 ( Z ) ( Δ 2 ) = { z C : z = ( w 1 ) 2 , w T } .
 (iv) 
The group generated by a is
e z a ( n ) = ( i z ) n H n ( i z ) ( n ) ! e z χ N 0 ( n ) , z C ,
n Z and e t a 1 e 4 t for t > 0 .
 (v) 
The cosine function C o s ( z , a ) = z n ( n ) ! e z + ( 1 ) n e z 2 χ N 0 ( n ) .
Proof. Items (i), (ii) and (iii) are similar to items (i), (ii) and (iii) in Theorem 3. Now, we define e z 1 ( Z ) con z C by e z ( n ) ( i z ) n H n ( i z ) ( n ) ! e z χ N 0 ( n ) where H n is the Hermite polynomial. First, we check that e z 1 ( Z ) ,
| | e z | | 1 = e ( z ) j = 0 | ( i z ) j H j ( i z ) j ! | e | z | j = 0 | z | j j ! | H j ( i z ) | = e | z | j = 0 | z | j j ! | 2 j π | i n ( z + t ) n e t 2 d t e | z | j = 0 | z | j j ! 2 j π ( | z | + | t | ) j e t 2 d t = 2 e | z | π j = 0 ( 2 | z | ) j j ! 0 ( | z | + t ) j e t 2 d t = 2 e | z | π 0 e t 2 j = 0 ( 2 | z | ( | z | + t ) ) j j ! d t = 2 e | z | π 0 e t 2 e 2 | z | ( | z | + t ) d t = 2 e 4 | z | π 0 e ( t | z | ) 2 d t = 2 e 4 | z | π | z | e u 2 d u < 2 e 4 | z | π e u 2 d u = 2 e 4 | z | < .
We calculate the discrete Fourier transform of e z
F ( e z ) ( w ) = n Z ( i z ) n H n ( i z ) ( n ) ! e z χ N 0 ( n ) w n = e z j = 0 H j ( i z ) j ! ( w i z ) j = e z e ( w i z ) 2 + 2 ( w i z ) ( i z ) = e z ( w 1 ) 2 = e z F ( a ) ( w ) = F ( e z a ) ( w ) ,
where we have applied ([14], Formula 8.951), for w T . As the discrete Fourier transform is injective, we conclude that e z ( n ) = e z a ( n ) for n Z .
To calculate the cosine function generated by a, we apply ([10], Example 3.14.15)
C o s ( z , a ) = 1 2 ( e z ( δ 0 δ 1 ) + e z ( δ 0 δ 1 ) ) = 1 2 e z z n ( n ) ! χ N 0 ( n ) + e z ( z ) n ( n ) ! χ N 0 ( n ) = z n ( n ) ! e z + ( 1 ) n e z 2 χ N 0 ( n ) ,
and we conclude the proof.□

3.3. The Operator 2

In this subsection, we treat the finite difference operator of order 2, defined by 2 f ( n ) = f ( n ) 2 f ( n 1 ) + f ( n 2 ) . Note that e 2 = 2 . We present some basic properties of operator 2 .
Theorem 5.
The operator 2 f = a f with a = δ 0 2 δ 1 + δ 2 has the following properties.
 (i) 
The norm of a in 1 ( Z ) is equal to 4, a 1 = 4 .
 (ii) 
The discrete Fourier transform F ( a ) ( z ) = ( 1 z 1 ) 2 for z T .
 (iii) 
The spectrum of a in 1 ( Z ) is equal to σ 1 ( Z ) ( 2 ) = { z C : z = ( 1 w 1 ) 2 , w T } .
 (iv) 
The group generated by a is given by
e z a ( n ) = ( i z ) n H n ( i z ) ( n ) ! e z χ N 0 ( n ) , z C ,
n Z and e t a 1 e 4 t , for t > 0 .
 (v) 
The cosine function C o s ( z , a ) ( n ) = z n n ! e z + ( 1 ) n e z 2 χ N 0 ( n ) .
Proof. We skip the proof of items (i), (ii) and (iii). To show (iv), we define e z ( n ) = ( i z ) n H n ( i z ) n ! e z χ N 0 ( n ) for z C . As in Theorem 4, e z 1 ( Z ) and
F ( e z ) ( w ) = n Z ( i z ) n H n ( i z ) n ! e z χ N 0 ( n ) w n = e z n = 0 H n ( i z ) n ! i z w n = e z e ( i z w ) 2 + 2 ( i z w ) ( i z ) = e z ( 1 w 1 ) 2 = e z F ( a ) ( w ) = F ( e z a ) ( w ) ,
for w T . We conclude that e z ( n ) = e z a ( n ) for z C . As 2 = ( ) 2 , we apply again ([10], Example 3.14.15), to obtain the cosine function generated by a. □
Remark 1.
The connection between semigroups and cosine functions is given by the so-called Weierstrass formula,
e a t = 1 π t 0 e s 2 4 t C o s ( s ) d s , t > 0 ,
for a 1 ( Z ) , see, for example, ([10], Theorem 3.14.17). In the case that we apply the Weierstrass formula in the conditions of Theorems 4 and 5 to obtain the well-known formula
( 2 i ) n π H n ( i ω ) = x n e ( x ω ) 2 d x , ω C ,
see, for example, ([14], Formula 3.462(4)).

4. Fundamental Solution for Semidiscrete Evolution Equations

In this section, we consider the operator B f ( n ) : = ( b f ) ( n ) , with b 1 ( Z ) , f p ( N 0 ) , p [ 1 , ] and n N 0 . We obtain an explicit representation of solutions for the following time/space fractional evolution equation:
D t β u ( n , t ) = B u ( n , t ) + g ( n , t ) , n N 0 , t > 0 . u ( n , 0 ) = φ ( n ) , u t ( n , 0 ) = ϕ ( n ) n N 0 ,
Here, β ( 0 , 2 ] is real number. We recall that D t β denotes the Caputo fractional derivative given by
D t β v ( t ) = 1 Γ ( 1 β ) 0 t ( t s ) β v ( s ) d s = ( g 1 β v ) ( t ) , t > 0 ,
for 0 < β < 1 and
D t β v ( t ) = 1 Γ ( 2 β ) 0 t ( t s ) 1 β v ( s ) d s = ( g 2 β v ) ( t ) , t > 0 ,
for 1 < β < 2 . For β = 1 and β = 2 , we consider the usual first- and second-order derivation. Note that
lim β 1 D t β v ( t ) = v ( t ) , lim β 2 D t β v ( t ) = v ( t ) , t > 0 ,
however,
lim β 0 + D t β v ( t ) = v ( t ) v ( 0 ) lim β 1 + D t β v ( t ) = v ( t ) v ( 0 ) , t > 0 ,
see, for example, [15,16].
The main result in this section is the following Theorem. The function E α , β ( b ) (with b 1 ( Z ) ) is the vector-valued Mittag-Leffler function given in Definition 1. A similar result was stated in the Banach space p ( Z ) in ([8], Theorem 5.1).
Theorem 6.
Let φ , ϕ p ( N 0 ) , and g : N 0 × R + C be such that, for each t R + , g ( · , t ) p ( N 0 ) and sup s [ 0 , t ] | | g ( · , s ) | | p < with 1 p .
 (i) 
For 0 < β < 1 , the function
u ( n , t ) = ( E β , 1 ( t β b ) φ ) ( n ) + 0 t ( t s ) β 1 E β , β ( ( t s ) β b ) g ( · , s ) ( n ) d s ,
is the unique solution of the initial value problem (4) with u ( · , t ) belonging to p ( N 0 ) for t > 0 and n N 0 .
 (ii) 
For 1 < β < 2 , the function
u ( n , t ) = ( E β , 1 ( t β b ) φ ) ( n ) + t ( E β , 2 ( t β b ) ϕ ) ( n ) + 0 t ( t s ) β 1 E β , β ( ( t s ) β b ) g ( · , s ) ( n ) d s ,
is the unique solution of the initial value problem (5) with u ( · , t ) belonging to p ( N 0 ) for t > 0 and n N 0 .
Proof. We show part (ii) in the case p = . Part (i) or 1 p < are proved in a similar way. We prove the result in several steps.
Step 1. First, we show the explicit solution for (5). Taking the Z -transform of (5), we obtain that
D t β u ( z , t ) = B u ( z , t ) + g ( z , t ) u ( z , 0 ) = φ ( z ) ; u t ( z , 0 ) = ϕ ( z ) .
Now, taking Laplace’s transformation to (13), we have:
u ^ ( z , λ ) = λ β 1 λ β F ( b ) ( θ ) φ ( z ) + λ β 2 λ β F ( b ) ( θ ) ϕ ( z ) + 1 λ β F ( b ) ( θ ) g ^ ( z , λ ) .
By inverse Laplace transform, see identity (9), we obtain
u ( z , t ) = E β , 1 ( F ( b ) ( θ ) t β ) φ ( z ) + t E β , 2 ( F ( b ) ( θ ) t β ) ϕ ( z ) + 0 t ( t s ) β 1 E β , β ( ( t s ) β F ( b ) ( θ ) ) ( z ) g ( z , s ) d s .
We apply Proposition 1 (ii) to obtain
u ( n , t ) = ( E β , 1 ( t β b ) φ ) ( n ) + t ( E β , 2 ( t β b ) ϕ ) ( n ) + 0 t ( t s ) β 1 E β , β ( ( t s ) β b ) g ( · , s ) ( n ) d s ,
for n N 0 and t > 0 .
Step 2. Now, we prove uniqueness. Suppose that the system (5) has two solutions u 1 and u 2 with the same initial values φ , ϕ , and write v : = u 1 u 2 . Then, v is a solution of the following ODE
D t β v ( z , t ) = B v ( z , t ) , v ( z , 0 ) = 0 , v t ( z , 0 ) = 0 .
Due to the above ordinary differential equation having its unique solution and that the function zero is a solution, we conclude that v ( z , t ) = 0 . As the Z -transform is injective, we conclude that v ( n , t ) = 0 for every n N 0 and t 0 . Hence, u 1 = u 2 .
Step 3. By Proposition 1 (i), we obtain that
E β , 1 ( t β b ) 1 E β , 1 ( t β b 1 ) E β , 2 ( t β b ) 1 E β , 2 ( t β b 1 ) E β , β ( t s ) β b ) 1 E β , 2 ( ( t s ) β b 1 )
for 0 < s < t and b 1 ( Z ) . Because sup s [ 0 , t ] | | g ( · , s ) | | < , we conclude that
0 t ( t s ) β 1 E β , β ( ( t s ) β b ) g ( · , s ) ( n ) d s ( N 0 ) ,
and the solution u ( · , t ) ( N 0 ) for t > 0 . □
Remark 2.
Now, we may shortly treat the behavior of the solution when β tends to natural parameter, i.e., β = 1 , 2 . For simplicity, we only present the homogeneous case, g = 0 . In the case that β 1 , the solution of Equation (4) tends to semigroup family operators E 1 , 1 ( t b ) , and when β 2 , the solution of Equation (5),
u ( · , t ) = E β , 1 ( t β b ) φ + t E β , 2 ( t β b ) ϕ , t > 0
tends to the well-known solution of second-order Cauchy problem, expressed in terms of cosine function and sine function generated by b, see ([10], Corollary 3.14.8).
However, in the case that β 1 + , the solution of Equation (5) converges to
u ( · , t ) = E 1 , 1 ( b t ) + t E 1 , 2 ( t b ) , t > 0 ,
as in the scalar case. We remark that this function is the solution of the following first-order modified Cauchy problem
v ( n , t ) = B v ( n , t ) + ϕ ( n ) , n N 0 , t > 0 , v ( n , 0 ) = φ ( n ) , n N 0 ,
for ϕ , φ p ( N 0 ) . This natural fact is connected with the interpolation property of the Caputo fractional differentiation, see (12).
The fundamental solution u β , 1 for Equations (4) and (5) is given by taking as initial values ψ = δ 0 and ϕ = 0 . In the case 1 < β 2 ( β = 2 , is the wave equation), a second fundamental solution u β , 2 is given by ψ = 0 and ϕ = δ 0 , see ([17], Remark 3.2).
As a corollary of Theorems 2 and 6 is the following subordination theorem for fundamental solutions. This result extends ([17], Corollary 3.5) in the space p ( N 0 ) .
We denote by Φ α the entire Wright function defined by
Φ α ( z ) : = n = 0 ( z ) n n ! Γ ( α n + 1 α ) = 1 2 π i γ μ α 1 e μ z μ α d μ , 0 < α < 1 ,
where γ is a complex path which starts and ends at and rounds the origin once counterclockwise. The Wright function is a known special function which has appeared in a wide variety of different contexts, for example, it is used for models in stochastic processes, see [18]. The proof of the next corollary is similar to ([8], Corollary 5.3), and we leave for the reader.
Corollary 1.
Let u β , 1 and u β , 2 be the fundamental solutions of problems (4) and (5) and Ψ α the Wright function defined by (15).
 (i) 
Let 0 < β < 1 . Then,
u β , 1 ( n , t ) = 0 Φ β ( τ ) u 1 , 1 ( n , τ t β ) d τ , n N 0 , t > 0 .
 (ii) 
Let 1 < β < 2 . Then,
u β , 1 ( n , t ) = 0 Φ β 2 ( τ ) u 2 , 1 ( n , τ t β 2 ) d τ , u β , 2 ( n , t ) = 0 t ( t u ) β 2 Γ ( 1 β 2 ) 0 Φ β 2 ( τ ) u 2 , 2 ( n , τ u β 2 ) d τ d u ,
for n N 0 and t > 0 .

5. Conclusions and Future Work

In this paper, we have considered Caputo fractional differential equations in the sequence Lebesgue spaces p ( N 0 ) with p 1 . The associate operator is given by a convolution with a sequence in the Banach Algebra which belongs to 1 ( Z ) . We use techniques from the Functional Analysis, such as the Guelfand theory in Banach algebra, to obtain more information about the problem. We calculate the explicit solution in terms of Mittag-Leffer functions. Some particular cases (sequences of compact support) of finite difference operators are treated in detail.
An interesting problem to address in the future is to consider these techniques in the continuous case. We may study these Caputo fractional differential equations in L p ( R + ) and L p ( R ) .

Author Contributions

Conceptualization, A.M. and P.J.M.; methodology, A.M. and P.J.M.; software, A.M. and P.J.M.; investigation, A.M. and P.J.M.; writing, A.M. and P.J.M.; funding acquisition, A.M. and P.J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study did not require ethical approval.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors have been partially supported by ProjectID2019-105979GBI00, DGI-FEDER, of the MCEI and Project E48-20R, Gobierno de Aragón, Spain.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Friesl, M.; Slavik, A.; Stehlik, P. Discrete-space partial dynamic equations on time scales and applications to stochastic processes. Appl. Math. Lett. 2014, 37, 86–90. [Google Scholar] [CrossRef]
  2. Slavik, A. Mixing problems with many tanks. Am. Math. Mon. 2013, 120, 806–821. [Google Scholar] [CrossRef] [Green Version]
  3. Feintuch, A.; Francis, B. Infinite chains of kinematic points. Autom. J. IFAC 2012, 48, 901–908. [Google Scholar] [CrossRef]
  4. Slavik, A. Asymptotic behavior of solutions to the semidiscrete diffusion equation. Appl. Math. Lett. 2020, 106, 106392. [Google Scholar] [CrossRef]
  5. Abadias, L.; de León-Contreras, M.; Torrea, J.L. Non-local fractional derivatives. Discrete and continuous. J. Math. Anal. Appl. 2017, 449, 734–755. [Google Scholar] [CrossRef] [Green Version]
  6. Ciaurri, O.; Lizama, C.; Roncal, L.; Varona, J.L. On a connection between the discrete fractional Laplacian and superdiffusion. Appl. Math. Lett. 2015, 49, 119–125. [Google Scholar] [CrossRef]
  7. Lizama, C.; Roncal, L. Hölder-Lebesgue regularity and almost periodicity for semidiscrete equations with a fractional Laplacian. Discr. Cont. Dyn. Syst. Ser. A 2018, 38, 1365–1403. [Google Scholar] [CrossRef] [Green Version]
  8. González-Camus, J.; Lizama, C.; Miana, P.J. Fundamental solutions for semidiscrete evolution equations via Banach algebras. Adv. Differ. Equ. 2021, 2021, 35. [Google Scholar] [CrossRef] [PubMed]
  9. Larsen, R. Banach Algebras: An Introduction; Marcel Dekker: New York, NY, USA, 1973. [Google Scholar]
  10. Arendt, W.; Batty, C.; Hieber, M.; Neubrander, F. Vector-Valued Laplace Transforms and Cauchy Problems; Monographs in Mathematics; Birkhäuser-Verlag: Basel, Switzerland, 2001; Volume 96. [Google Scholar]
  11. Sinclair, A.M. Continuous Semigroups in Banach Algebras; London Mathematical Society, Lecture Note Series; Cambridge University Press: Cambridge, UK, 1982; Volume 63. [Google Scholar]
  12. Bateman, H. Some simple differential difference equations and the related functions. Bull. Am. Math. Soc. 1943, 49, 494–512. [Google Scholar] [CrossRef] [Green Version]
  13. Ciaurri, O.; Gillespie, T.A.; Roncal, L.; Torrea, J.L.; Varona, J.L. Harmonic analysis associated with a discrete Laplacian. J. d’Anal. MathéMatique 2017, 132, 109–131. [Google Scholar] [CrossRef]
  14. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products, 7th ed.; Elsevier: Amsterdam, The Netherlands; Academic Press: Burlington, MA, USA; San Diego, CA, USA; London, UK, 2007. [Google Scholar]
  15. Bazlekova, E.G. Fractional Evolution Equations in Banach Spaces; Technische Universiteit: Eindhoven, The Netherlands, 2001. [Google Scholar]
  16. Gorenflo, R.; Mainardi, F. Parametric Subordination in Fractional Diffusion Processes. In Fractional Dynamics, Recent Advances; Klafter, J., Lim, S.C., Metzler, R., Eds.; World Scientific: Singapore, 2012; pp. 229–263. [Google Scholar]
  17. González-Camus, J.; Keyantuo, V.; Lizama, C.; Warma, M. Fundamental solutions for discrete dynamical systems involving the fractional Laplacian. Math. Meth. Appl. Sci. 2019, 42, 1–24. [Google Scholar] [CrossRef]
  18. Gorenflo, R.; Luchko, Y.; Mainardi, F. Analytical properties and applications of the Wright function. Fract. Calc. Appl. Anal. 1999, 2, 383–414. [Google Scholar]
Figure 1. Spectrum of finite difference operators in 1 ( Z ) .
Figure 1. Spectrum of finite difference operators in 1 ( Z ) .
Foundations 02 00059 g001
Table 1. Finite difference operators in 1 ( Z ) .
Table 1. Finite difference operators in 1 ( Z ) .
A A ( f ) = a f a a 1 F ( a ) ( z )
Δ f ( n ) f ( n + 1 ) δ 0 δ 1 2 1 z
f ( n ) f ( n 1 ) δ 0 δ 1 2 1 1 z
D f ( n + 1 ) f ( n 1 ) δ 1 δ 1 2 z 1 z
Δ d f ( n + 1 ) 2 f ( n ) + f ( n 1 ) δ 1 2 δ 0 + δ 1 4 z + 1 z 2
Δ 2 f ( n + 2 ) 2 f ( n + 1 ) + f ( n ) δ 0 2 δ 1 + δ 2 4 ( z 1 ) 2
2 f ( n ) 2 f ( n 1 ) + f ( n 2 ) δ 0 2 δ 1 + δ 2 4 ( 1 z 1 ) 2
Δ d d f ( n + 2 ) 2 f ( n ) + f ( n 2 ) δ 2 2 δ 0 + δ 2 4 ( z 1 z ) 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mahillo, A.; Miana, P.J. Caputo Fractional Evolution Equations in Discrete Sequences Spaces. Foundations 2022, 2, 872-884. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2040059

AMA Style

Mahillo A, Miana PJ. Caputo Fractional Evolution Equations in Discrete Sequences Spaces. Foundations. 2022; 2(4):872-884. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2040059

Chicago/Turabian Style

Mahillo, Alejandro, and Pedro J. Miana. 2022. "Caputo Fractional Evolution Equations in Discrete Sequences Spaces" Foundations 2, no. 4: 872-884. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2040059

Article Metrics

Back to TopTop