[论文解读] A First Look at First-Passage Processes

S. Redner|arXiv (Cornell University)|Jan 25, 2022

Diffusion and Search Dynamics被引用 6

一句话总结

本文全面介绍了随机动力学中首达过程的理论，通过生成函数推导了首达概率与占据概率之间的基本联系。研究建立了受限几何结构（如半直线、有限区间和吸收楔形）中扩散过程的关键结果，并将其应用于物理和生物系统，包括反应速率、神经元放电和种群灭绝，揭示了普遍存在的标度律，例如在生灭过程中存活概率按 $ t^{-1} $ 衰减，而在扩散过程中按 $ t^{-1/2} $ 衰减。

ABSTRACT

These notes are based on the lectures that I gave (virtually) at the Bruneck Summer School in 2021 on first-passage processes and some applications of the basic theory. I begin by defining what is a first-passage process and presenting the connection between the first-passage probability and the familiar occupation probability. Some basic features of first passage on the semi-infinite line and a finite interval are then discussed, such as splitting probabilities and first-passage times. I also treat the fundamental connection between first passage and electrostatics. A number of applications of first-passage processes are then presented, including the hitting probability for a sphere in greater than two dimensions, reaction rate theory and its extension to receptors on a cell surface, first-passage inside an infinite absorbing wedge in two dimensions, stochastic hunting processes in one dimension, the survival of a diffusing particle in an expanding interval, and finally the dynamics of the classic birth-death process.

研究动机与目标

• 建立利用生成函数将首达概率与占据概率联系起来的理论框架。
• 分析受限几何结构（如半直线和有限区间）中的首达现象，重点关注分裂概率与首达时间。
• 探讨首达过程中的静电类比，从而为复杂边界条件提供解析解。
• 将首达理论应用于真实世界系统，包括细胞表面的反应动力学、随机捕食行为以及种群动力学。
• 展示不同随机过程中首达标度律的普适性，例如生灭过程中的 $ S(t) \sim t^{-1} $ 与扩散过程中的 $ S(t) \sim t^{-1/2} $ 的对比。

提出的方法

• 推导占据概率 $ P(\mathbf{r},t) $ 与首达概率 $ F(\mathbf{r},t) $ 之间的卷积关系，表达式为 $ P(\mathbf{r},t) = \delta_{\mathbf{r},0}\delta_{t,0} + \sum_{t' \leq t} F(\mathbf{r},t') P(0,t-t') $。
• 使用生成函数 $ P(\mathbf{r},z) $ 和 $ F(\mathbf{r},z) $ 将卷积关系转化为代数关系：对于 $ \mathbf{r} \neq 0 $，有 $ F(\mathbf{r},z) = P(\mathbf{r},z)/P(0,z) $。
• 对带有吸收边界条件的扩散方程应用拉普拉斯变换，以计算连续空间与时间中的首达时间分布。
• 利用静电类比，通过格林函数和吸收边界处的通量，将首达问题映射到势论。
• 通过生成函数 $ g(z,t) = \sum_n P_n(t) z^n $ 求解生灭过程的主方程，将其转化为一阶偏微分方程：$ g_t = (1-z)^2 g_z $。
• 引入变量替换 $ y = 1/(1-z) $，将该偏微分方程化为波动方程 $ g_t = g_y $，并利用特征线法求解。

实验结果

研究问题

• RQ1 在随机游走中，首达概率与占据概率之间的确切关系是什么？
• RQ2 在半无限直线上和有限区间上，首达时间与分裂概率的行为如何？
• RQ3 首达过程与静电学之间有何联系？该联系如何用于求解边界值问题？
• RQ4 当粒子被限制在扩展区间中时，其存活概率如何随时间衰减？
• RQ5 生灭过程的首达特性是什么？其与简单扩散的首达特性相比有何异同？

主要发现

• 对于 $ \mathbf{r} \neq 0 $，首达概率 $ F(\mathbf{r},z) $ 完全由占据概率 $ P(\mathbf{r},z) $ 通过关系式 $ F(\mathbf{r},z) = P(\mathbf{r},z)/P(0,z) $ 确定。
• 在半直线上，通过拉普拉斯变换求解在 $ x=0 $ 处具有吸收边界条件的扩散方程，推导出从 $ x_0 > 0 $ 出发的粒子的首达时间分布。
• 对于在扩展区间中扩散的粒子，其存活概率按 $ t^{-1/2} $ 衰减，与标准扩散类似，但边界条件随时间变化。
• 在线性出生率的生灭过程中，存活概率 $ S(t) = 1 - P_0(t) $ 按 $ 1/(1+t) $ 衰减，表明尽管灭绝是必然的，但平均灭绝时间却为无穷大。
• 生灭过程的生成函数解给出了精确表达式：$ P_0(t) = t/(1+t) $，$ P_n(t) = t^{n-1}/(1+t)^{n+1} $（$ n \geq 1 $）。
• 静电类比使得首达问题可映射到势论，从而能够对复杂几何结构（如吸收楔形和球形受体）给出精确解。