[论文解读] Random walks and Brownian motion: A method of computation for first-passage times and related quantities in confined geometries

S. Condamin, Olivier Bénichou+1 more|arXiv (Cornell University)|Oct 9, 2006

Diffusion and Search Dynamics被引用 5

一句话总结

本文提出了一套统一的计算框架，用于在受限几何结构中计算离散随机游走和连续布朗运动的平均首达时间（MFPT）、分裂概率以及完整的首达时间分布。通过利用拉普拉斯算子在适当边界条件下的伪格林函数，该方法能够在任意有界区域中实现MFPT的精确解析解与近似解，特别针对球形、矩形及圆盘状几何结构给出了显式结果，显著拓展了以往在复杂环境中首达过程研究的成果。

ABSTRACT

In this paper we present a computation of the mean first-passage times both for a random walk in a discrete bounded lattice, between a starting site and a target site, and for a Brownian motion in a bounded domain, where the target is a sphere. In both cases, we also discuss the case of two targets, including splitting probabilities, and conditional mean first-passage times. In addition, we study the higher-order moments and the full distribution of the first-passage time. These results significantly extend our earlier contribution [Phys. Rev. Lett. 95, 260601].

研究动机与目标

• 开发一种适用于离散随机游走和连续布朗运动的有界区域中首达时间的通用计算方法。
• 通过计算不仅包括平均首达时间（MFPT），还涵盖高阶矩和完整的首达时间分布，以扩展先前的研究成果。
• 分析几何结构（特别是目标与源相对于边界的相对位置）对MFPT的影响，包括存在两个竞争目标的情况。
• 为对称区域（如球体、矩形）提供MFPT的显式解析表达式，并为一般形状提供实用的近似方法。
• 实现对真实生物和物理系统中首达时间的定量估计，例如扩散控制的反应以及细胞内搜索过程。

提出的方法

• 该方法基于Kac和Aldous-Fill针对不可约马尔可夫链的公式，利用平稳分布和转移概率在离散格点上计算MFPT。
• 对于连续布朗运动，该方法采用拉普拉斯算子的伪格林函数，并在目标处施加吸收边界条件，在区域边界施加反射条件。
• MFPT表示为伪格林函数的泛函形式，该函数在对称区域（如球体、矩形、圆盘）中可解析求解，其他情况下则进行近似。
• 在存在两个目标的情况下，使用相同的格林函数形式计算分裂概率和条件MFPT。
• 通过镜像法（如反射原理）将边界效应纳入考虑，以处理靠近平面或曲面边界的情况。
• 对于非对称区域，提出了伪格林函数的初步近似方法，经数值检验和物理一致性验证。

实验结果

研究问题

• RQ1 随机游走者的平均首达时间（MFPT）如何依赖于有界区域的几何形状，特别是当目标靠近边界时？
• RQ2 在受限几何结构中，首达时间的完整概率分布是什么？高阶矩（如方差）的行为如何？
• RQ3 当有兩個競爭性吸收目标存在于有界区域中时，分裂概率和条件MFPT的行为如何？
• RQ4 在二维和三维空间中，源与目标相对于边界的相对位置对MFPT的定量影响是什么？
• RQ5 如何计算或近似一般有界区域的拉普拉斯算子的伪格林函数，以实现首达时间的计算？

主要发现

• 在三维球形区域中，MFPT通过伪格林函数显式计算得出，边界接近引起的修正项与目标到边界的距离成反比。
• 在二维情况下，MFPT对源-目标距离呈对数依赖关系，边界效应更强，当目标靠近平面边界时，MFPT可能翻倍。
• 当源靠近边界时，MFPT降低，但该效应弱于目标靠近边界时引起的增加，尤其是在二维情况下。
• 在三维平面边界情况下，MFPT的修正项近似为 (1/(4π)) × (1/a − 1/|rS−rT| + 1/|rT−s(rT)| − 1/|s(rS)−rT|)，其中 s(r) 表示镜像点。
• 在二维情况下，边界修正项为 (1/(2π)) × (ln(|rS−rT|/a) + ln(|s(rS)−rT|/|s(rT)−rT|))，显示出对几何结构的对数依赖性。
• 该方法在离散和连续设置下均能对MFPT、分裂概率及完整的首达时间分布提供高精度估计，并为球形和矩形区域提供了显式公式。