计算化学公社

标题: BSSE计算单点能l502报错 [打印本页]

作者
Author: jnnn16 时间: 2019-3-19 15:21
标题: BSSE计算单点能l502报错
我的过渡态已经优化好了，在使用BSSE校正计算单点能时，总是出现l502报错，不知道是输入文件写错了还是怎么回事。我的输入及输出文件见附件 (, 下载次数 Times of downloads: 8)

作者
Author: sobereva 时间: 2019-3-19 16:57
当前片段划分完全莫名其妙，算弱相互作用也不是，算化学键键能也不是。先把片段设好了再说

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

如果要算化学键键能，根本不需要也不应该用CP考虑BSSE，仔细看
计算化学键键能时若基组合理则考虑BSSE是完全多余的
http://sobereva.com/381（http://bbs.keinsci.com/thread-6047-1-1.html）

作者
Author: jnnn16 时间: 2019-3-19 18:15

sobereva 发表于 2019-3-19 16:57
当前片段划分完全莫名其妙，算弱相互作用也不是，算化学键键能也不是。先把片段设好了再说

Sob老师，非常感谢您的回复。请问我这个该如何划分片段呢？我的体系由于弱相互作用影响较高，请问有什么更好的方法去处理吗？谢谢

作者
Author: sobereva 时间: 2019-3-19 18:28

jnnn16 发表于 2019-3-19 18:15
Sob老师，非常感谢您的回复。请问我这个该如何划分片段呢？我的体系由于弱相互作用影响较高，请问有什么 ...

怎么划分完全取决于你到底想研究什么
必须让能量差体现出来你要研究的量

作者
Author: l1234 时间: 2025-6-3 10:43

sobereva 发表于 2019-3-19 18:28
怎么划分完全取决于你到底想研究什么
必须让能量差体现出来你要研究的量

老师您好，我也遇到了类似的情况，我用oniom对二聚体进行优化和频率分析正常收敛高层分子没有虚频，但用优化好的结构做BSSE却出现了l502报错，这是什么原因呢？需要在原来结构基础上重新进行优化吗？（两个片段分别定义为二聚体中的两个单分子）下面是输入文件
%chk=2c-1-T-bsse.chk
%nprocshared=8
%mem=12000MB
#p B3LYP/6-31G** EmpiricalDispersion=GD3 counterpoise=2 test
2c-1-T-bsse
0 3
B 29.7630845246 20.4873380575 29.3669098956 1
……

作者
Author: wzkchem5 时间: 2025-6-3 11:56

l1234 发表于 2025-6-3 10:43
老师您好，我也遇到了类似的情况，我用oniom对二聚体进行优化和频率分析正常收敛高层分子没有虚频，但用 ...

与其用这么小的基组做计算然后做counterpoise校正，不如用def2-TZVP做计算然后加gCP校正

作者
Author: l1234 时间: 2025-6-3 15:13

wzkchem5 发表于 2025-6-3 11:56
与其用这么小的基组做计算然后做counterpoise校正，不如用def2-TZVP做计算然后加gCP校正

非常感谢老师的回复，老师那为什么结构优化收敛且没有虚频的情况下但算这个相互作用能却不收敛了呢？除了改变基组外还有别的解决方法吗？

作者
Author: wzkchem5 时间: 2025-6-3 18:29

l1234 发表于 2025-6-3 15:13
非常感谢老师的回复，老师那为什么结构优化收敛且没有虚频的情况下但算这个相互作用能却不收敛了呢？除了 ...

SCF是否收敛和结构优化是否收敛没有任何关系，“结构优化收敛但SCF却不收敛”这个“却”字完全不成立。
Counterpoise计算因为引入了大量并非局域在原子上的基组，对积分格点等要求高。所以核心在于避免用counterpoise。一定要注意counterpoise不是校正BSSE的唯一方法，而且校正BSSE也不是在所有情况下都是必须的。当基组足够大（如quadrupole zeta）时，不用校正BSSE；当基组比较大但没那么大（triple zeta）且精度要求不是特别高时，可以用近似的方法校正BSSE，例如gCP。如果基组只有double zeta，作用能即使校正了BSSE，精度也不够，结果不能发表

作者
Author: l1234 时间: 2025-6-4 08:54

wzkchem5 发表于 2025-6-3 18:29
SCF是否收敛和结构优化是否收敛没有任何关系，“结构优化收敛但SCF却不收敛”这个“却”字完全不成立。
...

万分感谢老师的回复，学生明白了

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