On a Class of Positive Deﬁnite Operators and Their Application in Fractional Calculus

: This paper is devoted to the spectral analysis of one class of integral operators, associated with the boundary value problems for differential equations of fractional order. Approximation matrices are also investigated. In particular, the positive deﬁniteness of the studied operators is shown, which makes it possible to select areas in the complex plane where there are no eigenvalues of these operators.


Introduction
As is known in fractional calculus and in the theory of mixed-type equations, an important role is played by the potential function 1 2 1 0 u(t)dt |x − t| 1/ρ   with density u(t) and with a power kernel 1 |x−t| 1/ρ , which is positive definite for 0 < 1 ρ < 1 according to a fact established by Tricomi [1][2][3][4][5][6][7][8].There are papers in which various generalizations of this result were given.First of all, we should note the paper of Gellerstedt [9], where an operator of the following form was investigated for positive definiteness Another direction was started in [10,11] where, in particular, it was shown that the operator A ρ : L 2 → L 2 , for 0 < 1 ρ < 1 (where A ρ = x 0 (x − t) 1/ρ−1 u(t)dt) is sectorial and also that the values of the form ( A ρ u, u), for 1 < 1 ρ < ∞, fill the whole complex plane [10].This manuscript is devoted to the study of the positive definiteness of operators of the form which are finite-dimensional perturbations (finite-dimensional perturbations of a special kind) of a fractional integration operator of a special kind.We suggest a principally new wide class of positive definite operators, which play an important role in fractional calculus and in applications.The obtained results are used to study some very important properties of functions of the Mittag-Leffler type.

On the Positive Definiteness of Operators of the Kind A [α,β] γ
Let us consider the operator This operator arises in the process of determining the solutions of boundary value problems for fractional differential equations [12].
Let us show that this operator (for specific α, γ, ρ) is positive definite.To highlight the main ideas, let us consider the simplest case.Let us consider in space L 2 (0, 1) the operator A [α,β] γ for α = β = γ = ρ, 0 < ρ < 2, i.e., we consider the operator A ρ [13].The more important case is for 0 < ρ < 1, as in this case, the operator A ρ corresponds to the differential equations of order more than 1.The case for 1 < ρ < ∞ in fractional calculus is not so interesting but to complete our investigation, we will to consider some results for this case too.
First of all we, note that the first term of operator A ρ is a fractional integral J 1/ρ of order 1/ρ.
Let us designate Obviously, the operator A ρ is different from the operator J 1/ρ by the positive constant.However, in the following, for easy reading, we will use the operator A ρ and do not pay attention to this difference.As is known in fractional calculus and in the theory of mixed-type equations, an important role is played by the potential with density u(t) and with a power kernel 1 |x−t| (1/ρ) which is positive definite, for 0 < 1/ρ < 1, and this fact was established by F. Tricomi [8].
F. Tricomi [8] showed that the symmetric component of the operator A ρ , i.e., It should be noted that the operator A ρ is strictly definite ( A ρ u, u) > 0 (the equality sign holds if and only if u = 0).A little later, Matsaev and Palant [11] showed that the operator A ρ is sectorial (0 < 1/ρ < 1), that is, the values of the form ( A ρ u, u) lie in the angle |argλ| < π 2ρ .
Further, Gokhberg and Krein [10] showed that the values of the form ( A ρ u, u) for (1 < 1/ρ < ∞) fill the whole complex plane.This paper provides further analysis of these operators.
There is another well-known result.The operator A ρ , for 1 < ρ < 2 for (A ρ : Let us now formulate and prove the following theorem Proof.It is known that the operator A is called positive definite if (Au, u) > 0, (u = 0).However, it is very difficult to verify this condition directly.Therefore, we will use the matrix approximation of the operator A ρ [6].As in [6], we denote the corresponding matrix by The system of eigenfunctions of the K is complete in the domain of values of the integral operator K f if and only if both sides in this inequality are equal.
It is known that the operator A is called positive definite if (Au, u) > 0, (u = 0).However, it is very difficult to verify this condition directly.Therefore, we will use the matrix approximation of the operator A ρ [6].As in [6], we denote the corresponding matrix by .
The matrix T n−1 (µ) has many useful properties.In particular, this matrix is positive, persymmetric, indecomposable, etc.It is known [6] that one of the necessary conditions for the positive definiteness of a matrix is the positivity of all its lead main minors.The fact that these minors are positive was shown in [6].Next, we need the following lemmas.

Lemma 1. The minors
of the matrix (here µ > 1) Proof.Let us consider the minor For i k ≤ j k , 1 ≤ k ≥ r, we may represent M r as follows To calculate the determinant M r , we consider Here, It is clear that which proves the lemma.
To prove that the matrix T n−1 (µ) is positive definite, we have studied the real component of this matrix Using the high-level mathematical package MATLAB, the eigenvalues of the matrix T R were considered for various values of µ and the dimension of the matrix N. It was shown that all eigenvalues of the matrix T R , for any N ≤ 3000 and µ > 0, are positive, that is, the above calculations confirm the hypothesis that the matrix T R is positive definite.This became the basis for us to assume that the matrix T (n−1) (µ) under study is positive definite.It is natural that the operator A ρ corresponding to the matrix T (n−1) (µ) will also be positive definite.
We give a strong proof of the positive definiteness of the matrix T R (µ).First, let us write the matrices T 6 (1/2), T * 6 (1/2), T R (1/2) using the MATLAB package That is, the following statements hold.
Proof.We write the formula for the general element of the matrix Obviously, the elements under the main diagonal are calculated as follows, µ > 1 and the elements under the main diagonal are So, these formulas show that the elements located above the main diagonal decrease.To consider the elements under the main diagonal, we introduce the generating function Obviously, the derivative of this function is positive on the segment x ∈ [1, N], which means that the function ϕ(x) increases on the segment x ∈ [1, N].This completes the proof.Lemma 3.For any fixed j 0 < i, the following relations a i,j 0 ≥ a i+1,j 0 , i ≥ j 0 ; a i,j 0 < a i+1,j 0 , i < j 0 hold.
Proof.The proof of Lemma 2 is similar to the proof of Lemma 1.
Lemma 4. The statements of Lemmas 1 and 2 are valid for the matrices T T n (µ) (T T n (µ) is the transposed matrix, and it shall be reminded that T(µ) for µ = 1 was studied by Krein [10], for µ = 2, it was studied in [6]).
Proof.It is obvious that all main lead minors of the matrix T R (0) are non-negative.In the same way, all main lead minors of the matrix T R (1) are positive.
Let us show that for µ ∈ (1, 2), all main lead minors of the matrix are T R (µ) = 0. To do this, it is enough to prove that all the rows (columns) of the leading main lead minors of the matrix are linearly independent.In proving this statement, without loss of generality, for definiteness, we consider rows with numbers k and k + 1.Then, it suffices to note that, by Lemma 4, that a k,1 a k+1,1 < 1 and a k,k+1 a k,k > 1, which proves the linear independence of these rows.
From the above-provided statements follows that the matrix T n−1 (µ) is positive definite.So, the operator A ρ is positive definite too, and this proves the lemma.
From this very important theorem, it follows that the operator A ρ is positive definite for 1/2 < ρ < ∞.

Conclusions
So, our spectral analysis of the operators generated by boundary value problems for fractional differential equations and boundary conditions of Sturm-Liouville type, using matrix calculus, shows that the spectral structure of these operators can be studied by the matrices we studied above.

Discussion
Operators generated (induced) by a differential expression of a fractional order and boundary conditions of the Sturm-Liouville type are non-self-adjoint and their spectral structure is almost not studied.The methods proposed by the authors are fundamentally new.They allow to study the completeness of systems of eigenfunctions and associated functions of these operators.