[論文レビュー] Frozen Gaussian approximation with surface hopping for mixed quantum-classical dynamics: A mathematical justification of fewest switches surface hopping algorithms

Jianfeng Lu, Zhennan Zhou|arXiv (Cornell University)|Feb 20, 2016

Spectroscopy and Quantum Chemical Studies参考文献 28被引用数 6

ひとこと要約

本稿では、半古典的領域における混合量子古典的ダイナミクスに対して、表面遷移付き凍結ガウス近似（FGA-SH）という表面遷移アルゴリズムの数学的に厳密な導出を提示する。経路積分の確率的表現を用いて解を定式化することで、複数のアディアバティック状態を伴うシュレーディンガー方程式のモンテカルロ近似が得られ、一貫性と安定性の両面で誤差がそれぞれ $\mathcal{O}(\varepsilon)$ および $\mathcal{O}(\varepsilon^2)$ のスケーリングを示し、漸近的に正しい結果が得られる。

ABSTRACT

We develop a surface hopping algorithm based on frozen Gaussian approximation for semiclassical matrix Schrödinger equations, in the spirit of Tully's fewest switches surface hopping method. The algorithm is asymptotically derived from the Schrödinger equation with rigorous approximation error analysis. The resulting algorithm can be viewed as a path integral stochastic representation of the semiclassical matrix Schrödinger equations. Our results provide mathematical understanding to and shed new light on the important class of surface hopping methods in theoretical and computational chemistry.

研究の動機と目的

• 非アディアバティック量子ダイナミクスにおける広く用いられている最少スイッチ表面遷移（FSSH）アルゴリズムの、数学的に厳密な基礎を提供すること。
• 計算化学分野で広く用いられているにもかかわらず、表面遷移法に理論的裏付けが欠けている問題に対処すること。
• 冷凍ガウス近似と経路積分の確率的表現に基づいて、半古典的極限において漸近的に正しい結果を得られる新しいアルゴリズム FGA-SH を導出すること。
• FGA-SH 解と行列シュレーディンガー方程式の真の解との間の定量的誤差推定を確立すること。
• 既存の FSSH 方法における概念的曖昧さを解消し、解の明確な確率的表現に結びつけること。

提案手法

• 半古典的行列シュレーディンガー方程式の解に対して、経路積分の確率的表現を導出し、位相空間内での区分的決定的軌道としてのダイナミクスを扱う。
• 経路積分の評価に用いるモンテカルロ法として、FGA-SH アルゴリズムを構築し、異なるポテンシャルエネルギー面に沿って古典的ハミルトニアンの流れに従って軌道が進化するようにする。
• 非アディアバティック結合点で表面遷移イベントを実装し、ハミルトニアン内の導出結合項から得られる遷移確率を用いる。
• 波動関数のアンザッツとして冷凍ガウス近似（ハーマン＝クラークプロパゲーター）を用い、位相空間積分による計算の効率化を実現する。
• 近似シュレーディンガー方程式における剰余項を用いて解の誤差を評価する、新規の安定性補題を用いて一貫性および安定性の推定を確立する。
• FGA-SH 近似における総誤差が、一貫性に関しては $\mathcal{O}(\varepsilon)$、安定性に関しては $\mathcal{O}(\varepsilon^2)$ で有界であり、時間および系のパrameterに依存する定数を伴うことを証明する。

実験結果

リサーチクエスチョン

• RQ1 最少スイッチ表面遷移アルゴリズムは、半古典的領域におけるシュレーディンガー方程式から数学的に厳密に導出可能か？
• RQ2 非アディアバティック遷移を伴う混合量子古典的ダイナミクスの経路積分表現の背後にある数学的構造は何か？
• RQ3 冷凍ガウス近似は、正確な誤差制御が可能な多状態系に表面遷移を含めて拡張可能か？
• RQ4 FGA-SH アルゴリズムと量子古典的リウヴィル方程式の関係は何か？ そして、前者は後者から正当化可能か？
• RQ5 半古典的パrameter $\varepsilon$ に従って、FGA-SH 近似の正確な誤差境界は何か？

主な発見

• FGA-SH アルゴリズムは、半古典的行列シュレーディンガー方程式の解に対する経路積分の確率的表現として導出され、表面遷移に数学的に堅固な根拠を提供する。
• FGA-SH 近似の一貫性誤差は、$C\varepsilon \frac{(C_t)^n}{n!} + C\varepsilon^2 e^{C_t}$ で有界であり、$n \to \infty$ のとき 0 に収束するため、漸近的に正しい結果が保証される。
• 安定性補題により、近似解と真の解との間の誤差は、$\frac{1}{\varepsilon} \int_0^t \mu^\varepsilon(s) ds$ のオーダーでしか増大しないことが示され、$\mu^\varepsilon(s)$ は近似方程式における剰余項を制御する。
• FGA-SH 近似における総誤差は、一貫性に関して $\mathcal{O}(\varepsilon)$、安定性に関して $\mathcal{O}(\varepsilon^2)$ で有界であり、時間および系のパrameterに依存する定数を伴う。
• 本手法により、最少スイッチ表面遷移アルゴリズムに対する理論的裏付けが明確にされ、非アディアバティックダイナミクスにおける長年の理論的懸念が解消された。
• 本アプローチにより、表面遷移が高次元の経路積分の直接的なモンテカルロ評価として解釈可能であることが明らかになり、アルゴリズムの理解に新たな概念的枠組みが提供された。