为何用orac算单点能在有无counterpoise下相差大

calvados80

近来算一个金属配体的单点能。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会高估相互作用能，这里得出的结论是相反的，为什么？烦请大佬指出我的错误






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

calvados80

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

冰释之川

标记冒号的原子都是鬼原子，只提供基函数

calvados80

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

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

sobereva

注意是ORCA不是orac

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

calvados80

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

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

谢谢老师提醒！