[논문 리뷰] Recent achievements in ab initio modelling of liquid water

Rustam Z. Khaliullin, Thomas D. Kühne|arXiv (Cornell University)|2013. 03. 08.

Spectroscopy and Quantum Chemical Studies참고 문헌 198인용 수 43

한 줄 요약

이 논문은 제2세대 Car-Parrinello 분자역학와 절대적으로 국소화된 분자 궤도(ALMO) 에너지 분해 분석을 조합한 방법을 제시하여, 이론적 분자역학 시뮬레이션을 통해 액체 물을 이전까지는 이루어지지 않았던 시간 및 길이 스케일에서 수행할 수 있게 한다. 이 방법은 열역학적 변동에 의해 유도되는 수소결합 네트워크 내 일시적인 전자 비대칭성이 X선 흡수 스펙트럼의 프리엣지 피크를 설명함으로써, 물의 구조에 대한 상반된 모델을 조율한다.

ABSTRACT

The application of newly developed first-principle modeling techniques to liquid water deepens our understanding of the microscopic origins of its unusual macroscopic properties and behaviour. Here, we review two novel ab initio computational methods: second-generation Car-Parrinello molecular dynamics and decomposition analysis based on absolutely localized molecular orbitals. We show that these two methods in combination not only enable ab initio molecular dynamics simulations on previously inaccessible time and length scales, but also provide unprecedented insights into the nature of hydrogen bonding between water molecules. We discuss recent applications of these methods to water clusters and bulk water.

연구 동기 및 목표

• 장시간 체계와 큰 체계에서 액체 물을 시뮬레이션하는 데 있어 전통적인 이론적 분자역학(AIMD)의 계산적 제약를 극복하기 위해.
• 높은 정확도로 수소결합의 전자 기원을 수소결합 클러스터와 액체 상태 물에서 조사하기 위해.
• 특히 프리엣지 피크를 포함한 액체 물의 X선 흡수 스펙트럼 해석에 대한 논란을 해결하기 위해.
• 기존의 대칭적인 정사면체 모델과 최근 제안된 비대칭적인 '사슬과 고리' 모델을 조율하기 위해.

제안 방법

• 제2세대 Car-Parrinello 분자역학(2nd-gen CPMD)을 사용하여, 큰 통합 시간 단계를 활용하면서도 장시간 이론적 시뮬레이션(최대 1 ns)을 가능하게 하여 CPMD의 효율성과 BOMD의 정확성을 결합함.
• 절대적으로 국소화된 분자 궤도(AlMOS)를 사용하여 에너지 분해 분석(EDA)을 수행하여 수소결합 상호작용에 기여하는 전자 기여도를 정량화함.
• ALMO-EDA 방법은 물 분자 간의 상호작용 에너지를 전기적 상호작용, 분극, 전하 이동 등의 물리적으로 의미 있는 구성요소로 분해함.
• 기여자와 수용자 수소결합의 상대 강도에 기반하여 대칭성 파라미터 Υ를 정의하여, 순간적인 비대칭성에 따라 물 분자를 네 그룹으로 동일한 크기로 분류함.
• 각 그룹의 비대칭성에 따라 X선 흡수(XA) 스펙트럼을 계산하여 실험 데이터와 직접 비교할 수 있도록 함.
• 이 방법은 피코초 시간 스케일에서 수소결합 네트워크의 일시적인 구조적 및 전자적 비대칭성을 분석할 수 있음.

실험 결과

연구 질문

• RQ1 액체 물의 X선 흡수 스펙트럼에서 두드러진 프리엣지 피크는 무엇에 기인하는가? 그리고 이것이 장기간 유지되는 구조적 모티프가 아니라 전자 비대칭성에 의해 설명될 수 있는가?
• RQ2 열역학적 변동은 평균적으로 정사면체 네트워크를 유지하고 있음에도 불구하고, 액체 물의 수소결합 네트워크 내에서 일시적인 비대칭성을 얼마나 유도하는가?
• RQ3 기여자와 수용자 수소결합 환경의 비대칭성에 따라 물 분자의 전자적 성질은 어떻게 변화하는가?
• RQ4 향상된 샘플링과 정확도를 갖춘 이론적 시뮬레이션을 통해 오랫동안 끊어진 대칭 모델과 비대칭 모델 간의 논란을 해결할 수 있는가?

주요 결과

• 실험적 액체 물의 X선 흡수 스펙트럼에서 관찰되는 프리엣지 피크는 '고리와 사슬'을 이루는 다수의 분자 때문이 아니라, 매우 비대칭적인 수소결합 환경에 있는 소수의 분자(~25%)에 기인함.
• 평균적으로 정사면체 네트워크이지만, 열역학적 변동으로 인해 어떤 순간에는 한 기여자 또는 수용자 수소결합이 다른 것보다 현저히 강해지며, 이로 인해 일시적인 전자 비대칭성이 발생함.
• 수소결합 강도의 순간적인 비대칭성은 분자 간 진동과 진동각운동으로 인해 수백 피코초 내에 회복됨.
• 비대칭성 그룹별로 분할된 계산된 X선 흡수 스펙트럼은, 프리엣지 피크가 가장 비대칭적인 25%의 분자에서 주로 기여하며, 여기서 가장 강한 상호작용은 두 번째로 강한 것보다 6배 이상 강함.
• 2nd-gen CPMD와 ALMO-EDA의 조합을 통해 나노초 스케일의 궤적과 정확한 전자 구조를 갖춘 액체 물의 시뮬레이션을 가능하게 하여, 유한한 크기와 샘플링 오차를 감소시킴.
• 결과는 물의 이국적인 성질이 지속적인 구조적 모티프 때문이 아니라, 변동하는 수소결합 네트워크 내에서의 동적이고 일시적인 전자 비대칭성에 기인한다는 것을 지지함.