计算化学公社

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

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

  [复制链接 Copy URL]

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

31#
发表于 Post on 2016-7-5 08:35:50 | 只看该作者 Only view this author
最近在做激发态优化,用的Opt=Vtight, Int=UltraFine,有几个结构总是有一个虚频,怎么都消不掉,试了maxstep=1也不行(Sob也不推荐这个法子)。但是TD计算没有解析二阶导数,无法使用Opt=CalcAll,不知道该如何处理比较合适?是换优化算法吗?

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

32#
 楼主 Author| 发表于 Post on 2016-7-5 09:39:23 | 只看该作者 Only view this author
ter20 发表于 2016-7-5 08:35
最近在做激发态优化,用的Opt=Vtight, Int=UltraFine,有几个结构总是有一个虚频,怎么都消不掉,试了maxst ...

注意做振动分析的时候此时也得用ultrafine

拿最后结构,opt里写readfc读取freq产生的精确Hessian再优化几步
也可以尝试gdiis。
或按照虚频模式位移一点点再用calcfc并继续优化
北京科音自然科学研究中心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

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

33#
发表于 Post on 2016-7-5 11:04:04 | 只看该作者 Only view this author
sobereva 发表于 2016-7-5 09:39
注意做振动分析的时候此时也得用ultrafine

拿最后结构,opt里写readfc读取freq产生的精确Hessian再优 ...

感谢Sob老师的回复,我挨个试试。还有些问题:
1. 我试了CalcFC,会报错,说TD没有解析梯度,我想应该是您的笔误吧
2. 我看了结果文件,里面有Using GEDIIS/GDIIS optimiser这样的语句,应该是已经用了GDIIS算法了
3. 如果关键字opt=(vtight,maxstep=1) Int=UltraFine freq这样肯定是没问题的对吧?你说的Freq时也要用Int=UltraFine是指Freq和Opt分开做是吧

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

34#
 楼主 Author| 发表于 Post on 2016-7-5 12:54:57 | 只看该作者 Only view this author
ter20 发表于 2016-7-5 11:04
感谢Sob老师的回复,我挨个试试。还有些问题:
1. 我试了CalcFC,会报错,说TD没有解析梯度,我想应该是 ...


1 你就freq之后readfc就完了。应该也有选项强行让高斯用数值Hessian。
2 不是,高斯默认是自动把GEDIIS、RFO通过一定规则组合在一起用的,和直接用GDIIS不同。有时候默认的优化方式不收敛,改用GDIIS能解决。
3 没问题,但opt里还应该增加notrust,否则优化过程中程序还是会自动根据情况自动改变步长上限。通常步长上限用3、5就够了,1太费时间。是。
北京科音自然科学研究中心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

177

帖子

0

威望

1567

eV
积分
1744

Level 5 (御坂)

35#
发表于 Post on 2016-7-5 13:15:48 | 只看该作者 Only view this author
sobereva 发表于 2016-7-5 12:54
1 你就freq之后readfc就完了。应该也有选项强行让高斯用数值Hessian。
2 不是,高斯默认是自动把GEDII ...

多谢Sob老师!又学习了不少

113

帖子

0

威望

903

eV
积分
1016

Level 4 (黑子)

36#
发表于 Post on 2016-7-11 16:29:08 | 只看该作者 Only view this author
sobereva 发表于 2015-1-28 14:14
你的几何优化和频率计算明显不是同一个级别进行的,出现一堆虚频是必然的

sob老师,我也遇到了类似的情况:
我的输入文件如下:
%chk=cation172.chk
%nprocshared=24
%mem=60GB
# opt=tight freq b3lyp/6-311G(d)

cation172

1 1
C                  2.14795413    1.19505188    0.11503324
C                  1.42130080   -0.06121699   -0.40010118
N                  0.03601540   -0.06227507    0.09171507
C                 -0.65757445    1.13905339   -0.39474465
C                 -0.65959441    1.14242752   -1.93473963
C                 -1.38621227    2.40096214   -2.44436410
C                 -1.38823228    2.40433624   -3.98435908
C                 -2.11485005    3.66287090   -4.49398356
C                  0.03794348   -0.06549587    1.56171027
C                  0.76456142   -1.32403044    2.07133474
C                 -0.65760824   -1.26144080   -0.40000417
C                 -2.10885962   -1.26254927    0.11523189
H                  1.41989737   -0.05887260   -1.47009770
H                  1.92615872   -0.93565334   -0.04601145
H                 -1.66591143    1.13828326   -0.03675591
H                 -0.15269190    2.01191552   -0.03682661
H                 -0.65901166   -1.25909641   -1.47000068
H                 -0.15275032   -2.13587716   -0.04591444
H                 -0.97039352   -0.06626604    1.91969897
H                  0.54282598    0.80736626    1.91962835
H                 -2.61374211   -2.13541141   -0.24268618
H                 -2.10745620   -1.26489365    1.18522841
H                 -2.61371754   -0.38811292   -0.23885785
H                  3.15629112    1.19582207   -0.24295548
H                  1.64309619    2.06948824   -0.23905647
H                  2.14935758    1.19270748    1.18502976
H                  0.76596484   -1.32637484    3.14133125
H                  0.25967891   -2.19689258    1.71341665
H                  1.77289841   -1.32326027    1.71334604
H                  0.34874256    1.14319759   -2.29272838
H                 -1.16447702    0.26956542   -2.29265766
H                 -2.39454924    2.40019209   -2.08637533
H                 -0.88132964    3.27382424   -2.08644610
H                 -0.37989531    2.40510622   -4.34234786
H                 -1.89311496    1.53147417   -4.34227708
C                 -2.11687009    3.66624499   -6.03397854
H                 -3.12318701    3.66210093   -4.13599477
H                 -1.60996736    4.53573296   -4.13606558
C                 -2.84348780    4.92477969   -6.54360301
H                 -1.10853312    3.66701491   -6.39196733
H                 -2.62175282    2.79338294   -6.39189652
H                 -2.84489134    4.92712402   -7.61359952
H                 -3.85182476    4.92400977   -6.18561421
H                 -2.33860506    5.79764173   -6.18568503

得到的log文件中,最后一个的最大受力、均方根位移等是4个YES:
         Item               Value     Threshold  Converged?
Maximum Force            0.000004     0.000015     YES
RMS     Force            0.000000     0.000010     YES
Maximum Displacement     0.000041     0.000060     YES
RMS     Displacement     0.000007     0.000040     YES
而最后显示有一个一个虚频,这个是因为做freq时上面的4项其实并不都是YES的么?只是没有写入到log文件?

cation172.gjf

2.68 KB, 下载次数 Times of downloads: 2

cation172.log

1.03 MB, 下载次数 Times of downloads: 6

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

37#
 楼主 Author| 发表于 Post on 2016-7-11 17:35:31 | 只看该作者 Only view this author
whtu 发表于 2016-7-11 16:29
sob老师,我也遇到了类似的情况:
我的输入文件如下:
%chk=cation172.chk


帖子里已经对原理分析讨论得很详细了,再仔细看。
关键是Hessian不同所导致的,优化时候是近似的,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

121

帖子

1

威望

519

eV
积分
660

Level 4 (黑子)

38#
发表于 Post on 2016-7-11 19:19:58 | 只看该作者 Only view this author
好贴,感谢!
233

3

帖子

0

威望

81

eV
积分
84

Level 2 能力者

39#
发表于 Post on 2017-7-14 09:56:55 | 只看该作者 Only view this author
老师好。结构优化不收敛后采用了帖子里的方法 SCF(novaracc,noincfock) int=acc2e=12  opt freq ,出现opt收敛freq后两项不收敛还差距较大的情况。http://gaussian.com/faq3/建议读取chk文件再优化,这个输入文件该如何写呢?本身计算花费了四天,如果calcall的话耗时会更多吧?感谢。

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

40#
 楼主 Author| 发表于 Post on 2017-7-14 10:07:34 | 只看该作者 Only view this author
小丫姐姐 发表于 2017-7-14 09:56
老师好。结构优化不收敛后采用了帖子里的方法 SCF(novaracc,noincfock) int=acc2e=12  opt freq ,出现opt ...

SCF(novaracc,noincfock) int=acc2e=12这是处理SCF不收敛的,和这里说的帮助结构收敛完全是两码事:
量子化学计算中帮助几何优化收敛的常用方法
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

3

帖子

0

威望

81

eV
积分
84

Level 2 能力者

41#
发表于 Post on 2017-7-14 10:27:03 | 只看该作者 Only view this author
sobereva 发表于 2017-7-14 10:07
SCF(novaracc,noincfock) int=acc2e=12这是处理SCF不收敛的,和这里说的帮助结构收敛完全是两码事:
量 ...

嗯嗯,理解错误。那么我现在的freq不收敛用calcall的话可能耗时太久,我尝试一下用opt=tight freq,int=ultrafine。另外,如果按照http://gaussian.com/faq3/说的读取chk文件的话输入文件该怎么写呢?菜鸟一枚,想着多试几种方法,看哪一种比较合适。

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

42#
 楼主 Author| 发表于 Post on 2017-7-14 10:41:50 | 只看该作者 Only view this author
小丫姐姐 发表于 2017-7-14 10:27
嗯嗯,理解错误。那么我现在的freq不收敛用calcall的话可能耗时太久,我尝试一下用opt=tight freq,int=u ...


甭管高斯页面上的,特别是对于没有判断能力的初学者来说
光从chk里读取结构什么的根本不可能解决问题
北京科音自然科学研究中心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

3

帖子

0

威望

81

eV
积分
84

Level 2 能力者

43#
发表于 Post on 2017-7-14 10:44:22 | 只看该作者 Only view this author
sobereva 发表于 2017-7-14 10:41
甭管高斯页面上的,特别是对于没有判断能力的初学者来说
光从chk里读取结构什么的根本不可能解决问题

好的,那我试试opt=tight freq,int=ultrafine.谢谢老师。

47

帖子

0

威望

131

eV
积分
178

Level 3 能力者

44#
发表于 Post on 2017-9-12 18:24:16 | 只看该作者 Only view this author
sob老师,请问优化分子结构的时候出现了一个虚频,请问是提取出最后一步加opt=tight呢,还是在刚开始的输入文件中加opt=tight 呢?

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

45#
 楼主 Author| 发表于 Post on 2017-9-12 18:51:27 | 只看该作者 Only view this author
曼竹w 发表于 2017-9-12 18:24
sob老师,请问优化分子结构的时候出现了一个虚频,请问是提取出最后一步加opt=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

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 12:13 , Processed in 0.460218 second(s), 28 queries , Gzip On.

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