计算化学公社

标题: 为什么同样的方法在gaussian和orca中算出来的CHELPG电荷不同？ [打印本页]

作者
Author: C41014007 时间: 2016-7-20 19:42
标题: 为什么同样的方法在gaussian和orca中算出来的CHELPG电荷不同？
在gaussian中使用B3LYP/6-31+g*对分子优化，然后用chelpg法计算电荷。同样在orca中用B3LYP/G 6-31+g*对分子优化，再用chelpg计算原子电荷。两个软件优化的几何构型差不多，但是得出的电荷相差比较大。请教各位知道这是什么原因吗？谢谢！

作者
Author: beefly 时间: 2016-7-20 20:03
gaussian版本的b3lyp在orca里是b3lyp/g

作者
Author: C41014007 时间: 2016-7-20 20:24

beefly 发表于 2016-7-20 20:03
gaussian版本的b3lyp在orca里是b3lyp/g

嗯，在ORCA中我对应用的就是B3LYP/G,是不是我在ORCA中计算电荷的时候不能像在gaussian中那样简单输入CHELPG，而是还要定义其它东西呢？

作者
Author: 冰释之川 时间: 2016-7-20 22:15
另外注意高斯和ORCA的默认精度问题，SCF的收敛限，计算的格点密度都会影响到结果的比较

作者
Author: sobereva 时间: 2016-7-20 23:59
可能拟合CHELPG用的积分格点不一样

你用完全相同的结构，在Gaussian里用B3LYP/6-31+G*得到波函数文件(.wfn或.wfx或.fch），再在ORCA里用B3LYP/G 6-31+G*得到.molden文件，放到Multiwfn计算CHELPG电荷，比较差异，以此弄清楚原因。如果还有差异，也许是泛函积分格点精度设定差异，把ORCA里的格点精度弄大点；如果没什么差异，说明是ORCA计算CHELPG电荷的模块不合理。

作者
Author: C41014007 时间: 2016-7-21 09:28
谢谢各位，我用Multiwfn计算试试，比较下结果

作者
Author: beefly 时间: 2016-7-21 16:40

冰释之川发表于 2016-7-20 22:15
另外注意高斯和ORCA的默认精度问题，SCF的收敛限，计算的格点密度都会影响到结果的比较

有可能是收敛精度造成的。orca单点计算的收敛精度比较低，大致相当于gaussian03单点默认的Sleazy。对能量影响不大，但是性质及波函分析在有些情况下会有定性差别

作者
Author: C41014007 时间: 2016-7-21 21:36
用相同的结构、相同的方法分别在gaussian和orca中得到.fchk和.molden文件，再放到Multiwfn中用chelpg去算电荷，感觉结果差异不大。图1是用orca中的文件的结果，图2是用g09中文件的结果。

作者
Author: sobereva 时间: 2016-7-22 11:29

C41014007 发表于 2016-7-21 21:36
用相同的结构、相同的方法分别在gaussian和orca中得到.fchk和.molden文件，再放到Multiwfn中用chelpg去算电 ...

还是有一点差异，诸如5C。
可以把ORCA的收敛限再设严一些、积分格点精度再设高一些，泛函换成没有含糊性的PBE（但不要开RI），基组换成def2SVP，如果在Multiwfn下算的电荷没有区别了，就是数值设定的原因。
如果ORCA之前输出的CHELPG电荷和当前Multiwfn及Gaussian输出的差很多，那就是ORCA自身计算CHELPG的模块不可靠了。

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