Fractional kinetics equation from a Markovian system of interacting Bouchaud trap models

技術要約：マルコフ的な相互作用ブシャード・トラップ・モデルから導かれる分数階キネティクス方程式

問題設定
本論文は、ランダムな環境下で進化する相互作用粒子系のマクロな流体力学的極限（HL）の導出に取り組んでいる。ランダムなコンダクタンスにおける単純対称排斥過程（SSEP）の流体力学的極限は確立されているが（熱方程式または準拡散をもたらす）、時間分数階劣拡散方程式（具体的には分数階キネティクス方程式、FKE）へとスケーリングされる厳密な微視的モデルは不足している。著者らは、重い裾を持つランダムなトラッピング環境において、適切な時空スケーリングによって、$d \ge 2$ では決定論的な測量（決定論的な測度）を満たす分数階偏微分方程式へ、$d=1$ ではランダムな準拡散へと収束する、マルコフ的な相互作用粒子系を構築することを目的としている。

モデルの定義
著者らは、格子 $\mathbb{Z}^d$ 上の**相互作用ブシャード・トラップ・モデル（BTM）**系を導入する。

• 環境: 環境は、重い裾を持つ分布に従う独立同一分布（i.i.d.）のランダム変数 $\alpha = (\alpha_x)_{x \in \mathbb{Z}^d}$ で定義される。ここで、$Q(\alpha_0 > u) \sim u^{-\beta}$ であり、$\beta \in (0, 1)$ である。
• 粒子の動力学: この系は、部分的な排斥規則に従って $\mathbb{Z}^d$ 上で進化する粒子からなる。各サイト $x$ における最大占有数は $\alpha_x$ である。サイト $x$ から隣接サイト $y$ へのジャンプ率は、$\alpha_x^{a-1} \alpha_y^a (1 - \eta(y)/\alpha_y)$ に比例する。ここで $\eta(y)$ はサイト $y$ における粒子数であり、$a \in [0, 1]$ は非対称性を制御するパラメータである。
  • 項 $(1 - \eta(y)/\alpha_y)$ は排斥規則を強制する。すなわち、占有されたサイトへのジャンプは抑制され、サイトの「容量」はトラップサイズ $\alpha_x$ によって決定される。
  • $a=0$ のとき、モデルは対称的であり、$a \in (0, 1]$ のとき、非対称的となる。
• 主要な特徴: 独立なBTMとは異なり、これらの粒子は排斥規則を通じて相互作用するが、系は確率的自己双対性という性質を保持しており、これが解析において極めて重要となる。

手法
流体力学的極限の証明は、相互作用系の確率的自己双対性に依拠している。これにより、従来のマルチンゲール問題のアプローチ（およびそのための1ブロック・エスティメイトや2ブロック・エスティメイト）を回避している。著者らは、これら従来の手法を時間分数階の極限へ拡張することは困難であると指摘している。

1. 双対関係: 著者らは、相互作用粒子系と、独立なBTM（および相互作用するBTMのペア）との間の双対関係を利用する。これにより、経験場の期待値を、単一粒子のスケーリング極限の研究へと還元できる。
2. 経験場: 2つの異なる経験場が解析される：
   • 経験密度場 ($X^n_t$): $n^{-d/\beta}$ でスケーリングされる。
   • 経験頻度場 ($Z^n_t$): $n^{-d}$ でスケーリングされ、粒子を局所的な容量 $\alpha_x$（すなわち $\eta(x)/\alpha_x$）で重み付けする。
3. 単一粒子のスケーリング極限: 解析は、既知の単一粒子BTMのスケーリング極限を活用する：
   • $d \ge 2$ の場合、リスケールされたBTMは分数階キネティクス過程（時間変更されたブラウン運動）に収束する。
   • $d = 1$ の場合、それはFIN拡散（特異な速度測度を持つ時間変更されたブラウン運動）に収束する。
4. 分解: 証明は、経験場を、動的ノイズ、初期条件の分散、初期条件の平均、および決定論的極限項の4つの項に分解する。著者らは、ノイズ項と分散項が消失すること、そして平均項が目標とするPDEの解に収束することを示す。

主要な結果

本論文は、特定の時間スケーリング $\theta_n$ の下での経験頻度場 $Z^n_t$ および密度場 $X^n_t$ の流体力学的極限を確立している：

• ケース $d \ge 2$:

  • 時間スケーリング: $\theta_n \sim n^{2/\beta}$ （$d=2$ では対数補正を伴う）。
  • 頻度場 ($Z^n_t$): 密度 $\rho_t$ が以下の分数階キネティクス方程式を解く決定論的な測量へと収束する：
    $$\frac{\partial^\beta \rho_t}{\partial t^\beta} = D_{\text{eff}} \Delta \rho_t$$
    ここで $\partial^\beta/\partial t^\beta$ はカプート（Caputo）型分数階微分である。環境のランダム性は均質化され、決定論的なマクロ極限が生じる。
  • 密度場 ($X^n_t$): $\rho_0 W$（ここで $W$ はポアソン点過程上に支持されるランダムな測量）というランダムな測量へと収束する。マクロな動力学は、環境のランダム性が極限においても残存するという意味で「凍結」されている。

• ケース $d = 1$:

  • 時間スケーリング: $\theta_n \sim n^{1+1/\beta}$。
  • 頻度場 ($Z^n_t$): 密度 $\rho_t$ がランダムな準拡散方程式を解むランダムな測量へと収束する：
    $$\frac{\partial \rho_t}{\partial t} = \frac{\partial}{\partial W} \frac{\partial \rho_t}{\partial x}$$
    ここで、微分はランダムな測量 $W$（FIN拡散の速度測量）の意味において解釈される。
  • 密度場 ($X^n_t$): これもまた、ランダムな測量 $\rho_t W$ へと収束する。

意義と主張
著者らは、本研究が、時間分数階の劣拡散方程式へとスケーリングされる最初の例（相互作用するBTM）を提供すると主張している。

• 手法の新規性: 本論文は、トラップモデルに関する先行研究がマルチンゲール問題に依存していたことを強調している。それらの手法は、非マルコフ的（分数階）な極限を扱うのが困難である。確率的双対性を利用することで、著者らは、前段階の相互作用系と流体力学的極限との間の直接的な接続を確立し、時間分数階の設定を成功裏に扱うことができた。
• 均質化 vs ランダム性: 2つの経験場の間の区別は重要な発見である。$d \ge 2$ では、頻度場は均質化を示し（ランダム性が消失して決定論的な分数階偏微分方程式となる）、一方で密度場は環境のランダム性を保持する。これは、$d=1$ のケース（両方の場がランダムなままとなる）とは対照的である。
• 分数階キネティクスの微視的起源: 本研究は、純粋にマルコフ的な相互作用粒子系の微視的システムが、重い裾を持つトラッピング時間と部分的な排斥のメカニズムを通じて、いかにしてマクロな非マルコフ過程（分数階キネティクス過程）を生み出し得るかを実証している。

結論として、本論文は、導出された分数階キネティクス方程式が、特定の微視的相互作用動力学の厳密な流体力学的極限であることを示し、微視的なトラップモデルとマクロな分数階偏微分方程式の間の溝を埋めている。