Enhanced diffusion through surface excursion: a master-equation approach to the narrow-escape-time problem

Diffusion-Controlled Reactions: An Overview

We review the milestones in the century-long development of the theory of diffusion-controlled reactions. Starting from the seminal work by von Smoluchowski, who recognized the importance of diffusion in chemical reactions, we discuss perfect and imperfect surface reactions, their microscopic origins, and the underlying mathematical framework. Single-molecule reaction schemes, anomalous bulk diffusions, reversible binding/unbinding kinetics, and many other extensions are presented. An alternative encounter-based approach to diffusion-controlled reactions is introduced, with emphasis on its advantages and potential applications. Some open problems and future perspectives are outlined.

Brownian dynamics simulations for the narrow escape problem in the unit sphere

The narrow escape problem is a first-passage problem that concerns the calculation of the time needed for a Brownian particle to leave a domain with localized absorbing boundary traps, such that the measure of these traps is asymptotically small compared to the domain size. A common model for the mean first-passage time (MFPT) as a function of particle's starting location in a given domain with constant diffusivity is given by a Poisson partial differential equation subject to mixed Dirichlet-Neumann boundary conditions. The primary objective of this work is to perform direct numerical simulations of multiple particles undergoing Brownian motion in a three-dimensional spherical domain with boundary traps, compute MFPT values by averaging Brownian escape times, and compare these with explicit asymptotic results obtained previously by approximate solution of the Poisson problem. A close agreement of MFPT values is observed already at 10^{4} particle runs from a single starting point, providing a computational validation of the Poisson equation-based continuum model. Direct Brownian dynamics simulations are also used to study additional features of particle dynamics in narrow escape problems that cannot be captured in a continuum approach, such as average times spent by particles in a thin layer near the domain boundary, and effects of isotropic vs anisotropic near-boundary diffusion.

Paradigm Shift in Diffusion-Mediated Surface Phenomena

Diffusion-mediated surface phenomena are crucial for human life and industry, with examples ranging from oxygen capture by lung alveolar surface to heterogeneous catalysis, gene regulation, membrane permeation, and filtration processes. Their current description via diffusion equations with mixed boundary conditions is limited to simple surface reactions with infinite or constant reactivity. In this Letter, we propose a probabilistic approach based on the concept of boundary local time to investigate the intricate dynamics of diffusing particles near a reactive surface. Reformulating surface-particle interactions in terms of stopping conditions, we obtain in a unified way major diffusion-reaction characteristics such as the propagator, the survival probability, the first-passage time distribution, and the reaction rate. This general formalism allows us to describe new surface reaction mechanisms such as for instance surface reactivity depending on the number of encounters with the diffusing particle that can model the effects of catalyst fooling or membrane degradation. The disentanglement of the geometric structure of the medium from surface reactivity opens far-reaching perspectives for modeling, optimization, and control of diffusion-mediated surface phenomena.

Time-averaged mean square displacement for switching diffusion

We consider a classic two-state switching diffusion model from a single-particle tracking perspective. The mean and the variance of the time-averaged mean square displacement (TAMSD) are computed exactly. When the measurement time (i.e., the trajectory duration) is comparable to or smaller than the mean residence times in each state, the ergodicity breaking parameter is shown to take arbitrarily large values, suggesting an apparent weak ergodicity breaking for this ergodic model. In this regime, individual random trajectories are not representative while the related TAMSD curves exhibit a broad spread, in agreement with experimental observations in living cells and complex fluids. Switching diffusions can thus present, in some cases, an ergodic alternative to commonly used and inherently non-ergodic continuous-time random walks that capture similar features.

The escape problem for mortal walkers

We introduce and investigate the escape problem for random walkers that may eventually die, decay, bleach, or lose activity during their diffusion towards an escape or reactive region on the boundary of a confining domain. In the case of a first-order kinetics (i.e., exponentially distributed lifetimes), we study the effect of the associated death rate onto the survival probability, the exit probability, and the mean first passage time. We derive the upper and lower bounds and some approximations for these quantities. We reveal three asymptotic regimes of small, intermediate, and large death rates. General estimates and asymptotics are compared to several explicit solutions for simple domains and to numerical simulations. These results allow one to account for stochastic photobleaching of fluorescent tracers in bio-imaging, degradation of mRNA molecules in genetic translation mechanisms, or high mortality rates of spermatozoa in the fertilization process. Our findings provide a mathematical ground for optimizing storage containers and materials to reduce the risk of leakage of dangerous chemicals or nuclear wastes.

Diffusive escape through a narrow opening: new insights into a classic problem

We study the mean first exit time (T_ε) of a particle diffusing in a circular or a spherical micro-domain with an impenetrable confining boundary containing a small escape window (EW) of an angular size ε.

Accelerating search kinetics by following boundaries

We derive exact expressions of the mean first-passage time to a bulk target for a random searcher that performs boundary-mediated diffusion in a circular domain. Although nonintuitive for bulk targets, it is found that boundary excursions, if fast enough, can minimize the search time. A scaling analysis generalizes these findings to domains of arbitrary shapes and underlines their robustness. Overall, these results provide a generic mechanism of optimization of search kinetics in interfacial systems, which could have important implications in chemical physics. In the context of animal behavior sciences, it shows that following the boundaries of a domain can accelerate a search process, and therefore suggests that thigmotactism could be a kinetically efficient behavior.

Splitting probabilities and interfacial territory covered by two-dimensional and three-dimensional surface-mediated diffusion

We consider the mean territory covered by a particle that performs surface-mediated diffusion inside a spherical confining domain (in two and three dimensions) before exit through an opening on the surface. This quantity can be expressed in terms of the splitting probability between two targets on the surface. We derive a general formula that relates this splitting probability to the mean first passage time to a single target that has been recently calculated for such a surface-mediated diffusion process. This formula is exact for pointlike targets and is shown to be accurate for extended targets. The mean covered territory is then found and analyzed for an arbitrary extension of the exit region in both two- and three-dimensional spherical domains.

Enhanced transport through desorption-mediated diffusion

We present a master equation approach to the study of the bulk-mediated surface diffusion mechanism in a three-dimensional confined domain. The proposed scheme allowed us to evaluate analytically a number of magnitudes that were used to characterize the efficiency of the bulk-mediated surface transport mechanism, for instance, the mean escape time from the domain, and the mean number of distinct visited sites on the confined domain boundary.

Exact mean exit time for surface-mediated diffusion

We present an exact expression for the mean exit time through the cap of a confining sphere for particles alternating phases of surface and of bulk diffusion. The present approach is based on an integral equation which can be solved analytically. In contrast to the statement of Berezhkovskii and Barzykin [J. Chem. Phys. 136, 54115 (2012)], we show that the mean exit time can be optimized with respect to the desorption rate, under analytically determined criteria.

Narrow-escape-time problem: the imperfect trapping case

We present a master equation approach to the narrow escape time (NET) problem, i.e., the time needed for a particle contained in a confining domain with a single narrow opening to exit the domain for the first time. We introduce a finite transition probability, ν, at the narrow escape window, allowing the study of the imperfect trapping case. Ranging from 0 to ∞, ν allowed the study of both extremes of the trapping process: that of a highly deficient capture and situations where escape is certain ("perfect trapping" case). We have obtained analytic results for the basic quantity studied in the NET problem, the mean escape time, and we have studied its dependence in terms of the transition (desorption) probability over (from) the surface boundary, the confining domain dimensions, and the finite transition probability at the escape window. Particularly we show that the existence of a global minimum in the NET depends on the "imperfection" of the trapping process. In addition to our analytical approach, we have implemented Monte Carlo simulations, finding excellent agreement between the theoretical results and simulations.

Interfacial territory covered by surface-mediated diffusion

We consider a minimal model of heterogeneous catalysis in which a molecule performs surface-mediated diffusion inside a confining domain whose boundary contains catalytic sites. We explicitly take into account the combination of surface and bulk diffusion, and we obtain exact results for the mean and variance of the territory covered on the boundary by the particle before its exit in the case of a two-dimensional spherical domain. Depending on the relative positions of the entrance and exit points, very different behaviors with respect to the mean adsorption time of the molecule on the surface are found. We also determine both exact lower and upper bounds and an approximate expression of the probability of reacting with catalytic sites before exiting the domain. These results provide a quantitative measure of the efficiency of an idealized catalyst.