[論文レビュー] On the Breakdown of the Ehrenfest Method for Molecular Dynamics on Surfaces

Ignacio Loaiza, Artur F. Izmaylov|arXiv (Cornell University)|Sep 11, 2018

Spectroscopy and Quantum Chemical Studies参考文献 26被引用数 10

ひとこと要約

本稿は、Landau-Zener式およびMasseyパラメータに基づく基準を導出することで、表面における分子動力学におけるEhrenfest（EH）法の破綻を同定し、定量的に評価する。非断熱結合が強く、核の速度が大きい場合にはEH法が失敗することを示し、表面ダイナミクスの複数電子状態を正確にシミュレートするため、より精度の高い手法（例えば、表面遷移法や集団変数法）を採用する必要がある状況を定量的に示す指標を提供する。

ABSTRACT

Due to a continuum of electronic states present in periodic systems, the description of molecular dynamics on surfaces poses a serious computational challenge. One of the most used families of approaches in these settings are friction theories, which are based on the Ehrenfest (EH) approach. Yet, a mean-field treatment of electronic degrees of freedom in the EH method makes this approach inaccurate in some cases. Our aim is to clarify when EH breaks down for molecular dynamics on surfaces. Answering this question provides limits of applicability for more approximate friction theories derived from EH. We assess the EH method on one-dimensional, numerically exactly solvable models with a large but finite number of electronic states. Using the Landau-Zener formula and the Massey parameter, an expression that determines when EH breaks down is deduced.

研究の動機と目的

• 非断熱効果による表面におけるEhrenfest法の失敗条件を明確化すること。
• 非断熱結合および核の運動の両方に基づいた、EH法の破綻の定量的基準を導出すること。
• EH法に基づく摩擦理論に代わって、より精度の高い手法（表面遷移法や集団変数法など）を採用すべき時期を実用的な指標として提供すること。
• 有限な電子状態を有する一次元モデルを用いて、現実的な表面ダイナミクスを模擬することで、EH法の妥当性を評価すること。
• EH法の破綻条件を、断熱的MasseyパラメータおよびLandau-Zener遷移確率に関連付けること。

提案手法

• 表面ダイナミクスを模擬するため、有限だが大きな数の電子状態を有する一次元数値的可解モデルを用いる。
• 核の運動を伴う非断熱遷移確率を計算するために、Landau-Zener式を適用する。
• 電子状態の断熱的ポテンシャルエネルギー曲線間の力の差から、断熱的Masseyパラメータを導出する。
• MasseyパラメータをLandau-Zener遷移率および核の速度に関して再定式化し、EH法の破綻条件を特定する。
• 時間に依存するシュレーディンガー方程式を断熱状態に射影することで、EH法および表面遷移法の運動方程式を導出する。
• クロスするポテンシャルエネルギー曲線を有するモデル系においてEH法と表面遷移法のダイナミクスを比較し、破綻基準の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 Ehrenfest法は、表面における分子動力学でどのような条件下で破綻するか？
• RQ2 非断熱結合および核の運動を用いて、Ehrenfest法の破綻を定量的に予測するにはどうすればよいか？
• RQ3 核の速度は、複数の電子状態を有する系においてEhrenfest法の失敗に果たす役割は何か？
• RQ4 断熱的Masseyパラメータは、表面ダイナミクスにおけるEhrenfest法の妥当性とどのように関係するか？
• RQ5 Landau-Zener式を用いて、Ehrenfest法を放棄してより精度の高い非断熱ダイナミクス手法に移行すべき時期を信頼性を持って導出できるか？

主な発見

• 断熱的Masseyパラメータが臨界閾値未満に下がると、Ehrenfest法は破綻する。これは強い非断熱結合を示している。
• 破綿基準は ξ = π∆² / δ として導出され、ここで ∆ は断熱的結合、δ はエネルギー差である。この式は、直接的にLandau-Zener遷移確率に関連している。
• 核の速度は極めて重要な役割を果たす。高速な核運動は、 avoided crossing に滞在する時間を短くし、非断熱遷移の適切な取り扱いが不十分になるため、EH法の失敗確率を高める。
• 導出された基準は、有限な電子状態を有する一次元モデルにおいて、失敗の発生を正しく特定でき、予測ツールとしての有効性を検証している。
• 複数の電子状態が、それぞれ異なる核依存性を有し、エネルギー的に近接した状態が同時に励起可能になると、EH法の失敗は最も顕著になる。特に、それらの状態の占有数が同程度である場合には顕著である。
• 本研究は、EH法に基づく摩擦近似が、これらの条件下では信頼できないことを確認しており、表面遷移法や集団電子変数を用いた高度な手法の採用が正当化されることを示している。