\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2018 (2018), No. 29, pp. 1--24.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2018 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2018/29\hfil Fractional differentiation operators]
{Spectral properties of fractional differentiation operators}

\author[M. V. Kukushkin \hfil EJDE-2018/29\hfilneg]
{Maksim V. Kukushkin}

\address{Maksim V. Kukushkin \newline
International Committee ``Continental'',
Geleznovodsk 357401, Russia}
\email{kukushkinmv@rambler.ru}

\dedicatory{Communicated by Ludmila S. Pulkina}

\thanks{Submitted October 10, 2017. Published January 29, 2018.}
\subjclass[2010]{47F05, 47F99, 46C05}
\keywords{Fractional derivative; fractional integral;  energetic space;
\hfill\break\indent sectorial operator; strong accretive operator}

\begin{abstract}
 We consider the fractional differentiation operators in a variety of
 senses. We show that the strong accretive property is the common property
 of fractional differentiation operators. Also we prove that the sectorial
 property holds for operators second order with fractional derivative in
 lower terms, we explore the location of spectrum and  resolvent sets and
 show that the generalized  spectrum is discrete. We prove that there is
 two-sided estimate for eigenvalues of real component of operators second
 order with fractional derivative in lower terms.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks

\section{Introduction}

The term accretive applicable to a linear operator $T$ acting in a Hilbert
space $H$ was introduced by Friedrichs in the work \cite{fridrichs1958},
and means that the operator has the following property:
The numeric domain of values $\Theta(T)$ (see \cite[p.335]{kato1966})
is a subset of the right half-plane i.e.
    $$
    \operatorname{Re}\langle Tu,u\rangle_{H}\geq0,\quad u\in \mathfrak{D}(T).
    $$
Accepting  a  notation   \cite{kipriyanov1960} we assume that $\Omega$
is convex domain of the $n$  dimensional Euclidean space $\mathbb{E}^{n}$,
$P$ is a fixed point of the boundary $\partial\Omega$,
$Q(r,\vec{e})$ is an arbitrary point of $\Omega;$ we denote by
$\vec{e}$ as a unit vector having the direction from $P$ to $Q$,
denote by $r=|P-Q|$ as a Euclidean distance between points $P$ and $Q$.
We will consider classes of Lebesgue $L_p(\Omega)$, $1\leq p<\infty $
complex valued functions. In polar coordinates, the summability $f$ on $\Omega$
of degree $p$ means that
\begin{equation}\label{1}
\int_{\Omega}|f(Q)|^pdQ=\int_{\omega}d\chi\int_0^{d(\vec{e})}
|f(Q)|^pr^{n-1}dr<\infty,
\end{equation}
where $d\chi$  is the element of the solid angle
the surface of a unit sphere in $\mathbb{E}^{n}$ and $\omega$
 is a  surface of this sphere,   $d:=d(\vec{e})$
 is the length of  segment of  ray going from point $P$ in the direction
$\vec{e}$ within the domain $\Omega$.
Without lose of   generality, we consider only those directions of
$\vec{e}$ for which the inner integral on the right side of equality
\eqref{1} exists and is finite, is  well known that  this is almost all directions.
Notation  $\operatorname{Lip} \mu$, $0<\mu\leq1 $ means the set of functions
satisfying the Holder-Lipschitz condition
$$
\operatorname{Lip} \mu:=\big\{\rho(Q):|\rho(Q)-\rho(P)|\leq M r^{\mu},\;
P,Q\in \bar{\Omega}\big\}.
$$
The operator of fractional differentiation in the  sense of Kipriyanov
 defined in \cite{1kipriyanov1960}  by formal expression
\begin{equation*}
\mathfrak{D}^{\alpha}(Q)=\frac{\alpha}{\Gamma(1-\alpha)}
\int_0^{r} \frac{[f(Q)-f(P+\vec{e}t)]}{(r - t)^{\alpha+1}}
\big(\frac{t}{r}big) ^{n-1} dt+
C^{(\alpha)}_n f(Q) r ^{ -\alpha},\quad P\in\partial\Omega,
\end{equation*}
where
$C^{(\alpha)}_n = (n-1)!/\Gamma(n-\alpha)$,
according to \cite[Theorem 2]{1kipriyanov1960} acting as follows
\begin{equation}\label{2}
\mathfrak{D}^{\alpha}:\mathring{W}_p ^l (\Omega)\to L_{q}(\Omega),\quad
lp\leq n,\quad 0<\alpha<l- \frac{n}{p} +\frac{n}{q},\quad
 p\leq q<\frac{np}{n-lp}.
\end{equation}
If in the condition  \eqref{2}  we have  the strict inequality $q>p$, then
for sufficiently small $\delta>0$  the next inequality holds
\begin{equation}\label{3}
\|\mathfrak{D}^{\alpha}f\|_{L_{q}(\Omega)}\leq \frac{K}{\delta^{\nu}}\|f\|_{L_p(\Omega)}+\delta^{1-\nu}\|f\|_{L^{l}_p(\Omega)},
\end{equation}
where
\begin{equation}\label{4}
 \nu=\frac{n}{l}\Big(\frac{1}{p}-\frac{1}{q} \Big)+\frac{\alpha+\beta}{l}.
\end{equation}
The constant  $K$ is independent of $\delta,f$, and point
$P\in\partial\Omega \beta$ is an arbitrarily small fixed positive number.
Further we assume that $(0<\alpha<1)$.
Using the terminology of \cite{samko1987}, the left-side,
 right-side   classes of functions representable by the
fractional integral on the segment  we will denote respectively by
$I_{a+}^{\alpha}(L_p(a,b))$,
$I_{b-}^{\alpha}(L_p(a,b))$, $1\leq p\leq\infty$.
 Denote $\operatorname{diam}\Omega =\mathfrak{d}$;
$C,C_i$ are constants for $i\in \mathbb{N}_0$. We use for inner product
of points $P=(P_1,P_2,\dots ,P_n) $ and $Q=(Q_1,Q_2,\dots ,Q_n)$ which
belong to  $\mathbb{E}^{n}$ a contracted notations
$P\cdot Q=P^{i}Q_i=\sum^{n}_{i=1}P_iQ_i$. As usually    $D_iu$ denotes
 the generalized derivative of function $u$ with respect to coordinate
variable with index  $1\leq i\leq n$.
We will assume that all functions has a zero extension outside  of
$\bar{\Omega}$. Symbols: $\operatorname{D} (L), \mathrm{R} (L)$
denote respectively the domain of definition, and range of values of operator $L$.
Everywhere,  if not stated otherwise we will use the notations of
\cite{kipriyanov1960}, \cite{1kipriyanov1960},
\cite{samko1987}.
Let us define the operators by the following integral constructions
\begin{gather*}
 (\mathfrak{I}^{\alpha}_{0+}g)(Q  ):=\frac{1}{\Gamma(\alpha)}
 \int^{r}_0\frac{g (P+t\vec{e})}{( r-t)^{1-\alpha}}
\big(\frac{t}{r}\big)^{n-1}dt, (\mathfrak{I}^{\alpha}_{d-}g)(Q  )
:=\frac{1}{\Gamma(\alpha)} \int_{r}^{d }
\frac{g (P+t\vec{e})}{(t-r)^{1-\alpha}}dt, \\
\;g\in L_p(\Omega),\;1\leq p\leq\infty.
\end{gather*}
In this way was defined operators we will call respectively  the left-sided,
 right-sided operator of fractional integration in the direction.
We introduce  the classes of functions representable by the fractional
integral in the direction of $\vec{e}$
 \begin{gather}\label{5}
  \mathfrak{I}^{\alpha}_{0+}(L_p)
:=\{ u:\,u(Q)=(\mathfrak{I}^{\alpha}_{0+}g)(Q  ),\, g\in L_p(\Omega),
\,1\leq p\leq\infty \}, \\
\label{6}
  \mathfrak{I}  ^{\alpha}_{ d -} (L_p  )
 =\{ u:\,u(Q)=(\mathfrak{I}^{\alpha}_{d-}g)(Q  ),\;g\in L_p(\Omega),
\;1\leq p\leq\infty \}.
 \end{gather}
We define the families  of operators $\psi^{+}_{\varepsilon }$,
$\psi^{-}_{\varepsilon}$, $\varepsilon>0$ as follows:
$\operatorname{D}(\psi^{+}_{  \varepsilon }),\operatorname{D}(\psi^{-}_{  \varepsilon })
\subset L_p(\Omega)$. In the left-side case
 \begin{equation}\label{7}
(\psi^{+}_{ \varepsilon }f)(Q)
=\begin{cases}
 \int_0^{r-\varepsilon }\frac{ f (P+\vec{e}r)r^{n-1}
- f(P+\vec{e}t)t^{n-1}}{(  r-t)^{\alpha +1}r^{n-1}}  dt, &\varepsilon\leq r\leq d ,
\\[4pt]
 \frac{ f(Q)}{\alpha} \big(\frac{1}{\varepsilon^{\alpha}}-\frac{1}{ r ^{\alpha} }
 \big),&  0\leq r <\varepsilon .
\end{cases}
\end{equation}
In the right-side case
\begin{equation*}
 (\psi^{-}_{  \varepsilon }f)(Q)
=  \begin{cases}
 \int_{r+\varepsilon }^{d }\frac{ f (P+\vec{e}r)
- f(P+\vec{e}t)}{( t-r)^{\alpha +1}} dt, & 0\leq r\leq d -\varepsilon,\\[4pt]
   \frac{ f(Q)}{\alpha} \big(\frac{1}{\varepsilon^{\alpha}}
-\frac{1}{(d -r)^{\alpha} }    \big), & d -\varepsilon <r \leq d .
\end{cases}
\end{equation*}
Following \cite[p.181]{samko1987}  we define a truncated
fractional derivative similarly the derivative in the sense of Marchaud,
in the left-side case
 \begin{equation}\label{8}
 ( \mathfrak{D} ^{\alpha}_{0+,\varepsilon}f)(Q)
=\frac{1}{\Gamma(1-\alpha)}f(Q) r ^{-\alpha}+\frac{\alpha}{\Gamma(1-\alpha)}
(\psi^{+}_{  \varepsilon }f)(Q),
 \end{equation}
in the right-side case
 \begin{equation*}
 ( \mathfrak{D }^{\alpha}_{d-,\varepsilon}f)(Q)
=\frac{1}{\Gamma(1-\alpha)}f(Q)(d-r)^{-\alpha}
+\frac{\alpha}{\Gamma(1-\alpha)}(\psi^{-}_{  \varepsilon }f)(Q).
 \end{equation*}
Left-side  and  right-side fractional derivatives  accordingly  will be
understood  as a limits in the sense of norm
$L_p(\Omega)$, $1\leq p<\infty$ of truncated fractional derivatives
 \begin{equation*}
 \mathfrak{D }^{\alpha}_{0+}f
=\lim_{\stackrel{\varepsilon\to 0}{ (L_p) }} \mathfrak{D }^{\alpha}_{0+,\varepsilon} f  ,\; \mathfrak{D }^{\alpha}_{d-}f=\lim_{\stackrel{\varepsilon\to 0}{ (L_p) }} \mathfrak{D }^{\alpha}_{d-,\varepsilon} f .
\end{equation*}
We need several auxiliary propositions, which we will present in the next section.

\section{Results}

We have the following theorem on the boundedness of operators fractional integration
in a direction.

 \begin{theorem}\label{T1}
Operators of fractional integration in the direction are bounded in
$L_p(\Omega)$, $1\leq p<\infty$, the following estimates holds
  \begin{equation}\label{9}
 \| \mathfrak{I}^{\alpha}_{0 +}u\|_{L_p(\Omega)}\leq C\|u \|_{L_p(\Omega)},\;\|   \mathfrak{I} ^{\alpha}_{d -}u\|_{L_p(\Omega)}\leq C\|u \|_{L_p(\Omega)},\;C=  \mathfrak{d} ^{\alpha}/ \Gamma(\alpha+1)  .
 \end{equation}
 \end{theorem}

 \begin{proof}
Let us prove the first estimate of \eqref{9}, the proof of the second estimate
is analogous. Using the generalized Minkowski inequality, we have
\begin{align*}
\| \mathfrak{I}^{\alpha}_{0 +}u\|_{L_p(\Omega)}
&=\frac{1}{\Gamma(\alpha)}\Big( \int_{\Omega}\Big|
\int^{r}_0\frac{g (P+t\vec{e})}{( r-t)^{1-\alpha}}\big(\frac{t}{r}\big)^{n-1}dt
  \Big|^pdQ\Big)^{1/p} \\
&=\frac{1}{\Gamma(\alpha)}\Big( \int_{\Omega}\Big|
\int^{r}_0\frac{g (Q-\tau  \vec{e})}{\tau^{1-\alpha}}
 \big(\frac{r-\tau}{r}\big)^{n-1}d\tau  \Big|^pdQ\Big)^{1/p}\\
&\leq \frac{1}{\Gamma(\alpha)}\Big( \int_{\Omega}\Big(  \int^{\mathfrak{d}}_0
\frac{|g (Q-\tau  \vec{e})|}{\tau^{1-\alpha}}d\tau  \Big)^pdQ\Big)^{1/p} \\
&\leq \frac{1}{\Gamma(\alpha)} \int^{\mathfrak{d}}_0\tau^{\alpha-1} d\tau
\Big( \int_{\Omega}  |g (Q-\tau  \vec{e})|^p dQ  \Big)^{1/p} \\
&\leq  \frac{ \mathfrak{d} ^{\alpha}}{\Gamma(\alpha+1)}\,  \| u\|_{L_p(\Omega)}.
\end{align*}
\end{proof}

\begin{theorem}\label{T2}
Assume $f\in L_p(\Omega)$ and exists
$\lim_{\varepsilon\to  0}\psi^{+}_{  \varepsilon }f$ or
$\lim_{\varepsilon\to  0}\psi^{-}_{  \varepsilon }f$ in the sense of norm
$L_p(\Omega),\;1\leq p<\infty$. Then respectively
$f\in \mathfrak{I} ^{\alpha}_{0 +}(L_p) $ or
$f\in \mathfrak{I }^{\alpha}_{d -}(L_p)$.
\end{theorem}
\begin{proof}
 Let $f\in L_p(\Omega)$ and
$\lim_{\stackrel{\varepsilon\to 0}{ (L_p) }}\psi^{+}_{  \varepsilon }f=\psi$.
Consider the function
$$
(\varphi^{+}_{ \varepsilon}f)(Q)=\frac{1}{\Gamma(1-\alpha)}
\big\{\frac{f(Q)}{ r ^{\alpha}}+\alpha (\psi^{+}_{ \varepsilon }f)(Q) \big\}.
 $$
 Note    \eqref{7}, we can   easily see  that
$\varphi^{+}_{  \varepsilon }f\in L_p(\Omega)$. From fundamental property
family of functions  $\{\varphi^{+}_{  \varepsilon }f\}$ follows, that there is
exists a limit
 $\varphi^{+}_{  \varepsilon }f\to \varphi\in L_p(\Omega)$.
In consequence of proved in theorem \ref{T1} continuous property of operator
$\mathfrak{I}^{\alpha}_{0+}$ in the space $ L_p(\Omega)$, for completing
this theorem sufficient to show that
 $\lim_{\varepsilon\to 0,\, (L_p)} \mathfrak{I}^{\alpha}_{0+}
\varphi^{+}_{ \varepsilon }f=f$.
For $\varepsilon\leq r\leq d$, we have
\begin{align*}
&(\mathfrak{I}^{\alpha}_{0 +}\varphi^{+}_{ \varepsilon }f)(Q)
\frac{\pi r^{n-1}}{\sin\alpha\pi} \\
&=  \int_{\varepsilon}^{r}\frac{f (P+y\vec{e})y ^{n-1-\alpha}}{( r-y)^{1-\alpha} } dy\\
&\quad +\alpha\int_{\varepsilon}^{r}( r-y)^{\alpha-1}  dy
\int_0^{y-\varepsilon }\frac{ f (P+y\vec{e})y^{n-1}
 - f(P+t\vec{e})t^{n-1}}{( y-t )^{\alpha +1}} dt \\
&\quad  +\frac{1}{\varepsilon^{\alpha}}
 \int_0^{\varepsilon  }f (P+y\vec{e})( r-y)^{\alpha-1} y^{n-1}  dy =I.
\end{align*}
After  conversion in the second summand, we have
 \begin{equation}\label{10}
\begin{aligned}
I&=   \frac{1}{\varepsilon^{\alpha}}\int_0^{r  }f (P+y\vec{e})
 ( r-y)^{\alpha-1}y^{n-1}   dy \\
&\quad  - \alpha\int_{\varepsilon}^{r}( r-y)^{\alpha-1} dy
 \int_0^{y-\varepsilon }\frac{  f(P+t\vec{e})}{(  y-t)^{\alpha +1}}t^{n-1} dt .
\end{aligned}
\end{equation}
 Making a change  variable in the second integral, changing the order of
integration and going back to the previous variable, we obtain
\begin{equation}\label{11}
\begin{aligned}
&\alpha\int_{\varepsilon}^{r}( r-y)^{\alpha-1} dy
\int_0^{y-\varepsilon }\frac{  f(P+t\vec{e})}{(  y-t)^{\alpha +1}}t^{n-1} dt \\
& =\alpha\int_0^{r-\varepsilon}( r-y-\varepsilon)^{\alpha-1} dy
 \int_0^{y  }\frac{  f(P+t\vec{e})}{(  y+\varepsilon-t)^{\alpha +1}}t^{n-1} dt \\
&=\alpha\int_0^{r-\varepsilon}f(P+t\vec{e})t^{n-1} dt
 \int_{t }^{r-\varepsilon  }\frac{( r-y-\varepsilon)^{\alpha-1}   }
 {(  y+\varepsilon-t)^{\alpha +1}}dy \\
&=\alpha\int_0^{r-\varepsilon}f(P+t\vec{e})t^{n-1} dt
 \int_{t +\varepsilon}^{r  } ( r-y )^{\alpha-1}   (  y -t)^{-\alpha -1} dy .
\end{aligned}
\end{equation}
Applying   \cite[(13.18) p.184]{samko1987} we have
\begin{equation}\label{12}
 \int_{t +\varepsilon}^{r  } ( r-y )^{\alpha-1}   (  y -t)^{-\alpha -1} dy=
  \frac{1}{\alpha \varepsilon^{\alpha}}\frac{(r-t-\varepsilon)^{\alpha}}{ r-t }.
\end{equation}
Rewrite  \eqref{10} taking into account the relations \eqref{11}, \eqref{12},
after that  make change a variable  $t=r-\varepsilon\tau$, we obtain
\begin{equation}\label{13}
\begin{aligned}
&(\mathfrak{I}^{\alpha}_{0 +}\varphi^{+}_{ \varepsilon }f)(Q)
\frac{\pi r^{n-1}}{\sin\alpha\pi} \\
&=\frac{1}{\varepsilon^{\alpha}}
\Big\{ \int_0^{r  }f (P+y\vec{e})(r-y)^{\alpha-1}y^{n-1}   dy \\
&\quad - \int_0^{r-\varepsilon  } \frac{f(P+t\vec{e})(r-t-\varepsilon)^{\alpha}}{ r-t }
t^{n-1}dt \Big\} \\
&=\frac{1}{ \varepsilon^{\alpha}} \int_0^{r  }
\frac{f(P+t\vec{e})\left[(r -t)^{\alpha}-(r-t-\varepsilon)_{+}^{\alpha}\right]}
{ r-t }t^{n-1}dt \\
&=\int_0^{r/\varepsilon   }\frac{\tau^{\alpha}-(\tau-1)_{+}^{\alpha}}{ \tau }
f(P+[r-\varepsilon \tau ]\vec{e})(r-\varepsilon \tau)^{n-1} d\tau,
\end{aligned}
 \end{equation}
where $\tau_{+}=\begin{cases}
\tau, & \tau\geq0;\\
 0, &\tau<0.
\end{cases}$

Consider the auxiliary function $ \mathcal{K}$ defined in
 \cite[p.105]{samko1987} and having the next properties
 \begin{equation}\label{14}
 \mathcal{K}(t)= \frac{\sin\alpha\pi}{\pi }
\frac{ t_{+}^{\alpha}-(t-1)_{+}^{\alpha}}{ t }\in L_p(\mathbb{R}^{1}),\quad
\int_0^{\infty  }\mathcal{K}(t)dt=1,\quad \mathcal{K}(t)>0.
\end{equation}
From \eqref{13}, \eqref{14}, since  $f$ has a zero extension outside of
$\bar{\Omega}$, we have
 \begin{equation}\label{15}
 (\mathfrak{I}^{\alpha}_{0+}\varphi^{+}_{ \varepsilon }f)(Q)-f(Q)
=  \int_0^{\infty  }\mathcal{K}(t)
\{f(P+[r-\varepsilon t]\vec{e})(1-\varepsilon t/r)_{+}^{n-1}-f(P+ r \vec{e})\}dt.
\end{equation}
If $0\leq  r <\varepsilon$, then in accordance with \eqref{7}
 after the changing a variable, we obtain
 \begin{equation}\label{16}
\begin{aligned}
&(\mathfrak{I}^{\alpha}_{0+}\varphi^{+}_{ \varepsilon }f)(Q)-f (Q)\\
& = \frac{\sin\alpha\pi}{\pi\varepsilon^{\alpha}}
 \int_0^{r }\frac{f (P+t\vec{e})}{(r-t)^{1-\alpha} }
 \big(\frac{t}{r}big)^{n-1} dt-f (Q)\\
&=\frac{\sin\alpha\pi}{\pi\varepsilon^{\alpha}}
\int_0^{r }\frac{f (P+[r-t]\vec{e})}{t^{1-\alpha} }
\big(\frac{r-t}{r} \big)^{n-1} dt-f (Q).
\end{aligned}
\end{equation}
Consider the domains
 \begin{equation}\label{17}
 \Omega_{\varepsilon} :=\{Q\in\Omega:d(\vec{e})\geq\varepsilon \},\quad
\Omega_{-\varepsilon}=\Omega\setminus \Omega_{\varepsilon}.
 \end{equation}
Accordingly with this definition  we can  divide the surface $\omega$
into two  parts $ \omega'$ and  $\omega''$, where $\omega'$ is
the subset of  $\omega$ for which  $d(\vec{e})\geq\varepsilon$,
$\omega''$ is the subset of  $\omega$ for which $d(\vec{e}) <\varepsilon$.
Taking into account \eqref{15}, \eqref{16}, we obtain
\begin{equation}\label{18}
\begin{aligned}
&\|(\mathfrak{I}^{\alpha}_{0+}\varphi^{+}_{ \varepsilon }f) -f\|^p_{L_p(\Omega)} \\
&=   \int_{\omega'}d\chi\int_{\varepsilon}^{ d}
  \Big|\int_0^{\infty  }\mathcal{K}(t)[f(Q- \varepsilon t \vec{e})
(1-\varepsilon t/r)_{+}^{n-1}-f(Q)]dt\Big|^pr^{n-1}dr  \\
&\quad + \int_{\omega'}d\chi\int_0^{\varepsilon }
\Big|  \frac{\sin\alpha\pi}{\pi\varepsilon^{\alpha}}
 \int_0^{r }\frac{f (P+[r-t]\vec{e})}{t^{1-\alpha} }
 \big(\frac{r-t}{r}\big)^{n-1} dt-f (Q)\Big|^pr^{n-1}dr \\
&\quad +  \int_{\omega''}d\chi\int_0^{d }
\Big| \frac{\sin\alpha\pi}{\pi\varepsilon^{\alpha}}\int_0^{r }
\frac{f (P+[r-t]\vec{e})}{t^{1-\alpha} }\big(\frac{r-t}{r}\big)^{n-1} dt
-f (Q) \Big|^pr^{n-1}dr \\
&=I_1+I_2+I_3.
\end{aligned}
 \end{equation}
Consider $I_1$; using the generalized Minkovski's inequality  we obtain
$$
  I^{1/p}_1 \leq  \int_0^{\infty  }\mathcal{K}(t)
  \Big(\int_{\omega' }d\chi\int_{\varepsilon}^{ d}
  |f(Q- \varepsilon t \vec{e})(1-\varepsilon t/r)_{+}^{n-1}-f(Q)|^pr^{n-1}dr
\big)^{1/p} dt.
$$
Let us introduce the notation
$$
 \mathcal{K}(t)\Big(\int_{\omega' }d\chi\int_{\varepsilon}^{ d}
  |f(Q- \varepsilon t \vec{e})(1-\varepsilon t/r)_{+}^{n-1}-f(Q)|^pr^{n-1}dr \Big)^{1/p} dt=h(\varepsilon,t),
$$
we have the  inequality
\begin{equation}\label{19}
  |h(\varepsilon,t)|\leq 2\mathcal{K}(t) \| f\|_{L_p(\Omega)},\quad
\forall\varepsilon>0.
\end{equation}
Note that
\begin{align*}
|h(\varepsilon,t)|
&\leq  \Big(\int_{\omega' }d\chi\int_{\varepsilon}^{ d}
 \big|(1-\varepsilon t/r)_{+}^{n-1}
[f(Q- \varepsilon t \vec{e})-f(Q)]\big|^pr^{n-1}dr \Big)^{1/p} dt \\
&\quad +\Big(\int_{\omega' }d\chi\int_0^{ d}
  |f(Q)  [1-(1-\varepsilon t/r)_{+}^{n-1}]|^pr^{n-1}dr \Big)^{1/p} dt\\
&=I_{11}+I_{12}.
\end{align*}
In consequence of property continuity on average in space $L_p(\Omega)$,
for all fixed\\ $0<t<\infty$, we have $I_{11}\to 0$, $\varepsilon\to 0$.
Consider $I_{12}$, it is obvious that for all fixed  $0< t<\infty$,
almost everywhere in $\Omega$ the following relations holds
$$
h_1(\varepsilon,t,r)=\big|f(Q)  [1-(1-\varepsilon t/r)_{+}^{n-1}]\big|
\leq |f(Q)|,\;h_1(\varepsilon,t,r)\to 0,\quad \varepsilon \to 0.
$$
 Applying the  majorant  theorem of Lebesgue we obtain
 $I_{12}\to 0$, as $\varepsilon\to 0$.
It implies that for all fixed  $0< t<\infty$, we have
\begin{equation}\label{20}
\lim_{\varepsilon\to 0} h(\varepsilon,t)=0.
\end{equation}
Note \eqref{19}, \eqref{20}, again by applying   majorant  theorem of
 Lebesgue, we obtain
$$
 I_1\to 0,\quad\text{as }\varepsilon\to0  .
$$
Using Mincovski's inequality we can estimate   $I_2$,
\begin{align*}
I^{1/p}_2
&\leq \frac{\sin\alpha\pi}{\pi\varepsilon^{\alpha}}
\Big( \int_{\omega' }d\chi\int_0^{\varepsilon }
\Big| \int_0^{r }\frac{f (Q-t\vec{e})}{t^{1-\alpha} }
\big(\frac{r-t}{r}\big)^{n-1} dt\Big|^pr^{n-1}dr\Big)^{1/p}\\
&\quad +\Big(\int_{\omega' }d\chi\int_0^{\varepsilon}
| f (Q) |^pr^{n-1}dr \Big)^{1/p}=I_{21} +I_{22}.
\end{align*}
Applying the generalized  Mincovski's inequality, we obtain
\begin{align*}
&I_{21}\frac{\pi}{\sin\alpha\pi} \\
&=\frac{1}{\varepsilon^{\alpha}}\Big(\int_{\omega' }d\chi\int_0^{\varepsilon }
\Big|  \int_0^{r }\frac{f (Q-t\vec{e})}{t^{1-\alpha} }
\big(\frac{r-t}{r}\big)^{n-1} dt\Big|^pr^{n-1}dr \Big)^{1/p}\\
&\leq\frac{1}{\varepsilon^{\alpha}}
\Big\{\int_{\omega' }\Big[\int_0^{\varepsilon }t ^{\alpha-1 }
\Big( \int_{t}^{\varepsilon }|f (Q -t \vec{e})|^p
\big(\frac{r-t}{r}\big)^{(p-1)(n-1)}(r-t)^{n-1} dr  \Big)^{1/p}dt\Big]^pd\chi
 \Big\}^{1/p}\\
&\leq\frac{1}{\varepsilon^{\alpha}}
\Big\{\int_{\omega' }\Big[\int_0^{\varepsilon } t ^{\alpha-1 }
\Big( \int_{t}^{\varepsilon } |f (P+[r-t]\vec{e})|^p  (r-t)^{n-1}dr  \Big)^{1/p}dt
\Big]^pd\chi \Big\}^{1/p}\\
&\leq\frac{1}{\varepsilon^{\alpha}}\Big\{\int_{\omega' }
\Big[\int_0^{\varepsilon } t ^{\alpha-1 }
\Big( \int_0^{\varepsilon } |f (P +r \vec{e})|^p  r^{n-1}dr  \Big)^{1/p}dt\Big]^p
d\chi \Big\}^{1/p}\\
&=\frac{1}{\alpha}\| f\| _{L_p( \Delta_{\varepsilon})},
\end{align*}
where $\Delta_{\varepsilon} :=\{Q\in\Omega_{\varepsilon},\; r<\varepsilon \}$.
Note that  $\operatorname{meas} \Delta_{\varepsilon}\to 0$,
$\varepsilon\to 0$, hence  $I_{21},I_{22}\to 0 $. It implies that
$I_{2 }\to 0$. Applying analogous reasoning we can get  that
 $I_{3 }\to 0$, $\varepsilon\to 0$. According to the note  given above we
came to conclusion that
$ \mathfrak{I}^{\alpha}_{0+}\varphi^{+}_{ \varepsilon }f \to f$ in  $L_p$.
From the remark at the beginning of this proof, we  complete the
proof corresponding to the left-side case.
The proof for the right-side case is analogous. We have to show that
$\lim_{ \varepsilon\to 0} \mathfrak{I}^{\alpha}_{d-}\varphi^{-}_{ \varepsilon }f=f$
in the sense of $L_p$-norm, for this purpose  we must repeating the
previous reasoning with minor technical differences.
 \end{proof}

\begin{theorem}\label{T3}
Let $f=\mathfrak{I}^{\alpha}_{0+} \psi$ or
$f= \mathfrak{I} ^{\alpha}_{d -} \psi$,
$\psi\in L_p(\Omega)$, $1\leq p<\infty$. Then, respectively,
$ \mathfrak{D }^{\alpha}_{0+}f =\psi$
or $ \mathfrak{D} ^{\alpha}_{d-}f =\psi$,
in the sense of norm $L_p$.
\end{theorem}

\begin{proof}
Consider the difference
\begin{align*}
& r^{n-1}f(Q)-(r-\tau)^{n-1}f(Q-\tau\vec{e}) \\
&= \frac{1}{\Gamma(\alpha)} \int_0^{ r}
 \frac{\psi (Q-t\vec{e})}{t^{1-\alpha}}(r-t)^{n-1}dt
-\frac{1}{\Gamma(\alpha)} \int_{\tau}^{ r}
 \frac{\psi (Q-t\vec{e})}{(t-\tau)^{1-\alpha}}(r-t)^{n-1}dt\\
& =\tau^{\alpha-1}  \int_0^{ r} \psi (Q-t\vec{e})
 k(\frac{t}{\tau})(r-t)^{n-1}dt,\;k(t)\\
&= \frac{1}{\Gamma(\alpha)} \begin{cases}t^{\alpha-1}, &0<t<1;\\
t^{\alpha-1}-(t-1)^{\alpha-1}, & t>1.
\end{cases}
\end{align*}
Hence with the assumptions   $\varepsilon\leq r\leq d$, we have
\begin{align*}
(\psi^{+}_{  \varepsilon }f)(Q)
&=\int_{ \varepsilon }^{ r }\frac{r^{n-1} f(Q)
 -(r-\tau)^{n-1}f(Q-\tau\vec{e})}{ r^{n-1}\tau ^{\alpha +1}} d\tau \\
&=\int_{ \varepsilon }^{ r } \tau^{-2} d\tau\int_0^{ r}
 \psi (Q-t\vec{e}) k(\frac{t}{\tau})( 1-t/r )^{n-1}dt\\
&=\int_{ 0 }^{ r }\psi (Q-t\vec{e}) ( 1-t/r )^{n-1}  dt
 \int_{\varepsilon}^{ r}  k(\frac{t}{\tau}) \tau^{-2}d \tau \\
&=\int_{ 0 }^{ r }\psi (Q-t\vec{e}) ( 1-t/r )^{n-1}t^{-1} dt
 \int_{ t/ r  }^{t/\varepsilon}  k (s )ds.
\end{align*}
Applying  \cite[(6.12) p.106]{samko1987}, we obtain
$$
(\psi^{+}_{  \varepsilon }f)(Q)\cdot\frac{\alpha}{\Gamma(1-\alpha)}
=\int_{ 0 }^{ r }\psi (Q-t\vec{e})
 ( 1-t/r )^{n-1}\big[ \frac{1}{\varepsilon}\mathcal{K}(\frac{t}{\varepsilon})
- \frac{1}{ r}\mathcal{K}(\frac{t}{ r}) \big]dt,
$$
since in accordance with \eqref{14}  we have
 $$
\mathcal{K}(\frac{t}{ r})
=[\Gamma(1-\alpha)\Gamma(\alpha) ]^{-1} (\frac{t}{ r})^{\alpha-1},
$$
it follows that
$$
(\psi^{+}_{  \varepsilon }f)(Q)\cdot\frac{\alpha}{\Gamma(1-\alpha)}
=\int_{ 0 }^{ r /\varepsilon }\mathcal{K} ( t   ) \psi (Q-\varepsilon t\vec{e})
( 1-\varepsilon t/r )^{n-1}
  dt-\frac{f(Q)}{\Gamma(1-\alpha) r ^{ \alpha}}.
$$
Since  the function   $\psi(Q)$ extended by zero outside of  $\bar{\Omega}$,
 then if we note \eqref{8},\eqref{14}, we obtain
$$
 ( \mathfrak{D} ^{\alpha}_{0+,\varepsilon}f)(Q)-\psi(Q)
=\int_{ 0 }^{\infty} \mathcal{K} ( t   )
[\psi (Q - \varepsilon t \vec{e})(1-\varepsilon t/r)_{+}^{n-1}-\psi (Q) ]dt,
$$
for $\varepsilon\leq r\leq d $.
For values $r$ such that  $  0\leq r<\varepsilon$ by \eqref{7},  we have
$$
( \mathfrak{D} ^{\alpha}_{0+,\varepsilon}f)(Q)-\psi(Q)
=\frac{f(Q)}{\varepsilon^{\alpha}\Gamma(1-\alpha)}-\psi(Q).
$$
Using generalized Mincovski's inequality, we can get the estimate
\begin{align*}
\|( \mathfrak{D} ^{\alpha}_{0+,\varepsilon}f)(Q)-\psi(Q)\|_{L_p(\Omega)}
&\leq \int_{ 0 }^{\infty}  \mathcal{K}(t)   \| \psi(Q - \varepsilon t \vec{e})
(1-\varepsilon t/r)_{+}^{n-1}-\psi (Q)\|_{L_p(\Omega)}dt \\
&\quad +\frac{1}{ \Gamma(1-\alpha)\varepsilon^{\alpha}}
 \|f\|_{L_p(\Delta'_{\varepsilon})}+\|\psi\|_{L_p(\Delta'_{\varepsilon})},
\end{align*}
where $\Delta'_{\varepsilon}=\Delta_{\varepsilon}\cup \Omega_{-\varepsilon}$.
Note that as it shown for the right side of   \eqref{18}, we can conclude
 that  all tree summands of the right side  of last inequality tends to
zero as $\varepsilon\to 0$.
\end{proof}

\begin{theorem} \label{T4}
Let $\rho\in \operatorname{Lip}\lambda$, $\alpha<\lambda\leq 1$,
$f\in H_0^{1}(\Omega) $, then $\rho f\in  \mathfrak{I} ^{\alpha}_{0 +}(L_2 )
 \cap \mathfrak{I} ^{\alpha}_{d -}(L_2 )$.
\end{theorem}

\begin{proof}
We  provide a proof only for the left-side  case, because  the proof
corresponding to the right-side case is analogous.
Suppose that all functions have  a zero extension outside of $\bar{\Omega}$.
 At first assume $f\in C_0^{\infty}(\Omega) $
 and in terms of notations \eqref{17} let us consider domains
$\Omega' = \Omega_{\varepsilon_1} $, $\Omega''=\Omega_{-\varepsilon_1} $,
we have
\begin{equation}\label{21}
\|\psi^{+}_{ \varepsilon_1}f  - \psi^{+}_{ \varepsilon_2}f \|_{L_2(\Omega )}
\leq\|\psi^{+}_{ \varepsilon_1}f  - \psi^{+}_{ \varepsilon_2}f \|_{L_2(\Omega')}
+\|\psi^{+}_{ \varepsilon_1}f  - \psi^{+}_{ \varepsilon_2}f \|_{L_2(\Omega'')}.
\end{equation}
Denote:   $\rho(P+\vec{e}t)t^{n-1}=\sigma (P+\vec{e}t) $ and consider
\begin{align*}
&\|\psi^{+}_{ \varepsilon_1}f  - \psi^{+}_{ \varepsilon_2}f \|_{L_2(\Omega')}\\
&\leq  \Big(\int_{\omega'}d\chi\int_{\varepsilon_1}^{d }
\Big|\int_{r-\varepsilon_1}^{r-\varepsilon_2}\frac{(\sigma  f )(Q)
 -(\sigma  f )(P+\vec{e}t)}{r^{n-1}(r-t)^{\alpha +1}} dt \Big|^2r^{n-1} dr
\Big)^{1/2} \\
&\quad +\Big(\int_{\omega'}d\chi\int_{ \varepsilon_2}^{ \varepsilon_1}
\Big|\int_0^{r- \varepsilon_1}\frac{(\sigma  f )(Q) }{r^{n-1}(r-t)^{\alpha +1}} dt \\
&\quad - \int_0^{r-\varepsilon_2}\frac{(\sigma  f )(Q)-(\sigma  f )
(P+\vec{e}t)}{r^{n-1}(r-t)^{\alpha +1}} dt\Big|^2r^{n-1} dr  \Big)^{1/2} \\
&\quad  +\Big(\int_{\omega'}d\chi\int_0^{\varepsilon_2 }
\Big|\int_0^{r-\varepsilon_1}
 \frac{(\sigma  f )(Q) }{r^{n-1}(r-t)^{\alpha +1}} dt \\
&\quad - \int_0^{r- \varepsilon_2}\frac{(\sigma  f )(Q) }{r^{n-1}(r-t)^{\alpha +1}} dt
\Big|^2r^{n-1} dr  Big)^{1/2} \\
& =I_1+I_2+I_3.
\end{align*}
Note since  $f\in C_0^{\infty}(\Omega)$, then for sufficient small
 $\varepsilon_1: f(Q)=0$, $r<\varepsilon_1$. This implies that
$I_2+I_3=0$ and also implies  that second summand of  \eqref{21} equals zero.
Making the change of variable in $I_1$, we have
$$
 I_1  =\Big(\int_{\omega'}d\chi\int_{\varepsilon_1}^{d}
\Big|\int_{ \varepsilon_1}^{ \varepsilon_2}\frac{(\sigma f )(Q)
 -(\sigma  f )(Q-\vec{e} t )}{ r^{n-1}t ^{\alpha +1}} dt\Big|^2r^{n-1} dr \Big)^{1/2}.
$$
Using the generalized Minkowski's inequality we obtain
\begin{align*}
I_1
&\leq\int_{ \varepsilon_2}^{ \varepsilon_1}t^{-\alpha -1}
\Big(\int_{\omega'}d \chi
\int_{\varepsilon_1}^{d}\Big| (\rho  f )(Q)-(1-t/r)^{n-1}(\rho  f )(Q-\vec{e}t)
\Big|^2  r^{n-1} dr  \Big)^{1/2} dt\\
& \leq\int_{ \varepsilon_2}^{ \varepsilon_1}t^{-\alpha -1}\Big(\int_{\omega'}
d \chi\int_{\varepsilon_1}^{d}\Big| (\rho f )(Q)- (\rho f )(Q-\vec{e}t) \Big|^2
 r^{n-1} dr  \Big)^{1/2} dt \\
&\quad +\int_{ \varepsilon_2}^{ \varepsilon_1}t^{-\alpha -1}
\Big(\int_{\omega'}d \chi\int_{\varepsilon_1}^{d}\left[ 1-( 1-t/r )^{n-1}\right]
| (\rho f )(Q-\vec{e}t) |^2  r^{n-1} dr \Big)^{1/2}dt \\
&\leq C_1\int_{ \varepsilon_2}^{ \varepsilon_1}t^{\lambda-\alpha-1 }d t \\
&\quad +\int_{ \varepsilon_2}^{ \varepsilon_1}t^{-\alpha }
\Big(\int_{\omega'}d \chi\int_{\varepsilon_1}^{d}
\Big|\frac{1}{r} \sum_{i=0}^{n-2}\big( \frac{t}{r}\big)^{i}(\rho f )(Q-\vec{e}t)
\Big|^2  r^{n-1} dr  \Big)^{1/2}dt.
\end{align*}
In consequence of the boundendedness property of function $f$, exists constant
 $\delta$ such that $  f  (Q-\vec{e}t )=0,\;r<\delta$. Finally   we have
the estimate
 $$
 I_1 \leq  \frac{C_1   }{\lambda-\alpha  }
(\varepsilon^{  \lambda-\alpha}_1-\varepsilon^{\lambda-\alpha}_2)
+ \|f\|_{L_2(\Omega)}\frac{ (n-1)  }{\delta(1-\alpha)  }
(\varepsilon^{  1- \alpha}_1-\varepsilon^{1- \alpha}_2)  .
$$
By theorem \ref{T1} we have the inclusion
$ \rho f \in  \mathfrak{I} ^{\alpha}_{0+}(L_2 ),\;f\in C_0^{\infty}(\Omega)$.

Let $f\in H^1_0(\Omega)$, then there exists a sequence
$\{f_n\} \subset C_0^{\infty}(\Omega)$,
$\rho f_n\stackrel{ L_2}{\longrightarrow}\rho f$.
According to the part proved above, we have
$\rho f_n= \mathfrak{I} ^{\alpha}_{0+}\varphi_n$,
$\{\varphi_n\}\in L_2(\Omega)$, therefore,
\begin{equation}\label{22}
 \mathfrak{I} ^{\alpha}_{0+}\varphi_n\stackrel{L_2 }{\longrightarrow} \rho f.
\end{equation}
We will show that exists $\varphi\in L_2(\Omega)$ such that
$\varphi_n\stackrel{L_2 }{\longrightarrow}\varphi$.
Note, that  by theorem  \ref{T2} we have
$\mathfrak{ D} ^{\alpha}_{0+}\rho f_n=\varphi_n$.
Thus introducing the notation $f_{n+m}-f_n=c_{n,m}$, we obtain
\begin{align*}
\|\varphi_{n+m}-\varphi_n\|_{L_2(\Omega)}
&\leq\frac{\alpha}{\Gamma(1-\alpha)}\Big(\int_{\Omega}
\Big|\int_0^{r}\frac{(\sigma c_{n,m})(Q)-(\sigma c_{n,m})(P+\vec{e}t)}{r^{n-1}
( t-r)^{\alpha +1}} dt\Big|^2 dQ  \Big)^{1/2} \\
&\quad +\frac{1}{\Gamma(1-\alpha)}\Big(\int_{\Omega}
\Big| \frac{(\rho c_{n,m})(Q) }{ r ^{\alpha }} \Big|^2dQ \Big)^{1/2}
=I_1 +I_2.
\end{align*}
Let us estimate $I_1 $,
\begin{align*}
\frac{\Gamma(1-\alpha)}{\alpha} I_1
&\leq  \Big\{\int_{\Omega}\Big|\int_0^{ r}\frac{  (\rho c_{n,m}) (Q)
-  (\rho c_{n,m}) (Q-\vec{e}t) }{ t^{\alpha +1}} dt\Big|^2 dQ \Big\}^{1/2} \\
&\quad  + \Big\{\int_{\Omega} \Big|\int_0^{ r}\frac{(\rho c_{n,m})
 (Q-\vec{e}t)[1-(1- t/r)^{n-1} ]  }{ t^{1+\alpha   }} dt\Big|^2 dQ \Big\}^{1/2}\\
&=  I_{01}+I_{02}.
\end{align*}
Now we consider $I_{01}$. It is obvious that
\begin{align*}
I_{01}
&\leq \sup_{Q\in \Omega}|\rho(Q)|\Big\{\int_{\Omega}
\Big(\int_0^{ r}\frac{  |  c_{n,m} (Q)-    c_{n,m} (Q-\vec{e}t) |}{ t^{\alpha +1}} dt
\Big)^2 dQ    \Big\}^{1/2} \\
&\quad +\Big\{\int_{\Omega}\Big|\int_0^{ r}\frac{ c_{n,m} (Q-\vec{e}t)
 [  \rho (Q)-    \rho (Q-\vec{e}t)] }{ t^{\alpha +1}} dt\Big|^2 dQ
 \Big\}^{1/2} \\
&=  I_{11}+I_{21} .
\end{align*}
Applying  the generalized Minkowski's inequality, and representing the function under
the integral by the derivative in the direction of $\vec{e}$, we obtain
\begin{align*}
I_{11}
&\leq C_1\int_0^{ \mathfrak{d}}t^{-\alpha -1}
\Big(\int_{\Omega} \big| c_{n,m}(Q)-c_{n,m}(Q-\vec{e}t) \big|^2  dQ \Big)^{1/2} dt\\
&=C_1\int_0^{ \mathfrak{d}}t^{-\alpha -1}
\Big(\int_{\Omega} \Big| \int_0^{t} c'_{n,m} (Q-\vec{e}\tau) d\tau\Big|^2
dQ\Big)^{1/2} dt.
\end{align*}
Using the Cauchy-Schwarz inequality, and  Fubini's  theorem, we have
\begin{align*}
I_{11} &\leq C_1\int_0^{ \mathfrak{d}}t^{-\alpha -1}
\Big(\int_{\Omega}dQ \int_0^{t} \left|c'_{n,m} (Q-\vec{e}\tau)\right|^2d\tau
\int_0^{t}d\tau    \Big)^{1/2} dt\\
& =C_1\int_0^{ \mathfrak{d}}t^{-\alpha -1/2 }
\Big(\int_0^{t}d\tau \int_{\Omega}\left|c'_{n,m} (Q-\vec{e}\tau)\right|^2 dQ
 \Big)^{1/2} dt\\
&\leq C_1  \frac{ \mathfrak{d}^{1-\alpha}  }{1-\alpha  }
 \|c'_{n,m} \|_{L_2(\Omega)}.
\end{align*}
  Using Holder's property of function $\rho$
 analogously  to the previous reasoning, we have the following estimate
$$
I_{21}
\leq M \int_0^{ \mathfrak{d}}t^{\lambda-\alpha -1}
\Big(\int_{\Omega} \left|  c_{n,m}(Q-\vec{e}t) \right|^2  dQ Big)^{1/2} dt \\
\leq M\frac{ \mathfrak{d}^{\lambda-\alpha}  }{\lambda-\alpha }
  \|c _{n,m} \|_{L_2(\Omega)}.
$$
Applying trivial estimates, we have
\begin{align*}
I_{02}
&\leq C_1\Big\{\int_{\Omega} \Big|\int_0^{ r}  | c_{n,m}(Q-\vec{e}t)
| \sum_{i=0}^{n-2}\big(\frac{t}{r}\big)^{i}  r^{-1} t^{ -\alpha   }
 dt\Big|^2 dQ    \Big\}^{1/2}\\
&\leq C_2\Big\{\int_{\Omega} \Big|\int_0^{ r}  | c_{n,m}(Q-\vec{e}t)|
  r^{-1} t^{ -\alpha   }  dt\Big|^2 dQ    \Big\}^{1/2}\\
&\leq C_2\Big\{\int_{\Omega} \Big(\int_0^{ r} t^{ -\alpha   }dt \int_{t}^{r}
\left| c'_{n,m}(Q-\vec{e}\tau)\right| d\tau  \Big)^2r^{-2}  dQ    \Big\}^{1/2} \\
&= C_2\Big\{\int_{\Omega} (\int_0^{ r}\left| c'_{n,m}(Q-\vec{e}\tau)\right|
 d\tau  \int_0^{\tau}t^{ -\alpha   }dt      )^2r^{-2}  dQ    \Big\}^{1/2} \\
& =\frac{C_2}{1-\alpha}\Big\{\int_{\Omega} \Big(\int_0^{ r}
\left| c'_{n,m}(Q-\vec{e}\tau)\right|(\frac{\tau}{r}) \tau^{-\alpha} d\tau Big)^2
dQ    \Big\}^{1/2}\\
& \leq\frac{C_2}{1-\alpha}\Big\{\int_{\Omega} \Big(\int_0^{ r}
\left| c'_{n,m}(Q-\vec{e}\tau)\right| \tau^{ -\alpha} d\tau \Big)^2  dQ \Big\}^{1/2}.
\end{align*}
Using the generalized Minkowski's inequality, we have
$$
I_{02}\leq C_3\int_0^{ \mathfrak{d}} \tau^{ -\alpha} d\tau
\Big(\int_{\Omega} \left| c'_{n,m}(Q-\vec{e}\tau)\right|^2 dQ  \Big)^{1/2}
\leq C_3\frac{\mathfrak{d}^{1-\alpha}}{1-\alpha} \|c'_{n,m}\|_{L_2(\Omega)}.
$$
Consider   $I_2$, representing the function under the integral by the
 derivative in the direction of $\vec{e}$,
\begin{align*}
I_2
& \leq\frac{C_{ 1}}{\Gamma(1-\alpha)}
\Big(\int_{\Omega} \left| c_{n,m}(Q)\right|^2   r ^{-2\alpha  } dQ
\Big)^{ \frac{1}{2}} \\
&=\frac{C_{ 1}}{\Gamma(1-\alpha)} (\int_{\Omega}  r ^{-2\alpha  }
\Big|\int_0^{r} c'_{n,m}(Q-\vec{e}t)dt\Big|^2dQ \Big)^{1/2}\\
& \leq\frac{C_{ 1}}{\Gamma(1-\alpha)}
\Big(\int_{\Omega}    \Big|\int_0^{r} c'_{n,m}(Q-\vec{e}t)t^{-\alpha}dt\Big|^2dQ
\Big)^{1/2}.
\end{align*}
Using the generalized Minkowski's inequality, then applying  obvious estimates
 we have
\begin{align*}
 I_2
&\leq  C_4 \Big\{\int_{\omega}\Big[\int_0^{d} t^{-\alpha} dt
\Big(\int_{t}^{d}|c'_{n,m}(Q-\vec{e}t)|^2   r^{ n-1 }dr\Big)^{1/2} \Big]^2d\chi
\Big\}^{1/2} \\
& \leq  C_4 \Big\{\int_{\omega}\Big[\int_0^{\mathfrak{d}} t^{-\alpha} dt
 \Big(\int_0^{d}|c'_{n,m}(Q-\vec{e}t)|^2   r^{ n-1 }dr\Big)^{1/2} \Big]^2d\chi
\Big\}^{1/2} \\
& =  C_4\int_0^{\mathfrak{d}} t^{-\alpha}   dt
\Big(\int_{\omega}d\chi\int_0^{d}|c'_{n,m}(Q-\vec{e}t)|^2   r^{ n-1 }dr\Big)^{1/2}  \\
&\leq  C_4\frac{\mathfrak{d}^{1-\alpha}}{1-\alpha}\|c'_{n,m}\|_{L_2(\Omega)}.
\end{align*}
 From the  fundamental  property of the sequences
 $\{c_{n,m}\}$, $\{c'_{n,m}\}$, it follows that  $I_1,I_2\to 0$. Hence the sequence
 $\{\varphi_n\} $ is fundamental and in consequence of  completeness property
 of the space $L_2(\Omega)$ there is  a limit of sequence  $\{\varphi_n\}$,
a  some function $\varphi \in L_2(\Omega)$.
In consequence of  theorem \ref{T1} the   operator of fractional integration
in the direction  is boundary acting in the space $L_2(\Omega)$, hence
$$
 \mathfrak{I }^{\alpha}_{0 +}\varphi_n\stackrel{L_2 }{\longrightarrow}
\mathfrak{I} ^{\alpha}_{0+}\varphi.
$$
 Since  \eqref{22} holds, we have   $\rho f=  \mathfrak{I} ^{\alpha}_{0+}\varphi$.
\end{proof}

\begin{lemma}\label{L1}
The operator $\mathfrak{D}^{ \alpha }$ is a contraction of operator
$\mathfrak{D}^{ \alpha }_{0+}$, exactly
$\mathfrak{D}^{ \alpha }\subset \mathfrak{D}^{ \alpha }_{0+}$.
\end{lemma}

\begin{proof}
We  show that the  equality holds
\begin{equation}\label{23}
 (\mathfrak{D}^{ \alpha }f)(Q)= (\mathfrak{D}^{ \alpha }_{0+} f)(Q), \quad
f\in \stackrel{0}{W_p ^l} (\Omega).
  \end{equation}
It follows  from the next obvious conversions:
\begin{equation}
\begin{aligned}
&r^{n-1}\mathfrak{D}^{ \alpha }v \\
&=\frac{\alpha}{\Gamma(1-\alpha)}\int_0^{r}
 \frac{v(Q)-v(P+\vec{e}t)}{(r- t)^{\alpha+1}}  t  ^{n-1} dt
+ \frac{C_n^{(\alpha)}}{\Gamma(1-\alpha)}   v(Q)   r ^{ n-1-\alpha} \\
& = \frac{\alpha}{\Gamma(1-\alpha)}\int_0^{r} \frac{r^{n-1}v(Q)
-t^{n-1}v(P+\vec{e}t)}{(r- t)^{\alpha+1}}dt\\
&\quad -v(Q) \frac{\alpha}{\Gamma(1-\alpha)}\int_0^{r} \frac{ r^{n-1} -t^{n-1} }
{(r- t)^{\alpha+1}}   dt
 +\frac{(n-1)!}{\Gamma(n-\alpha)}   v(Q)   r ^{ n-1-\alpha} \\
&= ( \mathfrak{D} ^{\alpha}_{0+}t^{n-1}v)(Q)
 -\frac{\alpha v(Q) }{\Gamma(1-\alpha)} \sum_{i=0}^{n-2}  r ^{n-2-i}
\int_0^{r} \frac{t^{i}}{(r- t)^{\alpha }}   dt \\
&\quad +\frac{(n-1)! } {\Gamma(n-\alpha)}  v(Q)   r ^{ n-1-\alpha}
 -\frac{1}{\Gamma(1-\alpha)}   v(Q)   r ^{ n-1-\alpha} \\
&= ( \mathfrak{D} ^{\alpha}_{0+}t^{n-1}v)(Q)- I_1 +I_2 -I_3.
\end{aligned} \label{24}
\end{equation}
Let us conduct  the following conversions  using the formula of fractional
integration of the exponential function  \cite[(2.44) p.47]{samko1987}
\begin{align*}
I_1
&= \frac{\alpha v(Q) }{\Gamma(1-\alpha)}r ^{n-2 }\int_0^{r}
 \frac{1}{(r- t)^{\alpha }}   dt
 +\frac{\alpha v(Q) }{\Gamma(1-\alpha)} \sum_{i=1}^{n-2}
r ^{n-2-i}\int_0^{r} \frac{t^{i}}{(r- t)^{\alpha }}   dt \\
& =  v(Q)\frac{\alpha }{ \Gamma(2-\alpha)}r ^{n-1-\alpha }
 +  v(Q)  \alpha \sum_{i=1}^{n-2}
r ^{n-2-i} (I^{1-\alpha}_{0+}t^{i})(r)   \\
&=v(Q)\frac{\alpha }{ \Gamma(2-\alpha)}r ^{n-1-\alpha }
+  v(Q)  \alpha \sum_{i=1}^{n-2}r ^{n-1-\alpha}\frac{i!}{\Gamma(2-\alpha+i)}.
\end{align*}
Consequently
\begin{equation}\label{25}
\begin{aligned}
r ^{-n+1+\alpha }(I_1 +I_3)/v(Q)
&=\frac{1}{ \Gamma(2-\alpha)}
 + \alpha \sum_{i=1}^{n-2} \frac{i!}{\Gamma(2-\alpha+i)}\\
&= \frac{2}{ \Gamma(3-\alpha)}
 +   \alpha \sum_{i=2}^{n-2} \frac{i!}{\Gamma(2-\alpha+i)} \\
&= \frac{3!}{ \Gamma(4-\alpha)}
 +  \alpha \sum_{i=3}^{n-2} \frac{i!}{\Gamma(2-\alpha+i)} \\
&=\frac{(n-2)!}{ \Gamma(n-1-\alpha)}
 +  \alpha\frac{(n-2)!}{\Gamma(n-\alpha )}=\frac{(n-1)!}{ \Gamma(n-\alpha)}.
\end{aligned}
\end{equation}
Hence $I_2-I_1-I_3=0$,
and equality \eqref{23}  follows from  \eqref{24},\eqref{25}.
The proof of the fact  of difference   the operators $\mathfrak{D}^{ \alpha }$
 and $\mathfrak{D}^{ \alpha }_{0+}$ implies from the following reasoning.
Let $f\in \mathfrak{I}^{\alpha}_{0+}  \varphi,\;\varphi \in L_p(\Omega) $,
then in consequence of theorem \ref{T2} we have
$ \mathfrak{D} _{0+} ^{ \alpha }  \mathfrak{I}^{\alpha}_{0+}   \varphi  =\varphi$.
Hence  $\mathfrak{I}^{\alpha}_{0+} ( L_p)\subset\operatorname{D}
(\mathfrak{D}^{ \alpha }_{0+})$. Given the above remains to note that
there exists $f\in \mathfrak{I}^{\alpha}_{0+} ( L_p)$, such that
$$
f(\Lambda)\neq 0,\;\Lambda \subset\partial\Omega,\quad
\operatorname{meas} \Lambda\neq 0,
$$
at the same time
$$
f(\partial \Omega)= 0,\quad \forall f\in \operatorname{D}(\mathfrak{D}^{ \alpha }  ) .
$$
\end{proof}

\begin{lemma}\label{L2}
The following equality holds
\begin{equation*}
  \mathfrak{D} _{0+} ^{ \alpha ^{*}}  = \mathfrak{D }^{\alpha}_{d-},
\end{equation*}
where
$$
\operatorname{D}( \mathfrak{D }_{0+} ^{ \alpha } )=\mathfrak{I}^{\alpha}_{0+} ( L_2),\quad
\operatorname{D}(\mathfrak{D }^{\alpha}_{d-})=  \mathfrak{I } ^{\alpha}_{d-}(  L_2).
 $$
\end{lemma}

\begin{proof}
We  show  the next relation
\begin{gather}\label{26}
 (\mathfrak{D}^{ \alpha }_{0+} f ,  g  )_{L_2(\Omega)}
=(f , \mathfrak{D }^{\alpha}_{d-} g )_{L_2(\Omega)}, \\
 f\in \mathfrak{I}^{\alpha}_{0+} ( L_2),\,g\in  \mathfrak{I } ^{\alpha}_{d-}(  L_2).
\nonumber
\end{gather}
Note that as a consequence of theorem  \ref{T3} the next equalities  holds:
 $\mathfrak{D}^{ \alpha }_{0+} \mathfrak{I}^{\alpha}_{0+} \varphi
 =\varphi\in L_2(\Omega)$,
$\mathfrak{D }^{\alpha}_{d-}  \mathfrak{I } ^{\alpha}_{d-}   \psi
=\psi\in L_2(\Omega)$.
Given the above and in consequence of theorem \ref{T1}, we obtain that the
left and right side of \eqref{26} exist and are finite.
Using the Fubini's theorem we can  perform the following conversions
\begin{equation}\label{27}
\begin{aligned}
(\mathfrak{D}^{ \alpha }_{0+} f ,  g  )_{L_2(\Omega)}
&=\int_{\omega}d\chi \int_0^{d}\varphi(P+\vec{e}r)
\overline{(\mathfrak{I } ^{\alpha}_{d-}\psi)(Q)}r^{n-1}dr \\
&=\frac{1}{\Gamma(\alpha)}\int_{\omega}d\chi\int_0^{d}\varphi(P+\vec{e}r)r^{n-1}dr
 \int_{r}^{d}\frac{\overline{\psi(P+\vec{e}t)}}{(t-r)^{1-\alpha}}\,dt \\
&=\frac{1}{\Gamma(\alpha)}\int_{\omega}d\chi\int_0^{d}
 \overline{\psi(P+\vec{e}t)}t^{n-1}dt
 \int_0^{t}\frac{\varphi(P+\vec{e}r)}{(t-r)^{1-\alpha}}( \frac{r}{t})^{n-1}dr \\
&=\int_{\Omega} (\mathfrak{I} ^{\alpha}_{0+} \varphi)(Q)\,\overline{\psi(Q)} \, dQ
 =(f , \mathfrak{D }^{\alpha}_{d-} g )_{L_2(\Omega)}.
\end{aligned}
\end{equation}
Inequality \eqref{26} is proved. From  equality  \eqref{26} follows that
  $\operatorname{D}(\mathfrak{D }^{\alpha}_{d-})\subset\operatorname{D}
( \mathfrak{D}_{0+}^{ \alpha^{*}})$.
Since  $\mathrm{R}(\mathfrak{D }^{\alpha}_{d-})=L_2$, then
$\mathrm{R}(\mathfrak{D}_{0+}^{ \alpha^{*}})=L_2$.
We will show that
$\operatorname{D}(\mathfrak{D}_{0+}^{ \alpha^{*}})\subset\operatorname{D}
(\mathfrak{D }^{\alpha}_{d-})$ thus completing the proof.
In accordance with the definition of conjugate operator,
for all element   $f\in \operatorname{D}(\mathfrak{D}^{ \alpha }_{0+}) $
 and   pars  of elements
$g\in \operatorname{D} (\mathfrak{D}_{0+}^{ \alpha^{*}}),\,g^{*}\in
  \mathrm{R}(\mathfrak{D}_{0+}^{ \alpha^{*}})$, the integral equality holds
\[
\langle \mathfrak{D}^{ \alpha }_{0+} f,g \rangle_{L_2(\Omega)}
=\langle   f,g^{*} \rangle_{L_2(\Omega)}.
\]
Suppose $f=\mathfrak{I}^{\alpha}_{0+} \varphi,\;\varphi  \in L_2(\Omega)$.
 Using the Fubini's theorem  and  performing  the conversion similar to \eqref{27},
we have
$$
\langle \mathfrak{D}^{ \alpha }_{0+} f,g- \mathfrak{I } ^{\alpha}_{d-}g^{*}
\rangle_{L_2(\Omega)}=0.
$$
By theorem \ref{T3} the image of operator $\mathfrak{D}^{ \alpha }_{0+}$
coincides with the space $L_2(\Omega)$. Hence   the element
 $(g- \mathfrak{I } ^{\alpha}_{d-}g^{*})\in L_2$ equals zero.
It implies that $\operatorname{D}(\mathfrak{D}_{0+}^{ \alpha^{*}})
\subset\operatorname{D}(\mathfrak{D }^{\alpha}_{d-})$.
\end{proof}

\section{Strong accretiveness property}

 The  following theorem establishes the strong accretive  property
 (see \cite[p. 352]{kato1966}) for the operator of
 fractional differentiation in the sense of Kipriyanov
acting in the complex weight space of Lebesgue summable
with squared functions.

\begin{theorem}\label{T5}
  Let $n\geq2,\;\rho(Q)$ is non-negative real function in  class
$\operatorname{Lip} \mu$, $\mu>\alpha$.
Then for  the operator of fractional differentiation in the sense
of Kipriyanov the inequality of a strong accretiveness    holds
\begin{equation}\label{28}
 \operatorname{Re} \langle f,\mathfrak{D}^{\alpha}f\rangle_{L_2(\Omega,\rho)}
\geq\frac{1}{\lambda^2}\|f\|^2_{L_2(\Omega,\rho)},\quad
f\in H^{1}_0 (\Omega).
\end{equation}
\end{theorem}

\begin{proof}
 First we assume that $f$ is real.
For $f\in C_0^{\infty}(\Omega) $ consider the following difference
 in which the second summand exists due to theorem \ref{T4},
\begin{align*}
&\rho(Q)f(Q)( \mathfrak{D}^{\alpha}f )(Q)-\frac{1}{2}
 (\mathfrak{D}^{\alpha} \rho f^2 ) (Q) \\
&=\frac{\alpha}{2\Gamma(1-\alpha)} \int_0^{r}
 \frac{\rho(Q)[f (P+\vec{e} r)- f(P+\vec{e}t)]^2}{(r - t)^{\alpha+1}}
\big(\frac{t}{r}\big)^{n-1}dt \\
&\quad +\frac{C^{(\alpha)}_n}{2}  \rho(Q)|f(Q)|^2r^{ -\alpha}dr\geq0.
\end{align*}
Therefore,
\begin{equation}\label{29}
\rho(Q)f(Q)( \mathfrak{D}^{\alpha}f )(Q)
\geq\frac{1}{2}(\mathfrak{D}^{\alpha} \rho f^2  )(Q).
\end{equation}
Integrating the left and right sides of inequality \eqref{29},
 then using a Fubini's theorem we obtain
\begin{align*}
&\int_0^{d   }
 f(Q)( \mathfrak{D}^{\alpha}f )(Q)\rho(Q) r^{n-1}dr \\
&\geq \frac{1}{2}\int_0^{d  } ( \mathfrak{D}^{\alpha} \rho f^2 )(Q)r^{n-1}dr \\
&= \frac{\alpha}{2\Gamma(1-\alpha)}\int_0^{d}t^{n-1}dt\int_{t}^{d}
\frac{(\rho f^2)(Q) -(\rho f^2)(P+\vec{e}t) }{(r - t)^{\alpha+1}} dr \\
&\quad + \frac{C^{(\alpha)}_n}{2} \int_0^{d}(\rho f ^2)(Q)r^{n-1-\alpha} dr \\
&=-\frac{1}{2 }\int_0^{d  }(\mathfrak{D }^{\alpha}_{d-}\rho f^2 )(Q)r^{n-1}dr
+\frac{C^{(\alpha)}_n}{2}\int_0^{d   }(\rho f ^2)(Q)r^{n-1-\alpha}dr \\
&\quad +\frac{1}{2\Gamma(1-\alpha)}\int_0^{d  } (\rho f ^2)(Q)r^{n-1}
 (d  -r)^{-\alpha} dr=I  .
\end{align*}
Let us rewrite the first summand of the last sum using the formula of
fractional integration of   exponential function
 \cite[(2.44) p.47]{samko1987}, we obtain
$$
 \int_0^{d  }(\mathfrak{D }^{\alpha}_{d-}\rho f^2 )(Q)r^{n-1}dr
 =\frac{(n-1)!}{ \Gamma(n-\alpha)\Gamma(\alpha)} \int_0^{d}
 (\mathfrak{D }^{\alpha}_{d-}\rho f^2 )(Q) dr
\int_0^{r}  \frac{t^{n-1-\alpha}}{(r-t)^{1-\alpha}} dt .
$$
Note that by theorems  \ref{T3} and  \ref{T4} we have
$\rho f^2= \mathfrak{I}^{\alpha}_{d-}(\mathfrak{D }^{\alpha}_{d-}\rho f^2)$.
Using   Fubini's theorem, we obtain
\begin{align*}
&\frac{(n-1)!}{ \Gamma(n-\alpha)\Gamma(\alpha)} \int_0^{d}
 (\mathfrak{D }^{\alpha}_{d-}\rho f^2 )(Q) dr
\int_0^{r}  \frac{t^{n-1-\alpha}}{(r-t)^{1-\alpha}} dt \\
& = \frac{(n-1)!}{ \Gamma(n-\alpha)}\int_0^{d}
  \big[\mathfrak{I}^{\alpha}_{d-}(\mathfrak{D }^{\alpha}_{d-}
 \rho f^2)\big](P+\vec{e}t)    t^{n-1-\alpha} dt \\
&= C_n^{(\alpha)}\int_0^{d}
(\rho f^2) (Q)    r^{n-1-\alpha} dr.
\end{align*}
Therefore,
$$
I  =  \frac{1}{2\Gamma(1-\alpha)}\int_0^{d  } (\rho f^2)(Q)r^{n-1} (d  -r)^{-\alpha}
dr
\geq \frac{\mathfrak{d}^{-\alpha }}{2\Gamma(1-\alpha)}\int_0^{d  }
(\rho f^2)(Q) r^{n-1} dr.
$$
Finally  for any direction $\vec{e}$, we obtain the inequality
$$
\int_0^{d (\vec{e} ) }
 f(Q)( \mathfrak{D}^{\alpha}f )(Q)\rho(Q) r^{n-1}dr
\geq\frac{\mathfrak{d}^{ -\alpha }}{
2\Gamma(1-\alpha) }\int_0^{d (\vec{e} )}
(\rho f^2)(Q)r^{n-1} dr.
$$
Integrating the left and right sides of the last inequality we obtain
\begin{equation}\label{30}
\langle f,\mathfrak{D}^{\alpha}f\rangle_{L_2(\Omega,\rho)}
\geq\frac{1}{\lambda^2 } \|f\|^2_{L_2(\Omega,\rho )},\quad
f\in C^{\infty}_0(\Omega),\quad
\lambda^2 = 2\Gamma(1-\alpha) \mathfrak{d}^{\alpha}.
\end{equation}
 Suppose   that   $f\in H_0^{1} (\Omega)$.
There is a sequence $\{f_k\}\in C^{\infty}_0(\Omega)$ such that
$f_k\stackrel {H^1 }{\longrightarrow}f$.
The conditions imposed on the weight function $\rho$ implies the equivalence
of norms $L_2(\Omega)$ and $L_2(\Omega,\rho)$, hence
$f_k\stackrel {L_2(\Omega,\rho) }{\longrightarrow}f$.
Using the smoothness of weight function $\rho$, the embedding of spaces
$L_p( \Omega )$, $p\geq1$, and the inequality \eqref{3},
 we obtain the  estimate
\[
\|\mathfrak{D}^{\alpha} f \|_{L_2(\Omega,\rho)}
\leq C_1 \|\mathfrak{D}^{\alpha} f \|_{L_{q}(\Omega)}
\leq C_2 \| f \|^2_{H_0 ^1(\Omega)},
\]
where $2<q<2n/(2\alpha-2+n)$, $C_i>0$, $(i= 1,2 )$.
Therefore
$ \mathfrak{D}^{\alpha}f_k\stackrel {L_2(\Omega,\rho) }{\longrightarrow}
\mathfrak{D}^{\alpha}f$.
Hence from the continuity properties of the inner product in the Hilbert space,
we obtain
$$
\langle f_k,\mathfrak{D}^{\alpha}f_k\rangle_{L_2(\Omega,\rho)}
\to \langle f ,\mathfrak{D}^{\alpha}f \rangle_{L_2(\Omega,\rho)} .
$$
Passing to the limit in the left and right side of inequality \eqref{30},
we obtain the inequality \eqref{28} in the real case.

Now consider the case when $f$ is complex-valued. Note that
\begin{gather}\label{31}
\operatorname{Re} \langle f,\mathfrak{D}^{\alpha}f\rangle_{L_2(\Omega,\rho)}
= \langle u,\mathfrak{D}^{\alpha}u\rangle_{L_2(\Omega,\rho)}+
\langle v,\mathfrak{D}^{\alpha}v\rangle_{L_2(\Omega,\rho)},\\
\label{32}
 u=\operatorname{Re}f,\,v=\operatorname{Im}f.
\end{gather}
Obviously  inequality \eqref{28} follows from relations \eqref{30} and \eqref{31}.
 \end{proof}

 \section{Sectorial property}

Consider a uniformly elliptic operator with real-valued
coefficients and fractional derivative in the sense of Kipriyanov  in
 lower terms,  defined by the expression
\begin{gather}
 Lu:=-  D_{j} ( a^{ij} D_iu)  +p\, \mathfrak{D}^{ \alpha }u,\quad
 (i,j=\overline{1,n}) ,  \nonumber \\
 \operatorname{D}(L)=H^2(\Omega)\cap H^{1}_0(\Omega), \nonumber\\
\label{33}
 a^{ij}(Q)\in C^{1}(\bar{\Omega})  ,\quad
a^{ij}\xi _i  \xi _{j}  \geq a_0 |\xi|^2 ,\quad a_0>0, \\
\label{34}
p(Q)>0,\quad p(Q)\in {\rm Lip\,\mu},\quad (0<\alpha<\mu).
\end{gather}
We also will consider the formal conjugate  operator
\[
 L^{+}u:=-  D_i ( a^{ij} D_{j}u)  + \mathfrak{D}  ^{\alpha}_{ d -}pu  , \quad
\operatorname{D}(L^{+})=\operatorname{D}(L),
\]
 and  the    operator
\begin{equation*}
  H=\frac{1}{2}(L+L^{+}).
 \end{equation*}
We will use the special case of the Green's formula
\begin{equation*}
-\int_{\Omega}D_{j}(a^{ij}D_iu)\,\bar{v}\, dQ
=\int_{\Omega}a^{ij}D_iu\, \overline{D_{j}v}\,  dQ\,,\quad
u\in H^2(\Omega),\; v\in H_0^{1}(\Omega) .
\end{equation*}
The following Lemma establishes a property of the closure of operator $L$.

\begin{lemma}\label{L3}
  The operators $L,L^{+},H$ have  closure
$\tilde{L},\tilde{L}^{+},\tilde{H}$, the domains of definition of  this  operators
is included  in  $  H_0 ^{1} (\Omega)$.
\end{lemma}

\begin{proof}
 At first consider operator $L$. For a function
$f\in \operatorname{D} (L) $ using the Green's formula,
 we have with the notations \eqref{32}
 \begin{equation}\label{35}
\begin{aligned}
-\int_{\Omega} D_{j} ( a^{ij} D_if)\bar{f}dQ
&=\int_{\Omega}    a^{ij} D_if \overline{D_{j} f }dQ \\
&=\int_{\Omega}    a^{ij} (D_iu  D_{j} u  +D_iv  D_{j} v  )dQ \\
&\quad + \imath \int_{\Omega}    a^{ij} (D_iv  D_{j} u  -D_iu  D_{j} v  )dQ.
\end{aligned}
 \end{equation}
Applying condition \eqref{33}, we obtain
\begin{equation}\label{36}
\begin{aligned}
\operatorname{Re} (  a^{ij} D_if ,D_{j}f)_{L_2(\Omega)}
&=  (  a^{ij} D_iu ,D_{j}u)_{L_2(\Omega)}
+ (  a^{ij} D_iv ,D_{j}v)_{L_2(\Omega)}  \\
&\geq a_0(\|   u \|^2_{L^{1}_2(\Omega)}+ \|   v \|^2_{L^{1}_2(\Omega)}) \\
& =a_0\|   f  \|^2_{L^{1}_2(\Omega)},\quad f\in H_0^{1}(\Omega).
\end{aligned}
\end{equation}
 From \eqref{35}, \eqref{36} it follows that
\begin{equation}\label{37}
-\operatorname{Re} ( D_{j} [a^{ij} D_if] ,f)_{L_2(\Omega)}\geq a_0\
|   f  \|^2_{L^{1}_2(\Omega)},\quad f\in \operatorname{D}(L).
\end{equation}
Choose an arbitrary $\varepsilon>0$. From  \eqref{37},
 theorem \ref{T5}, \eqref{34}, using a Jung's inequality it is easy to
show that for sequence  $ \{ f_n\} \subset \operatorname{D}(L)$,
 the next two-sided estimate holds
\begin{equation}\label{38}
\begin{aligned}
 C_1\|   f_n \|^2_{L^{1}_2(\Omega)}+ C_2\| f_n \|^2_{L_2(\Omega)}
&\leq   2 \,\operatorname{Re}(f_n,Lf_n)_{L_2(\Omega)} \\
&\leq\frac{1}{\varepsilon} \|Lf_n \|^2_{L_2(\Omega)}
+ \varepsilon \| f_n \|^2_{L_2(\Omega)},
\end{aligned}
\end{equation}
where $C_i>0$, $i=1,2$.
In consequence of  \cite[theorem 3.4 p. 337]  {kato1966},
from  lower estimate \eqref{38},   follows  that   operator $L$ has a closure.
Let $u\in \operatorname{D}(\tilde{L})$, then by definition,   exists a
sequence $\{f_n\}\subset \operatorname{D}(L) $ such as that
$ f_n\stackrel {L_2  }{\longrightarrow}  f$,
$\{Lf_n\}$ is fundamental sequence  in the sense of the  norm  $L_2(\Omega)$.
Hence the inequality \eqref{38} implies that a sequence $\{ f_n\}$
is fundamental  in the sense of the norm $H^{1}_0(\Omega)$.
  In consequence of the completeness property of the space
$H^{1}_0  (\Omega)$, we have  $u\in H^{1}_0  (\Omega)$.
The proof corresponding to the case of operator $L$  completed.
 For proving this result for the case of operators $L^{+},H$ we must note that
\begin{equation}\label{39}
\begin{gathered}
( L f ,g)_{L_2(\Omega)}=( f ,L^{+}g )_{L_2(\Omega)},\quad
\operatorname{Re}( Lf ,f)_{L_2(\Omega)}=( Hf , f )_{L_2(\Omega)},\\
f,g  \in \operatorname{D}(L),
\end{gathered}
\end{equation}
 and then repeat the previous proof using this remark.
\end{proof}

We have the following theorem describing the spectral properties
of the closed operator $\tilde{L}$.

\begin{theorem} \label{T6}
Operators $\tilde{L},\,\tilde{L}^{+}$  is strongly accretive,
the numerical range of values  belongs to  the sector
\begin{equation*}
   \mathfrak{S}:= \{\zeta\in\mathbb{C}: \arg (\zeta-\gamma)|\leq\theta\},
\end{equation*}
 where $\theta$ and $\gamma$ defined by the coefficients of operator $L$.
\end{theorem}

\begin{proof}
Consider operator $L$. Applying estimate  \eqref{37} we can get the
inequality
\begin{equation}\label{40}
  \operatorname{Re}  \langle f_n, L f_n  \rangle_{L_2(\Omega)}
\geq a_0\|   f_n \|^2_{L^{1}_2(\Omega)}
+ \operatorname{Re}\langle f_n,\mathfrak{D}^{\alpha}f_n\rangle_{L_2(\Omega,p)} ,\quad
\{f_n\}\subset\operatorname{D}(L).
\end{equation}
Assume $ f\in \operatorname{D}(\tilde{L})$. Note that exists sequence
$f_n\xrightarrow[L]{}f$, $\{f_n\}\subset\operatorname{D}(L)$
and in consequences  of lemma \ref{L2}: $f\in H_0^{1}(\Omega)$.
Using the continuity property  of the inner product and passing to the
limit in the left and right side of inequality \eqref{40},  we obtain
\begin{equation}\label{41}
 \operatorname{Re}\langle f , \tilde{L} f   \rangle_{L_2(\Omega)}
\geq a_0\|   f  \|^2_{L^{1}_2(\Omega)}+ \operatorname{Re}
\langle f ,\mathfrak{D}^{\alpha}f \rangle_{L_2(\Omega,p)} ,\quad
 f  \in\operatorname{D}(\tilde{L}).
\end{equation}
By theorem \ref{T5} we can rewrite the previous inequality in the  form
 \begin{equation}\label{42}
  \operatorname{Re}\langle f,\tilde{L}f
 \rangle_{L_2(\Omega)}\geq a_0\|f\|^2_{L^{1}_2(\Omega)}
+\frac{1}{ \lambda^2}\|f\|^2_{L_2(\Omega,p)} ,\quad f\in \operatorname{D}(\tilde{L}).
\end{equation}
Applying the   inequality of   Friedrichs - Poincare to the first summand
of   the right side  \eqref{42}, we will get  the inequality of a
strong accretiveness for operator $\tilde{L}$,
  \begin{equation}\label{43}
  \operatorname{Re}  \langle f,\tilde{L}f  \rangle_{L_2(\Omega)}
\geq \frac{1}{\mu^2}\|f\|^2_{L_2(\Omega)},\quad
 f  \in\operatorname{D}(\tilde{L}),\quad
 \mu^{-2}=a_0+\lambda^{-2}\inf_{Q\in \Omega}|p(Q)| .
 \end{equation}
Consider imaginary component  of the form generated by the operator $L$.
For $f\in \operatorname{D}(L)$ we obtain
\begin{equation}\label{44}
\begin{aligned}
|\operatorname{Im} \langle f,Lf    \rangle_{L_2(\Omega)}|
&\leq   \Big|\int_{\Omega} (a^{ij}  D_iu D_{j}v-a^{ij}  D_iv D_{j}u)dQ\Big|\\
&\quad +\big|  \langle u,\mathfrak{D}^{\alpha}v \rangle_{L_2(\Omega,p)}
-\langle v,\mathfrak{D}^{\alpha}u   \rangle_{L_2(\Omega,p)}\big| \\
&= I_1+I_2.
\end{aligned}
 \end{equation}
 Using the Cauchy-Schwarz inequality for a sum, then the Jung's inequality,
 we have
\begin{equation}\label{45}
\begin{gathered}
a^{ij}  D_iu D_{j}v \leq a(Q) |Du| |Dv|
 \leq  \frac{1}{2} a(Q)(|Du|^2 +|Dv|^2 ),\\
a(Q)=\Big(\sum_{i,j=1}^{n}|a_{ij}(Q)|^2 \Big)^{1/2}.
\end{gathered}
\end{equation}
Hence
\begin{equation*}
 I_1\leq   a_1\|f\|^2_{L^{1}_2(\Omega)},\;a_1=\sup_{Q\in \Omega}|a(Q)| .
\end{equation*}
 Applying  inequality \eqref{3}, and  Jung's inequality,   we obtain
\begin{gather}\label{46}
\begin{aligned}
&|\langle u,\mathfrak{D}^{\alpha}v \rangle_{L_2(\Omega,p)}| \\
&\leq C\|u\|_{L_2(\Omega)}\|\mathfrak{D}^{\alpha}v\|_{L_{q}(\Omega)}
 \leq C  \|u\|_{L_2(\Omega)}\Big\{\frac{K}{\delta^{\nu}}
  \|v\|_{L_2(\Omega)}+\delta^{1-\nu}\|v\|_{L^{1}_2(\Omega)} \Big\} \\
&\leq\frac{1}{\varepsilon} \|u\|^2_{L_2(\Omega)}
 + \varepsilon\big(\frac{ KC }{\sqrt{2}\delta^{\nu}}\big)^2 \|v\|^2_{L_2(\Omega)}
 +  \frac{\varepsilon}{2}( C \delta^{1-\nu})^2 \|v\|^2_{L^{1}_2(\Omega)}  ,
\end{aligned} \\
\label{47}
 2<q<\frac{2n}{2\alpha-2+n},\quad
C = (\operatorname{meas}\Omega)^{\frac{q-2}{q}}\sup_{Q\in\Omega}p(Q).
 \end{gather}
Hence
\begin{equation}\label{48}
\begin{aligned}
 I_2
&\leq\big|\langle u,\mathfrak{D}^{\alpha}v \rangle_{L_2(\Omega,p)}\big|
+\big|\langle v,\mathfrak{D}^{\alpha}u \rangle_{L_2(\Omega,p)}\big| \\
&\leq \frac{1}{\varepsilon}\big(\|u\|^2_{L_2(\Omega)}+\|v\|^2_{L_2(\Omega)}\big)
 +\varepsilon\big(\frac{ KC }{\sqrt{2}\delta^{\nu}}\big)^2
(\|u\|^2_{L_2(\Omega)}+\|v\|^2_{L_2(\Omega)} ) \\
&\quad +\frac{\varepsilon}{2}( C \delta^{1-\nu})^2(\|u\|^2_{L^{1}_2(\Omega)}
 +\|v\|^2_{L^{1}_2(\Omega)} ) \\
&= \big(\varepsilon \delta^{-2\nu}C_1  +\frac{1}{\varepsilon}\big)
 \|f\|^2_{L_2(\Omega)}
 + \varepsilon \delta^{2-2\nu} C_2 \|f\|^2_{L^{1}_2(\Omega)}  .
\end{aligned}
\end{equation}
Using \eqref{48}, \eqref{44} and applying a reasoning,
analogous to the one in  the proof of  \eqref{41}, we have the  estimate
\begin{align*}
&|\operatorname{Im} \langle f,\tilde{L}f    \rangle_{L_2(\Omega)}| \\
&\leq  (\varepsilon \delta^{-2\nu}C_1  +\frac{1}{\varepsilon})
 \|f\|^2_{L_2(\Omega)} + (\varepsilon \delta^{2-2\nu} C_2+a_1)
 \|  f\|^2_{L_2^{1}(\Omega)}, f\in \operatorname{D}(\tilde{L}).
\end{align*}
Thus  in consequence of \eqref{43}  for arbitrary $k>0$, the next
inequality holds
\begin{align*}
&\operatorname{Re}\langle f,\tilde{L}f  \rangle_{L_2(\Omega)}
 -k |\operatorname{Im} \langle f,\tilde{L}f  \rangle_{L_2(\Omega)}| \\
&\geq  \big(a_0-k    [\varepsilon\delta^{2-2\nu} C_2+a_1] \big)\|  f\|^2_{L_2^{1}(\Omega)}
 +\big( \frac{1}{\mu^2}- k[\varepsilon \delta^{-2\nu}C_1
 +\frac{1}{\varepsilon}] \big)\|f\|^2_{L_2(\Omega)}.
\end{align*}
 Choose  $k= a_0(\varepsilon \delta^{2-2\nu} C_2+ a_1 )^{-1}$, we obtain
\begin{equation}\label{49}
\begin{gathered}
|\operatorname{Im} \langle f,(\tilde{L}-\gamma) f \rangle_{L_2(\Omega)}|
 \leq \frac{1}{k}\operatorname{Re}\langle f,(\tilde{L}-\gamma)f
\rangle_{L_2(\Omega)},\\
\gamma= \frac{1}{\mu^2}- k[\varepsilon \delta^{-2\nu}C_1  +\frac{1}{\varepsilon}].
\end{gathered}
\end{equation}
The last inequality implies that the numerical range of values
$\Theta(\tilde{L})$ belongs to the sector with top in $\gamma$ and half-angle
$\theta=\arctan(1/k)$.
The prove for the case corresponding to the operator $\tilde{L}^{+}$
is obvious if we note the first relation \eqref{39}.
 \end{proof}

We will not research detailed conditions for the coefficients of operator
 $L$ under which  $\gamma>0$ holds, just we note that it follows
from the second  relation of \eqref{49}, that can be easy formulated.
In further reasoning , we assume that the coefficients $\delta,\varepsilon$
of the  operator $L$ are chosen according to the second relation of \eqref{49} so that
 $\gamma>0$.

\begin{theorem}\label{T7}
   The operators $\tilde{L},\tilde{L}^{+},\tilde{H}$ are $m$-sectorial, and
 operator $\tilde{H}$ is self-adjoint.
\end{theorem}

\begin{proof}
Let us prove the theorem for the case corresponding to the operator $\tilde{L}$.
By theorem \ref{T6} we know that the operator $\tilde{L}$ is sectorial i.e.
the numerical range of values  of  $\tilde{L}$   belongs to the  sector
 $\mathfrak{S}$. By \cite[ Theorem 3.2 p. 336]{kato1966}
 we arrive to the conclusion  that $\mathrm{R}(\tilde{L}-\zeta)$
is a closed  space  for any $\zeta\in \mathbb{C}\setminus\mathfrak{S}$
and the next relation  holds
\begin{equation}\label{50}
  \operatorname{def}(\tilde{L}-\zeta)=\mu,\; \mu={\rm const} .
 \end{equation}
Since \eqref{43} holds, then on the subspace
$\mathrm{R}(\tilde{L}+\zeta)$, $\operatorname{Re}\zeta>0$  defined the inverse
operator.     Accordingly  condition  \cite[(3.38) p.350]{kato1966},
 we need to show  that
\begin{equation}\label{51}
\operatorname{def}(\tilde{L}+\zeta)=0,\quad
\|(\tilde{L}+\zeta)^{-1}\|\leq (\operatorname{Re}\zeta)^{-1},\quad
\operatorname{Re}\zeta>0 .
\end{equation}
Let us prove the first relation of \eqref{51}.
Suppose that the parameters: $\delta,\varepsilon$   are chosen from the second
relation  \eqref{49} so  that $\gamma>0$.
   Hence, the left half-plane has included in  to the complement of the sector
$\mathfrak{S}$ in the complex plain.

Let $\zeta_0\in\mathbb{C}\setminus \mathfrak{S}$, $\operatorname{Re}\zeta_0 <0$.
In consequence of inequality  \eqref{43}, we have
 \begin{equation}\label{52}
  \operatorname{Re}  \langle f,(\tilde{L}-\zeta )f  \rangle_{L_2(\Omega)}
\geq  (\mu^{-2}- \operatorname{Re} \zeta ) \|f\|^2_{L_2(\Omega)} .
 \end{equation}
Since the  operator $\tilde{L}-\zeta_0$ has a closed range of values
$\mathrm{R} (\tilde{L}-\zeta_0)$, it follows that
\begin{equation*}
 L_2=\mathrm{R} (\tilde{L}-\zeta_0)\oplus \mathrm{R} (\tilde{L}-\zeta_0)^{\perp} .
 \end{equation*}
Note that intersection of sets  $  C^{\infty}_0(\Omega)$ and
$\mathrm{R} (\tilde{L}-\zeta_0)^{\perp}$ is  empty, because if we assume
otherwise, then applying inequality  \eqref{52}  for any element
 $u\in C^{\infty}_0(\Omega)\cap \mathrm{R}  (\tilde{L}-\zeta_0)^{\perp}$, we obtain
 \begin{equation*}
(\mu^{-2}-\operatorname{Re}\zeta_0) \|u\|^2_{L_2(\Omega)}
 \leq \operatorname{Re} \langle u,(\tilde{L}-\zeta_0)u  \rangle_{L_2(\Omega)}=0,
 \end{equation*}
hence $u=0$. Thus  intersection of sets  $  C^{\infty}_0(\Omega)$ and
$\mathrm{R} (\tilde{L}-\zeta_0)^{\perp}$ is empty, it implies that
$$
(g,v)_{L_2(\Omega)}=0,\quad \forall g\in  \mathrm{R}  (\tilde{L}-\zeta_0)^{\perp},\;
\forall v\in C^{\infty}_0(\Omega).
$$
Since $  C^{\infty}_0(\Omega)$  is dense set in $L_2(\Omega)$,
it follows that  $\mathrm{R}  (\tilde{L}-\zeta_0)^{\perp}=0$.
It implies that  $\operatorname{def} (\tilde{L}-\zeta_0) =0$ and if we
note \eqref{50} we arrive  to the conclusion that
$\operatorname{def} (\tilde{L}-\zeta )=0$,
$\zeta\in \mathbb{C}\setminus\mathfrak{S}$. Hence
$\operatorname{def} (\tilde{L}+\zeta )=0$,
$\operatorname{Re}\zeta>0$ and proof the first relation of  \eqref{51} is complete.

For proving the second relation \eqref{51} we  note that
 \begin{gather*}
(\mu^{-2}+\operatorname{Re}\zeta ) \|f\|^2_{L_2(\Omega)}
\leq \operatorname{Re} \langle f,(\tilde{L}+\zeta )f  \rangle_{L_2(\Omega)}
\leq \|f\|_{L_2(\Omega)}\|(\tilde{L}+\zeta )\|_{L_2(\Omega)}, \\
f\in \operatorname{D}(\tilde{L}),\quad \operatorname{Re}\zeta>0 .
\end{gather*}
Using the first relation of  \eqref{51} we have
 \begin{equation*}
\|(\tilde{L}+\zeta )^{-1}g\|_{L_2(\Omega)}
\leq(\mu^{-2}+\operatorname{Re}\zeta ) ^{-1} \|g\|_{L_2(\Omega)}
\leq ( \operatorname{Re}\zeta ) ^{-1} \|g\|_{L_2(\Omega)},\quad
g\in L_2(\Omega).
 \end{equation*}
This implies
$$
\|(\tilde{L}+\zeta )^{-1} \| \leq( \operatorname{Re}\,\zeta ) ^{-1},\quad
\operatorname{Re}\zeta>0.
$$
The proof of the case corresponding to  $\tilde{L}$ is complete.
The proof for the case corresponding to $\tilde{L}^{+}$ is analogous
 if we note the first relation of \eqref{39}.

Consider operator $\tilde{H}$. Obviously that $ \tilde{H}  $ is symmetric
operator. Hence the numerical range of values of operator $ \tilde{H}  $
belongs to the real axis. Note \eqref{39},\eqref{40}, using $H$-convergence
and  passing to the limit, analogously  way of obtaining  \eqref{43}
 we can get the  inequality
 \begin{equation*}
    \langle f,\tilde{H}f  \rangle_{L_2(\Omega)}
\geq \frac{1}{\mu^2}\|f\|^2_{L_2(\Omega)} .
 \end{equation*}
Reasoning as in the proof corresponding to the case of operator $\tilde{L} $
and applying  \cite[Theorem 3.2 p.336]{kato1966},
we conclude that  $\operatorname{def}(\tilde{H}-\zeta )=0$,
$\operatorname{Im}\zeta\neq 0 $ and
\begin{equation*}
\operatorname{def}(\tilde{H}+\zeta)=0,\quad
\|(\tilde{H}+\zeta)^{-1}\|\leq (\operatorname{Re}\zeta)^{-1},\quad
\operatorname{Re}\zeta>0 .
\end{equation*}
The last relations implies that operator $\tilde{H}$ is $m$-accretive.
as well as sectorial, then $m$-sectorial.
In consequence of   \cite[Theorem 3.16 p.340]{kato1966}
operator $\tilde{H}$ is self-adjoint.
 \end{proof}


\section{Main theorems}

Further reasoning would require the use  of  theory  sesquilinear sectorial forms.
If not stated, otherwise, we use the definitions and  notation from
\cite {kato1966}.
Consider the forms
\begin{gather}\label{53}
\mathbf{t}[u,v] = \int_{\Omega} a^{ij}D_iu \overline{D_{j}v}dQ
+ \int_{\Omega}  p\,\mathfrak{D}^{\alpha }u \, \bar{v } dQ, \;
 u,v\in H^{1}_0(\Omega) , \nonumber \\
\mathbf{t}^{*}[u,v] =\overline{\mathbf{t} [v,u]}
 =\int_{\Omega} a^{ij}D_{j}u \overline{D_iv}dQ
+ \int_{\Omega}   u  p\,  \overline{\mathfrak{D}^{\alpha }v }dQ, \nonumber \\
h= \frac{1}{2}(\mathbf{t}+\mathbf{t}^{*})=\mathfrak{Re} \mathbf{t}  . \nonumber
\end{gather}


 \begin{lemma}\label{L4}
The form $\mathbf{t}$ is a closed sectorial form, moreover
$\mathbf{t}=\mathfrak{\tilde{f}}$, where
 \begin{equation*}
 \mathfrak{ f }[u,v]=(\tilde{L}u,v)_{L_2},\quad
u,v\in \operatorname{D}(\tilde{L}).
 \end{equation*}
\end{lemma}

 \begin{proof}
At first we will prove that the following  two-sided  inequality   holds
\begin{equation}\label{54}
C_0\|f\|^2 _{H^{1}_0}\leq  \left|\mathbf{t}[f ]\right|\leq C_1\|f\|^2 _{H^{1}_0},\quad
f\in H^{1}_0(\Omega).
\end{equation}
Note that from \eqref{36}, theorem \ref{T5}, \eqref{34},
 we obtain a lower estimate \eqref{54}
\begin{equation}\label{55}
C_0\|f\|^2 _{H^{1}_0}\leq   \operatorname{Re} \mathbf{t}[f ]
\leq\left|\mathbf{t}[f ]\right|,\quad f\in H^{1}_0 (\Omega).
\end{equation}
Applying \eqref{45},\eqref{46},  we obtain the upper estimate \eqref{54},
\begin{equation}\label{56}
|\mathbf{t}[f]|\leq\big|(a^{ij}D_if,D_jf)_{L_2(\Omega)}\big|
+\big|(p \,\mathfrak{D}^{\alpha} f, f)_{L_2(\Omega)}\big|
 \leq C_1\|f\|^2_{H^{1}_0(\Omega)},\quad f\in H^{1}_0 (\Omega).
\end{equation}
The proof of the estimates \eqref{54} is complete.
By definition of form  $\mathbf{t}$ \eqref{53},
we have $ H^{1}_0(\Omega)=\operatorname{D}(\mathbf{ t })\subset \operatorname{D}
( \mathbf{\tilde{t}}) $. If $f\in \operatorname{D}( \mathbf{\tilde{t}} )$
then there is a sequence  $  f_n\xrightarrow[\mathbf{ t } ]{ }f  $,
applying lower estimate  \eqref{54} we can conclude that
$  f_n\xrightarrow[  ]{ H^{1}_0}f $. Hence $f\in H^{1}_0(\Omega)$,
and $\operatorname{D}( \mathbf{\tilde{t}}) \subset H^{1}_0(\Omega)$.
It implies that $\operatorname{D}( \mathbf{\tilde{t}})
=\operatorname{D}( \mathbf{ t })$ and $\mathbf{t}$ is closed form.
Proof of the sectorial property is contained in the proof of theorem \ref{T6}.

Let us prove that $\mathbf{t}=\mathfrak{\tilde{f}}$.
At first we need to show that
\begin{equation}\label{57}
\mathfrak{ \mathfrak{f} }[u,v]=\mathbf{t}[u,v],\quad
u,v\in \operatorname{D}(\mathfrak{f}).
\end{equation}
 Using the Green's formula, we have
\begin{equation}\label{58}
( L u ,v )_{L_2}=\mathbf{t}[u ,v ],\quad u,v\in \operatorname{D}(L).
 \end{equation}
 Hence  we can rewrite the  relation \eqref{54} as
\begin{equation}\label{59}
C_0\|f\|^2 _{H^{1}_0}\leq  \left|( L f,f)_{L_2}\right|
\leq C_1\|f\|^2 _{H^{1}_0},\quad f\in \operatorname{D}(L).
\end{equation}
 Suppose $f \in \operatorname{D}(\tilde{L})$. Since conditions of $L$-convergence
holds then exists sequence  $\{f_n\}\subset \operatorname{D}(L)$
such that  we have  $f_n\xrightarrow[L]{}f $. From  \eqref{59}, \eqref{54}
follows that
 $f_n\xrightarrow[\mathbf{t}]{}f$. Now  consider the elements
$u,v\in \operatorname{D}(\tilde{L})$, as shown above   exists a sequences
$\{u_n\},\{v_n\}$ including in  $\operatorname{D}(L)$ and  $L$-converging,
$\mathbf{t}$-converging  to $u,v$ respectively. Passing to the limit in
the left and right side of \eqref{58}, we obtain \eqref{57}.
 Hence   from \eqref{57},\eqref{54} follows
\begin{equation}\label{60}
C_0\|f\|^2 _{H^{1}_0}
\leq  |\mathfrak{ f }[f ]|
\leq C_1\|f\|^2 _{H^{1}_0},\quad f\in \operatorname{D}(\mathfrak{f}).
\end{equation}

Now we are ready to prove that $\mathbf{t}=\mathfrak{\tilde{f}}$.
 Note that in consequence of theorem \ref{T6} operator $\tilde{L}$
is sectorial, hence by   theorem 1.27 \cite[p. 399] {kato1966},
 the form $\mathfrak{f}$ has a closure.
If $f_0\in\operatorname{D}(\mathfrak{\tilde{f}})$ then  we have
$f_n\xrightarrow[\mathfrak{ f }]{}f_0,\;\{f_n\}\subset \operatorname{D}(\tilde{L})$.
Applying the  estimates  \eqref{60},\eqref{54}  we can easy to see that
 $f_n\xrightarrow[  \mathbf{t} ]{}f_0 $ and since $\mathbf{t}$  is closed form,
it follows that  $f_0\in \operatorname{D}(\mathbf{t})$.
 On the other hand let $f_0\in \operatorname{D}(\mathbf{t})$.
Note that  the set $ C_0^{\infty}(\Omega)=\operatorname{D}(\mathbf{t})$
 is a kernel of form $\mathbf{t}$ in the sense of definition
\cite[p. 397] {kato1966}, it follows from \eqref{54}.
Hence exists sequence $\{f_n\}\subset \operatorname{D}(\mathbf{t})$,
$f_n\xrightarrow[  \mathfrak{  \mathbf{t}  } ]{}f_0$. Applying  estimates
\eqref{54}, \eqref{60}, we have $f_n\xrightarrow[  \mathfrak{  f  } ]{}f_0$,
it implies that  $f_0\in \operatorname{D}( \mathfrak{ \tilde{ f}  })$.
Now let $u,v\in \operatorname{D}(\mathfrak{\tilde{f}})$, using given above,
\cite[Theorem 1.17 p. 395] {kato1966} and   passing to the
 limits in the left and right side of inequality \eqref{57}, we obtain
\begin{equation*}
\mathfrak{ \mathfrak{\tilde{f}} }[u,v]=\mathbf{t}[u,v],\quad
u,v\in \operatorname{D}(\mathfrak{\tilde{f}}).
\end{equation*}
On the other hand  if $u,v\in  \operatorname{D}( \mathbf{t} )$, by analogous way,
 we obtain
\begin{equation*}
\mathbf{t}[u,v]=\mathfrak{ \mathfrak{\tilde{f}} }[u,v],\quad
u,v\in \operatorname{D}( \mathbf{t} ).
\end{equation*}
Hence  $\mathbf{t}=\mathfrak{ \tilde{f} }$.
\end{proof}

\begin{lemma}\label{L5}
The form h is a closed symmetric sectorial  form, moreover
$h=\mathfrak{\tilde{k}}$, where
 \begin{equation*}
 \mathfrak{ k }[u,v]=(\tilde{H}u,v)_{L_2},\quad
u,v\in \operatorname{D}(\tilde{H}).
 \end{equation*}
\end{lemma}

\begin{proof}
The implementation of the  symmetric  property  in the sense of
\cite[(1.5) p.387] {kato1966}   follows  from the definition of
form $h$. It is sufficient to note that
$$
h[u,v]=\frac{1}{2}(t[u,v]+ \overline{t[v,u]}  )
=\frac{1}{2}\overline{( t[v,u] +\overline{t[u,v]})}
=\overline{h[v,u]},\quad u,v\in \operatorname{D}(h).
$$
It is obvious  that
$h[f]= \operatorname{Re}\mathbf{t}[f]$.
Hence applying the lower  estimate of \eqref{55}, estimate \eqref{56},   we have
\begin{equation}\label{61}
C_0\|f\| _{H^{1}_0}\leq  h[f ] \leq C_1\|f\| _{H^{1}_0},\quad f\in H^{1}_0(\Omega).
\end{equation}
Using $h$-convergence, it is easy to see that from lower estimate of \eqref{61}
consequences  that form $h$ is closed. Proof of the sectorial property
contains in the proof of theorem \ref{T6}.

Let us prove that $h=\mathfrak{\tilde{k}}$. At first we need to show that
\begin{equation}\label{62}
\mathfrak{ \mathfrak{k} }[u,v]=h[u,v],\quad u,v\in \operatorname{D}(\mathfrak{k}).
\end{equation}
Note the definition of operator $H$, applying  the Green's formula,
lemmas \ref{L1} and \ref{L2}, we have
\begin{equation}\label{63}
( H u,v)_{L_2}=h[u,v],\quad u,v\in \operatorname{D}(H).
\end{equation}
Using the  previous equality we can rewrite estimate \eqref{61}  as follows
\begin{equation}\label{64}
C_0\|f\| _{H^{1}_0}\leq    ( H f,f)_{L_2} \leq C_1\|f\| _{H^{1}_0},\quad
f\in \operatorname{D}(H).
\end{equation}
Note that in consequence of lemma  \ref{L3} operator $H$ has a closure and
$\operatorname{D}(\mathfrak{k})\subset H_0^{1}(\Omega)$.
Suppose $f \in \operatorname{D}(\tilde{H})$, since conditions of
$H$-convergence holds, then exists $\{f_n\}\subset \operatorname{D}(H)$,
such that  we have  $f_n\xrightarrow[H]{}f $. From  \eqref{64} and\eqref{61} it
follows that  $f_n\xrightarrow[h]{}f$.
Now  consider the elements $u,v\in \operatorname{D}(\tilde{H})$, as shown
above exists sequences  $\{u_n\},\{v_n\}$ including in
$\operatorname{D}(H)$ and  $H$-converging, $h$-converging to $u,v$ respectively.
Passing to the limit on the left and right side of \eqref{63} we obtain \eqref{62}.
 From inequalities  \eqref{62} and \eqref{61} it follows
\begin{equation}\label{65}
C_0\|f\| _{H^{1}_0}\leq   \mathfrak{ k }[f ] \leq C_1\|f\| _{H^{1}_0},\quad
f\in \operatorname{D}(\tilde{H}).
\end{equation}
Now we are ready to prove that $h=\mathfrak{\tilde{   k  }}$.
Note that by theorem \ref{T6} operator $\tilde{H}$ is sectorial,
hence by    \cite[Theorem 1.27 p. 399] {kato1966} the form
$\mathfrak{k}$ has  a closure.
If $f_0\in\operatorname{D}(\mathfrak{\tilde{k}})$ then  we have
$f_n\xrightarrow[\mathfrak{ k }]{}f_0,\;\{f_n\}\subset \operatorname{D}(\tilde{H})$.
Applying the  estimates  \eqref{65} and \eqref{61}  we can easy to see that
$f_n\xrightarrow[  h ]{}f_0 $ and since $h$ is a closed form,
then $f_0\in \operatorname{D}(h)$. On the other hand let
 $f_0\in \operatorname{D}(h)$. Note that  the set
$ C_0^{\infty}(\Omega)=\operatorname{D}(h)$  is the  kernel of form $h$ in the sense
of definition  \cite[p. 397] {kato1966},
it follows from \eqref{61}. Hence exists sequence
$\{f_n\}\subset \operatorname{D}(h)$, $f_n\xrightarrow[     h    ]{}f_0$.
Applying  estimates \eqref{61},\eqref{65}, we have
$f_n\xrightarrow[  \mathfrak{  k  } ]{}f_0$, it implies that
$f_0\in \operatorname{D}( \mathfrak{ \tilde{ k}  })$.
Now let $u,v\in \operatorname{D}(\mathfrak{\tilde{k}})$, using given above,
  \cite[Theorem 1.17 p. 395] {kato1966} and   passing to the limits in the left and right side of inequality \eqref{62}, we obtain
\begin{equation*}
\mathfrak{ \mathfrak{\tilde{k}} }[u,v]=h[u,v],\;u,v\in \operatorname{D}(\mathfrak{\tilde{k}}).
\end{equation*}
On the other hand  if $u,v\in  \operatorname{D}( h )$, by an analogous way, we obtain
\begin{equation*}
h[u,v]=\mathfrak{ \mathfrak{\tilde{k}} }[u,v],\quad u,v\in \operatorname{D}( h ).
\end{equation*}
This implies $h=\mathfrak{ \tilde{k} }$.
\end{proof}

\begin{theorem}\label{T8}
Operator $\tilde{H}$ has a compact resolvent, discrete  spectrum,  and
the following estimate for  eigenvalues of operator $ \tilde{H}$ holds
\begin{equation}\label{66}
\lambda_n(L_0)\leq\lambda_n(\tilde{H})\leq\lambda_n(L_1),\,n\in\mathbb{N},
\end{equation}
where $\lambda_n(L_{k})$  are  respectively eigenvalues of  operators with
constant coefficients  defined by   operator  $L$
\begin{equation}\label{67}
L_{k}f=-  a_{ k }^{ij} D _{j}D_if +p_{k}f,\quad f\in \operatorname{D}(L),\quad k=0,1.
\end{equation}
\end{theorem}

\begin{proof}
First we need to prove three propositions.
\begin{itemize}
\item[(i)] Operator $\tilde{H}$ is positive defined. For proving this
fact note that operator $\tilde{H} $   has domain of definition is  dense in
$L_2(\Omega);$ the symmetric property of operator $\tilde{H}$ follows
from the definition; in consequently of theorem  \ref{T5}, $\tilde{H} $
is positive and  bounded below. Hence one is positive defined.

\item[(ii)] The space $H^{1}_0(\Omega)$ coincides as a set of elements,
with energetic spaces  $\mathfrak{H}_{\tilde{H}},\mathfrak{H}_{L_{k}}$,
$k=0,1$. We must note that
 \begin{equation}\label{68}
 \|f\| _{\mathfrak{H}_{\tilde{H}}}=h[f],\quad f\in H_0^{1}(\Omega),
\end{equation}
it follows from reasoning of  lemma \ref{L5}.

\item[(iii)] We have a following estimates for energetic   norms
\begin{equation}\label{69}
 \|f\| _{\mathfrak{H} _{L_0}}\leq \|f\| _{\mathfrak{H}_{\tilde{H}}}
\leq  \|f\| _{\mathfrak{H} _{L_1}},\,f\in H^{1}_0(\Omega).
\end{equation}
\end{itemize}
 Applying reasoning that was used to obtain the inequality \eqref{54},
we have equivalence of norms $H_0^{1}$ and $\mathfrak{H} _{L_{k}}$, $k=0,1$.
 In particular   we can easy to see that  exists   operators \eqref{67} $L_0,L_1$,
such  that the next  inequalities holds
\begin{equation}\label{70}
\|f\| _{\mathfrak{H} _{L_0}} \leq C_0\|f\| _{H^{1}_0},\quad
C_1\|f\| _{H^{1}_0}\leq\|f\| _{\mathfrak{H} _{L_1}},\quad
f\in H^{1}_0(\Omega).
\end{equation}
Hence  from \eqref{61}, \eqref{68} and \eqref{70}, it follows  \eqref{69}.

 Now we can prove the main statements of this theorem.
In (i)  we proved that the operators $ L_0 ,L_1,\tilde{H}$ is positive defined.
Note  given above, it is easy to see that the  norms
$H_0^{1}$, $\mathfrak{H}_{\tilde{H}}$, $\mathfrak{H} _{L_{k}}$, $k=0,1$
 are equivalence.   Applying  the  Rellich-Kondrashov  theorem we have
that the energetic  spaces:
  $\mathfrak{H}_{\tilde{H}},\;\mathfrak{H} _{L_{k}}$,  $k=0,1$
are compactly embedded in $L_2(\Omega)$. Using
\cite[Theorem 3 p.216]{mihlin1970} we have that  the operators
$ L_0 ,L_1,\tilde{H}$ has a discrete   spectrum.
In consequence of theorem \ref{T5} operator $\tilde{H}$
is semibounded from  below with positive constant, in consequence of
theorem \ref{T7} one is self-adjoint. Hence a zero belongs to a resolvent
set of operator $\tilde{H}$. In consequence of
 \cite[Theorem 5 p.222]{mihlin1970}  operator
$\tilde{H}$ has a compact resolvent at point zero.
Hence  in accordance with    \cite[Theorem 6.29 p.237]{kato1966}
operator $\tilde{H}$ has a compact resolvent on the resolvent set.

Using  (i), (ii), (iii),   we have
$$
L_0\leq \tilde{H} \leq L_1,
$$
where order relation is understood in terms of \cite[p.225]{mihlin1970}.
Since as mentioned above operators $ L_0 ,L_1,\tilde{H}$ has a discrete
 spectrum,  then  using    theorem 1 \cite[p.225]{mihlin1970}
we obtain \eqref{66}.
\end{proof}

\begin{theorem}
Operator $\tilde{L}$ has a compact resolvent, and  discrete spectrum.
\end{theorem}

\begin{proof}
Note that by theorem \ref{T7} operators $\tilde{L},\tilde{H}$ is $m$-sectorial,
operator $\tilde{H}$  is self-adjoint. Applying lemma \ref{L4}, lemma \ref{L5},
 \cite[Theorem 2.9, p.409]{kato1966} we have that
 $T_{\mathbf{t}}=\tilde{L}$, $T_{h}=\tilde{H}$, where $T_{\mathbf{t}}$, $T_{h}$
are respectively extensions by Fridrichs of operators $\tilde{L},\;\tilde{H}$
(see \cite[p.409]{kato1966}). Since in accordance
with   definition \cite[p.424]{kato1966} operator
$\tilde{H}$ is a real part of operator $\tilde{L}$, then in consequences
of Theorem \ref{T8} and  \cite[Theorem 3.3 p.424]{kato1966}
operator $\tilde{L}$ has a compact resolvent, hence by
  \cite[Theorem 6.29, p.237]{kato1966}  spectrum of
operator $\tilde{L}$  is discrete.
\end{proof}

\begin{remark} \rm
It is  easy to see that the Kypriaynov operator  in the one-dimensional
case reduces to the Marchaud operator, along with that  the results of
this work are true only  for    dimensions   $2\leq n<\infty$,
it follows from from the conditions of theorem \ref{T5}. However we can
 apply obtained  technique for the cases corresponding to the Marchaud
 operator  and Riemann-Liouville operator by using
 \cite[Corollary 1]{1kukushkin2017}, which  establishes the
strong accretive property  for these operators.
\end{remark}


\subsection*{Conclusions}

This paper presents results cocerning the field of spectral theory of
operators of fractional differentiation. It  proves a number of
propositions which  represent  independent interest
in the theory of  fractional calculus.
It also introduces a new construction of multidimensional fractional integral
 in a direction.
 Having formulated the sufficient conditions of representability
functions by the fractional integral in  a  direction,
in particular it proves the embedding of a Sobolev space in classes
of functions representable by the fractional integral in direction.

Note that the technique of proof borrowed from the one-dimensional
case is of particular interest. It should be noted that  the
constructed extension of Kipriyanov operator, was  found  as a
conjugate operator. These all create   a complete picture reflecting
the qualitative properties of fractional differential operators.
 It should be noted that in as main new results, there
were proven the following  theorems:
the theorem establishing a  strong accretive  property  for the operator
of fractional differentiation  in the Kyprianov sense,
the theorem establishing a sectorial  property for
differential operator  second order with   operator of fractional differentiation
 in   lower terms.
For this operator it was proven the theorem establishing the
 maximum accretive  property. Proven a theorem on the discreteness of the
spectrum of real part of operator, obtained two-sided estimate of its
eigenvalues. As the main result  a theorem on the discreteness of the spectrum of
 differential operator  second order with fractional derivative in the lower
terms was proven.
With the help of the theory of bilinear forms we
obtained general theoretical results for differential operators  second order
with fractional derivative in  lower terms.
In this paper we consider  the proof corresponding to the multidimensional case,
 however, is possible to have reduction to the one-dimensional case.
For example the one-dimensional  case was described in
\cite{2kukushkin2017}.
It should also be noted that the results in this direction can be treated
on the real axis. It is noteworthy that the use of bilinear forms
as a tool to study differential operator  second order with fractional
derivative in  lower terms, this  gives    the opportunity to see a dominant
 of senior term in the manifestation  of the functional properties of operator.
This technique is new and can  be used for  study the spectrum of the perturbed
operator of fractional  differentiation.  Therefore, regardless of the results
 the idea of the  proof may be of interest.

\subsection*{Acknowledgments}
The author wish to thank Professor Alexander L. Skubachev-skii for the
valuable remarks and comments made during the presentation of this  work
that took place 31/10/2017 at Peoples' Friendship University of Russia,
 Moscow.

\begin{thebibliography}{99}


\bibitem{1Aleroev1984}  T. S. Aleroev;
\emph{Spectral analysis of one class of nonselfadjoint operators},
  Differential Equations,  \textbf{20}, No. 1 (1984),  171-172.


 \bibitem{Aleroev1989} T. S.  Aleroev, B. I. Aleroev;
 \emph{On eigenfunctions and eigenvalues of one non-selfadjoint operator},
  Differential Equations,  \textbf{25}, No. 11 (1989),  1996-1997.


\bibitem{1Aleroev1994} T. S.  Aleroev;
 \emph{On eigenvalues of the one  class of nonselfadjoint operators},
  Differential Equations,  \textbf{30}, No. 1 (1994),  169-171.

\bibitem{3Aleroev2000} T. S.  Aleroev;
\emph{On the completeness of system of eigenfunctions of one dierential operator
of fractional order}, Differential Equations,  \textbf{36}, No. 6 (2000), 829-830.


\bibitem{fridrichs1958}  K. Friedrichs;
\emph{Symmetric positive linear differential equations},
 Comm. Pure Appl. Math.,   \textbf{11} (1958),  238-241.


\bibitem{1Jrbashyan1970}  M. M. Jrbashyan;
\emph{Boundary value problem for the differential operator  fractional order
type of  Sturm-Liouville},
  Proceedings of the Academy of Sciences of the Armenian SSR,
\textbf{5}, No.2 (1970), 37-47.

\bibitem{kato1961}  T. Kato;
\emph{Fractional powers of dissipative operators},
J.Math.Soc.Japan, \textbf{13}, No.3 (1961), 246-274.

\bibitem{kato1966} T. Kato;
\emph{Perturbation theory for linear operators},
Springer-Verlag Berlin, Heidelberg, New York, 1966.


 \bibitem{kipriyanov1960}  I. A. Kipriyanov;
  \emph{On spaces of fractionally differentiable functions},
Proceedings of the Academy Of Sciences. USSR,  \textbf{24} (1960),  665-882.

\bibitem{1kipriyanov1960} I. A. Kipriyanov;
\emph{The operator of fractional differentiation and the degree of elliptic
 operators},   Proceedings of the Academy of Sciences. USSR,  \textbf{131} (1960),     238-241.

 \bibitem{2.1kipriyanov1960} I. A. Kipriyanov;
\emph{On some properties of fractional derivative in the direction},
Proceedings of the universities. Math., USSR,  No.2 ( 1960),  32-40.

 \bibitem{2.2kipriyanov1960}  I. A. Kipriyanov;
\emph{On the complete continuity of embedding operators  in the  spaces of
fractionally differentiable functions},  Russian Mathematical Surveys,
\textbf{17} (1962),  183-189.

\bibitem{2kukushkin2017} M.V.Kukushkin;
\emph{Evaluation of the eigenvalues of the Sturm-Liouville problem for a
differential operator with fractional derivative in the lower terms},
  Belgorod State University Scientific Bulletin, Math. Physics.,
\textbf{46}, No.6 (2017),  29-35.

\bibitem{1kukushkin2017} M. V. Kukushkin;
\emph{Theorem about bounded embedding of the energetic space generated by the
operator of fractional differentition in the sense of Marchaud on the axis},
 Belgorod State University Scientific Bulletin, Math. Physics.,
\textbf{48}, No.20 (2017),  24-30.

 \bibitem{1.1kukushkin2017}   M. V. Kukushkin;
\emph{On some qulitative properties of the operator fractional differentiation
in Kipriyanov sense},  Vestnik of Samara  University, Natural Science Series,
 Math.,  \textbf{23}, No.2  (2017),  32-43.

\bibitem{Lyakhov2014}  L. N. Lyakhov, I. P. Polovinkin, E. L. Shishkina;
\emph{On a one problem of I. A. Kipriyanov  for a singular ultrahyperbolic equation},
  Differential Equations,  \textbf{50}, No.4 (2014)  513-525.

\bibitem{mihlin1970} S. G. Mihlin;
\emph{Variational methods in mathematical physics},  Moscow Science, 1970.

\bibitem{1Nakhushev1977}  A. M. Nakhushev;
\emph{The Sturm-Liouville problem for ordinary differential equation
second order with fractional derivatives in the lower terms},
  Proceedings of the Academy of Sciences. USSR,  \textbf{234}, No. 2 (1977),
 308-311.

\bibitem{1Reutskiy2016} S. Yu. Reutskiy;
\emph{A novel method for solving second order fractional eigenvalue problems},
 Journal of Computational and Applied Mathematics,  \textbf{306}, (2016),  133-153.

\bibitem{samko1987} S. G. Samko, A. A. Kilbas, O. I. Marichev;
\emph{Integrals and derivatives of fractional order and some of their applications},
 Minsk Science and technology, 1987.

\end{thebibliography}


\end{document}

