计算化学公社

标题: 为何用orac算单点能在有无counterpoise下相差大 [打印本页]

作者
Author: calvados80 时间: 2020-3-31 19:32
标题: 为何用orac算单点能在有无counterpoise下相差大
本帖最后由 calvados80 于 2020-3-31 20:30 编辑

近来算一个金属配体的单点能。A为金属，B为配体。一种是用counterpoise做，计算级别为 RI-PWPB95-D3(BJ)/def2-TZVPP，得到三个能量分别为FINAL SINGLE POINT ENERGY    -3175.250830114327  E1（A+B）
FINAL SINGLE POINT ENERGY    -1778.444051727978 E2（A）
FINAL SINGLE POINT ENERGY    -1395.449367733336  E3    (B)       这里得到的相互作用能=E1-E2-E3=-1.3574
这个是我用mutiwfn导出的orac执行文件，三个单点任务电荷都是+2，自旋多重度为1。但其实后期我要上SMD考虑溶液情况下的相互作用能，这样就不好用counterpoise，于是我在计算级别完全相同的条件下，分别计算了A+B（+2,1）；A（+2,1）,B(0,1)的单点能，其中有机配体B是优化好的A+B直接砍掉A原子，然后不结构优化直接拿RI-PWPB95-D3(BJ)/def2-TZVPP算的。结果如下
E1=-3175.25083
E2=-1778.44102
E3=-1396.106295
E1 E2差距较小，但发现E3显著下降，这样相互作用能只有-0.70（暂不考虑 frequency thermal correction ），差距很大了
我想问E3为什么数值会下降？BSSE会高估相互作用能，这里得出的结论是相反的，为什么？烦请大佬指出我的错误

最后我想问如何判断过渡金属是高自旋还是低自旋，求一篇综述性文献学习一下，谢谢！

作者
Author: calvados80 时间: 2020-3-31 20:24
我观察到http://sobereva.com/542中sob老师提供文件的不同片段的原子坐标上打了冒号，这是屏蔽掉这个原子的意思吗？那样片段1和片段2就要分开设定电子和多重度了，不知道我的猜想对不对

作者
Author: 冰释之川 时间: 2020-3-31 21:02
标记冒号的原子都是鬼原子，只提供基函数

(, 下载次数 Times of downloads: 12)

作者
Author: calvados80 时间: 2020-3-31 21:27

冰释之川发表于 2020-3-31 21:02
标记冒号的原子都是鬼原子，只提供基函数

谢谢老师！鬼原子不带电子和核电荷，那显然我counterpoise片段设定的电荷数错了

作者
Author: sobereva 时间: 2020-4-1 02:49
注意是ORCA不是orac

各种自旋多重度都算算，对比能量便知该用什么自旋多重度。这种常见问题应当先搜公社论坛首页google搜索框，我回答过无数次

作者
Author: calvados80 时间: 2020-4-1 13:22

sobereva 发表于 2020-4-1 02:49
注意是ORCA不是orac

各种自旋多重度都算算，对比能量便知该用什么自旋多重度。这种常见问题应当先搜公社 ...

谢谢老师提醒！

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