[论文解读] Anomalous Stochastic Transport of Particles with Self-Reinforcement and Mittag–Leffler Distributed Rest Times

Daniel Han|arXiv (Cornell University)|Nov 15, 2021

Diffusion and Search Dynamics参考文献 65被引用 11

一句话总结

本文提出了一种具有自增强机制和Mittag-Leffler分布休息时间的持续随机游走模型，通过带有Riemann-Liouville分数阶导数的双曲型偏微分方程进行建模。研究通过解析和数值方法证明，尽管由于自增强作用在中间时间尺度出现瞬时超扩散，但长期行为由于重尾休息时间而恢复为亚扩散，从而为细胞内系统中异常输运提供了生物上合理的机制。

ABSTRACT

We introduce a persistent random walk model for the stochastic transport of particles involving self-reinforcement and a rest state with Mittag–Leffler distributed residence times. The model involves a system of hyperbolic partial differential equations with a non-local switching term described by the Riemann–Liouville derivative. From Monte Carlo simulations, we found that this model generates superdiffusion at intermediate times but reverts to subdiffusion in the long time asymptotic limit. To confirm this result, we derived the equation for the second moment and find that it is subdiffusive in the long time limit. Analyses of two simpler models are also included, which demonstrate the dominance of the Mittag–Leffler rest state leading to subdiffusion. The observation that transient superdiffusion occurs in an eventually subdiffusive system is a useful feature for applications in stochastic biological transport.

研究动机与目标

• 为具有自增强方向性和幂律分布休息时间的生物系统中的异常随机输运建模。
• 研究自增强与非马氏休息状态之间的相互作用，以确定长期输运行为。
• 推导并验证一个捕捉从超扩散到亚扩散转变的分数阶偏微分方程系统。
• 提供一种与细胞器运动中瞬时超扩散的实证观测一致的机制。

提出的方法

• 构建一个包含两个活跃状态（±ν）和一个休息状态的三状态随机模型，其演化由带有非局部切换机制的双曲型偏微分方程控制，切换过程通过Riemann-Liouville分数阶导数描述。
• 通过时变转移概率 r± = 1/3 ± α₀x/(2νt) 引入自增强机制，其中 α₀ 为自增强参数。
• 利用 Riemann-Liouville 分数阶导数 D₁₋ᵝₜ 对 Mittag-Leffler 分布的驻留时间进行建模。
• 在拉普拉斯空间中推导二阶矩方程，并执行逆拉普拉斯变换以分析长期标度行为。
• 使用 JIT 编译的 Python 并行化实现蒙特卡洛模拟，以验证解析预测结果。
• 分析粒子位置的概率密度函数（PDF），以确认存在瞬时双峰偏态及长期拉普拉斯衰减。

实验结果

研究问题

• RQ1 在具有 Mittag-Leffler 分布休息时间的持续随机游走中，自增强是否会引起超扩散？
• RQ2 在重尾休息时间存在的情况下，二阶矩的长期渐近行为如何？
• RQ3 在具有自增强的非马氏系统中，瞬时超扩散能否与最终的亚扩散共存？
• RQ4 幂律休息分布的时间尺度 τ₀ 如何影响瞬时超扩散的出现？
• RQ5 在长时间极限下，该系统在多大程度上可近似为分数阶扩散方程？

主要发现

• 当 α₀ ≠ 0 时，在中间时间尺度，二阶矩呈现 µ₂(t) ∼ t³|α₀| 的标度行为，表明存在瞬时超扩散。
• 在长时间极限下，二阶矩呈现 µ₂(t) ∼ tᵝ 的标度行为，其中 0 < β < 1，证实了由于 Mittag-Leffler 分布休息时间导致的亚扩散。
• 仅在第二个示例中存在分数阶扩散极限，其结果为具有 Dβ = ν²/(λ²Γ(β+1)τ₀ᵝ) 的分数阶扩散方程。
• 当 τ₀ 较小时，PDF 在中间时间呈现偏斜且双峰，表明自增强导致的方向持久性。
• 当 τ₀ 较大时，PDF 在长时间极限下恢复为拉普拉斯形式，与亚扩散行为一致。
• 非马氏休息状态在长时间极限下完全抑制了超扩散，无论自增强强度（α₀）如何。