Numerical analysis of first-order mean field games under displacement monotonicity

技術要約：変位単調性（Displacement Monotonicity）の下での一次平均場ゲームの数値解析

問題設定
本研究は、時間依存の一次平均場ゲーム（MFG）系の厳密な数値近似および解析に取り組むものである。具体的な焦点は、非分離型のハミルトニアンを持ち、かつ $P_2(\mathbb{R}^d)$（有限の二次モーメントを持つ確率測度の空間）における潜在的に特異な初期プレイヤー分布を特徴とするシステムである。連続的なMFG系は、値関数 $u$ に対するハミルトン–ヤコビ–ベルマン（HJB）方程式と、プレイヤー分布 $\rho$ の進化に対する連続の式（コルモゴロフ–フォッカー–プランク方程式の一次版）からなる、結合された前向き・後ろ向き偏微分方程式（PDE）系として記述される。

中心となる課題は、ハミルトニアン $H(x, p, \rho)$ の非分離構造、すなわち一般化された運動量 $p$ が分布 $\rho$ と加法的に分離できない点にある。さらに、本解析は、任意の時間ホライゾンおよび（ディラック測量のような）特異な初期分布を許容しており、これは分離構造、正則化ノイズ、あるいは絶対連続な初期測量を要求することが多かった先行研究の制限を越えるものである。

手法
著者らは、MFG系の確率的表現に基づいた粒子ベースの数値解法を提案している。このアプローチの核となるのは以下のプロセスである：

1. ハミルトニアン系の定式化: 連続的なMFG系は、非原子的な確率空間 $(\Omega, \mathcal{F}, P)$ 上で定義された決定論的な前向き・後ろ向きハミルトン系に関連付けられる。この系は、プレイヤーの軌跡 $X(t)$ に関する前向き方程式と、運動量のフロー $Y(t)$ に関する後ろ向き方程式で構成され、これらは $X(t)$ の法則（$\rho(t) = X(t)\#P$）を通じて結合されている。ハミルトニアン $H$ および終端コスト $g$ に対する変位単調性の仮定の下で、この系は一意のグローバル解を持つ。
2. 離散化スキーム: 数値スキームは、時間におけるインプリシット・オイラー離散化と空間におけるサンプリングを採用している。
   • 時間離散化: 時間区間 $[0, T]$ をステップサイズ $\tau_k$ で分割する。スキームは、左連続（$V^+_k$）および右連続（$V^-_k$）な空間を通じて定義される、状態および運動量のフローの区分的に定数な近似を用いる。
   • 空間離散化: 初期分布 $\rho_0$ は、$N$ 個の粒子を持つ決定論的で有限支持を持つ一連の測度 $\rho^N_0$ によって近似される。対応する初期変数 $X^N_0$ は、$N$ 個の異なる値を取る単純な確率変数である。
   • 離散系: 離散的なハミルトン系を解くために、積分をインプリシット・オイラー法を用いて近似する。離散解 $(X^{k,N}, Y^{k,N}, \rho^{k,N})$ は、離散的な軌跡の法則が連続的なフローを近似するように構築される。
3. 反復解法: 結果として得られる非線形な離散前向き・後ろ向き系を解くために、著者らは2つのアルゴリズムを提案している：
   • アルゴリズム A: 標準的なグローバル時間におけるピカール反復。
   • アルゴリズム B: Chassagneux–Crisan–Delarueに触発されたハイブリッド・アプローチであり、グローバルな不動点スキーム内でのローカル時間におけるピカール反復を利用する。これは、グローバルなピカール反復が失敗する可能性のある大きな時間ホライゾンにおいて、安定性と収束を確保するために設計されている。

主要な貢献と結果

• 離散スキームのウェルポーズドネス（適切設定性）: 著者らは、十分小さな時間ステップ $\tau_k$ および任意の空間サンプル数 $N$ に対して、離散数値スキームの解の存在と一意性を証明した。これは、アフィン摂動下での一般化された離散ハミルトン系の連続依存性境界を確立することによって達成される。決定的なことは、離散スキームが連続問題から変位単調性の性質を継承していることであり、これがハミルトニアンの分離性を必要とせずに安定性と一意性を保証している。
• 収束解析: 論文は、時間ステップ $\tau_k \to 0$ および空間サンプル数 $N \to \infty$ の両方における数値近似の強収束を確立している。
  • 誤差境界: プレイヤー分布については $L^\infty(W_2)$ メトリックにおいて、また最適軌跡に沿った値関数の勾配については $L^\infty(L^2)$ ノルムにおいて、漸近的な誤差境界を導出している。境界は $C(\|X_0 - X^N_0\|_{L^2} + \tau_k)$ の形式をとる。
  • 値関数の近似: ラグランジアンおよび終端コストに関する追加的な局所リプシッツ連続性の仮定の下で、離散値関数に対する $L^\infty(L^1)$ ノルムにおける誤差境界が確立されている。
• 明示的な収束率: 著者らは、空間近似の選択に基づく明示的な収束率を導出している：
  • 半離散最適輸送写像によって生成される近似の場合、収束率は次元 $d$ と初期測度のモーメント $q$ に依存する（例：$O(N^{-1/d} + \tau_k)$）。
  • ベクトル量子化による近似の場合、初期測量が非自明な絶対連続部分を持つならば、収束率は $O(N^{-1/d} + \tau_k)$ である。
• 特異な初期データの処理: 本手法は、既存の多くの一次MFG向け数値手法が苦戦したり失敗したりするシナリオである、特異な初期分布（例：ディラック質量）を扱うのに十分な堅牢性を備えている。

意義と主張
著者らは、本研究が以下の事項を同時に扱う、一次・非分離型MFGの最初の厳密な数値解析であると主張している：

1. 非分離型ハミルトニアン（$H$ が $H_0(x,p) - F(x,\rho)$ と分解できない場合）。
2. 潜在的に特異な初期エージェント分布。
3. 任意の長い時間ホライゾン。

その意義は、変位単調な数値スキームの開発にある。ハミルトニアンの小ささに関する条件（変位単調性が欠如している場合に通常必要とされる条件）を回避するために、連続系の変位単調性を活用して離散系の安定性と一意性を証明することで、この成果を得ている。本研究は、これらの複雑な系の解の存在に関する理論と、実用的な数値近似との間の溝を埋めるものである。これは、変位単調性が欠如している場合には解析的・数値的に困難とされる、一次・非分離の設定における一連の問題に対して、理論的に予測された収束率を提供するものである。

数値実験は、提案されたアルゴリズム（特にアルゴリズム B）が、特異な初期データを含む様々な空間次元（$d \in \{1, \dots, 6\}$）および大きな時間ホライゾンに対しても、安定性と予測された収束率を維持することを実証し、理論的知見を検証している。