如何解决高斯几何优化后计算的吸附能与psi4分解后总能量相差较大

slope · 发表于 Post on 2023-3-21 21:57:11

使用高斯wb97xb def2-svp计算吸附能为11.23kJ/mol，但是使用psi4的sapt0 def2-svp计算的总能量为18.70kJ/mol。这两个数值好像差得有点多，如何解决这种情况。psi4不使用推荐的sSAPT0/jun-cc-pVDZ（铜），因为后面的计算中有K，jun-cc-pVDZ好像没有定义，这样的计算结果是否可用，是否有价值

1.inp.dat (97.76 KB, 下载次数 Times of downloads: 8)

sobereva · 发表于 Post on 2023-3-21 23:15:49

这种对比没任何意义
DFT算弱相互作用能用def2-SVP太烂，不仅本身质量low，BSSE问题还极其严重。SAPT0结合def2-SVP也是极其糟糕的选择，没任何意义。要比也是两个靠谱的级别之间对比。K可以用def2-QZVP
不要把相互自由能叫做总能量

zjxitcc · 发表于 Post on 2023-3-22 01:13:05

def2SVP过于垃圾。本帖附件文件中无K元素，因此建议使用铜级别默认基组jun-cc-pVDZ。你发的另一帖子里有K元素，K可用def2TZVP之类的基组，其余元素仍可使用jun-cc-pVDZ。

七尺贱 · 发表于 Post on 2023-3-22 10:10:53

SAPT0结合def2-SVP我也见过一些文章，例如Phys Chem Chem Phys, 2019, 21, 16665-16675 和Phys. Chem. Chem. Phys., 2018, 20, 11347-11358，但这些都是巨大体系，原子数一般都200+，体系不大非常不推荐用这个组合

slope · 发表于 Post on 2023-3-22 15:50:43

sobereva 发表于 2023-3-21 23:15
这种对比没任何意义
DFT算弱相互作用能用def2-SVP太烂，不仅本身质量low，BSSE问题还极其严重。SAPT0结合d ...

谢谢sob老师的回复，那我将DFT与psi的计算换精度更高的基组重新进行计算

slope · 发表于 Post on 2023-3-22 15:53:02

zjxitcc 发表于 2023-3-22 01:13
def2SVP过于垃圾。本帖附件文件中无K元素，因此建议使用铜级别默认基组jun-cc-pVDZ。你发的另一帖子里有K元 ...

谢谢谢谢，我换一个重新算

slope · 发表于 Post on 2023-3-22 15:53:31

七尺贱发表于 2023-3-22 10:10
SAPT0结合def2-SVP我也见过一些文章，例如Phys Chem Chem Phys, 2019, 21, 16665-16675 和Phys. Chem. Chem ...

谢谢回复，我再去看一下

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

[新手求助] 如何解决高斯几何优化后计算的吸附能与psi4分解后总能量相差较大

评分 Rate

浏览过的版块