计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

16#
 楼主 Author| 发表于 Post on 2015-1-30 00:24:33 | 只看该作者 Only view this author
虎王 发表于 2015-1-29 22:14
如果优化加了int=ultrafine,是否计算能量、轨道、跃迁等都需要加int=ultrafine?

没这个必要
北京科音自然科学研究中心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

67

帖子

0

威望

1047

eV
积分
1114

Level 4 (黑子)

17#
发表于 Post on 2015-6-3 20:43:10 | 只看该作者 Only view this author
本帖最后由 nkallwar 于 2015-6-3 20:44 编辑

最近总出现这样的问题

两个displacement判据早早就变成了0
然后两个force判据就不动了,最后link 9999 错误


         Item               Value     Threshold  Converged?
Maximum Force            0.022910     0.000450     NO
RMS     Force            0.002915     0.000300     NO
Maximum Displacement     0.000000     0.001800     YES
RMS     Displacement     0.000000     0.001200     YES

****
Item               Value     Threshold  Converged?
Maximum Force            0.022910     0.000450     NO
RMS     Force            0.002915     0.000300     NO
Maximum Displacement     0.000000     0.001800     YES
RMS     Displacement     0.000000     0.001200     YES

*****
后面都是这样

这种情况是什么原因造成的?应该用什么方法帮助收敛呢?

opt.JPG (58.31 KB, 下载次数 Times of downloads: 142)

opt.JPG

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

18#
 楼主 Author| 发表于 Post on 2015-6-3 20:53:37 | 只看该作者 Only view this author
nkallwar 发表于 2015-6-3 20:43
最近总出现这样的问题

两个displacement判据早早就变成了0

gdiis和calcall都试试。
注意考察优化过程结构的变化
北京科音自然科学研究中心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

24

帖子

0

威望

54

eV
积分
78

Level 2 能力者

19#
发表于 Post on 2016-2-25 09:06:35 | 只看该作者 Only view this author
本帖最后由 dxyan 于 2016-2-25 18:16 编辑

我的体系就是这样,优化算完后4个YES,然后算频率第三项就变成了NO,而且奇怪的是我每个体系算下来都是这样的。昨天按照楼主说的方式用了opt=tight int=ultrafine由在之前的最后结构上重新做了优化,优化后依然4个YES,但是算出的能量比之前稍微高了一些。然后用这个优化的结果重新做了光计算频率的输入文件,此时没选tight,结果频率计算下来反而后两项都变成了NO,而且还出现了一个虚频。不明白为什么啊?

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

20#
 楼主 Author| 发表于 Post on 2016-2-25 18:28:20 | 只看该作者 Only view this author
dxyan 发表于 2016-2-25 09:06
我的体系就是这样,优化算完后4个YES,然后算频率第三项就变成了NO,而且奇怪的是我每个体系算下来都是这样 ...

你优化时候用opt=tight int=ultrafine,在freq的时候就必须也得用int=ultrafine
北京科音自然科学研究中心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

24

帖子

0

威望

54

eV
积分
78

Level 2 能力者

21#
发表于 Post on 2016-2-25 20:38:46 | 只看该作者 Only view this author
sobereva 发表于 2016-2-25 18:28
你优化时候用opt=tight int=ultrafine,在freq的时候就必须也得用int=ultrafine

注意:如果写opt=tight freq,则freq最后也会用tight标准来判断是否收敛。实际上只要能满足默认判断标准即可,不需要非得满足tight的判断标准
------------------------------------------那前面说算频率不小用tight啊,这算频率为什么要加int=ultrafine呢?

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

22#
 楼主 Author| 发表于 Post on 2016-2-25 22:15:29 | 只看该作者 Only view this author
dxyan 发表于 2016-2-25 20:38
注意:如果写opt=tight freq,则freq最后也会用tight标准来判断是否收敛。实际上只要能满足默认判断标准 ...


你几何优化用什么积分格点,频率计算就必须用什么积分格点,否则不是同一个级别
“算频率不小用tight啊”是病句,理解不能
北京科音自然科学研究中心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

24

帖子

0

威望

54

eV
积分
78

Level 2 能力者

23#
发表于 Post on 2016-2-26 08:43:53 | 只看该作者 Only view this author
sobereva 发表于 2016-2-25 22:15
你几何优化用什么积分格点,频率计算就必须用什么积分格点,否则不是同一个级别
“算频率不小用tight ...

打错字了,不好意思。。我意思是之前帖子说算频率不需要非得满足tight的判断标准,我理解就是算频率不需要加上tight。因为之前帖子说opt=tight,必须加上int=ultrafine。所以后来算频率这2个都没加。

24

帖子

0

威望

54

eV
积分
78

Level 2 能力者

24#
发表于 Post on 2016-3-5 21:35:14 | 只看该作者 Only view this author
Shannon 发表于 2015-1-17 05:21
有没办法让opt=calcall每一步都输出Frequency呢? 如果让gaussian每步都算hessian却不计算出每一步的Freq, ...

请问如何在gaussianview 里面看每一步虚频的方向,并且人工矫正啊?

24

帖子

0

威望

54

eV
积分
78

Level 2 能力者

25#
发表于 Post on 2016-3-5 21:52:43 | 只看该作者 Only view this author
sob老师,请教一下,比如我想优化A,B两个分子,A由于结构优化4个YES,再算频率时出现了2个NO,然后用opt=tight int=ultrafine重新优化了结构,频率计算的时候只加了int=ultrafine来解决这个问题。这样算下来,能量的绝对值比不加opt=tight int=ultrafine时大。那我要比较A,B结构的能量,如果B没有出现出现opt收敛而freq不收敛的情况,是不是B也要用opt=tight int=ultrafine重新优化,并且加int=ultrafine再做频率计算。

另外我要算A和C,B和C的相互作用能,那么C也要用opt=tight int=ultrafine优化吗?AC和BC的复合物在优化的时候是否也要用opt=tight int=ultrafine?谢谢!

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

26#
 楼主 Author| 发表于 Post on 2016-3-6 02:13:23 | 只看该作者 Only view this author
dxyan 发表于 2016-3-5 21:35
请问如何在gaussianview 里面看每一步虚频的方向,并且人工矫正啊?


没法获得每一步虚频方向。只有freq任务处理的那个结构才能观看虚频模式。
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

27#
 楼主 Author| 发表于 Post on 2016-3-6 02:17:00 | 只看该作者 Only view this author
dxyan 发表于 2016-3-5 21:52
sob老师,请教一下,比如我想优化A,B两个分子,A由于结构优化4个YES,再算频率时出现了2个NO,然后用opt=ti ...


只要比较相互的能量、求差值,积分格点必须完全一致,int=ultrafine要加就都加。

其它分子倒不一定非得opt=tight,当前情况是否用tight收敛限对结构和能量的影响只是极为轻微的。不过若你是完美主义,还是都用tight为宜。
北京科音自然科学研究中心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

24

帖子

0

威望

54

eV
积分
78

Level 2 能力者

28#
发表于 Post on 2016-3-6 17:38:41 | 只看该作者 Only view this author
sobereva 发表于 2016-3-6 02:17
只要比较相互的能量、求差值,积分格点必须完全一致,int=ultrafine要加就都加。

其它分子倒不一定 ...

哦哦,那就是出现了一个用了int=ultrafine,其他都要加。谢谢SOb老师。

39

帖子

0

威望

247

eV
积分
286

Level 3 能力者

29#
发表于 Post on 2016-4-20 22:14:49 | 只看该作者 Only view this author
学习了,谢谢sob老师

3

帖子

0

威望

14

eV
积分
17

Level 1 能力者

30#
发表于 Post on 2016-6-29 01:59:49 | 只看该作者 Only view this author
学习了

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 15:19 , Processed in 0.215305 second(s), 28 queries , Gzip On.

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