[論文レビュー] Split-step Fourier methods for the Gross-Pitaevskii equation

Juha Javanainen, Janne Ruostekoski|arXiv (Cornell University)|Nov 5, 2004

Photonic and Optical Devices参考文献 2被引用数 28

ひとこと要約

この論文は、非線形項に最新の波動関数を用いる場合、Gross-Pitaevskii方程式（GPE）に対する分割ステップフーリエ法が、線形シュレーディンガー方程式に対するそれらと同等の時間刻み精度を達成することを確立している。記号的計算を用いて、著者らは、この単純なルールが、高次手法を含むすべてのテスト済みの分割ステップスキームで、アルゴリズムの精度順序を保持することを証明し、任意の空間次元における最小分割ステップ法に対して普遍的に成り立つと仮説を立てている。

ABSTRACT

We perform a systematic study of the accuracy of split-step Fourier transform methods for the time dependent Gross-Pitaevskii equation using symbolic calculation. Provided the most recent approximation for the wave function is always used in the nonlinear atom-atom interaction potential energy, every split-step algorithm we have tried has the same-order time stepping error for the Gross-Pitaevskii equation and the Schroedinger equation.

研究の動機と目的

• 時間に依存するGross-Pitaevskii方程式（GPE）に対する分割ステップフーリエ法の精度を体系的に分析すること。GPEは、超低温原子物理学における重要な方程式である。
• GPEの数値実装における長年の曖昧さを解消すること。多くの研究グループが、収束保証が明確でない、恣意的な未検証手法を用いている。
• 線形シュレーディンガー方程式に対して優れた精度特性を示す分割ステップフーリエ法の性質が、非線形GPEに拡張された場合にも保たれるかどうかを特定すること。
• 特定のアルゴリズムに依存せず、任意の分割ステップGPEソルバーで高次時間精度を保証する、普遍的で計算コストが低いルールを同定すること。

提案手法

• 時間刻み $ h $ における記号的べき級数展開を用いて、GPEの正確な時間発展演算子と分割ステップ近似を体系的に比較した。
• 非可換な運動エネルギー演算子とポテンシャルエネルギー演算子の和の指数関数を、$ h $ ごとの順序で交換子項を追跡しながら演算子代数を用いて展開した。
• 誤差項を所望の順序まで打ち消すために、分割ステップ列の係数 $ \alpha_i $ と $ \beta_i $ に関する多変数多項式方程式を定式化した。
• 指数関数をべき級数に展開し、波動関数に逐次的に適用することで分割ステップアルゴリズムを実装した。記号的変換にはMathematicaを用いた。
• 位置ステップと運動量ステップを初期ステップとして用いた $ \mathcal{O}(h^3) $、$ \mathcal{O}(h^4) $、$ \mathcal{O}(h^5) $ の分割ステップスキームを複数テストした。
• 1次元の2階微分をラプラシアンに置き換えることで、多次元GPEに同様の手法を適用し、次元間での一貫性を確認した。

実験結果

リサーチクエスチョン

• RQ1 分割ステップフーリエ法を非線形Gross-Pitaevskii方程式に適用した場合の時間刻み誤差の次数は何か？
• RQ2 同じアルゴリズムが適用された場合、線形シュレーディンガー方程式に対して優れた精度特性を示す分割ステップ法の性質が、非線形GPEに対しても保たれるか？
• RQ3 GPEの非線形項 $ g|\psi|^2 $ を処理するための、線形の場合と同等の精度順序を保証する特定のルールは存在するか？
• RQ4 精度はステップの順序（例えば、位置または運動量から始めるか）や分割ステップスキームの係数の選択に依存するか？
• RQ5 観察された精度の保持は、任意の空間次元および高次分割ステップ法へ一般化可能か？

主な発見

• 分割ステップフーリエ法がGPEに適用された場合、非線形ポテンシャル項 $ g|\psi|^2 $ に利用可能な最新の波動関数を用いる限り、線形シュレーディンガー方程式の場合と同等の時間刻み誤差次数を達成する。
• テスト済みのすべての分割ステップスキーム（$ \mathcal{O}(h^3) $、$ \mathcal{O}(h^4) $、$ \mathcal{O}(h^5) $ を含む）において、非線形項で $ c_0 = 0 $、$ |c_1| = 1 $ とし、$ |\psi|^2 $ に最新の $ \psi $ を使用することで、$ \mathcal{O}(h^3) $ の精度が保証される。
• 最新の波動関数を使用するという要件は、非線形ケースにおける誤差次数を線形ケースと同等に保つために、単なる十分条件ではなく、テスト済みのスキームに基づくと必要条件であるように思われる。
• 著者らは、このルール（$ |\psi|^2 $ に最新の $ \psi $ を使用する）が、次数や次元にかかわらず、すべての最小分割ステップ法に普遍的に成り立つと仮説を立てている。
• 2次元および3次元でも、この方法の精度は保たれ、$ \mathcal{O}(h^3) $ の3指数スキームについて明示的な検証で確認された。
• この頑健性の背後にある理由は代数的であり、まだ説明がつかない。これは、ノルムを保存する非線形微分方程式に、より深い数学的構造が存在する可能性を示唆している。