当前位置: 首页 > 量子化学 >高斯不收敛,如何优化

高斯不收敛,如何优化

作者 sky423
来源: 小木虫 950 19 举报帖子
+关注

#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 返回小木虫查看更多

今日热帖
  • 精华评论
  • paramecium86

    引用回帖:
    9楼: Originally posted by sky423 at 2021-06-22 08:33:03
    仔细检查完之后发现结果的确有问题,所以重新构建了分子,但是新的分子却压根无法导入gaussian进行计算,报错,如下图分别是03和09的报错


    ...

    看起来是高斯的Scratch文件夹里内容被占用了。你修改好之后保存成另外一个名字的输入文件。并且把高斯的scratch文件夹内容清空 再计算就可以了。

  • paramecium86

    引用回帖:
    8楼: Originally posted by sky423 at 2021-06-22 05:40:46
    是含有重金属铼
    ...

    包含重金属,也不一定非要所有原子都用LanL2DZ 可以用混合基组 给重原子用LanL2DZ(这个2021年的今天已经很不推荐了精度太差,建议SDD) 给其他轻原子(前三周期)用普通的基组比如 6-31G(D)之类的。会合理很多。

  • paramecium86

    引用回帖:
    13楼: Originally posted by sky423 at 2021-06-25 01:53:14
    第六周期的重金属,直接把LanL2DZ换成SDD就可以了吧?别的什么关键词修饰还需要吗?
    ...

    之前那个任务直接改不太行的。要改成混合基组需要特别的写法。见我这个教程1.1节 https://muchong.com/t-14242735-1

  • paramecium86

    引用回帖:
    15楼: originally posted by sky423 at 2021-06-28 00:44:00
    你好,我能再麻烦您问您一个问题吗?就是我输入文件结构是正常的,每个键都连接很好,并且中间也没有死点,完成了,但是结构完成之后部分键就断开了,这属于正常现象吗?我就单纯的是普通优化,这样完全有背我原始 ...

    如果基组设置没有问题。电荷和多重度也设置的没有什么问题。结构本身原子间的距离并没有太大异常而仅仅是gaussview不显示某些键那一般不会有什么问题。
    见我这个教程3.3节 https://muchong.com/t-14242735-1

    3.3 化学键的显示问题。
    诸如 为什么优化完成后“缺了什么什么键“,“哪个键变成了单/双/三键”,“那个键不见了”等问题。
    这类问题的回答是:没有关系。gaussview怎么显示的,都没什么所谓。键显示的符合预期也好或者不符合预期也罢都不能说明计算结果本身是对还是错。
    用高斯做的量化计算需要注意的是量化计算中,或者说薛定谔方程中根本没有“键”这种东西。一切都是波函数。也就是一个一个的分子轨道。优化出来的就是原子间的相对位置,因为薛定谔方程需要的只是原子的位置信息。
    而gaussview在显示分子的时候会根据一个内在的键长标准(gaussview内置了一套原子的半径信息,会拿原子间的距离相应的两个原子的半径之和作比较,比如键长超过了原子半径之和的百分之多少就属于单键,之类之类的。),完全是根据原子间的距离来显示键的种类的。跟电子真正在分子中的排布一毛钱关系都没有。所以要分析原子间的成键需要进行波函数分析。算算键级之类的去判断。gaussview显示的完全就是为了好看,但是没有任何化学意义的。并且也是可以自己手动更改键的种类的。

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓