[論文レビュー] Universal iso-density polarizable continuum model for molecular solvents

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

Molecular spectroscopy and chirality被引用数 3

ひとこと要約

本論文は、分子性溶媒用の普遍的等密度PCM（極性誘電率モデル）を導入し、溶媒固有のパラメータフィッティングを必要とせず、体積平均熱力学的データと単一分子ab initio計算のみを用いて、高精度な溶媒化自由エネルギー予測を可能にする。この手法は、非線形誘電応答と溶媒の四極モーメント効果の補正を組み込むことで、プロトニックおよびアプロトニック溶媒、特にエチレングリコールを含む多様な溶媒において高い精度を達成する。

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の適用性を拡大すること。
• 容易に入手可能な体積平均熱力学的性質と単一分子量子化学的計算にのみ依存する、パラメータフリーまたは最小限のパrameter化を実現するアプローチの開発。
• 非電気的および非局所的効果（例：溶媒の四極モーメント）の取り扱いを改善し、標準モデルで不十分に捉えられるこれらの効果を補正すること。

提案手法

• 溶質の電子密度を用いてキャビティを定義することで、平面波DFTコードとの互換性を確保する。
• 極性溶媒における固定双極子の回転飽和を考慮する非線形誘電応答関数を組み込む。
• 溶媒固有のパrameterは、体積平均熱力学的データ（例：誘電率、モル体積）と、双極子および四極モーメントの単一分子ab initio計算から導出する。
• ネット双極モーメントがゼロだが四極モーメントが顕著な溶媒（例：エチレングリコール）に対しては、全分子ではなく断片ベースの双極子密度を用いることで補正を施す。
• 水、クロロホルム、エチレングリコールなど多様な溶媒における実験的溶媒化エネルギーと比較してモデルを検証する。
• 計算実装にはJDFTxを用い、revTPSS準GGA関数とノルム保存ポテンシャルを採用する。

実験結果

リサーチクエスチョン

• RQ1 各溶媒について実験的溶媒化エネルギーデータベースへのフィッティングを必要とせず、普遍的等密度PCMを構築できるか？
• RQ2 体積平均熱力学的データと単一分子ab initio計算のみを用いて、溶媒化自由エネルギーをどの程度の精度で予測できるか？
• 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の誤差を示すが（悪化）、補正済みバージョンでは異なる溶質間で一貫性が向上する。
• 断片ベースの双極子密度の導入により、ネット双極モーメントがゼロだが局所的双極子相互作用が強いエチレングリコールのような溶媒の予測精度が顕著に向上する。
• 溶媒固有のパラメータフィッティングを一切必要とせず、多様な溶媒において高精度な溶媒化エネルギー予測が可能になる。
• 計算コストが低く、平面波DFTと互換性があるため、固体-液体界面や溶液中における物質の応用が可能になる。