| This article is an orphan, as few or no other articles link to it. Please introduce links to this page from related articles; suggestions may be available. (January 2011) |
|
|
This article includes a list of references, related reading or external links, but its sources remain unclear because it lacks inline citations. Please improve this article by introducing more precise citations. (April 2011) |
The narrow escape problem is a ubiquitous problem in biology, biophysics and cellular biology.
The formulation is the following: a Brownian particle (ion, molecule, or protein) is confined to a bounded domain (a compartment or a cell) by a reflecting boundary, except for a small window through which it can escape. The narrow escape problem is that of calculating the mean escape time. This time diverges as the window shrinks, thus rendering the calculation a singular perturbation problem.
|
Contents
|
The motion of a particle is described by the Smoluchowski limit of the Langevin equation:

where
is the diffusion coefficient of the particle,
is the friction coefficient per unit of mass,
the force per unit of mass, and
is a Brownian motion.
A common question is to estimate the mean sojourn time of a particle diffusing in a bounded domain
before it escapes through a small absorbing window
in its boundary
. The time is estimated asymptotically in the limit 
The probability density function (pdf)
is the probability of finding the particle at position
at time
.
The pdf satisfies the Fokker–Planck equation

with initial condition

and mixed Dirichlet–Neumann boundary conditions (
)


The function

represents the mean sojourn time of particle, conditioned on the initial position
. It is the solution of the boundary value problem



The solution depends on the dimension of the domain. For a particle diffusing on a disk

where
is the surface of the domain. The function
does not depend on the initial position
, except for a small boundary layer near due to the asymptotic form. The first order term matters in dimension 2. For a circular disk of radius
, the mean escape time of a particle starting in the center is

The escape time averaged with respect to a uniform initial distribution of the particle is given by

The geometry of the small opening can affect the escape time: if the absorbing window is located at a corner of angle
, then
![E\tau = \frac{|\Omega|}{\alpha D} \left[\log \frac{1}{\varepsilon} +O(1)\right].](http://wpcontent.answcdn.com/wikipedia/en/math/c/5/4/c542fd806b0cace1a6d614746c53f05a.png)
More surprising, near a cusp in a two dimensional domain, the escape time
grows algebraically, rather than logarithmically: in the domain bounded between two tangent circles, the escape time is

where d > 1 is the ratio of the radii. Finally, when the domain is an annulus, the escape time to a small opening located on the inner circle involves a second parameter which is
the ratio of the inner to the outer radii, the escape time, averaged with respect to a uniform initial distribution, is
![E\tau = \frac{(R_2^2-R_1^2)}D\left[\log \frac{1}{\varepsilon} +
\log 2 + 2\beta^2 \right] +\frac{1}{2}\frac{R_2^2}{1-\beta^2}\log\frac{1}{\beta}- \frac{1}{4}R_2^2 +
O(\varepsilon,\beta^4)R_2^2.](http://wpcontent.answcdn.com/wikipedia/en/math/0/a/0/0a02c659296cf47a6469172e5fad1d4a.png)
This equation contains two terms of the asymptotic expansion of
and
is the angle of the absorbing boundary. The case
close to 1 remains open, and for general domains, the asymptotic expansion of the escape time remains an open problem. So does the problem of computing the escape time near a cusp point in three-dimensional domains. For Brownian motion in a field of force

the gap in the spectrum is not necessarily between the first and the second eigenvalues, depending on the relative size of the small hole and force barriers the particle has to overcome in order to escape. The escape stream is not necessarily Poissonian.
The forward rate of chemical reactions is the reciprocal of the narrow escape time, which generalizes the classical Smoluchowski formula for Brownian particles located in an infinite medium. A Markov description can be used to estimate the binding and unbinding to a small number of sites.
This entry is from Wikipedia, the leading user-contributed encyclopedia. It may not have been reviewed by professional editors (see full disclaimer)