[논문 리뷰] Universal iso-density polarizable continuum model for molecular solvents

Deniz Gunceler, T. A. Arias|arXiv (Cornell University)|2014. 03. 25.

Molecular spectroscopy and chirality인용 수 3

한 줄 요약

이 논문은 분자 용매를 위한 유니버설 등밀도 폴라라이즈러블 연속체 모델(PCM)을 소개하며, 용매별로 실험적 용해 자유 에너지 데이터베이스에 대한 피팅이 필요 없이, 단지 거시적 열역학적 데이터와 단일 분자 아비-아이티 계산만을 사용하여 정확한 용해 자유 에너지 예측을 가능하게 한다. 이 방법은 비극성 및 비극성 용매를 포함한 다양한 용매에서 높은 정확도를 달성하며, 비선형 유전율 반응과 용매의 쿼드루플 전하 효과를 보정함으로써 성능을 향상시킨다.

ABSTRACT

Implicit electron-density solvation models based on joint density-functional theory offer a computationally efficient solution to the problem of calculating thermodynamic quantities of solvated systems from firstprinciples quantum mechanics. However, despite much recent interest in such models, to date the applicability of such models in the plane-wave context to non-aqueous solvents has been limited because the determination of the model parameters requires fitting to a large database of experimental solvation energies for each new solvent considered. This work presents an alternate approach which allows development of new iso-density models for a large class of protic and aprotic solvents from only simple, single-molecule ab initio calculations and readily available bulk thermodynamic data.

연구 동기 및 목표

• 기존의 등밀도 PCM이 각 새로운 용매에 대해 실험적 용해 자유 에너지 데이터베이스에 광범위하게 피팅이 필요하다는 한계를 극복하기 위해.
• 등밀도 PCM의 적용 범위를 극성 및 비극성 용매 전반에 넓히며, 실험적 용해 데이터가 제한된 용매들까지 포함하기 위해.
• 사용 가능한 거시적 열역학적 성질과 단일 분자 양자화학 계산만을 기반으로 한 매개변수 없음 또는 최소 매개변수화된 접근법을 개발하기 위해.
• 비전자적 및 국소적이지 않은 효과(예: 용매의 쿼드루플 전하 모멘트)의 처리를 향상시키기 위해, 기존 모델에서 잘 포착되지 않는 요소들을 보완하기 위해.

제안 방법

• 모델은 용질의 전자 밀도를 기반으로 캐비티를 정의하여 평면파 DFT 코드와의 호환성을 확보한다.
• 극성 용매에서 영구 전기 dipole의 회전 포화를 고려한 비선형 유전율 반응 함수를 포함한다.
• 용매 특성 매개변수들은 거시적 열역학적 데이터(예: 유전율, 몰 부피)와 분자 단위의 아비-아이티 계산을 통해 유도된 전기 dipole 및 쿼드루플 모멘트로부터 유도된다.
• 순수 전기 dipole가 0이지만 쿼드루플 모멘트가 뚜렷한 용매(예: 에틸렌 글리콜)의 경우, 전체 분자 대신 프ragment 기반 전기 dipole 밀도를 사용하여 보정을 적용한다.
• 수소결합 및 비극성 용매를 포함한 다양한 용매에서의 실험적 용해 자유 에너지와의 비교를 통해 모델을 검증한다.
• 계산 구현은 JDFTx를 사용하며, revTPSS 메타-GGA 함수와 노름-보존 허위파면을 적용한다.

실험 결과

연구 질문

• RQ1 각 용매에 대해 실험적 용해 자유 에너지 데이터베이스에 피팅이 필요 없이 유니버설 등밀도 PCM을 구축할 수 있는가?
• RQ2 단지 거시적 열역학적 데이터와 단일 분자 아비-아이티 계산만을 사용하여 용해 자유 에너지를 얼마나 정확하게 예측할 수 있는가?
• RQ3 용매의 쿼드루플 모멘트가 용해 에너지 예측에 미치는 영향은 무엇이며, 모델에서 이를 어떻게 보정할 수 있는가?
• RQ4 모델은 극성 및 비극성 용매, 특히 복잡한 분자 대칭성을 가진 용매에서도 용해를 정확하게 기술할 수 있는가?
• RQ5 비선형 유전율 반응을 포함함으로써 선형 모델 대비 예측 성능이 어떻게 향상되는가?

주요 결과

• 에틸렌 글리콜에서 벤젠의 용해에 대해 평균 절대 오차가 1.41 kcal/mol이며, 쿼드루플 보정을 적용하면 0.56 kcal/mol으로 향상된다.
• 플루오로벤젠의 경우, 보정 후 실험값과 0.40 kcal/mol 이내로 예측되며, 보정 없이 약 1.09 kcal/mol의 오차를 보인다.
• 나프탈렌의 경우 보정 없이 오차가 1.80에서 3.95 kcal/mol로 악화되지만, 보정된 버전은 다양한 용질에 대해 보다 일관된 예측 성능을 보인다.
• 프래그먼트 기반 전기 dipole 밀도의 도입은 에틸렌 글리콜과 같이 순수 전기 dipole가 0이지만 局소 전기 dipole 상호작용이 강한 용매의 경우 예측 정확도를 크게 향상시킨다.
• 모델은 용매별 매개변수 최적화 없이도 다양한 용매에서 정확한 용해 에너지 예측이 가능하다.
• 이 방법은 계산적으로 효율적이며 평면파 DFT와 호환되어 고체-액체 인터페이스 및 용액 내 물질의 응용이 가능하다.