计算化学公社

标题: 溶剂化模型下，作用能Eint计算问题 [打印本页]

作者
Author: 枕上卢生 时间: 2021-7-11 23:30
标题: 溶剂化模型下，作用能Eint计算问题
老师，您好！现在遇到一个i棘手的问题，就是在使用溶剂化模型结构优化后能够正常结束，可是计算作用能时，不能再溶剂化模型下计算，一直显示L1错误。当去掉溶剂化关键词时能够正常运行。可是计算出来的作用能很高，想请教一下老师，这种问题如何解决？如何在溶剂化环境下计算作用能？相应的输入与输出文件见附件。期待老师的解答！

作者
Author: 喵星大佬 时间: 2021-7-12 00:16
不用考虑BSSE，因为虚原子没有范德华半径所以没法定义溶剂空腔。

而且溶剂模型的误差比BSSE大，算单点的时候基组够大就行了，直接就正常算结合能就好了不用管counterpoise

作者
Author: 枕上卢生 时间: 2021-7-12 09:36
老师，您的意思是计算溶剂化环境下化合物的作用能不需要加入溶剂化关键词吗？

作者
Author: 喵星大佬 时间: 2021-7-12 23:15
我的意思是，让你不要写counterpoise=2，因为在这个关键词下定义了带基组的虚原子，但是在做隐式溶剂模型的时候需要原子的范德华半径来生成溶剂化空腔，这些虚原子是没有范德华半径的自然也生成不了。
而且要考虑溶剂结合能的话显然要考虑溶剂效应，尤其是非极性部分(需要溶剂化表面的)，所以直接正常在溶剂模型下做片段相减就好了，不用counterpoise=2，基组选大些，比如def2qzvp之类的或者aug-cc-pVTZ之类的，此时BSSE已经可以忽略了，其他的处理BSSE的方法，比如gcp在此类基组下已经没有定义了。

作者
Author: 枕上卢生 时间: 2021-7-13 09:21
好的老师，明白了，谢谢您的解答。非常感谢！

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