280 likes | 307 Views
New Insights Into the (Old) Narrow Escape Problem. by Denis Grebenkov. Laboratoire de Physique de la Mati ère Condensée , CNRS – Ecole Polytechnique, 91128 Palaiseau France. Saint Petersburg, Russia, 22/09/2017. Outline. Overview. Part I: Role of geometry.
E N D
New Insights Into the (Old) Narrow Escape Problem by Denis Grebenkov Laboratoire de Physique de la Matière Condensée, CNRS – Ecole Polytechnique, 91128 Palaiseau France Saint Petersburg, Russia, 22/09/2017
Outline Overview Part I: Role of geometry Part II: Role of interactions and barrier All related publications can be found on my webpage: http://pmc.polytechnique.fr/pagesperso/dg
First passage and exit problems S. Redner, A guide to first-passage processes (2001) R. Metzler, et al., First-passage phenomena and their applications (2014) D. Holcman, Z. Schuss, SIAM Rev. 56, 213-257 (2014)
First passage and exit problems Target size/shape/location/type Domain shape Starting point location S. Redner, A guide to first-passage processes (2001) R. Metzler, et al., First-passage phenomena and their applications (2014) D. Holcman, Z. Schuss, SIAM Rev. 56, 213-257 (2014) Type of diffusive motion
First passage time (FPT) Survival probability MFPT: Mixed boundary value problem Global MFPT: on target on the rest
Escape problems Lord Rayleigh (1877): exit from a sphere through a circular hole of radius a Singer et al. (2006): Exact solution for the disk: 2D manifolds: Asymptotic behavior: Rayleigh, The theory of sound (1877); Singer et al., J. Stat. Phys. (2006)
Various extensions Multiple targets Partial adsorption on target/remaining part Gating problem Interaction with the surface Anomalous and intermittent diffusion Probability distribution Splitting probability First passage to targets inside the domain Optimality of the MFPT R. Metzler, et al., First-passage phenomena and their applications (2014) D. Holcman, Z. Schuss, SIAM Rev. 56, 213-257 (2014)
Summary: current paradigm Narrow escape limit: Global MFPT (2D) is proportional to and diverges as ln(1/) We will show that this paradigm is incomplete and can be misleading
Part I Arbitrary planar domains DG, PRL 117, 260201 (2016)
Two basic counter-examples Common picture: L h (x,y) (x,y) is independent of y and h and does not diverge as T is large while is large T is small while is small These are not “regular” domains:
Conformal mapping D x0 is the harmonic measure of seen from x0 on on on the rest on the rest
Integral formula for the MFPT eventual reflections MFPT from x0 to the whole boundary
Asymptotic behavior of the MFPT O(1) We retrieve the leading logarithmic behavior, with the perimeter replaced by the harmonic measure
Two basic counter-examples New picture: L h (x,y) (x,y) Makarov’s theorem: Makarov, Proc. London Math. Soc. 51, 369 (1985)
Illustration for a nontrivial domain We fix x0 and and move the region 2 3 5 15 4 MFPT 1 W p w x0 | |/ log(1/ ) 0 10 W p e | |/ log(1/ ) 7 10 6 T 8 9 5 0 0 2 4 6 8 10 s
Conclusions Integral formula for the MFPT in planar domains Dependence on the starting point Shape of the domain through conformal mapping Accounting for space-dependent diffusivity Fundamental relation between the MFPT and the HM, which substitutes the normalized perimeter Simple logarithmic term with the harmonic measure Scaling with the area is not universal
Part II Long-range interactions with boundary and escape barrier G. Oshanin Laboratory of Condensed Matter Theoretical Physics, University Paris-6 DG, G. Oshanin, PCCP 19, 2723 (2017)
Circular and spherical domains Mixed boundary value problem on target on the rest on target Two ingredients: entropic/energetic barrier long-range radial potential on the rest
Self-consistent approximation Mixed Robin-Neumann boundary condition on target on the rest on target on the rest Inhomogeneous Neumann boundary condition Note: solution of the modified problem is not unique! But one can get the leading (divergent) term
Main result MFPT to fully absorbing whole boundary Contribution due to partial escape region Contribution due to partial reflections
Accuracy of the approximation SCA asympt 1 10 FEM MC 2 /R 0 e 10 D T -1 10 0 0.5 1 1.5 2 2.5 3 e (no long-range interaction case)
Accuracy of the approximation 150 e SCA ( = 0.1) e FEM ( = 0.1) e SCA ( = 0.2) 100 e 2 FEM ( = 0.2) /R e p SCA ( = /4) e e p D T FEM ( = /4) 50 0 0 1 2 3 10 10 10 10 k R/D (no long-range interaction case)
Narrow escape limit MFPT to fully absorbing whole boundary Contribution due to partial escape region Contribution due to partial reflections Conventional scaling In presence of entropic/energetic barrier, the escape process is “barrier-limited”, not diffusion-limited
Role of long-range interactions Toy model: linear radial potential with a finite extent For this model, one can explicitly find the radial functions satisfying Attractive potential (U0<0) Repulsive potential (U0>0) MFPT is decreased MFPT is increased
Optimal range of the potential 60 50 40 2 /R (3) 30 e DT 20 U = -1 0 U = -2 10 0 U = -5 0 0 0 0.2 0.4 0.6 0.8 1 r /R 0 The optimal range lies near the boundary
Conclusions We proposed an accurate approximation to MFPT for rotation-invariant domains The escape process is “barrier-limited”, the asymptotics being , not The narrow escape behavior is only determined by the potential and its derivative at the boundary The range of attractive potentials can be optimized to minimize the MFPT
Theory and Modeling of Complex Systems in Life Sciences We thank all the participants!!! Scientific Organizers Local Organizers Denis Grebenkov Sergey Nechaev Stanislav Smirnov Nadia Zalesskaya Tatiana Vinogradova Natalia Kirshner Liza Kruykova http://inadilic.fr/conference-2017/ Leonhard Euler (1707-1783) We acknowledge the financial support by