计算NPA，所得电荷Total与初始设置不一致

guo7212c · 发表于 Post on 2015-6-25 10:00:53

本帖最后由 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。原始文件参见附件。

LA-npa.log (116.56 KB, 下载次数 Times of downloads: 27)

sobereva · 发表于 Post on 2015-6-25 10:33:22

把弥散函数去掉再试，肯定就没问题了。
NBO分析对基组很不敏感，弥散函数更是不需要。而有弥散函数时则容易引入各种数值问题。

guo7212c · 发表于 Post on 2015-6-25 22:09:12

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

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

sobereva · 发表于 Post on 2015-6-26 05:04:45

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

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

guo7212c · 发表于 Post on 2015-6-26 09:43:24

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[NBO相关] 计算NPA，所得电荷Total与初始设置不一致

浏览过的版块