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

yflchx

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

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

KiritsuguPapa

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

kekexili_08

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

liyuanhe211

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

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

liyuanhe211

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

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

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

498746012

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

liyuanhe211

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

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

498746012

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

mengzi

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

huangtf

请问各位老师，SCF=NoVarAcc和SCF=XQC可以同时使用吗？在#p添加SCF(NoVarAcc,XQC,MaxConventionalCyc=N)。谢谢！

huangtf

卡开发发 发表于 2016-5-25 07:09
其实如果知道Helleman-Feynman定理的话，这个也很容易想通。

不得不吐个槽，vasp默认就是类似IOp(5/13 ...

卡老师，您好！请问怎样关闭“vasp默认就是类似IOp(5/13=1)这样的效果”？  谢谢！

卡开发发

huangtf 发表于 2024-12-10 00:47
卡老师，您好！请问怎样关闭“vasp默认就是类似IOp(5/13=1)这样的效果”？  谢谢！

较新的版本你或许可以使用stop_on='Alert'来达到这个目的，当然程序其实还是会输出力等信息。早期版本我们这边是通过自己在main部分更改一些代码来实现lstop或labort的。