Narrow escape problem

Last updated January 15, 2024

The narrow escape problem^[1]^[2] is a ubiquitous problem in biology, biophysics and cellular biology.

The mathematical 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.^[3]^[4]^[5]^[6]^[7]^[8]^[9]

When escape is even more stringent due to severe geometrical restrictions at the place of escape, the narrow escape problem becomes the dire strait problem.^[10]^[11]

The narrow escape problem was proposed in the context of biology and biophysics by D. Holcman and Z. Schuss,^[12] and later on with A.Singer and led to the narrow escape theory in applied mathematics and computational biology.^[13]^[14]^[15]

Formulation

The motion of a particle is described by the Smoluchowski limit of the Langevin equation:^[16]^[17]

dX_{t}={\sqrt {2D}}\,dB_{t}+{\frac {1}{\gamma }}F(x)\,dt,

where $D$ is the diffusion coefficient of the particle, $\gamma$ is the friction coefficient per unit of mass, $F(x)$ the force per unit of mass, and $B_{t}$ is a Brownian motion.

Mean first passage time and the Fokker-Planck equation

A common question is to estimate the mean sojourn time of a particle diffusing in a bounded domain $\Omega$ before it escapes through a small absorbing window $\partial \Omega _{a}$ in its boundary $\partial \Omega$ . The time is estimated asymptotically in the limit ${\textstyle \varepsilon ={\frac {|\partial \Omega _{a}|}{|\partial \Omega |}}\ll 1}$

The probability density function (pdf) $p_{\varepsilon }(x,t)$ is the probability of finding the particle at position $x$ at time $t$ .

The pdf satisfies the Fokker–Planck equation:

{\frac {\partial }{\partial t}}p_{\varepsilon }(x,t)=D\Delta p_{\varepsilon }(x,t)-{\frac {1}{\gamma }}\nabla (p_{\varepsilon }(x,t)F(x))

with initial condition

p_{\varepsilon }(x,0)=\rho _{0}(x)\,

and mixed Dirichlet–Neumann boundary conditions ( $t>0$ )

p_{\varepsilon }(x,t)=0{\text{ for }}x\in \partial \Omega _{a}

D{\frac {\partial }{\partial n}}p_{\varepsilon }(x,t)-{\frac {p_{\varepsilon }(x,t)}{\gamma }}F(x)\cdot n(x)=0{\text{ for }}x\in \partial \Omega -\partial \Omega _{a}

The function

u_{\varepsilon }(y)=\int _{\Omega }\int _{0}^{\infty }p_{\varepsilon }(x,ty)\,dt\,dx

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

D\Delta u_{\varepsilon }(y)+{\frac {1}{\gamma }}F(y)\cdot \nabla u_{\varepsilon }(y)=-1

u_{\varepsilon }(y)=0{\text{ for }}y\in \partial \Omega _{a}

{\frac {\partial u_{\varepsilon }(y)}{\partial n}}=0{\text{ for }}y\in \partial \Omega _{r}

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

u_{\varepsilon }(y)={\frac {A}{\pi D}}\ln {\frac {1}{\varepsilon }}+O(1),

where $A$ is the surface of the domain. The function $u_{\epsilon }(y)$ does not depend on the initial position $y$ , except for a small boundary layer near the absorbing boundary due to the asymptotic form.

The first order term matters in dimension 2: for a circular disk of radius $R$ , the mean escape time of a particle starting in the center is

E(\tau |x(0)=0)={\frac {R^{2}}{D}}\left(\log \left({\frac {1}{\varepsilon }}\right)+\log 2+{\frac {1}{4}}+O(\varepsilon )\right).

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

E(\tau )={\frac {R^{2}}{D}}\left(\log \left({\frac {1}{\varepsilon }}\right)+\log 2+{\frac {1}{8}}+O(\varepsilon )\right).

The geometry of the small opening can affect the escape time: if the absorbing window is located at a corner of angle $\alpha$ , then:

E\tau ={\frac {|\Omega |}{\alpha D}}\left[\log {\frac {1}{\varepsilon }}+O(1)\right].

More surprising, near a cusp in a two dimensional domain, the escape time $E\tau$ grows algebraically, rather than logarithmically: in the domain bounded between two tangent circles, the escape time is:

E\tau ={\frac {|\Omega |}{(d-1)D}}\left({\frac {1}{\varepsilon }}+O(1)\right),

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 $\beta ={\frac {R_{1}}{R_{2}}}<1,$ 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}.

This equation contains two terms of the asymptotic expansion of $E\tau$ and $2\epsilon$ is the angle of the absorbing boundary. The case $\beta$ 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

F(x)\neq 0

the gap in the spectrum is not necessarily small between the first and the second eigenvalues, depending on the relative size of the small hole and the force barriers, the particle has to overcome in order to escape. The escape stream is not necessarily Poissonian.

Analytical results

A theorem that relates the Brownian motion escape problem to a (deterministic) partial differential equation problem is the following.

Theorem — Let $\Omega$ be a bounded domain with smooth boundary $\partial \Omega$ and $\Gamma$ be a closed subset of $\partial \Omega$ . For each $x\in \Omega$ , let $\tau _{x}$ be the first time of a particle hitting $\Gamma$ , assuming that the particle starts from $x$ , is subject to the Brownian motion in $\Omega$ , and reflects from $\partial \Omega$ . Then, the mean first passage time, $T(x):=\mathbb {E} [\tau _{x}]$ , and its variance, $v(x):=\mathbb {E} [(\tau _{x}-T(x))^{2}]$ , are solutions of the following boundary value problems:

-\Delta T=2{\text{ in }}\Omega ,~T=0{\text{ on }}\Gamma ,~\partial _{n}T=0{\text{ on }}\partial \Omega \setminus \Gamma

-\Delta v=2\vert \nabla T\vert ^{2}{\text{ in }}\Omega ,~v=0{\text{ on }}\Gamma ,~\partial _{n}v=0{\text{ on }}\partial \Omega \setminus \Gamma

Here $\partial _{n}:=n\cdot \nabla$ is the derivative in the direction $n$ , the exterior normal to $\partial \Omega .$ Moreover, the average of the variance can be calculated from the formula

{\bar {v}}:={\frac {1}{\vert \Omega \vert }}\int _{\Omega }v(x)dx={\frac {1}{\vert \Omega \vert }}\int _{\Omega }T^{2}(x)dx=:T^{2}

The first part of the theorem is a classical result, while the average variance was proved in 2011 by Carey Caginalp and Xinfu Chen.^[18]^[19]^[20]

The escape time has been the subject of a number of studies using the small gate as an asymptotically small parameter. The following closed form result^[18]^[19]^[20] gives an exact solution that confirms these asymptotic formulae and extends them to gates that are not necessarily small.

Theorem (Carey Caginalp and Xinfu Chen Closed Formula) — In 2-D, with points identified by complex numbers, let

{\displaystyle \Omega

Then the mean first passage time $T(z)$ , for $z\in {\bar {\Omega }}$ , is given by

T(z)={\frac {1-\vert z\vert ^{2}}{2}}+2\log {\left|{\frac {1-z+{\sqrt {(1-ze^{-i\varepsilon })(1-ze^{i\varepsilon })}}}{2\sin {\frac {\varepsilon }{2}}}}\right|}

Another set of results concerns the probability density of the location of exit.^[19]

Theorem (Carey Caginalp and Xinfu Chen Probability Density) — The probability density of the location of a particle at time of its exit is given by

{\bar {j}}(e^{i\theta }):=-{\frac {1}{2\pi }}{\frac {\partial }{\partial r}}T(e^{i\theta })={\begin{cases}0,&{\text{if }}\varepsilon <\theta <2\pi -\varepsilon \\{\frac {1}{2\pi }}{\frac {\cos {\frac {\theta }{2}}}{\sqrt {\sin ^{2}{\frac {\varepsilon }{2}}-\sin ^{2}{\frac {\theta }{2}}}}},&{\text{if }}\vert \theta \vert <\varepsilon \end{cases}}

That is, for any (Borel set) $\gamma \subset \partial \Omega$ , the probability that a particle, starting either at the origin or uniformly distributed in $\Omega$ , exhibiting Brownian motion in $\Omega$ , reflecting when it hits $\partial \Omega \setminus \Gamma$ , and escaping once it hits $\Gamma$ , ends up escaping from $\gamma$ is

P(\gamma )=\int _{\gamma }{\bar {j}}(y)dS_{y}

where $dS_{y}$ is the surface element of $\partial \Omega$ at $y\in \partial \Omega$ .

Simulations of Brownian motion escape

In simulation there is a random error due to the statistical sampling process. This error can be limited by appealing to the central limit theorem and using a large number of samples. There is also a discretization error due to the finite size approximation of the step size in approximating the Brownian motion. One can then obtain empirical results as step size and gate size vary. Using the exact result quoted above for the particular case of the circle, it is possible to make a careful comparison of the exact solution with the numerical solution.^[21]^[22] This illuminates the distinction between finite steps and continuous diffusion. A distribution of exit locations was also obtained through simulations for this problem.

Biological applications

Stochastic chemical reactions in microdomains

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.^[23]

Related Research Articles

The propagation constant of a sinusoidal electromagnetic wave is a measure of the change undergone by the amplitude and phase of the wave as it propagates in a given direction. The quantity being measured can be the voltage, the current in a circuit, or a field vector such as electric field strength or flux density. The propagation constant itself measures the dimensionless change in magnitude or phase per unit length. In the context of two-port networks and their cascades, propagation constant measures the change undergone by the source quantity as it propagates from one port to the next.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

<span class="mw-page-title-main">Fermi gas</span> Physical model of gases composed of many non-interacting identical fermions

A Fermi gas is an idealized model, an ensemble of many non-interacting fermions. Fermions are particles that obey Fermi–Dirac statistics, like electrons, protons, and neutrons, and, in general, particles with half-integer spin. These statistics determine the energy distribution of fermions in a Fermi gas in thermal equilibrium, and is characterized by their number density, temperature, and the set of available energy states. The model is named after the Italian physicist Enrico Fermi.

In continuum mechanics, the infinitesimal strain theory is a mathematical approach to the description of the deformation of a solid body in which the displacements of the material particles are assumed to be much smaller than any relevant dimension of the body; so that its geometry and the constitutive properties of the material at each point of space can be assumed to be unchanged by the deformation.

In vector calculus, Green's theorem relates a line integral around a simple closed curve $C$ to a double integral over the plane region $D$ bounded by $C$ . It is the two-dimensional special case of Stokes' theorem.

In quantum mechanics and quantum field theory, the propagator is a function that specifies the probability amplitude for a particle to travel from one place to another in a given period of time, or to travel with a certain energy and momentum. In Feynman diagrams, which serve to calculate the rate of collisions in quantum field theory, virtual particles contribute their propagator to the rate of the scattering event described by the respective diagram. These may also be viewed as the inverse of the wave operator appropriate to the particle, and are, therefore, often called (causal) Green's functions.

Geometrical optics, or ray optics, is a model of optics that describes light propagation in terms of rays. The ray in geometrical optics is an abstraction useful for approximating the paths along which light propagates under certain circumstances.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

In quantum physics, Fermi's golden rule is a formula that describes the transition rate from one energy eigenstate of a quantum system to a group of energy eigenstates in a continuum, as a result of a weak perturbation. This transition rate is effectively independent of time and is proportional to the strength of the coupling between the initial and final states of the system as well as the density of states. It is also applicable when the final state is discrete, i.e. it is not part of a continuum, if there is some decoherence in the process, like relaxation or collision of the atoms, or like noise in the perturbation, in which case the density of states is replaced by the reciprocal of the decoherence bandwidth.

In quantum mechanics, a two-state system is a quantum system that can exist in any quantum superposition of two independent quantum states. The Hilbert space describing such a system is two-dimensional. Therefore, a complete basis spanning the space will consist of two independent states. Any two-state system can also be seen as a qubit.

The Duffing equation, named after Georg Duffing (1861–1944), is a non-linear second-order differential equation used to model certain damped and driven oscillators. The equation is given by

In mathematics, in particular in algebraic geometry and differential geometry, Dolbeault cohomology (named after Pierre Dolbeault) is an analog of de Rham cohomology for complex manifolds. Let M be a complex manifold. Then the Dolbeault cohomology groups $depend on a pair of integers p and q and are realized as a subquotient of the space of complex differential forms of degree (p, q).$

In mathematics, the Möbius energy of a knot is a particular knot energy, i.e., a functional on the space of knots. It was discovered by Jun O'Hara, who demonstrated that the energy blows up as the knot's strands get close to one another. This is a useful property because it prevents self-intersection and ensures the result under gradient descent is of the same knot type.

In the fields of nonlinear optics and fluid dynamics, modulational instability or sideband instability is a phenomenon whereby deviations from a periodic waveform are reinforced by nonlinearity, leading to the generation of spectral-sidebands and the eventual breakup of the waveform into a train of pulses.

In continuum mechanics, a compatible deformation tensor field in a body is that unique tensor field that is obtained when the body is subjected to a continuous, single-valued, displacement field. Compatibility is the study of the conditions under which such a displacement field can be guaranteed. Compatibility conditions are particular cases of integrability conditions and were first derived for linear elasticity by Barré de Saint-Venant in 1864 and proved rigorously by Beltrami in 1886.

In mathematical physics and the theory of partial differential equations, the solitary wave solution of the form $is said to be orbitally stable if any solution with the initial data sufficiently close to forever remains in a given small neighborhood of the trajectory of$

In mathematics, the Kodaira–Spencer map, introduced by Kunihiko Kodaira and Donald C. Spencer, is a map associated to a deformation of a scheme or complex manifold X, taking a tangent space of a point of the deformation space to the first cohomology group of the sheaf of vector fields on X.

Phase reduction is a method used to reduce a multi-dimensional dynamical equation describing a nonlinear limit cycle oscillator into a one-dimensional phase equation. Many phenomena in our world such as chemical reactions, electric circuits, mechanical vibrations, cardiac cells, and spiking neurons are examples of rhythmic phenomena, and can be considered as nonlinear limit cycle oscillators.

The redundancy principle in biology expresses the need of many copies of the same entity to fulfill a biological function. Examples are numerous: disproportionate numbers of spermatozoa during fertilization compared to one egg, large number of neurotransmitters released during neuronal communication compared to the number of receptors, large numbers of released calcium ions during transient in cells, and many more in molecular and cellular transduction or gene activation and cell signaling. This redundancy is particularly relevant when the sites of activation are physically separated from the initial position of the molecular messengers. The redundancy is often generated for the purpose of resolving the time constraint of fast-activating pathways. It can be expressed in terms of the theory of extreme statistics to determine its laws and quantify how the shortest paths are selected. The main goal is to estimate these large numbers from physical principles and mathematical derivations.

Hybrid stochastic simulations are a sub-class of stochastic simulations. These simulations combine existing stochastic simulations with other stochastic simulations or algorithms. Generally they are used for physics and physics-related research. The goal of a hybrid stochastic simulation varies based on context, however they typically aim to either improve accuracy or reduce computational complexity. The first hybrid stochastic simulation was developed in 1985.

References

↑ Schuss, Z.; Singer, A.; Holcman, D. (2007-09-27). "The narrow escape problem for diffusion in cellular microdomains". Proceedings of the National Academy of Sciences. Proceedings of the National Academy of Sciences USA. 104 (41): 16098–16103. Bibcode:2007PNAS..10416098S. doi: 10.1073/pnas.0706599104 . ISSN 0027-8424. PMC 1994903 . PMID 17901203.
↑ D Holcman, Z Schuss, The narrow escape problem SIAM Review 56 (2), 213-257 (2014)
↑ Singer, A.; Schuss, Z.; Holcman, D. (2008-11-14). "Narrow escape and leakage of Brownian particles". Physical Review E. American Physical Society (APS). 78 (5): 051111. arXiv: 0808.2288 . Bibcode:2008PhRvE..78e1111S. doi:10.1103/physreve.78.051111. ISSN 1539-3755. PMID 19113099. S2CID 8739640.
↑ M. J. Ward, S. Pillay, A. Peirce, and T. Kolokolnikov An Asymptotic Analysis of the Mean First Passage Time for Narrow Escape Problems: Part I: Two-Dimensional Domains
↑ Holcman, D; Schuss, Z (2008-04-02). "Diffusion escape through a cluster of small absorbing windows". Journal of Physics A: Mathematical and Theoretical. IOP Publishing. 41 (15): 155001. Bibcode:2008JPhA...41o5001H. doi:10.1088/1751-8113/41/15/155001. ISSN 1751-8113. S2CID 4179599.
↑ Holcman, D., & Schuss, Z. (2015). Stochastic Narrow Escape in Molecular and Cellular Biology: Analysis and Applications. Springer.
↑ Cheviakov, Alexei F.; Ward, Michael J.; Straube, Ronny (2010). "An Asymptotic Analysis of the Mean First Passage Time for Narrow Escape Problems: Part II: The Sphere". Multiscale Modeling & Simulation. Society for Industrial & Applied Mathematics (SIAM). 8 (3): 836–870. doi:10.1137/100782620. hdl: 11858/00-001M-0000-0013-908F-6 . ISSN 1540-3459.
↑ Cheviakov, Alexei F.; Zawada, Daniel (2013-04-22). "Narrow-escape problem for the unit sphere: Homogenization limit, optimal arrangements of large numbers of traps, and the N² conjecture". Physical Review E. American Physical Society (APS). 87 (4): 042118. Bibcode:2013PhRvE..87d2118C. doi:10.1103/physreve.87.042118. ISSN 1539-3755. PMID 23679384.
↑ Coombs, Daniel; Straube, Ronny; Ward, Michael (2009). "Diffusion on a Sphere with Localized Traps: Mean First Passage Time, Eigenvalue Asymptotics, and Fekete Points". SIAM Journal on Applied Mathematics. Society for Industrial & Applied Mathematics (SIAM). 70 (1): 302–332. doi:10.1137/080733280. hdl: 11858/00-001M-0000-0013-9335-3 . ISSN 0036-1399.
↑ D. Holcman Z. Schuss, The dire strait time, SIAM Multiscale Modeling and simulations, 10(4), 1204–1231.
↑ Holcman, D; Schuss, Z (2013-06-20). "Control of flux by narrow passages and hidden targets in cellular biology". Reports on Progress in Physics. IOP Publishing. 76 (7): 074601. Bibcode:2013RPPh...76g4601H. doi:10.1088/0034-4885/76/7/074601. ISSN 0034-4885. PMID 23787818. S2CID 2102724.
↑ Holcman, D.; Schuss, Z. (2004). "Escape Through a Small Opening: Receptor Trafficking in a Synaptic Membrane". Journal of Statistical Physics. Springer Science and Business Media LLC. 117 (5–6): 975–1014. Bibcode:2004JSP...117..975H. doi:10.1007/s10955-004-5712-8. ISSN 0022-4715. S2CID 6324415.
↑ Singer, A.; Schuss, Z.; Holcman, D.; Eisenberg, R. S. (2006-01-20). "Narrow Escape, Part I". Journal of Statistical Physics. Springer Science and Business Media LLC. 122 (3): 437–463. arXiv: math-ph/0412048 . Bibcode:2006JSP...122..437S. doi:10.1007/s10955-005-8026-6. ISSN 0022-4715. S2CID 14014727.
↑ Singer, A.; Schuss, Z.; Holcman, D. (2006-01-20). "Narrow Escape, Part II: The Circular Disk". Journal of Statistical Physics. Springer Science and Business Media LLC. 122 (3): 465–489. arXiv: math-ph/0412050 . Bibcode:2006JSP...122..465S. doi:10.1007/s10955-005-8027-5. ISSN 0022-4715. S2CID 15765954.
↑ Singer, A.; Schuss, Z.; Holcman, D. (2006-01-20). "Narrow Escape, Part III: Non-Smooth Domains and Riemann Surfaces". Journal of Statistical Physics. Springer Science and Business Media LLC. 122 (3): 491–509. Bibcode:2006JSP...122..491S. doi:10.1007/s10955-005-8028-4. ISSN 0022-4715. S2CID 12317568.
↑ Z. Schuss, Theory and Applications of Stochastic Differential Equations (Wiley Series in Probability and Statistics - (1980)
↑ Z. Schuss, Theory and Applications of Stochastic Processes. An Analytical Approach. Series: Applied Mathematical Sciences, Vol. 170.
1 2 Caginalp, Carey; Chen, Xinfu (2011-02-01). "Analytical and numerical results for first escape time in 2D". Comptes Rendus Mathématique. 349 (3–4): 191–194. doi:10.1016/j.crma.2010.11.024. ISSN 1631-073X.
1 2 3 Chen, Xinfu; Caginalp, Carey (2012-01-01). "Analytical and Numerical Results for an Escape Problem". Archive for Rational Mechanics and Analysis. 203 (1): 329–342. Bibcode:2012ArRMA.203..329C. doi:10.1007/s00205-011-0455-6. ISSN 1432-0673. S2CID 32394342.
1 2 Caginalp, Carey (2011). Analytical and numerical results on escape (B. Phil. Thesis). University of Pittsburgh.
↑ Hughes, Aoife; Morris, Richard; Tomkins, Melissa (2020-03-31). "PyEscape: A narrow escape problem simulator package for Python". Journal of Open Source Software. 5 (47): 2072. Bibcode:2020JOSS....5.2072H. doi: 10.21105/joss.02072 . ISSN 2475-9066.
↑ Hughes, Aoife; Faulkner, Christine; Morris, Richard J.; Tomkins, Melissa (2021). "Intercellular communication as a series of narrow escape problems". IEEE Transactions on Molecular, Biological and Multi-Scale Communications. 7 (2): 89–93. doi: 10.1109/TMBMC.2021.3083719 . ISSN 2332-7804.
↑ Holcman, D.; Schuss, Z. (2005-03-15). "Stochastic chemical reactions in microdomains". The Journal of Chemical Physics. AIP Publishing. 122 (11): 114710. arXiv: math-ph/0412089 . Bibcode:2005JChPh.122k4710H. doi:10.1063/1.1849155. ISSN 0021-9606. PMID 15836246. S2CID 845444.

External links

Applied Mathematics and Computational Biology in Ecole Normale Superieure, Paris
Carey Caginalp publications and lectures http://www.pitt.edu/~careycag/
Comptes Rendus paper http://www.pitt.edu/~careycag/paper1.pdf
ARMA paper http://www.pitt.edu/~careycag/paper2.pdf

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Schuss, Z.; Singer, A.; Holcman, D. (2007-09-27). "The narrow escape problem for diffusion in cellular microdomains". Proceedings of the National Academy of Sciences. Proceedings of the National Academy of Sciences USA. 104 (41): 16098–16103. Bibcode:2007PNAS..10416098S. doi: 10.1073/pnas.0706599104 . ISSN 0027-8424. PMC 1994903 . PMID 17901203.

[2] D Holcman, Z Schuss, The narrow escape problem SIAM Review 56 (2), 213-257 (2014)

[Singer2008-3] Singer, A.; Schuss, Z.; Holcman, D. (2008-11-14). "Narrow escape and leakage of Brownian particles". Physical Review E. American Physical Society (APS). 78 (5): 051111. arXiv: 0808.2288 . Bibcode:2008PhRvE..78e1111S. doi:10.1103/physreve.78.051111. ISSN 1539-3755. PMID 19113099. S2CID 8739640.

[4] M. J. Ward, S. Pillay, A. Peirce, and T. Kolokolnikov An Asymptotic Analysis of the Mean First Passage Time for Narrow Escape Problems: Part I: Two-Dimensional Domains

[5] Holcman, D; Schuss, Z (2008-04-02). "Diffusion escape through a cluster of small absorbing windows". Journal of Physics A: Mathematical and Theoretical. IOP Publishing. 41 (15): 155001. Bibcode:2008JPhA...41o5001H. doi:10.1088/1751-8113/41/15/155001. ISSN 1751-8113. S2CID 4179599.

[6] Holcman, D., & Schuss, Z. (2015). Stochastic Narrow Escape in Molecular and Cellular Biology: Analysis and Applications. Springer.

[7] Cheviakov, Alexei F.; Ward, Michael J.; Straube, Ronny (2010). "An Asymptotic Analysis of the Mean First Passage Time for Narrow Escape Problems: Part II: The Sphere". Multiscale Modeling & Simulation. Society for Industrial & Applied Mathematics (SIAM). 8 (3): 836–870. doi:10.1137/100782620. hdl: 11858/00-001M-0000-0013-908F-6 . ISSN 1540-3459.

[8] Cheviakov, Alexei F.; Zawada, Daniel (2013-04-22). "Narrow-escape problem for the unit sphere: Homogenization limit, optimal arrangements of large numbers of traps, and the N² conjecture". Physical Review E. American Physical Society (APS). 87 (4): 042118. Bibcode:2013PhRvE..87d2118C. doi:10.1103/physreve.87.042118. ISSN 1539-3755. PMID 23679384.

[9] Coombs, Daniel; Straube, Ronny; Ward, Michael (2009). "Diffusion on a Sphere with Localized Traps: Mean First Passage Time, Eigenvalue Asymptotics, and Fekete Points". SIAM Journal on Applied Mathematics. Society for Industrial & Applied Mathematics (SIAM). 70 (1): 302–332. doi:10.1137/080733280. hdl: 11858/00-001M-0000-0013-9335-3 . ISSN 0036-1399.

[10] D. Holcman Z. Schuss, The dire strait time, SIAM Multiscale Modeling and simulations, 10(4), 1204–1231.

[11] Holcman, D; Schuss, Z (2013-06-20). "Control of flux by narrow passages and hidden targets in cellular biology". Reports on Progress in Physics. IOP Publishing. 76 (7): 074601. Bibcode:2013RPPh...76g4601H. doi:10.1088/0034-4885/76/7/074601. ISSN 0034-4885. PMID 23787818. S2CID 2102724.

[12] Holcman, D.; Schuss, Z. (2004). "Escape Through a Small Opening: Receptor Trafficking in a Synaptic Membrane". Journal of Statistical Physics. Springer Science and Business Media LLC. 117 (5–6): 975–1014. Bibcode:2004JSP...117..975H. doi:10.1007/s10955-004-5712-8. ISSN 0022-4715. S2CID 6324415.

[13] Singer, A.; Schuss, Z.; Holcman, D.; Eisenberg, R. S. (2006-01-20). "Narrow Escape, Part I". Journal of Statistical Physics. Springer Science and Business Media LLC. 122 (3): 437–463. arXiv: math-ph/0412048 . Bibcode:2006JSP...122..437S. doi:10.1007/s10955-005-8026-6. ISSN 0022-4715. S2CID 14014727.

[14] Singer, A.; Schuss, Z.; Holcman, D. (2006-01-20). "Narrow Escape, Part II: The Circular Disk". Journal of Statistical Physics. Springer Science and Business Media LLC. 122 (3): 465–489. arXiv: math-ph/0412050 . Bibcode:2006JSP...122..465S. doi:10.1007/s10955-005-8027-5. ISSN 0022-4715. S2CID 15765954.

[15] Singer, A.; Schuss, Z.; Holcman, D. (2006-01-20). "Narrow Escape, Part III: Non-Smooth Domains and Riemann Surfaces". Journal of Statistical Physics. Springer Science and Business Media LLC. 122 (3): 491–509. Bibcode:2006JSP...122..491S. doi:10.1007/s10955-005-8028-4. ISSN 0022-4715. S2CID 12317568.

[16] Z. Schuss, Theory and Applications of Stochastic Differential Equations (Wiley Series in Probability and Statistics - (1980)

[17] Z. Schuss, Theory and Applications of Stochastic Processes. An Analytical Approach. Series: Applied Mathematical Sciences, Vol. 170.

[:0-18] 1 2 Caginalp, Carey; Chen, Xinfu (2011-02-01). "Analytical and numerical results for first escape time in 2D". Comptes Rendus Mathématique. 349 (3–4): 191–194. doi:10.1016/j.crma.2010.11.024. ISSN 1631-073X.

[:1-19] 1 2 3 Chen, Xinfu; Caginalp, Carey (2012-01-01). "Analytical and Numerical Results for an Escape Problem". Archive for Rational Mechanics and Analysis. 203 (1): 329–342. Bibcode:2012ArRMA.203..329C. doi:10.1007/s00205-011-0455-6. ISSN 1432-0673. S2CID 32394342.

[:2-20] 1 2 Caginalp, Carey (2011). Analytical and numerical results on escape (B. Phil. Thesis). University of Pittsburgh.

[21] Hughes, Aoife; Morris, Richard; Tomkins, Melissa (2020-03-31). "PyEscape: A narrow escape problem simulator package for Python". Journal of Open Source Software. 5 (47): 2072. Bibcode:2020JOSS....5.2072H. doi: 10.21105/joss.02072 . ISSN 2475-9066.

[22] Hughes, Aoife; Faulkner, Christine; Morris, Richard J.; Tomkins, Melissa (2021). "Intercellular communication as a series of narrow escape problems". IEEE Transactions on Molecular, Biological and Multi-Scale Communications. 7 (2): 89–93. doi: 10.1109/TMBMC.2021.3083719 . ISSN 2332-7804.

[23] Holcman, D.; Schuss, Z. (2005-03-15). "Stochastic chemical reactions in microdomains". The Journal of Chemical Physics. AIP Publishing. 122 (11): 114710. arXiv: math-ph/0412089 . Bibcode:2005JChPh.122k4710H. doi:10.1063/1.1849155. ISSN 0021-9606. PMID 15836246. S2CID 845444.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]

[21]

[22]

[23]