Next Article in Journal
Map of a Bending Problem for Self-Similar Beams into the Fractal Continuum Using the Euler–Bernoulli Principle
Previous Article in Journal
Stability Analysis of the Nabla Distributed-Order Nonlinear Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution Equations in Hilbert Spaces via the Lacunae Method

by
Maksim V. Kukushkin
Moscow State University of Civil Engineering, 129337 Moscow, Russia
Fractal Fract. 2022, 6(5), 229; https://doi.org/10.3390/fractalfract6050229
Submission received: 11 March 2022 / Revised: 11 April 2022 / Accepted: 16 April 2022 / Published: 20 April 2022

Abstract

:
In this paper, we consider evolution equations in the abstract Hilbert space under the special conditions imposed on the operator at the right-hand side of the equation. We establish the method that allows us to formulate the existence and uniqueness theorem and find a solution in the form of a series on the root vectors of the right-hand side. We consider fractional differential equations of various kinds as an application. Such operators as the Riemann-Liouville fractional differential operator, the Riesz potential, the difference operator have been involved.

1. Introduction

In the paper [1], we obtained the clarification of the results by Lidskii [2] on the decomposition on the root vector system of the non-selfadjoint operator. We used a technique of the entire function theory and introduce a so-called Schatten-von Neumann class of the convergence exponent. Considering strictly accretive operators satisfying special conditions formulated in terms of the norm, we constructed a sequence of contours of the power type on the contrary to the results by Lidskii [2], where a sequence of contours of the exponential type was used.

1.1. Motivation for the Current Study

In this paper, we produce the application of the mentioned method to evolution equations in the abstract Hilbert space with the right-hand side of the special type. Here, to show our motivation, we should appeal to a plenty of applications to concrete differential equations connected with modeling various physical-chemical processes: filtration of liquid and gas in the highly porous fractal medium; heat exchange processes in a medium with fractal structure and memory; casual walks of a point particle that starts moving from the origin by self-similar fractal set; oscillator motion under the action of elastic forces which is characteristic for a viscoelastic medium, etc. In particular, we would like to study the existence and uniqueness theorems for evolution equations with the right-hand side—a differential operator with a fractional derivative in final terms. In this connection such operators as a Riemann-Liouville fractional differential operator, Kipriyanov operator, Riesz potential, difference operator are involved. Nowadays, the concept of the fractional integro-differentiation is efficiently used in the study of problems related with the approximate controllability of systems with damping [3], approximate controllability of semi-linear stochastic systems [4,5], partial-approximate controllability of semi-linear systems [6], asymptotic stability of stochastic differential equations in Banach spaces [7].

1.2. Novelty for the Current Study

Note that analysis of the required conditions imposed upon the right-hand side of the studied class of evolution equations deserves to be mentioned. In this regard, we should note a well-known fact (see for instance [8]) that a particular interest appears in the case when a senior term of the operator (see [9]) is not selfadjoint at least, for in the contrary case there is a plenty of results devoted to the topic, within the framework of which the following papers are well-known [8,10,11,12,13]. Indeed, most of them deal with a decomposition of the operator to a sum, where the senior term must be either a selfadjoint or normal operator. In other cases, the methods of the papers [9,14] become relevant and allow us to study spectral properties of operators, whether we have the mentioned above representation or not. Here, we ought to stress that the results of the papers [12,15] can be also applied to study non-selfadjoin operators but based on the sufficiently strong assumption regarding the numerical range of values of the operator (the numerical range belongs to a parabolic domain). The methods of [9] that are applicable to study non-selfadjoint operators can be used in the natural way, if we deal with a more abstract construction—the infinitesimal generator of a semigroup of contraction [16]. The central challenge of the latter paper is how to create a model representing a composition of fractional differential operators in terms of the semigroup theory. Here, we should note that motivation arouse in connection with the fact that a second-order differential operator can be presented as a some kind of a transform of the infinitesimal generator of a shift semigroup and stress that the eigenvalue problem for the operator was previously studied by methods of the theory of functions [17,18]. Having been inspired by novelty of the idea, we generalize a differential operator with a fractional integro-differential composition in the final terms to some transform of the corresponding infinitesimal generator of the shift semigroup. By virtue of the methods obtained in the paper [9], we managed to study spectral properties of the infinitesimal generator transform and obtained an outstanding result–asymptotic equivalence between the real component of the resolvent and the resolvent of the real component of the operator. The relevance is based on the fact that the asymptotic formula for the operator real component can be established in most cases due to well-known asymptotic relations for the regular differential operators as well as for the singular ones [19]. Thus, we have theorems establishing spectral properties of some class of non-selfadjoint operators which allow us, jointly with the results [1], to study the Cauchy problem for the evolution equation by the functional analysis methods. Note that the abstract approach to the Cauchy problem for the fractional evolution equation was previously implemented [20,21]. However, the main advantage of this paper is the obtained formula for the solution of the evolution equation with the relatively wide conditions imposed upon the right-hand side, wherein the derivative at the left-hand side is supposed to be of the real order. We consider the evolution equations with the right-hand side—an operator function of the power type. This problem appeals to many ones that lie in the framework of the theory of differential equations. For instance, in the paper [22] the solution of the evolution equation modeling the switching kinetics of ferroelectrics in the injection mode can be obtained in the analytical way, if we impose the conditions upon the right-hand side. The following papers deal with equations which can be studied by the obtained in this paper abstract method [23,24,25,26,27]. Thus, we can claim that the offered approach is undoubtedly novel and relevant.

2. Preliminaries

2.1. Basic Notations

Let C , C i , i N 0 be positive constants. We assume that a value of C can be different in various formulas and parts of formulas but values of C i are certain. Denote by Fr M the set of boundary points of the set M . Everywhere further, if the contrary is not stated, we consider linear densely defined operators acting on a separable complex Hilbert space H . Denote by B ( H ) the set of linear bounded operators on H . Denote by L ˜ the closure of an operator L . We establish the following agreement on using symbols L ˜ i : = ( L ˜ ) i , where i is an arbitrary symbol. Denote by D ( L ) , R ( L ) , N ( L ) the domain of definition, the range, and the kernel or null space of an operator L , respectively. The deficiency (codimension) of R ( L ) , dimension of N ( L ) are denoted by def L , nul L respectively. Assume that L is a closed operator acting on H , N ( L ) = 0 , let us define a Hilbert space H L : = f , g D ( L ) , ( f , g ) H L = ( L f , L g ) H . Consider a pair of complex Hilbert spaces H , H + , the notation H + H means that H + is dense in H as a set of elements and we have a bounded embedding provided by the inequality f H C 0 f H + , C 0 > 0 , f H + , moreover any bounded set with respect to the norm H + is compact with respect to the norm H . Let L be a closed operator, for any closable operator S such that S ˜ = L , its domain D ( S ) will be called a core of L . Denote by D 0 ( L ) a core of a closeable operator L . Let P ( L ) be the resolvent set of an operator L and R L ( ζ ) , ζ P ( L ) , [ R L : = R L ( 0 ) ] denotes the resolvent of an operator L . Denote by λ i ( L ) , i N the eigenvalues of an operator L . Suppose L is a compact operator and N : = ( L L ) 1 / 2 , r ( N ) : = dim R ( N ) ; then the eigenvalues of the operator N are called the singular numbers (s-numbers) of the operator L and are denoted by s i ( L ) , i = 1 , 2 , , r ( N ) . If r ( N ) < , then we put by definition s i = 0 , i = r ( N ) + 1 , 2 , . Let ν ( L ) denotes the sum of all algebraic multiplicities of an operator L . Denote by n ( r ) a function equals to the quantity of the elements of the sequence { a n } 1 , | a n within the circle | z < r . Let A be a compact operator, denote by n A ( r ) , counting function a function n ( r ) corresponding to the sequence { s i 1 ( A ) } 1 . Let S p ( H ) , 0 < p < be a Schatten-von Neumann class and S ( H ) be the set of compact operators. Denote by S ˜ ρ ( H ) the class of the operators such that A S ˜ ρ ( H ) { A S ρ + ε , A ¯ S ρ ε , ε > 0 } . In accordance with [1], we will call it Schatten-von Neumann class of the convergence exponent. Suppose L is an operator with a compact resolvent and s n ( R L ) C n μ , n N , 0 μ < ; then we denote by μ ( L ) order of the operator L (see [8]). Denote by Re L : = ( L + L ) / 2 , Im L : = ( L L ) / 2 i the real and imaginary components of an operator L respectively. In accordance with the terminology of the monograph [28], the set Θ ( L ) : = { z C : z = ( L f , f ) H , f D ( L ) , f H = 1 } is called the numerical range of an operator L . An operator L is called sectorial if its numerical range belongs to a closed sector L ι ( θ ) : = { ζ : | arg ( ζ ι ) θ < π / 2 } , where ι is the vertex and θ is the semi-angle of the sector L ι ( θ ) . If we want to stress the correspondence between ι and θ , then we will write θ ι . An operator L is called bounded from below if the following relation holds Re ( L f , f ) H γ L f H 2 , f D ( L ) , γ L R , where γ L is called a lower bound of L . An operator L is called accretive if γ L = 0 . An operator L is called strictly accretive if γ L > 0 . An operator L is called m-accretive if the next relation holds ( A + ζ ) 1 B ( H ) , ( A + ζ ) 1 ( Re ζ ) 1 , Re ζ > 0 . An operator L is called symmetric if one is densely defined and the following equality holds ( L f , g ) H = ( f , L g ) H , f , g D ( L ) . Consider a sesquilinear form t [ · , · ] (see [28]) defined on a linear manifold of the Hilbert space H . Let h = ( t + t ) / 2 , k = ( t t ) / 2 i be a real and imaginary component of the form t respectively, where t [ u , v ] = t [ v , u ] ¯ , D ( t ) = D ( t ) . Denote by t [ · ] the quadratic form corresponding to the sesquilinear form t [ · , · ] . According to these definitions, we have h [ · ] = Re t [ · ] , k [ · ] = Im t [ · ] . Denote by t ˜ the closure of a form t . The range of a quadratic form t [ f ] , f D ( t ) , f H = 1 is called the range of the sesquilinear form t and is denoted by Θ ( t ) . A form t is called sectorial if its range belongs to a sector having a vertex ι situated at the real axis and a semi-angle 0 θ < π / 2 . Due to Theorem 2.7 ([28], p. 323) there exist unique m-sectorial operators T t , T h associated with the closed sectorial forms t , h respectively. The operator T h is called a real part of the operator T t and is denoted by Re T t . Considering an arbitrary compact operator B : H H such that Θ ( B ) L 0 ( θ ) , π < θ < π , we put the following contour in correspondence to the operator
γ ( B ) : = λ : | λ = r > 0 , | arg λ θ + ε λ : | λ > r , | arg λ = θ + ε ,
where ε > 0 is an arbitrary small number, the number r is chosen so that the operator ( I λ B ) 1 is regular within the corresponding closed circle. Everywhere further, unless otherwise stated, we use notations of the papers [28,29,30,31,32]. Consider the following hypotheses regarding a non-selfadjoint operator.
(H1)
There exists a Hilbert space H + H and a linear manifold M that is dense in H + . The operator L is defined on M .
(H2)
( L f , g ) H C 1 f H + g H + , Re ( L f , f ) H C 2 f H + 2 , f , g M , C 1 , C 2 > 0 .
Throughout the paper we consider a restriction W of the operator L on the set M . We also use the short-hand notations A : = R W ˜ , μ : = μ ( H ) , where H : = Re W ˜ .

2.2. Way of Approach

Firstly, we consider general statements proved in [1] with the made refinement related to the involved notion of the convergence exponent as well as newly constructed sequence of contours allowing to arrange the eigenvectors in the power type way, we used this expression following the literary style of the monograph [2]. We implement the approach that refers us to the notion—operator order, it gives us an opportunity to reformulate results of the spectral theory in the more convenient and applicable way. Recall that in the paper [2] there was considered a sequence of contours of the exponential type, the condition α > ρ (here and further ρ denotes the index of the Schatten-von Neumann class of the convergence exponent) is imposed (see [1,2]). We improved this result in the paper [1] in the following sense, we produced a sequence of the power type contours what gives us the opportunity to obtain a solution of the problem in the case α = ρ . Moreover, the used method allows us to claim that the results are well-applicable to the class of non-selfadjoint operators with the numerical range of values contained in a parabolic domain. Such a significant achievement is obtained by virtue of the way of choosing a sector, containing the parabolic domain, with the infinitely small semi-angle. However, let us consider more general assumptions regarding the numerical range of values. Note that using the condition H 2 , it is not hard to prove that
Re ( W ˜ f , f ) H k | Im ( W ˜ f , f ) H ( C 2 k C 1 ) f H + 2 , k > 0 ,
from what follows the fact Θ ( A ) L 0 ( θ ) , θ = arctan ( C 1 / C 2 ) . Generally, the last relation gives us a range of the semi-angle π / 4 < θ < π / 2 , hence the conditions H 1 , H 2 are not sufficient to guaranty a value of the semi-angle less than π / 4 . At the same time some relevant results can be obtained in the case of sufficiently small values of the semi-angle, it gives us a motivation to consider an additional assumption.
(H3)
| Im ( L ˜ f , g ) H C 3 f H + g H , , f , g M , C 3 > 0 .
In this case, we have Re ( W ˜ f , f ) H k | Im ( W ˜ f , f ) H C 2 f H + k C 3 { ε f H + 2 / 2 + f H 2 / 2 ε } ( C 2 k ε C 3 ) | f H + 2 / 2 + ( C 2 / C 0 k C 3 / ε ) f H 2 / 2 , k > 0 . Thus, choosing ε = C 2 / k C 3 , we get Θ ( W ˜ ) L ι ( θ ι ) , where ι = C 2 / 2 C 0 ( k C 3 ) 2 / 2 C 2 , θ ι = arctan ( 1 / k ) . This relation guarantees that we can choose a sufficiently small value of the semi-angle θ ι . We put the following contour in correspondence to an operator L satisfying the additional condition H 3
Γ ( A ) : = Fr ( L ι ( θ ι + ε ) L 0 ( θ 0 + ε ) ) C r , C r : = λ : | λ < r , | arg λ θ 0 ,
where r is chosen so that the operator ( E λ A ) 1 is regular within the corresponding closed circle, ε > 0 is sufficiently small.

2.3. Basic Results, Theorems from the Referred Articles

The auxiliary theorems given bellow (see [1]) give us a tool to study the existence and uniqueness theorems in the abstract Hilbert space. Moreover, we obtain a solution that can be presented by the series on the operator A root vectors e q ξ + i with the coefficients c q ξ + i , where the index q relates to the eigenvalue, the index ξ relates to the geometrical multiplicity, the index i relates to the algebraic multiplicity, the convergence is understood in the Abel-Lidskii sense (see [2]). The following relation is a central aim of the Abel-Lidslii concept of decomposition on the root vector series
1 2 π i ϖ e λ α t A ( I λ A ) 1 f d λ = ν = 0 q = N ν + 1 N ν + 1 ξ = 1 m ( q ) i = 0 k ( q ξ ) e q ξ + i c q ξ + i ( t ) , f H ,
where ϖ is a contour in the complex plane, α > 0 , { N ν } 0 is a sequence of natural numbers,
ν = 0 q = N ν + 1 N ν + 1 ξ = 1 m ( q ) i = 0 k ( q ξ ) e q ξ + i c q ξ + i ( t ) H < .
The theorems given bellow devoted to the proof of the latter relation, they are formulated under the assumption that either the condition Θ ( A ) L 0 ( θ ) , θ < π / 2 α or the condition H 3 holds. In accordance with such an alternative, we put in correspondence ϖ : = γ ( A ) or ϖ : = Γ ( A ) respectively. The idea of reasonings belongs to Lidskii, however, we produce the proofs in [1], since the made refinement corresponding to the case, when ρ does not equal the index of the Schatten-von Neumann class, deserves to be considered itself.
Theorem 1.
Assume that Θ ( A ) L 0 ( θ ) , θ < π / 2 α , A S ˜ ρ , 0 < ρ α . Moreover in the case A S ˜ ρ S ρ the additional condition holds
n A m + 1 ( r m + 1 ) r ρ 0 , m = [ ρ ] .
Then a sequence of natural numbers { N ν } 0 can be chosen so that relation (1) holds.
Theorem 2.
Assume that the operator W ˜ satisfies the condition α > 2 / μ , μ ( 0 , 1 ] and α > 1 , μ ( 1 , ) . Then a sequence of natural numbers { N ν } 0 can be chosen so that relation (1) holds. Moreover, the following inequality holds for the eigenvalues
| λ N ν + k | λ N ν + k 1 | C λ N ν + k 1 1 / τ , k = 2 , 3 , , N ν + 1 N ν , 0 < τ < μ .
Theorem 3.
Assume that a normal operator W ˜ satisfies the condition ( ln 1 + 1 / α x ) λ i ( H ) = o ( i 1 / α ) , α > 1 . Then a sequence of natural numbers { N ν } 0 can be chosen so that relation (1) holds. Moreover, the following inequality holds for the eigenvalues
| λ N ν + k | λ N ν + k 1 | C λ N ν + k 1 1 / τ , k = 2 , 3 , , N ν + 1 N ν , 0 < τ < 1 / α .

3. Main Results

In this section, we consider evolution equations in the abstract Hilbert space. Having used an abstract theorem formulated in terms of the operator order, we produce an example of the class of differential equations for which the made refinement regarding the convergence exponent is relevant. More precisely, under the assumption ρ = α the sequence of contours may be chosen in a concrete—power type way, what provides a peculiar validity of the statement. We prove the existence and uniqueness theorem and supply it with a plenty of applications. We consider applications to the differential equations in the concrete Hilbert spaces and involve such operators as Riemann-Liouville operator, Kipriyanov operator, Riesz potential, difference operator. Moreover, we produce the artificially constructed normal operator for which the clarification of the Lidskii results relevantly works. Further, we will consider a Hilbert space H which consists of element-functions u : R + H , u : = u ( t ) , t 0 and we will assume that if u belongs to H then the fact holds for all values of the variable t . Notice that under such an assumption all standard topological properties as completeness, compactness etc. remain correctly defined. We understand such operations as differentiation and integration in the generalized sense that is caused by the topology of the Hilbert space H , more detailed information can be found in the Chapter 4 [33]. Consider a Cauchy problem
d u d t = W ˜ n u , u ( 0 ) = h D ( W ˜ ) , n = 1 , 2 , .
In the case when W ˜ n is accretive, we can suppose h H (here we should note that the case n = 1 was considered by Lidskii [2]). Bellow, we formulate a result which follows from the auxiliary Theorems 1–3 given above.
Theorem 4.
Assume that the conditions of one of the Theorems 1–3 hold under the assumption α = n , then there exists a solution of the Cauchy problem (2) in the form (1), a sequence of natural numbers { N ν } 0 can be chosen in accordance with the claim of the corresponding theorem. Moreover, the existing solution is unique, if the operator W ˜ n is accretive.
Proof. 
Let us find a solution of the problem (2) in the form (1). We need prove that the following integral converges, i.e.,
1 2 π i ϖ e λ n t ( I λ A ) 1 h d λ H , h H .
For this purpose consider a contour γ k : = { λ ϖ , | λ < R k , k N } , R k . Using simple estimating, then applying Lemmas 6 and 9 [1] respectively, we get
γ k e λ n t ( I λ A ) 1 h d λ H γ k e λ n t | · ( I λ A ) 1 h H | d λ C h H γ k e t Re λ n | d λ .
It is clear that
γ k e t Re λ n | d λ C , k .
The last fact gives us the desired result. Since A is bounded, then we have
W ˜ u ( t ) = W ˜ 1 2 π i ϖ e λ n t A ( I λ A ) 1 h d λ = W ˜ A 1 2 π i ϖ e λ n t ( I λ A ) 1 h d λ =
= 1 2 π i ϖ e λ n t ( I λ A ) 1 h d λ .
Combining the last relation with (3), we obtain u D ( W ˜ ) . Analogously to the above, using Lemmas 6 and 9 [1] respectively, we can show that the the following derivative exists, i.e.,
d u d t = 1 2 π i ϖ e λ n t λ n A ( I λ A ) 1 h d λ H .
Notice that λ n A n ( I λ A ) 1 = ( I λ A ) 1 ( I + λ A + + λ n 1 A n 1 ) , substituting this relation to the above formula, we obtain
A n 1 d u d t = 1 2 π i ϖ e λ n t ( I λ A ) 1 h d λ + 1 2 π i ϖ e λ n t k = 0 n 1 λ k A k h d λ .
The second integral equals to zero, by virtue of the fact that the function under the integral is analytical inside the intersection of the domain G with the circle of the arbitrary radius R , and it decreases sufficiently fast on the arch of the radius R , when R . Here we denote by G the interior of the contour ϖ . Thus, we have come to the relation
A n 1 d u d t = 1 2 π i ϖ e λ n t ( I λ A ) 1 h d λ .
Since the left-hand side of the last relation belongs to D ( W ˜ n 1 ) , we can claim that it is so for the right-hand side also. It follows that u D ( W n ) . Now, if we recall the expression for u , we get A n 1 u t + W ˜ u = 0 . Applying the operator W ˜ n 1 to both sides of the last relation, we obtain the fact that u is a solution of the Equation (2). Let us show that the initial condition holds in the sense u ( t ) H h , t + 0 . It becomes clear in the case h D ( W ˜ ) , for in this case it suffices to apply Lemmas 7 and 10 [1] respectively, what gives us the desired result, i.e., we can put u ( 0 ) = h . Consider a case when h is an arbitrary element of the Hilbert space H and let us involve the accretive property of the operator W ˜ n . Consider an operator
S t h = 1 2 π i ϖ e λ n t A ( I λ A ) 1 h d λ , t > 0 .
In accordance with the above, it is clear that S t : H H . Let us prove that S t H H 1 , t > 0 . Firstly, assume that h D ( W ˜ ) . Let us multiply the both sides of the relation (2) on u in the sense of the inner product, we get u t , u H + ( W ˜ n u , u ) H = 0 . Consider a real part of the last relation, we have Re u t , u H + Re ( W ˜ n u , u ) H = u t , u H / 2 + u , u t H / 2 + Re ( W ˜ n u , u ) H . Therefore ( u ( t ) H 2 ) t = 2 Re ( W ˜ n u , u ) H 0 . Integrating both sides, we get
u ( τ ) H 2 u ( 0 ) H 2 = 0 τ d d t u ( t ) H 2 d t 0 .
The last relation can be rewritten in the form S t h H h H , h D ( W ˜ ) . Since D ( W ˜ ) is a dense set in H , we obviously obtain the desired result, i.e., S t H H 1 . Now, having assumed that h n H h , n , { h n } D ( W ˜ ) , h H , consider the following reasonings u ( t ) h H = S t h h H = S t h S t h n + S t h n h n + h n h H S t · h h n H + S t h n h n H + h n h H . Note that S t h n H h n , t + 0 . It is clear that if we choose n so that h h n H < ε / 3 and after that choose t so that S t h n h n H < ε / 3 , then we obtain ε > 0 , δ ( ε ) : u ( t ) h H < ε , t < δ . Thus, we can put u ( 0 ) = h and claim that the initial condition holds in the case h H . The decomposition on the series of the root vectors (1) is given by virtue of Theorems 1–3 respectively. The uniqueness follows easily from the fact that W ˜ n is accretive. In this case, repeating the previous reasonings, we come to
ϕ ( τ ) H 2 ϕ ( 0 ) H 2 = 0 τ d d t ϕ ( t ) 2 d t 0 ,
where ϕ is a sum of solutions u 1 and u 2 . Notice that by virtue of the initial conditions, we have ϕ ( 0 ) = 0 , thus relation (4) can hold only if ϕ = 0 . The proof is complete. □
Remark 1.
Note that the assumption n > ρ , that is additional to the ones of the above theorems, guaranties an opportunity to choose a sequence of numbers { N ν } 0 so that in the used terms, we have
| λ N ν + k | λ N ν + k 1 e | λ N ν τ , 0 < τ < n ρ .
This fact follows from the results by Lidskii ([2], p. 23).

3.1. Evolution Equations with the Quasi-Polynomial Right-Hand Side

Let J : = ( a , b ) R , Ω : = [ 0 , ) , consider functions u ( t , x ) , t Ω , x J ¯ . In accordance with the said above, we will consider functional spaces with respect to the variable x and we will assume that if u belongs to a functional space then the fact holds for all values of the variable t . Consider a Cauchy problem
d u d t = k = s s Q k D a + k ϑ u = : P s , ϑ u , ϑ > 0 , u ( 0 , x ) = h ( x ) L 2 ( J ) ,
where Q k = const , at the right-hand side we have a linear combination of the Riemann-Liouville fractional integro-differential operators acting in L 2 ( J ) with respect to the variable x , see [32]
D a + ψ f ( x ) : = d d x [ ψ ] + 1 I a + 1 { ψ } φ ( x ) , ψ > 0 , D a + ψ : = I a + ψ ,
I a + ψ φ ( x ) : = 1 Γ ( ψ ) a x ( x t ) ψ 1 φ ( t ) d t , φ L 1 ( J ) , ψ > 0 .
we will call (analogously to the theory of ordinary differential equations) the right-hand side of Equation (5) quasi-polynomial. Consider a case when the right-hand side of the Equation (5) can be represented as follows
P s , ϑ u = ( η D 2 + ξ D a + β ) n u , 0 < β < 1 / n , η < 0 , ξ > 0 .
Let us show that under such assumptions the problem (5) has a unique solution which can be found due to a certain formula. Consider the operator W : = η D 2 + ξ D a + β , D ( W ) = C 0 ( J ) . Note that the hypotheses H 1 , H 2 hold regarding the operator, if we assume that H : = L 2 ( J ) , H + : = H 0 1 ( J ) , where L 2 ( J ) is a Lebesque space, H 0 k ( J ) , k N is a subspace of the Sobolev space H k ( J ) defined to be the closure of C 0 ( J ) . It follows from the strictly accretive property of the fractional differential operator (see [14]) and the estimate
D a + β f L 2 C f H 0 1 , f C 0 ( J ) , β ( 0 , 1 ) .
Let us prove that
C ( D 2 f , f ) L 2 Re ( W f , f ) L 2 C ( D 2 f , f ) L 2 , f C 0 ( J ) .
Using relation (7) and the Friedrichs inequality, we obtain Re ( D a + β f , f ) L 2 f L 2 D a + β f L 2 C f H 0 1 2 = C ( D 2 f , f ) L 2 , f C 0 ( J ) , what gives us the upper estimate. The lower estimate follows easily from the accretive property of the fractional differential operator of the order less than one. Using relation (8), the corollary of the minimax principle, we get λ j ( H ) λ j ( D 2 ) , where H is a real part of the operator W ˜ . Therefore, taking into account the well-known fact λ j ( D 2 ) = π 2 j 2 / ( b a ) 2 , we get λ j ( H ) j 2 , it follows that μ ( H ) > 1 . Now we can study the Cauchy problem (5) by restricting the one to the problem
d u d t = ( η D 2 + ξ D a + β ) n u , u ( 0 ) = h D ( W ˜ ) .
Note that in the case α = n = 1 the condition θ < π / 2 α is fulfilled by virtue of the conditions H 1 , H 2 . In the case n > 1 , it can be verified directly that the condition H 3 holds, we should apply (7) and the Cauchy-Schwartz inequality. In accordance with Theorem 4, we are able to present a solution of the problem (5), with the restricted assumptions regarding h , as follows
u ( t ) = 1 2 π i ϖ e λ n t A ( I λ A ) 1 h d λ ,
where the used terms relate to the operator W ˜ . Thus, we have in the reminder a question how to weaken conditions imposed upon the function h as well as whether the representation (6) holds. To answer the questions consider the following reasonings. Further, for the sake of the simplicity, we consider a case when η = 1 , ξ = 1 . This assumption does not restrict generality of reasonings. Let us show that D ( W ˜ ) H 0 2 ( J ) . Using H 2 , we have the implication
f k W f f k H 0 1 f , { f k } 1 C 0 ( J ) .
Applying (7), we get D a + β f k L 2 D a + β f . The following fact can be obtained easily, we have omitted the proof
{ f k W f , D a + β f k L 2 D a + β f } D 2 f k L 2 D 2 f .
Combining the above implications we obtain the desired result, i.e., D ( W ˜ ) H 0 2 ( J ) . Consider a set H 0 + s ( J ) : = { f : f H s ( J ) , f ( k ) ( a ) = 0 , k = 0 , 1 , , s 1 } , s N . It is clear that H 0 s ( J ) H 0 + s ( J ) , thus we can define the operator W + as the extension of the operator W ˜ on the set H 0 + 2 ( J ) , we have W ˜ W + . Let us show that D ( W + n ) = H 0 + 2 n ( J ) . Assume that f H 0 + 2 n ( J ) , then f D ( W + n ) , it can be verified directly. If f D ( W + n ) , then in accordance with the definition, we have W + n 1 f H 0 + 2 ( J ) . It follows that W + g 1 H 0 + 2 ( J ) , where g 1 = W + n 2 f H 0 + 2 ( J ) . Hence D 2 g 1 + D a + β g 1 H 0 + 2 ( J ) . Applying the operator I a + 2 to the both sides of the last relation, we easily get
g 1 + I a + 2 β g 1 H 0 + 4 ( J ) .
Using the constructive features of relation (10), we can conclude firstly g 1 H 0 + 3 ( J ) and due to the same reasonings establish the fact g 1 H 0 + 4 ( J ) secondly, what gives us W + n 2 f H 0 + 4 ( J ) . Using the absolutely analogous reasonings we prove that W + n k f H 0 + 2 k ( J ) , k = 1 , 2 , , n . Thus, we obtain the desired result. Let us show that
W + n f = k = 0 n ( 1 ) n k C n k D a + β k + 2 ( n k ) f , f H 0 + 2 n ( J ) .
We need establish the formula D a + β k D 2 ( n k ) f = D 2 ( n k ) D a + β k f = D a + β k + 2 ( n k ) f , f H 0 + 2 n ( J ) , ( k = 1 , 2 , , n ) for this purpose, in accordance with the Theorem 2.5 ([32], p. 46), we should prove that f I a + β k + 2 ( n k ) ( L 1 ) or f = I a + 2 ( n k ) φ a . e . , φ I a + β k ( L 1 ) . It is clear that almost everywhere, we have f = I a + 2 n D 2 n f = I a + 2 ( n k ) D 2 ( n k ) f = I a + 2 ( n k ) φ , where φ : = D 2 ( n k ) f . Note that the conditions of the Theorem 13.2 ([32], p. 229) hold, i.e., the Marchaud derivative of the function φ belongs to L 1 ( J ) . Hence φ I a + β k ( L 1 ) and we obtain the required formula. Using the well-known formulas for linear operators ( A + B ) C A C + B C , C ( A + B ) = C A + C B , applying the Leibniz formula, we obtain (11). Now, combining the obvious inclusion W ˜ n W + n with (11), we get
W ˜ n k = 0 n ( 1 ) n k C n k D a + β k + 2 ( n k ) ,
where C n k are binomial coefficients. The next question is whether the operator W ˜ n is accretive. By direct calculation, we have
Re ( W ˜ n f , f ) L 2 = k = 0 n C n k Re D a + β k + ( n k ) f , D n k f L 2 = k = 0 n C n k Re D a + β k g k , g k L 2 0 ,
where g k : = D n k f , f D ( W ˜ n ) . Note that the last inequality holds by virtue of the strictly accretive property of the fractional differential operator of the order less than one (see [14]). Thus, the uniqueness and the opportunity to weaken conditions imposed on h follow from Theorem 4. Here, we should remark that the last theorem gives us the fact that the existing solution is unique in the set D ( W ˜ n ) , we should recall the fact W ˜ n W + n . However, using the same method, we can establish the uniqueness of the solution of the problem (5). Having known the root vectors of the operator A , applying (1), we can represent the obtained solution as a series. Note that using the same reasonings, we can solve the problem (5), with the right-hand side P s , ϑ = ( η D 2 + ξ I a + β ) n . Using the ordinary properties of the homogenous equations, combining the obtained results, we can consider the problem (5) with the right-hand side of the mixed type, i.e., it contains fractional integrals as well as fractional derivatives. Here, we should remark that such technical peculiarities (we left them to the reader) do not lie in the scope of this paper devoted mostly to methods and their applications.

3.2. Kipriyanov Operator

Using notations of the paper [30] we assume that Ω is a convex domain of the m-dimensional Euclidean space E m , P is a fixed point of the boundary Ω ,   Q ( r , e ) is an arbitrary point of Ω ; we denote by e a unit vector having a direction from P to Q , denote by r = | P Q the Euclidean distance between the points P , Q , and use the shorthand notation T : = P + e t , t R . We consider the Lebesgue classes L p ( Ω ) , 1 p < of complex valued functions. For the function f L p ( Ω ) , we have
Ω | f ( Q ) p d Q = ω d χ 0 d ( e ) | f ( Q ) p r m 1 d r < ,
where d χ is an element of solid angle of the unit sphere surface (the unit sphere belongs to E m ) and ω is a surface of this sphere, d : = d ( e ) is the length of the segment of the ray going from the point P in the direction e within the domain Ω . Without loss of generality, we consider only those directions of e for which the inner integral on the right-hand side of the equality (12) exists and is finite. It is the well-known fact that these are almost all directions. Denote by D i f a distributional derivative of the function f with respect to a coordinate variable with index 1 i m . For convenience, we use the Einstein convention a i b i : = i = 1 m a i b i .
Here, we study the case β ( 0 , 1 ) . Assume that Ω E m is a convex domain, with a sufficient smooth boundary ( C 3 class) of the m-dimensional Euclidian space. For the sake of the simplicity, we consider that Ω is bounded, but the results can be extended to some type of unbounded domains. In accordance with the definition given in the paper [16], we consider the directional fractional integrals. By definition, put
( I 0 + β f ) ( Q ) : = 1 Γ ( β ) 0 r f ( P + t e ) ( r t ) 1 β t r m 1 d t , ( I d β f ) ( Q ) : = 1 Γ ( β ) r d f ( P + t e ) ( t r ) 1 β d t ,
f L p ( Ω ) , 1 p ,
where Γ ( β ) is the gamma function. Also, we consider auxiliary operators, the so-called truncated directional fractional derivatives (see [34]). By definition, put
( D 0 + , ε β f ) ( Q ) = β Γ ( 1 β ) 0 r ε f ( Q ) r m 1 f ( P + e t ) t m 1 ( r t ) β + 1 r m 1 d t + f ( Q ) Γ ( 1 β ) r β , ε r d , ( D 0 + , ε β f ) ( Q ) = f ( Q ) ε β , 0 r < ε ;
( D d , ε β f ) ( Q ) = β Γ ( 1 β ) r + ε d f ( Q ) f ( P + e t ) ( t r ) β + 1 d t + f ( Q ) Γ ( 1 β ) ( d r ) β , 0 r d ε , ( D d , ε β f ) ( Q ) = f ( Q ) β 1 ε β 1 ( d r ) β , d ε < r d .
Now, we can define the directional fractional derivatives as follows
D 0 + β f = lim ( L p ) ε 0 D 0 + , ε β f , D d β f = lim ( L p ) ε 0 D d , ε β f , 1 p .
The properties of these operators are described in detail in the papers [14,34]. We suppose I 0 + 0 = I . Nevertheless, this fact can be easily proved dy virtue of the reasonings corresponding to the one-dimensional case, and is given in [32]. We also consider integral operators with a weighted factor (see [32], p. 175) defined by the following formal construction
I 0 + β ϕ f ( Q ) : = 1 Γ ( β ) 0 r ( ϕ f ) ( P + t e ) ( r t ) 1 β t r m 1 d t ,
where ϕ is a real-valued function. Consider a linear combination of the uniformly elliptic operator, which is written in the divergence form, and a composition of a fractional integro-differential operator, where the fractional differential operator is understood as the adjoint operator regarding the Kipriyanov operator (see [14,30,31])
L : = T + I 0 + σ ϕ D d β , σ [ 0 , 1 ) , D ( L ) = H 2 ( Ω ) H 0 1 ( Ω ) ,
where T : = D j ( a i j D i · ) , i , j = 1 , 2 , , m , under the following assumptions regarding coefficients
a i j ( Q ) C 2 ( Ω ¯ ) , Re a i j ξ i ξ j γ a | ξ 2 , γ a > 0 , Im a i j = 0 ( m 2 ) , ϕ L ( Ω ) .
Note that in the one-dimensional case, the operator I 0 + σ ϕ D d β is reduced to a weighted fractional integro-differential operator composition, which was studied properly by many researchers [35,36,37,38], more detailed historical review see in ([32], p. 175). In accordance with the Theorem 9 [16], we claim that the conditions H 1 , H 2 are fulfilled, if γ a is sufficiently large in comparison with ϕ , where we put M : = C 0 ( Ω ) . Note that the order μ of the operator H can be evaluated easily through the order of the regular differential operator and since the latter can be found by the methods described in [19]. More precisely, we have
C ( Re T f , f ) H ( H f , f ) H C ( Re T f , f ) H , f C 0 ( Ω ) .
Applying the minimax principle, we get λ j ( H ) λ j ( Re T ) . Using the well-known formula for regular differential operators λ j ( Re T ) j 2 / m (see [19]), we get μ ( H ) = 2 / m . Note that the case m = 1 leads us to the previous paragraph, a special interest may appear by virtue of the fact that the composition of a fractional integro-differential operator has been involved. In the case m 2 , we are forced to assume that n m . The last assumption is based on the opportunity to fulfill the condition n = 2 / μ corresponding to Theorem 1 (since ρ 2 / μ , see the Theorem 3 [16]) and condition n > 2 / μ corresponding to Theorem 2. Since the hypothesis H 1 , H 2 guarantee the fulfilment of the condition Θ ( A ) L 0 ( θ ) , θ < π / 2 n in the case only n = 1 , then to consider higher powers, we should verify the fulfilment of the condition H 3 . Firstly, we should assume that a i j = a j i . Secondly, we should use relation (27) [16] in accordance with which, we have
I 0 + σ ϕ D d α f , g L 2 C f H 0 1 g L 2 , f , g C 0 ( Ω ) .
Thus, the additional condition H 3 holds. Then in accordance with Theorem 4, we can claim that there exists a solution of the problem (2), where W is a restriction of L on the set C 0 ( Ω ) , the coefficients (13) are sufficiently smooth to guarantee the fact the right-hand side of (2) has a sense. Note that the solvability of the uniqueness problem as well as the opportunity to extend the initial condition depends on the accretive property of the operator W ˜ n . The last problem can be studied by the methods similar to the ones used in the previous paragraph. Indeed, we have established the accretive property in the one-dimensional case.

3.3. Riesz Potential

Consider a space L 2 ( Ω ) , Ω : = ( , ) . We denote by H 0 2 , ς ( Ω ) the completion of the set C 0 ( Ω ) with the norm
f H 0 2 , ς ( Ω ) = f L 2 ( Ω ) 2 + D 2 f L 2 ( Ω , ω ς ) 2 1 / 2 , ς R ,
where
L 2 ( Ω , ω ς ) : = f : f L 2 ( Ω , ω ς ) < , f L 2 ( Ω , ω ς ) : = Ω | f ( t ) 2 ω ς ( t ) d t 1 / 2 , ω ( x ) : = 1 + | x ,
the above integral is understood in the Lebesque sense. Let us notice the following fact (see Theorem 1 [39]), if ς > 4 , then H 0 2 , ς ( Ω ) L 2 ( Ω ) . Consider a Riesz potential
I β f ( x ) = B β f ( s ) | s x β 1 d s , B β = 1 2 Γ ( β ) cos ( β π / 2 ) , β ( 0 , 1 ) ,
where f is in L p ( Ω ) , 1 p < 1 / β . It is obvious that I β f = B β Γ ( β ) ( I + β f + I β f ) , where
I ± β f ( x ) = 1 Γ ( β ) 0 f ( s x ) s β 1 d s ,
these operators are known as fractional integrals on the whole real axis (see [32], p. 94). Assume that the following condition holds σ / 2 + 3 / 4 < β < 1 , where σ is a non-negative constant. Following the idea of the monograph ([32], p. 176), consider a sum of a differential operator and a composition of fractional integro-differential operators
W : = D 2 a D 2 + I + σ ξ I 2 ( 1 β ) D 2 + δ I , D ( W ) = C 0 ( Ω ) ,
where ξ ( x ) L ( Ω ) , a ( x ) L ( Ω ) C 2 ( Ω ) , Re a ( x ) > γ a ( 1 + | x ) 5 . Let Ω : = [ 0 , ) , consider the functions u ( t , x ) , t Ω , x Ω . Similarly to the previous paragraph, we will consider functional spaces with respect to the variable x and we will assume that if u belongs to a functional space then this fact holds for all values of the variable t , wherewith all standard topological properties of a space as completeness, compactness etc. remain correctly defined. Consider a Cauchy problem (2) in the corresponding terms. Notice that in accordance with the results [16], we can claim that if min { γ a , δ } is sufficiently large, then the hypothesis H 1 , H 2 hold regarding: the operator W ˜ , the set C 0 ( Ω ) , the spaces L 2 ( Ω ) , H 0 2 , 5 ( Ω ) , more precisely we should put H : = L 2 ( Ω ) , H + : = H 0 2 , 5 ( Ω ) , M : = C 0 ( Ω ) . Thus, in accordance with the condition H 2 , we have
( H f , f ) L 2 = Re ( W f , f ) L 2 γ a f H 0 2 , 5 2 = C ( D 2 w D 2 f , f ) L 2 + C ( f , f ) L 2 , f C 0 ( Ω ) .
where w ( x ) = ( 1 + | x | ) 5 . Let us consider the operator B = D 2 w D 2 + δ I , D ( B ) = C 0 ( Ω ) . It is clear that by virtue of the minimax principle, we can estimate the eigenvalues of the operator W ˜ via estimating the eigenvalues of the operator B ˜ . Hence, we have come to the problem of estimating the eigenvalues of the singular operator. Here, we should point out that there exists the Fefferman concept that covers such a kind of problems. For instance, the Rozenblyum result is presented in the monograph ([19], p. 47), in accordance with which we can choose such an unbounded subset of R that the relation λ j ( B ˜ ) j 4 holds. Thus, we left this question to the reader for a more detailed study and reasonably allow ourselves to assume that the condition μ ( H ) = 4 holds. Note that in the case n = 1 , in accordance with the Theorem 3 [16], Theorem 1, Theorem 4, we are able to present a solution of the problem (2) in the form (1). Now, assume additionally that Im a = 0 , then D 2 a D 2 is selfadjoint. It follows that the operator W ˜ is selfadjoint. Applying Formula (40) [16], we get
I + σ ξ I 2 ( 1 α ) D 2 f L 2 C f H 0 2 , 5 , f H 0 2 , 5 ( Ω ) .
Using this fact, applying the Cauchy-Schwartz inequality, we conclude that the condition H 3 holds. Moreover, since W ˜ is selfadjoint, then we can easily prove that Re ( W ˜ n f , f ) H 0 , n = 1 , 2 , . Therefore, applying Theorem 4, we can claim that there exists a solution of the problem (2), where the coefficients of the operator W are sufficiently smooth to guaranty the fact the right-hand side of the Formula (2) has a sense. Moreover, we can assume that h H and claim that the existing solution is unique.

3.4. Difference Operator

The approach implemented in studying the difference operator is remarkable due to the appeared opportunity to set the problem within the framework of the created theory having constructed a suitable perturbation of the operator composition. Consider a space L 2 ( Ω ) , Ω : = ( , ) , define a family of operators
T t f ( x ) : = e c t k = 0 ( c t ) k k ! f ( x d μ ) , f L 2 ( Ω ) , c , d > 0 , t 0 ,
where convergence is understood in the sense of L 2 ( Ω ) norm. In accordance with the Lemma 6 [16], we know that T t is a C 0 semigroup of contractions, the corresponding infinitesimal generator and its adjoint operator are defined by the following expressions
Y f ( x ) = c [ f ( x ) f ( x d ) ] , Y f ( x ) = c [ f ( x ) f ( x + d ) ] , f L 2 ( Ω ) .
Let us find a representation for fractional powers of the operator Y . Using Formula (45) [16], we get
Y β f = k = 0 M k f ( x k d ) , f L 2 ( Ω ) , M k = β Γ ( k β ) k ! Γ ( 1 β ) c β , β ( 0 , 1 ) .
we need the following theorem (see the Theorem 14 [16]).
Theorem 5.
Assume that Q is a closed operator acting in L 2 ( Ω ) , Q 1 S ( L 2 ) , the operator N is strictly accretive, bounded, R ( Q ) D ( N ) . Then a perturbation
L : = Y a Y + b Y β + Q N Q , a , b L ( Ω ) ,
satisfies conditions H1 and H2, if γ N > σ Q 1 2 , where we put M : = D 0 ( Q ) ,
σ = 4 c a L + b L β c β Γ ( 1 β ) k = 0 Γ ( k β ) k ! .
Observe that by virtue of the made assumptions regarding Q , we have H Q L 2 ( Ω ) . We have chosen the space L 2 ( Ω ) as a space H and the space H Q as a space H + . Applying the condition H 2 , we get
C ( Q N Q f , f ) H ( H f , f ) H C ( Q N Q f , f ) H , f D 0 ( Q ) ,
where H is a real part of W ˜ . Therefore, by virtue of the minimax principle, we get λ j ( H ) λ j ( Q N Q ) . Hence μ ( H ) = μ ( Q N Q ) . Thus, we have naturally come to the significance of the operator Q and the remarkable fact that we can fulfill the conditions of Theorem 4 choosing the operator Q in the artificial way. Applying Theorem 4, we can claim that there exists a solution of the problem (2), where W is a restriction of L on the set M (see Introduction), functions a , b are sufficiently smooth to guaranty the fact the right-hand side of (2) has a sense. The extension of the initial conditions on the whole space H , as well as solvability of the uniqueness problem can be implemented in the case when the operator W ˜ n is accretive. In its own turn, it is clear that the particular methods, to establish the accretive property, can differ and may depend on the concrete form of the operator Q .

3.5. Artificially Constructed Normal Operator

In this paragraph we consider an operator class which cannot be completely studied by methods [2], at the same time Theorem 3 gives us a rather relevant result. Our aim is to construct a normal operator N being satisfied the Theorem 3 conditions, such that N S ˜ α \ S α . Let us consider the following example as a perquisite for the further reasonings.
Example 1.
Here we would like to produce an example of the sequence { μ n } 1 that satisfies the condition ( ln κ + 1 x ) μ n = o ( n κ ) , ( 0 < κ < 1 ) , and at the same time
n = 1 1 | μ n 1 / κ = .
Consider a sequence μ n = n κ ln κ ( n + q ) · ln κ ln ( n + q ) , q > e e 1 , then using the integral test for convergence we can easily see that the previous series diverges. At the same time substituting, we get easily
ln κ μ n μ n C ln κ ( n + q ) n κ ln κ ( n + q ) · ln κ ln ( n + q ) = C n κ ln κ ln ( n + q ) , n = 1 , 2 , ,
what gives us the fulfilment of the first condition.
Consider the abstract separable Hilbert space H and an operator N acting in the space as follows
N f = n = 1 λ n f n e n , f n = ( f , e n ) H , λ n = μ n + i η n ,
where { e n } 1 H is an orthonormal basis, the sequence { μ n } 1 is defined in Example 1, the sequence { η n } 1 is satisfied the following condition | η n | λ n 1 / 2 , n = 1 , 2 , . Define the space H + as follows
H + : = f H : f H + 2 : = n = 1 | λ n | f n 2 < .
It is clear that H + is dense in H , since { e n } 1 H + . Let us show that embedding of the spaces H + H is compact. Consider the operator B : H H defined as follows
B f = n = 1 | λ n 1 / 2 f n e n .
The fact that B is a compact operator can be proved easily due to the well-known criterion of compactness in the Banach space endowed with a basis (we left the prove to the reader). Notice that if f H + , then g H , where g is defined by its Fourier coefficients g n = | λ n 1 / 2 | f n . By virtue of such a correspondence, we can consider any bounded set in the space H + as a bounded set in the space H . Applying the operator B to the element g , we get the element f . Since B is a compact operator, then we conclude that the image of the bounded set in the sense of the norm H + is a compact set in the sense of the norm H . Define the set M as a linear manifold generated by the basis vectors. Thus, we have obtained the relation H + H and established the fulfilment of the condition H 1 . The first relation of the condition H 2 can be obtained easily due to the application of the Cauchy-Schwarz inequality. To obtain the second one, we should note that η n 2 + μ n 2 | λ n + μ n 2 ; μ n | λ n | ( 1 λ n 1 ) 1 / 2 M | λ n | , M = ( 1 μ 1 1 ) 1 / 2 , n = 1 , 2 , . Therefore
Re ( N f , f ) H = n = 1 Re λ n | f n 2 M n = 1 | λ n | f n 2 = M f H + 2 .
Now, we conclude that hypotheses H 1 , H 2 hold. Let us show that the condition H 3 is satisfied, we have
Im ( N f , g ) H = n = 1 Im λ n f n g ¯ n n = 1 | λ n 1 / 2 | f n | g n f H + g H ,
thus we have obtained the desired result. Consider a Cauchy problem
d u d t = N 1 / κ u , u ( 0 ) = h D ( N ) , κ = 1 / 2 , 1 / 3 , ,
where h is supposed to be an arbitrary element, if the operator N 1 / κ is accretive. In accordance with Theorem 4, we conclude that there exists a solution of problem (14) presented by the series (1). Let us establish conditions under which being imposed the right-hand side of (14) is an accretive operator. Note that using the relation between the real and imaginary parts of an eigenvalue, we get μ n ( 1 | λ n 1 ) 1 / 2 | λ n | ( 1 λ n 1 ) 1 / 2 η n 2 ( 1 | μ 1 1 ) 1 / 2 η n 2 > ( 1 e κ ) 1 / 2 η n 2 . Hence, the eigenvalues lie in a parabolic domain W : = { z C , Re z > ( 1 e κ ) 1 / 2 Im 2 z } . Notice that a condition
| arg λ n π κ / 2 , n = 1 , 2 ,
guarantees the following fact
Re λ n 1 / κ = | λ n 1 / κ cos arg λ n κ 0 , n = 1 , 2 ,
from what follows the desired accretive property
Re ( N 1 / κ f , f ) H = n = 1 Re λ n 1 / κ | f n 2 0 .
It is clear that the condition (15) holds for the eigenvalues with sufficiently large numbers of indexes, since they lie in the parabolic domain W . Thus, we see that a finite number of eigenvalues do not satisfy the condition (15). Here, we should note that using simple reasonings (we left them to the reader), we can find the initial number, starting from which the condition (15) holds. Now, it is clear that if we additionally assume that | η n / μ n tan ( π κ / 2 ) , n = 1 , 2 , , N , where N is a sufficiently large number, then we obtain (16). Therefore, we can assume that h H and claim that the existing solution is unique. The constructed normal operator indicates the significance of the made in Theorem 3 clarification of the results [2]. However, we produce one more relevant application of the mentioned theorem in the following paragraph.

3.6. Evolution Equations with the Fractional Derivative at the Left-Hand Side

In this paragraph, we still consider a Hilbert space H consists of element-functions u : R + H , u : = u ( t ) , t 0 assuming that if u belongs to H then the fact holds for all values of the variable t . We understand such operations as differentiation and integration in the generalized sense that is caused by the topology of the Hilbert space H . The derivative is understood as a limit
u ( t + Δ t ) u ( t ) Δ t H d u d t , Δ t 0 .
Let t J : = [ a , b ] , 0 < a < b < . The following integral is understood in the Riemann sense as a limit of partial sums
i = 0 n u ( ξ i ) Δ t i H J u ( t ) d t , ζ 0 ,
where ( a = t 0 < t 1 < < t n = b ) is an arbitrary splitting of the segment J , ζ : = max i ( t i + 1 t i ) , ξ i is an arbitrary point belonging to [ t i , t i + 1 ] . The sufficient condition of the last integral existence is a continuous property (see [33], p. 248), i.e. u ( t ) H u ( t 0 ) , t t 0 , t 0 J . The improper integral is understood as a limit
a b u ( t ) d t H a c u ( t ) d t , b c , c [ 0 , ] .
Combining the operations we can consider a generalized fractional derivative in the Riemann-Liouville sense (see [32]), in the formal form, we have
D 1 / α f ( t ) : = 1 Γ ( 1 1 / α ) d d t 0 f ( t + x ) x 1 / α d x , α > 1 .
Let us study a Cauchy problem
D 1 / α u = W ˜ u , u ( 0 ) = h D ( W ˜ ) ,
where in the case when the operator composition D 1 1 / α W ˜ is accretive, we assume that h H .
Theorem 6.
Assume that the Theorem 3 conditions hold, then there exists a solution of the Cauchy problem (17) in the form (1), a sequence of natural numbers { N ν } 0 can be chosen in accordance with the claim of Theorem 3. Moreover, the existing solution is unique if the operator composition D 1 1 / α W ˜ is accretive.
Proof. 
Let us find a solution of problem (17) in the form (1). Analogously to the reasonings of Theorem 4, using Lemmas 6 and 9 [1] respectively, it is not hard to prove that the following integrals exist, i.e.
1 2 π i ϖ e λ α t ( I λ A ) 1 h d λ H ; d u d t = 1 2 π i ϖ e λ α t λ α A ( I λ A ) 1 h d λ H .
Note that the first relation gives us the fact u ( t ) D ( W ˜ ) . Applying the scheme of the proof corresponding to the ordinary integral calculus, using the contour γ k (see reasonings of Theorem 4), applying Lemmas 6 and 9 [1] respectively, we can establish a formula
0 x σ d x ϖ e λ ζ ( t + x ) λ m A ( I λ A ) 1 h d λ = ϖ e λ ζ t λ m A ( I λ A ) 1 h d λ 0 x σ e λ ζ x d x ,
where σ ( 0 , 1 ) , ζ > 0 , m N 0 . In the same way, we get
d d t ϖ λ 1 α e λ α t A ( I λ A ) 1 h d λ = ϖ λ e λ α t A ( I λ A ) 1 h d λ .
Applying the obtained formulas, taking into account a relation
0 x 1 / α e λ α x d x = Γ ( 1 1 / α ) λ 1 α ,
we get
D 1 / α u = 1 2 π i ϖ e λ α t λ A ( I λ A ) 1 h d λ .
Making the substitution using the formula λ A ( I λ A ) 1 = ( I λ A ) 1 E , we obtain
D 1 / α u = 1 2 π i ϖ e λ α t ( I λ A ) 1 h d λ 1 2 π i ϖ e λ α t h d λ = I 1 + I 2 .
The second integral equals zero by virtue of the fact that the function under the integral is analytical inside the intersection of the domain G with the circle of an arbitrary radius R and it decreases sufficiently fast on the arch of the radius R , when R , here we denote by G the interior of the contour ϖ . Now, if we consider the expression for u , we obtain the fact that u is a solution of the equation i.e. D 1 / α u = W ˜ u . By virtue of Theorem 3, we obtain the decomposition on the series of the root vectors (1). Let us show that the initial condition holds in the sense u ( t ) H h , t + 0 . If h D ( W ˜ ) , then it becomes obvious, for in this case it suffices to apply Lemmas 7 and 10 [1] respectively, what gives us the desired result. To establish the fact in the case h H , we should involve the accretive property of the operator composition D 1 1 / α W ˜ . Following the scheme of the Theorem 4 reasonings, consider an operator
S t h = 1 2 π i ϖ e λ α t A ( I λ A ) 1 h d λ , t > 0 .
Using the first relation (18), the fact A B ( H ) , we can easily prove that S t : H H . Let us prove that S t H H 1 , t > 0 . We need to establish the following formula
D 1 1 / α D 1 / α u = u .
Using relation (19), we get
0 x 1 / α 1 d x ϖ e λ α ( t + x ) λ A ( I λ A ) 1 h d λ = ϖ e λ α t λ A ( I λ A ) 1 h d λ 0 x 1 / α 1 e λ α x d x =
= Γ ( 1 / α ) ϖ e λ α t A ( I λ A ) 1 h d λ = 2 π i Γ ( 1 / α ) u ( t ) .
Substituting (20) to the first term of the last relation, using the given above definition of the generalized fractional derivative, we obtain (21). Now, assume that h D ( W ˜ ) . Applying the operator D 1 1 / α to both sides of relation (17), using Formula (21), we get u + D 1 1 / α W ˜ u = 0 . Let us multiply the both sides of the last relation on u in the sense of the inner product, we get u t , u H + ( D 1 1 / α W ˜ u , u ) H = 0 . Consider a real part of the last relation, we have Re u t , u H + Re ( D 1 1 / α W ˜ u , u ) H = u t , u H / 2 + u , u t H / 2 + Re ( D 1 1 / α W ˜ u , u ) H . Therefore
u ( t ) H 2 t = 2 Re ( D 1 1 / α W ˜ u , u ) H 0 .
Integrating both sides, we get
u ( τ ) H 2 u ( 0 ) H 2 = 0 τ d d t u ( t ) H 2 d t 0 .
The last relation can be rewritten in the form S t h H h H , h D ( W ˜ ) . Since D ( W ˜ ) is a dense set in H , then we obviously obtain the desired result, i.e., S t H H 1 . The rest part of the proof is analogous to the proof of Theorem 4. The decomposition on the series of the root vectors (1) is given by virtue of Theorem 3. The uniqueness follows from the accretive property of the operator composition D 1 1 / α W ˜ .
Remark 2.
Eventually, we want to note that Theorem 6 is applicable to the operators studied in the previous paragraphs. The question that may appear is related to the description of the operator class for which the considered above operator composition is accretive.

4. Conclusions

In this paper, we studied Cauchy problems for the evolution equation in the abstract Hilbert space. The made approach allows us to obtain a solution analytically for the right-hand side belonging to a sufficiently wide class of operators. In this regard, such operators as the Riemann-Liouville fractional differential operator, the Kipriyanov operator, the Riesz potential, the difference operator have been involved. Moreover, we produced the artificially constructed normal operator for which the clarification of the Lidskii results relevantly works. It is remarkable that Theorem 6 covers many previously obtained results in the framework of the theory of factional differential equations. However, the main advantage is the obtained abstract formula for the solution. Moreover, the norm-convergence of the series representing the solution allows us to apply the methods of the approximation theory. It is also worth noticing that Theorem 6 jointly with Theorem 3 give us the opportunity to minimize conditions imposed on the fractional order at the left-hand side of the equation, and the artificially constructed normal operator provides the relevance and significance of such an achievement. Here, we should admit that the analogs of Theorem 6 can be obtained by Lidskii methods [2], but the used in the paper [1] variant of the natural lacunae method allows us to formulate the optimal conditions in comparison with the Lidskii results.
The interesting problem that may appear is how to consider the evolution equations of the arbitrary real order. We should note that the applied technique can be relevant in the case of the integer order higher than one, for the operator function under the integral presenting the solution is sufficiently “good”, what allows us to differentiate it. The latter opportunity becomes more valuable due to the appeared formula connecting the derivative of the integer order and the corresponding power of the operator. This creates a prerequisite to consider a generalized Taylor series and raise a question on convergence at least and representing the solution at most. It is remarkable that the supposed correspondence between derivatives and the operator powers leads us to the study of series of unbounded operators. Apparently, the hypotheses on the representing the solution by the Taylor series is equivalent to existing the exponent of an unbounded operator. We hope that this idea creates a prerequisite for further study in the direction.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Dedicated to the blessed memory of my teacher, Isai I. Mikaelyan.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kukushkin, M.V. Natural lacunae method and Schatten-von Neumann classes of the convergence exponent. arXiv 2021, arXiv:2112.10396v1. [Google Scholar]
  2. Lidskii, V.B. Summability of series in terms of the principal vectors of non-selfadjoint operators. Tr. Mosk. Mat. Obs. 1962, 11, 3–35. [Google Scholar]
  3. Haq, A.; Sukavanam, N. Existence and approximate controllability of Riemann-Liouville fractional integrodifferential systems with damping. Chaos Solitons Fractals 2020, 139, 110043. [Google Scholar] [CrossRef]
  4. Shukla, A.; Sukavanam, N.; Pandey, D.N. Approximate controllability of semilinear fractional control systems of order α ∈ (1,2]. J. Dyn. Control Syst. 2017, 23, 679–691. [Google Scholar] [CrossRef]
  5. Shukla, A.; Sukavanam, N.; Pandey, D.N. Approximate controllability of semilinear fractional stochastic control system. Asian-Eur. J. Math. 2018, 11, 1850088. [Google Scholar] [CrossRef]
  6. Haq, A. Partial-approximate controllability of semi-linear systems involving two Riemann-Liouville fractional derivatives. Chaos Solitons Fractals 2022, 157, 111923. [Google Scholar] [CrossRef]
  7. Singh, A.; Shukla, A.; Vijayakumar, V.; Udhayakumar, R. Asymptotic stability of fractional order (1, 2] stochastic delay differential equations in Banach spaces. Chaos Solitons Fractals 2021, 150, 111095. [Google Scholar] [CrossRef]
  8. Shkalikov, A.A. Perturbations of selfadjoint and normal operators with a discrete spectrum. Russ. Math. Surv. 2016, 71, 113–174. [Google Scholar] [CrossRef]
  9. Kukushkin, M.V. On one Method of Studying Spectral Properties of Non-selfadjoint Operators. Abstr. Appl. Anal. Hindawi Lond. UK 2020, 2020, 1461647. [Google Scholar] [CrossRef]
  10. Katsnelson, V.E. Conditions under which systems of eigenvectors of some classes of operators form a basis. Funct. Anal. Appl. 1967, 1, 122–132. [Google Scholar] [CrossRef]
  11. Krein, M.G. Criteria for completeness of the system of root vectors of a dissipative operator. Am. Math. Soc. Transl. Ser. Am. Math. Soc. Provid. RI 1963, 26, 221–229. [Google Scholar]
  12. Markus, A.S.; Matsaev, V.I. Operators generated by sesquilinear forms and their spectral asymptotics. Linear Oper. Integral Equ. Mat. Issled. Stiintsa Kishinev 1981, 61, 86–103. [Google Scholar]
  13. Markus, A.S. Expansion in root vectors of a slightly perturbed selfadjoint operator. Sov. Math. Dokl. 1962, 3, 104–108. [Google Scholar]
  14. Kukushkin, M.V. Asymptotics of eigenvalues for differential operators of fractional order. Fract. Calc. Appl. Anal. 2019, 22, 658–681. [Google Scholar] [CrossRef] [Green Version]
  15. Agranovich, M.S. On series with respect to root vectors of operators associated with forms having symmetric principal part. Funct. Anal. Its Appl. 1994, 28, 151–167. [Google Scholar] [CrossRef]
  16. Kukushkin, M.V. Abstract fractional calculus for m-accretive operators. Int. J. Appl. Math. 2021, 34, 1–41. [Google Scholar] [CrossRef]
  17. Nakhushev, A.M. The Sturm-Liouville problem for an ordinary differential equation of the second order with fractional derivatives in lower terms. Proc. Acad. Sciences. USSR 1977, 234, 308–311. [Google Scholar]
  18. Aleroev, T.S. On eigenvalues of one class of non-selfadjoint operators. Differ. Equ. 1994, 30, 169–171. [Google Scholar]
  19. Rozenblyum, G.V.; Solomyak, M.Z.; Shubin, M.A. Spectral theory of differential operators. Results Sci. Technol. Ser. Mod. Probl. Math. Dir. 1989, 64, 5–242. [Google Scholar]
  20. Bazhlekova, E. The abstract Cauchy problem for the fractional evolution equation. Fract. Calc. Appl. Anal. 1998, 1, 255–270. [Google Scholar]
  21. Clément, P.; Gripenberg, G.; Londen, S.-O. Hölder regularity for a linear fractional evolution equation. In Topics in Nonlinear Analysis: The Herbert Amann Anniversary Volume; Birkhäuser: Basel, Switzerland, 1998; pp. 69–82. [Google Scholar]
  22. Moroz, L.; Maslovskaya, A.G. Hybrid stochastic fractal-based approach to modeling the switching kinetics of ferroelectrics in the injection mode. Math. Model. Comput. Simulations 2020, 12, 348–356. [Google Scholar] [CrossRef]
  23. Mainardi, F. The fundamental solutions for the fractional diffusion-wave equation. Appl. Math. Lett. 1966, 9, 23–28. [Google Scholar] [CrossRef] [Green Version]
  24. Mamchuev, M.O. Boundary value problem for the time-fractional telegraph equation with Caputo derivatives Mathematical Modelling of Natural Phenomena. Spec. Funct. Anal. PDEs 2017, 12, 82–94. [Google Scholar] [CrossRef] [Green Version]
  25. Mamchuev, M.O. Solutions of the main boundary value problems for the time-fractional telegraph equation by the Green function method. Fract. Calc. Appl. Anal. 2017, 20, 190–211. [Google Scholar] [CrossRef] [Green Version]
  26. Pskhu, A.V. The fundamental solution of a diffusion-wave equation of fractional order. Izv. Math. 2009, 73, 351–392. [Google Scholar] [CrossRef]
  27. Wyss, W. The fractional diffusion equation. J. Math. Phys. 1986, 27, 2782–2785. [Google Scholar] [CrossRef]
  28. Kato, T. Perturbation Theory for Linear Operators; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1980. [Google Scholar]
  29. Gohberg, I.C.; Krein, M.G. Introduction to the Theory of Linear Non-Selfadjoint Operators in a Hilbert Space; Nauka, Fizmatlit: Moscow, Russia, 1965. [Google Scholar]
  30. Kipriyanov, I.A. On spaces of fractionally differentiable functions. Proc. Acad. Sci. USSR 1960, 24, 665–882. [Google Scholar]
  31. Kipriyanov, I.A. The operator of fractional differentiation and powers of the elliptic operators. Proc. Acad. Sci. USSR 1960, 131, 238–241. [Google Scholar]
  32. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science: Philadelphia, PA, USA, 1993. [Google Scholar]
  33. Krasnoselskii, M.A.; Zabreiko, P.P.; Pustylnik, E.I.; Sobolevskii, P.E. Integral Operators in the Spaces of Summable Functions; Science, FIZMATLIT: Moscow, Russia, 1966. [Google Scholar]
  34. Kukushkin, M.V. Spectral properties of fractional differentiation operators. Electron. J. Differ. Equ. 2018, 2018, 1–24. [Google Scholar]
  35. Dimovski, I.H.; Kiryakova, V.S. Transmutations, convolutions and fractional powers of Bessel-type operators via Maijer’s G-function. In Proceedings of the Complex Analysis and Applications, Varna, Bulgaria, 13 September 1983; Publishing House Bulgarian Academy of Sciences: Sofia, Bulgaria, 1986; pp. 45–66. [Google Scholar]
  36. Erdelyi, A. Fractional integrals of generalized functions. J. Austral. Math. Soc. 1972, 14, 30–37. [Google Scholar] [CrossRef] [Green Version]
  37. McBride, A. A note of the index laws of fractional calculus. J. Austral. Math. Soc. A 1983, 34, 356–363. [Google Scholar] [CrossRef] [Green Version]
  38. Nakhushev, A.M. Fractional Calculus and Its Application; Fizmatlit: Moscow, Russia, 2003. [Google Scholar]
  39. Adams, R.A. Compact lmbeddings of Weighted Sobolev Spaces on Unbounded Domains. J. Differ. Equ. 1971, 9, 325–334. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kukushkin, M.V. Evolution Equations in Hilbert Spaces via the Lacunae Method. Fractal Fract. 2022, 6, 229. https://doi.org/10.3390/fractalfract6050229

AMA Style

Kukushkin MV. Evolution Equations in Hilbert Spaces via the Lacunae Method. Fractal and Fractional. 2022; 6(5):229. https://doi.org/10.3390/fractalfract6050229

Chicago/Turabian Style

Kukushkin, Maksim V. 2022. "Evolution Equations in Hilbert Spaces via the Lacunae Method" Fractal and Fractional 6, no. 5: 229. https://doi.org/10.3390/fractalfract6050229

APA Style

Kukushkin, M. V. (2022). Evolution Equations in Hilbert Spaces via the Lacunae Method. Fractal and Fractional, 6(5), 229. https://doi.org/10.3390/fractalfract6050229

Article Metrics

Back to TopTop