Uniform asymptotic approximation of diffusion to a small target: generalized reaction models
arXiv:1605.01279 · doi:10.1103/PhysRevE.94.042414
Abstract
The diffusion of a reactant to a binding target plays a key role in many biological processes. The reaction-radius at which the reactant and target may interact is often a small parameter relative to the diameter of the domain in which the reactant diffuses. We develop uniform in time asymptotic expansions in the reaction-radius of the full solution to the corresponding diffusion equations for two separate reactant-target interaction mechanisms: the Doi or volume reactivity model, and the Smoluchowski-Collins-Kimball partial absorption surface reactivity model. In the former, the reactant and target react with a fixed probability per unit time when within a specified separation. In the latter, upon reaching a fixed separation, they probabilistically react or the reactant reflects away from the target. Expansions of the solution to each model are constructed by projecting out the contribution of the first eigenvalue and eigenfunction to the solution of the diffusion equation, and then developing matched asymptotic expansions in Laplace-transform space. Our approach offers an equivalent, but alternative, method to the pseudo-potential approach we previously employed in [Isaacson and Newby 2013] for the simpler Smoluchowski pure absorption reaction mechanism. We find that the resulting asymptotic expansions of the diffusion equation solutions are identical with the exception of one parameter: the diffusion limited reaction rates of the Doi and partial absorption models. This demonstrates that for biological systems in which the reaction-radius is a small parameter, properly calibrated Doi and partial absorption models may be functionally equivalent.
References in corpus (7)
- First-passage times in complex scale-invariant media
- Mean first-passage times of non-Markovian random walkers in confinement
- Uniform asymptotic approximation of diffusion to a small target
- A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems
- A First-Passage Kinetic Monte Carlo Method for Reaction-Drift-Diffusion Processes
- The Influence of Spatial Variation in Chromatin Density Determined by X-ray Tomograms on the Time to Find DNA Binding Sites
- Reactive Boundary Conditions as Limits of Interaction Potentials for Brownian and Langevin Dynamics