| 查看: 977 | 回复: 2 | ||
[求助]
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的人吗？ 先感谢。 |
回复此楼» 猜你喜欢
半透性容器进行高湿或者低湿研究
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有294人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2012-02-16 22:46:12
mophyworld
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1282.1
- 散金: 10
- 红花: 4
- 帖子: 459
- 在线: 88.2小时
- 虫号: 631894
- 注册: 2008-10-20
- 性别: MM
- 专业: 理论和计算化学
【答案】应助回帖
感谢参与，应助指数 +1
minatomirai(金币+2): ★有帮助 好的，非常感谢这条建议！ 2012-02-19 06:29:34
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），结果在继续优化后不久结构优化即完成了。 以上是一点经验，希望对您有所帮助。 |
3楼2012-02-18 10:31:41
30