[論文レビュー] Anomalous diffusion: Fractional Brownian motion vs fractional Ito motion

Iddo Eliazar, Tal Kachman|arXiv (Cornell University)|Jan 1, 2022

Fractional Differential Equations Solutions参考文献 133被引用数 14

ひとこと要約

本稿では、異常拡散のモデリングのための分数カリオ運動（FIM）を導入し、分数ブラウン運動（FBM）の扱いやすい代替手段として提示する。FBMとは異なり、FIMはマルコフ的かつマルティングールである。また、非ガウス的かつ非定常な速度を有するため、シミュレーションが容易で解析的扱いも可能である。主な貢献は、FIMが対数ポテンシャル内での拡散と正確に等価であることにあり、これは、べき乗則に従う平均二乗変位を示す、非ガウス的かつ解析的に解けるサブ拡散およびスーパー拡散のモデルを提供する。

ABSTRACT

Contains fulltext : 247529.pdf (Publisher’s version ) (Open Access)

研究の動機と目的

• 分数ブラウン運動（FBM）の限界を解決すること。FBMは、広く異常拡散モデリングに用いられているが、非マルコフ的かつ非マルティングール的であり、解析的に扱いにくい。
• 分数カリオ運動（FIM）を、自己相似的かつ連続的軌道を持つプロセスとして提案すること。FIMは、FBMと同様に、ハースト指数に依存する拡散挙動（サブ拡散、通常拡散、スーパー拡散）を維持する。
• FIMがマルコフ過程かつマルティングールであることを示し、FBMの複雑さとは対照的に、解析的取り扱いと効率的なシミュレーションが可能であることを明らかにすること。
• FIMと対数ポテンシャル内での確率的力学との間の厳密な数学的関係を確立し、プロセスに物理的解釈を与えること。
• FIMとFBMを、非ガウス性、速度の定常性、有限時間分布のカルバック・ライブラー発散といった、主な統計的および動的特性において比較すること。

提案手法

• 逆ガウス分布族を用いてFIMの確率密度関数（PDF）を導出し、それが非ガウス的で、重い尾を持つことを示す。
• 自己相似性を用いて、FIMを法則的に IH(t) = t^H · IH(1) として表現し、有限時間分布のスケーリングを可能にする。
• 伊藤積分を用いて、FIMの確率微分方程式（SDE）を導出する：dIH(t) = |IH(t)|^{1 - 1/(2H)} dB(t)。これはべき乗則に従うボラティリティを持つ。
• 非線形写像 ϕ(x) = 2H |x|^{1/(2H)} sign(x) を用いてFIMを変換し、等価なランジュバンSDE dξH(t)/dt = (1/2 - H)/ξH(t) + dB(t) を得る。これにより、FIMは対数ポテンシャル内での拡散と関連づけられる。
• FBMおよびFIMの有限時間分布と初期時間分布との間のカルバック・ライブラー（KL）発散を計算し、異なるスケーリング行動を明らかにする。
• FIMの増分の分散を分析し、長時間において Var[ΔIH] ∝ t^{2H-1} と成り立つことを示し、非定常な速度および非マルコフ的スケーリングに類似した振る舞いを示す。

実験結果

リサーチクエスチョン

• RQ1 分数カリオ運動（FIM）の統計的および動的挙動は、分数ブラウン運動（FBM）と比べて、マルコフ性、マルティングール性、ガウス性の観点でどのように異なるか？
• RQ2 FIMの正確な確率的力学は何か？ また、対数ポテンシャルを持つランジュバン方程式にどのように変換できるか？
• RQ3 FIMとFBMのカルバック・ライブラー発散は時間経過とともにどのように変化するか？ これにより、プロセス内の情報損失または増加はどのように明らかになるか？
• RQ4 FIMの増分分散の挙動は何か？ また、これは速度プロセスの非定常性をどのように反映するか？
• RQ5 FIMは、束縛ポテンシャル内での物理的拡散プロセスとして解釈可能か？ もしそうなら、そのポテンシャルの性質は何か？

主な発見

• FIMはマルコフ過程かつマルティングールである。これに対してFBMは、両方とも満たさない。これにより、マルティングール理論からの強力な解析的ツールが利用可能になる。
• FIMはガウス過程ではない。時間tにおける周辺分布は非ガウス的であり、|x|^{1/H - 2} のように大きな|x|で減衰するPDFを持つ。
• FIMの時間tにおける分布と時間1における分布とのカルバック・ライブラー発散は、漸近的に (1 - H)(t - 1 - ln t) とスケーリングし、FBMとは異なり、対数的でない発散行動を示すことが証明された。
• FIMの増分分散は、大規模なtにおいて Var[ΔIH] ≈ v₁Δ²ᴴ t^{2H-1} とスケーリングし、非定常な速度および時間依存の拡散係数を示している。
• FIMは数学的に、V(x) = (1/2 - H) ln|x| の対数ポテンシャル内を拡散する粒子と正確に等価であり、SDE dξH/dt = (1/2 - H)/ξH + dB(t) で記述される。
• FIMプロセスは、FBMと同様にべき乗則に従う平均二乗変位 φ(t) = c t^{2H} を示し、サブ拡散、通常拡散、スーパー拡散の各領域において、異常拡散のモデリングとしての有効性が確認された。