온도 적분을 통한 자유 에너지 차이 효율적 계산법

초록

본 논문은 동일한 자유도와 상이한 에너지 지형을 가진 두 시스템 사이의 자유 에너지 차이를 효율적으로 구하기 위한 ‘온도 적분(Temperature Integration)’ 방법을 제안한다. 각 시스템에 대해 병렬 온도 교환(Parallel Tempering) 시뮬레이션을 수행해 두 온도에서의 로그 분배함수(ln Z) 차이를 계산하고, 고온에서 두 시스템이 동일한 위상공간 부피를 갖는 경우 이 차이를 이용해 자유 에너지 차이를 얻는다. 또한 위상공간 부피가 알려진 경우 절대 자유 에너지 계산도 가능하다. 저자는 1차원 고리 위의 강체 막대 모델을 통해 방법의 효율성을 입증한다.

상세 분석

온도 적분 방법은 전통적인 자유 에너지 차이 계산 기법인 Thermodynamic Integration(TI)이나 Free Energy Perturbation(FEP)과는 근본적인 차이를 보인다. TI는 두 시스템 사이의 해밀토니안 변화를 연속적인 λ 파라미터에 따라 샘플링하지만, 복잡한 에너지 지형에서는 충분한 샘플링이 어려워 편향이 발생한다. 반면 온도 적분은 각 시스템을 독립적으로 고온에서 저온으로 연결하는 온도 경로를 이용한다. 병렬 온도 교환(Parallel Tempering) 알고리즘을 적용하면 고온에서의 자유도 탐색이 용이해져, 에너지 장벽을 넘는 샘플링 효율이 크게 향상된다.

핵심 아이디어는 ln Z(T)라는 온도 의존적인 함수의 차이를 직접 측정하는 것이다. 두 온도 T₁, T₂에 대해 ln Z(T₁)−ln Z(T₂)=∫_{β₂}^{β₁}⟨E⟩_β dβ(β=1/k_BT)라는 관계를 이용한다. 여기서 ⟨E⟩_β는 해당 β에서의 평균 에너지이며, Parallel Tempering을 통해 다양한 β값에서 안정적인 ⟨E⟩를 얻을 수 있다. 두 시스템이 동일한 자유도와 동일한 위상공간 부피를 가질 경우, 충분히 높은 온도 T_high에서 ln Z는 거의 동일하게 수렴한다. 따라서 ΔF=F_B−F_A=−k_BT_low