NWChem COSMO模块的最新改进

《Journal of Chemical Theory and Computation》：Recent Improvements to the NWChem COSMO Module

高分辨率图像下载地址：/cms/10.1021/acs.jctc.5c01368/asset/images/large/ct5c01368_0016.jpeg

　　本研究介绍了一项在NWChem软件中对COSMO溶剂化模型的最新改进，旨在提升其在计算溶剂化能方面的准确性和稳定性。COSMO（Conductor-like Screening Model）是一种基于连续介质模型的溶剂效应计算方法，其核心思想是将溶剂视为一种均匀的、理想导体的连续介质，并通过计算溶质在溶剂边界上的表面电荷分布来模拟溶剂对溶质的极化作用。在COSMO模型中，表面电荷分布由一系列离散的点电荷来近似，这些点电荷被放置在分子表面的各个小区域中心，从而能够更精确地描述溶质与溶剂之间的相互作用。

传统的COSMO模型在构建分子腔体（cavity）时通常使用范德华（vdW）表面，这种表面由原子半径定义，每个原子周围被看作一个球体，这些球体相互重叠以形成分子表面。然而，这种方法在某些情况下会导致分子腔体的表面过于粗糙，甚至可能留下一些不合理的空隙，从而影响计算结果的准确性。为了更好地模拟溶剂分子在溶质表面的分布，本研究引入了一种基于溶剂排除表面（SES）的新型腔体构建方法，这种方法能够更精确地反映分子的实际形状，同时避免了传统方法中可能产生的数值问题。

SES表面的构建依赖于GEPOL93算法，这是一种经典的表面生成算法，能够通过添加额外的球体来填充那些不被溶剂分子访问的区域，从而形成一个更为精确的表面模型。然而，这种方法并非完全精确，因为它会在某些区域引入不真实的曲率平滑现象。为了提高计算效率，研究团队对GEPOL93算法进行了优化，使得用户可以根据需要选择不同的细分程度（coarse、fine或xfine），从而在精度与计算效率之间取得平衡。此外，为了应对某些情况下表面区域过于接近所导致的近似奇异性（near-singularity），研究团队引入了一种简单的合并策略，即将距离过近的表面区域合并为一个，并将它们的面积相加，以减少计算过程中可能出现的数值不稳定现象。

除了腔体构建的改进，本研究还对COSMO模型中常见的“外露电荷”（outlying charges）问题进行了深入探讨。在量子化学计算中，由于分子表面的离散化，部分电子密度会分布在腔体之外，这种现象被称为外露电荷。传统的COSMO方法通常通过缩放因子（scaling factors）或拉格朗日乘数（Lagrange multipliers）来修正这一问题，但这些方法往往无法完全消除外露电荷带来的误差，特别是在处理带电分子时。因此，研究团队提出了一种新的修正策略，即在计算完表面电荷后，同时对电势进行修正，以确保电势在修正后的电荷分布下满足导体边界条件。这一改进使得修正后的电势与原始电势之间的偏差大幅降低，从而显著提高了溶剂化能计算的准确性。

为了验证新方法的有效性，研究团队对大约100种化学结构各异的分子进行了计算，包括中性分子、小离子和常见离子液体成分。这些分子的溶剂化能被计算后，与GAMESS软件中采用的双腔体（double-cavity）方法进行对比。尽管双腔体方法被认为在理论上更为精确，但本研究发现，通过适当的电势修正，NWChem中基于缩放因子或拉格朗日乘数的简单修正方法也可以达到与双腔体方法相近的精度。在未进行任何修正的情况下，平均绝对偏差约为0.6 kcal/mol，而在进行修正后，偏差降至0.15 kcal/mol左右，这表明新的修正策略在实际应用中具有很高的可靠性。

此外，研究团队还评估了不同修正方法在COSMO-RS和COSMO-SAC模型中的表现。COSMO-RS是一种基于表面电荷密度的溶剂化模型，能够预测混合物的热力学性质，而COSMO-SAC则进一步引入了方程状态模型，以考虑压力对溶剂化过程的影响。为了确保COSMO模型能够与这些热力学模型兼容，研究团队对NWChem中的COSMO输出格式进行了调整，使其能够直接用于COSMO-RS/SAC的后续计算。这一改进不仅提升了计算效率，还增强了模型在预测混合物相平衡方面的实用性。

在实际应用中，COSMO模型被广泛用于预测分子在不同溶剂中的溶解行为以及相平衡特性。本研究中，基于NWChem的SES腔体和修正后的电势，研究团队开发了一种改进的COSMO-SAC模型（称为CS25），并利用其预测了多种混合物的气液平衡（VLE）和液液平衡（LLE）。结果表明，该模型在预测某些典型混合物的相平衡时表现优异，如醋酸和环己烷的混合物，以及氯仿和乙醚的混合物。然而，在某些情况下，如含有羧酸或胺类分子的混合物，模型仍然存在一定的局限性，特别是在描述分子间缔合现象时。这表明，尽管COSMO模型在许多情况下表现良好，但在某些复杂体系中，仍然需要更精确的修正方法或引入额外的物理模型来提高预测的准确性。

总体而言，本研究的成果为COSMO模型在NWChem中的应用提供了重要的改进。通过引入SES腔体构建方法和更精确的电势修正策略，研究团队显著提高了溶剂化能计算的精度和稳定性。这些改进不仅适用于中性分子，也适用于带电分子，如阴离子。此外，新的COSMO-SAC模型在预测混合物相平衡方面也表现出良好的性能，特别是在处理氢键供体和受体的混合物时，其预测结果与实验数据高度一致。然而，研究团队也指出，目前的修正方法在某些情况下仍然无法完全消除外露电荷带来的影响，尤其是在处理非球形分子时，不同修正方法在空间分布上存在一定的差异，这可能会对高阶偶极矩和多极矩产生一定的扰动。因此，未来的研究方向可能包括开发更加精确的修正方案，以消除外露电荷的影响，同时保持高阶偶极矩和多极矩的完整性。

本研究的成果不仅提升了NWChem在计算溶剂效应方面的能力，也为其他量子化学软件提供了有价值的参考。通过将COSMO模型与NWChem的计算功能紧密结合，研究团队为化学家提供了一种更加高效和准确的工具，以模拟分子在不同溶剂环境下的行为。此外，这些改进也为进一步研究分子间相互作用、混合物的热力学性质以及溶剂化过程的微观机制奠定了基础。随着计算化学在各个领域的广泛应用，COSMO模型的改进将有助于更精确地理解和预测分子在不同环境下的行为，从而推动相关研究的发展。