计算化学公社

标题: iop(1/8=10) 这个关键词是否有利于收敛？ [打印本页]

作者
Author: 15whyhq 时间: 2017-1-11 21:48
标题: iop(1/8=10) 这个关键词是否有利于收敛？
我看了 iop(1/8=k) 这个关键词的相关资料但是里面有很多词不太明白什么意思，我不知道这个关键词在优化的时候是否有利于收敛？谢谢！

作者
Author: liyuanhe211 时间: 2017-1-12 00:45
不要用iop(1/8=x)，用opt=maxstep=x代替。有关键词就不要用iop

在优化的时候一般不必动它，将其缩小通常会使得几何优化步骤增加；但如果遇到优化震荡，使其缩小至3~10可能打破震荡使得几何优化收敛。

我个人觉得增大maxstep的数值一般也没有必要，我在我跑过的文件里搜了一下（包括很多极小点和过渡态的优化），往往超过trustradius，需要scale down的情况，在每个优化里最多也就5步左右，以1~3步居多，所以即使调大了maxstep，也最多不过省个一两步而已，还增加了不收敛的风险。

对TS优化默认就是notrustradius我还是第一次知道，感谢提及

作者
Author: sobereva 时间: 2017-1-12 02:43
这里专门批了有关键词却还写IOp的行为
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

靠恰当利用maxstep帮助解决不收敛这里明说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

想了解更多几何优化知识看
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=105

不要图省事在优化一开始用这个减小步长上限，否则收敛所需步数往往多几甚至十几倍。一定要等到发现震荡了再用这个减小步长上限来解决震荡。

PS：应当对txt文档的内容注明出处：
Gaussian中有用的IOp一览
http://sobereva.com/93

作者
Author: wswrpd 时间: 2017-1-12 05:03

sobereva 发表于 2017-1-12 02:43
这里专门批了有关键词却还写IOp的行为
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

我一个电子计算领域小白，还是得多听sob老大的一点人生经验比较好

作者
Author: 15whyhq 时间: 2017-1-12 17:03
"用gview检查优化过程的轨迹，看是否最后开始震荡了" 这句话在gview具体是怎么操作的呀？谢谢！

作者
Author: bomsaude 时间: 2017-1-12 21:17
本帖最后由 bomsaude 于 2017-1-12 13:18 编辑

15whyhq 发表于 2017-1-12 09:03
"用gview检查优化过程的轨迹，看是否最后开始震荡了" 这句话在gview具体是怎么操作的呀？谢谢！

打开没有优化结束的文件，记着选上读取中间结构的选项，然后用"Results"中的"Optimization"查看。

作者
Author: 15whyhq 时间: 2017-1-13 09:26

bomsaude 发表于 2017-1-12 21:17
打开没有优化结束的文件，记着选上读取中间结构的选项，然后用"Results"中的"Optimization"查看。

“选上读取中间结构的选项” 是怎么操作的呢？谢谢！

作者
Author: 15whyhq 时间: 2017-1-13 10:37

sobereva 发表于 2017-1-12 02:43
这里专门批了有关键词却还写IOp的行为
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

#p opt=(calcfc,ts,noeigentest) oniom(b3pw91/6-31g(d,p):hf) freq(noraman) 老师，我不知道我的输入文件里面的关键词是不是应该加上trustupdate这个关键词？还是不太明白那段话得意思，谢谢！

作者
Author: sobereva 时间: 2017-1-13 14:05

15whyhq 发表于 2017-1-13 10:37
#p opt=(calcfc,ts,noeigentest) oniom(b3pw91/6-31g(d,p):hf) freq(noraman) 老师，我不知道我的输入 ...

trustupdate毫无意义
要加都是加notrust，而对于找TS时本来就是默认的

作者
Author: bomsaude 时间: 2017-1-13 17:35
本帖最后由 bomsaude 于 2017-1-13 09:48 编辑

15whyhq 发表于 2017-1-13 01:26
“选上读取中间结构的选项” 是怎么操作的呢？谢谢！

作者
Author: yjcmwgk 时间: 2017-1-14 09:25
本帖最后由 yjcmwgk 于 2017-1-14 09:26 编辑

看到楼主和sob君、老李等人的“追问型互动”，就觉得楼主真该通读一遍Exp Chem with Elect Struct Methods再来提问

这本书的第二版，盗版满天飞，随便下载

作者
Author: sobereva 时间: 2017-1-14 16:39

yjcmwgk 发表于 2017-1-14 09:25
看到楼主和sob君、老李等人的“追问型互动”，就觉得楼主真该通读一遍Exp Chem with Elect Struct Methods ...

那本书不解决问题，群里各种日经、半周经、周经、半月经、月经...问题在书里都不会去介绍怎么以最好的方式解决
其实他问的>90%的问题都是科音量化初级班里都会讲的内容，不系统地好好学一次，诸如查看收敛趋势这种初级问题碰见一个问一个真是很麻烦，效率极低

欢迎光临计算化学公社 (http://bbs.keinsci.com/)