[论文解读] Continuous time random walk, Mittag-Leffler waiting time and fractional diffusion: mathematical aspects

Rudolf Gorenflo, Francesco Mainardi|ArXiv.org|May 6, 2007

Fractional Differential Equations Solutions参考文献 21被引用 11

一句话总结

本文建立了连续时间随机游走（CTRW）中幂律等待时间分布与Mittag-Leffler等待时间分布之间的渐近等价性，证明在适当的时间与空间重标度下，系统收敛于分数阶扩散过程。主要贡献在于严格推导出空间-时间分数阶扩散方程作为具有重尾等待时间的CTRW的标度极限，其中Mittag-Leffler函数作为基本核。

ABSTRACT

We show the asymptotic long-time equivalence of a generic power law waiting time distribution to the Mittag-Leffler waiting time distribution, characteristic for a time fractional CTRW. This asymptotic equivalence is effected by a combination of "rescaling" time and "respeeding" the relevant renewal process followed by a passage to a limit for which we need a suitable relation between the parameters of rescaling and respeeding. Turning our attention to spatially 1-D CTRWs with a generic power law jump distribution, "rescaling" space can be interpreted as a second kind of "respeeding" which then, again under a proper relation between the relevant parameters leads in the limit to the space-time fractional diffusion equation. Finally, we treat the `time fractional drift" process as a properly scaled limit of the counting number of a Mittag-Leffler renewal process.

研究动机与目标

• 在长时间极限下，建立幂律等待时间分布与Mittag-Leffler等待时间分布之间的数学等价性。
• 证明在时间与空间适当重标度下，具有幂律等待时间的连续时间随机游走收敛于时间分数阶CTRW。
• 将空间-时间分数阶扩散方程作为此类CTRW模型的标度极限进行推导。
• 通过拉普拉斯变换与傅里叶变换，形式化Mittag-Leffler函数与分数阶扩散过程之间的联系。
• 通过Gnedenko与Kovalenko提出的泊松过程薄化理论，将泊松过程推广为使用Mittag-Leffler分布的分数阶更新过程。

提出的方法

• 通过因子 $\tau$ 对时间进行重标度，并通过因子 $a$ 对更新过程进行重速度调整，通过 $\tau$ 与 $a$ 之间的特定关系，确保分布收敛。
• 利用拉普拉斯变换与傅里叶变换分析特征函数的收敛性，从而实现对极限过程的弱收敛性证明。
• 采用广义Mittag-Leffler函数 $E_{\beta,\beta}(-t^\beta)$ 作为等待时间密度，其等价于 $-\frac{d}{dt}E_\beta(-t^\beta)$，用于表征Mittag-Leffler更新过程。
• 通过因子 $h$ 进行空间重标度以建模空间中的加速，结合时间重标度，形成第二次重速度调整，最终导出空间-时间分数阶扩散方程。
• 应用Gnedenko与Kovalenko最初提出的更新过程薄化理论，以证明收敛至Mittag-Leffler等待时间分布。
• 使用Riesz-Feller分数阶导数表示空间分数阶扩散算子，其符号由严格稳定分布的特征函数导出。

实验结果

研究问题

• RQ1 如何通过重标度与重速度调整，使幂律等待时间分布渐近等价于Mittag-Leffler等待时间分布？
• RQ2 当时间与空间均经过适当重标度时，具有幂律等待时间的连续时间随机游走的极限行为是什么？
• RQ3 时间与空间重标度的结合如何导致空间-时间分数阶扩散方程的出现？
• RQ4 Mittag-Leffler函数在表征分数阶更新过程的等待时间分布中起什么作用？
• RQ5 时间分数阶漂移过程能否作为Mittag-Leffler更新过程的适当标度极限导出？

主要发现

• 在适当的重标度与重速度调整下，Mittag-Leffler等待时间密度 $\phi^{\text{ML}}(t) = -\frac{d}{dt}E_\beta(-t^\beta)$ 成为幂律等待时间分布的普遍极限，其拉普拉斯变换为 $\widetilde{\phi}(s) = \frac{1}{1+s^\beta}$。
• 生存函数 $\Psi(t) = E_\beta(-t^\beta)$ 在 $0 < \beta < 1$ 时完全单调，保持了Mittag-Leffler过程作为泊松过程分数阶推广的概率解释。
• 当 $t$ 较大时，$\phi(t)$ 的渐近行为为 $\sim \frac{\sin(\beta\pi)}{\pi} \frac{\Gamma(\beta+1)}{t^{\beta+1}}$，证实其具有重尾特性。
• 在时间重标度 $\tau$ 与空间重标度 $h$ 之间存在特定关系时，具有幂律等待时间与幂律跳跃的CTRW的标度极限即为包含Riesz-Feller空间分数阶导数的分数阶扩散方程。
• Mittag-Leffler更新过程的计数函数 $N(t)$ 在分布上收敛于时间分数阶漂移过程，后者由分数阶Fokker-Planck方程表征。
• Mittag-Leffler函数 $E_\beta(z)$ 及其导数提供了生存函数与等待时间密度的精确表达式，可通过拉普拉斯变换与谱测度给出积分表示。