计算化学公社

标题: bsse矫正总是出现l502错误 [打印本页]

作者
Author: xiaoshan 时间: 2021-12-30 10:42
标题: bsse矫正总是出现l502错误
优化好无虚頻的二体另存为gjf以后再进行bsse计算，总显示l502，其他二体都没出现这种情况，不知道这个是怎么回事
(, 下载次数 Times of downloads: 23)

作者
Author: snljty 时间: 2021-12-30 11:01
和BSSE没直接关系。
看卢老师博文《解决SCF不收敛问题的方法》http://sobereva.com/61 。
不要说Lxxxx报错，看《Gaussian 常见报错及解决方法；新手求助报错时的注意事项》http://bbs.keinsci.com/thread-4829-1-1.html

作者
Author: xiaoshan 时间: 2021-12-30 11:27
好的，谢谢您，我学习一下

作者
Author: xiaoshan 时间: 2021-12-30 11:30

snljty 发表于 2021-12-30 11:01
和BSSE没直接关系。
看卢老师博文《解决SCF不收敛问题的方法》http://sobereva.com/61 。
不要说Lxxxx报 ...

我想问一下，优化好的构型，再加bsse矫正，怎么又出现不收敛的情况呢，之前没遇到过这种情况

作者
Author: snljty 时间: 2021-12-30 12:57

xiaoshan 发表于 2021-12-30 11:30
我想问一下，优化好的构型，再加bsse矫正，怎么又出现不收敛的情况呢，之前没遇到过这种情况

1.不知道你说的是复合物基组算单体时候还是复合物时候SCF不收敛。
2.几何优化每一步读取上一步结构收敛的波函数做初猜，BSSE里面的算单点任务每步重新构造初猜。
3.一般几何优化相对不需要很大的基组，但是计算BSSE往往需要更大的基组精度才比较高。由于大概就是复合物单点能三倍左右耗时，一般用计算BSSE的基组会比几何优化更大，而更大基组尤其有很大弥散的时候自洽场就是会更难收敛。

作者
Author: 喵星大佬 时间: 2021-12-30 13:02

snljty 发表于 2021-12-30 12:57
1.不知道你说的是复合物基组算单体时候还是复合物时候SCF不收敛。
2.几何优化每一步读取上一步结构收敛 ...

先算小基组单体用鬼原子不带基组，加基组读波函数重算，最后手工BSSE校正可以完美解决

作者
Author: snljty 时间: 2021-12-30 13:21

喵星大佬发表于 2021-12-30 13:02
先算小基组单体用鬼原子不带基组，加基组读波函数重算，最后手工BSSE校正可以完美解决

Gaussian的话，第一步要加鬼原子么？

作者
Author: 冰释之川 时间: 2021-12-30 13:44

snljty 发表于 2021-12-30 13:21
Gaussian的话，第一步要加鬼原子么？

手工算BSSE要设置鬼原子的

作者
Author: snljty 时间: 2021-12-30 14:40

冰释之川发表于 2021-12-30 13:44
手工算BSSE要设置鬼原子的

我理解的喵哥的意思是要算复合物基组下的单体波函数的时候，为了降低SCF不收敛的概率，用单体基组下的单体波函数做初猜，再读这个初猜做复合物基组下的单体波函数计算，是这个意思吧？那第一步没必要加不带基组的鬼原子吧

作者
Author: 喵星大佬 时间: 2021-12-30 15:15

snljty 发表于 2021-12-30 14:40
我理解的喵哥的意思是要算复合物基组下的单体波函数的时候，为了降低SCF不收敛的概率，用单体基组下的单 ...

理论上没必要，不过没试过读原子数不一样的波函数，不知道有没有问题，这样肯定没问题就是了。

而且BSSE不管自动的还是手动的，用DFT-D3的话D3校正部分都要手动处理吧

作者
Author: snljty 时间: 2021-12-30 15:20

喵星大佬发表于 2021-12-30 15:15
理论上没必要，不过没试过读原子数不一样的波函数，不知道有没有问题，这样肯定没问题就是了。

而且BS ...

前者我之前试过好像没有问题，不过确实看起来按照您说的这样更稳妥。
后者我没明白为什么，是因为Gaussian没有对Bq原子指定D3系数所以会报错？不过带有基函数的Bq原子比如C-Bq好像不会有这个报错，之前看到卢老师的博文里用法也没提到要特殊处理D3...还请赐教

作者
Author: 喵星大佬 时间: 2021-12-30 15:31
本帖最后由喵星大佬于 2021-12-30 15:33 编辑

snljty 发表于 2021-12-30 15:20
前者我之前试过好像没有问题，不过确实看起来按照您说的这样更稳妥。
后者我没明白为什么，是因为Gaussi ...

我就是说不带基函数的Bq没有D3参数，不过讲道理用DFT-D这个级别计算即使是弱相互作用计算，基组大点，比如大到aug-cc-pVTZ基本也不存在什么BSSE问题了，就算是在瓜丝里稍微好点的计算条件算个单点也不是什么大问题，这种时候再过多考虑BSSE感觉多此一举了，或者说这个时候额外考虑BSSE的影响未必能有改进。

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