[论文解读] Mean first-passage time to a small absorbing target in three-dimensional elongated domains

Denis S. Grebenkov, Alexei T. Skvortsov|arXiv (Cornell University)|Feb 15, 2022

Diffusion and Search Dynamics参考文献 68被引用 12

一句话总结

本文推导出在横截面缓慢变化的三维轴对称长形域中，对任意形状小吸收目标的平均首达时间（MFPT）的近似解析公式。通过将三维福克-普朗克问题简化为在目标位置具有半透膜的1D有效扩散方程，该方法通过基于静电电容的等效反应性，捕捉了MFPT对目标位置、径向分布及目标形状的依赖关系，并通过蒙特卡洛模拟得到验证。

ABSTRACT

We derive an approximate formula for the mean first-passage time (MFPT) to a small absorbing target of arbitrary shape inside an elongated domain of a slowly varying axisymmetric profile. For this purpose, the original Poisson equation in three dimensions is reduced an effective one-dimensional problem on an interval with a semi-permeable semi-absorbing membrane. The approximate formula captures correctly the dependence of the MFPT on the distance to the target, the radial profile of the domain, and the size and the shape of the target. This approximation is validated by Monte Carlo simulations.

研究动机与目标

• 为具有任意轴对称截面的三维长形域中，小吸收目标的平均首达时间（MFPT）发展一种通用近似方法。
• 解决窄逃逸理论在高度各向异性域中的局限性，其中目标在横向尺寸小但沿长度方向不可忽略。
• 通过基于静电电容推导的等效捕获系数，将目标形状与径向位置对MFPT的影响纳入模型。
• 通过在多种域几何形状（圆柱形、锥形、振荡形）中与蒙特卡洛模拟对比，验证该方法的有效性。
• 通过考虑径向限制与熵漂移效应，将二维平面方法推广至三维。

提出的方法

• 将任意形状的小吸收目标替换为具有相同捕获系数K的等效吸收圆盘，K由其静电电容C通过K = 4πDC推导得出。
• 在目标纵向位置zT处引入一个具有有效反应性κ = K(rT)/S(zT)的半透半吸收膜，其中S(z)为横截面积。
• 将MFPT建模为受后向福克-普朗克方程控制的1D问题：d/dz [S(z) dT/dz] = -S(z)/D，边界条件为z=0和z=ℓ处的诺伊曼条件。
• 基于轴向与壁面相邻极限之间的插值，对位置相关的捕获系数K(rT)建立关于目标尺寸ν = a/r(zT)与径向偏移η = rT/r(zT)的解析近似。
• 在两个区间（0,zT）与（zT,ℓ）中分别求解1D MFPT方程，并在zT处施加匹配条件，以获得完整的MFPT分布。
• 通过在具有不同截面形状（圆柱形、锥形、振荡形）的域中进行大量蒙特卡洛模拟，验证解析结果。

实验结果

研究问题

• RQ1 在三维长形域中，小目标的MFPT如何依赖于其径向位置与形状，而不仅限于窄逃逸极限？
• RQ2 有效1D膜模型能否准确捕捉具有任意目标几何形状的各向异性、缓慢变化的3D域中的MFPT？
• RQ3 用于替代小目标的膜的正确有效反应性κ是什么？其与目标电容及位置的关系如何？
• RQ4 径向限制与熵漂移如何影响长形域中的MFPT？该效应能否在简化的1D模型中被准确捕捉？
• RQ5 该近似方法在非球形目标与复杂域截面下的适用范围有多大？

主要发现

• MFPT对目标的径向位置rT有显著依赖，这是由于位置相关的捕获系数K(rT)考虑了约束域曲率的影响。
• 有效反应性κ = K(rT)/S(zT)在1D模型中成功捕捉了目标的吸收效率，其中K(rT)由插值轴向与壁面相邻极限的校准函数Ψ(ν,η)推导得出。
• 通过Fick-Jacobs型方程d/dz [S(z) dT/dz] = -S(z)/D实现的1D降维方法，能准确再现具有缓慢变化截面的长形域中的MFPT。
• 该方法在多种域几何形状（圆柱形、锥形、振荡形）中均表现出高精度，与蒙特卡洛模拟结果一致。
• 即使目标在横向尺寸小但沿长度方向不可忽略时，该近似仍保持有效，克服了经典窄逃逸理论的局限性。
• 该方法可实现对涉及管状或丝状结构中扩散搜索的复杂生物与物理系统，进行快速、解析的MFPT估算。