Resonance-based integrators for stochastic Schrödinger equations. Convergence and long-time error bounds

技術的サマリー：確率シュレーディンガー方程式のための共鳴に基づく積分法

問題定義
本研究は、加法的な $Q$-ウィーナー雑音に駆動される確率シュレーディンガー方程式（SSE）の数値近似を取り扱う。研究は、粗いポテンシャルを持つ線形方程式（$F(\xi) = V(x)\xi$）と、3 次非線形方程式（$F(\xi) = |\xi|^2\xi$）の 2 つの具体的なケースに焦点を当てている。この分野における主要な課題は、確率的畳み込みによって引き起こされる解の低い時間的正則性であり、これは通常、標準的な数値スキームの収束次数を制限する。さらに、そのような方程式に対する既存の理論的解析は、ほとんどが短時間スケールに限定されている。本論文は、低空間的正則性の仮定のもとでより高い収束次数を達成する共鳴に基づく低正則性積分法を開発・解析し、決定的に、小さなポテンシャルと雑音の領域において一様な長時間誤差 bound を確立することを目的としている。

手法
提案された手法は、フーリエ解析と共鳴に基づく数値積分法の構築に依存している。中核的な戦略は以下の通りである：

1. ねじれた変数（Twisted Variables）： 解 $\xi(t)$ をねじれた変数 $v(t) = e^{-it\Delta}\xi(t)$ に変換し、支配的な振動項を分離する。
2. 共鳴項の厳密積分： 積分法は、非線形相互作用またはポテンシャル相互作用のフーリエ展開に現れる支配的な項を厳密に積分することで構築される。これには、振動関数の時間積分（例：$\int_0^\tau e^{is(\ell_1^2 + 2\ell_1\ell_2)} ds$）を近似することが含まれる。ここで、共鳴の場合（位相が消滅する場合）には厳密な積分を保持し、非共鳴の場合には近似を行う。
3. 正則性補償振動（RCO）技法： 決定論的設定で以前に使用された RCO 技法の新たな拡張が、確率的枠組みに適応される。この技法は、確率的畳み込みによって引き起こされる時間的正則性の損失（通常 $O(\sqrt{\tau})$）を克服するために不可欠である。非共鳴項における振動の打ち消しを利用し、高周波数に対するカットオフパラメータを用いることで、著者らは改善された誤差評価を導き出す。
4. 切断戦略： 指数積分性の欠如により強収束が妨げられる 3 次非線形の場合において、著者らは経路ごとの収束結果を導出するために、低正則性技法と組み合わせた切断ベースの戦略を採用する。

主要な貢献と結果

• 粗いポテンシャルを持つ線形方程式：

  • 収束： 共鳴型低正則性積分法（SLR1）が構築される。これは、$H^{\sigma+\gamma}$ 内の解に対して、$L^2(\Omega; H^\sigma)$ において次数 $\gamma \in (0, 1]$ の強収束を達成する。注目すべきは、$H^\sigma$ における 1 次収束のために、この手法が初期データとポテンシャルに対して $H^{\sigma+1}$ の正則性のみを要求し、従来の $H^{\sigma+2}$ の要件を改善している点である。
  • 長時間 bound： $O(\varepsilon^2)$ のポテンシャルと $O(\varepsilon)$ の雑音を持つ領域において、厳密解は $O(\varepsilon^{-2})$ までの時間に対して一様なモーメント bound を持つことが示される。この領域向けに非共鳴積分法（SNRLR1）が設計され、$L^2(\Omega; L^2)$ において $O(\varepsilon^2 \tau)$ の一様な長時間誤差 bound を示し、従来の手法よりも大幅な改善が見られる。

• 3 次非線形方程式：

  • 収束： 共鳴型（SLR2）と非共鳴型（SNRLR2）の両方の積分法が導入される。低正則性において経路ごとの収束結果が導出され、$H^{\sigma+\gamma}$ 内の解に対して $H^\sigma$ において次数 $\gamma \in (0, 1)$ が達成される。
  • 長時間 bound： 弱非線形領域（$O(\varepsilon^2)$ の非線形性、$O(\varepsilon)$ の雑音）において、非共鳴積分法（SNRLR2）は、$O(\varepsilon^{-2})$ までの時間に対して、任意の $\delta < 1$ に対して $O(\varepsilon^2 \tau^\delta)$ の一様な長時間経路ごとの誤差 bound を満たすことが証明される。

• 理論的革新：

  • 本論文は、RCO 技法を確率分散方程式に拡張した最初の研究である。この拡張により、著者らは確率的雑音によって課される標準的な $\sqrt{\tau}$ の制限を克服する改善された誤差評価を回復でき、他の低正則性積分法と比較して長時間誤差を $\varepsilon^2$ 倍に低減している。

意義と主張
著者らは、この研究が確率分散方程式に適用される低正則性積分法に対する長時間誤差 bound の確立を初めて行ったと主張している。その意義は以下の点にある：

1. 低減された正則性要件： 1 次収束を達成するために必要な空間的正則性が低減され（$H^{\sigma+1}$ 対 $H^{\sigma+2}$）、より粗いポテンシャルや初期データに適用可能となっている。
2. 長時間安定性： $O(\varepsilon^{-2})$ の時間スケールにわたる一様な誤差 bound の導出は、以前の結果が短時間スケールに限定されていた文献のギャップを埋めるものである。
3. 確率的 RCO 拡張： 確率的設定への RCO 技法の成功した適応は、確率的畳み込み固有の時間的不規則性を処理するための新たな解析ツールを提供する。

数値的検証
理論的結果は数値実験によって確認されている。これらの実験は以下を実証する：

• $H^1$ 正則性のもとでの線形ケースにおける、次数 1 の強収束率および経路ごとの収束率。
• 共鳴型（SLR1）および古典的指数法で観察される誤差の増大とは対照的に、長時間区間（$O(\varepsilon^{-2})$）における非共鳴積分法（SNRLR1）の誤差の有界性。
• 小さなポテンシャル/雑音領域において、非共鳴スキームの誤差の大きさにおける $O(\varepsilon^2)$ の低減。
• 提案された積分法を用いた 3 次非線形ケースにおける同様の収束および長時間挙動。