[论文解读] 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），仅通过宏观热力学数据和单分子从头算计算即可实现溶剂化自由能的精确预测，无需针对每种溶剂进行实验溶剂化能数据库的拟合。该方法在包括乙二醇在内的多种质子性与非质子性溶剂中均表现出高精度，通过引入非线性介电响应并校正溶剂四极矩效应实现。

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代码。
• 引入非线性介电响应函数，以考虑极性溶剂中永久偶极子的取向饱和效应。
• 溶剂特异性参数由宏观热力学数据（如介电常数、摩尔体积）和单分子从头算计算得到的偶极矩与四极矩确定。
• 对于净偶极矩为零但具有显著四极矩的溶剂（如乙二醇），通过使用片段基偶极密度而非完整分子进行校正。
• 模型通过与多种溶剂（包括水、氯仿和乙二醇）的实验溶剂化能对比进行验证。
• 计算实现采用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（更差），但校正后版本在不同溶质间表现出更一致的预测性能。
• 引入片段基偶极密度显著提升了对乙二醇等溶剂的预测精度，尽管其净偶极矩为零，但局域偶极相互作用强烈。
• 该方法可在无需溶剂特异性参数化的情况下，实现对多种溶剂的精确溶剂化能预测。
• 该方法计算效率高，兼容平面波DFT，适用于固液界面及溶液中材料的模拟。