A General Non-Markovian Electronic Friction and Langevin Dynamics Framework for Nonequilibrium Systems

技術要約：非平衡系のための一般的な非マルコフ的電子摩擦およびランジュバン動力学フレームワーク

問題提起
金属表面と相互作用する分子のビブロン動力学をシミュレーションすることは、分子散乱、表面反応、非平衡電荷輸送などの現象を理解するために不可欠である。しかし、これらの系はしばしば、動力学が高度に非マルコフ的かつ非摂動的となる強結合領域で動作する。階層型方程式（HEOM）やマルチレイヤー多配置時間依存ハートリー法（ML-MCTDH）のような完全量子手法は正確なベンチマークを提供するが、多くの振動モードや大きな励起を持つ系に対しては計算コストが極めて高い。

その結果、電子摩擦およびランジュバン動力学（EFLD）に特化した混合量子・古典アプローチが広く用いられている。これらの手法は通常、電子ノイズをホワイトノイズとして扱い、電子摩擦が瞬時に減衰すると仮定するマルコフ近似に依存している。この近似は、特定の極限（例：重い原子核や強結合）では成功しているが、多体相互作用やバイアス駆動による非弾性遷移が存在する系では、マルコフ近似が不適切であるという証拠が増えている。具体的には、マルコフ近似は負の電子摩擦といった物理的に不自然な不安定性を引き起こす可能性があり、非平衡領域における電子のバックアクションの完全な動的構造を捉えることができない。非マルコフ的EFLDを拡張する際の主要な課題は、座標依存性を持つ量子カラーノイズ、特に揺動散逸定理（FDT）が非平衡下で破れている場合のサンプリングの困難さである。

手法
著者らは、非マルコフ的電子摩擦およびランジュバン動力学（NM-EFLD）を用いた非平衡ビブロン動力学をシミュレートするための一般的なフレームワークを提示している。本手法の核心は、メモリカーネルを持つ一般化ランジュバン方程式を、一連の結合されたマルコフ的確率微分方程式へと変換するマルコフ埋め込みスキームである。

主要な手法構成要素は以下の通りである：

1. 一般化ランジュバン方程式： 振動動力学は、平均断熱力、非マルコフ的摩擦項、およびカラーノイズを持つガウス型確率力によって記述される一般化ランジュバン方程式に従う。
2. 独立極展開： 座標と時間の因子分解を仮定することが多い従来の埋め込みスキームとは異なり、本フレームワークでは、電子摩擦カーネル $\gamma(x, t)$ と確率力の相関関数 $D(x, t)$ の両方に対して独立した指数関数展開を行う。
   • $\gamma(x, t) = \text{Re} \left\{ \sum a_j(x) e^{-b_j(x)t} \right\}$
   • $D(x, t) = \text{Re} \left\{ \sum c_j(x) e^{-d_j(x)t} \right\}$
     これにより、量子FDTが破れており、$\gamma$ と $D$ が単純な関係にない非平衡シナリオを扱うことが可能となる。
3. 補助変数： メモリカーネルは、摩擦の履歴を表す決定論的な補助変数 ($u_j$) と、カラーノイズを表す確率的補助変数 ($w_j$) にマッピングされる。これらの変数は、電子力の完全な座標依存性を保持しながら、結合されたオルンシュタイン＝ウーレンベック過程に従って進化する。
4. 伝播アルゴリズム： 著者らは、結合された運動方程式を積分するために、更新されたABOBAアルゴリズム（確率プロパゲーターの2次対称分解）を採用している。このアルゴリズムには、物理的な振動自由度とともに補助変数の時間発展を扱うための、特定の演算子分割が含まれている。

主な貢献

• 一般的な非平衡フレームワーク： 電子力の完全な座標依存性を保持し、古典的なFDTを必要とせずに量子カラーノイズを扱う手法を導入しており、バイアスのかかった非平衡輸送への適用が可能である。
• 独立展開： 摩擦カーネルとノイズ相関関数の展開を分離することで、非平衡下ではこれらの量に異なる数の指数項が必要であり、単純な平衡関係によって結びつけられないことを認めている。
• 安定性と精度： 本手法は、従来のマルコフ的EFLD（M-EFLD）が失敗するパラメータ領域、具体的にはM-EFLDが物理的に不自然な負の摩擦を予測して不安定になる領域においても、数値的に安定していることが示されている。

結果
本手法は、振動結合したドナー・アクセプター分子接合モデルを用いて検証された。結果は、数値的に厳密なHEOMシミュレーションと比較された：

• 平衡状態の検証： ゼロバイアスにおいて、NM-EFLDアプローチは量子熱期待値および確率力の正確な相関関数を正しく再現しており、補助変数表現の妥当性を確認している。
• 非平衡性能：
  • 正のエネルギー分裂 ($\Delta > 0$) を持つモデルでは、NM-EFLDアプローチは高バイアス電圧における定常状態の振動励起を正確に予測し、HEOMのベンチマークと密接に一致した。対照的に、M-EFLDは励起を過小評価した。
  • 負のエネルギー分裂 ($\Delta < 0$) を持つモデルでは、M-EFLDアプローチは中程度のバイアス電圧 ($\Phi \approx 0.6$ V) で不安定になり、厳密な量子結果よりも数桁大きい振動励起を予測した。NM-EFLDアプローチは全バイアス範囲にわたって安定しており、HEOMの結果を定量的に再現した。
• 不安定性のメカニズム： 電子力のスペクトル解析により、ゼロ周波数（マルコフ的）の摩擦係数が負（駆動を示唆）であっても、完全な周波数依存の摩擦は正であり支配的であることが明らかになった。マルコフ近似は、有限周波数成分を無視することでエネルギー交換を定性的に誤認し、不安定性を招いている。NM-EFLDフレームワークは、これらの有限周波数のメモリ効果を捉えることで、動力学を正しく安定化させている。

意義と主張
本論文は、非マルコフ的な電子力が、特にマルコフ近似が失敗する領域において、非平衡振動動力学を決定づける役割を果たすことを主張している。著者らは、これらの系における電子緩和は、分子と金属の結合だけでなく、電子のエネルギー分裂によっても制御され、複数のタイムスケールにわたるメモリ効果を生み出すと論じている。

本研究の意義は、電流駆動型の分子輸送、振動エネルギー転移、および表面散乱を研究するための、大規模なシミュレーションに組み込める実用的かつ一般的なフレームワークを提供した点にある。独立した極展開を用いたマルコフ埋め込みを利用することで、本手法は、計算コストの高い完全量子手法と、広く用いられているがしばしば不正確なマルコフ的EFLDとの間の溝を埋め、複雑な分子・金属界面を非平衡下で研究するための堅牢なツールを提供している。