Next Article in Journal
Some Remarks about Exponential and Semi-Exponential Operators
Next Article in Special Issue
Bessel Collocation Method for Solving Fredholm–Volterra Integro-Fractional Differential Equations of Multi-High Order in the Caputo Sense
Previous Article in Journal
Analysis and Prediction for Confirmed COVID-19 Cases in Czech Republic with Uncertain Logistic Growth Model
Previous Article in Special Issue
Analysis and Optimal Control of φ-Hilfer Fractional Semilinear Equations Involving Nonlocal Impulsive Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectrum of Fractional and Fractional Prabhakar Sturm–Liouville Problems with Homogeneous Dirichlet Boundary Conditions

by
Malgorzata Klimek
Department of Mathematics, Faculty of Mechanical Engineering and Computer Science, Czestochowa University of Technology, 42-201 Czestochowa, Poland
Symmetry 2021, 13(12), 2265; https://doi.org/10.3390/sym13122265
Submission received: 26 October 2021 / Revised: 15 November 2021 / Accepted: 18 November 2021 / Published: 28 November 2021
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus)

Abstract

:
In this study, we consider regular eigenvalue problems formulated by using the left and right standard fractional derivatives and extend the notion of a fractional Sturm–Liouville problem to the regular Prabhakar eigenvalue problem, which includes the left and right Prabhakar derivatives. In both cases, we study the spectral properties of Sturm–Liouville operators on function space restricted by homogeneous Dirichlet boundary conditions. Fractional and fractional Prabhakar Sturm–Liouville problems are converted into the equivalent integral ones. Afterwards, the integral Sturm–Liouville operators are rewritten as Hilbert–Schmidt operators determined by kernels, which are continuous under the corresponding assumptions. In particular, the range of fractional order is here restricted to interval ( 1 / 2 , 1 ] . Applying the spectral Hilbert–Schmidt theorem, we prove that the spectrum of integral Sturm–Liouville operators is discrete and the system of eigenfunctions forms a basis in the corresponding Hilbert space. Then, equivalence results for integral and differential versions of respective eigenvalue problems lead to the main theorems on the discrete spectrum of differential fractional and fractional Prabhakar Sturm–Liouville operators.

1. Introduction

The aim of this paper is to study the fundamental properties of fractional eigenvalue problems developed by the construction of the Sturm–Liouville operator (SLO) with left and right fractional derivatives. In classical differential equations theory, this is a linear differential operator of the second order and yields an eigenvalue problem of the form (here, x [ 0 , b ] in the case when we consider the problem on a finite interval):
L q y ( x )   =   d d x p ( x ) d y ( x ) d x + q ( x ) y ( x )   =   λ w ( x ) y ( x )
with boundary conditions appearing as follows:
c 1 y ( 0 ) + c 2 d y ( 0 ) d x   =   0 , d 1 y ( b ) + d 2 d y ( b ) d x   =   0 .
Let us point out that, depending on the choice of coefficient functions and boundary conditions, such problems provide various systems of orthogonal eigenfunctions, orthogonal polynomials and families of special functions. Orthogonal systems of the solutions of classical Sturm–Liouville problems are widely applied in the analysis and solving of fundamental differential equations of mathematics, physics, mechanics , and economics.
In most of the FSLPs presented at the beginning of fractional Sturm–Liouville theory, first-order derivatives in a standard Sturm–Liouville problem were replaced with fractional order derivatives. The resulting equations were solved using some numerical schemes [1,2,3,4]. However, in these works, the essential properties, such as the orthogonality of the eigenfunctions of the fractional operator, were not investigated. In addition, the question of whether the associated eigenvalues are real or not is not addressed. Some results concerning these properties have been obtained in papers [5,6], where the discussed equations contain a classical SLO extended by including a sum of the left and the right derivatives. Then, in paper [7], we proposed the construction of a fractional Sturm–Liouville operator which preserves the orthogonality of the eigenfunctions corresponding to distinct eigenvalues and provides real eigenvalues. The FSLO contains both the left and right derivatives and is a symmetric operator on function space restricted by fractional boundary conditions which generalize conditions (1).
A fractional version of Bessel SLO has been developed and applied to anomalous diffusion in [8], where the space-fractional differential operator has a form analogous to the FSLO proposed in a general form in [7]. Some special cases of singular fractional Sturm–Liouville problems were also studied in [9,10], where exact solutions and eigenvalues were calculated.
In our earlier works [7,11,12,13,14], we focused on the construction of a fractional version of operator L q , which includes standard fractional derivatives. The characteristic feature of the proposed approach is the mixture of the left and right fractional derivatives in the fractional Sturm–Liouville operator (FSLO). This construction provides eigenvalue problems with orthogonal eigenfunctions and discrete spectra under the appropriate homogeneous boundary conditions.
In recent years, fractional eigenvalue problems have also been discussed within the framework of tempered and conformable fractional calculus. In the papers [15,16], a fractional Sturm–Liouville operator is built by using the left and right tempered derivatives. Next, in [17,18], an FSLO is constructed as a composition of conformable fractional derivatives. In addition, in paper [19], the authors show how to build an FSLO with composite fractional derivatives.
Here, we add the generalization of fractional eigenvalue problems to problems with operators, including Prabhakar derivatives. The regular fractional and fractional Prabhakar Sturm–Liouville operators considered here include the left and the right derivatives, and the derived equations are in fact of a variational nature; i.e., they are Euler–Lagrange equations for respective actions (compare [11,20] and the references therein for FSLE). The properties of the spectra and eigenfunctions’ systems of FSLP can be studied by applying the variational method [12,21]. Here, we shall develop the transformation method for FSLP and PSLP with Dirichlet boundary conditions, which means that we rewrite the FSLP/PSLP as the equivalent integral eigenvalue problem.
The paper is organized as follows. In the next section, we present the necessary definitions and properties of fractional and fractional Prabhakar operators, as well as the formulation of a regular fractional Sturm–Liouville problem with its generalization to the Prabhakar Sturm–Liouville problem. In Section 3, we define the problems with homogeneous Dirichlet boundary conditions and derive equivalence results for both types of fractional eigenvalue problems. It appears that by applying composition rules for derivatives and integrals, they can be converted into the equivalent integral ones. Spectral properties of integral versions of fractional and fractional Prabhakar Sturm–Liouville operators are discussed in Section 4. We shall prove that these operators are Hilbert–Schmidt integral operators, which are compact and self-adjoint on the L w 2 ( 0 , b ) space. Applying the spectral Hilbert–Schmidt theorem, we derive results on discrete spectra both for fractional and fractional Prabhakar Sturm–Liouville operators. The equivalence of differential and integral versions of eigenvalue problems leads to the corresponding spectral results for differential operators.
The paper closes with a brief discussion of results and future investigations. The Appendix A contains two parts. First, we present results on Hölder continuity of kernels defining integral Sturm–Liouville operators. Then, we prove a useful theorem on the convergence of convolutions’ series in a general case, which is applied in the construction of integral Sturm–Liouville operators.

2. Preliminaries

We start with a summary of definitions and properties of fractional integrals and derivatives which shall be applied in the construction of fractional and fractional Prabhakar eigenvalue problems. First, we recall the left and right Riemann–Liouville fractional derivatives of order α 0 , 1 [22,23]:
D 0 + α y x   : =   d d x I 0 + 1 α y x , D b α y x   : =   d d x I b 1 α y x ,
where the operators I 0 + α and I b α are respectively the left and the right fractional Riemann–Liouville integrals of order α > 0 defined by the following formulas
I 0 + α y x   : =   0 x ( x t ) α 1 y ( t ) Γ α d t , x > 0 ,
I b α y x   : =     x b ( t x ) α 1 y ( t ) Γ α d t , x < b .
Next, we have Caputo fractional derivatives:
c D 0 + α y ( x )   =   D 0 + α y ( x ) y ( 0 ) , c D b α y ( x )   =   D b α y ( x ) y ( b )
and we note that when y ( 0 )   =   y ( b )   =   0 , both types of derivatives coincide, i.e.,
c D 0 + α y ( x )   =   D 0 + α y ( x ) , c D b α y ( x )   =   D b α y ( x ) .
We also recall some of the composition rules of fractional operators for the case of order α 0 , 1 ; namely, for the left-sided Caputo derivative and left-sided fractional integral, we have
I 0 + α c D 0 + α y x     =     y x y 0 ,
c D 0 + α I 0 + α y x   =   y x ,
while for the right-sided Riemann–Liouville derivatives, the following relations are valid
I b α D b α y x   =   y x I b 1 α y b · ( b x ) α 1 Γ ( α ) ,
D b α I b α y x   =   y x .
All of the above rules are fulfilled for all points x [ 0 , b ] when function y is a continuous one. Let us note that for the continuous function fulfilling condition y ( 0 )   =   0 , rules (6) and (8) look as follows:
I 0 + α c D 0 + α y x   =   y x , I b α D b α y x   =   y x .
The fractional operators, described above, are generalized to Prabhakar integrals and derivatives. They are defined using a three-parameter Mittag–Leffler function [22,24]:
E ρ , μ γ ( z )   : =   1 Γ ( γ ) k   =   0 Γ ( γ + k ) Γ ( ρ k + μ ) · z k k !
and Prabhakar function [24,25]:
e ρ , μ γ ( ω z ρ )   : =   z μ 1 E ρ , μ γ ( ω z ρ ) ,
both defined on the complex space when R e ( ρ ) > 0 and R e ( μ ) > 0 .
These functions lead to the left and right Prabhakar derivatives [24]:
D ρ , γ , ω , 0 + α y x   : =   d d x E ρ , γ , ω , 0 + 1 α y x , D ρ , γ , ω , b α y x   : =   d d x E ρ , γ , ω , b 1 α y x ,
where operators E ρ , γ , ω , 0 + α and E ρ , γ , ω , b α are respectively the left and the right fractional Prabhakar integrals:
E ρ , γ , ω , 0 + α y x   : =   0 x e ρ , α γ ( ω ( x t ) ρ ) y ( t ) d t , x > 0 ,
E ρ , γ , ω , b α y x   : =   x b e ρ , α γ ( ω ( x t ) ρ ) y ( t ) d t , x < b .
Similar to Caputo derivatives, given in (5), we have Caputo-type Prabhakar derivatives defined as follows
c D ρ , γ , ω , 0 + α y x   =   D ρ , γ , ω , 0 + α y ( x ) y ( 0 ) ,
c D ρ , γ , ω , b α y x   =   D ρ , γ , ω , b α y ( x ) y ( b )
coinciding with Prabhakar derivatives (13) when y ( 0 )   =   0 or y ( b )   =   0 , respectively. Restricting function space to continuous functions fulfilling condition y ( 0 )   =   0 , we arrive at composition rules of Prabhakar operators analogous to (7), (9), and (10):
c D ρ , γ , ω , 0 + α E ρ , γ , ω , 0 + α y ( x )   =   y ( x ) ,
E ρ , γ , ω , 0 + α c D ρ , γ , ω , 0 + α y ( x )   =   y ( x ) ,
D ρ , γ , ω , b α E ρ , γ , ω , b α y ( x )   =   y ( x ) ,
E ρ , γ , ω , b α D ρ , γ , ω , b α y ( x )   =   y ( x ) .
Now, we shall quote the general formulation of the fractional eigenvalue problem, introduced and investigated in papers [7,11,12,13,14,21].
Definition 1
(compare Definition 5 in [7]). Let α ( 0 , 1 ] . With the notation
L q   : =   D b α p ( x ) c D 0 + α + q ( x ) ,
consider the fractional Sturm–Liouville equation (FSLE)
L q y λ ( x )   =   λ w ( x ) y λ ( x ) ,
where p ( x ) 0 , w ( x ) > 0 x [ 0 , b ] , functions p , q , w are real-valued continuous functions in [ 0 , b ] and boundary conditions are:
c 1 y λ ( 0 ) + c 2 I b 1 α p ( x ) D 0 + α y λ ( x ) x = 0 = 0 ,
d 1 y λ ( b ) + d 2 I b 1 α p ( x ) D 0 + α y λ ( x ) x = b = 0
with c 1 2 + c 2 2 0 and d 1 2 + d 2 2 0 . The problem of finding number λ (eigenvalue) such that the BVP has a non-trivial solution, y λ (eigenfunction) will be called the regular fractional Sturm–Liouville eigenvalue problem (FSLP).
We include Prabhakar derivatives into the construction of FSLO and formulate below the Prabhakar Sturm–Liouville problem.
Definition 2.
Let α ( 0 , 1 ] . With the notation
L q   : =   D ρ , γ , ω , b α p ( x ) c D ρ , γ , ω , 0 + α + q ( x ) ,
consider the fractional Prabhakar Sturm–Liouville equation (PSLE)
L q y λ ( x )   =   λ w ( x ) y λ ( x ) ,
where p ( x ) 0 , w ( x ) > 0 x [ 0 , b ] , functions p , q , w are real-valued continuous functions in [ 0 , b ] and boundary conditions are:
c 1 y λ ( 0 ) + c 2 E ρ , γ , ω , b 1 α p ( x ) D ρ , γ , ω , 0 + α y λ ( x ) x = 0 = 0 ,
d 1 y λ ( b ) + d 2 E ρ , γ , ω , b 1 α p ( x ) D ρ , γ , ω , 0 + α y λ ( x ) x = b = 0
with c 1 2 + c 2 2 0 and d 1 2 + d 2 2 0 . The problem of finding number λ (eigenvalue) such that the BVP has a non-trivial solution, y λ (eigenfunction) will be called the regular fractional Prabhakar Sturm–Liouville eigenvalue problem (PSLP).

3. Formulation of the Problem and Methods

In this section, we shall focus on fractional eigenvalue problems subjected to the homogeneous Dirichlet boundary conditions. We choose values c 2   =   d 2   =   0 in Definitions 1 and 2 and formulate the corresponding definitions of FSLP and PSLP. First, we have the fractional Sturm–Liouville problem with Dirichlet boundary conditions.
Definition 3.
Let α ( 0 , 1 ] . With the notation
L q   : =   D b α p ( x ) c D 0 + α + q ( x ) ,
consider the fractional Sturm–Liouville Equation (23), where p ( x ) 0 , w ( x ) > 0 x [ 0 , b ] , functions p , q , w are real-valued continuous functions in [ 0 , b ] and the boundary conditions are:
y λ ( 0 ) = y λ ( b ) = 0 .
The problem of finding number λ (eigenvalue) such that the BVP has a non-trivial solution, y λ (eigenfunction) will be called the regular fractional Sturm–Liouville eigenvalue problem (FSLP) with homogeneous Dirichlet boundary conditions.
Next, we formulate the definition of the Prabhakar Sturm–Liouville problem with Dirichlet boundary conditions.
Definition 4.
Let α ( 0 , 1 ] . With the notation
L q   : =   D ρ , γ , ω , b α p ( x ) c D ρ , γ , ω , 0 + α + q ( x ) ,
consider the fractional Prabhakar Sturm–Liouville Equation (27), where p ( x ) 0 , w ( x ) > 0 x [ 0 , b ] , functions p , q , w are real-valued continuous functions in [ 0 , b ] and the boundary conditions are:
y λ ( 0 )   =   y λ ( b )   =   0 .
The problem of finding number λ (eigenvalue) such that the BVP has a non-trivial solution, y λ (eigenfunction) is the regular fractional Prabhakar Sturm–Liouville eigenvalue problem (PSLP) with homogeneous Dirichlet boundary conditions.
We shall study the spectral properties of the eigenvalue problems described in the above definitions. Let us point out that an FSLP with a Dirichlet boundary condition spectrum was investigated in papers [12,21] using variational methods. Here, we extend the study to the Prabhakar Sturm–Liouville problem and develop the results by transforming both differential fractional problems into the respective equivalent integral ones. Then, we analyse properties of the integral versions of fractional Sturm–Liouville operators (22) and (26) and apply the Hilbert–Schmidt spectral theorem to prove that their spectrum is purely discrete. Equivalence of the respective differential and integral fractional eigenvalue problems yields the theorems on spectra of the differential fractional and fractional Prabhakar eigenvalue problems given by Definitions 3 and 4. We begin our considerations with the case when q   =   0 .

3.1. Equivalence Results for Differential and Integral FSLP, PSLP: Case q   =   0

Here, we shall prove equivalence results for the FSLP/PSLP with an equation containing the fractional differential operators (22) and (26) and investigate the properties of the integral eigenvalue problem connected to the FSLE/PSLE in the case of order α fulfilling condition 1 α > 1 / 2 and solutions’ space restricted by the homogeneous Dirichlet boundary conditions.
In the first part, we transformed the differential fractional Sturm–Liouville problem (Definition 3) into the integral one on the subspace of the continuous functions defined below:
C D [ 0 , b ]   : =   { y C [ 0 , b ] ; y ( 0 )   =   y ( b )   =   0 } .
Let us note that the composition rules of fractional operators (7) and (9) allow us a to write a fractional Sturm–Liouville Equation (23) on the C D [ 0 , b ] space in the case of q   =   0 as follows:
L 0 1 λ I 0 + α 1 p I b α w ( x ) y ( x )   =   0
which leads to the integral equation
1 λ I 0 + α 1 p I b α w ( x ) y ( x )   =   C 1 w + C 2 w I 0 + α ( b x ) α 1 p ( x ) .
Constants C 1 w and C 2 w are determined by the homogeneous Dirichlet boundary conditions
C 1 w   =   0 , C 2 w   =   λ I 0 + α 1 p I b α w ( x ) y ( x ) | x = b I 0 + α ( b x ) α 1 p ( x ) | x = b .
The above calculations lead to the integral form of FSLE (23) with q   =   0
1 λ y ( x )   =   T w y ( x ) ,
where linear integral operator T w is built using the left and right Riemann–Liouville integrals and acts as follows:
T w y ( x )   =   I 0 + α 1 p I b α w ( x ) y ( x ) I 0 + α 1 p I b α w ( x ) y ( x ) | x = b I 0 + α ( b x ) α 1 p ( x ) | x = b · I 0 + α ( b x ) α 1 p ( x ) .
Similar considerations yield the integral form of PSLE (27) when q   =   0
1 λ y ( x )   =   T w y ( x ) ,
where linear integral operator T w is constructed using the left and right Prabhakar integrals and acts as follows
T w y ( x )   =   E ρ , γ , ω , 0 + α 1 p E ρ , γ , ω , b α w ( x ) y ( x )
E ρ , γ , ω , 0 + α 1 p E ρ , γ , ω , b α w ( x ) y ( x ) | x = b E ρ , γ , ω , 0 + α e ρ , α γ ( ω ( b x ) ρ ) p ( x ) | x = b · E ρ , γ , ω , 0 + α e ρ , α γ ( ω ( b x ) ρ ) p ( x ) .
We note that the above integral operators (35) and (37) can be rewritten as operators indexed by the arbitrary continuous function r (here, r   =   w ) and determined by the corresponding kernels— G 1 for FSLP and G 2 for PSLP:
T r y ( x )   : =   0 b G j ( x , s ) r ( s ) y ( s ) d x , j   =   1 , 2 ,
where kernels are of the form:
G 1 ( x , s )   : =   K 1 ( x , s ) K 1 ( b , x ) K 1 ( b , s ) K 1 ( b , b ) ,
G 2 ( x , s )   : =   K 1 P ( x , s ) K 1 P ( b , x ) K 1 P ( b , s ) K 1 P ( b , b ) ,
K 1 ( x , s )   =   0 min { x , s } ( x t ) α 1 Γ ( α ) · ( s t ) α 1 Γ ( α ) 1 p ( t ) d t ,
K 1 P ( x , s )   =   0 min { x , s } e ρ , α γ ( ω ( x t ) ρ ) e ρ , α γ ( ω ( s t ) ρ ) p ( t ) d t .
It is easy to check the following properties of kernels. First, they are symmetric functions on square Δ   =   [ 0 , b ] × [ 0 , b ]
K 1 ( x , s )   =   K 1 ( s , x ) , K 1 P ( x , s )   =   K 1 P ( s , x ) , G j ( x , s )   =   G j ( s , x )
and, in addition, we have
K 1 ( 0 , s )   =   K 1 ( b , 0 )   =   0 , K 1 P ( 0 , s )   =   K 1 P ( b , 0 )   =   0 , G j ( 0 , s )   =   G j ( b , s )   =   0 .
In our results developed in this paper, we apply two types of assumptions.
Hypothesis 1 (H1).
1 α > 1 / 2 , 1 p C [ 0 , b ] and function 1 p be positive on [ 0 , b ] or negative.
Hypothesis 2 (H2).
1 α > 1 / 2 , 1 p C [ 0 , b ] and function 1 p be positive on [ 0 , b ] or negative. In addition, let the real parameters α , ρ , γ , ω fulfil the conditions:
min { ρ , γ } > 0 , ω < 0 , α ρ γ , ρ < 1 .
Proposition 1.
If (H1) is fulfilled and function y L 2 ( 0 , b ) , then its image T r y C D [ 0 , b ] for any function r C [ 0 , b ] and operator defined by kernel (39).
If (H2) is fulfilled and function y L 2 ( 0 , b ) , then its image T r y C D [ 0 , b ] for any function r C [ 0 , b ] and operator defined by kernel (40).
Proof 
We sketch here the proof of the first part of the discussed proposition and omit the proof of the second one as it is analogous. By Corollary A1, kernel G 1 fulfills the Hölder condition; therefore, we find
| T r y ( x ) T r y ( x ) |   0 b | G 1 ( x , s ) G 1 ( x , s ) | · | r ( s ) y ( s ) | d s
M 1 | x x | β 0 b | r ( s ) y ( s ) | d s M 1 b · | | r | | · | | y | | L 2 · | x x | β
and we infer that image T r y is a continuous function and is even uniformly continuous on interval [ 0 , b ] .
We check that it obeys the homogeneous Dirichlet boundary conditions as well, because kernel G 1 fulfils the conditions (44):
T r y ( 0 )   =   0 b G 1 ( 0 , s ) y ( s ) r ( s ) d s   =   0 ,
T r y ( b )   =   0 b G 1 ( b , s ) y ( s ) r ( s ) d s   =   0 .
For functions belonging to the C D [ 0 , b ] space, we can prove the equivalence of the differential and integral form of the FSLP and PSLP, respectively. That is, the following two propositions are valid when q   =   0 . The first one concerns differential and integral fractional Sturm–Liouville problems.
Proposition 2.
If (H1) is fulfilled and w C [ 0 , b ] , then the following equivalence is valid on the C D [ 0 , b ] space
L 0 y ( x )   =   λ w ( x ) y ( x ) T w y ( x )   =   1 λ y ( x ) ,
where operator L 0 is defined in (22) and operator T w contains kernel (39).
Proof 
Assuming that y C D [ 0 , b ] is an eigenfunction corresponding to eigenvalue λ :
1 w ( x ) L 0 y ( x )   =   λ y ( x )
we act with the T w operator on both sides of this equation:
T w 1 w ( x ) L 0 y ( x )   =   λ T w y ( x )
and by applying composition rules (10), we obtain the integral eigenvalue equation
1 λ y ( x )   =   T w y ( x ) .
Next, we assume that function y L 2 ( 0 , b ) is an eigenfunction of the integral FSLP, i.e., Equation (47) is fulfilled. According to Proposition 1, eigenfunction y is a continuous one and belongs to the C D [ 0 , b ] space. Then, we calculate composition L 0 T w using the composition rules (7) and (9)
L 0 T w y ( x )   =   w ( x ) y ( x )
and by applying Equation (47), we arrive at the implication
L 0 T w y ( x )   =   w ( x ) y ( x )   =   1 λ L 0 y ( x ) L 0 y ( x )   =   λ w ( x ) y ( x ) .
Therefore, we conclude that on the C D [ 0 , b ] space, the equivalence of the differential and integral FSLP is valid. □
Below, we formulate the extended version of Proposition 2, where we describe the appropriate equivalence for Prabhakar Sturm–Liouville operators. Its proof is analogous to that presented above.
Proposition 3.
If (H2) is fulfilled and w C [ 0 , b ] , then the following equivalence
L 0 f ( x )   =   λ w ( x ) f ( x ) T w f ( x )   =   1 λ f ( x ) ,
is valid on the C D [ 0 , b ] space, where the L 0 operator is defined in (26) and the T w operator contains kernel (40).
Equivalence of the integral and differential fractional and fractional Prabhakar eigenvalue problems is an important step in deriving results on the spectrum for the problems described in Definitions 3 and 4. In the next section, we shall extend the equivalence results to the case where q 0 .

3.2. Equivalence Results for Differential and Integral FSLP, PSLP: General Case q 0

We begin our discussion with the fractional Sturm–Liouville problem. We write Equation (23) in the following form
1 w L q λ y ( x )   =   0
and apply composition rules for fractional operators (7) and (9)
1 w L 0 1 + I 0 + α 1 p I b α q ( x ) λ I 0 + α 1 p I b α w ( x ) y ( x )   =   0 .
The fractional differential Sturm–Liouville Equation (23) now takes the form of integral equation
y ( x ) + I 0 + α 1 p I b α q ( x ) + C 1 q + C 2 q I 0 + α ( b x ) α 1 p ( x )
= λ I 0 + α 1 p I b α w ( x ) y ( x ) + C 1 w + C 2 w I 0 + α ( b x ) α 1 p ( x ) ,
where constants are determined by the homogeneous Dirichlet boundary conditions; namely, C 1 w , C 2 w are given by (33) and for C 1 q , C 2 q , we have
C 1 q   =   0 , C 2 q   =   I 0 + α 1 p I b α q ( x ) y ( x ) | x   =   b I 0 + α ( b x ) α 1 p ( x ) | x = b .
To conclude, Equation (23) is now an integral equation
( 1 + T q ) y ( x )   =   λ T w y ( x ) ,
where the T w operator is given in (35) and the T q operator is given by the formula below
T q y ( x )   =   I 0 + α 1 p I b α q ( x ) y ( x ) I 0 + α 1 p I b α q ( x ) y ( x ) | x = b I 0 + α ( b x ) α 1 p ( x ) | x = b · I 0 + α ( b x ) α 1 p ( x ) .
Let us point out that, similar to the calculations presented in the previous part, both of the above integral operators can also be rewritten as integral operators (38) with kernel (39) for r   =   w and r   =   q , respectively.
Our aim is to reformulate the intermediate integral Equation (50) to the form of an eigenvalue equation. We apply Theorem A1 to invert the operator on the left-hand side. First, we check the assumption of Theorem A1, particularly when condition (H1) is fulfilled and w C [ 0 , b ] . We then apply Corollary A1, denoting K ( x , s )   =   G 1 ( x , s ) , and obtain:
G w ( · , s ) )   =   sup v [ 0 , b ] | G w ( v , s ) |
=   sup v [ 0 , b ] | G 1 ( v , s ) w ( s ) | w sup v [ 0 , b ] | G 1 ( v , s ) |
w sup v [ 0 , b ] | G 1 ( v , s ) G 1 ( 0 , s ) | + | G 1 ( 0 , s ) |
w · M 1 sup v [ 0 , b ] v α 1 / 2   =   w · M 1 · b α 1 / 2 < .
Next, we write condition (A8) in the explicit form:
ξ   =   sup x [ 0 , b ] 0 b | q ( s ) G 1 ( x , s ) | d s
=   sup x [ 0 , b ] 0 b | q ( s ) | · K 1 ( x , s ) K 1 ( b , x ) K 1 ( b , s ) K 1 ( b , b ) d s < 1 .
All the above considerations lead to the proposition on convergence of the series associated with the intermediate fractional integral eigenvalue problem given in (50) and (A5). Analogous convolutions’ series were also studied on the C [ a , b ] and L 2 ( a , b ) function spaces for FSLPs with homogeneous mixed and Robin boundary conditions, respectively [13,14].
Proposition 4.
Let (H1) be fulfilled, w , q C [ 0 , b ] and function w be positive. If condition (52) is fulfilled, then for any function y L 2 ( 0 , b ) series on the right-hand side of the formula below is uniformly convergent on interval [ 0 , b ] :
T y ( x )   : =   ( 1 + T q ) 1 T w y ( x )   =   T w y ( x ) + n   =   1 ( T q ) n T w y ( x ) ,
where operators T q , T w are defined in (A6) and (A7) with K ( x , s )   =   G 1 ( x , s ) . In addition, series (A9) determining the kernel of integral operator T in (53) is uniformly convergent on square Δ and kernel G is continuous on Δ.
Proof 
Let us observe that the composition of operators T q T w is an integral operator
T q T w y ( x )   =   0 b d s G q ( x , s ) 0 b G w ( s , u ) y ( u ) d u
=   0 b d u y ( u ) 0 b G q ( x , s ) G w ( s , u ) d s   =   0 b G q * G w ( x , u ) y ( u ) d u ,
where the kernel is defined by the following convolution:
A * B ( x , u )   : =   0 b A ( x , s ) B ( s , u ) d s .
We shall prove that the compositions ( T q ) n T w are also defined by convolutions of kernels G q and G w . We start with the induction hypothesis:
( T q ) n T w y ( x )   =   0 b ( G q * n ) * G w ( x , u ) y ( u ) d u
and we prove that this formula is valid for the next step n + 1 as well:
( T q ) n + 1 T w y ( x )   =   0 b ( G q * ( n + 1 ) ) * G w ( x , u ) y ( u ) d u .
We begin with the left-hand side, applying the induction hypothesis and associativity property of the convolutions of continuous functions:
( T q ) n + 1 T w y ( x )   =   0 b d s G q ( x , s ) ( T q ) n T w y ( s )
=   0 b d s G q ( x , s ) 0 b G q * n * G w ( s , u ) y ( u ) d u
=   0 b d u y ( u ) G q * G q * n * G w ( x , u )
=   0 b G q * ( n + 1 ) * G w ( x , u ) y ( u ) d u .
As inductive hypothesis (54) leads to the validity of the next step n + 1 ; we infer that formula (54) holds for any natural number n 1 .
Now, we apply Theorem A1 and calculate kernel G for integral operator T   : =   ( 1 + T q ) 1 T w :
T y ( x )   =   T w y ( x ) + n   =   1 ( T q ) n T w y ( x )
=   0 b G w ( x , s ) y ( s ) d s + n   =   1 ( 1 ) n 0 b G q * n * G w ( x , s ) y ( s ) d s
=   0 b G w ( x , s ) + n   =   1 ( 1 ) n G q * n * G w ( x , s ) y ( s ) d s   =   0 b G ( x , s ) y ( s ) d s .
The above calculations lead to the thesis of Proposition 4; namely, operator T, defined by series (53), is correctly defined on space L w 2 ( 0 , b )   =   L 2 ( 0 , b ) as an integral operator with a continuous kernel G:
T y ( x )   =   0 b G ( x , s ) y ( s ) d s .
Having constructed operator T, we now prove the equivalence result, connecting the differential and integral fractional Sturm–Liouville problems in the general case.
Proposition 5.
If (H1) and condition (52) are fulfilled, w , q C [ 0 , b ] and function w is positive, then the following equivalence is valid on the C D [ 0 , b ] space
L q y ( x )   =   λ w ( x ) y ( x ) T y ( x )   =   1 λ y ( x ) ,
where the L q operator is defined in (22) and the T operator is given in (53) with a kernel determined by series (A9) with K ( x , s )   =   G 1 ( x , s ) .
Proof 
We recall that for any function y C D [ 0 , b ] , we have (proof of Proposition 2)
T w 1 w ( x ) L 0 y ( x )   =   y ( x ) ,
and we extend this equality to the analogous formula for operators T and L q
T 1 w ( x ) L q y ( x )   =   T 1 w ( x ) L 0 y ( x ) + T q ( x ) w ( x ) y ( x )
=   T w + n   =   1 ( T q ) n T w 1 w ( x ) L 0 y ( x ) + T q y ( x ) + n   =   1 ( T q ) n T q y ( x )
=   y ( x ) + n   =   1 ( T q ) n y ( x ) + T q y ( x ) + n   =   1 ( T q ) n T q y ( x )
=   y ( x ) ,
where we calculate the corresponding formulas for series by using the fact that operator T is a uniformly convergent series (Proposition 4) when acting on the C D [ 0 , b ] space. For differential FSLE,
1 w ( x ) L q y ( x )   =   λ y ( x )
after calculating the image of the T operator of functions on both sides of FSLE
T 1 w ( x ) L q y ( x )   =   y ( x )   =   λ T y ( x ) ,
we obtain the integral fractional Sturm–Liouville equation in the form of
T y ( x )   =   1 λ y ( x ) .
In the next step, we assume that the above integral FSLE is fulfilled. Then, function y C D [ 0 , b ] . We apply the differential operator L q to both sides of the integral FSLE
L q T y ( x )   =   1 λ L q y ( x ) .
For the composition of operators on the left-hand side, we get for continuous functions f , y C D [ 0 , b ]
L 0 T w f ( x )   =   w ( x ) f ( x ) , L 0 T q f ( x )   =   q ( x ) f ( x ) ,
L 0 ( T q ) n T w y ( x )   =   q ( x ) ( T q ) n 1 T w y ( x ) .
Applying Proposition 4 again, we obtain the following result for the composition of the L q and T operators
L q T y ( x )   =   ( q ( x ) + L 0 ) T w y ( x ) + n   =   1 ( T q ) n T w y ( x )
=   q ( x ) T w y ( x ) + q ( x ) n   =   1 ( T q ) n T w y ( x ) + w ( x ) y ( x ) q ( x ) n   =   1 ( T q ) n 1 T w y ( x )
=   w ( x ) y ( x ) .
From this relation, we derive the differential fractional eigenvalue equation
w ( x ) y ( x )   =   1 λ L q y ( x )
which leads to the differential fractional Sturm–Liouville equation:
L q y ( x )   =   λ w ( x ) y ( x )
and this ends the proof of equivalence (55). □
Now, we generalize the Sturm–Liouville operator L q by introducing Prabhakar derivatives and we move on to the Prabhakar Sturm–Liouville problem (PSLP) determined in Definitions 2 and 4 and discussed in [26] in the case when the solutions’ space is restricted by the mixed homogeneous boundary conditions.
1 w L q λ y ( x )   =   0
We obtain the intermediate form of the integral fractional Prabhakar eigenvalue equation applying composition rules (18)–(21)
( 1 + T q ) y ( x )   =   T w y ( x ) ,
where integral operator T w is given in Formula (37) and operator T q looks as follows:
T q y ( x )   =   E ρ , γ , ω , 0 + α 1 p E ρ , γ , ω , b α q ( x ) y ( x )
E ρ , γ , ω , 0 + α 1 p E ρ , γ , ω , b α q ( x ) y ( x ) | x   =   b E ρ , γ , ω , 0 + α e ρ , α γ ( ω ( b x ) ρ ) p ( x ) | x = b · E ρ , γ , ω , 0 + α e ρ , α γ ( ω ( b x ) ρ ) p ( x ) .
Similar to the previous calculations for FSLP, operators (37) and (57) can be rewritten as integral operators (38), with kernel G 2 given in (40) for r   =   w and r   =   q , respectively. Again, we apply Theorem A1 to invert operator 1 + T q . First, we check the assumption of Theorem A1, assuming that (H2) is fulfilled and applying Corollary A1 with K ( x , s )   =   G 2 ( x , s ) :
G w ( · , s ) )   =   sup v [ 0 , b ] | G w ( v , s ) |
=   sup v [ 0 , b ] | G 2 ( v , s ) w ( s ) | w sup v [ 0 , b ] | G 2 ( v , s ) |
w sup v [ 0 , b ] | G 2 ( v , s ) G 2 ( 0 , s ) | + | G 2 ( 0 , s ) |
w · M 2 sup v [ 0 , b ] | v | β   =   w · M 2 · b β < .
Next, we write condition (A8) in the explicit form:
ξ   =   sup x [ 0 , b ] 0 b | q ( s ) G 2 ( x , s ) | d s
=   sup x [ 0 , b ] 0 b | q ( s ) | · K 1 P ( x , s ) K 1 P ( b , x ) K 1 P ( b , s ) K 1 P ( b , b ) d s < 1 .
In the proposition below, we describe the inverse operator ( 1 + T q ) 1 connected to the intermediate Equation (56). We omit the proof as it is a straightforward corollary of Theorem A1, and the full proof is analogous to that of Proposition 4.
Proposition 6.
Let (H2) be fulfilled, w , q C [ 0 , b ] and function w be positive. If condition (58) is fulfilled, then for any function y L 2 ( 0 , b ) series on the right-hand side of the formula below is uniformly convergent on interval [ 0 , b ] :
T y ( x )   : =   ( 1 + T q ) 1 T w y ( x )   =   T w y ( x ) + n   =   1 ( T q ) n T w y ( x ) ,
where operators T q , T w are defined in (A6) and (A7) with K ( x , s )   =   G 2 ( x , s ) . In addition, series (A9) determining kernel of integral operator T in (59) is uniformly convergent on square Δ and kernel G is continuous on Δ.
Similar to Proposition 5, we formulate the equivalence result for integral and differential version of eigenvalue equations corresponding to PSLP. The proof is based on the composition rules (18) and (19) and on Proposition 6, which describes inverse integral operator ( 1 + T q ) 1 . We omit the proof as it is analogous to the proof of Proposition 5.
Proposition 7.
If (H2) and condition (58) are fulfilled, w , q C [ 0 , b ] and function w is positive, then the following equivalence is valid on the C D [ 0 , b ] space
L q y ( x )   =   λ w ( x ) y ( x ) T y ( x )   =   1 λ y ( x ) ,
where the L q operator is defined in (26), T operator is given in (59) with the kernel determined by the series (A9) and K ( x , s )   =   G 2 ( x , s ) .

4. Results on the Spectrum of Integral and Differential Fractional and Fractional Prabhakar Sturm–Liouville Problems

In the previous section, we discussed and proved the results on the equivalence of differential and integral forms of fractional eigenvalue problems. First, Propositions 2 and 3 describe the equivalence for fractional and fractional Prabhakar Sturm–Liouville problems when fractional differential operators are respectively L 0 and L 0 , i.e., q   =   0 . In this case, the corresponding integral operators are T w with kernels G 1 and G 2 . We prove the spectral results for these operators by applying the Hilbert–Schmidt theorem.

4.1. Case: q   =   0

Theorem 1.
If (H1) is fulfilled and w C [ 0 , b ] is a positive function, then the spectrum of operator T w defined by (38) and (39) is a discrete one, enclosed in the interval ( 1 , 1 ) , with 0 being its only limit point. Eigenfunctions belong to the C D [ 0 , b ] space and form an orthogonal basis in the L w 2 ( 0 , b ) space.
If (H2) is fulfilled and w C [ 0 , b ] is a positive function, then the spectrum of operator T w defined by (38) and (39) is a discrete one, enclosed in the interval ( 1 , 1 ) , with 0 being its only limit point. Eigenfunctions belong to the C D [ 0 , b ] space and form an orthogonal basis in the L w 2 ( 0 , b ) space.
Proof 
Let us observe that when weight function fulfils the assumptions of the theorem, we have for functions spaces
L 2 ( 0 , b )   =   L w 2 ( 0 , b ) , L 2 ( Δ )   =   L w w 2 ( Δ ) .
The integral Hilbert–Schmidt operator T w , defined by kernel G 1 , is a compact one, as this kernel is a function continuous on square Δ and G 1 L w w 2 ( Δ ) .
It is also a self-adjoint operator on L w 2 ( 0 , b ) , because kernel G 1 is a symmetric function on square Δ , and for an arbitrary pair of functions f , g L w 2 ( 0 , b ) , we obtain:
g , T w f w   =   0 b d x w ( x ) g ( x ) ¯ 0 b G 1 ( x , s ) f ( s ) w ( s ) d s
=   0 b d s w ( s ) f ( s ) 0 b G 1 ( s , x ) g ( x ) ¯ w ( x ) d x
=   f , T w g w ¯   =   T w g , f w .
The thesis is a straightforward result of the Hilbert–Schmidt spectral theorem. We omit the proof of the second part as it is analogous to the one presented above. □
The spectral theorem for integral fractional and Prabhakar Sturm–Liouville operators together with the equivalence results, included in Propositions 2 and 3, lead to the theorem on the spectrum of differential fractional eigenvalue problems subjected to the homogeneous Dirichlet boundary conditions in the case when q   =   0 .
Theorem 2.
If (H1) is fulfilled and w C [ 0 , b ] is a positive function, then the spectrum of operator L 0 defined by (22) and considered on the C D [ 0 , b ] space is a discrete one, and | λ n | . Eigenfunctions belonging to the C D [ 0 , b ] space form an orthogonal basis in the L w 2 ( 0 , b ) space.
If (H2) is fulfilled and w C [ 0 , b ] is a positive function, then the spectrum of operator L 0 , defined by (26) and considered on the C D [ 0 , b ] space is a discrete one and | λ n | . Eigenfunctions belonging to the C D [ 0 , b ] space form an orthogonal basis in the L w 2 ( 0 , b ) space.

4.2. General Case q 0

We observe that the analogous equivalence of differential and integral FSLP holds in the general case q 0 as well. This result is given by Proposition 5. Analogously, Proposition 7 gives the equivalence relation of both versions of the fractional Prabhakar Sturm–Liouville problem. The results, included in the mentioned propositions, allow us to rewrite eigenvalue equations, replacing the differential FSLO and PSLO with the corresponding integral operators T. These operators, first determined as operator series with convergence described in Propositions 4 and 6, are in fact integral Hilbert–Schmidt operators. Their kernels—sums of a uniformly convergent series of convolutions—are continuous functions on square Δ . The theorem below describes the spectrum of fractional integral operators T with kernel G, determined by kernels G 1 and G 2 , respectively.
Theorem 3.
If (H1) and condition (52) are fulfilled, w , q C [ 0 , b ] and w is a positive function; then the spectrum of operator T defined by (53) with kernel G given in (A9) with K ( x , s )   =   G 1 ( x , s ) is a discrete one, enclosed in interval ( 1 , 1 ) , with 0 being its only limit point. Eigenfunctions belong to the C D [ 0 , b ] space and form an orthogonal basis in the L w 2 ( 0 , b ) space.
If (H2) and condition (58) are fulfilled, w , q C [ 0 , b ] and w is a positive function, then the spectrum of operator T is defined by (59), with kernel G given in (A9) and with K ( x , s )   =   G 2 ( x , s ) is a discrete one, enclosed in interval ( 1 , 1 ) , with 0 being its only limit point. Eigenfunctions belong to the C D [ 0 , b ] space and form an orthogonal basis in the L w 2 ( 0 , b ) space.
Proof 
Let us again observe that when the weight function fulfils assumptions of the theorem; we have for spaces considered as sets of functions
L 2 ( 0 , b )   =   L w 2 ( 0 , b ) , L 2 ( Δ )   =   L w w 2 ( Δ ) .
Integral Hilbert–Schmidt operator T, defined by kernel G, is a compact one as this kernel is a continuous function on square Δ and G L w w 2 ( Δ ) .
We recall (proof of Theorem 1) that on the L w 2 ( 0 , b ) space, the following equality holds for the arbitrary pair of functions f , g L w 2 ( 0 , b ) :
g , T w f w   =   T w g , f w
because kernel G 1 is a symmetric function on square Δ . Next, for the composition of operators T q T w , we obtain the relation
g , T q T w f w   =   g , T w q w T w f w   =   q w T w g , T w f w   =   T q T w g , f w .
Now, we apply the mathematical induction principle to prove that such relations hold for arbitrary n > 1 natural. We formulate an induction hypothesis in the form of
g , ( T q ) n T w f w   =   ( T q ) n T w g , f w
and for step n + 1 , we achieve
g , ( T q ) n + 1 T w f w   =   g , T q ( T q ) n T w f w   =   q w T w g , ( T q ) n T w f w
=   ( T q ) n T w q w T w g , f w   =   ( T q ) n + 1 T w g , f w .
Applying the mathematical induction principle, we infer that Formula (61) is valid for all natural numbers n 1 . We use this formula in the proof of the fact that integral operator T is a self-adjoint one. Remembering that it is represented by a series, uniformly convergent on the Hilbert space (Proposition 4), we calculate the scalar product term by term
g , T f w   =   g , T w + n   =   1 ( T q ) n T w f w
=   T w g , f w + n   =   1 ( 1 ) n ( T q ) n T w g , f w
=   T w + n   =   1 ( T q ) n T w g , f w
=   T g , f w .
To conclude, the integral operator T with a kernel G given in (A9) with K ( x , s )   =   G 1 ( x , s ) is a compact and self-adjoint operator on Hilbert space L w 2 ( 0 , b ) . Therefore, the thesis of the first part of the theorem holds by the Hilbert–Schmidt spectral theorem.
Proof of the second part for operator T, associated with the integral PSLP with homogeneous Dirichlet boundary conditions, is analogous. □
Now, we apply the above spectral theorem for integral fractional eigenvalue problems, with equivalence results enclosed in Propositions 5 and 7 to formulate a theorem on discrete spectra for differential fractional and fractional Prabhakar Sturm–Liouville problems.
Theorem 4.
If (H1) and condition (52 ) are fulfilled, w , q C [ 0 , b ] and w is a positive function, then the spectrum of operator L q defined by (22) and considered on the C D [ 0 , b ] space is a discrete one, and | λ n | . Eigenfunctions, belonging to the C D [ 0 , b ] space, form an orthogonal basis in the L w 2 ( 0 , b ) space.
If (H2) and condition (58) are fulfilled, w , q C [ 0 , b ] and w is a positive function; then the spectrum of operator L q defined by (26) and considered on the C D [ 0 , b ] space, is a discrete one and | λ n | . Eigenfunctions, belonging to the C D [ 0 , b ] space, form an orthogonal basis in the L w 2 ( 0 , b ) space.

5. Discussion

In this paper, we presented results on the discrete spectrum of fractional and fractional Prabhakar Sturm–Liouville problems in a case when eigenfunctions’ space is subjected to the homogenous Dirichlet boundary conditions. First, we extended the idea of the fractional to the fractional Prabhakar eigenvalue problem, where the Sturm–Liouville operator was constructed by using the left and right Prabhakar derivatives.
Prabhakar derivatives, with respect to time, were recently applied in anomalous diffusion models [27,28]. The derived spectral results for regular PSLP with Dirichlet boundary conditions will be used in developing equations with fractional partial derivatives with respect to the space–variable.
It appears that the method of converting the differential eigenvalue problem into the equivalent integral one can be applied to both types of Sturm–Liouville operator. This approach, developed in [13,14] for fractional eigenvalue problems subject to the homogeneous mixed and Robin boundary conditions, is extended to the case of FSLP with Dirichlet boundary conditions and generalized to PSLP with the same type of conditions.
Let us point out that the spectrum and eigenfunctions of fractional eigenvalue problems with Dirichlet boundary conditions were also studied in [11,16] by applying variational methods. The first of these papers describes the spectrum of FSLP for a fractional order in the range ( 1 / 2 , 1 ] , and the spectral result was extended to range ( 0 , 1 / 2 ] in [16]. Comparing both of the methods—the variational one and the transformation into integral FSLP/PSLP—we observe that in the case of Dirichlet boundary conditions, the range of order is wider in the variational method. Nevertheless, the approach proposed here has an advantage of providing the spectral results for regular PSLP as well. Simultaneously, we obtain eigenfunctions’ systems for both types of eigenvalue problems, which provide orthogonal bases in the corresponding Hilbert spaces. Such bases are a meaningful tool in applications in constructing and solving partial differential fractional equations, for example, space-fractional diffusion equations in the finite domain, as well as fractional equations governing control systems (compare references and examples in [29]).

6. Conclusions

The results developed in this paper describe the spectrum and eigenfunctions properties for FSLP and PSLP subjected to homogeneous Dirichlet boundary conditions. It seems that the conversion method can also be easily applied to other Prabhakar Sturm–Liouville problems; in particular, we shall construct the corresponding mixed, Robin, and Neumann boundary conditions and develop the equivalence results. Then, we will construct the integral PSLO with kernels analogous to those from the papers [13,14] and study the spectral properties, both for the integral and differential PSLPs.
Regarding the extension of the range of fractional order for the conversion method, we observe that so far we proved equivalence results on the space of continuous solutions. This restriction is connected to the version of Hölder condition for kernels, as discussed in Lemma A1 and Corollary A1. Thus, the aim of our future work will be to weaken this condition and to extend the range of fractional order.
Further, our investigations will include numerical simulations in order to derive approximate values of eigenvalues and eigenfunctions. As was shown in the papers [13,14], the integral form of the fractional Sturm–Liouville eigenvalue equation is particularly useful as a first step of the numerical method of solving FSLP. Thus, our aim will be to discretize integral eigenvalue problems and apply the equivalence results, enclosed in Propositions 2 and 3 for the case q   =   0 , and in Propositions 5 and 7, when q 0 . In this way, we shall arrive at numerical solutions of differential FSLP and PSLP with Dirichlet boundary conditions.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Appendix A.1

Let us point out that the three-parameter Mittag–Leffler function (11), which appears in the definition of the Prabhakar function (12), is a completely monotone function [30], and this property leads to the following two inequalities. First, when parameters α , ρ , γ , ω are real and obey conditions
α ( 0 , 1 ] , min { ρ , γ } > 0 , ω < 0 , α ρ γ , ρ 1 ,
the three-parameter Mittag–Leffler function is bounded on any interval [ 0 , b ] ( M e is a constant)
| E ρ , α γ ( ω x ρ ) | M e
and it fulfills the Lipschitz condition on this interval ( M L is a constant)
| E ρ , α γ ( ω ( x ) ρ ) E ρ , α γ ( ω x ρ ) | M L | ( x ) ρ x ρ | .
In addition, we remember that the power function obeys the Hölder condition on interval [ 0 , 1 ] when ρ 1 ( M ρ is a constant)
| ( y ) ρ y ρ | M ρ | y y | ρ , y , y [ 0 , 1 ] .
All the above inequalities will be applied to derive properties of a fractional integral operator associated to the differential Prabhakar Sturm–Liouville operator (PSLO). In particular, they are important in the study of Hölder continuity and the continuity of kernels of integral versions of Prabhakar Sturm–Liouville operators. The lemma below summarizes the Hölder continuity properties of kernels K 1 , K 1 P and was proven in [26] (compare Properties 3.2 and 3.3).
Lemma A1.
If (H1) is fulfilled, then kernel K 1 , given by (41), obeys the Hölder-type condition, i.e., there exists coefficient β ( 0 , 1 ] and function m L 2 ( 0 , b ) such that
| K 1 ( x , s ) K 1 ( x , s ) | m ( s ) | x x | β ,
where β   =   α 1 / 2 and
m ( s )   =   2 b α 1 / 2 1 / p ( Γ ( α ) ) 2 ( α 1 / 2 )
is a constant function.
If (H2) is fulfilled, then kernel K 1 P , given by (42), obeys the Hölder-type condition, i.e., there exists coefficient β ( 0 , 1 ] and function m L 2 ( 0 , b ) such that
| K 1 P ( x , s ) K 1 P ( x , s ) | m ( s ) | x x | β ,
where β   =   min { α 1 / 2 , ρ } and
m ( s )   =   max { b α 1 / 2 , b 2 α 1 ρ } · 1 / p · M e α 1 / 2 · ( 2 M e + M L M ρ b ρ )
is a constant function.
Analyzing the construction of kernels G 1 and G 2 , we obtain the following corollary.
Corollary A1.
If (H1) is fulfilled, then kernel G 1 , defined by Formulas (39) and (41), obeys the Hölder-type condition, i.e., there exists coefficient β ( 0 , 1 ] and constant M 1 such that
| G 1 ( x , s ) G 1 ( x , s ) | M 1 | x x | β ,
where β   =   α 1 / 2 and
M 1   =   2 b α 1 / 2 | | 1 / p 1 + 1 / p · p ( Γ ( α ) ) 2 ( α 1 / 2 ) .
If (H2) is fulfilled, then kernel G 2 , defined by Formulas (40) and (42), obeys the Hölder-type condition, i.e., there exists coefficient β ( 0 , 1 ] and constant M 2 such that
| G 2 ( x , s ) G 2 ( x , s ) | M 2 | x x | β ,
where β   =   min { α 1 / 2 , ρ } and
M 2   =   max { b α 1 / 2 , b 2 α 1 ρ } · 1 / p 1 + 1 / p · p · M e α 1 / 2 · ( 2 M e + M L M ρ b ρ ) .
Proof 
We prove the Hölder-type condition for kernel G 1 by applying Lemma A1 and the symmetry property of kernel K 1 given in (43). We begin by estimating values K 1 ( b , b ) and K 1 ( b , s ) :
| K 1 ( b , b ) | = a b ( b t ) 2 α 2 ( Γ ( α ) ) 2 p ( t ) d t   =   a b ( b t ) 2 α 2 ( Γ ( α ) ) 2 | p ( t ) | d t
( b a ) 2 α 1 ( Γ ( α ) ) 2 ( 2 α 1 ) p ,
| K 1 ( b , s ) | = a s ( b t ) α 1 Γ ( α ) · ( s t ) α 1 Γ ( α ) 1 p ( t ) d t
1 p · ( b a ) 2 α 1 ( Γ ( α ) ) 2 ( 2 α 1 ) .
Now, we apply the derived inequalities and condition (A1)
| G 1 ( x , s ) G 1 ( x , s ) | =
=   K 1 ( x , s ) K 1 ( b , x ) K 1 ( b , s ) K 1 ( b , b ) K 1 ( x , s ) + K 1 ( b , x ) K 1 ( b , s ) K 1 ( b , b )
| K 1 ( x , s ) K 1 ( x , s ) | + | K 1 ( x , b ) K 1 ( x , b ) | · K 1 ( b , s ) K 1 ( b , b )
| x x | β m ( s ) 1 + K 1 ( b , s ) K 1 ( b , b )
m ( s ) 1 + 1 / p · p | x x | β
=   M 1 | x x | β ,
where
M 1   =   2 b α 1 / 2 1 / p 1 + 1 / p · p ( Γ ( α ) ) 2 ( α 1 / 2 ) .
The proof of the Hölder condition for kernel G 2 is analogous. □
The next corollary results from the Hölder conditions (A3) and (A4) and symmetry properties of kernels G 1 , G 2 given in (43) and yields continuity of both kernels on square Δ   =   [ 0 , b ] × [ 0 , b ] .
Corollary A2.
If (H1) is fulfilled, then kernel G 1 , defined by Formulas (39) and (41), is continuous on square Δ   =   [ 0 , b ] × [ 0 , b ] .
If (H2) is fulfilled, then kernel G 2 , defined by Formulas (40) and (42) is continuous on square Δ   =   [ 0 , b ] × [ 0 , b ] .
Proof 
Let us note that the symmetry of kernel G 1 allows us to write condition (A3) in the following form
| G 1 ( x , s ) G 1 ( x , s ) | M 1 | x x | β + | s s | β .
To prove continuity of the kernel, we apply the Cauchy definition of continuous function, i.e., we take arbitrary ϵ > 0 and assume that the distance between points ( x , s ) , ( x , s ) is smaller than δ ( ϵ )   =   ϵ 2 M 1 1 / β , which means
| ( x , s ) ( x , s ) |   =   ( x x ) 2 + ( s s ) 2 < δ ( ϵ ) .
We observe that the following inequalities are then valid
| x x | < δ ( ϵ ) , | s s | < δ ( ϵ ) .
Now, we check the distance between the values of function G 1 :
| G 1 ( x , s ) G 1 ( x , s ) | M 1 | x x | β + | s s | β
M 1 · 2 ( δ ( ϵ ) ) β   =   ϵ .
We see that for arbitrary ϵ > 0 , bound δ ( ϵ ) for the distance of points exists, such that the implication below is valid
| ( x , s ) ( x , s ) | < δ ( ϵ ) | G 1 ( x , s ) G 1 ( x , s ) | < ϵ .
Thus, kernel G 1 is a continuous function on square Δ by the Cauchy definition of continuity.
Proof for kernel G 2 is analogous. □

Appendix A.2

We shall study properties of integral equations of the form:
( 1 + T q ) y ( x )   =   λ T w y ( x )
determined on the L w 2 ( a , b ) function space. Such an equation is the intermediate stage of transformation of the fractional differential eigenvalue problems into the equivalent integral ones (see examples in papers [13,14]). In cases where the integral operator on the left-hand side of (A5) is invertible, we can convert fractional differential Sturm–Liouville operator into an integral one. Then, we can study spectral properties of the integral operator and derive results for the spectrum and eigenfunctions of the fractional differential Sturm–Liouville problems connected to various homogeneous boundary conditions.
Operators T q and T w are integral ones, with kernels given in the form of
T q y ( x )   : =   a b G q ( x , s ) y ( s ) d s , G q ( x , s )   =   K ( x , s ) q ( s ) ,
T w y ( x )   : =   a b G w ( x , s ) y ( s ) d s , G w ( x , s )   =   K ( x , s ) w ( s ) .
We formulate below a theorem which we shall apply to analyse integral eigenvalue problems associated with the fractional differential ones.
Theorem A1.
Let function q C [ a , b ] and function | | G w ( · , s ) | | : =   sup v [ a , b ] | G w ( v , s ) | be bounded on interval [ a , b ] . If condition
ξ   : =   sup x [ a , b ] a b | G q ( x , v ) | d v < 1
is fulfilled, then the series
G ( x , s )   : =   G w ( x , s ) + n   =   1 ( 1 ) n G q * n * G w ( x , s )
is uniformly convergent on square Δ; i.e., the sum of this series G is determined for all points ( x , s ) Δ .
If, in addition, kernels G q , G w C ( Δ ) , then sum G C ( Δ ) .
Proof 
We shall apply the mathematical induction principle to estimate all terms of series (A9). First, we estimate the absolute value of the first convolution term:
| G q * G w ( x , s ) |   =   a b G q ( x , v ) G w ( v , s ) d v
a b | G q ( x , v ) G w ( v , s ) | d v G w ( · , s ) · sup x [ a , b ] a b | G q ( x , v ) | d v   =   ξ · | | G w ( · , s ) | |
and for the second term, we obtain
| G q * G q * G w ( x , s ) | ξ sup v [ a , b ] | G q * G w ( v , s ) |
ξ 2 · | | G w ( · , s ) | | .
Now, we formulate the induction hypothesis (here, n > 2 is a natural number):
| G q * n * G w ( x , s ) | ξ n · | | G w ( · , s ) | |
and we shall prove that it holds for the next step n + 1
| G q * ( n + 1 ) * G w ( x , s ) | ξ n + 1 · | | G w ( · , s ) | | .
We begin from the left-hand side of the above inequality and we find
| G q * ( n + 1 ) * G w ( x , s ) |   =   | G q * G q * n * G w ( x , s ) |
ξ sup v [ a , b ] | G q * n * G w ( v , s ) |
ξ n + 1 sup v [ a , b ] | G w ( v , s ) | ξ n + 1 · | | G w ( · , s ) | | .
The induction hypothesis (A12) implies the validity of the next step for n + 1 ; therefore, we infer that estimation (A12) is valid for all terms indexed by n 1 . Now, we are ready to consider the convergence of the function series (A9) by using the Weierstrass convergence test and inequality (A12). We observe that the majorant number series (a geometric one) is absolutely convergent under the assumption (A8). Thence, the function series (A9) is absolutely and uniformly convergent, as we achieve for any point ( x , s ) Δ
G w ( x , s ) + n   =   1 ( 1 ) n G q * n * G w ( x , s )
| G w ( x , s ) | + n   =   1 ξ n | | G w ( · , s ) | |   =   | G w ( x , s ) | + | | G w ( · , s ) | | · ξ 1 ξ .
In the second part of Theorem 2, we note that continuity of kernels G q , G w implies that all terms of the series (A9) are continuous as convolutions of continuous functions. The absolutely and uniformly convergent series (A9) leads to sum G, which is also continuous on Δ . □

References

  1. Al-Mdallal, Q.M. An efficient method for solving fractional Sturm-Liouville problems. Chaos Solitons Fractals 2009, 40, 183–189. [Google Scholar] [CrossRef]
  2. Al-Mdallal, Q.M. On the numerical solution of fractional Sturm-Liouville problems. Int. J. Comput. Math. 2010, 87, 2837–2845. [Google Scholar] [CrossRef]
  3. Erturk, V.S. Computing eigenelements of Sturm-Liouville Problems of fractional order via fractional differential transform method. Math. Comput. Appl. 2011, 16, 712–720. [Google Scholar] [CrossRef]
  4. Duan, J.-S.; Wang, Z.; Liu, Y.-L.; Qiu, X. Eigenvalue problems for fractional ordinary differential equations. Chaos Solitons Fractals 2013, 46, 46–53. [Google Scholar] [CrossRef]
  5. Qi, J.; Chen, S. Eigenvalue problems of the model from nonlocal continuum mechanics. J. Math. Phys. 2011, 52, 073516. [Google Scholar] [CrossRef]
  6. Atanackovic, T.M.; Stankovic, S. Generalized wave equation in nonlocal elasticity. Acta Mech. 2009, 208, 1–10. [Google Scholar] [CrossRef]
  7. Klimek, M.; Agrawal, O.P. On a Regular Fractional Sturm-Liouville Problem with Derivatives of Order in (0,1). In Proceedings of the 13th International Carpathian Control Conference, (ICCC 20212), High Tatras, Slovakia, 28–21 May 201; pp. 284–289. [CrossRef]
  8. d’Ovidio, M. From Sturm-Liouville problems to to fractional and anomalous diffusions. Stoch. Process. Appl. 2012, 122, 3513–3544. [Google Scholar] [CrossRef] [Green Version]
  9. Rivero, M.; Trujillo, J.J.; Velasco, M.P. A fractional approach to the Sturm-Liouville problem. Cent. Eur. J. Phys. 2013, 11, 1246–1254. [Google Scholar] [CrossRef]
  10. Zayernouri, M.; Karniadakis, G.E. Fractional Sturm–Liouville eigen-problems: Theory and numerical approximation. J. Comput. Phys. 2013, 252, 495–517. [Google Scholar] [CrossRef]
  11. Klimek, M.; Agrawal, O.P. Fractional Sturm-Liouville Problem. Comput. Math. Appl. 2013, 66, 795–812. [Google Scholar] [CrossRef]
  12. Klimek, M.; Odzijewicz, T.; Malinowska, A.B. Variational methods for the fractional Sturm-Liouville problem. J. Math. Anal. Appl. 2014, 416, 402–426. [Google Scholar] [CrossRef]
  13. Klimek, M.; Ciesielski, M.; Blaszczyk, T. Exact and numerical solutions of the fractional Sturm-Liouville problem. Fract. Calc. Appl. Anal. 2018, 21, 45–71. [Google Scholar] [CrossRef] [Green Version]
  14. Klimek, M. Homogeneous Robin boundary conditions and discrete spectrum of fractional eigenvalue problem. Fract. Calc. Appl. Anal. 2019, 22, 78–94. [Google Scholar] [CrossRef]
  15. Zayernouri, M.; Ainsworth, M.; Karniadakis, G.E. Tempered fractional Sturm–Liouville eigenproblems. SIAM J. Sci. Comput. 2015, 37, A1777–A1800. [Google Scholar] [CrossRef] [Green Version]
  16. Pandey, P.; Pandey, R.; Yadav, S.; Agrawal, O.P. Variational Approach for Tempered Fractional Sturm–Liouville Problem. Int. J. Appl. Comput. Math. 2021, 7, 51. [Google Scholar] [CrossRef]
  17. Al-Refai, M.; Abdeljawad, T. Fundamental results of conformable Sturm-Liouville eigenvalue problems. Complexity 2017, 2017, 3720471. [Google Scholar] [CrossRef] [Green Version]
  18. Mortazaasl, H.; Jodayree Akbarfam, A.A. Two classes of conformable fractional Sturm-Liouville problems: Theory and applications. Math. Meth. Appl. Sci. 2021, 44, 166–195. [Google Scholar] [CrossRef]
  19. Li, J.; Qi, J. On a nonlocal Sturm–Liouville problem with composite fractional derivatives. Math. Meth. Appl. Sci. 2021, 44, 1931–1941. [Google Scholar] [CrossRef]
  20. Malinowska, A.B.; Odzijewicz, T.; Torres, D.F.M. Advanced Methods in the Fractional Calculus of Variations; Springer International Publishing: London, UK, 2015. [Google Scholar]
  21. Pandey, P.K.; Pandey, R.K.; Agrawal, O.P. Variational approximation for fractional Sturm–Liouville problem. Fract. Calc. Appl. Anal. 2020, 23, 861–874. [Google Scholar] [CrossRef]
  22. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  23. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  24. Giusti, A.; Colombaro, I.; Garra, R.; Garappa, R.; Polito, F.; Popolizio, M.; Mainardi, F. A practical guide to Prabhakar fractional calculus. Fract. Calc. Appl. Anal. 2020, 23, 9–54. [Google Scholar] [CrossRef] [Green Version]
  25. Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 1–11. [Google Scholar]
  26. Klimek, M. On properties of exact and numerical solutions to integral eigenvalue problems associated to fractional differential ones. In Selected Topics in Contemporary Mathematical Modeling; Czestochowa University of Technology Press: Czestochowa, Poland, 2021; pp. 49–65. [Google Scholar]
  27. Stanislavsky, A.; Weron, A. Transient anomalous diffusion with Prabhakar-type memory. J. Chem. Phys. 2018, 149, 044107. [Google Scholar] [CrossRef] [PubMed]
  28. Dos Santos, M.A.F. Fractional Prabhakar derivative in diffusion equation with non-static stochastic resetting. Physics 2019, 1, 40–58. [Google Scholar] [CrossRef] [Green Version]
  29. Shukla, A.; Sukavanam, N.; Pandey, D.N. Approximate controllabillity of semilinear fractional control systems of order α∈(1,2] with infinite delay. Mediterr. J. Math. 2016, 13, 2539–2550. [Google Scholar] [CrossRef]
  30. Górska, K.; Horzela, A.; Lattanzi, A.; Pogany, T.K. On complete monotonicity of three parameter Mittag-Leffler function. Appl. Anal. Discret. Math. 2021, 15, 118–128. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Klimek, M. Spectrum of Fractional and Fractional Prabhakar Sturm–Liouville Problems with Homogeneous Dirichlet Boundary Conditions. Symmetry 2021, 13, 2265. https://doi.org/10.3390/sym13122265

AMA Style

Klimek M. Spectrum of Fractional and Fractional Prabhakar Sturm–Liouville Problems with Homogeneous Dirichlet Boundary Conditions. Symmetry. 2021; 13(12):2265. https://doi.org/10.3390/sym13122265

Chicago/Turabian Style

Klimek, Malgorzata. 2021. "Spectrum of Fractional and Fractional Prabhakar Sturm–Liouville Problems with Homogeneous Dirichlet Boundary Conditions" Symmetry 13, no. 12: 2265. https://doi.org/10.3390/sym13122265

APA Style

Klimek, M. (2021). Spectrum of Fractional and Fractional Prabhakar Sturm–Liouville Problems with Homogeneous Dirichlet Boundary Conditions. Symmetry, 13(12), 2265. https://doi.org/10.3390/sym13122265

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