Gaussian官方: IOp(5/13=1) 不应在任何计算中使用

yflchx · 发表于 Post on 2016-5-27 15:58:20

liyuanhe211 发表于 2016-5-25 19:14
我还见过国内某知名大学的SI里光明正大的写着自己用了IOp(5/13=1)

应该学到一定程度后再开始钻研IOP。

KiritsuguPapa · 发表于 Post on 2016-6-10 15:42:53

建议此贴链接收入sob那篇文章末尾的批判部分

kekexili_08 · 发表于 Post on 2016-6-13 11:49:42

本帖最后由 kekexili_08 于 2016-6-13 11:52 编辑

我觉得这个也不能说太绝对了。其实第一次SCF迭代循环的时候不能正确收敛，用IOP（5/13=1) 跳过，以后的几何下， SCF正确收敛，多数情况是没有问题的。其实关键是，用了以后要检查收敛到的态是不是想要的。还有，用了IOP（5/13=1)，没正常收敛也是有提示的，有经验的人看输出文件就知道。
还有IOP（1/8=x），用这个是因为以前的某个版本 opt(maxstep)有bug， gaussiian官方回复建议用的IOP（1/8=x），或许后续版本改过了吧。习惯而已，又不是不对，谈不上不好。

liyuanhe211 · 发表于 Post on 2016-6-13 12:59:09

本帖最后由 liyuanhe211 于 2016-6-14 04:21 编辑

kekexili_08 发表于 2016-6-13 11:49
我觉得这个也不能说太绝对了。其实第一次SCF迭代循环的时候不能正确收敛，用IOP（5/13=1) 跳过，以后的 ...

我个人觉得，

1. 明确清楚IOp(5/13=1)在做什么，可能产生什么后果，并且还知道用了之后即使Normal termination也要去检查优化的最后一步是否收敛、收敛到的电子态对不对的用户，都是资深的明白人。而对于原理清楚的资深用户来说，不仅是各种IOp，就算直接去修改代码也显然是可以的。对非资深的、脑子里对原理一团浆糊甚至一无所知的糊涂虫，直接禁止他们这么做是很合理的。

2. 优化的时候非要用IOp跳过一步/几步才能开始收敛多半说明几何结构的初猜是不合理的，在SCF不收敛的情况下用及其不靠谱甚至毫无物理意义的force去猜测下一步的几何结构包含明显的碰大运成分，如果我上一条说的合理，作为“获准使用IOp(5/13=1)的高级用户”也显然应该知道通过调整几何初猜、或者调整SCF的收敛方法等措施来达成SCF收敛，这些方法解决起来优雅、简单的多（考虑到IOp(5/13=1)算完还需要去检查收敛没有、由于中间步骤force极不准确导致的结构合不合理）。

3. maxstep的bug在D.01里已经修复了。本文邮件里也说过是“能用关键词的就用关键词”，有Bug的时候iop凑合一下，没有了还用它干啥，不论从易读性还是出bug的可能性上这个论点都没有任何问题，否则怎么不把整条路径连带B3LYP/6-31G(d)都用IOp写呢，又不是不对，谈不上不好是不？

kekexili_08 · 发表于 Post on 2016-6-14 16:39:01

实践用证明，对于不涉及到高对称性，高自旋态，绝大多数情况IOp(5/13=1)跳过第一SCF迭代循环最后得到的结果是正确的。对于初学者来说，很多时候先到正确的结果，比知道为什么更重要。结构稍微复杂点的，让初学者凋整出好的初始结构，是很难的，让机器从一个不好的初始结构多算点时间得到好的最终结构，未尝不可。
过多争论这个没有太多意义。不过我想，做的时间长了，应该看一些IOP，以便遇到问题的时候知道自己去找出解决的方法。

liyuanhe211 · 发表于 Post on 2016-6-14 16:46:45

kekexili_08 发表于 2016-6-14 16:39
实践用证明，对于不涉及到高对称性，高自旋态，绝大多数情况IOp(5/13=1)跳过第一SCF迭代循环最后得到的结 ...

实践也证明，初学者拿IOp(5/13=1)优化导致优化不收敛甚至最后拿着一个没收敛的能量用的并不少见，知道去检查的十分之少。实践同时证明，只要构型不离谱、调整SCF关键词解决能解决相当大部分的不收敛问题、不一定非要把构型调的多完美。实践同时证明，对比较正常的分子用力场、半经验预优化一下获得合理的构型几乎没有任何技术门槛。

至于最后一点前面说过了。

498746012 · 发表于 Post on 2016-11-29 21:59:10

liyuanhe211 老师您好！Increase integral accuracy 会不会影响能量？一个体系的其他结构是不是也得用Increase integral accuracy

liyuanhe211 · 发表于 Post on 2016-11-29 23:20:12

498746012 发表于 2016-11-29 21:59
liyuanhe211 老师您好！Increase integral accuracy 会不会影响能量？一个体系的其他结构是不是也得用Incre ...

会。有一个用了int=ultrafine其他的都得用

498746012 · 发表于 Post on 2016-11-30 12:30:22

liyuanhe211 老师您好!能不能麻烦您一下：Integral=(Acc2E=N)会不会影响能量？一个体系的其他结构是不是也得用Integral=(Acc2E=N)？还有opt=tight是不是也是同样的情况？

mengzi · 发表于 Post on 2016-11-30 14:09:33

谢谢Li老师热心而又细致的讲解，受益匪浅。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] Gaussian官方: IOp(5/13=1) 不应在任何计算中使用

评分 Rate