[论文解读] An encounter-based approach for restricted diffusion with a gradient drift

Denis S. Grebenkov|arXiv (Cornell University)|Oct 23, 2021

Diffusion and Search Dynamics参考文献 74被引用 47

一句话总结

本文提出了一种基于相遇的受限扩散方法，结合梯度漂移，通过将狄利克雷-诺伊曼算子扩展以推导全传播器的谱分解——即粒子位置与边界局部时间的联合概率密度。该方法可精确计算反应动力学、存活概率及首次通过时间统计，揭示了漂移如何增强边界相遇并改变受限几何结构中局部时间的统计特性。

ABSTRACT

We develop an encounter-based approach for describing restricted diffusion with a gradient drift towards a partially reactive boundary. For this purpose, we introduce an extension of the Dirichlet-to-Neumann operator and use its eigenbasis to derive a spectral decomposition for the full propagator, i.e., the joint probability density function for the particle position and its boundary local time. This is the central quantity that determines various characteristics of diffusion-influenced reactions such as conventional propagators, survival probability, first-passage time distribution, boundary local time distribution, and reaction rate. As an illustration, we investigate the impact of a constant drift onto the boundary local time for restricted diffusion on an interval. More generally, this approach accesses how external forces may influence the statistics of encounters of a diffusing particle with the reactive boundary.

研究动机与目标

• 开发一个系统化的框架，用于建模朝向部分反应边界具有梯度漂移的受限扩散。
• 将边界局部时间作为核心随机变量，以描述粒子与反应表面的相遇。
• 将狄利克雷-诺伊曼算子扩展，以实现漂移存在下全传播器的谱分解。
• 实现对关键反应-扩散量（如存活概率、首次通过时间与反应速率）的精确计算。
• 分析外力（漂移）如何定量影响受限几何结构中边界相遇的统计特性。

提出的方法

• 引入一个扩展的狄利克雷-诺伊曼算子，以同时包含位置与边界局部时间的动力学。
• 利用扩展算子的本征函数，推导全传播器的谱分解。
• 将位置与边界局部时间的联合概率密度作为所有反应-扩散特性的核心量。
• 将谱方法应用于具有常数漂移的罗宾型边界条件下的福克-普朗克方程。
• 通过一维区间上的随机模拟与解析解验证该方法。
• 通过将边界局部时间作为动力学变量，统一斯科罗霍德随机微分方程与福克-普朗克形式化。

实验结果

研究问题

• RQ1 常数漂移如何影响区间上受限扩散中边界局部时间的分布？
• RQ2 能否通过扩展的狄利克雷-诺伊曼算子对带漂移的扩散全传播器实现谱分解？
• RQ3 外力如何改变粒子与部分反应边界相遇的统计特性？
• RQ4 漂移对受限几何结构中存活概率与首次通过时间分布有何影响？
• RQ5 基于相遇的框架能否推广至罗宾边界条件之外，以建模与相遇相关的反应性？

主要发现

• 通过扩展的狄利克雷-诺伊曼算子的本征函数，推导出全传播器的谱分解，从而实现对位置-局部时间联合分布的精确求解。
• 常数漂移显著改变边界局部时间的分布：正漂移（朝向边界）增加局部时间，负漂移（远离）则减少之，且在区间上导出了显式公式。
• 存活概率与首次通过时间分布被精确计算，并显示其对漂移强度与边界反应性的非平凡依赖关系。
• 该方法揭示，漂移可增强边界相遇的频率与持续时间，其量化通过边界局部时间的均值与方差实现。
• 该方法可推广至罗宾边界条件之外，支持对与相遇相关的反应性及非马尔可夫表面反应的建模。
• 在强漂移极限下，系统从扩散主导过渡为漂移主导动力学，边界局部时间呈现ℓt ∼ t的标度行为（正漂移），而强负漂移下则趋近于零。