请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 547|回复: 11

[Gaussian/gview] iop(1/8=10) 这个关键词是否有利于收敛?

[复制链接]

424

帖子

0

威望

951

eV
积分
1375

Level 4 (黑子)

发表于 2017-1-11 21:48:24 | 显示全部楼层 |阅读模式
我看了 iop(1/8=k) 这个关键词的相关资料但是里面有很多词不太明白什么意思,我不知道这个关键词在优化的时候是否有利于收敛?谢谢!
图1.png

IOP.txt

14.77 KB, 下载次数: 7

评分

参与人数 1eV +5 收起 理由
zsu007 + 5 谢谢

查看全部评分

1727

帖子

19

威望

5516

eV
积分
7623

Level 6 (一方通行)

发表于 2017-1-12 00:45:59 | 显示全部楼层
不要用iop(1/8=x),用opt=maxstep=x代替。有关键词就不要用iop

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

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

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

9462

帖子

16

威望

1万

eV
积分
21596

管理员

公社社长

发表于 2017-1-12 02:43:58 | 显示全部楼层
这里专门批了有关键词却还写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

思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。共计近4000人,讨论计算化学为主,加入时必须注明研究方向,否则一律不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

101

帖子

0

威望

255

eV
积分
356

Level 3 能力者

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

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

424

帖子

0

威望

951

eV
积分
1375

Level 4 (黑子)

 楼主| 发表于 2017-1-12 17:03:27 | 显示全部楼层
"用gview检查优化过程的轨迹,看是否最后开始震荡了"  这句话在gview具体是怎么操作的呀?谢谢!

125

帖子

1

威望

1789

eV
积分
1934

Level 5 (御坂)

发表于 2017-1-12 21:17:30 | 显示全部楼层
本帖最后由 bomsaude 于 2017-1-12 13:18 编辑
15whyhq 发表于 2017-1-12 09:03
"用gview检查优化过程的轨迹,看是否最后开始震荡了"  这句话在gview具体是怎么操作的呀?谢谢!


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

评分

参与人数 1eV +1 收起 理由
sobereva + 1

查看全部评分

424

帖子

0

威望

951

eV
积分
1375

Level 4 (黑子)

 楼主| 发表于 2017-1-13 09:26:43 | 显示全部楼层
bomsaude 发表于 2017-1-12 21:17
打开没有优化结束的文件,记着选上读取中间结构的选项,然后用"Results"中的"Optimization"查看。

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

424

帖子

0

威望

951

eV
积分
1375

Level 4 (黑子)

 楼主| 发表于 2017-1-13 10:37:10 | 显示全部楼层
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这个关键词?还是不太明白那段话得意思,谢谢!

9462

帖子

16

威望

1万

eV
积分
21596

管理员

公社社长

发表于 2017-1-13 14:05:00 | 显示全部楼层
15whyhq 发表于 2017-1-13 10:37
#p opt=(calcfc,ts,noeigentest) oniom(b3pw91/6-31g(d,p):hf) freq(noraman)   老师,我不知道我的输入 ...

trustupdate毫无意义
要加都是加notrust,而对于找TS时本来就是默认的
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。共计近4000人,讨论计算化学为主,加入时必须注明研究方向,否则一律不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

125

帖子

1

威望

1789

eV
积分
1934

Level 5 (御坂)

发表于 2017-1-13 17:35:04 | 显示全部楼层
本帖最后由 bomsaude 于 2017-1-13 09:48 编辑
15whyhq 发表于 2017-1-13 01:26
“选上读取中间结构的选项”  是怎么操作的呢?谢谢!


426

帖子

11

威望

2694

eV
积分
3340

Level 5 (御坂)

密度泛函·小卒

发表于 2017-1-14 09:25:14 | 显示全部楼层
本帖最后由 yjcmwgk 于 2017-1-14 09:26 编辑

看到楼主和sob君、老李等人的“追问型互动”,就觉得楼主真该通读一遍Exp Chem with Elect Struct Methods再来提问
这本书的第二版,盗版满天飞,随便下载
东风快递,使命必达!六分钟覆盖全球,哪里不服点哪里!

9462

帖子

16

威望

1万

eV
积分
21596

管理员

公社社长

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


那本书不解决问题,群里各种日经、半周经、周经、半月经、月经...问题在书里都不会去介绍怎么以最好的方式解决
其实他问的>90%的问题都是科音量化初级班里都会讲的内容,不系统地好好学一次,诸如查看收敛趋势这种初级问题碰见一个问一个真是很麻烦,效率极低
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。共计近4000人,讨论计算化学为主,加入时必须注明研究方向,否则一律不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

Archiver|手机版|小黑屋|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2017-6-26 16:53 , Processed in 0.124930 second(s), 26 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表