YuYuqi的个人博客分享 http://blog.sciencenet.cn/u/YuYuqi

博文

ONIOM: Microiterations cycle 26 out of a maximum of 25

已有 3013 次阅读 2018-2-4 12:39 |个人分类:量化|系统分类:科研笔记

用ONIOM优化时遇到:

Microiterations cycle 26 out of a maximum of 25 ,maximum number of microiterations cycles exceeded!

解决方法:


1.

可以去一部分水再试。

另外整体结构可以先在UFF下优化,然后再尝试改用ONIOM;或者先尝试不用电子嵌入,优化成功了再开启

(来自计算化学公社)

2.

Actually I have a solution for this problem. Strangely, it seems there is no keyword to change the maximum number

of microiteration cycles. But there is a way to get around it. Use the keywords: " opt=quadmacro iop(1/98=66,1/19=7)"

Useful options might be iop(1/8=3,1/13=6). The deal is the following: "quadmacro" stands for quadratic macrostep.

This is an awesome feature in g09, I heavily encourage everybody to use it in all QM/MM calculations, especially for

TS searches. Iop(1/98=66) makes a twist in the optimization workflow, it disables microiterations, only second order steps

are made everytime in both the QM and MM regions. Keep in mind that the MM region needs much more steps to

converge than the QM region because of its sheer size. If you run an optimization like this, you will notice that it will take

steps faster than the standard protocol, and eventually the QM region will reach convergence, but the MM region will

most probably be far from convergence. So after like a few tens of steps (for example 50-100) stop, get the latest

structure, and start an optimization with microiterations. Mine looks like:

"opt=(quadmacro,readfc) iop(1/8=3,1/128=1,1/110=2,1/19=7,1/13=6)". 1/128 and 1/110 are microiteration realted options

I think they are useful too in a lot of cases. 1/19=7 is a must (at least for my calculations) with quadmacro because of a

bug. If you get still the same error, than get another round of pire quadmacro optimization. And try again microiterations

afterwards. You can reach full convergence without microiterations, I suppose, but it would be quite inefficient, I think.  

(from Researchgate Gergely Bánóczi




https://blog.sciencenet.cn/blog-3375268-1098402.html

上一篇:选择小分子附近的完整残基
下一篇:word2016合并单元格不显示内容也不能编辑的问题
收藏 IP: 202.127.19.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-25 18:14

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部