The Fokker-Planck equation, also known as Kolmogorov forward equation, is a hyperbolic partial differential equation, the evolution equation of the probability density of a certain class of stochastic differential equations.
has an associated Fokker-Planck equation (as physicists say) aka Kolmogorov forward equation (as mathematicians say), that describes the time evolution of the probability distribution of the Langevin equation: From the physicist’s POV this probability distribution tells us what the probability is to find the particle at a certain time in a certain place.
(The Smoluchowski Equation is a special kind of Fokker-Planck equation.)
For a one dimensional Ito stochastic differential equation of the form
the probability density is a solution of this one-dimensional Fokker-Planck equation:
The mathematical theorem making all of this precise is the Feynman-Kac formula.
The differential operator on the right side is sometimes called the Fokker-Planck operator.
While the Fokker-Planck equation describes the evolution of the probability distribution forwards in time, the backward Fokker-Planck equation, or Kolmogorov backward equation, describes the evolution of the probability distribution backwards in time. The equation is an evolution equation with the formal adjoint of the Fokker-Planck operator:
The Fokker-Planck equation describes the time evolution of a probability distribution, and therefore satisfies a conservation law that describes that no probability is lost, the probability of the trivial event is still 1. For the sake of simplicity we state it in one spatial dimension only, we have:
with the probability current
A stationary solution is a solution that does not depend on time and is therefore an eigenvector of the Fokker-Planck operator to the eigenvalue . In one dimension without explicit time dependency of the Fokker-Planck operator, it can be calculated in closed form:
is a normalization constant to ensure that
This formula is obviously valid only if the function can be normalized in this way, of course.
A particle in one dimension in a potential that is described by a Langevin equation of the form
has the following Fokker-Planck equation:
Therefore, the stationary solution is according to our general formula, this:
This is true for all potential functions that are differentiable and such that .
In our toy example for stochastic resonance, the stationary solution to the time independent system with the bistable potential
is, according to the preceding paragraph, simply
The following picture shows both the potential and the stationary solution (not normalized though).
For applications it is often of interest, when a system or a particle leaves a certain area and crosses a threshold. This kind of problem is often called a First-hitting-time problem, and these can often be handled using techniques involving Fokker-Planck equations.
Here we will consider an example discussed here:
The noise matrix of the SDE system is
In the multidimensional case we have to calculate the diffusion matrix which is
The Fokker-Planck equation for the probability distribution is therefore
Note that the Fokker-Planck equation has no explicit dependency on the angle : We expected this, since the original system is rotationally invariant.
Some handwaving tells us that there should be a stationary equillibirium solution, and that this should be rotationally invariant, too. If we assume that actually depends on only, the Fokker-Planck equation simplifies to
If we further impose the conditions that the function value and its first derivative should vanish in , we get as solution
The normalization constant N has to be determined numerically.
We’ll consider two examples, and and . For the first case we get , for the second .
(Green is and red is .)
The following book is a classic reference for the practicioner, Risken explains model building and approximate solution methods like linear response theory? using a minimum of mathematical background:
The mathematically inclined reader will find a nice exposition in the following classic textbook, chapter 4 “Brownian Motion and Partial Differential Equations”: