[論文レビュー] Numerical computation of soliton dynamics for NLS equations in a driving potential

Marco Caliari, Marco Squassina|ArXiv.org|Aug 25, 2009

Advanced Mathematical Physics Problems参考文献 23被引用数 4

ひとこと要約

本稿では、時間に依存しない外部ポテンシャルを伴う非線形シュレーディンガー方程式におけるソリトンダイナミクスの数値的検証を提示している。半古典的極限（$\varepsilon \to 0$）において、解の重心がニュートンの運動方程式 $\ddot{x} = -\nabla V(x)$ に従う古典的軌道に従うことを示している。スペクトル法と指数型積分法を用いて、ソリトンダイナミクスの挙動を検証し、2次元の調和ポテンシャルにおいて $\mathbb{H}_\varepsilon$ ノルムで $\mathcal{O}(\varepsilon)$ の誤差境界が確認された。

ABSTRACT

We provide some numerical computations for the soliton dynamics of the nonlinear Schrödinger equation with an external potential. After computing the ground state solution $r$ of a related elliptic equation we show that, in the semi-classical regime, the center of mass of the solution with initial datum modelled on $r$ is driven by the solution of a Newtonian type law. Finally, we provide some examples and analyze the numerical errors in the two dimensional case when $V$ is an harmonic potential.

研究の動機と目的

• 外部ポテンシャルを伴う非線形シュレーディンガー方程式におけるソリトンダイナミクスを、半古典的領域で数値的に調査すること。
• 固定された $L^2$ ノルム制約の下でエネルギー汎関数の変分的最小化により、関連する楕円型方程式の基底状態解 $r$ を計算すること。
• 解の重心が古典的軌道 $\ddot{x} = -\nabla V(x)$ に従うという解析的予測を検証すること。
• 調和ポテンシャルを用いた2次元シミュレーションにおける数値誤差を分析すること。

提案手法

• 固定された $L^2$ ノルム制約の下でエネルギー汎関数の変分的最小化により、基底状態 $r$ の数値的計算を行う。
• 空間方向にスペクトルコロケーションを用い、高精度を得るためのフーリエスペクトル法を採用する。
• 時間積分に $\varphi_1$ および $\varphi_2$ 関数に基づく指数型積分法を適用し、2次精度を確保する。
• $\hat{R}_{n+1} - \hat{R}_{n2}$ からの誤差推定値を用いて、時間ステップを適応的に制御することで計算効率を向上させる。
• 数値計算の安定性を確保するため、$\varepsilon \to 0$ の極限において問題をスケーリングする変数変換を実施する。
• $\mathbb{H}_\varepsilon$ ノルム（$\|\phi\|_{\mathbb{H}_\varepsilon}^2 = \varepsilon^{2-N}\|\nabla\phi\|_{L^2}^2 + \varepsilon^{-N}\|\phi\|_{L^2}^2$）を用いた誤差解析により、収束性を定量的に評価する。

実験結果

リサーチクエスチョン

• RQ1 外部ポテンシャルを伴うNLS方程式の解の重心は、半古典的極限において、古典的軌道 $\ddot{x} = -\nabla V(x)$ に従うか？
• RQ2 関連する楕円型方程式の基底状態解は、数値的にどの程度の精度で計算可能か？
• RQ3 解をソリトンプロファイル $r((x - x(t))/\varepsilon)$ で近似する際、$\mathbb{H}_\varepsilon$ ノルムにおける数値誤差の大きさはどの程度か？
• RQ4 ポテンシャルの選択（例えば、有理数比・無理数比の調和ポテンシャル）は、ソリトン重心の長期的ダイナミクスにどのように影響するか？
• RQ5 残差誤差推定値に基づく適応的時間ステップ法は、長時間シミュレーションにおける計算効率を向上させることができるか？

主な発見

• 初期データが基底状態 $r$ に従う解の重心は、半古典的領域において、古典的軌道 $\ddot{x} = -\nabla V(x)$ に従うことが確認された。
• 数値的シミュレーションにより、真の解とソリトンプロファイル $r((x - x(t))/\varepsilon)$ の間の $\mathbb{H}_\varepsilon$ ノルム誤差が、解析理論で予測された $\mathcal{O}(\varepsilon)$ で抑えられていることが確認された。
• 有理数比の2次元調和ポテンシャルでは、ソリトン軌道は周期的である。一方、無理数比では、軌道はエルゴディックである。
• ダブルバンプ初期条件ではソリトン衝突が発生し、数値スキームがそのダイナミクスを正確に捉えていることが示された。
• 適応的時間ステップを用いた指数型積分法は、空間方向でスペクトル精度、時間方向で2次精度を達成しており、誤差推定値を用いた効率的な計算が可能である。