bsse矫正总是出现l502错误

xiaoshan

优化好无虚頻的二体另存为gjf以后再进行bsse计算，总显示l502，其他二体都没出现这种情况，不知道这个是怎么回事

snljty

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

xiaoshan

好的，谢谢您，我学习一下

xiaoshan

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

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

snljty

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

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

喵星大佬

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

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

snljty

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

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

冰释之川

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

手工算BSSE要设置鬼原子的

snljty

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

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

喵星大佬

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

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

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

snljty

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

而且BS ...

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

喵星大佬

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

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