[論文レビュー] Tunneling-splittings from path-integral molecular dynamics using a Langevin thermostat

Christophe L. Vaillant, David J. Wales+1 more|arXiv (Cornell University)|Mar 12, 2018

Spectroscopy and Quantum Chemical Studies被引用数 6

ひとこと要約

本稿では、半古典的インスタントン経路に沿った熱力学的統合とランジュバン thermostat を用いた改良型経路積分分子動力学（PIMD）法を提案し、分子系におけるトンネル分裂を計算する。この手法により、水ダイマーのような複雑な系においても、大きな分裂と小さな分裂の両方を精度よく確率的計算可能となり、バリエーショナル基準結果と10–20%以内の一致を達成する。

ABSTRACT

We report an improved method for the calculation of tunneling splittings between degenerate configurations in molecules and clusters using path-integral molecular dynamics (PIMD). Starting from an expression involving a ratio of thermodynamic density matrices at the bottom of the symmetric wells, we use thermodynamic integration with molecular dynamics simulations and a Langevin thermostat to compute the splittings stochastically. The thermodynamic integration is performed by sampling along the semiclassical instanton path, which provides an efficient reaction coordinate as well as being physically well-motivated. This approach allows us to carry out PIMD calculations of the multi-well tunnelling splitting pattern in water dimer, and to refine previous PIMD calculations for one-dimensional models and malonaldehyde. The large (acceptor) splitting in water dimer agrees to within 20% of benchmark variational results, and the smaller splittings are within 10%.

研究の動機と目的

• 縮退した井戸を持つ分子系におけるトンネル分裂をより正確かつ効率的に計算するための手法の開発、特にバリエーショナル法が失敗するような系を対象とする。
• リングポリマーインスタントン（RPI）法などの半古典的近似に起因する限界を克服すること、特に高い非調和性や複雑なノード表面を持つ系で失敗する可能性がある点を改善する。
• 固定ノードDFTMが収束問題により失敗する、水ヘキサマーのような小さな分裂を示す系へのPIMDの適用範囲を拡大すること。
• インスタントン経路に基づく物理的に妥当な反応座標を導入することで、PIMDにおけるサンプリング効率と精度を向上させること。
• 一様な井戸モデル、マロンアルデヒド、および水ダイマーを含むベンチマーク系を用いて手法の妥当性を検証し、分裂の大きさが異なる状況下でも安定性を示す。

提案手法

• トンネル分裂は、対称配置における熱力学的密度行列比を用いて計算され、量子統計力学から導かれた式 Iν(β) = tanh[Δν/2 (β − β̄ν)] を用いる。
• 熱力学的統合は、サンプリングに効率的で物理的に妥当な反応座標としての半古典的インスタントン経路に沿って実行される。
• PIMDシミュレーションにおけるサンプリング効率を向上させるために、ランジュバル thermostat が導入されている。
• トンネル分裂 Δν は、式 (4) の解析的形を用いて、2つの逆温度 β(1) と β(2) における Iν(β) のフィッティングにより抽出され、エネルギー準位の直接的収束を回避する。
• 多井戸系では、射影演算子とキャラクターの直交性を用いて、η±ν(β) 関数を通じて個々の振動準位の寄与を分離する。
• 固定ノードDFTMとは異なり、ノード表面の事前知識が不要であり、極めて小さな分裂の計算が可能となる。

実験結果

リサーチクエスチョン

• RQ1 ランジュバル thermostat とインスタントン経路に沿った熱力学的統合を備えたPIMDベースの手法は、大きな分裂と小さな分裂を示す両方の系において、トンネル分裂を正確に計算できるか？
• RQ2 特に受容体のトンネル運動に注目した場合、水ダイマーにおいて、本手法はバリエーショナル法およびRPI法と比較して、精度と効率の面で優れているか？
• RQ3 インスタントン経路を反応座標として用いることで、従来のPIMD実装と比較して、サンプリング効率と収束性がどの程度向上するか？
• RQ4 ノード表面の複雑さのため、固定ノードDFTMが到達できない極めて小さな分裂（例：∼10−5 cm−1）を、本手法が信頼性を持って計算できるか？
• RQ5 一様な井戸モデルおよびマロンアルデヒドにおいて、本手法はベンチマーク結果をどの程度再現できるか？また、その結果から大規模系への応用にどのような示唆が得られるか？

主な発見

• 本手法は、水ダイマーにおける大きな受容体トンネル分裂を、バリエーショナル基準結果と比較して20%の精度で計算した。
• 水ダイマーにおけるより小さな分裂も、バリエーショナル基準値と10%以内の一致を示し、小さな分裂に対する高い精度を実証した。
• ランジュバル thermostat の導入により、従来の実装と比較して、PIMDシミュレーションにおけるサンプリング効率と収束性が顕著に向上した。
• インスタントン経路に沿った熱力学的統合は、物理的に妥当で効率的な反応座標を提供し、分裂計算の安定性と精度を向上させた。
• 本手法は、水ダイマーにおける全トンネル分裂パターン（大きな分裂と小さな分裂を含む）を成功裏に計算でき、DFTMおよびRPI法が困難とする系に対しても有効であることを示した。
• 一様な井戸モデルおよびマロンアルデヒドを用いた検証では、一貫性があり正確な結果が得られ、ノード表面が複雑な系や極めて小さな分裂を示す系へのPIMDの適用範囲を拡大した。