高斯不收敛,如何优化
#p b3lyp/lanl2dz opt pop=(full,npa)
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(UB+HF-LYP) = -1714.45115901 A.U. after 129 cycles
Convg = 0.8068D-05 -V/T = 2.0529
S**2 = 0.9172
KE= 1.628346096785D+03 PE=-1.275059909267D+04 EE= 5.005782636124D+03
Annihilation of the first spin contaminant:
S**2 before annihilation 0.9172, after 0.7568
Convergence failure -- run terminated.
Error termination via Lnk1e in d:program filesg03l502.exe at Fri Jun 18 18:59:12 2021.
Job cpu time: 1 days 5 hours 55 minutes 14.0 seconds.
File lengths (MBytes): RWF= 206 Int= 0 D2E= 0 Chk= 18 Scr= 1
请问高手,我这个从报错看应该是没有收敛,但是应该如何优化或者修改?救救孩子吧。。。。@jiewei
返回小木虫查看更多
京公网安备 11010802022153号
看起来是高斯的Scratch文件夹里内容被占用了。你修改好之后保存成另外一个名字的输入文件。并且把高斯的scratch文件夹内容清空 再计算就可以了。
包含重金属,也不一定非要所有原子都用LanL2DZ 可以用混合基组 给重原子用LanL2DZ(这个2021年的今天已经很不推荐了精度太差,建议SDD) 给其他轻原子(前三周期)用普通的基组比如 6-31G(D)之类的。会合理很多。
之前那个任务直接改不太行的。要改成混合基组需要特别的写法。见我这个教程1.1节 https://muchong.com/t-14242735-1
如果基组设置没有问题。电荷和多重度也设置的没有什么问题。结构本身原子间的距离并没有太大异常而仅仅是gaussview不显示某些键那一般不会有什么问题。
见我这个教程3.3节 https://muchong.com/t-14242735-1
3.3 化学键的显示问题。
诸如 为什么优化完成后“缺了什么什么键“,“哪个键变成了单/双/三键”,“那个键不见了”等问题。
这类问题的回答是:没有关系。gaussview怎么显示的,都没什么所谓。键显示的符合预期也好或者不符合预期也罢都不能说明计算结果本身是对还是错。
用高斯做的量化计算需要注意的是量化计算中,或者说薛定谔方程中根本没有“键”这种东西。一切都是波函数。也就是一个一个的分子轨道。优化出来的就是原子间的相对位置,因为薛定谔方程需要的只是原子的位置信息。
而gaussview在显示分子的时候会根据一个内在的键长标准(gaussview内置了一套原子的半径信息,会拿原子间的距离相应的两个原子的半径之和作比较,比如键长超过了原子半径之和的百分之多少就属于单键,之类之类的。),完全是根据原子间的距离来显示键的种类的。跟电子真正在分子中的排布一毛钱关系都没有。所以要分析原子间的成键需要进行波函数分析。算算键级之类的去判断。gaussview显示的完全就是为了好看,但是没有任何化学意义的。并且也是可以自己手动更改键的种类的。
,