[論文レビュー] Suppression of reflection from the grid boundary in solving the time-dependent Schroedinger equation by split-step technique with fast Fourier transform

Arkady Gonoskov, Ivan Gonoskov|ArXiv.org|Jul 12, 2006

Electromagnetic Scattering and Analysis参考文献 1被引用数 3

ひとこと要約

本稿では、時間に依存するシュレーディンガー方程式を解くために、高速フーリエ変換を用いた新しい分割ステップフーリエ法を提案し、境界処理を組み込むことで、グリッド端縁からの反射を任意の精度で抑える。この手法により、波長が長い波動関数を含むシミュレーションにおいて特に重要な、長波長波の効果的吸収が可能になる。吸収領域の特性を入射波のスケールに合わせて調整することで実現される。

ABSTRACT

We present an approach to numerically solving the time-dependent Schroedinger equation and other parabolic equations by the split-step technique with fast Fourier transform, which suppresses the backreflection of waves from the grid boundaries with any specified accuracy. Most importantly, all known methods work well only for a narrow region of incident waves spectrum, and the proposed method provides absorption of any wave whose length is large enough in comparison with the size of absorption region.

研究の動機と目的

• 時間に依存するシュレーディンガー方程式の数値解法における計算グリッド境界での波の反射という、長年の課題に取り組む。
• 従来の吸収境界条件が、入射波のスペクトル範囲に限定してしか有効でないという制限を克服する。
• 吸収領域のサイズよりもはるかに大きな波長を持つ波に対して、反射を抑えることができる手法を開発する。
• 不自然な反射が解を歪めることのない、高精度な量子ダイナミクスの数値シミュレーションを保証する。

提案手法

• 実空間で時間に依存するシュレーディンガー方程式を解くために、高速フーリエ変換を用いた分割ステップフーリエ法を適用する。
• 境界領域における波動関数の振る舞いを変更する、特許的な吸収境界条件を導入する。
• 期待される入射波の波長に基づいて、パrameterを調整可能な複素数ポテンシャルまたは減衰プロファイルを境界領域に導入する。
• 波動パケットの支配的波長に合わせて、吸収層の厚さと強度を最適化し、反射を最小限に抑える。
• スペクトル精度を保ちながら吸収を適用するために、実空間とフーリエ空間の操作を交互に実行する分割ステップ手法を用いる。
• 吸収プロファイルのパrameterを調整することで、反射抑制を任意の精度に達成できることを示す。

実験結果

リサーチクエスチョン

• RQ1 分割ステップFFT法を用いた時間に依存するシュレーディンガー方程式ソルバーにおいて、グリッド境界からの数値的反射をどのように抑制できるか。
• RQ2 標準的な実装では、なぜ従来の吸収境界条件が長波長入射波に対して失敗するのか。
• RQ3 広い波数範囲にわたり反射を効果的に抑えるために、吸収層が満たすべき条件は何か。
• RQ4 境界層のパrameterを調整することで、吸収効率を入射波の波長に依存させないようにできるか。
• RQ5 吸収領域のサイズと吸収される波の波長との関係は何か。

主な発見

• 提案手法は、グリッド境界からの反射を任意の指定精度で抑えることができ、高精度なシミュレーションを可能にする。
• 従来手法とは異なり、吸収領域サイズよりもはるかに大きな波長の波に対しても有効である。
• 分割ステップフーリエ法と高速フーリエ変換を組み合わせた場合でも、数値的安定性とスペクトル精度を維持する。
• 入射波パケットの特徴的な波長に合わせて、吸収プロファイルを体系的に調整可能である。
• 時間に依存するシュレーディンガー方程式を超える、他の放物型偏微分方程式に対しても適用可能である。
• 数値結果により、境界層の設計を適切に行うことで、反射係数を無視できる程度にまで低下させられることを確認した。