Linear integro-differential operator and Nonlocal minimal surfaces: Difference between pages

From nonlocal pde
(Difference between pages)
Jump to navigation Jump to search
imported>Luis
No edit summary
 
imported>Nestor
No edit summary
 
Line 1: Line 1:
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
In broad and vague terms, these surfaces arise as the boundaries of domains $E \subset \mathbb{R}^n$ that minimize (within a class of given admissible configurations) the energy functional:


\[ 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) \]
\[ J_s(E)= C_{n,s}\int_{E}\int_{E^c}\frac{1}{|x-y|^{n+s}}dxdy,\;\; s \in (0,1) \]
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. Many theorems, and in particular regularity theorems, require extra assumptions in the kernels $K$. These assumptions restrict the study to certain sub-classes of linear operators. The simplest of all is the [[fractional Laplacian]]. We list below several extra assumptions that are usually made.
It can be checked easily that this agrees (save for a factor of $2$) with  norm of the characteristic function $\chi_E$ in the homogenous Sobolev space  $\dot{H}^{\frac{s}{2}}$. The dimensional constant $C_{n,s}$ blows up as $s \to 1^-$, in which case (at least when the boundary of $E$ is smooth enough) one can check that $J_s(E)$ converges to the perimeter of $E$.  


== Absolutely continuous measure ==
Classically,  [[minimal surfaces]] (or generally [[surfaces of constant mean curvature]] ) arise in physical situations where one has two phases interacting (eg. water-air, water-ice ) and the energy of interaction is proportional to the area of the interface, which is due to the interaction between particles/agents in both phases being negligible when they are far apart.


In most cases, the nonnegative measure $\mu$ is assumed to be absolutely continuous: $\mathrm{d} \mu_x(y) = K(x,y) \mathrm{d}y$.
Nonlocal minimal surfaces then, describe physical phenomena where the interaction potential does not decay fast enough as particles are apart, so that two particles on different phases and far from the interface still contribute a non-trivial amount to the total interaction energy, in particular, one may consider much more general energy functionals corresponding to different interaction potentials


We keep this assumption in all the examples below.
\[ J_K(E)= \int_{E}\int_{E^c}K(x,y) dxdy \]


== Purely integro-differential operator ==
== Definition ==  


In this case we neglect the local part of the operator
Following the most accepted convention for [[minimal surfaces]],  a nonlocal minimal surface is the boundary $\Sigma$ of an open set $E \subset \mathbb{R}^n$ such that $\chi_E \in \dot{H}^{s/2}$ and whose [[Nonlocal mean curvature]] $H_s$ is identically zero, that is
\[ 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 ==
H_s(x) = \[ \int_{\mathbb{R}^n} \frac{\chi_E(y)-\chi_{E^c}(y)}{|x-y|^{n+s}}dy=0 \;\;\forall\; x \in \Sigma\]
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
In this case we say that $\Sigma$ is a nonlocal minimal surface in $\Omega$.
\[ 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))$.
Note that for this definition to make sense, $\Sigma$ must be the boundary of some open set $E$, in this article, we will often refer to the set $E$ itself as the minimal surface, and no confusion should arise from this.  


== Translation invariant operators ==
Example: Suppose that $E$ and $\Omega$ is such that for any other domain $F$ such that $F \Delta E \subset \subset \Omega$ (i.e. $F$ agrees with $E$ outside $\Omega$) we have
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 ==
$J_s(E) \leq J_s(F)$


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.
Then, if it is the case that $E$ has a smooth enough boundary, one can check that the nonlocal mean curvature of $\Sigma$ is identically zero in $\Omega$.


\[ -(-\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. \]
== Nonlocal mean curvature ==


== Stable operators ==
== Surfaces minimizing non-local energy functionals ==


These are the operators whose kernel is homogeneous in $y$
== The Caffarelli-Roquejoffre-Savin Regularity Theorem==
\[ K(x,y)=\frac{a(x,y/|y|)}{|y|^{n+s}}.\]
They are the generators of stable Lévy processes.
 
== 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
\[ |\partial_y^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.
 
== More singular/irregular kernels ==
 
The concept of uniform ellipticity can be relaxed in various ways. The following family of kernels was considered in the paper of Silvestre and Schwab.<ref name="schwab2014regularity" />
 
== 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}}. \]
This family of kernels has the outstanding property that the corresponding linear operator $L$ coincides with $Lu(x) = a_{ij} \partial_{ij}\left[(-\Delta)^{(s-2)/2} u \right] (x)$ for some coefficients $a_{ij}$.
 
== Second order elliptic operators as limits of purely integro-differential ones ==
 
Given any bounded, even, positive function $a: \mathbb{R}^n\to \mathbb{R}$, the family of operators
 
\[L_\sigma u(x) = (2-\sigma) \int_{\mathbb{R}^n} (u(x+y)+u(x-y)-2u(x))\frac{a(y)}{|y|^{n+\sigma}}dy,\;\; \sigma \in (0,2), \]
 
define in the limit $\sigma \to 2^-$ a second order linear elliptic operator (possibly degenerate). This can be checked for any fixed $C^2$ function $u$ by a straightforward computation using the second order Taylor expansion. A class of kernels that is big enough to recover all translation invariant elliptic operators of the form  $Lu(x) = Tr ( A \cdot D^2u(x) )$ is given by the kernels
 
\[ K_A(y) = (2-\sigma) \frac{1}{|Ay|^{n+\sigma}},\]
 
where $A$ is an invertible symmetric matrix.
 
== Characterization via global maximum principle ==
 
A bounded linear operator
 
\[ L: C^2_0(\mathbb{R}^n) \to C(\mathbb{R}^n) \]
 
is said to satisfy the global maximum principle if given any $u \in C^2_0(\mathbb{R}^n)$ with a global maximum at some point $x_0$ we have
 
\[ (Lu)(x_0) \leq 0 \]
 
It turns out this property imposes strong restrictions on the operator $L$, and we have the following theorem due to Courrège <ref name="C65"/> <ref name="C64"/>: if $L$ satisfies the global maximum principle then it has the form
 
\[ Lu(x) = \mathrm{tr} \, A(x) \cdot D^2 u + b(x) \cdot \nabla u + c(x) u + \int_{\R^n} (u(x+y) - u(x) - y \cdot \nabla u(x) \chi_{B_1}(y)) \, \mathrm{d} \mu_x(y) \]
 
where again $A(x)$ is a nonnegative matrix for all $x$, $c(x)\leq 0$ and $\mu_x$ is a nonnegative measure for all $x$ satisfying
 
\[ \int_{\R^n} \min(y^2 , 1) \mathrm{d} \mu_x(y) < +\infty. \]
 
and $A(x),c(x)$ and $b(x)$ are bounded.
 
 
== See also ==
 
* [[Fractional Laplacian]]
* [[Levy processes]]
* [[Dirichlet form]]
 
 
== References ==
{{reflist|refs=
<ref name="C64">{{Citation | last1=Courrège | first1=Philippe | title=Générateur infinitésimal d'un semi-groupe de convolution sur $R^n$, et formule de Lévy-Khinchine | year=1964 | journal=Bulletin des Sciences Mathématiques. 2e Série | issn=0007-4497 | volume=88 | pages=3–30}}</ref>
<ref name="C65">{{Citation | last1=Courrège | first1=P. | title=Sur la forme intégro-différentielle des opéateurs de  $C_k^\infty(\mathbb{R}^n)$  dans $C(\mathbb{R}^n)$ satisfaisant au principe du maximum | journal=Sém. Théorie du potentiel (1965/66) Exposé | volume=2}}</ref>
<ref name="schwab2014regularity">{{Citation | last1=Schwab | first1= Russell W | last2=Silvestre | first2= Luis | title=Regularity for parabolic integro-differential equations with very irregular kernels | journal=arXiv preprint arXiv:1412.3790}}</ref>
}}

Revision as of 13:10, 31 May 2011

In broad and vague terms, these surfaces arise as the boundaries of domains $E \subset \mathbb{R}^n$ that minimize (within a class of given admissible configurations) the energy functional:

\[ J_s(E)= C_{n,s}\int_{E}\int_{E^c}\frac{1}{|x-y|^{n+s}}dxdy,\;\; s \in (0,1) \]

It can be checked easily that this agrees (save for a factor of $2$) with norm of the characteristic function $\chi_E$ in the homogenous Sobolev space $\dot{H}^{\frac{s}{2}}$. The dimensional constant $C_{n,s}$ blows up as $s \to 1^-$, in which case (at least when the boundary of $E$ is smooth enough) one can check that $J_s(E)$ converges to the perimeter of $E$.

Classically, minimal surfaces (or generally surfaces of constant mean curvature ) arise in physical situations where one has two phases interacting (eg. water-air, water-ice ) and the energy of interaction is proportional to the area of the interface, which is due to the interaction between particles/agents in both phases being negligible when they are far apart.

Nonlocal minimal surfaces then, describe physical phenomena where the interaction potential does not decay fast enough as particles are apart, so that two particles on different phases and far from the interface still contribute a non-trivial amount to the total interaction energy, in particular, one may consider much more general energy functionals corresponding to different interaction potentials

\[ J_K(E)= \int_{E}\int_{E^c}K(x,y) dxdy \]

Definition

Following the most accepted convention for minimal surfaces, a nonlocal minimal surface is the boundary $\Sigma$ of an open set $E \subset \mathbb{R}^n$ such that $\chi_E \in \dot{H}^{s/2}$ and whose Nonlocal mean curvature $H_s$ is identically zero, that is

H_s(x) = \[ \int_{\mathbb{R}^n} \frac{\chi_E(y)-\chi_{E^c}(y)}{|x-y|^{n+s}}dy=0 \;\;\forall\; x \in \Sigma\]

In this case we say that $\Sigma$ is a nonlocal minimal surface in $\Omega$.

Note that for this definition to make sense, $\Sigma$ must be the boundary of some open set $E$, in this article, we will often refer to the set $E$ itself as the minimal surface, and no confusion should arise from this.

Example: Suppose that $E$ and $\Omega$ is such that for any other domain $F$ such that $F \Delta E \subset \subset \Omega$ (i.e. $F$ agrees with $E$ outside $\Omega$) we have

$J_s(E) \leq J_s(F)$

Then, if it is the case that $E$ has a smooth enough boundary, one can check that the nonlocal mean curvature of $\Sigma$ is identically zero in $\Omega$.

Nonlocal mean curvature

Surfaces minimizing non-local energy functionals

The Caffarelli-Roquejoffre-Savin Regularity Theorem