In studying the diffusion equation we have neglected the possible presence of any external force. We could ask, however, what happens if we include them; as we will see here, we will obtain the so called Fokker-Planck equation.
So, let's go back to the beginning and start with the same assumptions we made in The diffusion equation, with the difference that this time we also include an external force . Therefore if is the position of the particle at time , after a little time interval it will be:
is the so called mobility
of the particles.
Depending on the nature of the system, the expression of the mobility can vary. For example, if our particles are diffusing through a dense fluid, the mobility will be essentially related to the viscosity of the fluid itself; in fact, it is a known fact that the equation of motion of a body in a fluid is (in one dimension):
is the velocity of the body,
the viscosity of the fluid and
a geometric factor (for example, from Stokes's law in the case of spherical bodies we have
the radius of the sphere). In this case if
is constant the body will rapidly acquire a drift velocity:
so in the case of particles diffusing through a viscous fluids we have
On the other hand, in the case of a dilute fluid, supposing again that
is constant (or equivalently that the free mean path of the particles is much smaller than the scales over which
varies) we will have:
where we are only considering the drift due to
(we have neglected the possible term
is the velocity of the particle after a scattering, since it is isotropic and so it does not contribute to the net displacement of the particles). Now, since
and remembering also that
, so that the mean displacement after a scattering is
(the mean square velocity of the particles) we have:
where of course
is the kinetic energy of the particles.
Now that we also know what the mobility is, we can proceed to determine how the diffusion equation changes with the presence of the external and static force .
Consider an arbitrary function ; then:
where the presence of and the integration over are due to the fact that at time the coordinate will be equal to with probability (just like we have seen in the derivation of the diffusion equation in The diffusion equation).
Since and are small, we can expand in the right hand side:
Neglecting all the higher order terms and remembering that
, we get:
Integrating by parts:
vanish sufficiently fast at infinity, the last two terms are null so that we are left with:
is arbitrary, this equality holds if the remaining part of the integrands are equal; therefore, dividing also by
and taking the limits
, we have:
This is the so called Fokker-Planck equation
is Langevin's equation
Let us note that this is still compatible with a continuity equation, but this time the flow of particles is given (in three dimensions) by:
where the first term is the flow due to the external force, while the second one is the usual diffusive term due to the density inhomogeneities of the system.
Just to see the Fokker-Planck equation at work in a very simple case, let us determine its stationary solution when is a uniform gravitational force, with being the vertical coordinate (which is the only significant one; we assume that the positive direction is that pointing downwards); since (with the mass of the particles), then if:
is the gravitational potential.
- ↑ Of course, we could have equivalently solved the equation from a direct application of , namely from .