Fully nonlinear integro-differential equations and Linear integro-differential operator: Difference between pages

From nonlocal pde
(Difference between pages)
Jump to navigation Jump to search
imported>Luis
No edit summary
 
imported>Luis
 
Line 1: Line 1:
Fully nonlinear integro-differential equations are a nonlocal version of fully nonlinear elliptic equations of the form $F(D^2 u, Du, u, x)=0$. The main examples are the integro-differential [[Bellman equation]] from optimal control, and the [[Isaacs equation]] from stochastic games.
The linear integro-differential operators that we consider ''in this wiki'' are the generators of [[Levy processes]]. According to the Levy-Kintchine formula, they have the general form


The general definition of ellipticity provided below does not require a specific form of the equation. However, the main two applications are the two above.
\[ Lu(x) = \mathrm{tr} \, A(x) \cdot D^2 u + b(x) \cdot \nabla u + c(x) u + d(x) + \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, \mathrm{d} \mu_x(y) \]
where $A(x)$ is a nonnegative matrix for all $x$, and $\mu_x$ is a nonnegative measure for all $x$ satisfying
\[ \int_{\R^n} \min(y^2 , 1) \mathrm{d} \mu_x(y) < +\infty. \]


== Definition <ref name="CS"/><ref name="CS2"/> ==
The above definition is very general. In most cases we are interested in some subclass of linear operators. The simplest of all is the [[fractional Laplacian]]. We list below several extra assumptions that are usually made.


Given a family of [[linear integro-differential operators]] $\mathcal{L}$, we define the [[extremal operators]] $M^+_\mathcal{L}$ and $M^-_\mathcal{L}$:
== Absolutely continuous measure ==
\begin{align*}
M^+_\mathcal{L} u(x) &= \sup_{L \in \mathcal{L}} \, L u(x) \\
M^-_\mathcal{L} u(x) &= \inf_{L \in \mathcal{L}} \, L u(x)
\end{align*}


We define a nonlinear operator $I$ to be elliptic in a domain $\Omega$ with respect to the class $\mathcal{L}$ if it assigns a continuous function $Iu$ to every function $u \in L^\infty(\R^n) \cap C^2(\Omega)$, and moreover for any two such functions $u$ and $v$:
In most cases, the nonnegative measure $\mu$ is assumed to be absolutely continuous: $\mathrm{d} \mu_x(y) = K(x,y) \mathrm{d}y$.
\[M^-_\mathcal{L} [u-v](x)\leq Iu(x) - Iv(x) \leq M^+_\mathcal{L} [u-v] (x), \]
for any $x \in \Omega$.


A fully nonlinear elliptic equation with respect to $\mathcal{L}$ is an equation of the form $Iu=0$ in $\Omega$ with $I$ uniformly elliptic respect to $\mathcal{L}$.
We keep this assumption in all the examples below.


<div style="background:#DDEEFF;">
== Purely integro-differential operator ==
<blockquote>
'''Note'''. If $\mathcal L$ consists of purely second order operators of the form $\mathrm{tr} \, A \cdot D^2 u$ with $\lambda I \leq A \leq \Lambda I$, then $M^+_{\mathcal L}$ and $M^-_{\mathcal L}$ denote the extremal Pucci operators. It is ''folklore'' statement that then nonlinear operator $I$ elliptic respect to $\mathcal L$ in the sense described above must coincide with a fully nonlinear elliptic operator of the form $Iu = F(D^2u,x)$. However, this proof may have never been written anywhere.
</blockquote>
</div>


== Examples ==
In this case we neglect the local part of the operator
\[ Lu(x) = \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, K(x,y) \mathrm d y. \]


The two main examples are the following.
== Symmetric kernels ==
If the kernel is symmetric $K(x,y) = K(x,-y)$, then we can remove the gradient term from the integral and replace the difference by a second order quotient.


* The [[Bellman equation]] is the equality
In the purely integro-differentiable case, it reads as
\[ \sup_{a \in \mathcal{A}} \, L_a u(x) = f(x), \]
\[ Lu(x) = \frac 12 \int_{\R^n} (u(x+y)+u(x+y)-2u(x)) \, K(x,y) \mathrm d y. \]
where $L_a$ is some family of linear integro-differential operators indexed by an arbitrary set $\mathcal{A}$.


The equation appears naturally in problems of stochastic control with [[Levy processes]].
The second order incremental quotient is sometimes abbreviated by $\delta u(x,y) := (u(x+y)+u(x+y)-2u(x))$.


The equation is uniformly elliptic with respect to any class $\mathcal{L}$ that contains all the operators $L_a$.
== Translation invariant operators ==
In this case, all coefficients are independent of $x$.
\[ Lu(x) = \mathrm{tr} \, A \cdot D^2 u + b \cdot \nabla u + c u + d + \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, K(y) \mathrm{d}y. \]


* The [[Isaacs equation]] is the equality
== The fractional Laplacian ==
\[ \sup_{a \in \mathcal{A}} \ \inf_{b \in \mathcal{B}} \ L_{ab} u(x) = f(x), \]
where $L_{ab}$ is some family of linear integro-differential operators with two indices $a \in \mathcal A$ and $b \in \mathcal B$.


The equation appears naturally in zero sum stochastic games with [[Levy processes]].
The [[fractional laplacian]] is the simplest and most common purely integro-differential operator. It corresponds to a translation invariant operator for which $K(y)$ is radially symmetric and homogeneous.


The equation is uniformly elliptic with respect to any class $\mathcal{L}$ that contains all the operators $L_{ab}$.
\[ -(-\Delta)^{s/2} u(x) = C_{n,s} \int_{\R^n} (u(x+y)+u(x+y)-2u(x)) \frac{1}{|y|^{n+s}} \mathrm d y. \]


<div style="background:#DDEEFF;">
== Uniformly elliptic of order $s$ ==
<blockquote>
'''Note'''. In several articles <ref name="BI"/><ref name="BIC2"/><ref name="BIC"/>, fully nonlinear integro-differential equations of the form $F(D^2 u, Du, u, x, Lu)=f(x)$ are analyzed, where $L$ is a [[linear integro-differential operator]]. This is a rigid structure for an equation because if, for example, an equation is purely integro-differential (it does not depend on $D^2u$, $Du$ or $u$) then it is forced to be linear: $Lu(x) = [F(x,\cdot)^{-1}f(x)]$.


On the other hand, the results in these papers apply to the more general definitions of fully nonlinear integro-differential equations as well. The reason for that restriction seems to be just to have an equation that is short to write down.
This corresponds to the assumption that the kernel is comparable to the one of the fractional Laplacian of the same order.
</blockquote>
\[ \frac {(2-s)\lambda}{|y|^{n+s}} \leq K(x,y) \leq \frac {(2-s)\Lambda}{|y|^{n+s}}. \]
</div>


== References ==
The normalizing factor $(2-s)$ is a normalizing factor which is only important when $s$ approaches two.
{{reflist|refs=
 
<ref name="CS2">{{Citation | last1=Caffarelli | first1=Luis | last2=Silvestre | first2=Luis | title=Regularity results for nonlocal equations by approximation | publisher=[[Springer-Verlag]] | location=Berlin, New York | journal=Archive for Rational Mechanics and Analysis | issn=0003-9527 | pages=1–30}}</ref>
An operator of variable order can be either one for which $s$ depends on $x$, or one for which there are two values $s_1<s_2$, one for the left hand side and another for the right hand side.
<ref name="CS">{{Citation | last1=Caffarelli | first1=Luis | last2=Silvestre | first2=Luis | title=Regularity theory for fully nonlinear integro-differential equations | url=http://dx.doi.org/10.1002/cpa.20274 | doi=10.1002/cpa.20274 | year=2009 | journal=[[Communications on Pure and Applied Mathematics]] | issn=0010-3640 | volume=62 | issue=5 | pages=597–638}}</ref>
 
<ref name="BIC">
== Smoothness class $k$ of order $s$ ==
{{Citation | last1=Barles | first1=Guy | last2=Chasseigne | first2=Emmanuel | last3=Imbert | first3=Cyril | title=Hölder continuity of solutions of second-order non-linear elliptic integro-differential equations | url=http://dx.doi.org/10.4171/JEMS/242 | doi=10.4171/JEMS/242 | year=2011 | journal=Journal of the European Mathematical Society (JEMS) | issn=1435-9855 | volume=13 | issue=1 | pages=1–26}}</ref>
This class (sometimes denoted as $\mathcal L_k^s$) corresponds to kernels that are uniformly elliptic of order $s$ and, moreover, their derivatives are also bounded
<ref name="BIC2">{{Citation | last1=Barles | first1=G. | last2=Chasseigne | first2=Emmanuel | last3=Imbert | first3=Cyril | title=On the Dirichlet problem for second-order elliptic integro-differential equations | url=http://dx.doi.org/10.1512/iumj.2008.57.3315 | doi=10.1512/iumj.2008.57.3315 | year=2008 | journal=Indiana University Mathematics Journal | issn=0022-2518 | volume=57 | issue=1 | pages=213–246}}</ref>
\[ D^r K(x,y) \leq \frac {\Lambda}{|y|^{n+s+r}} \ \ \text{for all } r\leq k. \]
<ref name="BI">{{Citation | last1=Barles | first1=Guy | last2=Imbert | first2=Cyril | title=Second-order elliptic integro-differential equations: viscosity solutions' theory revisited | url=http://dx.doi.org/10.1016/j.anihpc.2007.02.007 | doi=10.1016/j.anihpc.2007.02.007 | year=2008 | journal=Annales de l'Institut Henri Poincaré. Analyse Non Linéaire | issn=0294-1449 | volume=25 | issue=3 | pages=567–585}}</ref>
 
}}
== Order strictly below one ==
 
If a non symmetric kernel $K$ satisfies the extra local integrability assumption
\[ \int_{\R^n} \min(|y|,1) K(x,y) \mathrm d y < +\infty, \]
then the extra gradient term is not necessary in order to define the operator.
 
\[ Lu(x) = \int_{\R^n} (u(x+y) - u(x)) \, K(x,y) \mathrm d y. \]
 
The modification in the integro-differential part of the operator becomes an extra drift term.
 
A uniformly elliptic operator of order $s<1$ satisfies this condition.
 
== Order strictly above one ==
 
If a non symmetric kernel $K$ satisfies the extra integrability assumption on its tail.
\[ \int_{\R^n} \min(|y|^2,|y|) K(x,y) \mathrm d y < +\infty, \]
then the gradient term in the integral can be taken global instead of being cut off in the unit ball.
 
\[ Lu(x) = \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x)) \, K(x,y) \mathrm d y. \]
 
The modification in the integro-differential part of the operator becomes an extra drift term.
 
A uniformly elliptic operator of order $s>1$ satisfies this condition.
 
== Indexed by a matrix ==
In some cases, it is interesting to study a family of kernels $K$ that are indexed by a matrix. For example, given the matrix $A$, one can consider the kernel of order $s$:
\[ K_A(y) = \frac{(2-s) \langle y , Ay \rangle}{|y|^{n+2+s}}. \]

Revision as of 14:39, 27 May 2011

The linear integro-differential operators that we consider in this wiki are the generators of Levy processes. According to the Levy-Kintchine formula, they have the general form

\[ Lu(x) = \mathrm{tr} \, A(x) \cdot D^2 u + b(x) \cdot \nabla u + c(x) u + d(x) + \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, \mathrm{d} \mu_x(y) \] where $A(x)$ is a nonnegative matrix for all $x$, and $\mu_x$ is a nonnegative measure for all $x$ satisfying \[ \int_{\R^n} \min(y^2 , 1) \mathrm{d} \mu_x(y) < +\infty. \]

The above definition is very general. In most cases we are interested in some subclass of linear operators. The simplest of all is the fractional Laplacian. We list below several extra assumptions that are usually made.

Absolutely continuous measure

In most cases, the nonnegative measure $\mu$ is assumed to be absolutely continuous: $\mathrm{d} \mu_x(y) = K(x,y) \mathrm{d}y$.

We keep this assumption in all the examples below.

Purely integro-differential operator

In this case we neglect the local part of the operator \[ Lu(x) = \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, K(x,y) \mathrm d y. \]

Symmetric kernels

If the kernel is symmetric $K(x,y) = K(x,-y)$, then we can remove the gradient term from the integral and replace the difference by a second order quotient.

In the purely integro-differentiable case, it reads as \[ Lu(x) = \frac 12 \int_{\R^n} (u(x+y)+u(x+y)-2u(x)) \, K(x,y) \mathrm d y. \]

The second order incremental quotient is sometimes abbreviated by $\delta u(x,y) := (u(x+y)+u(x+y)-2u(x))$.

Translation invariant operators

In this case, all coefficients are independent of $x$. \[ Lu(x) = \mathrm{tr} \, A \cdot D^2 u + b \cdot \nabla u + c u + d + \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, K(y) \mathrm{d}y. \]

The fractional Laplacian

The fractional laplacian is the simplest and most common purely integro-differential operator. It corresponds to a translation invariant operator for which $K(y)$ is radially symmetric and homogeneous.

\[ -(-\Delta)^{s/2} u(x) = C_{n,s} \int_{\R^n} (u(x+y)+u(x+y)-2u(x)) \frac{1}{|y|^{n+s}} \mathrm d y. \]

Uniformly elliptic of order $s$

This corresponds to the assumption that the kernel is comparable to the one of the fractional Laplacian of the same order. \[ \frac {(2-s)\lambda}{|y|^{n+s}} \leq K(x,y) \leq \frac {(2-s)\Lambda}{|y|^{n+s}}. \]

The normalizing factor $(2-s)$ is a normalizing factor which is only important when $s$ approaches two.

An operator of variable order can be either one for which $s$ depends on $x$, or one for which there are two values $s_1<s_2$, one for the left hand side and another for the right hand side.

Smoothness class $k$ of order $s$

This class (sometimes denoted as $\mathcal L_k^s$) corresponds to kernels that are uniformly elliptic of order $s$ and, moreover, their derivatives are also bounded \[ D^r K(x,y) \leq \frac {\Lambda}{|y|^{n+s+r}} \ \ \text{for all } r\leq k. \]

Order strictly below one

If a non symmetric kernel $K$ satisfies the extra local integrability assumption \[ \int_{\R^n} \min(|y|,1) K(x,y) \mathrm d y < +\infty, \] then the extra gradient term is not necessary in order to define the operator.

\[ Lu(x) = \int_{\R^n} (u(x+y) - u(x)) \, K(x,y) \mathrm d y. \]

The modification in the integro-differential part of the operator becomes an extra drift term.

A uniformly elliptic operator of order $s<1$ satisfies this condition.

Order strictly above one

If a non symmetric kernel $K$ satisfies the extra integrability assumption on its tail. \[ \int_{\R^n} \min(|y|^2,|y|) K(x,y) \mathrm d y < +\infty, \] then the gradient term in the integral can be taken global instead of being cut off in the unit ball.

\[ Lu(x) = \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x)) \, K(x,y) \mathrm d y. \]

The modification in the integro-differential part of the operator becomes an extra drift term.

A uniformly elliptic operator of order $s>1$ satisfies this condition.

Indexed by a matrix

In some cases, it is interesting to study a family of kernels $K$ that are indexed by a matrix. For example, given the matrix $A$, one can consider the kernel of order $s$: \[ K_A(y) = \frac{(2-s) \langle y , Ay \rangle}{|y|^{n+2+s}}. \]