1 Introduction
Let d be some fixed positive integer number. By ${\mathbb{R}^{d}}$ we denote the real d-dimensional Euclidean space. As usual, we denote by $(\cdot ,\cdot )$ the inner product and by $|\cdot |$ the norm in ${\mathbb{R}^{d}}$ (we use the last notation for denoting the absolute value of a real number or a multi-index and the module of a complex number).
Let us consider a family of pseudo-differential operators ${(A(t,x))_{t\ge 0,x\in {\mathbb{R}^{d}}}}$ defined by the symbols ${(a(t,x,\lambda ))_{\lambda \in {\mathbb{R}^{d}}}}$ for every $t\ge 0$, $x\in {\mathbb{R}^{d}}$. That is,
\[ A(t,x)f(x)=\frac{1}{{(2\pi )^{d}}}{\int _{{\mathbb{R}^{d}}}}a(t,x,\lambda )F(\lambda ){e^{i(\lambda ,x)}}\hspace{0.1667em}\mathrm{d}\lambda ,\hspace{1em}t\ge 0,\hspace{2.5pt}x\in {\mathbb{R}^{d}},\]
where F is the Fourier transform of the function f: $F(\lambda )={\textstyle\int _{{\mathbb{R}^{d}}}}f(x){e^{-i(\lambda ,x)}}\hspace{0.1667em}\mathrm{d}x$, $\lambda \in {\mathbb{R}^{d}}$. Note that the function a may be complex-valued.We assume that the following conditions from [6, Ch. 4] about the function a and the operator A are satisfied (we denote by $Lf(r,\cdot )(x)$ the result of the operator L used on $f(r,x)$ as the function of x).
Theorem 4.3 from [6] (see also Theorem 4.1 there) states that there exists a bounded nonterminating strict Markov process ${(x(t))_{t\ge 0}}$ without second kind discontinuities, and the fundamental solution ${(g(s,x,t,y))_{0\le s\lt t,x,y\in {\mathbb{R}^{d}}}}$ to the equation
-
A1) The function a is homogeneous of a degree $\alpha \in (1,2)$ with respect to the variable λ and $\operatorname{Re}a(t,x,\lambda )\ge {a_{0}}\gt 0$ for all $t\ge 0$, $x,\lambda \in {\mathbb{R}^{d}}$, $|\lambda |=1$.
-
A2) The function a has $N\ge 2d+3$ (which is some natural number) continuous derivatives in λ for all $\lambda \ne 0$, and\[ |{\partial ^{\varkappa }}a(t,x,\cdot )(\lambda )|\le {C_{N}}|\lambda {|^{\alpha -|\varkappa |}},\]\[ |{\partial ^{\varkappa }}[a(t,x,\cdot )-a(s,y,\cdot )](\lambda )|\le {C_{N}}\left(|x-y{|^{\gamma }}+|t-s{|^{\gamma /\alpha }}\right)|\lambda {|^{\alpha -|\varkappa |}}\]for all multi-indexes ϰ with $|\varkappa |\le N$, $x,y,\lambda \in {\mathbb{R}^{d}}$, $\lambda \ne 0$, $s,t\in [0,+\infty )$. Here $\gamma \in (0,1)$ is some constant.
\[ \frac{\partial }{\partial t}u(t,x)+A(t,x)u(t,\cdot )(x)=0,\hspace{1em}t\gt 0,\hspace{2.5pt}x\in {\mathbb{R}^{d}},\]
is its transition probability density. The function g can be constructed by the parametrix method (see [6, Sec. 4.1.4]).If the function a is defined by the equality
\[ a(t,x,\lambda )=c|\lambda {|^{\alpha }},\hspace{1em}t\ge 0,\hspace{2.5pt}x\in {\mathbb{R}^{d}},\hspace{2.5pt}\lambda \in {\mathbb{R}^{d}},\]
where $c\gt 0$ is some constant, the corresponding Markov process is an isotropic α-stable process. The function g can be presented by the equality
\[ g(s,x,t,y)=\frac{1}{{(2\pi )^{d}}}{\int _{{\mathbb{R}^{d}}}}\exp \{i(x-y,\lambda )-c(t-s)|\lambda {|^{\alpha }}\}\mathrm{d}\lambda \]
in this case. The operator A (which does not depend on t and x) is the generator of this process. This is the simplest example of the processes considered in this paper. Therefore we name the process ${(x(t))_{t\ge 0}}$ the α-stable type process.Choosing some $0\lt \beta \le 1$, let us denote the β-order “pseudo-gradient” by ${\nabla _{\beta }}$, i.e. it is the operator which is defined by the symbol ${(i\lambda |\lambda {|^{\beta -1}})_{\lambda \in {\mathbb{R}^{d}}}}$. Note that ${\nabla _{\beta }}$ is a vector-valued operator. In the case of $\beta =1$, it is the ordinary gradient. We restrict ourselves to the case $\beta \ne 1$ in this paper.
Our goal is to construct a perturbation of the considered process using the operator $\left(b,{\nabla _{\beta }}\right)$, where b is some measurable ${\mathbb{R}^{d}}$-valued function. Under the perturbation of the Markov process ${(x(t))_{t\ge 0}}$ we understand the construction of a two-parameter family of operators ${({𝕋_{st}})_{0\le s\lt t}}$, defined on the set ${C_{b}}({\mathbb{R}^{d}})$ of real-valued continuous bounded functions defined on ${\mathbb{R}^{d}}$, such that for every $\varphi \in {C_{b}}({\mathbb{R}^{d}})$ and $t\gt 0$, the function $u(s,x,t)={𝕋_{st}}\varphi (x)$ satisfies in some sense the following Cauchy problem:
\[ \frac{\partial }{\partial s}u(s,x,t)+L(s,x)u(s,\cdot ,t)(x)=0,\hspace{1em}0\le s\lt t,\hspace{2.5pt}x\in {\mathbb{R}^{d}},\]
where $L(s,x)=A(s,x)+(b(s,x),{\nabla _{\beta }})$.The problem of perturbation of Markov processes was and remains in the focus of attention of many researchers. The diffusion case, that is, if $\alpha =2$, was considered by M. I. Portenko (see [14, 15] and the references there). The perturbation operators there are of the form $(b,\nabla )$, where the function b belongs to some ${L_{p}}$-space of functions or is some generalized function of the delta function type. In [19], the parabolic equation of the type $\nabla (A\nabla u)+B\nabla u-{u_{t}}=0$ was considered. The case $\alpha \lt 2$ was considered in the works of K. Bogdan and T. Jakubowski [2], who studied such a perturbation with the function b from a Kato class. In [7], T. Jakubowski considered a fundamental solution of the fractional diffusion equation (like written above) with a singular drift, i.e. the case of $\beta =1$ and some singular function b. In [8], an α-stable process was perturbed by the gradient operator with a multiplier b satisfying a certain integral space-time condition. S. I. Podolynny, M. I. Portenko [13] and M. I. Portenko [16–18] investigated such perturbations with the function b from ${L_{p}}$. Y. Maekawa and H. Miura were concerned with nonlocal parabolic equations in the presence of a divergence free drift term in [10]. J.-U. Loebus and M. I. Portenko [9] perturbed the infinitesimal operator of a one-dimensional symmetric α-stable process using the operator $(q{\delta _{0}},{\partial _{\alpha -1}})$, where ${\partial _{\alpha -1}}$ is some pseudo-differential operator of the order $\alpha -1$. Perturbation of an α-stable process by a fractional Laplacian operator was considered in [5]. The results for $\alpha \in (1,2)$ and perturbation operators of the type $(b,{\nabla _{\alpha -1}})$ can be found in [11, 12]. We studied the case of an α-stable process and perturbation operators $(b,{\nabla _{\beta }})$ with a ${\mathbb{R}^{d}}$-valued time independent function b from ${L_{p}}({\mathbb{R}^{d}})$ and $0\lt \beta \lt \alpha $ in [3]. In [4], the case of a generalized coefficient b is considered.
This article is structured as follows. The next section is devoted to some auxiliary facts. The perturbation equation is solved in Section 3. In Section 4, we study some properties of the corresponding two-parameter evolutionary family of operators. The last section is devoted to the construction of a solution to a corresponding Cauchy problem.
We will use the notation C for different constants, if their values are not important. If we need to emphasize the dependence of the constant C on parameter π, we will write ${C_{\pi }}$. We will use the notation $(a\vee b)$ for $\max (a,b)$ and $(a\wedge b)$ for $\min (a,b)$.
2 Auxiliaries
First of all let us note that if $0\lt \beta \lt 1$ the operator ${\nabla _{\beta }}$ can be represented in the integral form:
where ${n_{\beta }}=-{2^{-1}}{\pi ^{-(d-1)/2}}\Gamma \left(-\beta /2\right)\Gamma \left((d+\beta +1)/2\right)/\Gamma (-\beta )$ (here one have to use the equality $\Gamma (1+x)=x\Gamma (x)$ to expand the Euler gamma function Γ to negative noninteger arguments). Representation (1) is true if a function f is at least Lipschitz continuous and bounded. To obtain the normalizing factor ${n_{\beta }}$, it is sufficient to apply the operator ${\nabla _{\beta }}$ to the function ${f_{\lambda }}(x)={e^{i(\lambda ,x)}}$, $x\in {\mathbb{R}^{d}}$, for any $\lambda \in {\mathbb{R}^{d}}$. Moreover, this is the way to prove (1).
(1)
\[ {\nabla _{\beta }}f(x)={n_{\beta }}{\int _{{\mathbb{R}^{d}}}}(f(x+y)-f(x))\frac{y}{|y{|^{d+\beta +1}}}\hspace{0.1667em}\mathrm{d}y,\]Next, we will need some auxiliary statements. The following lemma is proved in [6] (see Lemma 1.11 there). It will be used frequently in this article.
Lemma 1.
The inequality ($0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$, remind that $1\lt \alpha \lt 2$)
holds with some constant $C\gt 0$ that depends only on d, α, k and l for all ϰ, λ, k, l, satisfying the inequalities $0\lt k\lt \alpha +\varkappa $, $0\lt l\lt \alpha +\lambda $. Here $B(\cdot ,\cdot )$ is the Euler beta function.
(2)
\[\begin{array}{r}\displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{\lambda /\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+l}}}\frac{{(t-\tau )^{\varkappa /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+k}}}\hspace{0.1667em}\mathrm{d}z\le \\ {} \displaystyle C\left[B\left(1+\frac{\varkappa -k}{\alpha },1+\frac{\lambda }{\alpha }\right)\frac{{(t-s)^{1+(\varkappa +\lambda -k)/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+l}}}\right.+\\ {} \displaystyle \left.B\left(1+\frac{\varkappa }{\alpha },1+\frac{\lambda -l}{\alpha }\right)\frac{{(t-s)^{1+(\varkappa +\lambda -l)/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+k}}}\right]\end{array}\]Remark 2.
Analyzing the proof of Lemma 1, one can see that the following inequality holds for all $0\le s\lt \tau \lt t$, $x,y\in {\mathbb{R}^{d}}$:
(3)
\[\begin{array}{r}\displaystyle {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{\lambda /\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+l}}}\frac{{(t-\tau )^{\varkappa /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+k}}}\hspace{0.1667em}\mathrm{d}z\le \\ {} \displaystyle C\left[\frac{{(\tau -s)^{\lambda /\alpha }}{(t-\tau )^{(\varkappa -k)/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+l}}}+\frac{{(\tau -s)^{(\lambda -l)/\alpha }}{(t-\tau )^{\varkappa /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+k}}}\right].\end{array}\]The following lemma can be obtained from the results of [6, Ch. 4]).
Lemma 2.
If the assumptions A1–A3 hold, then the transition probability density of the process ${(x(t))_{t\ge 0}}$ has the following properties ($0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$): the function g is continuously differentiable with respect to $x\in {\mathbb{R}^{d}}$ and
where ${\partial ^{k}}$ means some derivative of the order k;
Here, positive constants ${N_{k,T}}$ and ${N_{\beta ,T}}$ can be depended on T.
(4)
\[ |{\partial ^{k}}g(s,\cdot ,t,y)(x)|\le {N_{k,T}}\frac{{(t-s)^{1-(\gamma +k)/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}},\hspace{1em}k=0,1,\](5)
\[\begin{array}{r}\displaystyle |{\nabla _{\beta }}g(s,\cdot ,t,y)(x)|\le {N_{\beta ,T}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}\right.+\\ {} \displaystyle \left.\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right).\end{array}\]Proof.
Let us note that the function g can be constructed by the parametrix method (see [6, Th. 4.1]). That is, it can be presented by the equality
where
and the function Φ is determined by some integral equation.
Theorem 4.1 cited above states that the following inequalities hold:
The function ${g_{0}}$ satisfies the inequality (see [6, eq. (4.1.25)])
in which ${\partial ^{k}}$ means some derivative of the integer order $0\le k\le N-2d-1$ (the constant N is defined in assumption A2), ${N_{k}}\gt 0$ are some constants. Note that k can be at least 0, 1, or 2.
(7)
\[\begin{array}{r}\displaystyle |h(s,x,t,y)|\le \\ {} \displaystyle C\left(\frac{{(t-s)^{1+\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha }}}+\frac{t-s}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right),\end{array}\](9)
\[ |{\partial ^{k}}{g_{0}}(s,\cdot ,t,y)(x)|\le {N_{k}}\frac{t-s}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha +k}}}\]To prove inequality (4) for the case of $k=1$, we have to show the equality
(10)
\[\begin{array}{r}\displaystyle \partial {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{g_{0}}(s,\cdot ,\tau ,z)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)=\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\partial {g_{0}}(s,\cdot ,\tau ,z)(x)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z.\end{array}\]Due to (8) and (9),
\[\begin{array}{r}\displaystyle \left|\partial {g_{0}}(s,\cdot ,\tau ,z)(x)\Phi (\tau ,z,t,y)\right|\le \\ {} \displaystyle \frac{{N_{1}}(\tau -s)}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha +1}}}\frac{C}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\le \\ {} \displaystyle {N_{1}}{(\tau -s)^{-1-(d+1)/\alpha }}\frac{C}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}.\end{array}\]
The right hand side of the last inequality is integrable with respect to z on ${\mathbb{R}^{d}}$ for every fixed $0\le s\lt \tau \lt t$, $y\in {\mathbb{R}^{d}}$ and we can write down the equality
\[ \partial {\int _{{\mathbb{R}^{d}}}}{g_{0}}(s,\cdot ,\tau ,z)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)={\int _{{\mathbb{R}^{d}}}}\partial {g_{0}}(s,\cdot ,\tau ,z)(x)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z.\]
Inequality (3) leads us to the following ($0\le s\lt \tau \lt t\le T$, $x,y\in {\mathbb{R}^{d}}$) inequalities:
\[\begin{array}{r}\displaystyle \left|\partial {\int _{{\mathbb{R}^{d}}}}{g_{0}}(s,\cdot ,\tau ,z)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)\right|\le \\ {} \displaystyle C\left[\frac{(\tau -s){(t-\tau )^{-1+\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha +1}}}+\frac{{(\tau -s)^{-1/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right]\le \\ {} \displaystyle {C_{T}}\left((\tau -s){(t-\tau )^{-1+\gamma /\alpha }}+{(\tau -s)^{-1/\alpha }}\right){(t-s)^{-1-(d+1)/\alpha }}.\end{array}\]
The last expression is integrable with respect to τ on $(s,t)$, so we have (10).Using (2) we can obtain the following inequalities valid for $k=0,1$:
\[\begin{array}{r}\displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\left|{\partial ^{k}}{g_{0}}(s,\cdot ,\tau ,z)(x)\Phi (\tau ,z,t,y)\right|\mathrm{d}z\le \\ {} \displaystyle {C_{k}}{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{\tau -s}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha +k}}}\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\hspace{0.1667em}\mathrm{d}z\le \\ {} \displaystyle {C_{k}}\left(\frac{{(t-s)^{1-k/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha +k}}}+\frac{{(t-s)^{1+\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right)\le \\ {} \displaystyle {C_{k}}\frac{{(t-s)^{1-k/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}},\end{array}\]
where ${C_{k}}\gt 0$ are some constants. Therefore (use (9)),
\[\begin{array}{r}\displaystyle |{\partial ^{k}}g(s,\cdot ,t,y)(x)|\le ({N_{k}}\vee {C_{k}})\frac{{(t-s)^{1-k/\alpha }}}{{({(t-s)^{1-(k+\gamma )/\alpha }}+|y-x|)^{d+\alpha +-\gamma }}}\times \\ {} \displaystyle \left[\frac{{(t-s)^{(k+\gamma )/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{k+\gamma }}}+{(t-s)^{\gamma /\alpha }}\right]\le \\ {} \displaystyle {C_{k,T}}\frac{{(t-s)^{1-(\gamma +k)/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}},\end{array}\]
for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and every $T\gt 0$. Here the constants ${C_{k,T}}\gt 0$ depend on T.For proving inequality (5), we use representation (1) of the operator ${\nabla _{\beta }}$ (remind that $0\lt \beta \lt 1$). So, since $a(t,x,\cdot )$ is a homogeneous function and
valid for all $0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$ with some positive constant ${C_{\beta }}$.
\[\begin{array}{r}\displaystyle {\nabla _{\beta }}{g_{0}}(s,\cdot ,t,y)(x)=\\ {} \displaystyle \frac{i}{{(2\pi )^{d}}}{\int _{{\mathbb{R}^{d}}}}\lambda |\lambda {|^{\beta -1}}\exp \{i(x-y,\lambda )-a(t,y,{(t-s)^{1/\alpha }}\lambda )\}\mathrm{d}\lambda =\\ {} \displaystyle \frac{i}{{(2\pi )^{d}}{(t-s)^{(d+\beta )/\alpha }}}{\int _{{\mathbb{R}^{d}}}}\lambda |\lambda {|^{\beta -1}}\exp \left\{i\left(\frac{x-y}{{(t-s)^{1/\alpha }}},\lambda \right)-a(t,y,\lambda )\right\}\mathrm{d}\lambda ,\end{array}\]
the result of [6, Lemma 4.2] leads us to the inequality
(11)
\[ |{\nabla _{\beta }}{g_{0}}(s,\cdot ,t,y)(x)|\le {C_{\beta }}\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}\]Let us consider the function
\[ f(x)={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{g_{0}}(s,x,\tau ,z)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z,\hspace{1em}x\in {\mathbb{R}^{d}},\]
for every fixed $0\le s\lt t$, $y\in {\mathbb{R}^{d}}$. Using representation (1) we can write down the equality
\[ {\nabla _{\beta }}f(x)={k_{\beta }}{\int _{{\mathbb{R}^{d}}}}\frac{u}{|u{|^{d+\beta +1}}}(f(x+u)-f(x))\hspace{0.1667em}\mathrm{d}u={I_{1}}+{I_{2}},\]
where ${I_{1}}$ and ${I_{2}}$ are the integrals of the same integrand as on the left hand side but over ${B_{\varepsilon }}(0)=\{u\in {\mathbb{R}^{d}}:|u|\le \varepsilon \}$ and ${B_{\varepsilon }}{(0)^{C}}$, respectively.Since, for $u\in {B_{\varepsilon }}(0)$ and fixed $x\in {\mathbb{R}^{d}}$, $f(x+u)-f(x)=(\nabla f(x+\theta u),u)$, where $\theta =\theta (u)\in (0,1)$,
\[\begin{array}{r}\displaystyle {I_{1}}=\\ {} \displaystyle {k_{\beta }}{\int _{{B_{\varepsilon }}(0)}}\frac{u}{|u{|^{d+\beta +1}}}\left(u,\nabla {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{g_{0}}(s,\cdot ,\tau ,z)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x+\theta u)\right)\hspace{0.1667em}\mathrm{d}u=\\ {} \displaystyle {k_{\beta }}{\int _{{B_{\varepsilon }}(0)}}\frac{u}{|u{|^{d+\beta +1}}}{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\left(u,\nabla {g_{0}}(s,\cdot ,\tau ,z)(x+\theta u)\right)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}\hspace{0.1667em}\mathrm{d}u.\end{array}\]
Due to (8) and (9), the integrand here can be estimated as follows:
\[\begin{array}{r}\displaystyle \left|\frac{u}{|u{|^{d+\beta +1}}}\left(u,\nabla {g_{0}}(s,\cdot ,\tau ,z)(x+\theta u)\right)\Phi (\tau ,z,t,y)\right|\le \\ {} \displaystyle \frac{C}{|u{|^{d+\beta -1}}}\frac{\tau -s}{{({(\tau -s)^{1/\alpha }}+|z-x-\theta u|)^{d+\alpha +1}}}\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}.\end{array}\]
For fixed u the last function is integrable with respect to $(\tau ,z)$ on $(s,t)\times {\mathbb{R}^{d}}$ and the integral is not greater then
\[ \frac{C}{|u{|^{d+\beta -1}}}\frac{{(t-s)^{1-1/\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x-\theta u|)^{d+\alpha -\gamma }}}.\]
Let us take $\varepsilon \lt \frac{1}{2}{(t-s)^{1/2}}$. The obtained estimator is not greater then the following function (for fixed $0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$) integrable with respect to u on ${B_{\varepsilon }}(0)$:
Here we have used the obvious inequalities
The integrand in ${I_{2}}$ can be estimated as follows ($|u|\ge \varepsilon $):
\[\begin{array}{r}\displaystyle \left|\frac{u}{|u{|^{d+\beta +1}}}({g_{0}}(s,x+u,\tau ,z)-{g_{0}}(s,x,\tau ,z))\Phi (\tau ,z,t,y)\right|\le \\ {} \displaystyle \frac{C}{|u{|^{d+\beta }}}\left(\frac{\tau -s}{{({(\tau -s)^{1/\alpha }}+|z-x-u|)^{d+\alpha }}}+\frac{\tau -s}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha }}}\right)\times \\ {} \displaystyle \frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}.\end{array}\]
Integrating this function with respect to $(\tau ,z)$ on $(s,t)\times {\mathbb{R}^{d}}$ for every fixed $u\in {B_{\varepsilon }}{(0)^{C}}$ one can obtain the estimator
\[\begin{array}{r}\displaystyle \frac{C}{|u{|^{d+\beta }}}\left(\frac{t-s}{{({(t-s)^{1/\alpha }}+|y-x-u|)^{d+\alpha -\gamma }}}+\frac{t-s}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right)\le \\ {} \displaystyle \frac{C}{|u{|^{d+\beta }}}{(t-s)^{-(d-\gamma )/\alpha }},\end{array}\]
which is integrable with respect to u on ${B_{\varepsilon }}{(0)^{C}}$ for all fixed $0\le s\lt t$.Therefore, integrals ${I_{1}}$ and ${I_{2}}$ are absolutely convergent. So, using the Fubini theorem we can state that
\[\begin{array}{r}\displaystyle {\nabla _{\beta }}{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{g_{0}}(s,\cdot ,\tau ,z)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)=\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{\nabla _{\beta }}{g_{0}}(s,\cdot ,\tau ,z)(x)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z.\end{array}\]
Now, using (8), (11) and Lemma 1, one can obtain the estimate
\[\begin{array}{r}\displaystyle |{\nabla _{\beta }}h(s,\cdot ,t,y)(x)|\le C\left(\frac{{(t-s)^{\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}\right.\\ {} \displaystyle +\left.\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right),\end{array}\]
and inequality (5). The lemma is proved. □3 Solving of perturbation equation
Let us consider the perturbation equation ($0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$)
Formally applying the operator ${\nabla _{\beta }}$ to both sides of equation (12) with respect to the variable x we can search for its solution in the form
where the function ${(v(s,x,t,y))_{0\le s\lt t,x,y\in {\mathbb{R}^{d}}}}$ is a solution to the equation
in which
(12)
\[\begin{array}{r}\displaystyle G(s,x,t,y)=g(s,x,t,y)+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),{\nabla _{\beta }}G(\tau ,\cdot ,t,y)(z))\hspace{0.1667em}\mathrm{d}z.\end{array}\](13)
\[ G(s,x,t,y)=g(s,x,t,y)+{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z,\]Theorem 1.
Let assumptions A1–A3 hold and the ${\mathbb{R}^{d}}$-valued function ${(b(t,x))_{t\gt 0,x\in {\mathbb{R}^{d}}}}$ belong to ${L_{p}}([0,T]\times {\mathbb{R}^{d}})$ for every $T\gt 0$ with some $p\gt \frac{d+\alpha }{\alpha -1}$, including $p=+\infty $. Then the following statements are true: there exists a unique solution to equation (14) in the class of functions satisfying the estimate
the function (13) is a solution to equation (12) and satisfies the estimate
for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and every $T\gt 0$, where a positive constant ${C_{T}}$ can depend on T.
(16)
\[\begin{array}{r}\displaystyle |v(s,x,t,y)|\le \\ {} \displaystyle {C_{T}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right);\end{array}\](17)
\[\begin{array}{r}\displaystyle |G(s,x,t,y)|\le \\ {} \displaystyle {C_{T}}\left(\frac{{(t-s)^{\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{1-\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right)\end{array}\]Before proving this theorem, let us prove the following auxiliary statement.
Lemma 3.
If the assumptions of Theorem 1 hold and a ${\mathbb{R}^{d}}$-valued function ${(v(s,x,t,y))_{0\le s\lt t,x,y\in {\mathbb{R}^{d}}}}$ satisfies inequality (16), then the equality
is valid.
(18)
\[\begin{array}{r}\displaystyle {\nabla _{\beta }}\left({\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,\cdot ,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\right)(x)=\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{\nabla _{\beta }}g(s,\cdot ,\tau ,z)(x)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\end{array}\]Proof.
To prove this lemma we use calculations similar to those used in the proof of Lemma 2. First of all, let us prove the differentiability of the function
\[ f(x)={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z,\hspace{1em}x\in {\mathbb{R}^{d}},\]
for every fixed $0\le s\lt t\le T$, $y\in {\mathbb{R}^{d}}$. Since
\[\begin{array}{r}\displaystyle \left|\partial g(s,\cdot ,\tau ,z)(x)(b(\tau ,z),v(\tau ,z,t,y))\right|\le {C_{T}}\frac{{(\tau -s)^{1-(\gamma +1)/\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha -\gamma }}}\times \\ {} \displaystyle \left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right)|b(\tau ,z)|\end{array}\]
for all fixed $0\le s\lt \tau \lt t\le T$, $x,y,z\in {\mathbb{R}^{d}}$ we obtain that the last expression is not greater then
\[\begin{array}{r}\displaystyle {C_{T}}{(\tau -s)^{-(d+1)/\alpha }}|b(\tau ,z)|\times \\ {} \displaystyle \left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right).\end{array}\]
Moreover1 (we take $q=p/(p-1)$),
\[\begin{array}{r}\displaystyle {\int _{{\mathbb{R}^{d}}}}|b(\tau ,z)|\left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right)\hspace{0.1667em}\mathrm{d}z\le \\ {} \displaystyle C{\left({\int _{{\mathbb{R}^{d}}}}|b(\tau ,z){|^{p}}\hspace{0.1667em}\mathrm{d}z\right)^{1/p}}\left({\int _{{\mathbb{R}^{d}}}}\left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{q(d+\beta )}}}+\right.\right.\\ {} \displaystyle {\left.\left.\frac{{(t-\tau )^{q(1-\beta /\alpha )}}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{q(d+\alpha -\gamma )}}}\right)\mathrm{d}z\right)^{1/q}}\le \\ {} \displaystyle C{\left({\int _{{\mathbb{R}^{d}}}}|b(\tau ,z){|^{p}}\hspace{0.1667em}\mathrm{d}z\right)^{1/p}}{(t-\tau )^{(d-1)/q\alpha -(d+\beta )/\alpha }}\lt +\infty \end{array}\]
almost everywhere on $(s,t)$. Therefore, we have
\[\begin{array}{r}\displaystyle \partial {\int _{{\mathbb{R}^{d}}}}g(s,\cdot ,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)=\\ {} \displaystyle {\int _{{\mathbb{R}^{d}}}}\partial g(s,\cdot ,\tau ,z)(x)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\end{array}\]
and (see Remark 2)
\[\begin{array}{r}\displaystyle \left|\partial {\int _{{\mathbb{R}^{d}}}}g(s,\cdot ,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)\right|\le \\ {} \displaystyle C{\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{1-(\gamma +1)/\alpha }}|b(\tau ,z)|}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha -\gamma }}}\times \end{array}\]
\[\begin{array}{r}\displaystyle \left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right)\hspace{0.1667em}\mathrm{d}z\le \\ {} \displaystyle {C_{T}}{\left({\int _{{\mathbb{R}^{d}}}}|b(\tau ,z){|^{p}}\hspace{0.1667em}\mathrm{d}z\right)^{1/p}}{(\tau -s)^{-(d+p)/\alpha p}}{(t-\tau )^{-(d+\beta p)/\alpha p}}\times \\ {} \displaystyle \left(\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\vee \frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}\right)\le \\ {} \displaystyle {C_{T}}{\left({\int _{{\mathbb{R}^{d}}}}|b(\tau ,z){|^{p}}\hspace{0.1667em}\mathrm{d}z\right)^{1/p}}\times \\ {} \displaystyle {(\tau -s)^{-(d+p)/\alpha p}}{(t-\tau )^{-(d+\beta p)/\alpha p}}\left({(t-s)^{-(d+\alpha -\gamma )/\alpha }}\vee {(t-s)^{-(d+\beta )/\alpha }}\right).\end{array}\]
The last function is integrable with respect to τ on $(s,t)$ for every fixed $0\le s\lt t\le T$ in the assumptions of the lemma.Therefore,
\[\begin{array}{r}\displaystyle \partial {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,\cdot ,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)=\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau \hspace{0.1667em}\partial {\int _{{\mathbb{R}^{d}}}}g(s,\cdot ,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}(x)=\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\partial g(s,\cdot ,\tau ,z)(x)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z.\end{array}\]
Consequently, we can use representation (1) of the operator ${\nabla _{\beta }}$:
\[ {\nabla _{\beta }}f(x)={k_{\beta }}{\int _{{\mathbb{R}^{d}}}}\frac{u}{|u{|^{d+\beta +1}}}(f(x+u)-f(x))\hspace{0.1667em}\mathrm{d}u={I_{1}}+{I_{2}},\]
where ${I_{1}}$ and ${I_{2}}$ are the integrals of the same integrand as on the left hand side but over ${B_{\varepsilon }}(0)$ and ${B_{\varepsilon }}{(0)^{C}}$, respectively. Namely,
\[\begin{array}{r}\displaystyle {I_{1}}={k_{\beta }}{\int _{{B_{\varepsilon }}(0)}}\frac{u\hspace{0.1667em}\mathrm{d}u}{|u{|^{d+\beta +1}}}\times \\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}(g(s,x+u,\tau ,z)-g(s,x,\tau ,z))(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z=\\ {} \displaystyle {k_{\beta }}{\int _{{B_{\varepsilon }}(0)}}\frac{u\hspace{0.1667em}\mathrm{d}u}{|u{|^{d+\beta +1}}}{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}(\nabla g(s,\cdot ,\tau ,z)(x+\theta u),u)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z,\end{array}\]
where $\theta \in (0,1)$ is some constant dependent of s, t, x, y and u only;
\[\begin{array}{r}\displaystyle {I_{2}}={k_{\beta }}{\int _{{B_{\varepsilon }}{(0)^{C}}}}\frac{u}{|u{|^{d+\beta +1}}}\times \\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}(g(s,x+u,\tau ,z)-g(s,x,\tau ,z))(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\hspace{0.1667em}\mathrm{d}u.\end{array}\]
The integrand of ${I_{1}}$ satisfies the inequality
For the integral ${I_{2}}$ we have
The right hand sides of (19) and (20) are integrable on the domains of integrals ${I_{1}}$ and ${I_{2}}$, respectively. Thus, using the Fubini theorem, we obtain equality (18) and the lemma is proved. □
(19)
\[\begin{array}{r}\displaystyle \left|\frac{u}{|u{|^{d+\beta +1}}}(\nabla g(s,\cdot ,\tau ,z)(x+\theta u),u)(b(\tau ,z),v(\tau ,z,t,y))\right|\le \\ {} \displaystyle \frac{{C_{T}}}{|u{|^{d+\beta -1}}}\frac{{(\tau -s)^{1-(\gamma +1)/\alpha }}|b(\tau ,z)|}{{({(\tau -s)^{1/\alpha }}+|z-x-\theta u|)^{d+\alpha -\gamma }}}\times \\ {} \displaystyle \left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right).\end{array}\](20)
\[\begin{array}{r}\displaystyle \left|\frac{u}{|u{|^{d+\beta +1}}}(g(s,x+u,\tau ,z)-g(s,x,\tau ,z))(b(\tau ,z),v(\tau ,z,t,y))\right|\le \\ {} \displaystyle \frac{{C_{T}}|b(\tau ,z)|}{|u{|^{d+\beta }}}\times \\ {} \displaystyle \left(\frac{{(\tau -s)^{1-(\gamma +1)/\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x-u|)^{d+\alpha -\gamma }}}+\frac{{(\tau -s)^{1-(\gamma +1)/\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha -\gamma }}}\right)\times \\ {} \displaystyle \left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right).\end{array}\]Proof of Theorem 1.
Let us solve equation (14) for all fixed $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$, using the method of successive approximations. Namely, consider the sequence of functions ${({v_{k}}(s,x,t,y))_{0\le s\lt t\le T,x,y\in {\mathbb{R}^{d}}}}$, $k=0,1,2,\dots \hspace{0.1667em}$, given by the recurrence relation
\[ {v_{k+1}}(s,x,t,y)={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,\tau ,z)(b(\tau ,z),{v_{k}}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z,\]
where ${v_{0}}$ is defined by (15).Using inequality (5) one can obtain the relation ($0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$)
Iterating relation (21) by using inequality (2), we come to the estimate
where $\theta =1-((d+\alpha )/p+\beta )/\alpha $ (note that $\theta \gt (1-\beta )/\alpha $) and the sequence $\{{R_{k,T}}:k=0,1,2,\dots \hspace{0.1667em}\}$ satisfies the relations ${R_{0}}={N_{\beta ,T}}$ and
for all $k=1,2,\dots \hspace{0.1667em}$. Here C is the maximum of a finite number of constants taken from inequality (2). Indeed, when $k=0$, inequality (22) coincides with inequality (5). If (22) is correct for some $k\in \mathbb{N}$, we obtain
(21)
\[\begin{array}{r}\displaystyle |{v_{k+1}}(s,x,t,y)|\le \\ {} \displaystyle 2{N_{\beta ,T}}\| b{\| _{p}^{T}}\left({\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{|{v_{k}}(\tau ,z,t,y){|^{q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\beta )q}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle {\left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{|{v_{k}}(\tau ,z,t,y){|^{q}}{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\alpha -\gamma )q}}}\hspace{0.1667em}\mathrm{d}z\right)^{1/q}}.\end{array}\](22)
\[\begin{array}{r}\displaystyle |{v_{k}}(s,x,t,y)|\le {R_{k,T}}\times \\ {} \displaystyle \left(\frac{{(t-s)^{k\theta }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{k\theta +1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right),\end{array}\](23)
\[\begin{array}{r}\displaystyle {R_{k+1,T}}={R_{k,T}}{8^{1+1/q}}{N_{\beta ,T}}\| b{\| _{p}^{T}}{C^{1/q}}(1+{T^{\gamma /\alpha }})\times \\ {} \displaystyle \left(B{\left((k+1)\theta q,1\right)^{1/q}}\vee B{\left(1+k\theta q,\theta q\right)^{1/q}}\right).\end{array}\]
\[\begin{array}{r}\displaystyle |{v_{k+1}}(s,x,t,y)|\le 8{N_{\beta ,T}}\| b{\| _{p}^{T}}{R_{k,T}}\times \\ {} \displaystyle \left[{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\alpha -\gamma )q}}}\frac{{(t-\tau )^{k\theta q}}}{{({(t-\tau )^{1/\alpha }}+|z-y|)^{(d+\beta )q}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle \left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\alpha -\gamma )q}}}\frac{{(t-\tau )^{k\theta q+(1-\beta /\alpha )q}}}{{({(t-\tau )^{1/\alpha }}+|z-y|)^{(d+\alpha -\gamma )q}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle \left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{1}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\beta )q}}}\frac{{(t-\tau )^{k\theta q+(1-\beta /\alpha )q}}}{{({(t-\tau )^{1/\alpha }}+|z-y|)^{(d+\alpha -\gamma )q}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle {\left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{1}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\beta )q}}}\frac{{(t-\tau )^{k\theta q}}}{{({(t-\tau )^{1/\alpha }}+|z-y|)^{(d+\beta )q}}}\hspace{0.1667em}\mathrm{d}z\right]^{1/q}}.\end{array}\]
Using inequality (2), we can write
\[\begin{array}{r}\displaystyle |{v_{k+1}}(s,x,t,y)|\le 8{N_{\beta ,T}}\| b{\| _{p}^{T}}{R_{k,T}}{C^{1/q}}\times \\ {} \displaystyle \left[B\left((k+1)\theta q,1+\left(1-\frac{\beta }{\alpha }\right)q\right)\frac{{(t-s)^{(k+1)\theta q+(1-\beta /\alpha )q}}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\alpha -\gamma )q}}}+\right.\end{array}\]
\[\begin{array}{r}\displaystyle \left.B\left(1+k\theta q,\theta q+\frac{\gamma }{\alpha }q\right)\frac{{(t-s)^{(k+1)\theta q+q\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\beta )q}}}+\right.\\ {} \displaystyle \left.B\left((k+1)\theta q+\frac{\gamma }{\alpha }q,1+\left(1-\frac{\beta }{\alpha }\right)q\right)\frac{{(t-s)^{(k+1)\theta q+(1-\beta /\alpha +\gamma /\alpha )q}}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\alpha -\gamma )q}}}+\right.\end{array}\]
\[\begin{array}{r}\displaystyle \left.B\left(1+k\theta q+\left(1-\frac{\beta }{\alpha }\right),\theta q+\frac{\gamma }{\alpha }q\right)\frac{{(t-s)^{(k+1)\theta q+(1-\beta /\alpha +\gamma /\alpha )q}}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\alpha -\gamma )q}}}+\right.\\ {} \displaystyle \left.B\left((k+1)\theta q+\frac{\gamma }{\alpha }q,1\right)\frac{{(t-s)^{(k+1)\theta q+q\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\beta )q}}}+\right.\end{array}\]
\[\begin{array}{r}\displaystyle \left.B\left(1+k\theta q+\left(1-\frac{\beta }{\alpha }\right)q,\theta q\right)\frac{{(t-s)^{(k+1)\theta q+(1-\beta /\alpha )q}}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\alpha -\gamma )q}}}+\right.\\ {} \displaystyle \left.B\left((k+1)\theta q,1\right)\frac{{(t-s)^{(k+1)\theta q}}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\beta )q}}}+\right.\\ {} \displaystyle {\left.B\left(1+k\theta q,\theta q\right)\frac{{(t-s)^{(k+1)\theta q}}}{{({(t-s)^{1/\alpha }}+|x-y|)^{(d+\beta )q}}}\right]^{1/q}}.\end{array}\]
Since the function ${(B(x,y))_{x\gt 0,y\gt 0}}$ decreases with respect to each of its arguments, we have the inequality
\[\begin{array}{r}\displaystyle |{v_{k+1}}(s,x,t,y)|\le 8{N_{\beta ,T}}\| b{\| _{p}^{T}}{R_{k,T}}{C^{1/q}}\times \\ {} \displaystyle \left(B{\left((k+1)\theta q,1\right)^{1/q}}\vee B{\left(1+k\theta q,\theta q\right)^{1/q}}\right){8^{1/q}}(1+{T^{\gamma /\alpha }})\times \\ {} \displaystyle \left(\frac{{(t-s)^{(k+1)\theta }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{(k+1)\theta +1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right),\end{array}\]
that is, (22) holds for $k+1$. Therefore, it is correct for all $k=0,1,2,\dots \hspace{0.1667em}$, and relation (23) is true.Since ${\lim \nolimits_{k\to \infty }}B(1+k\theta q,\theta q)={\lim \nolimits_{k\to \infty }}B((k+1)\theta q,1)=0$, the series ${\textstyle\sum _{k=0}^{\infty }}{v_{k}}(s,x,t,y)$ converges uniformly with respect to $x,y\in {\mathbb{R}^{d}}$ and locally uniformly with respect to $0\le s\lt t$. Let ${(v(s,x,t,y))_{0\le s\lt t,x,y\in {\mathbb{R}^{d}}}}$ be the sum of this series. The function v is a solution to equation (14) and the estimate
\[ |v(s,x,t,y)|\le {C_{T}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right)\]
holds for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and every $T\gt 0$ with some constant ${C_{T}}\gt 0$ depended on T.The uniqueness of this solution follows from the fact that the difference ${v^{\ast }}$ of every two such solutions satisfies the equation
\[ {v^{\ast }}(s,x,t,y)={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,\tau ,z)(b(\tau ,z),{v^{\ast }}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z.\]
Therefore, since ${v^{\ast }}$ satisfies estimate (16), using inequality (2), one can obtain the inequality ($0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$)
\[\begin{array}{r}\displaystyle |{v^{\ast }}(s,x,t,y)|\le {R_{k,T}}\times \\ {} \displaystyle \left(\frac{{(t-s)^{k\theta }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{k\theta +1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right),\end{array}\]
for all $k\in \mathbb{N}$, $T\gt 0$, where ${R_{k,T}}$ is defined by (23). Since ${\lim \nolimits_{k\to \infty }}{R_{k,T}}{T^{k\theta }}=0$ for all $T\gt 0$, this means that ${v^{\ast }}(s,x,t,y)\equiv 0$.Let us define the function G by equality (13) with the function v just constructed. Using (2), (4) and Hölder’s inequality we can write the following chain of inequalities:
\[\begin{array}{r}\displaystyle |G(s,x,t,y)|\le {N_{0,T}}\frac{{(t-s)^{1-\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}+{2^{1+1/q}}{N_{0,T}}{C_{T}}\| b{\| _{p}^{T}}\times \\ {} \displaystyle \left[\left({\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\alpha -\gamma )q}}}\times \right.\right.\\ {} \displaystyle \left.{\left.\frac{1}{{({(t-\tau )^{1/\alpha }}+|z-y|)^{(d+\beta )q}}}\hspace{0.1667em}\mathrm{d}z\right)^{1/q}}+\right.\\ {} \displaystyle \left.\left({\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\alpha -\gamma )q}}}\times \right.\right.\\ {} \displaystyle \left.{\left.\frac{{(t-\tau )^{(1-\beta /\alpha )q}}}{{({(t-\tau )^{1/\alpha }}+|z-y|)^{(d+\beta )q}}}\hspace{0.1667em}\mathrm{d}z\right)^{1/q}}\right]\le \\ {} \displaystyle {C_{T}}\left[\frac{{(t-s)^{\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{1-\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right],\end{array}\]
where ${C_{T}}$ (in the last expression) is some positive constant, which might be depended on T.Let us prove that ${\nabla _{\beta }}G(s,\cdot ,t,y)(x)\equiv v(s,x,t,y)$. Using the statement of Lemma 3, we obtain the equality
\[\begin{array}{r}\displaystyle {\nabla _{\beta }}G(s,\cdot ,t,y)(x)={v_{0}}(s,x,t,y)+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,\tau ,z)(b(\tau ,z),v(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z=v(s,x,t,y),\end{array}\]
which is true for all $0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$.4 Family of evolution operators
Let us define the two-parameter family of operators $\{{𝕋_{st}}:0\le s\lt t\}$ in the space of continuous bounded functions ${C_{b}}({\mathbb{R}^{d}})$ defined by the equality
Similarly to the proof of Theorem 1 one can prove the following statement.
(24)
\[ {𝕋_{st}}\varphi (x)={\int _{{\mathbb{R}^{d}}}}G(s,x,t,y)\varphi (y)\hspace{0.1667em}\mathrm{d}y,\hspace{1em}\varphi \in {C_{b}}({\mathbb{R}^{d}}),\hspace{2.5pt}x\in {\mathbb{R}^{d}}.\]Lemma 4.
The function $w(s,x,t,\varphi )={\textstyle\int _{{\mathbb{R}^{d}}}}v(s,x,t,y)\varphi (y)\hspace{0.1667em}\mathrm{d}y$, $0\le s\lt t$, $x\in {\mathbb{R}^{d}}$ and $\varphi \in {C_{b}}({\mathbb{R}^{d}})$ (v is defined in Theorem 1), is a unique (in the class of functions, which satisfy the inequality $|w(s,x,t,\varphi )|\le {C_{T}}{(t-s)^{-\beta /\alpha }}$) solution to the equation
for all $0\le s\lt t\le T$, $x\in {\mathbb{R}^{d}}$ and every $T\gt 0$. Here ${v_{0}}$ is defined by (15) and ${w_{0}}(s,x,t,\varphi )={\textstyle\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,t,y)\varphi (y)\hspace{0.1667em}\mathrm{d}y$.
(25)
\[\begin{array}{r}\displaystyle w(s,x,t,\varphi )={w_{0}}(s,x,t,\varphi )+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z,\end{array}\]Proof.
Equation (25) is obtained from equation (14) multiplying it by the function φ and using the Fubini theorem. The justification of the usage of the Fubini theorem is based on estimates (5), (16) (see also (15)). Indeed,
\[\begin{array}{r}\displaystyle {\int _{{\mathbb{R}^{d}}}}\hspace{0.1667em}\mathrm{d}y{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}|{v_{0}}(s,x,\tau ,z)||b(\tau ,z)||v(\tau ,z,t,y)||\varphi (y)|\hspace{0.1667em}\mathrm{d}z\le \\ {} \displaystyle 4\| \varphi \| \| b{\| _{p}^{T}}{N_{\beta ,T}}{C_{T}}\times \\ {} \displaystyle {\int _{{\mathbb{R}^{d}}}}\hspace{0.1667em}\mathrm{d}y\left({\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{1}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{(d+\beta )q}}}\times \right.\\ {} \displaystyle \left.\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{(d+\beta )q}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle \left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{1}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{(d+\beta )q}}}\frac{{(t-\tau )^{(1-\beta /\alpha )q}}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{(d+\alpha -\gamma )q}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle \left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{(d+\alpha -\gamma )q}}}\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{q(d+\beta )}}}\hspace{0.1667em}\mathrm{d}z+\right.\\ {} \displaystyle \left.{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\beta /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{(d+\alpha -\gamma )q}}}\times \right.\\ {} \displaystyle {\left.\frac{{(t-\tau )^{(1-\beta /\alpha )q}}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{(d+\alpha -\gamma )q}}}\hspace{0.1667em}\mathrm{d}z\right)^{1/q}}\le \\ {} \displaystyle {C_{T}}{\int _{{\mathbb{R}^{d}}}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right)\hspace{0.1667em}\mathrm{d}y\le \\ {} \displaystyle {C_{T}}{(t-s)^{-\beta /\alpha }}\end{array}\]
for each $T\gt 0$ and every $0\le s\lt t\le T$, $x\in {\mathbb{R}^{d}}$ with some constants ${C_{T}}\gt 0$. Here we used the well-known formula2 ${\textstyle\int _{{\mathbb{R}^{d}}}}{(a+|x|)^{-d-\varkappa }}\hspace{0.1667em}\mathrm{d}x={a^{-\varkappa }}B(d,\varkappa )\frac{2{\pi ^{d/2}}}{\Gamma (d/2)}$, which is valid for all $a\gt 0$ and $\varkappa \gt 0$. □The next theorem contains the properties of the family of operators (24).
Theorem 2.
Let the assumptions of Theorem 1 hold. Then the following statements are true:
-
• the operators ${𝕋_{st}}$, $0\le s\lt t$, are linear and bounded on ${C_{b}}({\mathbb{R}^{d}})$;
-
• if $\varphi (x)\equiv 1$ then ${𝕋_{st}}\varphi (x)\equiv 1$;
-
• the family of operators $\{{𝕋_{st}}:0\le s\lt t\}$ has an evolution property, that is, ${𝕋_{s\tau }}{𝕋_{\tau t}}={𝕋_{st}}$ for all $0\le s\lt \tau \lt t$;
-
• ${\textit{w-}\text{lim}_{s\uparrow t}}{𝕋_{st}}=I$, where I is the identical operator, i.e. ${\lim \nolimits_{s\uparrow t}}{𝕋_{st}}\varphi (x)=\varphi (x)$, $x\in {\mathbb{R}^{d}}$ for all $\varphi \in {C_{b}}({\mathbb{R}^{d}})$.
Proof.
The linearity of operator ${𝕋_{st}}$ is evident. Let us prove its boundedness. If $\varphi \in {C_{b}}({\mathbb{R}^{d}})$, then using inequality (17) we can write ($\| \varphi \| ={\max _{x\in {\mathbb{R}^{d}}}}|\varphi (x)|$)
\[\begin{array}{r}\displaystyle |{𝕋_{st}}\varphi (x)|\le \| \varphi \| {\int _{{\mathbb{R}^{d}}}}|G(s,x,t,y)|\hspace{0.1667em}\mathrm{d}y\le \\ {} \displaystyle {C_{T}}\| \varphi \| {\int _{{\mathbb{R}^{d}}}}\left(\frac{{(t-s)^{\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta }}}+\frac{{(t-s)^{1-\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\alpha -\gamma }}}\right)\hspace{0.1667em}\mathrm{d}y\le \\ {} \displaystyle {C_{T}}\| \varphi \| \end{array}\]
for all $0\le s\lt t\le T$ and each $T\gt 0$. Therefore the operators ${𝕋_{st}}$ are bounded.Next, if $\varphi (x)\equiv 1$, then
\[\begin{array}{r}\displaystyle {𝕋_{st}}\varphi (x)={\int _{{\mathbb{R}^{d}}}}g(s,x,t,y)\hspace{0.1667em}\mathrm{d}y+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),{\int _{{\mathbb{R}^{d}}}}v(\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}y)\hspace{0.1667em}\mathrm{d}z.\end{array}\]
The function $w(s,x,t,1)={\textstyle\int _{{\mathbb{R}^{d}}}}v(\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}y$, $0\le s\lt t$, $x\in {\mathbb{R}^{d}}$, is a solution to equation (25) with ${w_{0}}(s,x,t,1)\equiv 0$. So, $w(s,x,t,1)\equiv 0$ and ${𝕋_{st}}\varphi (x)={\textstyle\int _{{\mathbb{R}^{d}}}}g(s,x,t,y)\hspace{0.1667em}\mathrm{d}y\equiv 1$.Although the evolution property can be proved in a standard way (see, for example, [11, 12, 14, 16, 17]), we will provide it here. For this, let us choose arbitrary $0\le s\lt u\lt t$, $\varphi \in {C_{b}}({\mathbb{R}^{d}})$, $x\in {\mathbb{R}^{d}}$ and consider
Using (25), one can obtain (by changing the order of integration)
\[\begin{array}{r}\displaystyle {𝕋_{st}}\varphi (x)={𝕋_{st}^{0}}\varphi (x)+{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z=\\ {} \displaystyle {𝕋_{su}^{0}}({𝕋_{ut}^{0}}\varphi )(x)+{\int _{s}^{u}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z+\\ {} \displaystyle {\int _{{\mathbb{R}^{d}}}}g(s,x,u,y)\hspace{0.1667em}\mathrm{d}y{\int _{u}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(u,y,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z=\\ {} \displaystyle {𝕋_{su}^{0}}({𝕋_{ut}}\varphi )(x)+{\int _{s}^{u}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z,\end{array}\]
where
(26)
\[ {𝕋_{st}^{0}}\varphi (x)={\int _{{\mathbb{R}^{d}}}}g(s,x,t,y)\varphi (y)\hspace{0.1667em}\mathrm{d}y.\]
\[ w(s,x,t,\varphi )={\int _{{\mathbb{R}^{d}}}}\varphi (y)\hspace{0.1667em}\mathrm{d}y{\int _{{\mathbb{R}^{d}}}}{\nabla _{\beta }}g(s,\cdot ,u,z)(x)g(u,z,t,y)\hspace{0.1667em}\mathrm{d}z+\]
\[\begin{array}{r}\displaystyle {\int _{s}^{u}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z+\\ {} \displaystyle {\int _{u}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,u,y)g(u,y,\tau ,z)\hspace{0.1667em}\mathrm{d}y(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z=\\ {} \displaystyle {w_{0}}(s,x,u,{𝕋_{ut}}\varphi )+{\int _{s}^{u}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{v_{0}}(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z.\end{array}\]
Here we used the well-known relation (the Chapman–Kolmogorov equation)
which is true for all $0\le s\lt u\lt t$, $x,y\in {\mathbb{R}^{d}}$.Since equation (25) has a unique solution, we state that
for all $\varphi \in {C_{b}}({\mathbb{R}^{d}})$, $x\in {\mathbb{R}^{d}}$, $0\le s\lt u\lt t$. Therefore ${𝕋_{st}}\varphi (x)={𝕋_{su}}({𝕋_{ut}}\varphi )(x)$ and the evolution property is proved.
The last statement of this theorem follows from these two facts: first, using (26) we have the equality ${\textit{w}\text{-lim}_{s\uparrow t}}{𝕋_{st}^{0}}=I$, and second (note that $1\le q\lt \frac{d+\alpha }{d+1}$),
\[\begin{array}{r}\displaystyle \left|{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z\right|\le \\ {} \displaystyle \| b{\| _{p}^{T}}{N_{0,T}}{C_{T}}{\left({\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{{(\tau -s)^{(1-\gamma /\alpha )q}}}{{({(\tau -s)^{1/\alpha }}+|x-z|)^{(d+\alpha -\gamma )q}}}{(t-\tau )^{-q\beta /\alpha }}\hspace{0.1667em}\mathrm{d}z\right)^{1/q}}\le \\ {} \displaystyle {C_{T}}{\left({\int _{s}^{t}}{(\tau -s)^{-(q-1)d/\alpha }}{(t-\tau )^{-q\beta /\alpha }}\hspace{0.1667em}\mathrm{d}\tau \right)^{1/q}}=\\ {} \displaystyle {C_{T}}{(t-s)^{1-q\beta /\alpha -(q-1)d/\alpha }}\to 0,\hspace{2.5pt}\text{as}\hspace{2.5pt}s\uparrow t.\end{array}\]
This completes the proof of the theorem. □Remark 3.
We cannot state that the operators ${𝕋_{st}}$ preserve the cone of nonnegative functions. We have no proof that the function ${𝕋_{st}}\varphi (x)$ can have negative values if $\varphi (x)\ge 0$, $x\in {\mathbb{R}^{d}}$. But the example of an α-stable process and $b(t,x)\equiv b\in {\mathbb{R}^{d}}$ confirms this fact. Exactly analogously to how it is done in [1] for the case $\beta =\alpha -1$, we can obtain the equality
in our case. The function G here has negative values as the Fourier transform of a nonpositive definite function.
Thus, the family of evolution operators $\{{𝕋_{st}}:0\le s\lt t\}$ does not define any of the Markov processes but only a pseudo-process possessing a “Markov property”.
5 Cauchy problem
In this section we fix some $T\gt 0$ and prove that the function $u(s,x,t)={𝕋_{st}}\varphi (x)$, $0\le s\lt t\le T$, $x\in {\mathbb{R}^{d}}$, is a solution (in some sense) to the following Cauchy problem:
for every $t\in (0,T]$, where $L(s,x)=A(s,x)+(b(s,x),{\nabla _{\beta }})$ and $\varphi \in {C_{b}}({\mathbb{R}^{d}})$.
Let a sequence $\{{b_{n}}:n\in \mathbb{N}\}\subset {L_{p}}([0,T]\times {\mathbb{R}^{d}})$ be such that ${L_{p}}\text{-}{\lim \nolimits_{n\to \infty }}{b_{n}}=b$ and the functions ${b_{n}}$ satisfy the assumptions of Theorem 1. Let ${u_{n}}(s,x,t)={𝕋_{st}^{(n)}}\varphi (x)$, where the operators ${𝕋_{st}^{(n)}}$ are constructed using Theorems 1, 2 and functions ${b_{n}}$ instead of b. If the functions ${u_{n}}$ satisfy the Cauchy problem (27), (28) with the function ${b_{n}}$ instead of b and the pointwise limit $u={\lim \nolimits_{n\to \infty }}{u_{n}}$ exists, we will call the function u by the generalized solution to the Cauchy problem (27), (28).
The following auxiliary statement will be useful in constructing this generalized solution.
Lemma 5.
Let functions ${(\tilde{b}(s,x))_{s\ge 0,x\in {\mathbb{R}^{d}}}}$ and ${(\hat{b}(s,x))_{s\ge 0,x\in {\mathbb{R}^{d}}}}$ satisfy the assumptions of Theorem 1 with the same p. Then the corresponding functions $\tilde{G}$ and $\hat{G}$ satisfy the inequality
for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and each $T\gt 0$ with some positive constant ${C_{T}}$.
(29)
\[\begin{array}{r}\displaystyle |\tilde{G}(s,x,t,y)-\hat{G}(s,x,t,y)|\le {C_{T}}\| \tilde{b}-\hat{b}{\| _{p}^{T}}(1+\| \tilde{b}{\| _{p}^{T}}+\| \hat{b}{\| _{p}^{T}})\times \\ {} \displaystyle \frac{1}{{({(t-s)^{1/\alpha }}+|x-y|)^{d+\beta -\gamma }}},\end{array}\]Proof.
Using (13), we can obtain the equality
for all $0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$, where $W(s,x,t,y)=(\hat{b}(s,x),\hat{v}(s,x,t,y))-(\tilde{b}(s,x),\tilde{v}(s,x,t,y))$, in which $\hat{v}$ and $\tilde{v}$ are solutions to equations obtained from (14) replacing the function b by the functions $\hat{b}$ and $\tilde{b}$, respectively. Relation (14) leads us to the equation
(30)
\[ \hat{G}(s,x,t,y)-\tilde{G}(s,x,t,y)={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}g(s,x,\tau ,z)W(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\]
\[\begin{array}{r}\displaystyle W(s,x,t,y)={W_{0}}(s,x,t,y)+{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}({v_{0}}(s,x,\tau ,z),\hat{b}(s,x))W(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z+\\ {} \displaystyle {\int _{{\mathbb{R}^{d}}}}{W_{0}}(s,x,\tau ,z)(\tilde{b}(\tau ,z),\tilde{v}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z,\end{array}\]
where ${W_{0}}(s,x,t,y)=(\hat{b}(s,x)-\tilde{b}(s,x),{v_{0}}(s,x,t,y))$. Considering (5), we obtain the inequality
\[\begin{array}{r}\displaystyle |{W_{0}}(s,x,t,y)\le |\hat{b}(s,x)-\tilde{b}(s,x)|{N_{\beta ,T}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}\right.+\\ {} \displaystyle \left.\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right),\end{array}\]
valid for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and each $T\gt 0$. Moreover,
\[\begin{array}{r}\displaystyle \left|{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{W_{0}}(s,x,\tau ,z)(\tilde{b}(\tau ,z),\tilde{v}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\right|\le \\ {} \displaystyle {N_{\beta ,T}}|\hat{b}(s,x)-\tilde{b}(s,x)|{C_{T}}{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}|\tilde{b}(\tau ,z)|\times \\ {} \displaystyle \left(\frac{1}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\beta }}}+\frac{{(\tau -s)^{1-\beta /\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha -\gamma }}}\right)\times \\ {} \displaystyle \left(\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}+\frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\right)\hspace{0.1667em}\mathrm{d}z=\\ {} \displaystyle {C_{T}}|\hat{b}(s,x)-\tilde{b}(s,x)|{\sum \limits_{k=1}^{4}}{I_{k}},\end{array}\]
where
\[ {I_{1}}={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}|\tilde{b}(\tau ,z)|\frac{1}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\beta }}}\frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}\hspace{0.1667em}\mathrm{d}z,\]
\[\begin{array}{r}\displaystyle {I_{2}}={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}|\tilde{b}(\tau ,z)|\frac{1}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\beta }}}\times \\ {} \displaystyle \frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\hspace{0.1667em}\mathrm{d}z,\end{array}\]
\[\begin{array}{r}\displaystyle {I_{3}}={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}|\tilde{b}(\tau ,z)|\frac{{(\tau -s)^{1-\beta /\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha -\gamma }}}\times \\ {} \displaystyle \frac{1}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\beta }}}\hspace{0.1667em}\mathrm{d}z,\end{array}\]
\[\begin{array}{r}\displaystyle {I_{4}}={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}|\tilde{b}(\tau ,z)|\frac{{(\tau -s)^{1-\beta /\alpha }}}{{({(\tau -s)^{1/\alpha }}+|z-x|)^{d+\alpha -\gamma }}}\times \\ {} \displaystyle \frac{{(t-\tau )^{1-\beta /\alpha }}}{{({(t-\tau )^{1/\alpha }}+|y-z|)^{d+\alpha -\gamma }}}\hspace{0.1667em}\mathrm{d}z.\end{array}\]
Since $\tilde{b}\in {L_{p}}([0,T]\times {\mathbb{R}^{d}})$ for each $T\gt 0$, using inequality (2) we obtain the estimates
\[ {I_{1}}\le \| \tilde{b}{\| _{p}^{T}}{(2CB(1,\theta q))^{1/q}}\frac{{(t-s)^{\theta }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}},\]
\[\begin{array}{r}\displaystyle {I_{2}}\le \| \tilde{b}{\| _{p}^{T}}{\left(2C\left(B\left(1,\theta q+q\frac{\gamma }{\alpha }\right)\vee B\left(1+\left(1-\frac{\beta }{\alpha }\right)q,\theta q\right)\right)\right)^{1/q}}\times \\ {} \displaystyle \left(\frac{{(t-s)^{\theta +\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}+\frac{{(t-s)^{\theta +1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right),\end{array}\]
\[\begin{array}{r}\displaystyle {I_{3}}\le \| \tilde{b}{\| _{p}^{T}}{\left(2C\left(B\left(1,\theta q+q\frac{\gamma }{\alpha }\right)\vee B\left(1+\left(1-\frac{\beta }{\alpha }\right)q,\theta q\right)\right)\right)^{1/q}}\times \\ {} \displaystyle \left(\frac{{(t-s)^{\theta +\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}+\frac{{(t-s)^{\theta +1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right),\end{array}\]
\[ {I_{4}}\le \| \tilde{b}{\| _{p}^{T}}{\left(2CB\left(1+\left(1-\frac{\beta }{\alpha }\right)q,\theta q+q\frac{\gamma }{\alpha }\right)\right)^{1/q}}\frac{{(t-s)^{\theta +1-\beta /\alpha +\gamma /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}},\]
where, as it was above, $\theta =1-((d+\alpha )/p+\beta )/\alpha $, $q=p/(p-1)$ and C is maximum of positive constants derived from inequality (2) in considered four cases. Therefore,
\[\begin{array}{r}\displaystyle \left|{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{W_{0}}(s,x,\tau ,z)(\tilde{b}(\tau ,z),\tilde{v}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z\right|\le {C_{T}}\| \tilde{b}{\| _{p}^{T}}|\hat{b}(s,x)-\tilde{b}(s,x)|\times \\ {} \displaystyle \left(\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}+\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right).\end{array}\]
Let us denote
This equation can be solved by the method of successive approximations, i.e. a solution of it can be found in the form $W(s,x,t,y)={\textstyle\sum _{k=0}^{\infty }}{W_{k}^{\ast }}(s,x,t,y)$. The terms of this series satisfy the relation ($k=1,2,\dots \hspace{0.1667em}$)
valid for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and each $T\gt 0$, where ${C_{T}}$ is some positive constant depended on T.
\[ {W_{0}^{\ast }}(s,x,t,y)={W_{0}}(s,x,t,y)+{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{W_{0}}(s,x,\tau ,z)(\tilde{b}(\tau ,z),\tilde{v}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z.\]
Then we can write the equality
(31)
\[\begin{array}{r}\displaystyle W(s,x,t,y)={W_{0}^{\ast }}(s,x,t,y)+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}({v_{0}}(s,x,\tau ,z),\hat{b}(s,x))W(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z.\end{array}\]
\[ {W_{k}^{\ast }}(s,x,t,y)={\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}({v_{0}}(s,x,\tau ,z),\hat{b}(s,x)){W_{k-1}^{\ast }}(\tau ,z,t,y))\hspace{0.1667em}\mathrm{d}z.\]
To justify this, note that
\[\begin{array}{r}\displaystyle |{W_{0}^{\ast }}(s,x,t,y)|\le {C_{T}}(1+\| \tilde{b}{\| _{p}^{T}})|\hat{b}(s,x)-\tilde{b}(s,x)|\times \\ {} \displaystyle \left(\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}+\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right).\end{array}\]
Moreover, by the method of mathematical induction one can prove the estimate
\[\begin{array}{r}\displaystyle |{W_{k}^{\ast }}(s,x,t,y)|\le {R_{k}}|\hat{b}(s,x)|\| \hat{b}-\tilde{b}{\| _{p}^{T}}\times \\ {} \displaystyle \left(\frac{{(t-s)^{k\theta }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}+\frac{{(t-s)^{k\theta +1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right)\end{array}\]
valid for all $k\ge 1$, where ${R_{1}}={C_{T}}(1+\| \tilde{b}{\| _{p}^{T}})$ and ${C_{T}}$ is some positive constant depended on T, and, for $k\ge 2$,
\[\begin{array}{r}\displaystyle {R_{k}}={R_{k-1}}{C_{T}}(1+{T^{\gamma /\alpha }})(1+\| \tilde{b}{\| _{p}^{T}})\| \hat{b}{\| _{p}^{T}}\times \\ {} \displaystyle {(8C(B(1,k\theta q)\vee B(\theta q,(k-1)\theta q+1)))^{1/q}}.\end{array}\]
Therefore, the series ${\textstyle\sum _{k=0}^{\infty }}{W_{k}^{\ast }}(s,x,t,y)$ converges uniformly with respect to $x,y\in {\mathbb{R}^{d}}$ and locally uniformly with respect to $0\le s\lt t$. So, its sum W is a solution to equation (31). In addition, we obtain the estimate
(32)
\[\begin{array}{r}\displaystyle |W(s,x,t,y)|\le {C_{T}}\left((1+\| \tilde{b}{\| _{p}^{T}})|\hat{b}(s,x)-\tilde{b}(s,x)|+|\hat{b}(s,x)|\| \hat{b}-\tilde{b}{\| _{p}^{T}}\right)\times \\ {} \displaystyle \left(\frac{1}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\beta }}}+\frac{{(t-s)^{1-\beta /\alpha }}}{{({(t-s)^{1/\alpha }}+|y-x|)^{d+\alpha -\gamma }}}\right),\end{array}\]Lemma 6.
Let the function w be defined in Lemma 4. For every $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$, the inequality
holds with some positive constant ${C_{T}}$ depended on T.
Proof.
Using (25) we can write down the relation ($0\le s\lt t$, $x,y\in {\mathbb{R}^{d}}$)
Let us remark that (see (6))
(33)
\[\begin{array}{r}\displaystyle w(s,x,t,\varphi )-w(s,y,t,\varphi )={w_{0}}(s,x,t,\varphi )-{w_{0}}(s,y,t,\varphi )+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}({v_{0}}(s,x,\tau ,z)-{v_{0}}(s,y,\tau ,z))(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z.\end{array}\]
\[ {v_{0}}(s,x,t,y)={\nabla _{\beta }}{g_{0}}(s,\cdot ,t,y)(x)+{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}{\nabla _{\beta }}{g_{0}}(s,\cdot ,\tau ,z)(x)\Phi (\tau ,z,t,y)\hspace{0.1667em}\mathrm{d}z.\]
Moreover, for all $0\le s\lt t$, $x,y,z\in {\mathbb{R}^{d}}$,
\[\begin{array}{r}\displaystyle {\nabla _{\beta }}{g_{0}}(s,\cdot ,t,z)(x)-{\nabla _{\beta }}{g_{0}}(s,\cdot ,t,z)(y)=\\ {} \displaystyle \frac{i}{{(2\pi )^{d}}}{\int _{{\mathbb{R}^{d}}}}\left({e^{i(x-z,\lambda )}}-{e^{i(y-z,\lambda )}}\right)\lambda |\lambda {|^{\beta -1}}\exp \left\{-a(t,z,\lambda )(t-s)\right\}\hspace{0.1667em}\mathrm{d}\lambda =\\ {} \displaystyle -\frac{1}{{(2\pi )^{d}}}{\int _{{\mathbb{R}^{d}}}}{e^{i(\eta -z,\lambda )}}(\lambda ,x-y)\lambda |\lambda {|^{\beta -1}}\exp \left\{-a(t,z,\lambda )(t-s)\right\}\hspace{0.1667em}\mathrm{d}\lambda =\\ {} \displaystyle |x-y|{D_{\beta +1}}{g_{0}}(s,\cdot ,t,z)(\eta ),\end{array}\]
where $\eta =\theta x+(1-\theta )y$ with some $\theta \in (0,1)$ and the operator ${D_{\beta +1}}$ is defined by the symbol ${\left(\left(\lambda ,\frac{x-y}{|x-y|}\right)\lambda |\lambda {|^{\beta -1}}\right)_{\lambda \in {\mathbb{R}^{d}}}}$. This symbol satisfies the assumptions of [6, Lemma 4.2], which leads us to the estimate
\[ |{D_{\beta +1}}{g_{0}}(s,\cdot ,t,z)(\eta )|\le \frac{C}{{({(t-s)^{1/\alpha }}+|z-\eta |)^{d+\beta +1}}}.\]
Thus, using (2), we obtain the inequality
\[\begin{array}{r}\displaystyle |{v_{0}}(s,x,t,z)-{v_{0}}(s,y,t,z)|\le \frac{C|x-y|}{{({(t-s)^{1/\alpha }}+|z-\eta |)^{d+\beta +1}}}+\\ {} \displaystyle {\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\frac{C|x-y|}{{({(\tau -s)^{1/\alpha }}+|u-\eta |)^{d+\beta +1}}}\frac{C}{{({(t-\tau )^{1/\alpha }}+|z-u|)^{d+\alpha -\gamma }}}\hspace{0.1667em}\mathrm{d}u\le \\ {} \displaystyle |x-y|{C_{T}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|z-\eta |)^{d+\beta +1}}}+\frac{{(t-s)^{1-(\beta +1)/\alpha }}}{{({(t-s)^{1/\alpha }}+|z-\eta |)^{d+\alpha -\gamma }}}\right),\end{array}\]
valid for all $0\le s\lt t\le T$, $x,y,z\in {\mathbb{R}^{d}}$. And, as consequence, the following inequalities are true for all $0\le s\lt t\le T$, $x,y\in {\mathbb{R}^{d}}$ and $\varphi \in {C_{b}}({\mathbb{R}^{d}})$:
\[ |{v_{0}}(s,x,t,\varphi )-{v_{0}}(s,y,t,\varphi )|\le {C_{T}}|x-y|\| \varphi \| {(t-s)^{-(\beta +1)/\alpha }};\]
\[\begin{array}{r}\displaystyle \left|{\int _{s}^{t}}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}({v_{0}}(s,x,\tau ,z)-{v_{0}}(s,y,\tau ,z))(b(\tau ,z),w(\tau ,z,t,\varphi ))\hspace{0.1667em}\mathrm{d}z\right|\le \\ {} \displaystyle {C_{T}}|x-y|\| b{\| _{p}^{T}}\left({\int _{s}^{t}}{(t-\tau )^{-\beta q/\alpha }}\mathrm{d}\tau {\int _{{\mathbb{R}^{d}}}}\left(\frac{1}{{({(t-s)^{1/\alpha }}+|z-\eta |)^{d+\beta +1}}}+\right.\right.\\ {} \displaystyle {\left.{\left.\frac{{(t-s)^{1-(\beta +1)/\alpha }}}{{({(t-s)^{1/\alpha }}+|z-\eta |)^{d+\alpha -\gamma }}}\right)^{q}}\hspace{0.1667em}\mathrm{d}z\right)^{1/q}}\le \\ {} \displaystyle {C_{T}}|x-y|{(t-s)^{1-(1+d/\alpha )/p-(2\beta +1)/\alpha }},\end{array}\]
where $q=p/(p-1)$.Therefore, using (33) and the fact that $1-(1+d/\alpha )/p-\beta /\alpha \gt (1-\beta )/\alpha $, we obtain the statement of the lemma. □
In the next theorem, we construct a generalized solution to the Cauchy problem formulated at the beginning of the section.
Theorem 3.
Let the assumptions of Theorem 1 hold and the function G be constructed there. Then the function
\[ u(s,x,t)={\int _{{\mathbb{R}^{d}}}}G(s,x,t,y)\varphi (y)\hspace{0.1667em}\mathrm{d}y,\hspace{1em}0\le s\lt t\le T,\hspace{2.5pt}x\in {\mathbb{R}^{d}},\]
is a generalized solution to the Cauchy problem (27), (28) for each $\varphi \in {C_{b}}({\mathbb{R}^{d}})$.
Proof.
Let us consider a sequence $\left\{{b_{n}}:n\in \mathbb{N}\right\}$ of ${\mathbb{R}^{d}}$-valued functions, which are infinitely continuously differentiable, have compact supports and belong to ${L_{p}}([0,T]\times {\mathbb{R}^{d}})$, where p is defined in Theorem 1. Assume that ${L_{p}}\text{-}{\lim \nolimits_{n\to \infty }}{b_{n}}=b$, where b is the function from Theorem 1.
We denote by ${v_{n}}$, ${w_{n}}$ and ${G_{n}}$ the objects that are defined as v, w and G, respectively, using ${b_{n}}$ instead of b.
Lemma 5 allows us to state that the sequence of corresponding functions ${G_{n}}$ constructed in Theorem 1 using the functions ${b_{n}}$ instead of b converges to the function G uniformly with respect to $y\in {\mathbb{R}^{d}}$ for each fixed $0\le s\lt t\le T$ and $x\in {\mathbb{R}^{d}}$.
Let us consider the function
\[ {f_{n}}(s,x,t)={\int _{{\mathbb{R}^{d}}}}({b_{n}}(s,x),{\nabla _{\beta }}{G_{n}}(s,\cdot ,t,y)(x))\varphi (y)\hspace{0.1667em}\mathrm{d}y,\hspace{1em}0\le s\lt t\le T,\hspace{2.5pt}x\in {\mathbb{R}^{d}}.\]
Remind that ${\nabla _{\beta }}{G_{n}}(s,\cdot ,t,y)(x)={v_{n}}(s,x,t,y)$. Moreover,
\[\begin{array}{r}\displaystyle |{f_{n}}(s,x,t)-{f_{n}}(s,y,t)|\le |{b_{n}}(s,x)||{w_{n}}(s,x,t,\varphi )-{w_{n}}(s,y,t,\varphi )|+\\ {} \displaystyle |{b_{n}}(s,x)-{b_{n}}(s,y)||{w_{n}}(s,y,t,\varphi )|.\end{array}\]
Since the function ${b_{n}}$ is Lipschitz continuous with respect to x, uniformly to s and has a compact support, then, as follows from Lemma 6, the function ${f_{n}}(s,x,t)$ is Lipschitz continuous with respect to x, uniformly to $s\in [0,t)$ for every fixed $t\gt 0$.