| 查看: 869 | 回复: 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的人吗？ 先感谢。 |
回复此楼» 猜你喜欢
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
物理化学论文润色/翻译怎么收费?
已经有180人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
调剂招生
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
2楼2012-02-16 22:46:12
