Dirichlet form and Levy processes: Difference between pages

From nonlocal pde
(Difference between pages)
Jump to navigation Jump to search
imported>Nestor
No edit summary
 
imported>Nestor
No edit summary
 
Line 1: Line 1:
$$
A Lévy process is an important type of [[stochastic process]] (namely, a family of $\mathbb{R}^d$ valued random variables each indexed by a positive number $t\geq 0$).  In the context of parabolic integro-differential equations they play the same role that Brownian motion and more general diffusions play in the theory of second order parabolic equations.
\newcommand{\dd}{\mathrm{d}}
\newcommand{\R}{\mathbb{R}}
$$


A Dirichlet form in $\mathbb{R}^n$ is a bilinear function
Informally speaking, a Lévy process is a Brownian motion which may jump, the times, length and direction of the jumps being random variables. A prototypical example would be $X(t)=B(t)+N(t)$ where $B(t)$ is the standard [[Brownian motion]] and $N(t)$ is a [[Compound Poisson process]], the trajectory described by typical sample path of this process would look like the union of several disconnected Brownian motion paths.


\begin{equation*}
== Definition ==
\mathcal{E}: D\times D \to \mathbb{R}
\end{equation*}


with the following properties
A stochastic process $X=\{X(t)\}_{t \geq 0}$ with values in $\mathbb{R}^d$ is said to be a Lévy process if


1) The domain $D$ is a dense subset of $L^2(\mathbb{R}^n)$
1.For any sequence $0 \leq t_1 < t_2 <...<t_n$ the random variables $X(t_0),X(t_1)-X(t_0),...,X(t_n)-X(t_{n-1})$ are independent.


2) $\mathcal{E}$ is symmetric, that is $\mathcal{E}(u,v)=\mathcal{E}(v,u)$ for any $u,v \in D$.
2.For any positive times $s\leq t$ the random variables $X(t-s)$ and $X(t)-X(s)$ have the same probability law.


3) $\mathcal{E}(u,u) \geq 0$ for any $u \in D$.
3.Almost surely, the trajectory of $X(t)$ is continuous from the right, with limit from the left also known as "càdlàg" for its acronym in french.  


4) The set $D$ equipped with the inner product defined by $(u,v)_{\mathcal{E}} := (u,v)_{L^2(\mathbb{R}^n)} + \mathcal{E}(u,v)$ is a real Hilbert space.
== Lévy-Khintchine Formula ==


5) For any $u \in D$ we have that $u_* = (u\vee 0) \wedge 1 \in D$ and $\mathcal{E}(u_*,u_*)\leq \mathcal{E}(u,u)$
It follows from the first two properties above that if $X$ is a Lévy process and we further assume $X(0)=0$ a.s. then for each fixed positive $t$ the random variable $X(t)$ is infinitely divisible, that is, it can be written as the sum of $n$ independent and identically distribued random variables, for all $n\in\mathbb{N}$. Indeed, let $h=\tfrac{t}{n}$, then


\[X(t) = \left( X(h)-X(0)\right)+\left( X(2h)-X(h)\right)+...+\left( X(t)-X((n-1)h)\right)\]


In other words, a Dirichlet form is nothing but a positive symmetric bilinear form defined in a dense subset of $L^2(\mathbb{R}^n)$ such that 4) and 5) hold. Alternatively, the quadratic form $u \to \mathcal{E}(u,u)$ itself is known as the Dirichlet form and it is still denoted by $\mathcal{E}$, so $\mathcal{E}(u):=\mathcal{E}(u,u)$.
and by the above definition the differences $X(kh)-X((k-1)h)$ are independent and distributed the same as $X(h)$. From the infinite divisibility of $X(t)$ it follows by a theorem of Lévy and Khintchine that  for any $\xi \in \mathbb{R}^d$ we have


The best known Dirichlet form is the Dirichlet energy
\[ \mathbb{E} \left [ e^{i\xi\cdot X_t}\right ] = e^{t\eta(\xi)}\]
\begin{equation*}
\mathcal{E}(u) = \int_{\mathbb{R}^n} |\nabla u|^2\;dx
\end{equation*}


which gives rise to the space $H^1(\mathbb{R}^n)$.  Another example of a Dirichlet form is given by  
the function $\eta(\xi)$ given by
\begin{equation*}
\mathcal{E}(u) = \iint_{\R^n \times \R^n} (u(y)-u(x))^2 k(x,y)\, \dd x \dd y
\end{equation*}
where $K$ is some non-negative symmetric kernel.


If the kernel $K$ satisfies the bound $K(x,y) \leq \Lambda |x-y|^{-n-s}$, then the quadratic form is bounded in $\dot H^{s/2}$ . If moreover, $\lambda |x-y|^{-n-s} \leq K(x,y)$, then the form is comparable to the norm in $\dot H^{s/2}$ squared and in that case the set $D \subset L^2(\mathbb{R}^n)$ defined above is given by  $H^{s/2}(\mathbb{R}^n)$
\[\eta(\xi)=y\cdot b i  -\tfrac{1}{2}(A\xi,\xi)+\int_{\mathbb{R}^d} \left ( e^{\xi\cdot y}-1-i\xi\cdot y \chi_{B_1}(y) \right ) d\mu(y) \]


Then we see that Dirichlet forms are natural generalizations of the Dirichlet integrals
where $b$ is a vector, $A$ is a positive matrix and $\mu$ is a Lévy measure, that is, a Borel measure in $\mathbb{R}^d$ such that
\[ \int a_{ij}(x) \partial_i u \partial_j u \dd x, \]
where $a_{ij}(x)$ is a positive matrix.


The Euler-Lagrange equation of a Dirichlet form is a non-local analogue  of an elliptic equations in divergence form. Equations of this type are studied using variational methods and they are expected to satisfy similar properties <ref name="BBCK"/><ref name="K"/><ref name="CCV"/>.
\[ \int_{\mathbb{R}^d}\frac{|y|^2}{1+|y|^2}d\mu(y) <+\infty \]


== References ==
== Connection to linear integro-differential operators ==
(There should be a lot more references here)
{{reflist|refs=
<ref name="CCV">{{Citation | last1=Caffarelli | first1=Luis | last2=Chan | first2=Chi Hin | last3=Vasseur | first3=Alexis | title= | doi=10.1090/S0894-0347-2011-00698-X | year=2011 | journal=[[Journal of the American Mathematical Society]] | issn=0894-0347 | issue=24 | pages=849–869}}</ref>
<ref name="BBCK">{{Citation | last1=Barlow | first1=Martin T. | last2=Bass | first2=Richard F. | last3=Chen | first3=Zhen-Qing | last4=Kassmann | first4=Moritz | title=Non-local Dirichlet forms and symmetric jump processes | url=http://dx.doi.org/10.1090/S0002-9947-08-04544-3 | doi=10.1090/S0002-9947-08-04544-3 | year=2009 | journal=[[Transactions of the American Mathematical Society]] | issn=0002-9947 | volume=361 | issue=4 | pages=1963–1999}}</ref>
<ref name="K">{{Citation | last1=Kassmann | first1=Moritz | title=A priori estimates for integro-differential operators with measurable kernels | url=http://dx.doi.org/10.1007/s00526-008-0173-6 | doi=10.1007/s00526-008-0173-6 | year=2009 | journal=Calculus of Variations and Partial Differential Equations | issn=0944-2669 | volume=34 | issue=1 | pages=1–21}}</ref>
}}


Any Lévy process $X(t)$ such that $X(0)=0$ almost surely defines a semigroup $\{U_t\}_{t\geq0}$ on the space  of continuous functions $f:\mathbb{R}^d\to\mathbb{R}^d$ as follows


{{stub}}
\[(U_tf)(x)= \mathbb{E}\left [ f(x+X(t) \right ]\]
 
Given the initial assumption on $X(0)$ it is clear that $U_0$ is the identity, and given that $X(t)-X(s)$ is distributed as $X(t-s)$ it follows that $U_t \circ U_s = U_{t+s}$.
 
As a semigroup, $U_t$ has an infinitesimal generator which turns out to be a [[Linear integro-differential operator]]. More precisely, if we let $f(x,t):=(U_tf)(x)$, then, assuming that $f(x,t)$ has enough regularity it can be checked that
 
\[\partial_t f = Lf \;\;\;\mbox{ for all } (x,t)\in\mathbb{R}^d\times \mathbb{R}_+\]
 
where for any smooth function $\phi$, we have
 
\[ L\phi(x) =  b \cdot \nabla \phi(x) +\mathrm{tr} \,( A\cdot D^2 \phi )+ \int_{\R^d} (\phi(x+y) - \phi(x) - y \cdot \nabla \phi(x) \chi_{B_1}(y)) \, \mathrm{d} \mu(y) \]

Revision as of 14:36, 22 January 2012

A Lévy process is an important type of stochastic process (namely, a family of $\mathbb{R}^d$ valued random variables each indexed by a positive number $t\geq 0$). In the context of parabolic integro-differential equations they play the same role that Brownian motion and more general diffusions play in the theory of second order parabolic equations.

Informally speaking, a Lévy process is a Brownian motion which may jump, the times, length and direction of the jumps being random variables. A prototypical example would be $X(t)=B(t)+N(t)$ where $B(t)$ is the standard Brownian motion and $N(t)$ is a Compound Poisson process, the trajectory described by typical sample path of this process would look like the union of several disconnected Brownian motion paths.

Definition

A stochastic process $X=\{X(t)\}_{t \geq 0}$ with values in $\mathbb{R}^d$ is said to be a Lévy process if

1.For any sequence $0 \leq t_1 < t_2 <...<t_n$ the random variables $X(t_0),X(t_1)-X(t_0),...,X(t_n)-X(t_{n-1})$ are independent.

2.For any positive times $s\leq t$ the random variables $X(t-s)$ and $X(t)-X(s)$ have the same probability law.

3.Almost surely, the trajectory of $X(t)$ is continuous from the right, with limit from the left also known as "càdlàg" for its acronym in french.

Lévy-Khintchine Formula

It follows from the first two properties above that if $X$ is a Lévy process and we further assume $X(0)=0$ a.s. then for each fixed positive $t$ the random variable $X(t)$ is infinitely divisible, that is, it can be written as the sum of $n$ independent and identically distribued random variables, for all $n\in\mathbb{N}$. Indeed, let $h=\tfrac{t}{n}$, then

\[X(t) = \left( X(h)-X(0)\right)+\left( X(2h)-X(h)\right)+...+\left( X(t)-X((n-1)h)\right)\]

and by the above definition the differences $X(kh)-X((k-1)h)$ are independent and distributed the same as $X(h)$. From the infinite divisibility of $X(t)$ it follows by a theorem of Lévy and Khintchine that for any $\xi \in \mathbb{R}^d$ we have

\[ \mathbb{E} \left [ e^{i\xi\cdot X_t}\right ] = e^{t\eta(\xi)}\]

the function $\eta(\xi)$ given by

\[\eta(\xi)=y\cdot b i -\tfrac{1}{2}(A\xi,\xi)+\int_{\mathbb{R}^d} \left ( e^{\xi\cdot y}-1-i\xi\cdot y \chi_{B_1}(y) \right ) d\mu(y) \]

where $b$ is a vector, $A$ is a positive matrix and $\mu$ is a Lévy measure, that is, a Borel measure in $\mathbb{R}^d$ such that

\[ \int_{\mathbb{R}^d}\frac{|y|^2}{1+|y|^2}d\mu(y) <+\infty \]

Connection to linear integro-differential operators

Any Lévy process $X(t)$ such that $X(0)=0$ almost surely defines a semigroup $\{U_t\}_{t\geq0}$ on the space of continuous functions $f:\mathbb{R}^d\to\mathbb{R}^d$ as follows

\[(U_tf)(x)= \mathbb{E}\left [ f(x+X(t) \right ]\]

Given the initial assumption on $X(0)$ it is clear that $U_0$ is the identity, and given that $X(t)-X(s)$ is distributed as $X(t-s)$ it follows that $U_t \circ U_s = U_{t+s}$.

As a semigroup, $U_t$ has an infinitesimal generator which turns out to be a Linear integro-differential operator. More precisely, if we let $f(x,t):=(U_tf)(x)$, then, assuming that $f(x,t)$ has enough regularity it can be checked that

\[\partial_t f = Lf \;\;\;\mbox{ for all } (x,t)\in\mathbb{R}^d\times \mathbb{R}_+\]

where for any smooth function $\phi$, we have

\[ L\phi(x) = b \cdot \nabla \phi(x) +\mathrm{tr} \,( A\cdot D^2 \phi )+ \int_{\R^d} (\phi(x+y) - \phi(x) - y \cdot \nabla \phi(x) \chi_{B_1}(y)) \, \mathrm{d} \mu(y) \]