24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4238  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-05 4/200 2026-04-05 20:29 by 啵啵啵0119
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 385分 生物学(071000)求调剂 +11 qf626 2026-04-01 11/550 2026-04-05 17:35 by Ecowxq666!
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 材料调剂 +18 一样YWY 2026-04-02 19/950 2026-04-04 22:14 by hemengdong
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 298求调剂 +5 zzz,,r 2026-04-02 8/400 2026-04-04 19:55 by 蓝云思雨
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 一志愿华东理工大学,080500学硕,317分,求调剂 +13 s1145 2026-03-31 15/750 2026-04-03 11:44 by msi123
[考研] 一志愿北京科技材料科学与工程288分,求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 085601一志愿西北工业大学初试346 +4 085601初试346 2026-03-30 4/200 2026-03-31 07:47 by jp9609
信息提示
请填处理意见