24小时热门版块排行榜    

查看: 2676  |  回复: 10

791618275

捐助贵宾 (正式写手)

[求助] 出现l502错误,虫友提示了我很多种做法,我不是特别懂做法的含义,望大神指教,谢谢 已有3人参与

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RM06L) =  -4577.14687025     A.U. after  129 cycles
            NFock=128  Conv=0.16D-05     -V/T= 2.1661
KE= 3.925093620919D+03 PE=-4.561870577694D+04 EE= 1.953551716892D+04
Convergence failure -- run terminated.
Error termination via Lnk1e in /share/apps/g09.D01/g09/l502.exe at Mon Sep 11 21:54:09 2017.
Job cpu time:       1 days 15 hours 28 minutes 32.4 seconds.
File lengths (MBytes):  RWF=    629 Int=      0 D2E=      0 Chk=     57 Scr=      1

虫友们给了很多解决方法,我分不太清楚,是不是很多的方法是一样的啊?最好给我说说关键词的含义!
opt=maxcycle=300
iop(1/8=5)或者设成1
optcyc=500,scfcyc=500
scf=maxcyc
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
791618275: 金币+15, ★★★★★最佳答案 2017-09-13 16:20:34
你给出的这几种方法,我不知道是谁建议你的,5种写法其实只有3位方法, 还有两种我认为是没什么用的。
其中optcyc=opt(maxcyc), scfcyc=scf(maxcyc),这两个方法为增加scf和opt的迭代步数,绝大多数情况下没有用,除非你通过观察优化的势能曲线发现能量确实是一直在呈现降低趋势。
iop(1/8=5)或者设成1的含义是改变步长上限至5或1,默认是30,单位是0.01Bohr,在势能面振荡,位移不大但受力较大的情况下可以考虑使用,防止优化直接跨过了极小点。
因l502的作用是求解SCF,所以发生l502错误时,更多应该考虑的是scf不收敛的问题,对于解决scf不收敛,具体方法具体分析,为避免误人子弟,这里贴出sob大神的神贴,LZ可以自己去细细琢磨
http://sobereva.com/61
2楼2017-09-12 10:06:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ra2ghgzh

木虫 (正式写手)

学痴

【答案】应助回帖

感谢参与,应助指数 +1
这个是典型的电子态不收敛问题。
SCF Done:  E(RM06L) =  -4577.14687025     A.U. after  129 cycles
SCF默认是128次迭代,但是实际跑了129次还不迭代。
NFock=128  Conv=0.16D-05     -V/T= 2.1661
说明收敛情况达到了Conv=0.16D-05,默认要到1D-07以下才收敛。
针对这种情况可以设置:
SCF=MaxCycle=256或者更大,这个是让SCF的迭代次数更多。
SCF=Conver=5或者更小,这个是让SCF的收敛条件放宽。如果下次还是Conv=0.16D-05,那就可以收敛了。
SCF=NoVarAcc,这个是让SCF计算时保持较高的积分精度。
上述三个关键词同时服用效果更佳。
我永远当不了老板
3楼2017-09-12 10:07:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

791618275

捐助贵宾 (正式写手)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2017-09-11 14:06:32
你给出的这几种方法,我不知道是谁建议你的,5种写法其实只有3位方法, 还有两种我认为是没什么用的。
其中optcyc=opt(maxcyc), scfcyc=scf(maxcyc),这两个方法为增加scf和opt的迭代步数,绝大多数情况下没有用,除 ...

谢谢,你说的这个贴我也看过,确实很有用......iop(1/8=5)或者设成1的含义是改变步长上限至5或1,默认是30,单位是0.01Bohr,在势能面振荡,位移不大但受力较大的情况下可以考虑使用,防止优化直接跨过了极小点,这里面怎么看位移大不大和受力大不大呢
4楼2017-09-13 16:20:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

791618275

捐助贵宾 (正式写手)

引用回帖:
3楼: Originally posted by ra2ghgzh at 2017-09-11 14:07:05
这个是典型的电子态不收敛问题。
SCF Done:  E(RM06L) =  -4577.14687025     A.U. after  129 cycles
SCF默认是128次迭代,但是实际跑了129次还不迭代。
NFock=128  Conv=0.16D-05     -V/T= 2.1661
说明收敛情 ...

SCF=Conver=5这个关键词不能随便设置吧,这是放宽收敛标准啊,SCF=NoVarAcc,保持较高积分精度是什么意思呢
5楼2017-09-13 16:22:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

铁杆木虫 (正式写手)

本帖仅楼主可见
6楼2017-09-13 17:44:30
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

ra2ghgzh

木虫 (正式写手)

学痴

引用回帖:
5楼: Originally posted by 791618275 at 2017-09-13 16:22:30
SCF=Conver=5这个关键词不能随便设置吧,这是放宽收敛标准啊,SCF=NoVarAcc,保持较高积分精度是什么意思呢...

如果你实在优化分子结构那我觉得SCF=Conver=5可以一试,因为有时分子结构不太好会使得SCF的收敛情况不好,导致一楼出现问题停止优化。
但是你用了这个关键词就有可能跳过这个不好的结构,继续优化,有可能最终的结构更有利于SCF的收敛。
SCF在迭代时要做一些积分计算,这些积分的精度是不断变化的,SCF=NoVarAcc可以让积分精度不变,维持高积分精度。
至于增大SCF=MaxCycle,sob老师说:“如果一个搞量化的碰见SCF不收敛时第一反应是使用这个办法解决,或者看到某人的输入文件里总是带着比如SCF=maxcyc=500这样的关键词,那么他一定是个菜鸟!”
这方面我不赞同,因为这个关键词是Michel Dupuis教授教我使用的,我想他作为国际量子分子科学院院士应该不是菜鸟吧
我永远当不了老板
7楼2017-09-13 18:42:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

791618275

捐助贵宾 (正式写手)

引用回帖:
7楼: Originally posted by ra2ghgzh at 2017-09-12 22:42:54
如果你实在优化分子结构那我觉得SCF=Conver=5可以一试,因为有时分子结构不太好会使得SCF的收敛情况不好,导致一楼出现问题停止优化。
但是你用了这个关键词就有可能跳过这个不好的结构,继续优化,有可能最终的结 ...

我觉得咱俩可以加个扣扣 ,详谈一下- -791618275
8楼2017-09-13 21:44:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

791618275

捐助贵宾 (正式写手)

引用回帖:
6楼: Originally posted by SYST at 2017-09-12 21:44:30
人生短短数十年,该珍惜,该努力岁月匆匆漫漫,你可以选择错误不断的人生,但青春只有一次,加油吧!!!...

研究马克思的,你怎么跑量化的专题来了- -还复制我的个性签名
9楼2017-09-13 21:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

704865589

铁虫 (初入文坛)

【答案】应助回帖

http://sobereva.com/164  直接看sob的帖子吧 这个是我认为目前最靠谱的帖子 一定能解决你的问题
10楼2017-09-18 10:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 791618275 的主题更新
信息提示
请填处理意见