[論文レビュー] Ergodicity of a Time-Reversibly Thermostated Harmonic Oscillator and the 2014 Ian Snook Prize

Wm.G. Hoover, C.G. Hoover|arXiv (Cornell University)|Jan 30, 2020

Advanced Physical and Chemical Molecular Interactions被引用数 3

ひとこと要約

本稿では、数値的シミュレーションと位相空間解析を通じて、時間反転可逆にトロスレートされた調和振動子のエルゴード性を、Marta-Klein-Tuckerman (MKT)、Hoover-Holian (HH)、Patra-Bhattacharya (PB) の3つの力学的モデルを用いて調査している。MKTおよびHHモデルは、カノニカルなガウス分布と整合するエルゴード的挙動を示すが、PBモデルはトーラス構造や密度のずれといった非エルゴード的特徴を示し、MKTの非エルゴード性に関する以前の主張に反論するとともに、これらの系におけるエルゴード性の決定的証明を求める2014年イアン・スヌーク賞を提供する。

ABSTRACT

Shuichi Nos\'e opened up a new world of atomistic simulation in 1984. He formulated a Hamiltonian tailored to generate Gibbs' canonical distribution dynamically. This clever idea bridged the gap between microcanonical molecular dynamics and canonical statistical mechanics. Until then the canonical distribution was explored with Monte Carlo sampling. Nos\'e's dynamical Hamiltonian bridge requires the "ergodic" support of a space-filling structure in order to reproduce the entire distribution. For sufficiently small systems, such as the harmonic oscillator, Nos\'e's dynamical approach failed to agree with Gibbs' sampling and instead showed a complex structure, partitioned into a chaotic sea, islands, and chains of islands, that is familiar textbook fare from investigations of Hamiltonian chaos. In trying to enhance small-system ergodicity several more complicated "thermostated" equations of motion were developed. All were consistent with the canonical Gaussian distribution for the oscillator coordinate and momentum. The ergodicity of the various approaches has undergone several investigations, with somewhat inconclusive ( contradictory ) results. Here we illustrate several ways to test ergodicity and challenge the reader to find even more convincing algorithms or an entirely new approach to this problem.

研究の動機と目的

• 時間反転可逆にトロスレートされた調和振動子のエルゴード性を、異なる力学的モデルのもとで検証すること。
• PatraとBhattacharyaがMKTが非エルゴード的であると主張したことに鑋いて、MKTおよびHHモデルに関する矛盾する主張を解消すること。
• 時間平均、投影、密度測定、周期軌道解析を通じた数値的および幾何的証拠を提供し、これらのモデルが全カノニカル分布を正しくサンプリングしているかどうかを検証すること。
• 研究コミュニティに、より厳密または新しいアルゴリズムを開発してエルゴード性を証明するよう挑戦し、2014年イアン・スヌーク賞をインcentiveとして提供すること。

提案手法

• MKT、HH、PBモデルの軌道をシミュレートするために、1000万ステップ、ステップサイズ0.001の4次ルンゲ・クッタ法を用いた数値積分。
• 時間平均による速度モーメント（p²、p⁴、p⁶）の計算を行い、カノニカルなガウス分布の期待値（1、3、15）と比較。
• 位相空間の投影を(q,p)、(q,ζ)、(q,ξ)、(p,ζ)、(p,ξ)、(ζ,ξ)平面に施し、分布構造の可視化と穴や対称性の有無の検出。
• 10億ステップを用いて、(q,p)円の面積を0.01πから0.01π/16まで縮小させた範囲で密度測定を行い、原点付近の均一性を評価。
• 初期条件を原点付近に設定し、特にPBモデルにおいて周期軌道を特定。これによりネストされたトーラス構造が明らかになった。
• ランダムまたはグリッドベースの初期条件を用いたアンサンブルテストにより、二峰性分布の有無を検出。非エルゴード性の兆候とみなされる。

実験結果

リサーチクエスチョン

• RQ1 MKT、HH、PBモデルは、時間平均された速度モーメントがカノニカルなガウス分布と整合的か？
• RQ2 これらのモデルの位相空間投影は、穴や非対称性がなく、ガウス的構造を示すか？ これはエルゴード性を示唆する。
• RQ3 各モデルについて、4次元位相空間における原点(q,p,ζ,ξ) = (0,0,0,0)付近の確率密度は一様か？
• RQ4 PB方程式は周期軌道やトーラス的構造を示し、非エルゴード的挙動を示唆するか？
• RQ5 MKTまたはHHモデルについて、エルゴード性の決定的数学的証明を構築できるか？ これにより、矛盾する主張を解消できる。

主な発見

• MKTおよびHHモデルは、速度モーメント（p²、p⁴、p⁶）をそれぞれ（1.000、3.000、15.00）として、4桁の有効数字でカノニカルなガウス分布と一致する。
• PBモデルはモーメントとして（1.000、3.800、23.50）を示し、ガウス的挙動から顕著に逸脱しており、非エルゴード的であると示唆される。
• MKTおよびHHモデルの位相空間投影は、原点付近にわずかな穴を伴う混沌とした領域を示すが、PBモデルはqとpの間に強い正の相関と明確なトーラス的構造を示す。
• 半径を小さくした(q,p)円内の密度測定では、MKTおよびHHモデルにおいては顕著な変化が認められず、原点付近での均一なサンプリングが示された。
• PBモデルの密度はMKT/HHの約半分であり、その位相空間測度のおよそ50％が準周期的で、混沌としたサンプリングから除外されていると示唆される。
• PBモデルは(q,p,ζ)空間において明確なトーラス的構造を示し、ネストされた周期軌道と混沌とした領域内の大きな空洞が確認され、非エルゴード的挙動が確認された。