24小时热门版块排行榜    

查看: 981  |  回复: 6

xiaowu787

木虫 (正式写手)


[交流] GAMESS 收敛失败

GAMESS 收敛失败
收敛的标准不算很严格吧,怎么就收敛失败呢?
CODE:
$contrl scftyp=uhf runtyp=optimize mult=2
         icharg=0 coord=zmt mplevl=2 nzvar=0 $end
$system timlim=2 memory=2000000 $end
$statpt opttol=1.0e-5 $end
$basis  gbasis=n21 ngauss=3 $end
$guess  guess=huckel $end
$data
Cu-O
C1
O
Cu 1 r

r=1.73
$end

CPU     0: STEP CPU TIME=     0.01 TOTAL CPU TIME=        0.2 (    0.0 MIN)
TOTAL WALL CLOCK TIME=        1.6 SECONDS, CPU UTILIZATION IS  11.18%

SCF DID NOT CONVERGE...NO MPLEVL=2 CALCULATION
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Tue Nov 29 20:11:28 2011
                580000  WORDS OF DYNAMIC MEMORY USED
CPU     0: STEP CPU TIME=     0.00 TOTAL CPU TIME=        0.2 (    0.0 MIN)
TOTAL WALL CLOCK TIME=        1.7 SECONDS, CPU UTILIZATION IS  10.53%
ddikick.x: application process 0 quit unexpectedly.
ddikick.x: Sending kill signal to DDI processes.
ddikick.x: Execution terminated due to error(s).
----- accounting info -----

[ Last edited by xiaowu787 on 2011-11-30 at 14:18 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787(金币+1):谢谢参与
如果总能量没有任何收敛趋势,说明程序不知道该算哪一种轨道占据和电子态。多参考体系和近简并体系经常会遇到这种情况,一般需要指定电子占据,或者做多参考计算。

如果趋向收敛,加大迭代次数,或者改用QC收敛方案。
2楼2011-11-29 23:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abdoman

木虫 (著名写手)


★ ★ ★ ★
xiaowu787(金币+1):谢谢参与
小红豆(金币+3): 谢谢回帖 2011-11-30 13:40:28
收敛标准opttol=1.0e-5  这个设置有点点严格了吧。默认的是1.0e-04
你的基组又太粗糙,建议降低收敛标准试试。比如用默认标准
3楼2011-11-30 11:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)


引用回帖:
2楼: Originally posted by beefly at 2011-11-29 23:19:03:
如果总能量没有任何收敛趋势,说明程序不知道该算哪一种轨道占据和电子态。多参考体系和近简并体系经常会遇到这种情况,一般需要指定电子占据,或者做多参考计算。

如果趋向收敛,加大迭代次数,或者改用QC收敛 ...

确实是没有收敛趋势,能量一直波动。多参考体系的计算还没做过,能否指点一下?我再查查相关文献
4楼2011-11-30 17:46:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
xiaowu787(金币+2): 谢谢指点 2011-12-03 14:34:44
引用回帖:
4楼: Originally posted by xiaowu787 at 2011-11-30 17:46:28:
确实是没有收敛趋势,能量一直波动。多参考体系的计算还没做过,能否指点一下?我再查查相关文献

可以先用gaussian程序,用同样的计算方法算一下。如果能收敛,可以对照gaussian的轨道占据,在gamess计算中交换轨道,调整为同样的占据。需要注意坐标方位必须一致,不过双原子分子坐标方位都是约定俗成的,不会遇到这种问题。

还有个办法,先用几个杂化泛函试试,或者做闭壳层离子的计算。如果能收敛,从产生的punch文件找到波函,复制到输入文件,改变初始猜测设置,做正常的计算也许就能收敛。
5楼2011-12-01 00:55:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)


用Gaussian的mp2是可以收敛的,对于新手,上面的操作感觉比较复杂呀,谢谢指点!

[ Last edited by xiaowu787 on 2011-12-3 at 14:43 ]
6楼2011-12-03 14:41:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nibelungen

新虫 (初入文坛)



xiaowu787(金币+1): 谢谢参与
statpt group 里的opttol不是SCF的收敛标准,应该是几何优化的收敛标准吧?
楼主可以试试放宽scf group里的conv试试,默认是1.0d-5。不过不知道降低scf的收敛标准会不会让结果偏差很大
7楼2012-06-09 20:23:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowu787 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见