计算化学公社

标题: 用Gaussian做BSSE修正时遇见 Overlap in AMove.报错 [打印本页]

作者
Author: poly-Z 时间: 2025-2-21 19:28
标题: 用Gaussian做BSSE修正时遇见 Overlap in AMove.报错
用Gaussian16算两个分子的结合能
关键词是#p cc-pvtz counterpoise=2 m062x int=grid=300590 SCF=novaracc guess=INDO
算了一个多小时后出现报错
不知道是什么原因造成的，附上out文件

作者
Author: sobereva 时间: 2025-2-22 08:57
cc-pVTZ根本不适合DFT算弱相互作用能，本来对DFT计算就不划算，而且弥散函数都没有（更何况当前显著参与相互作用的有不少显著带负电荷的原子），光靠counterpoise也救不回来

改成6-311+G(2d,p)结合counterpoise，对于DFT算弱相互作用又划算又好，看下文2.3节的表格的对比
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

用合理的级别，有什么问题到时候再说

作者
Author: poly-Z 时间: 2025-2-22 18:20

sobereva 发表于 2025-2-22 08:57
cc-pVTZ根本不适合DFT算弱相互作用能，本来对DFT计算就不划算，而且弥散函数都没有（更何况当前显著参与相 ...

感谢卢老师回复，但是我换成6-311+G（2d，p）之后，运行五十多分钟之后又出现了一样的报错

作者
Author: poly-Z 时间: 2025-2-22 19:53

poly-Z 发表于 2025-2-22 18:20
感谢卢老师回复，但是我换成6-311+G（2d，p）之后，运行五十多分钟之后又出现了一样的报错

补一句，在计算相互作用能之前先用M062X,6-31G(d)做过结构优化

作者
Author: sobereva 时间: 2025-2-22 20:47

poly-Z 发表于 2025-2-22 18:20
感谢卢老师回复，但是我换成6-311+G（2d，p）之后，运行五十多分钟之后又出现了一样的报错

去掉out=wfn guess=INDO

作者
Author: poly-Z 时间: 2025-2-23 02:01

sobereva 发表于 2025-2-22 20:47
去掉out=wfn guess=INDO

老师，还是不行，报错原因和之前一样

作者
Author: sobereva 时间: 2025-2-23 02:03

poly-Z 发表于 2025-2-23 02:01
老师，还是不行，报错原因和之前一样

按下文自己利用鬼原子手动做counterpoise
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46
或者尝试其它版本
当前是bug

作者
Author: poly-Z 时间: 2025-2-23 13:54

sobereva 发表于 2025-2-23 02:03
按下文自己利用鬼原子手动做counterpoise
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

谢谢卢老师的指导，用鬼原子解决了

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