计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12045|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 请教高斯优化过渡态未收敛,振动方式正确的结构算正确吗?

[复制链接 Copy URL]

4

帖子

0

威望

21

eV
积分
25

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
请教一下,用高斯优化过渡态的时候,优化后出现l9999报错,但是振动方式和猜测的相同,而且在log文件中也给出了能量,请问这个结构可以使用吗?
也用过其他方法尝试继续收敛,可是收敛之后优化出的结构和初猜结构相差很大,而且也不是想要的振动方式~
希望前辈们能够帮帮忙,十分感谢!!!!

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2016-12-7 00:21:44 | 只看该作者 Only view this author
如果你是指优化过渡态没收敛,然后用最后一帧的结构直接做振动分析,发现虚频模式是想要的,那么能否用要看freq任务最后对收敛的判断情况。如果受力和位移比默认收敛限大不超过2倍,勉强可以接受,否则必须重新优化过渡态,calcfc/calcall结合maxstep减小步长上限一般都能解决。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

4

帖子

0

威望

21

eV
积分
25

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2016-12-7 08:49:52 | 只看该作者 Only view this author
sobereva 发表于 2016-12-7 00:21
如果你是指优化过渡态没收敛,然后用最后一帧的结构直接做振动分析,发现虚频模式是想要的,那么能否用要看 ...

sob老师,我的是在l9999的时候虚频模式是想要的,而继续用最后优化的结构继续优化则出现了错误的振动模式,使用的是calcall,有两个结构都是这样,这是因为我的初猜错了吗?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2016-12-7 12:51:09 | 只看该作者 Only view this author
阿秋Autumn 发表于 2016-12-7 08:49
sob老师,我的是在l9999的时候虚频模式是想要的,而继续用最后优化的结构继续优化则出现了错误的振动模式 ...

用opt(TS,maxstep=5,calcall)基于之前最后的结构继续优化过渡态,如果确实跑到其它地方去了,那即便之前那个结构虚频正确也白搭

有唯一的虚频且振动方向正确也不代表结构就是过渡态结构,必须结合受力为0的条件(做振动分析时末尾输出的受力接近0)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

4

帖子

0

威望

21

eV
积分
25

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2016-12-8 08:48:25 | 只看该作者 Only view this author
sobereva 发表于 2016-12-7 12:51
用opt(TS,maxstep=5,calcall)基于之前最后的结构继续优化过渡态,如果确实跑到其它地方去了,那即便之前 ...

sob老师,我用您的方法对最后的结构继续进行了优化,可是最后还是出现了l9999的错误,信息如下:
Error termination request processed by link 9999.
Error termination via Lnk1e in /lustre/software/g09/l9999.exe at Thu Dec  8 00:32:58 2016.
Job cpu time:       0 days  7 hours 15 minutes 46.9 seconds.
File lengths (MBytes):  RWF=    285 Int=      0 D2E=      0 Chk=      8 Scr=      1

希望老师指点,感谢~

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

6#
发表于 Post on 2016-12-8 17:23:52 | 只看该作者 Only view this author
阿秋Autumn 发表于 2016-12-8 08:48
sob老师,我用您的方法对最后的结构继续进行了优化,可是最后还是出现了l9999的错误,信息如下:
Error ...

别光看l9999,要用gview观看优化过程结构和能量、受力变化趋势,判断怎么处理
凡是可能有用的做法都在这里提了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

4

帖子

0

威望

21

eV
积分
25

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2016-12-8 18:09:51 | 只看该作者 Only view this author
sobereva 发表于 2016-12-8 17:23
别光看l9999,要用gview观看优化过程结构和能量、受力变化趋势,判断怎么处理
凡是可能有用的做法都在这 ...

好的,谢谢sob老师~

142

帖子

3

威望

4379

eV
积分
4581

Level 6 (一方通行)

8#
发表于 Post on 2017-10-10 00:05:52 | 只看该作者 Only view this author
sobereva 发表于 2016-12-7 12:51
用opt(TS,maxstep=5,calcall)基于之前最后的结构继续优化过渡态,如果确实跑到其它地方去了,那即便之前 ...

Sobereva老师,请问什么情况下会出现TS优化出只有一个虚频且是反应方向上的虚频,而freq时的Forces远大于0的情况呢?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

9#
发表于 Post on 2017-10-10 00:15:54 | 只看该作者 Only view this author
laoman 发表于 2017-10-10 00:05
Sobereva老师,请问什么情况下会出现TS优化出只有一个虚频且是反应方向上的虚频,而freq时的Forces远大于 ...

不会出现这种情况。force远大于0的情况opt任务根本不会正常结束。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

142

帖子

3

威望

4379

eV
积分
4581

Level 6 (一方通行)

10#
发表于 Post on 2017-10-10 06:05:58 | 只看该作者 Only view this author
sobereva 发表于 2017-10-10 00:15
不会出现这种情况。force远大于0的情况opt任务根本不会正常结束。

明白了,感谢社长的回答。最近发现原来从freq的chk(来自PES扫描后的最高点,有反应方向的虚频,-1600到-1800的范围)读取geom和guess会比导出坐标算的快了一点(虽然也用了guess),读取chk的前提应该是freq和后面的ts都用的同样的基组吧。

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

11#
发表于 Post on 2017-10-10 14:01:43 | 只看该作者 Only view this author
laoman 发表于 2017-10-10 06:05
明白了,感谢社长的回答。最近发现原来从freq的chk(来自PES扫描后的最高点,有反应方向的虚频,-1600到- ...


快一点是因为从chk里读取了初猜,导致SCF部分基本不耗什么时间。
guess=read前后的基组可以不同,程序自动会做投影
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

142

帖子

3

威望

4379

eV
积分
4581

Level 6 (一方通行)

12#
发表于 Post on 2017-10-10 15:38:50 | 只看该作者 Only view this author
sobereva 发表于 2017-10-10 14:01
快一点是因为从chk里读取了初猜,导致SCF部分基本不耗什么时间。
guess=read前后的基组可以不同,程序 ...

好的,谢谢社长~

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 10:51 , Processed in 0.176883 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list