为什么同样的方法在gaussian和orca中算出来的CHELPG电荷不同？

C41014007 · 发表于 Post on 2016-7-20 19:42:43

在gaussian中使用B3LYP/6-31+g*对分子优化，然后用chelpg法计算电荷。同样在orca中用B3LYP/G 6-31+g*对分子优化，再用chelpg计算原子电荷。两个软件优化的几何构型差不多，但是得出的电荷相差比较大。请教各位知道这是什么原因吗？谢谢！

sobereva · 发表于 Post on 2016-7-22 11:29:01

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

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

C41014007 · 发表于 Post on 2016-7-21 21:36:09

用相同的结构、相同的方法分别在gaussian和orca中得到.fchk和.molden文件，再放到Multiwfn中用chelpg去算电荷，感觉结果差异不大。图1是用orca中的文件的结果，图2是用g09中文件的结果。

beefly · 发表于 Post on 2016-7-21 16:40:39

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

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

C41014007 · 发表于 Post on 2016-7-21 09:28:42

谢谢各位，我用Multiwfn计算试试，比较下结果

sobereva · 发表于 Post on 2016-7-20 23:59:23

可能拟合CHELPG用的积分格点不一样

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

冰释之川 · 发表于 Post on 2016-7-20 22:15:47

另外注意高斯和ORCA的默认精度问题，SCF的收敛限，计算的格点密度都会影响到结果的比较

C41014007 · 发表于 Post on 2016-7-20 20:24:51

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

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

beefly · 发表于 Post on 2016-7-20 20:03:52

gaussian版本的b3lyp在orca里是b3lyp/g

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

[新手求助] 为什么同样的方法在gaussian和orca中算出来的CHELPG电荷不同？

评分 Rate

浏览过的版块