[論文レビュー] Some Further Results for the Stationary Points and Dynamics of Supercooled Liquids

David J. Wales, Jonathan P. K. Doye|arXiv (Cornell University)|Sep 2, 2003

Advanced Thermodynamics and Statistical Mechanics参考文献 4被引用数 98

ひとこと要約

本稿では、超冷却液体のポテンシャルエネルギー面を、局所的エネルギー最小値だけでなく、すべての停留点に対して、ニュートン・ラプソン法と固有ベクトルフォローアップ法を用いた手法により、流域に分割する方法を提示する。弱い相互作用を示す部分系において、ヘッセ行列のインデックスごとの停留点数の解析的表現を導出し、低インデックスのヘッセ行列のサドル点（特にインデックス-1）がより局所的で空間的に独立した変位を示すことを明らかにした。これは、動的モデルにおける有効な遷移状態としての利用を支持する。

ABSTRACT

We present some new theoretical and computational results for the stationary points of bulk systems. First we demonstrate how the potential energy surface can be partitioned into catchment basins associated with every stationary point using a combination of Newton-Raphson and eigenvector-following techniques. Numerical results are presented for a 256-atom supercell representation of a binary Lennard-Jones system. We then derive analytical formulae for the number of stationary points as a function of both system size and the Hessian index, using a framework based upon weakly interacting subsystems. This analysis reveals a simple relation between the total number of stationary points, the number of local minima, and the number of transition states connected on average to each minimum. Finally we calculate two measures of localisation for the displacements corresponding to Hessian eigenvectors in samples of stationary points obtained from the Newton-Raphson-based geometry optimisation scheme. Systematic differences are found between the properties of eigenvectors corresponding to positive and negative Hessian eigenvalues, and localised character is most pronounced for stationary points with low values of the Hessian index.

研究の動機と目的

• すべての停留点（局所的エネルギー最小値に限らない）に関連する流域に、ポテンシャルエネルギー面を体系的かつ安定的に分割するための手法を開発すること。
• 弱い相互作用を示す部分系において、ヘッセ行列インデックスごとの停留点数の解析的表現を導出し、従来の統計力学フレームワークを拡張すること。
• ヘッセ行列固有ベクトルの空間的局在性と独立性を調査し、正の固有値と負の固有値を持つモードの違いを明確にすること。
• エネルギー最小値と真の遷移状態（インデックス-1）に基づく従来の動的理論が、超冷却液体において高インデックスのサドル点が支配的である場合にも有効であるかを評価すること。

提案手法

• ニュートン・ラプソンに基づく幾何最適化スキームと固有ベクトルフォローアップを組み合わせ、ポテンシャルエネルギー面上のすべての停留点を特定・分類する。
• 各停留点から勾配降下路を追跡することで流域の吸引域を特定し、流域分割を実施する。
• 弱い相互作用を示す部分系に基づく組み合わせ的モデルを構築し、系のサイズとヘッセ行列インデックスの関数として、停留点数を解析的に計算する。
• 空間的広がり（L）と空間的独立性（N）の2つの指標を用いて、ヘッセ行列固有ベクトルの局在性を定量化し、正の固有値と負の固有値を持つ固有ベクトルを比較する。
• 256原子のバイナリ・レナード・ジョーンズ系のシミュレーションから得られた停留点を用い、解析的モデルの妥当性を検証するとともに、固有ベクトルの性質を評価する。
• 低インデックスサドル点に注目し、異なるヘッセ行列インデックスの固有ベクトルの空間的特性を比較分析する。

実験結果

リサーチクエスチョン

• RQ1 すべての停留点（高インデックスサドル点を含む）に対して、体系的かつ数値的に安定した方法でポテンシャルエネルギー面を流域に分割することは可能か？
• RQ2 弱い相互作用系において、停留点総数、局所的エネルギー最小値の数、および最小値1つあたりの遷移状態平均数との間には、解析的な関係が存在するか？
• RQ3 正の固有値モードと負の固有値モードの間で、ヘッセ行列固有ベクトルの空間的局在性と独立性にはどのような違いがあるか？
• RQ4 負の固有値に対応する固有ベクトル（サドル点）は、正の固有値に対応するものよりもより局在的かつ空間的に独立しており、特にヘッセ行列の低インデックスでその傾向が顕著になるか？
• RQ5 超冷却液体において高インデックスサドル点が優勢である場合でも、エネルギー最小値とインデックス-1の真の遷移状態に基づく従来の動的理論は有効に保たれるか？

主な発見

• ニュートン・ラプソン法と固有ベクトルフォローアップ法を用いた手法により、すべての停留点に対してポテンシャルエネルギー面の流域分割が成功し、局所的エネルギー最小値を超えた体系的分析が可能となった。
• 弱い相互作用部分系に対する解析的モデルにより、停留点総数、局所的エネルギー最小値の数、および最小値1つあたりの遷移状態平均数との間の単純な関係が予測された。
• 負のヘッセ行列固有値に対応する固有ベクトルは、正の固有値に対応するものと比べて顕著に高い空間的局在性と独立性を示した。
• 固有ベクトルの局在性と空間的独立性は、ヘッセ行列のインデックスが低い（特に I=1）停留点で最も顕著に現れた。
• 本結果は、低インデックスサドル点が、部分系からの真の遷移状態の組み合わせとして効果的に記述可能であるという見解を支持しており、Shell らの先行研究とも整合的である。
• 超冷却液体において高インデックスサドル点が顕著に存在するにもかかわらず、動的挙動は依然としてエネルギー最小値と真の遷移状態（インデックス-1）に基づく理論で適切に記述可能である。これは、最低の障壁が依然としてインデックス-1サドル点を介して媒介されているからである。