[논문 리뷰] Driven Langevin dynamics: heat, work and pseudo-work

David A. Sivak, John D. Chodera외 1명|arXiv (Cornell University)|2011. 07. 14.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 9인용 수 2

한 줄 요약

이 논문은 유한 시간 간격 랑주뱅 동역학에서 수치 이산화로 인한 오차를 보정하기 위해 '섀도 워크(Shadow Work)'의 개념을 도입한다. 적분기(스텝)를 비평형 과정으로 간주함으로써, 비평형 일의 정리(Work Theorems)를 활용해 시간 간격 의존성 없는 자유 에너지 추정이 가능해지며, 섀도 워크를 통한 보정으로 시간 간격에 의존하는 편향이 제거된다.

ABSTRACT

Common algorithms for computationally simulating Langevin dynamics must discretize the stochastic differential equations of motion. These resulting finite time step integrators necessarily have several practical issues in common: Microscopic reversibility is violated, the sampled stationary distribution differs from the desired equilibrium distribution, and the work accumulated in nonequilibrium simulations is not directly usable in estimators based on nonequilibrium work theorems. Here, we show that even with a time-independent Hamiltonian, finite time step Langevin integrators can be thought of as a driven, nonequilibrium physical process. Once an appropriate work-like quantity is defined -- here called the shadow work -- recently developed nonequilibrium fluctuation theorems can be used to measure or correct for the errors introduced by the use of finite time steps. In particular, we demonstrate that amending estimators based on nonequilibrium work theorems to include this shadow work removes the time step dependent error from estimates of free energies. We also quantify, for the first time, the magnitude of deviations between the sampled stationary distribution and the desired equilibrium distribution for equilibrium Langevin simulations of solvated systems of varying size. While these deviations can be large, they can be eliminated altogether by Metropolization or greatly diminished by small reductions in the time step. Through this connection with driven processes, further developments in nonequilibrium fluctuation theorems can provide additional analytical tools for dealing with errors in finite time step integrators.

연구 동기 및 목표

• 유한 시간 간격 랑주뱅 적분기의 경우 마이크로스코픽 역행성( microscopic reversibility)을 위반하고 올바른 정상 분포를 샘플링하지 못하는 근본적인 문제를 해결하기 위해.
• 용액 시스템에서 샘플링된 정상 분포와 진정한 평형 분포 사이의 편차를 규명하고 정량화하기 위해.
• 비평형 일의 정리에서 시간 간격 오차를 보정할 수 있는 일 유사 양수인 '섀도 워크'를 개발하기 위해.
• 섀도 워크를 추정기(_estimator)에 포함시켜 자유 에너지 계산에서 시간 간격 의존성 오차가 제거됨을 보여주기 위해.
• 메트로폴리제이션 또는 시간 간격 감소가 평형 시뮬레이션에서 발생하는 샘플링 오차를 어떻게 제거하거나 최소화할 수 있는지 탐색하기 위해.

제안 방법

• 유한 시간 간격 이산화로 인해 랑주뱅 동역학에 도입되는 비평형 구동 효과를 정량화하는 새로운 일 유사 양수인 '섀도 워크'를 정의한다.
• 최근 개발된 비평형 열역학 정리( fluctuation theorems)를 섀도 워크에 적용하여 자유 에너지 추정기의 편향을 보정한다.
• 크기가 다른 용액 시스템의 수치 시뮬레이션을 통해 정상 분포의 평형에서의 편차 크기를 정량화한다.
• 표준 적분기와 메트로폴라이제이션된 버전의 성능을 비교하여 보정 프로토콜이 오차를 얼마나 줄이는지 평가한다.
• 시간 간격 이산화에 의해 유도되는 효과적인 비평형 동역학을 유도하고 분석하며, 이를 구동 과정으로 프레임화한다.
• 섀도 워크가 시간 간격 크기와 무관하게 일관된 자유 에너지 추정이 가능하게 함을 입증한다.

실험 결과

연구 질문

• RQ1 유한 시간 간격 이산화로 인해 시간에 독립적인 해밀토니안을 가진 랑주뱅 동역학에서도 비평형 효과가 어떻게 유도되는가?
• RQ2 용액 시스템에서 샘플링된 정상 분포와 진정한 평형 분포 사이의 편차 크기는 어느 정도인가?
• RQ3 비평형 일의 정리에서 시간 간격에 의한 오차를 보정할 수 있는 잘 정의된 일 유사 양수인 '섀도 워크'를 구성할 수 있는가?
• RQ4 섀도 워크는 자유 에너지 추정기에서 시간 간격 의존성 편향을 어느 정도 제거할 수 있는가?
• RQ5 메트로폴리제이션과 시간 간격 감소는 어떤 정도로 시간 간격 유한성으로 인한 샘플링 오차를 완화하는가?

주요 결과

• 유한 시간 간격 랑주뱅 적분기의 경우, 수치 이산화로 인해 시간에 독립적인 해밀토니안이더라도 비평형 구동 과정을 유도한다.
• 평형 시뮬레이션에서 샘플링된 정상 분포는 특히 작은 용액 시스템에서 진정한 평형 분포와 크게 다를 수 있다.
• 크기가 다른 용액 시스템에 대해 이러한 편차의 크기를 정량화한 결과, 보정이 없을 경우 상당한 오차가 있음을 확인했다.
• 메트로폴리제이션은 적분기의 상세 균형을 강제함으로써 샘플링 오차를 완전히 제거한다.
• 시간 간격을 줄이면 샘플링 오차가 크게 감소하지만, 완전히 제거되지는 않는다.
• 비평형 일의 정리에 섀도 워크를 포함시킴으로써 자유 에너지 추정에서 시간 간격 의존성 편향이 성공적으로 제거된다.