[論文レビュー] Fast numerical evaluation of time-derivative nonadiabatic couplings for mixed quantum-classical methods

Ilya G. Ryabinkin, Jayashree Nagesh+1 more|arXiv (Cornell University)|Aug 6, 2015

Spectroscopy and Quantum Chemical Studies被引用数 6

ひとこと要約

本論文は、スレーター行列式の代わりに分子軌道を用いて計算を再定式化することにより、混合量子古典的ダイナミクスにおける時間微分非断熱結合（TDNAC）を数値的に効率的に計算する手法を提示する。これにより、計算コストの高い重なり行列式の評価が不要となり、ベンチマーク系では最大1,372倍の高速化が達成された。この結果、TDNACの評価が非断熱ダイナミクスシミュレーションにおけるボトルネックから取り除かれ、実用的な計算が可能になった。

ABSTRACT

We have developed a numerical differentiation scheme which eliminates evaluation of overlap determinants in calculating the time-derivative non-adiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater-determinant representation. The central idea of our approach is, first, to reduce the analytic time derivatives of Slater determinants to time derivatives of molecular orbitals, and then to apply a finite-difference formula. Benchmark calculations prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calculation step for midsize molecules.

研究の動機と目的

• TDNAC計算中にスレーター行列式の重なり行列を評価することによって生じる混合量子古典的非断熱ダイナミクスにおける計算ボトルネックを解消すること。
• TDNACの数値微分スキームにおいて、重なり行列式の明示的評価を回避すること。
• 多数電子波動関数の重なりに依存する従来の有限差分法の代替として、数値的に安定かつ効率的な手法を開発すること。
• TDNAC評価のコストを著しく削減することで、中規模および大規模分子における非断熱ダイナミクスの実用的シミュレーションを可能にすること。
• 異なる分子系および基底関数集合において、本手法の精度と性能を検証すること。

提案手法

• 本手法は、電子波動関数をスレーター行列式ではなく分子軌道で表現することにより、時間微分非断熱結合（TDNAC）の計算を再定式化する。
• 一階前進差分および二階中央差分の有限差分式を、全多数電子波動関数の代わりに分子軌道の時間微分に直接適用する。
• 中心的イノベーションは、スレーター行列式の解析的時間微分を、その構成分子軌道の時間微分に還元することにより、非正規直交なスレーター行列式間の重なり行列式の計算を回避することにある。
• 本手法は、軌道係数の行列-行列乗算を活用し、スクリーニング手順を不要とし、分子軌道の直交性を活かして計算を簡略化する。
• 標準的な量子化学ソフトウェア（Gaussian）を用いて実装されており、スピン制限形式と互換性があり、αスピンおよびβスピンチャネルの重複計算を回避することで計算コストを削減する。
• TDNACは連鎖則を用いて τ_KJ = Ṙ · d_KJ として計算され、d_KJは軌道微分から得られ、最終的な結合は軌道レベルでの有限差分近似により得られる。

実験結果

リサーチクエスチョン

• RQ1 スレーター行列式の重なり計算を明示的に回避することで、時間微分非断熱結合の評価が高速化可能か？
• RQ2 TDNACの軌道ベースの有限差分スキームは、行列式ベースの手法と比較して十分な数値的精度を維持できるか？
• RQ3 提案手法が混合量子古典的非断熱ダイナミクスシミュレーションの計算コストをどの程度低減できるか？
• RQ4 システムサイズおよび基底関数サイズの増大に伴い、本手法の性能はどのようにスケーリングするか？
• RQ5 本手法は、CISなどの標準的電子構造手法を用いた中規模有機分子に効果的に適用可能か？

主な発見

• 提案された軌道ベース手法により、C18H14O や C25H18 といった中規模有機分子において、TDNAC評価で最大1,372倍の高速化が達成された。特にSTO-3Gなどの小規模基底関数において顕著な向上が見られた。
• 6-31G**基底関数を用いたC18H14Oの50 fs FSSH軌道において、合計シミュレーション時間は元のNewton-Xの99時間から本手法では45時間に短縮された。TDNAC計算時間は54時間からわずか12分にまで削減された。
• 本手法は高い数値的精度を維持しており、時間刻みが0.05 fs程度でもTDNACの誤差が1×10⁻⁷未満に抑えられ、二階中央差分スキームは一階前進差分スキームに比べて優れた収束性を示した。
• 基底関数サイズが増大しても、軌道ベース手法は従来の行列式ベース手法を2桁以上上回る性能を示し、特に占有軌道数が全基底サイズと同等に近い系で顕著であった。
• 本手法は異なる電子状態および分子系においても安定しており、C18H14OおよびC25H18の両者において、複数の基底関数セットで一貫した性能向上が確認された。
• 重なり行列式の評価を排除することで、スピン制限形式の効率的利用が可能となり、従来のスピン非制限実装に伴う二重計算コストを回避できた。