| 查看: 975 | 回复: 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的人吗？ 先感谢。 |
回复此楼» 猜你喜欢
申请2027年材料/化学类博士
已经有6人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有192人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有63人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
10