24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 3320  |  回复: 19
【悬赏金币】回答本帖问题,作者sky423将赠送您 50 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sky423

新虫 (著名写手)

[求助] 高斯不收敛,如何优化已有1人参与

#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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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显示的完全就是为了好看,但是没有任何化学意义的。并且也是可以自己手动更改键的种类的。
16楼2021-06-28 13:59:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

DX0929

至尊木虫 (著名写手)

增加循环次数

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

机会是留给有准备的人的。
2楼2021-06-21 16:38:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sky423

新虫 (著名写手)

送红花一朵
引用回帖:
2楼: Originally posted by DX0929 at 2021-06-21 16:38:33
增加循环次数

#p b3lyp/lanl2dz opt pop=(full,npa) SCF=(maxcyc=200) SCF=QC


类似于这样吗?谢谢,我太菜了,谢谢您

发自小木虫Android客户端
3楼2021-06-21 17:15:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

此楼不小心重复了
4楼2021-06-22 05:30:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见