[論文レビュー] 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

ひとこと要約

本稿では、自己強化とミタグ・レフラー分布の休息時間を持つ持続的ランダムウォークモデルを提案し、リーマン・リウヴィルの分数マージン導関数を用いた双曲型PDEで記述する。解析的および数値的に、自己強化により中間時間では一時的な超拡산が発生するが、長時間行動では重い尾を持つ休息時間のため劣拡散に戻ることを示しており、細胞内系における異常輸送の生物学的に妥当なメカニズムを提供する。

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.

研究の動機と目的

• 自己強化された方向性とべき乗則分布の休息時間を有する生物学的系における異常確率的輸送をモデル化すること。
• 自己強化と非マルコフ的休息状態の間の相互作用が、長時間輸送行動をどのように決定するかを調査すること。
• 超拡散から劣拡散への遷移を捉えることができる分数マージンPDE系を導出し、検証すること。
• 細胞小器官の移動における一時的超拡散の実験的観察と整合するメカニズムを提供すること。

提案手法

• アクティブ状態（±ν）と休息状態を有する三状態確率的モデルを定式化し、リーマン・リウヴィルの分数マージン導関数による非局所的遷移を伴う双曲型PDEで記述する。
• 時間依存の遷移確率 r± = 1/3 ± α₀x/(2νt) を用いて自己強化を導入する。ここで α₀ は自己強化パラメータである。
• リーマン・リウヴィルの分数マージン導関数 D₁₋ᵝₜ を用いて、ミタグ・レフラー分布の滞在時間を休息時間としてモデル化する。
• ラプラス空間における2番目のモーメント方程式を導出し、逆ラプラス変換を実行して長時間スケーリングを分析する。
• JITコンパイルされたPythonを用いた並列化によるモンテカルロシミュレーションを実施し、解析的予測の妥当性を検証する。
• 粒子位置の確率密度関数（PDF）を分析し、一時的な二峰性の歪みと長時間におけるラプラス型減衰を確認する。

実験結果

リサーチクエスチョン

• RQ1 自己強化は、ミタグ・レフラー分布の休息時間を有する持続的ランダムウォークにおいて超拡散を誘発するか？
• RQ2 重い尾を持つ休息時間の存在下で、2番目のモーメントの長時間漸近的挙動はいかなるものか？
• RQ3 非マルコフ的系に自己強化が存在する場合、一時的超拡散と最終的な劣拡散が共存しうるか？
• RQ4 べき乗則の休息分布の特徴的時間スケール τ₀ は、一時的超拡散の出現にどのように影響するか？
• RQ5 長時間極限において、この系はどの程度分数拡散方程式で近似可能か？

主な発見

• α₀ ≠ 0 のとき、中間時間では2番目のモーメントが µ₂(t) ∼ t³|α₀| とスケーリングし、一時的超拡散を示す。
• 長時間極限では、2番目のモーメントが µ₂(t) ∼ tᵝ とスケーリングし、0 < β < 1 となるため、ミタグ・レフラー分布の休息時間による劣拡散が確認される。
• 分数拡散極限は第2の例でのみ存在し、Dβ = ν²/(λ²Γ(β+1)τ₀ᵝ) を持つ分数拡散方程式を導く。
• τ₀ が小さい場合、中間時間におけるPDFは歪んでおり、二峰性を示し、自己強化による方向性の持続性を示す。
• τ₀ が大きい場合、長時間極限ではPDFはラプラス型に回帰し、劣拡散的挙動と整合的である。
• 非マルコフ的休息状態は、自己強化の強さ（α₀）にかかわらず、長時間極限において超拡散を完全に抑制する。