| 查看: 843 | 回复: 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人回复
-
七嗪类物质合成求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有210人回复
-
东方理-中科大联合博士生招聘
已经有0人回复
-
钴酸锂半电池小倍率容量上不去
已经有1人回复
-
吉林大学材料物理本科生求问调剂信息
已经有23人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
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