Diffusion process

For the marketing term, see Diffusion of innovations.

In probability theory and statistics, diffusion processes are a class of continuous-time Markov process with almost surely continuous sample paths. Diffusion process is stochastic in nature and hence is used to model many real-life stochastic systems. Brownian motion, reflected Brownian motion and Ornstein–Uhlenbeck processes are examples of diffusion processes. It is used heavily in statistical physics, statistical analysis, information theory, data science, neural networks, finance and marketing.

A sample path of a diffusion process models the trajectory of a particle embedded in a flowing fluid and subjected to random displacements due to collisions with other particles, which is called Brownian motion. The position of the particle is then random; its probability density function as a function of space and time is governed by a convection–diffusion equation.

Mathematical definition

A diffusion process is a Markov process with continuous sample paths for which the Kolmogorov forward equation is the Fokker–Planck equation.^[1]

A diffusion process is defined by the following properties. Let ${\displaystyle a^{ij}(x,t)}$ be uniformly continuous coefficients and ${\displaystyle b^{i}(x,t)}$ be bounded, Borel measurable drift terms. There is a unique family of probability measures ${\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}$ (for ${\displaystyle \tau \geq 0}$, ${\displaystyle \xi \in \mathbb {R} ^{d}}$) on the canonical space ${\displaystyle \Omega =C([0,\infty ),\mathbb {R} ^{d})}$, with its Borel ${\displaystyle \sigma }$-algebra, such that:

1. (Initial Condition) The process starts at ${\displaystyle \xi }$ at time ${\displaystyle \tau }$: ${\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }[\psi \in \Omega$ :\psi (t)=\xi {\text{ for }}0\leq t\leq \tau ]=1.}

2. (Local Martingale Property) For every ${\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times [\tau ,\infty ))}$, the process ${\displaystyle M_{t}^{[f]}=f(\psi (t),t)-f(\psi (\tau ),\tau )-\int _{\tau }^{t}{\bigl (}L_{a;b}+{\tfrac {\partial }{\partial s}}{\bigr )}f(\psi (s),s)\,ds}$ is a local martingale under ${\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}$ for ${\displaystyle t\geq \tau }$, with ${\displaystyle M_{t}^{[f]}=0}$ for ${\displaystyle t\leq \tau }$.

This family ${\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}$ is called the ${\displaystyle {\mathcal {L}}_{a;b}}$-diffusion.

SDE Construction and Infinitesimal Generator

It is clear that if we have an ${\displaystyle {\mathcal {L}}_{a;b}}$-diffusion, i.e. ${\displaystyle (X_{t})_{t\geq 0}}$ on ${\displaystyle (\Omega ,{\mathcal {F}},{\mathcal {F}}_{t},\mathbb {P} _{a;b}^{\xi ,\tau })}$, then ${\displaystyle X_{t}}$ satisfies the SDE ${\displaystyle dX_{t}^{i}={\frac {1}{2}}\,\sum _{k=1}^{d}\sigma _{k}^{i}(X_{t})\,dB_{t}^{k}+b^{i}(X_{t})\,dt}$. In contrast, one can construct this diffusion from that SDE if ${\displaystyle a^{ij}(x,t)=\sum _{k}\sigma _{i}^{k}(x,t)\,\sigma _{j}^{k}(x,t)}$ and ${\displaystyle \sigma ^{ij}(x,t)}$, ${\displaystyle b^{i}(x,t)}$ are Lipschitz continuous. To see this, let ${\displaystyle X_{t}}$ solve the SDE starting at ${\displaystyle X_{\tau }=\xi }$. For ${\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times [\tau ,\infty ))}$, apply Itô's formula: ${\displaystyle df(X_{t},t)={\bigl (}{\frac {\partial f}{\partial t}}+\sum _{i=1}^{d}b^{i}{\frac {\partial f}{\partial x_{i}}}+v\sum _{i,j=1}^{d}a^{ij}\,{\frac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}{\bigr )}\,dt+\sum _{i,k=1}^{d}{\frac {\partial f}{\partial x_{i}}}\,\sigma _{k}^{i}\,dB_{t}^{k}.}$ Rearranging gives ${\displaystyle f(X_{t},t)-f(X_{\tau },\tau )-\int _{\tau }^{t}{\bigl (}{\frac {\partial f}{\partial s}}+L_{a;b}f{\bigr )}\,ds=\int _{\tau }^{t}\sum _{i,k=1}^{d}{\frac {\partial f}{\partial x_{i}}}\,\sigma _{k}^{i}\,dB_{s}^{k},}$ whose right‐hand side is a local martingale, matching the local‐martingale property in the diffusion definition. The law of ${\displaystyle X_{t}}$ defines ${\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}$ on ${\displaystyle \Omega =C([0,\infty ),\mathbb {R} ^{d})}$ with the correct initial condition and local martingale property. Uniqueness follows from the Lipschitz continuity of ${\displaystyle \sigma \!,\!b}$. In fact, ${\displaystyle L_{a;b}+{\tfrac {\partial }{\partial s}}}$ coincides with the infinitesimal generator ${\displaystyle {\mathcal {A}}}$ of this process. If ${\displaystyle X_{t}}$ solves the SDE, then for ${\displaystyle f(\mathbf {x} ,t)\in C^{2}(\mathbb {R} ^{d}\times \mathbb {R} ^{+})}$, the generator ${\displaystyle {\mathcal {A}}}$ is ${\displaystyle {\mathcal {A}}f(\mathbf {x} ,t)=\sum _{i=1}^{d}b_{i}(\mathbf {x} ,t)\,{\frac {\partial f}{\partial x_{i}}}+v\sum _{i,j=1}^{d}a_{ij}(\mathbf {x} ,t)\,{\frac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}+{\frac {\partial f}{\partial t}}.}$

See also

• Stochastic differential equation
• Itô calculus
• Fokker–Planck equation
• Markov process
• Diffusion
• Itô diffusion
• Jump diffusion
• Sample-continuous process