[論文レビュー] Convergence of a Normalized Gradient Algorithm for Computing Ground States

Erwan Faou, Tiphaine Jézéquel|arXiv (Cornell University)|Mar 8, 2016

Advanced Mathematical Physics Problems参考文献 14被引用数 19

ひとこと要約

本稿は、1次元立方非線形シュレーディンガー方程式の基底状態を計算するための正規化勾配法（虚時間法）の収束性を証明する。中心差分による空間離散化と線形的陰的時間積分法を用いることで、離散ソリトンが空間メッシュサイズ $h$ および領域切断 $K$ に依存する誤差境界を伴い、指数関数的に収束する。主な貢献は、$H^1$ ノルムにおける明示的な誤差推定を含む厳密な収束解析である。

ABSTRACT

We consider the approximation of the ground state of the one-dimensional cubic nonlinear Schr{\\"o}dinger equation by a normalized gradient algorithm combined with linearly implicit time integrator, and finite difference space approximation. We show that this method, also called imaginary time evolution method in the physics literature, is con-vergent, and we provide error estimates: the algorithm converges exponentially towards a modified solitons that is a space discretization of the exact soliton, with error estimates depending on the discretization parameters.

研究の動機と目的

• 1次元立方非線形シュレーディンガー方程式の基底状態を計算するための正規化勾配（虚時間）法の収束性を厳密に確立すること。
• 空間的および時間的離散化の下で、計算された解と正確なソリトンとの間の誤差を分析すること。
• 線形的陰的時間スキームが離散基底状態を正確に保存することを示し、安定性および精度を保証すること。
• メッシュサイズ $h$、時間ステップ $\tau$、領域切断 $K$ に依存する完全離散アルゴリズムの定量的誤差推定を提供すること。

提案手法

• アルゴリズムは、エネルギー汎関数の $L^2$-勾配の負の方向に沿う放物型方程式 $\partial_t \psi = \frac{1}{2}\Delta\psi + |\psi|^2\psi$ を解くために線形的陰的時間積分法を用いる。
• スキームは $\psi_n^* = \psi_n - \tau \widehat{\nabla H}(\psi_n, \psi_n^*)$ で定義され、ここで $-\widehat{\nabla H}(\psi_n, \psi_n^*) = \frac{1}{2}\Delta_h \psi_n^* + |\psi_n|^2 \psi_n^*$ であり、無条件安定性と基底状態の正確な保存を保証する。
• 各時間ステップ後に $\psi_{n+1} = \psi_n^* / \|\psi_n^*\|_{L^2}$ により $L^2$ ノルムで正規化することで、単位質量を維持する。
• 空間離散化は、$|x| \leq Kh$ の領域で中央差分を用い、境界条件をゼロとする。収束は離散 $H^1$ ノルムで測定される。
• 離散基底状態 $\eta_{h,K}$ は、単位 $L^2$ ノルムの下でのエネルギー汎関数の唯一の最小化子として定義され、離散固有値方程式を満たす。
• 安定性および凸性の解析のため、$\eta_{h,K}$ の周囲に局所座標系を構築し、指数的収束の証明を可能にする。

実験結果

リサーチクエスチョン

• RQ1 線形的陰的時間積分法を用いた正規化勾配法は、空間的および時間的離散化のもとで、1次元立方NLS方程式の真の基底状態に収束するか？
• RQ2 このアルゴリズムの収束速度は何か。また、$h$、$\tau$、$K$ の離散化パラメータにどのように依存するか？
• RQ3 線形的陰的スキームは離散基底状態を正確に保存できるか。この性質は収束性および安定性を向上させるか？
• RQ4 離散ソリトン $\eta_{h,K}$ と正確な連続ソリトン $\eta$ の間の $H^1$ ノルムにおける近さはどの程度か？
• RQ5 代替スキーム（半陽的または完全陽的）についても収束するか。離散設定における誤差境界は何か？

主な発見

• アルゴリズムは時間に関して指数関数的に離散基底状態 $\eta_{h,K}$ に収束し、$\|\psi_n - \eta_{h,K}\|_h \leq C e^{-cn\tau}$ を満たす。ここで $C, c > 0$ は $h$、$K$、$=\tau$ に依存しない定数である。
• 離散基底状態 $\eta_{h,K}$ は離散固有値方程式 $\frac{1}{2}(\Delta_h \eta_{h,K})^\ell + |\eta_{h,K}^\ell|^2 \eta_{h,K}^\ell = \lambda_h \eta_{h,K}^\ell$ を満たし、スキームによって正確に保存される。
• 連続ソリトン $\eta$ と離散ソリトン $\eta_{h,K}$ の間の誤差は $\|i_h \eta_{h,K} - \eta\|_{H^1} \leq C(h + \frac{1}{h^2} e^{-C_1 h K})$ で抑えられ、$K$ に関してスペクトル収束、$h$ に関して代数的収束を示す。
• 完全離散解 $\psi_n$ は $\|i_h \psi_n - \eta\|_{H^1} \leq C(e^{-cn\tau} + h + \frac{1}{h^2} e^{-C_1 h K})$ を満たし、時間的指数減衰と空間離散化誤差が組み合わされる。
• 線形的陰的スキームは基底状態を正確に保存する。この性質は、半陽的または完全陽的スキームには共有されず、数値的安定性および精度を向上させる。
• 代替スキーム（半陽的または完全陽的）についても、修正された離散基底状態 $\eta_{\tau,h,K}$ に収束し、誤差 $\|i_h \eta_{\tau,h,K} - \eta\|_{H^1} \leq C(\tau + h + \frac{1}{h^2} e^{-C_1 h K})$ を満たす。時間離散化が $\tau$ に依存する追加誤差をもたらすことが示された。