[論文レビュー] Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schrödinger equation

Jiří Vaníček|arXiv (Cornell University)|Feb 20, 2023

Advanced Fiber Laser Technologies被引用数 4

ひとこと要約

本稿は、ヘラーの「解凍ガウス」、変分ガウス、および単一ヘッシアン近似などのさまざまなガウス波パッケージ動的法を、状態依存の二次ポテンシャルを有する非線形シュレーディンガー方程式の共通フレームに統一する。一般化された運動方程式を導出し、ノルム、エネルギー、シンプレクティック構造の保存を証明するとともに、局所的立方近似よりも精度を向上させつつ、重要な幾何的性質を低コストで保持する新しい「単一四次」変分法を提案する。

ABSTRACT

Many approximate solutions of the time-dependent Schrödinger equation can be formulated as exact solutions of a nonlinear Schrödinger equation with an effective Hamiltonian operator depending on the state of the system. We show that several well-known Gaussian wavepacket dynamics methods, such as Heller's original thawed Gaussian approximation or Coalson and Karplus's variational Gaussian approximation, fit into this framework if the effective potential is a quadratic polynomial with state-dependent coefficients. We study such a nonlinear Schrödinger equation in general: in particular, we derive general equations of motion for the Gaussian's parameters, demonstrate the time reversibility and norm conservation, and analyze conservation of the energy, effective energy, and symplectic structure. We also describe efficient geometric integrators of arbitrary even orders of accuracy in the time step for the numerical solution of this nonlinear Schrödinger equation. The general presentation is illustrated by examples of this family of Gaussian wavepacket dynamics, including the variational and nonvariational thawed and frozen Gaussian approximations, and special limits of these methods based on the global harmonic, local harmonic, single-Hessian, local cubic, and single quartic approximations for the potential energy. Without substantially increasing the cost, the proposed single quartic variational thawed Gaussian wavepacket dynamics improves the accuracy over the local cubic approximation and, at the same time, conserves both the effective energy and symplectic structure, in contrast to the much more expensive local quartic approximation. Most results are presented in both Heller's and Hagedorn's parametrizations of the Gaussian wavepacket.

研究の動機と目的

• 状態依存の二次ポテンシャルを有する非線形シュレーディンガー方程式に基づく、一貫した理論的枠組みにおいて、多様なガウス波パッケージ動的法を統一すること。
• ガウス波パッケージのパラメータ（中心、広がり行列、位相、運動量）の一般化された運動方程式を導出し、ノルム、有効エネルギー、シンプレクティック構造といった基本的保存則を証明すること。
• 非線形シュレーディンガー方程式の数値的解法に適した高次幾何的積分法を構築すること。
• 局所的立方近似よりも精度を向上させつつ、エネルギーおよびシンプレクティック保存を維持する新しい「単一四次」変分ガウス近似を提案・分析すること。
• 本フレームワークの一般性を示すために、変分的および非変分的解凍・凍結ガウス近似に適用し、グローバル調和、局所的調和、単一ヘッシアン、局所的立方、局所的四次近似の極限を含む。

提案手法

• ガウス波パッケージの動的法を、波パッケージの瞬時の状態に依存するパrameter（二次ポテンシャル係数）を有する有効ハミルトニアンをもつ非線形時間に依存するシュレーディンガー方程式の正確な解として定式化する。
• ガウスアンザッツに時間依存変分原理を適用し、ガウス波パッケージの中心、広がり行列、位相、運動量の一般化された運動方程式を導出する。
• 非線形系の背後にある幾何的・代数的性質を分析することで、時間反転対称性、ノルム保存、有効エネルギーおよびシンプレクティック構造の保存を証明する。
• シンプレクティック構造および有効エネルギーを保存する高次幾何的積分法を導入し、数値シミュレーションにおける長期的安定性と精度を確保する。
• 本フレームワークを、変分的および非変分的解凍・凍結ガウス近似、および調和、立方、四次ポテンシャル近似を用いた極限に適用する。
• 局所的立方近似に四次微分項を追加することで、計算コストを低く保ちつつ精度を向上させ、幾何的不変量を維持する新しい「単一四次」変分法を提案する。

実験結果

リサーチクエスチョン

• RQ1 ヘラーの解凍ガウスやコールソン＝カープルスの変分ガウスといった、さまざまなガウス波パッケージ動的法が、状態依存ポテンシャルを有する非線形シュレーディンガー方程式の枠組みで統一可能か。
• RQ2 この非線形枠組みにおいて、ガウス波パッケージパラメータの一般化された運動方程式は何か。また、ノルム、エネルギー、シンプレクティック構造といった幾何的性質は、どの程度保存されるか。
• RQ3 非線形シュレーディンガー方程式の数値的解法において、シンプレクティック構造および有効エネルギーを保存する高次幾何的積分法をどのように構築できるか。
• RQ4 局所的立方近似よりも精度を向上させつつ、エネルギーおよびシンプレクティック保存を損なわない新しい変分ガウス近似を設計可能か。
• RQ5 ポテンシャルの高次微分（例：四次微分）を含めることで、ガウス波パッケージ動的法における精度および保存性にどのような影響を与えるか。

主な発見

• 有効エネルギーは、非線形シュレーディンガー方程式の枠組みにおいて正確に保存され、複数のガウス波パッケージ動的法を一般化・統一する。
• 運動方程式によってシンプレクティック構造が保存され、数値的シミュレーションにおける長期的安定性と精度が保証される。
• 提案された「単一四次」変分ガウス近似は、局所的立方近似よりも精度を向上させつつ、有効エネルギーおよびシンプレクティック構造を両方とも保存する。
• フレームワークにより、単一ヘッシアン近似が単純であるにもかかわらず、有効エネルギーを正確に保存し、拡張された位相空間におけるハミルトニアン構造を持つことが明らかになった。
• 一般理論は、ヘラーのパrametrizationおよびハゲルドンのパラメータ化の両方に対して有効であり、異なる定式化に広く適用可能である。
• シンプレクティック構造および有効エネルギーを保存する高次幾何的積分法が導出され、ガウス波パッケージの効率的かつ高精度な長時間発展が可能になった。