Law of Large Numbers and Central Limit Theorem for random sets of solitons of the focusing nonlinear Schrödinger equation

技術的概要：集束型非線形シュレーディンガー方程式のソリトン集合に対する大数の法則と中心極限定理

問題提起
本論文は、1 次元空間における立方集束型非線形シュレーディンガー（fNLS）方程式の解の統計的挙動を取り扱います：
$$i\psi_t + \frac{1}{2}\psi_{xx} + |\psi|^2\psi = 0, \quad x \in \mathbb{R}, t \in \mathbb{R}^+.$$
線形方程式における無相関な波の重ね合わせを通じて、ランダムな初期データの進化はよく理解されていますが、非線形の場合には重大な課題が生じます。非線形波では、波場の確率分布が時間とともに著しく変形します。従来の研究は、主に波動運動方程式によって記述される弱非線形領域に焦点を当ててきました。しかし、ソリトンに代表されるような大振幅の根本的に非線形な解に対しては、厳密な統計理論が存在していません。

著者らは、スペクトルデータ（固有値と規格化定数）がランダム化された $N$ ソリトン解を調査します。具体的には、固有値 $\{\lambda_k\}_{k=1}^N$ は、上半平面 $\mathbb{C}_+$ 内のコンパクトな台を持つ確率分布からサンプリングされた独立同分布（i.i.d.）の確率変数です。規格化定数は、これらの固有値の滑らかな補間関数によって決定されます。中心的な問いは、$N \to \infty$ において、ランダムな $N$ ソリトン解が決定的な極限（「ソリトンガス」）に収束するかどうか、そしてこの極限からの揺らぎが中心極限定理（CLT）に従うかどうかです。

手法
解析は、逆散乱法（IST）による fNLS 方程式の可積分性に依存しています。解は、リーマン・ヒルベルト（RH）問題を用いて散乱データから再構成されます。

1. ランダム RH 定式化：著者らは、$N$ ソリトン解をランダムな双有理 RH 問題として定式化します。漸近解析を容易にするため、離散固有値に関連する極条件を、固有値分布の台を囲む滑らかな輪郭 $\gamma_+$ および $\gamma_-$ におけるジャンプ関係に変換します。これにより、ランダムなジャンプ行列 $J_N(z; x, t, \lambda)$ が得られます。
2. 平均化された決定的極限：ランダムなジャンプ行列の期待値 $J(z; x, t) = \mathbb{E}[J_N(z; x, t, \lambda)]$ を取ることで、決定的な「平均化された」RH 問題が構築されます。固有値の i.i.d. 性質と規格化定数の補間により、この平均化されたジャンプ行列は $N$ に依存しません。この決定的 RH 問題の解は、ソリトンガス解として解釈される極限関数 $\psi_\infty(x, t)$ を定義します。
3. 誤差解析：ランダムな解 $\psi_N$ と決定的極限 $\psi_\infty$ の差は、誤差行列 $E(z) = M_N(z)M(z)^{-1}$ を調べることで解析されます。ここで、$M_N$ と $M$ はそれぞれランダムおよび平均化された RH 問題の解です。
4. 小ノルム議論：誤差問題のジャンプ行列 $J_E$ は、単位行列の摂動 $J_E = I + W_N$ であることが示されます。項 $W_N$ はランダム固有値の線形統計量に依存します。確率的な評価（特に線形統計量 $X_N^f$ に関する境界）を用いて、著者らは $N \to \infty$ において、高い確率で $W_N$ のノルムが小さくなることを示します。これにより、誤差行列 $E$ を収束するノイマン級数として展開することが可能になります。
5. 確率的極限：
   • 大数の法則（LLN）：ノイマン級数の誤差項を境界付け、「悪い」構成（ノルムが小さくない場合）の確率を制御することで、著者らは $\psi_N$ および $|\psi_N|^2$ がそれぞれ $\psi_\infty$ および $|\psi_\infty|^2$ に平均収束（$L^1$）することを証明します。
   • 中心極限定理（CLT）：スケーリングされた差 $\sqrt{N}(\psi_N - \psi_\infty)$ の展開における主要項は、ランダム固有値の線形統計量として特定されます。古典的な CLT により、この項はガウス確率変数に収束します。ノイマン級数の高次項は確率的に無視できることが示されます。

主要な貢献と結果

• 決定的ソリトンガスの存在：本論文は、スペクトルデータの期待値から導かれる fNLS 方程式の一意な古典解 $\psi_\infty(x, t)$ の存在を確立します。この解は、ランダムなソリトンガスの巨視的極限として解釈されます。
• 大数の法則（定理 2.6）：著者らは、$\mathbb{R} \times \mathbb{R}^+$ のコンパクト集合内の $(x, t)$ に対して、ランダムな $N$ ソリトン解 $\psi_N(x, t; \lambda)$ およびそのモジュラスの二乗が、それぞれ決定的解 $\psi_\infty(x, t)$ および $|\psi_\infty(x, t)|^2$ に平均収束することを証明します：
  $$\lim_{N \to \infty} \mathbb{E}[|\psi_N - \psi_\infty|] = 0, \quad \lim_{N \to \infty} \mathbb{E}[||\psi_N|^2 - |\psi_\infty|^2|] = 0.$$
• 中心極限定理（定理 2.7）：本論文は、平均からの揺らぎが $\sqrt{N}$ でスケーリングされた場合、分布においてガウス確率変数に収束することを証明します。具体的には、$\sqrt{N}(\psi_N - \psi_\infty)$ は複素ガウス変数 $X_{G1}$ に収束し、$\sqrt{N}(|\psi_N|^2 - |\psi_\infty|^2)$ は実ガウス変数 $X_{G2}$ に収束します。両者とも平均はゼロであり、その分散は、平均化された RH 問題の解および補間関数を含む積分として明示的に計算されます。
• 相関関数（定理 2.8）：著者らは揺らぎの 2 点相関関数を計算し、それらが極限ガウス変数の共分散に収束することを示します。

意義と主張
本論文は、基礎となる非線形偏微分方程式の解析を通じて、集束型 NLS 方程式の文脈におけるランダムなソリトンガスに対する運動論的理論の結果を初めて厳密に導出したと主張しています。ソリトンガスに対しては、一般化流体力学（Generalized Hydrodynamics）などの運動論的理論が提案され、数値的に検証されてきましたが、この研究は、この設定における大数の法則と中心極限定理に対する厳密な数学的基盤を提供します。

著者らは、散乱データの線形設定にランダム性を導入するアプローチを採用しており、これは解が進化するにつれて無相関かつ分布を変えずに維持されるため、ランダムな大振幅非線形ダイナミクスの研究を可能にしていると強調しています。結果は、個々のソリトンの微視的記述とソリトンガスの巨視的統計的挙動の間のギャップを埋め、広大な空間・時間スケールにおける大振幅波に対する予測的な統計理論を提供します。この研究は、可積分な非線形分散系の統計力学を理解するための厳密な一歩として提示されています。