NARROW ESCAPE, PART I A. Singer ∗, Z. Schuss†, D. Holcman‡§, R. S. Eisenberg¶ August 12, 2018 Abstract A Brownian particle with diffusion coefficient D is confined to a bounded domain of volume V in R3 by a reflecting boundary, except for a small absorbing window. The mean time to absorption diverges as the win- dow shrinks, thus rendering the calculation of the mean escape time a singular perturbation problem. We construct an asymptotic approxi- mation for the case of an elliptical window of large semi axis a V 1/3 V ≪ and show that the mean escape time is Eτ K(e), where e is ∼ 2πDa the eccentricity of the ellipse; and K( ) is the complete elliptic integral · of the first kind. In the special case of a circular hole the result reduces V to Lord Rayleigh’s formula Eτ , which was derived by heuristic ∼ 4aD considerations. For the special case of a spherical domain, we obtain V a R a the asymptotic expansion Eτ = 1+ log + O . This 4aD R a R problem is important in understanding the flow of ions in and out of arXiv:math-ph/0412048v1 15 Dec 2004 narrow valves that control a wide range of biological and technological function. 1 Introduction We consider the exit problem of a Brownian motion from a bounded domain, whose boundary is reflecting, except for a small absorbing window. The mean ∗Department of Applied Mathematics, Tel-Aviv University, Ramat-Aviv, 69978 Tel-Aviv, Israel, e-mail: [email protected] †Department of Mathematics, Tel-Aviv University, Tel-Aviv 69978, Israel, e-mail: [email protected]. ‡Department of Mathematics, Weizmann Institute of Science, Rehovot 76100 Israel, e- mail [email protected]. §Keck Center, department of Physiology, UCSF, 513 Parnassus Ave, San Francisco 94143 USA, e-mail [email protected]. ¶Department of Molecular Biophysics and Physiology, Rush Medical Center, 1750 Har- rison St., Chicago, IL 60612, email: [email protected] 1 first passage time to the absorbing window (MFPT), Eτ, is the solution of a mixed Neumann-Dirichlet boundary value problem (BVP) for the Poisson equation, known as the corner problem, which has singularity at the boundary of the hole [1]-[3]. The MFPT grows to infinity as the window size shrinks to zero, thus rendering its calculation a singular perturbation problem for the mixed BVP, which we call the narrow escape problem. The narrow escape problem has been considered in the literature in only a few special cases, be- ginning with Lord Rayleigh (in the context of acoustics), who found the flux through a small hole by using a result of Helmholtz [6]. He stated [7] (p.176) “Among different kinds of channels an important place must be assigned to those consisting of simple apertures in unlimited plane walls of infinitesimal thickness. In practical applications it is sufficient that a wall be very thin in proportion to the dimensions of the aperture, and approximately plane within a distance from aperture large in proportion to the same quantity.” More recently, Rayleigh’s result was shown to fit the MFPT obtained from Brow- nian dynamics simulations [8]. Another result was presented in [9], where a two-dimensional narrow escape problem was considered and whose method is generalized here. A related problem is that of escape from a domain, whose boundary is reflecting, except for an absorbing sphere, disjoint from the reflect- ing part of the boundary [11, and references therein]. It differs from the narrow escape problem in that there is no singularity at the boundary and there is no boundary layer. The mixed boundary value problems of classical electrostatics (e.g., the electrified disk problem [12]), elasticity (punch problems), diffusion and conductance theory, hydrodynamics, and acoustics were solved, by and large, for special geometries by separation of variables. In axially symmetric geometries this method leads to a dual series or to integral equations that can be solved by special techniques [13]-[17]. The special case of asymptotic rep- resentation of the solution of the corner problem for small Dirichlet and large Neumann boundaries was not done for general domains. The first attempt in this direction seems to be [9]. The narrow escape problem does not seem to fall within the theory of large deviations [18]. It is different from Kolmogorov’s exit problem [19] of a diffusion process with small noise from an attractor of the drift (e.g., a stable equilibrium or limit cycle) in that the narrow escape problem has no large and small coefficients in the equation. The singularity of Kolmogorov’s problem is the degeneration of a second order elliptic operator into a first order operator in the limit of small noise, whereas the singularity of the narrow escape problem is the degeneration of the mixed BVP to a Neu- mann BVP on the entire boundary. There exist precise asymptotic expansions of Eτ for Kolmogorov’s exit problem, including error estimates (see, e.g., [4], [5]), which show that the MFPT grows exponentially with decreasing noise. In contrast, the narrow escape time grows algebraically rather than exponentially, as the window shrinks. Our first main result is a derivation of the leading or- der term in the expansion of the MFPT of a Brownian particle with diffusion 2 coefficient D, from a general domain of volume V to an elliptical hole of large semi axis a that is much smaller than V 1/3, V Eτ K(e), (1.1) ∼ 2πDa where e is the eccentricity of the ellipse, and K( ) is the complete elliptic integral of the first kind. In the special case of a circular· hole (1.1) reduces to V Eτ . (1.2) ∼ 4aD Equation (1.1) shows that the MFPT depends on the shape of the hole, and not just on its area. This result was known to Lord Rayleigh [7], who considered the problem of the electrified disk (which he knew was equivalent to finding the flow of an incompressible fluid through a channel and to the problem of finding the conductance of the channel), who reduced the problem to that of solving an integral equation for the flux density through the hole. The solution of the integral equation, which goes back to Helmholtz [6] and is discussed in [16], is proportional to (a2 ρ2)−1/2 in the circular case, where ρ is the distance from − the center of the hole [12]-[14]. Note that equations (1.1) and (1.2) are leading order approximations and do not contain an error estimate. We prove (1.1) by using the singularity properties of Neumann’s function for three-dimensional domains, in a manner similar to that used in [9] for two-dimensional problems. The leading order term is the solution of Helmholtz’s integral equation [6]. Our second main result is a derivation of the second term and error estimate for a ball of radius R with a small circular hole of radius a in the boundary, V a R a Eτ = 1+ log + O . (1.3) 4aD R a R Equation (1.3) contains both the second term in the asymptotic expansion of the MFPT and an error estimate. We use Collins’ method [20, 21] of solving dual series of equations and expand the resulting solutions for small ε = a/R. The estimate of the error term, which turns out to be O(ε log ε), seems to be a new result. An error estimate for eq.(1.1) for a general domain is still an open problem. We conjecture that it is O(ε log ε), as is the case for the ball. If the absorbing window touches a singular point of the boundary, such as a corner or cusp, the singularity of the Neumann function changes and so do the asymptotic results. In three dimensions the class of isolated singularities of the boundary is much richer than in the plane, so the methods of [22] cannot be generalized in a straightforward manner to three dimensions. We postpone the investigation of the MFPT to windows at isolated singular points in three dimensions to a future paper. In Section 2 we derive a leading order approximation to the MFPT for a general domain with a general small window. 3 The leading order term is expressed in terms of a solution to Helmholtz’s integral equation, which is solved explicitly for an elliptical window. In Section 3 we obtain two terms in the asymptotic expansion of the MFPT from a ball with a circular window and an error estimate. Finally, we present a summary and list some applications in Section 4. This is the first paper in a series of three, the second of which considers the narrow escape problem from a bounded simply connected planar domain, and the third of which considers the narrow escape problem from a bounded domain with boundary with corners and cusps on a two-dimensional Riemannian manifold. 2 General 3D bounded domain A Brownian particle diffuses freely in a bounded domain Ω R3, whose bound- ary ∂Ω is sufficiently smooth. The trajectory of the Brownian⊂ particle, denoted x(t), is reflected at the boundary, except for a small hole ∂Ωa, where it is ab- sorbed. The reflecting part of the boundary is ∂Ωr = ∂Ω ∂Ωa. The lifetime of the particle in Ω is the first passage time τ of the Brownian− particle from any point x Ω to the absorbing boundary ∂Ω . The MFPT, ∈ a v(x)= E[τ x(0) = x], | is finite under quite general conditions [24]. As the size (e.g., the diameter) of the absorbing hole decreases to zero, but that of the domain remains finite, we assume that the MFPT increases indefinitely.