24小时热门版块排行榜    

查看: 3990  |  回复: 12

calebby

捐助贵宾 (著名写手)

[求助] gaussian能量震荡的问题 已有1人参与

gaussian优化配合物出现L9999错误,输入文件如下:
%chk=EGME.chk
%mem=8GB
%nprocshared=8
# opt(cartesian) rb3lyp/genecp

3D Atomistic

0 1
*****(原子坐标省略)
C H O 0
6-31g(d)
****
Zn 0
Lanl2dz
****

Zn 0
Lanl2dz
查看log文件发现原因是不收敛,如下:
        Item               Value     Threshold  Converged?
Maximum Force            0.010194     0.000450     NO
RMS     Force            0.002456     0.000300     NO
Maximum Displacement     0.014784     0.001800     NO
RMS     Displacement     0.006932     0.001200     NO
Predicted change in Energy=-2.177814D-03
Optimization stopped.
    -- Number of steps exceeded,  NStep= 241
    -- Flag reset to prevent archiving.
之前也增加了循环次数,还是不收敛。所以现在怀疑是能量震荡,但是log文件用GV打开后,results--optimization按钮是灰色的,点不了。不知道有没有什么办法确认是能量震荡?以及这个错误如何处理呢?
按照论坛上的说法,直接提取最后一步的构型重新优化,试了一次还是不行,现在继续提取最后一步,接着优化中。。。
回复此楼

» 收录本帖的淘帖专辑推荐

量子化学和计算化学

» 猜你喜欢

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

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

绍敏郡主

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
gmy1990: 金币+1 2013-12-17 20:50:57
看log文件中每一个循环之后E的变化趋势
2楼2013-11-24 20:48:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
如LS所说,看下每步E的趋势,另外可用opt=gdiis,据说是解决震荡神器
3楼2013-11-25 07:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)

本帖仅楼主可见
4楼2013-12-09 16:15:23
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

javy

金虫 (正式写手)

【答案】应助回帖

★ ★
gmy1990: 金币+2 2013-12-17 20:51:08
可以用gausssum软件来看,这个比较方便看优化结果是否能量振荡。如果是能量在振荡的话,首先可以减小优化步长,IOP(1/8=n),默认是30,你可以选择更小的数,一般可以收敛。opt=GDIIS在减小步长不收敛的情况下用一般会管用,不过花的时间相对更长。~
5楼2013-12-09 17:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)

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

雪狼乖乖

木虫 (正式写手)


【答案】应助回帖

★ ★
gmy1990: 金币+2 2013-12-19 08:32:33
有两个问题:
1 你确定这个计算是闭壳层计算?
2 是否是结构不合理造成的。一般L9999能量出现震荡,gdiis未解决的 都是结构问题。
7楼2013-12-17 21:11:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marson

银虫 (正式写手)

引用回帖:
4楼: Originally posted by 四毛@ at 2013-12-09 16:15:23
G09好像默认的就是gdiis
试过了,木有用。。。。...

G09默认的是gediis吧,不是gdiis~
思想有多远,我们就能走多远~
8楼2013-12-18 15:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)

本帖仅楼主可见
9楼2013-12-18 17:31:24
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
gmy1990: 金币+2 2013-12-19 08:32:43
引用回帖:
9楼: Originally posted by 四毛@ at 2013-12-18 17:31:24
现在我后两个判据都yes了,但是前两个整死都是no,opt(maxcyc=500)了,还在算。不知道能不能行,做的是激发态结构的优化opt td。...

同意7楼的说法,先检查电荷、自旋是否正确,然后判断合理的结构。

如果上述都没问题,换一种方法或是泛函作为预优化,以及适当提升基组也有可能能够解决。
不一定挂在论坛,计算问题问题欢迎留言。
10楼2013-12-18 23:20:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 calebby 的主题更新
信息提示
请填处理意见