[論文レビュー] An Electronically Non-Adiabatic Generalization of Ring Polymer Molecular Dynamics

Timothy J. H. Hele|arXiv (Cornell University)|Aug 19, 2013

Advanced Physical and Chemical Molecular InteractionsChemistry参考文献 7被引用数 17

ひとこと要約

本稿では、電子的非断熱系における量子効果および再通過ダイナミクスの正確なモデル化を可能にする、リングポリマー分子動力学（RPMD）理論の非断熱一般化を提案する。この手法は、修正されたポテンシャルエネルギー表面を通じて電子状態の結合を組み込むことで、古典的RPMDを拡張し、深いつながりトンネル領域ですら正確な量子レートと2倍以内の一致を達成する。計算コストは系のサイズに線形に比例する。

ABSTRACT

In this thesis I generalize Ring Polymer Molecular Dynamics (RPMD) rate theory to electronically non-adiabatic systems, followed by application to two one-dimensional curve crossing models and a multidimensional spin-boson model.

研究の動機と目的

• 核運動が複数の電子的ポテンシャルエネルギー表面の上を進行する電子的非断熱系に、リングポリマー分子動力学（RPMD）速度理論を拡張すること。
• トンネル効果や再通過ダイナミクスといった量子効果を捉える計算的に効率的な手法を開発すること。
• 電子状態が1つである極限において従来の断熱RPMDに還元され、高温極限において物理的に整合性を保つようにすること。
• 正確な量子ベンチマークを用いて、1次元の曲線交差モデルおよび多次元スピンボソンモデルでの手法の妥当性を検証すること。
• 系のサイズに線形に比例する計算コストを維持しながら、他の近似手法を上回る精度を示すことを示すこと。

提案手法

• 量子系を調和スプリングで接続された古典的ビーズのリングポリマーとして表現する古典的同型性を一般化し、電子的自由度を含む拡張された位相空間を導入する。
• ハミルトニアンにおける修正された結合項を通じて、電子的非断熱結合を組み込んだ、リングポリマー用の新しい有効ポテンシャルエネルギー表面を導出する。
• ベンネット＝チャンドラ因子分解を適用し、RPMD速度を、分割表面における重心確率密度と生成物領域における正の運動量を持つ重心のフラックスの積に分解する。
• 非対称な曲線交差モデルにおける最適な分割表面を決定するために変分法を用いる。
• 拡張位相空間上で古典的分子動力学を実装し、リングポリマー重心の時間相関関数から反応速度を計算する。
• スピンボソンモデルにおける量子遷移状態理論（QTST）の発散するレートを克服するため、物理的に意味のある2つの有限な寄与に収束する因子分解を用いる。

実験結果

リサーチクエスチョン

• RQ1 RPMD速度理論は、正確性と計算効率を保ちつつ、電子的非断熱系に一般化可能か？
• RQ2 対称的および非対称的結合を有する1次元曲線交差モデルにおいて、非断熱RPMD速度は正確な量子速度とどのように比較されるか？
• RQ3 反応座標の質量がゼロに近づく極限においても、Debyeスペクトル密度を有するスピンボソンモデルを正確に記述できるか？
• RQ4 電子状態が結合していない場合に、非断熱RPMD速度は断熱RPMDの極限に正しく還元されるか？
• RQ5 標準的なQTSTが質量要因の発散や未定義のレートにより失敗する系において、RPMD速度は信頼性を持って計算可能か？

主な発見

• 対称な1次元曲線交差モデルでは、100 Kという低温でも古典的レートが桁違いに失敗する中で、RPMDレートは正確な量子レートと2倍以内の一致を示す。
• 非対称な曲線交差モデルでは、対称ケースよりもさらに正確に正確な量子レートを近似することができ、断熱極限における偏差と同等の程度のずれを示す。
• Debyeスペクトル密度を有するスピンボソンモデルにおいて、QTSTが質量要因の発散により未定義であるにもかかわらず、RPMDレートは収束する因子分解により計算可能である。
• 非断熱RPMD速度は、メタフェルド法や半古典的手法を含む、テストされたほぼすべての他の近似手法よりも正確な量子レートに近い。
• 電子結合をオフにした場合に、従来の断熱RPMDに正しく還元されることを確認し、既存理論との整合性が裏付けられた。
• 計算コストは系のサイズに線形に比例するため、電子的非断熱ダイナミクスのスケーラブルな大規模シミュレーションに適している。