Narrow capture problem: An encounter-based approach to partially reactive targets.

Abstract

A general topic of current interest is the analysis of diffusion problems in singularly perturbed domains with small interior targets or traps (the narrow capture problem). One major application is to intracellular diffusion, where the targets typically represent some form of reactive biochemical substrate. Most studies of the narrow capture problem treat the target boundaries as totally absorbing (Dirichlet), that is, the chemical reaction occurs immediately on first encounter between particle and target surface. In this paper, we analyze the three-dimensional narrow capture problem in the more realistic case of partially reactive target boundaries. We begin by considering classical Robin boundary conditions. Matching inner and outer solutions of the single-particle probability density, we derive an asymptotic expansion of the Laplace transformed flux into each reactive surface in powers of ε, where ερ is a given target size. In turn, the fluxes determine the splitting probabilities for target absorption. We then extend our analysis to more general types of reactive targets by combining matched asymptotic analysis with an encounter-based formulation of diffusion-mediated surface reactions. That is, we derive an asymptotic expansion of the joint probability density for particle position and the so-called boundary local time, which characterizes the amount of time that a Brownian particle spends in the neighborhood of a point on a totally reflecting boundary. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. Robin boundary conditions are recovered in the special case of an exponential law for the stopping local times. Finally, we illustrate the theory by exploring how the leading-order contributions to the splitting probabilities depend on the choice of surface reactions. In particular, we show that there is an effective renormalization of the target radius of the form ρ→ρ-Ψ[over ̃](1/ρ), where Ψ[over ̃] is the Laplace transform of the stopping local time distribution.