应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gaussian能量震荡的问题
131/2返回列表12下一页
查看: 4460  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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按钮是灰色的，点不了。不知道有没有什么办法确认是能量震荡？以及这个错误如何处理呢？
按照论坛上的说法，直接提取最后一步的构型重新优化，试了一次还是不行，现在继续提取最后一步，接着优化中。。。
回复此楼

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

量子化学和计算化学

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-24 20:45:54
已阅   回复此楼   关注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 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 oepdf5z972 2026-06-20 3/150 2026-06-21 21:01 by zo4669e5n4
[基金申请] E0414, 我的本子有没有希望？ +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[博后之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
关闭
请填处理意见
关闭