[논문 리뷰] Open Systems out of Equilibrium: Theory and Simulation

Roya Ebrahimi Viand, Felix Höfling외 2명|arXiv (Cornell University)|2020. 05. 15.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 51인용 수 5

한 줄 요약

이 논문은 비평형 상태의 개방계를 위한 이론적 프레임워크인 Bergmann-Lebowitz(BL) 이론과 일치시키는 방식으로 Adaptive Resolution Simulation(AdResS) 방법의 비평형 확장을 제안한다. 열기울기가 작용하는 개방 경계 조건 하에서 Lennard-Jones 유체를 시뮬레이션함으로써, 밀도, 온도, 열류의 정적인 프로파일이 전자기적 원자모형 시뮬레이션과 뛰어난 일치를 보이며, AdResS가 비평형 분자의 역학을 효율적이고 신뢰성 있게 시뮬레이션할 수 있음을 검증한다.

ABSTRACT

We consider the theoretical model of Bergmann and Lebowitz for open systems out of equilibrium and translate its principles in the adaptive resolution molecular dynamics technique (AdResS). We simulate Lennard-Jones fluids with open boundaries in a thermal gradient and find excellent agreement of the stationary responses with results obtained from the simulation of a larger, locally forced closed system. The encouraging results pave the way for a computational treatment of open systems far from equilibrium framed in a well-established theoretical model that avoids possible numerical artifacts and physical misinterpretations.

연구 동기 및 목표

• 비평형 상태에서 물리적으로 일관되며 계산적으로 효율적인 개방 분자계 시뮬레이션 방법을 개발하는 것.
• 개방계 이론 모델(Bergmann-Lebowitz)과 실용적 분자역학 시뮬레이션 간 격차를 메우는 것.
• 대규모 시스템의 전자기적 원자모형 시뮬레이션과의 비교를 통해 AdResS 방법의 비평형 조건 하에서의 타당성을 검증하는 것.
• 소프트 물질 및 나노스케일계에서의 복잡한 비평형 현상을 계산 비용을 줄여 연구할 수 있도록 하는 것.
• 입자-연속체 결합을 포함한 복잡한 환경에서 질량 및 에너지 유량을 시뮬레이션할 수 있는 이식 가능한 프레임워크를 제공하는 것.

제안 방법

• 열기울기와 저장소 기반 설정(트레이서, 열조절장치, 열역학적 힘 포함)을 도입하여 AdResS 방법을 비평형 조건으로 확장한다.
• 계를 세 영역으로 나눈다: 원자모형(AT), 열조절장치와 한 입자 힘을 갖는 델타(Δ), 그리고 저장소를 나타내는 트레이서(TR).
• 열역학적 힘 F_th(r)는 평형 AdResS 시뮬레이션에서 유도되며, AT 영역 내에서 적절한 통계역학을 보장한다.
• 비평형 평균은 D-NEMD 접근법을 통해 계산되며, 비평형 시뮬레이션에 사용되는 초기 조건으로 상관관계가 없는 평형 궤적을 사용한다. 이때 열기울기가 임포즈된다.
• 열류는 제어체적 내에서 입자 에너지, 응력 텐서 및 속도를 포함하는 표준 표현식을 통해 계산된다.
• 수치적 안정성과 장시간 통합을 위해 GPU 가속을 활용한 해밀토니안 역학을 활용한다.

실험 결과

연구 질문

• RQ1 Bergmann-Lebowitz 이론 모델과 일치시키는 AdResS 방법이 비평형 상태의 개방계를 정확하게 시뮬레이션할 수 있는가?
• RQ2 AdResS 설정이 전자기적 원자모형 시뮬레이션과 비교해 밀도, 온도, 열류의 정적인 프로파일을 얼마나 잘 재현하는가?
• RQ3 열역학적 힘과 열조절장치가 비평형 AdResS 설정에서 물리적 일관성을 유지하는 데 어떤 역할을 하는가?
• RQ4 이 방법은 초임계 유체와 같이 높은 압축성과 고압을 갖는 시스템을 신뢰성 있게 시뮬레이션할 수 있는가?
• RQ5 AdResS-AdResS 프레임워크는 질량 및 에너지 유량을 포함한 다른 비평형 현상으로 일반화될 수 있는 정도는 어느 정도인가?

주요 결과

• 열기울기가 작용하는 개방 경계 조건 하에서 AdResS 설정은 전자기적 원자모형 시뮬레이션과 정밀한 정량적 일치를 보이며, AT 영역 내에서 밀도, 온도, 열류의 정적인 프로파일을 재현한다.
• Δ 영역의 온도 프로파일은 Δ/TR 경계 근처에서 약간의 편차를 보이지만, 이는 Δ 영역 내에서 빠르게 회복되며 효과적인 열평형을 나타낸다.
• ρ̄ = 0.3 σ⁻³ 및 T̄ = 1.5 k_B/ε, ΔT = 0.2 k_B/ε 조건에서 초임계 Lennard-Jones 유체에 대해 이 방법은 기준 시뮬레이션과 높은 일치를 유지하며, 다양한 유체 상태에서의 강건성을 확인한다.
• 유도된 열역학적 힘 F_th(r)는 개방계 조건을 충족함을 보여주는 보조 자료에서 확인되었으며, 이론적 일관성을 보장한다.
• 이 방법은 잡음이 많은 고정 크기의 NEMD 시뮬레이션에서 흔히 발생하는 잡음 없이 비평형 현상을 정확하게 시뮬레이션할 수 있으며, 계산 비용을 줄일 수 있다.
• 이 프레임워크는 향후 입자-연속체 결합 및 복잡한 분자 환경에서의 응용에 적합한 이식 가능하고 일반적인 비평형 개방계 모델을 제공한다.