The modeling and analysis of the COVID-19 pandemic with vaccination and isolation: a case study of Italy

1.   Introduction

Fractional calculus is a notably attractive subject owing to having wide-ranging application areas of theoretical and applied sciences. Despite the fact that there are a large number of worthwhile mathematical works on the fractional differential calculus, there is no noteworthy parallel improvement of fractional difference calculus up to lately. This statement has shown that discrete fractional calculus has certain unforeseen hardship.

Fractional sums and differences were obtained firstly in Diaz-Osler ^[1], Miller-Ross ^[2] and Gray and Zhang ^[3] and they found discrete types of fractional integrals and derivatives. Later, several authors began to touch upon discrete fractional calculus; Goodrich-Peterson ^[4], Baleanu et al. ^[5], Ahrendt et al. ^[6]. Nevertheless, discrete fractional calculus is a rather novel area. The first studies have been done by Atıcı et al. ^{[7,8,9,10,11]}, Abdeljawad et al. ^[12,13,14], Mozyrska et al. ^[15,16,17], Anastassiou ^[18,19], Hein et al. ^[20] and Cheng et al. ^[21] and so forth ^{[22,23,24,25,26]}.

Self-adjoint operators have an important place in differential operators. Levitan and Sargsian ^[27] studied self-adjoint Sturm-Liouville differential operators and they obtained spectral properties based on self-adjointness. Also, they found representation of solutions and hence they obtained asymptotic formulas of eigenfunctions and eigenvalues. Similarly, Dehghan and Mingarelli ^[28,29] obtained for the first time representation of solution of fractional Sturm-Liouville problem and they obtained asymptotic formulas of eigenfunctions and eigenvalues of the problem. In this study, firstly we obtain self-adjointness of DFSL operator within nabla fractional Riemann-Liouville and delta fractional Grünwald-Letnikov operators. From this point of view, we obtain orthogonality of distinct eigenfunctions, reality of eigenvalues. In addition, we open a new gate by obtaining representation of solution of DFSL problem for researchers study in this area.

Self-adjointness of fractional Sturm-Liouville differential operators have been proven by Bas et al. ^[30,31], Klimek et al. ^[32,33]. Variational properties of fractional Sturm-Liouville problem has been studied in ^[34,35]. However, self-adjointness of conformable Sturm-Liouville and DFSL with Caputo-Fabrizio operator has been proven by ^[36,37]. Nowadays, several studies related to Atangana-Baleanu fractional derivative and its discrete version are done ^{[38,39,40,41,42,43,44,45]}.

In this study, we consider DFSL operators within Riemann-Liouville and Grünwald-Letnikov sense, and we prove the self-adjointness, orthogonality of distinct eigenfunctions, reality of eigenvalues of DFSL operator. However, we get sum representation of solutions for DFSL equation by means Laplace transform for nabla fractional difference equations. Finally, we compare the results for the solution of DFSL problem, discrete Sturm-Liouville (DSL) problem with the second order, fractional Sturm-Liouville (FSL) problem and classical Sturm-Liouville (CSL) problem with the second order. The aim of this paper is to contribute to the theory of DFSL operator.

We discuss DFSL equations in three different ways with;

$i)$ Self-adjoint (nabla left and right) Riemann-Liouville (R-L) fractional operator,

$ii)$ Self-adjoint (delta left and right) Grünwald-Letnikov (G-L) fractional operator,

$iii)$(nabla left) DFSL operator is defined by R-L fractional operator,

2.   Preliminaries

Definition 2.1. ^[4] Delta and nabla difference operators are defined by respectively

Definition 2.2. ^[46] Falling function is defined by, $\alpha \in \mathbb{R}$

where $\Gamma$ is Euler gamma function.

Definition 2.3. ^[46] Rising function is defined by, $\alpha \in \mathbb{R}$,

Remark 1. Delta and nabla operators have the following properties

Definition 2.4. ^[2,7] Fractional sum operators are defined by,

$\left(i\right)$ The left defined nabla fractional sum with order $\mu > 0$ is defined by

$\left(ii\right)$ The right defined nabla fractional sum with order $\mu > 0$ is defined by

where $\rho \left(t\right) = t-1$ is called backward jump operators, $\mathbb{N} _{a} = \left \{ a, a+1, ...\right \},$ $_{b} \mathbb{N} = \left \{ b, b-1, ...\right \}$.

Definition 2.5. ^[47] Fractional difference operators are defined by,

$\left(i\right)$ The nabla left fractional difference of order $\mu > 0$ is defined

$\left(ii\right)$ The nabla right fractional difference of order $\mu > 0$ is defined

Fractional differences in $\left(7-8\right)$ are called the Riemann-Liouville (R-L) definition of the $\mu$-th order nabla fractional difference.

Definition 2.6. ^[1,21,48] Fractional difference operators are defined by,

$\left(i\right)$ The left defined delta fractional difference of order $\mu,$ $0 < \mu \leq 1,$ is defined by

$\left(ii\right)$ The right defined delta fractional difference of order $\mu,$ $0 < \mu \leq 1,$ is defined by

Fractional differences in $\left(9-10\right)$ are called the Grünwald-Letnikov (G-L) definition of the $\mu$-th order delta fractional difference.

Theorem 2.7. ^[47] We define the summation by parts formula for R-L fractional nabla difference operator, $u$ is defined on $_{b} \mathbb{N}$ and $v$ is defined on $\mathbb{N} _{a}$, then

Theorem 2.8. ^[26,48] We define the summation by parts formula for G-L delta fractional difference operator, $u$, $v$ is defined on $\left \{ 0, 1, ..., n\right \}$, then

Definition 2.9. ^[20] $f: \mathbb{N} _{a}\rightarrow \mathbb{R}$, $s\in \Re,$ Laplace transform is defined as follows,

where $\Re = \mathbb{C} \backslash \left \{ 1\right \}$ and $\Re$ is called the set of regressive (complex) functions.

Definition 2.10. ^[20] Let $f, g: \mathbb{N}_{a}\rightarrow \mathbb{R}$, all $t\in \mathbb{N}_{a+1},$ convolution property of $f$ and $g$ is given by

where $\rho \left(s\right)$ is the backward jump function defined in ^[46] as

Theorem 2.11. ^[20] $f, g: \mathbb{N} _{a}\rightarrow \mathbb{R},$ convolution theorem is expressed as follows,

Lemma 2.12. ^[20] $f: \mathbb{N} _{a}\rightarrow \mathbb{R},$ the following property is valid,

Theorem 2.13. ^[20] $f: \mathbb{N} _{a}\rightarrow \mathbb{R},$ $0 < \mu < 1,$ Laplace transform of nabla fractional difference

Definition 2.14. ^[20] For $\left \vert p\right \vert < 1$, $\alpha > 0,$ $\beta \in \mathbb{R}$ and $t\in \mathbb{N} _{a}$, discrete Mittag-Leffler function is defined by

where $t^{\overline{n}} = \left \{ \begin{array}{cc} t\left(t+1\right) \cdots \left(t+n-1\right), & n\in \mathbb{Z} \\ \frac{\Gamma \left(t+n\right) }{\Gamma \left(t\right) }, & n\in \mathbb{R} \end{array} \right.$ is rising factorial function.

Theorem 2.15. ^[20] For $\left \vert p\right \vert < 1$, $\alpha > 0,$ $\beta \in \mathbb{R}$, $\left \vert 1-s\right \vert < 1$, and $\left \vert s\right \vert ^{\alpha } > p$, Laplace transform of discrete Mittag-Leffler function is as follows,

Definition 2.16. Laplace transform of $f\left(t\right) \in \mathbb{R} ^{+},$ $t\geq 0$ is defined as follows,

Theorem 2.17. For $z,$ $\theta \in \mathbb{C}, {Re}(\delta) > 0$, Mittag-Leffler function with two parameters is defined as follows

Theorem 2.18. Laplace transform of Mittag-Leffler function is as follows

Property 2.19. ^[28] $f: \mathbb{N} _{a}\rightarrow \mathbb{R},$ $0 < \mu < 1,$ Laplace transform of fractional derivative in Caputo sense is as follows, $0 < \alpha < 1,$

Property 2.20. ^[28] $f: \mathbb{N} _{a}\rightarrow \mathbb{R},$ $0 < \mu < 1,$ Laplace transform of left fractional derivative in Riemann-Liouville sense is as follows, $0 < \alpha < 1,$

here $I_{0^{+}}^{\alpha }$ is left fractional integral in Riemann-Liouville sense.

3.   Main results

3.1.   Discrete fractional Sturm-Liouville equations

We consider discrete fractional Sturm-Liouville equations in three different ways as follows:

$First$ $Case:$ Self-adjoint $L_{1}$ DFSL operator is defined by (nabla right and left) R-L fractional operator,

where $p\left(t\right) > 0,$ $r\left(t\right) > 0,$ $q\left(t\right)$ is a real valued function on $\left[a+1, b-1\right]$ and real valued, $\lambda$ is the spectral parameter, $t\in \left[a+1, b-1\right],$ $x\left(t\right) \in l^{2}\left[a+1, b-1\right].$ In $\ell ^{2}\left(a+1, b-1\right),$ the Hilbert space of sequences of complex numbers $u\left(a+1\right), ..., u\left(b-1\right)$ with the inner product is given by,

for every $u\in D_{L_{1}}$, let's define as follows

$Second$ $Case:$ Self-adjoint $L_{2}$ DFSL operator is defined by(delta left and right) G-L fractional operator,

where $p, r, \lambda$ is as defined above, $q\left(t\right)$ is a real valued function on $\left[0, n\right],$ $t\in \left[0, n\right],$ $x\left(t\right) \in l^{2}\left[0, n\right].$ In $\ell ^{2}\left(0, n\right),$ the Hilbert space of sequences of complex numbers $u\left(0\right), ..., u\left(n\right)$ with the inner product is given by, $n$ is a finite integer,

for every $u\in D_{L_{2}},$ let's define as follows

$Third$ $Case:$$L_{3}$ DFSL operator is defined by (nabla left) R-L fractional operator,

$p, r, \lambda$ is as defined above, $q\left(t\right)$ is a real valued function on $\left[a+1, b-1\right],$ $t\in \left[a+1, b-1\right].$

Firstly, we consider the first case and give the following theorems and proofs$;$

Theorem 3.1. DFSL operator $L_{1}$ is self-adjoint.

Proof.

If $\left(16-17\right)$ is subtracted from each other

and sum operator from $a+1$ to $b-1$ to both side of the last equality is applied, we get

If we apply the summation by parts formula in $\left(11\right)$ to right hand side of $\left(18\right),$ we have

Hence, the proof completes.

Theorem 3.2. Two eigenfunctions, $u(t, \lambda_{\alpha})$ and $v(t, \lambda_{\beta})$, of the equation $\left(13\right)$ are orthogonal as $\lambda_{\alpha}\neq\lambda_{\beta}$.

Proof. Let $\lambda _{\alpha }$ and $\lambda _{\beta }$ are two different eigenvalues corresponds to eigenfunctions $u\left(t\right)$ and $v\left(t\right)$ respectively for the the equation $\left(13\right)$,

If we multiply last two equations by $v\left(t\right)$ and $u\left(t\right)$ respectively, subtract from each other and apply definite sum operator, owing to the self-adjointness of the operator $L_{1},$ we have

since $\lambda _{\alpha }\neq \lambda _{\beta, }$

Hence, the proof completes.

Theorem 3.3. All eigenvalues of the equation $\left(13\right)$ are real.

Proof. Let $\lambda = \alpha +i\beta,$ owing to the self-adjointness of the operator $L_{1},$ we can write

Since $\left \langle u\left(t\right), u\left(t\right) \right \rangle _{r}\neq 0,$

and hence $\beta = 0.$ The proof completes.

Secondly, we consider the second case and give the following theorems and proofs$;$

Theorem 3.4. DFSL operator $L_{2}$ is self-adjoint.

Proof.

If $\left(19-20\right)$ is subtracted from each other

and definite sum operator from $0$ to $t$ to both side of the last equality is applied, we have

If we apply the summation by parts formula in $\left(12\right)$ to r.h.s. of $\left(21\right),$ we get

Hence, the proof completes.

Theorem 3.5. Two eigenfunctions, $u(t, \lambda_{\alpha})$ and $v(t, \lambda_{\beta})$, of the equation $\left(14\right)$ are orthogonal as $\lambda_{\alpha}\neq\lambda_{\beta}$. orthogonal.

Proof. Let $\lambda _{\alpha }$ and $\lambda _{\beta }$ are two different eigenvalues corresponds to eigenfunctions $u\left(t\right)$ and $v\left(t\right)$ respectively for the the equation $\left(14\right)$,

If we multiply last two equations to $v\left(t\right)$ and $u\left(t\right)$ respectively, subtract from each other and apply definite sum operator, owing to the self-adjointness of the operator $L_{2},$ we get

since $\lambda _{\alpha }\neq \lambda _{\beta, }$

So, the eigenfunctions are orthogonal. The proof completes.

Theorem 3.6. All eigenvalues of the equation $\left(14\right)$ are real.

Proof. Let $\lambda = \alpha +i\beta,$ owing to the self-adjointness of the operator $L_{2}$

Since $\left \langle u, u\right \rangle _{r}\neq 0,$

and hence $\beta = 0.$ The proof completes.

3.2.   Sum representation of solution of discrete fractional Sturm-Liouville problem

Now, we consider the third case and give the following theorem and proof$;$

Theorem 3.7.

where $p\left(t\right) > 0,$ $r\left(t\right) > 0,$ $q\left(t\right)$ is defined and real valued, $\lambda$ is the spectral parameter$.$ The sum representation of solution of the problem $\left(22\right) -\left(23\right)$ is found as follows,

where $\left \vert \lambda \right \vert < 1$, $\left \vert 1-s\right \vert < 1$, and $\left \vert s\right \vert ^{\alpha } > \lambda$ from Theorem 2.15.

Proof. Let's use the Laplace transform of both side of the equation $\left(22\right)$ by Theorem 2.13, and let $q\left(t\right) x\left(t\right) = g\left(t\right),$

from Lemma 2.12, we get

Applying inverse Laplace transform to the equation $\left(25 \right)$, then we get representation of solution of the problem $\left(22\right) -\left(23\right),$

4.   General discussions

Now, let us consider comparatively discrete fractional Sturm-Liouville (DFSL) problem, discrete Sturm-Liouville (DSL) problem, fractional Sturm-Liouville (FSL) problem and classical Sturm-Liouville (CSL) problem respectively as follows by taking $q\left(t\right) = 0,$

$DFSL$ problem:

and its analytic solution is as follows by the help of Laplace transform in Lemma 2.12

$DSL$ problem$:$

and its analytic solution is as follows

$FSL$ problem$:$

and its analytic solution is as follows by the help of Laplace transform in Property 2.19 and 2.20

$CSL$ problem$:$

and its analytic solution is as follows

where the domain and range of function $x\left(t\right)$ and Mittag-Leffler functions must be well defined. Note that we may show the solution of CSL problem can be obtained by taking $\mu \rightarrow 1$ in the solution of FSL problem and similarly, the solution of DSL problem can be obtained by taking $\mu \rightarrow 1$ in the solution of DFSL problem.

Firstly, we compare the solutions of DFSL and DSL problems and from here we show that the solutions of DFSL problem converge to the solutions of DSL problem as $\mu \rightarrow 1$ in Figure 1 for discrete Mittag-Leffler function $E_{p, \alpha, \beta }\left(t, a\right) = \sum \limits_{k = 0}^{1000}p^{k}\frac{\left(t-a\right) ^{\overline{\alpha k+\beta }}}{\Gamma \left(\alpha k+\beta +1\right) }$; let $\lambda = 0.01,$

Secondly, we compare the solutions of DFSL, DSL, FSL and CSL problems for discrete Mittag-Leffler function $E_{p, \alpha, \beta }\left(t, a\right) = \sum \limits_{k = 0}^{1000}p^{k}\frac{ \left(t-a\right) ^{\overline{\alpha k+\beta }}}{\Gamma \left(\alpha k+\beta +1\right) }$. At first view, we observe the solution of DSL and CSL problems almost coincide in any order $\mu$, and we observe the solutions of DFSL and FSL problem almost coincide in any order $\mu$. However, we observe that all of the solutions of DFSL, DSL, FSL and CSL problems almost coincide to each other as $\mu \rightarrow 1$ in Figure 2. Let $\lambda = 0.01,$

Thirdly, we compare the solutions of DFSL problem $\left(22-23\right)$ with different orders, different potential functions and different eigenvalues for discrete Mittag-Leffler function $E_{p, \alpha, \beta }\left(t, a\right) = \sum \limits_{k = 0}^{1000}p^{k}\frac{\left(t-a\right) ^{\overline{\alpha k+\beta }} }{\Gamma \left(\alpha k+\beta +1\right) }$ in the Figure 3;

Eigenvalues of DFSL problem $\left(22-23 \right),$ correspond to some specific eigenfunctions for numerical values of discrete Mittag-Leffler function $E_{p, \alpha, \beta }\left(t, a\right) = \sum \limits_{k = 0}^{i}p^{k}\frac{\left(t-a\right) ^{\overline{\alpha k+\beta }}}{\Gamma \left(\alpha k+\beta +1\right) },$ is given with different orders while $q\left(t\right) = 0$ in Table 1;

Finally, we give the solutions of DFSL problem $\left(22-23\right)$ with different orders, different potential functions and different eigenvalues for discrete Mittag-Leffler function $E_{p, \alpha, \beta }\left(t, a\right) = \sum \limits_{k = 0}^{100}p^{k}\frac{\left(t-a\right) ^{\overline{\alpha k+\beta }} }{\Gamma \left(\alpha k+\beta +1\right) }$ in Tables 2–4;

4.1.   Discussions on eigenvalues and eigenfunctions of DFSL, DSL, FSL, CSL problems

Now, let's consider the problems together DFSL $\left(26)-(27 \right),$ DSL $\left(29)-(30\right)$, FSL $\left(32)-(33\right)$ and CSL $\left(35)-(36\right).$ Eigenvalues of these problems are the roots of the following equation

Thus, if we apply the solutions $\left(28\right),$ $\left(31 \right),$ $\left(34\right)$ and $\left(37\right)$ of these four problems to the equation above respectively, we can find the eigenvalues of these problems for the orders $\mu = 0.9$ and $\mu = 0.99$ respectively in Table 5, and Table 6,

In here, we observe that these four problems have real eigenvalues under different orders $\mu = 0.9$ and $\mu = 0.99$, hence we can find eigenfunctions putting these eigenvalues into the four solutions. Furthermore, as the order changes, we can see that eigenvalues change for DFSL problems.

5.   Conclusion

We consider firstly discrete fractional Sturm-Liouville (DFSL) operators with nabla Riemann-Liouville and delta Grünwald-Letnikov fractional operators and we prove self-adjointness of the DFSL operator and fundamental spectral properties. However, we analyze DFSL problem, discrete Sturm-Liouville (DSL) problem, fractional Sturm-Liouville (FSL) problem and classical Sturm-Liouville (CSL) problem by taking $q\left(t\right) = 0$ in applications. Firstly, we compare the solutions of DFSL and DSL problems and we observe that the solutions of DFSL problem converge to the solutions of DSL problem when $\mu \rightarrow 1$ in Fig. 1. Secondly, we compare the solutions of DFSL, DSL, FSL and CSL problems in Fig. 2. At first view, we observe the solutions of DSL and CSL problems almost coincide with any order $\mu$, and we observe the solutions of DFSL and FSL problem almost coincide with any order $\mu$. However, we observe that all of solutions of DFSL, DSL, FSL and CSL problems almost coincide with each other as $\mu \rightarrow 1$. Thirdly, we compare the solutions of DFSL problem $\left(22-23\right)$ with different orders, different potential functions and different eigenvalues in Fig. 3.

Eigenvalues of DFSL problem $\left(22-23\right)$ corresponded to some specific eigenfunctions is given with different orders in Table 1. We give the eigenfunctions of DFSL problem $\left(22-23\right)$ with different orders, different potential functions and different eigenvalues in Table 2, Table 3 and Table 4.

In Section 4.1, we consider DFSL, DSL, FSL and CSL problems together and thus, we can compare the eigenvalues of these four problems in Table 5 and Table 6 for different values of $\mu$. We observe that these four problems have real eigenvalues under different values of $\mu$, from here we can find eigenfunctions corresponding eigenvalues. Moreover, when the order change, eigenvalues change for DFSL problems.

Consequently, important results in spectral theory are given for discrete Sturm-Liouville problems. These results will lead to open gates for the researchers studied in this area. Especially, representation of solution will be practicable for future studies. It worths noting that visual results both will enable to be understood clearly by readers and verify the results to the integer order discrete case while the order approaches to one.

Acknowledgment

This paper includes a part of Ph.D. thesis data of Ramazan OZARSLAN.

Conflict of interest

The authors declare no conflict of interest.