24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 867  |  回复: 2

minatomirai

铁虫 (初入文坛)

[求助] ONIOM优化,l9999后restart报错

ONIOM B3LYP/6-31G*优化一个大分子(约7000原子),两层,high layer(QM)约有70原子,是一个较大的体系。坐标用的是Z-Matrix。
使用的关键词如下

# oniom(b3lyp/6-31g(d):amber=hardfirst) nosymm opt=(ts,calcfc,noeigentest,maxcycles=5000) geom=connectivity iop(2/15=3) test

最终以l9999报错,原因是没有收敛。NStep=444(即做了444次),这对于很多做Gaussian优化来说好像步数已经太多了,不过做那么大体系的ONIOM,用的又是B3LYP,还真没底。查看能量,是呈一个下降的趋势,然后判断收敛的几个数值也是朝着减小的方向走。自己认为再加一些循环次数就能够收敛。

然后想restart,我用的Gaussian 09。我直接把上面那行关键词改成

# oniom(b3lyp/6-31g(d):amber=hardfirst) nosymm opt=restart geom=connectivity iop(2/15=3) test

当然,chk文件保持不变。
结果……就报错了:
All quantities printed in internal units (Hartrees-Bohrs-Radians)
MOMM, but no IMMCRS
Error termination via Lnk1e in /home/zjw/g09/l103.exe at Wed Feb 15 14:15:56 2012.

然后上网搜索这个错误……无果。也有人提问过,同是做ONIOM的,但就是没有人回答。
问问这里有人碰到过吗?有人知道是怎么回事吗?
或者,这里有做大分子体系ONIOM的人吗?

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

minatomirai

铁虫 (初入文坛)

木有人碰到过类似的错误?
有人做过restart么?
2楼2012-02-16 22:46:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mophyworld

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
minatomirai(金币+2): 有帮助 好的,非常感谢这条建议! 2012-02-19 06:29:34
没做过ONIOM,但是之前做优化时,用restart一开始计算就出错,在网上搜索过后,听说是G09的restart任务需要有.chk和rwf文件才可以进行。所以我在之后做优化时都会把rwf文件保存下来。但是,后来做任务时,出现优化步数已经达到最大值,使用opt(restart, maxcycle=***)关键字又出现RdWrOT: Data mismatch on MaxStp/MaxJob问题,即设定的maxcycle超过程序给出的maxcycle最大值的问题,后来的解决方式是提取之前优化出现错误时最后一次计算的结构,然后重新开始优化计算而不是使用opt(restart),结果在继续优化后不久结构优化即完成了。

以上是一点经验,希望对您有所帮助。
Ineedsomeinspirationandsomeconfidence
3楼2012-02-18 10:31:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minatomirai 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336求调剂 +5 kiyy 2026-04-01 5/250 2026-04-02 07:14 by chixmc
[考研] 085801 总分275 本科新能源 求调剂 +10 bradoner 2026-04-01 11/550 2026-04-01 22:57 by 来看流星雨10
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[考研] 085601材料工程找调剂 +20 oatmealR 2026-03-29 21/1050 2026-04-01 21:00 by lijunpoly
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 314求调剂 +4 溪云珂 2026-03-26 4/200 2026-04-01 17:00 by oooqiao
[考研] 348求调剂 +5 小懒虫不懒了 2026-03-27 6/300 2026-04-01 16:13 by yanflower7133
[考研] 085900土木水利336分求调剂 +3 Zhangjiangj 2026-03-31 5/250 2026-04-01 09:14 by Zhangjiangj
[考研] 311(085601)求调剂 +12 liziyeyeye 2026-03-28 13/650 2026-04-01 00:34 by fmesaito
[考研] 085602 307分 求调剂 +10 不知道叫什么! 2026-03-26 10/500 2026-03-31 19:53 by Dyhoer
[考研] 08开头看过来!!! +3 wwwwffffff 2026-03-31 5/250 2026-03-31 17:45 by 星光/
[考研] 求化学调剂 +12 wulanna 2026-03-28 12/600 2026-03-31 16:38 by 690616278
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
[考博] 材料专业申博 +5 杜雨婷dyt 2026-03-29 5/250 2026-03-31 11:19 by oooqiao
[考研] 085602化工求调剂(331分) +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
信息提示
请填处理意见