[논문 리뷰] Ergodicity of a Time-Reversibly Thermostated Harmonic Oscillator and the 2014 Ian Snook Prize

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

Advanced Physical and Chemical Molecular Interactions인용 수 3

한 줄 요약

이 논문은 수치 시뮬레이션과 위상공간 분석을 통해 세 가지 동역학 모델—Marta-Klein-Tuckerman(MKT), Hoover-Holian(HH), Patra-Bhattacharya(PB)—을 사용하여 시간 역전 가능하게 온도 조절된 조화 진동자의 에르고딕성(ergodicity)을 조사한다. MKT 및 HH 모델은 캐논리컬 가우시안 분포와 일치하는 에르고딕 행동를 보이며, PB 모델은 토러스 구조와 밀도의 편차를 보이며 비에르고딕 특성을 나타내어 이전에 MKT의 비에르고딕성을 주장한 바를 도전하고, 이러한 시스템에서의 에르고딕성에 대한 명확한 증명을 위한 2014년 Ian Snook 상을 제안한다.

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년 Ian Snook 상을 제안하는 것.

제안 방법

• 1000만 개의 스텝 동안 시간 간격 0.001로 4차 룬게-쿠타 방법을 사용한 수치적 적분을 통해 MKT, HH, PB 모델의 궤도를 시뮬레이션.
• 시간 평균을 통한 속도 모멘트(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 모델이 캐논리컬 가우시안 분포와 일치하는 시간 평균 속도 모멘트(p², p⁴, p⁶)를 생성하는가?
• 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)로 산출하여 캐논리컬 가우시안 분포와 네 자리 유효숫자까지 정확히 일치시킨다.
• PB 모델은 모멘트로 (1.000, 3.800, 23.50)를 산출하여 가우시안 행동에서의 상당한 편차를 보이며 비에르고딕성을 나타낸다.
• MKT 및 HH의 위상공간 투영은 원점 근처에 약간의 구멍이 있는 혼돈의 바다를 보이며, PB는 q와 p 사이에 강한 양의 상관관계와 뚜렷한 토러스 구조를 보인다.
• 반경이 감소하는 (q,p) 원형 영역(최소 1/16 원래 면적) 내의 밀도 측정 결과, MKT 및 HH에서는 변화가 없어 원점 근처의 균일한 샘플링을 나타낸다.
• PB 모델의 밀도는 MKT/HH의 약 절반으로, 약 50%의 위상공간 측도가 준주기적이고 혼돈적 샘플링에서 제외됨을 시사한다.
• PB 모델은 (q,p,ζ) 공간에서 명확한 토러스 구조를 보이며, 내재된 주기 궤도와 혼돈의 바다 내 큰 빈공간이 존재함을 확인하여 비에르고딕 행동를 확인한다.