计算化学公社

标题: 计算NPA，所得电荷Total与初始设置不一致 [打印本页]

作者
Author: guo7212c 时间: 2015-6-25 10:00
标题: 计算NPA，所得电荷Total与初始设置不一致
本帖最后由 guo7212c 于 2015-6-25 10:00 编辑

读取opt所得chk文件，设置计算：# b3lyp/6-311+G(d,p) guess=read scrf=(solvent=water) geom=check pop=npa。所得结果出现——row 1957: Total -0.98233，与初始设置（-1）不一致，Mulliken 电荷计算是 -1。原始文件参见附件。 (, 下载次数 Times of downloads: 27)

作者
Author: sobereva 时间: 2015-6-25 10:33
把弥散函数去掉再试，肯定就没问题了。
NBO分析对基组很不敏感，弥散函数更是不需要。而有弥散函数时则容易引入各种数值问题。

作者
Author: guo7212c 时间: 2015-6-25 22:09

sobereva 发表于 2015-6-25 10:33
把弥散函数去掉再试，肯定就没问题了。
NBO分析对基组很不敏感，弥散函数更是不需要。而有弥散函数时则容 ...

由于计算资源所限，现在才计算完成。如sob老师所说，去掉弥散函数后，NPA电荷显示正常。谢谢sob老师。
这里还有有个问题，看您的论文《原子电荷计算方法的对比》中提到，NPA并不适合分析静电相互作用。我的目的是分析表面活性剂A（阴离子）的电荷分布情况以及添加阳离子（如Na+）后对其电荷分布的影响，另A与水分子之间的作用情况。按照您所说，此处NPA并不适用，那么选用您后面提到的MK电荷或未提到的Chelpg如何？想确认，分析分子间作用时是否要选用同样的基组（6-311+G**）进行计算？

作者
Author: sobereva 时间: 2015-6-26 05:04

guo7212c 发表于 2015-6-25 22:09
由于计算资源所限，现在才计算完成。如sob老师所说，去掉弥散函数后，NPA电荷显示正常。谢谢sob老师。
...

此时用NPA问题不大。关键是不能用NPA电荷根据库仑公式去直接计算静电相互作用能，否则误差很大。
复合物和每个片段计算时用的基组应当一致。

作者
Author: guo7212c 时间: 2015-6-26 09:43

sobereva 发表于 2015-6-26 05:04
此时用NPA问题不大。关键是不能用NPA电荷根据库仑公式去直接计算静电相互作用能，否则误差很大。
复合物 ...

明白了，NPA计算所得的电荷也就不能用在分子模拟上，谢谢Sob老师，老师注意身体！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)