[논문 리뷰] The hydration free energy of water

Gerhard Hummer, Lawrence R. Pratt외 1명|ArXiv.org|1995. 07. 20.

Fluid dynamics and aerodynamics studies인용 수 14

한 줄 요약

이 연구는 변동 이론과 다수의 시뮬레이션 방법을 사용하여 전하의 함수로써 물의 수화 자유 에너지를 계산하며, 연속체 모델이 예측하는 2차 함수 형태와의 상당한 이탈을 입증한다. 저자들은 겹치는 히스토그램, 수락 비율, 기하 평균 보간 방법 간 일관된 결과를 도출하여 Rick과 Berne의 열역학적 통합 결과를 확인하면서도, 단순 선형 모델을 무효화하는 비정규 분포의 전기적 상호작용 변동을 드러낸다.

ABSTRACT

We study the chemical potential of water as a function of charge based on perturbation theory. By calculating the electrostatic-energy fluctuations of two states (fully charged and uncharged) we are able to determine accurate values for the dependence of the chemical potential on charge. We find identical results for the chemical-potential difference of fully charged and uncharged water from overlapping-histogram and acceptance-ratio methods and by smoothly connecting the curves of direct exponential averages. Our results agree with those of Rick and Berne (J. Am. Chem. Soc., 1994, 116, 3949) with respect to both the chemical-potential difference and its dependence on the charge coupling parameter. We observe significant deviations from simple Gaussian-fluctuation statistics. The dependence on the coupling parameter is not quadratic, as would be inferred from linear continuum models of electrostatics.

연구 동기 및 목표

• 균형 시뮬레이션을 사용하여 전하의 함수로써 물의 화학적 포텐셜에 기여하는 전기적 기여도를 정확히 결정하는 것.
• 두 번째 순서의 변동 이론과 정규 분포 통계가 물의 수화 자유 에너지를 모델링하는 데 있어 타당한지 테스트하는 것.
• 물의 전기적 상호작용 시뮬레이션에서 유한 체적 효과와 자기상호작용 보정의 영향을 평가하는 것.
• 동일한 시스템에 대해 겹치는 히스토그램, 수락 비율, 기하 평균 보간과 같은 다수의 자유 에너지 추정 기법을 비교하는 것.
• Rick과 Berne(1994)의 기존 열역학적 통합 데이터와의 일치를 검증하는 것.

제안 방법

• 전하가 있는 물과 없는 물 간의 화학적 포텐셜 차이를 전기적 에너지 차이의 지수 평균과 연결하기 위해 변동 이론을 사용한다.
• 고정된 전하 상태에서의 구성 평균을 통해 초과 화학적 포텐셜을 계산하기 위해 전위 분포 정리(potential-distribution theorem)를 적용한다.
• 장거리 전기적 상호작용에 대해 Ewald 합산과 일반화된 반응성장(GRF) 방법을 사용하며, 자기상호작용 보정을 일관되게 포함한다.
• 독립된 시뮬레이션에서 두 전하 상태에서의 자유 에너지 차이를 추정하기 위해 Bennett의 겹치는 히스토그램 및 수락 비율 방법을 사용한다.
• 전하 결합 매개변수 λ를 통해 직접적인 지수 평균을 부드럽게 연결하는 기하 방법을 사용한다.
• 유한 체적 효과를 평가하고 수렴성을 검증하기 위해 크기가 다른 시스템(N = 64 및 N = 256)에서 시뮬레이션을 수행한다.

실험 결과

연구 질문

• RQ1 물의 수화 자유 에너지가 Born 모델이 예측한 것처럼 전하 결합 매개변수 λ에 대해 2차 함수 형태를 띠는가?
• RQ2 겹치는 히스토그램 및 수락 비율과 같은 변동 기반 방법이 전하가 있는 물과 없는 물 간의 화학적 포텐셜 차이를 추정하는 데 얼마나 정확한가?
• RQ3 유한 체적 효과와 자기상호작용 항목의 생략이 물 시뮬레이션에서 자유 에너지 추정에 얼마나 큰 영향을 미치는가?
• RQ4 자기상호작용 보정이 포함된 경우, Ewald 및 GRF 전기적 모델 간에 일관된 결과를 도출할 수 있는가?
• RQ5 전기적 에너지의 비정규 분포 변동이 수용액계에서 두 번째 순서의 변동 이론의 타당성에 어떤 영향을 미치는가?

주요 결과

• 완전히 전하가 있는 물과 전하가 없는 물 간의 화학적 포텐셜 차이는 Δμ^ex = 35.56 kJ mol⁻¹로, Rick과 Berne의 열역학적 통합 결과인 35.52 kJ mol⁻¹와 뛰어난 일치를 보였다.
• 겹치는 히스토그램, 수락 비율, 기하 평균 보간 방법이 모두 동일한 화학적 포텐셜 차이를 도출하여 방법론적 일관성을 확인하였다.
• 전하 결합 매개변수 λ에 대한 2차 함수 형태에서의 상당한 이탈이 관찰되어 단순 선형 연속체 전기적 모델의 무효성을 입증하였다.
• 자기상호작용 항목을 포함함으로써 체적 의존성이 감소하고 Ewald 및 GRF 시뮬레이션 간의 일치성이 향상되었으며, 특히 작은 시스템(N = 64)에서 두드러졌다.
• 자기상호작용 보정이 없을 경우 Ewald와 GRF 결과 간의 격차는 약 1.7 kJ mol⁻¹에 이르며, 이러한 보정의 필요성을 강조한다.
• 두 잘 측정된 상태 사이의 보간이 고차 누산수 전개보다 더 정확하며, 후자는 분포 꼬리의 열악한 표본 추출으로 인해 문제가 있음을 입증하였다.