24小时热门版块排行榜    

查看: 4185  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
请填处理意见