计算化学公社

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

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

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

61#
 楼主 Author| 发表于 Post on 2020-12-21 00:39:23 | 只看该作者 Only view this author
li447fan 发表于 2020-12-21 00:30
社长,用了gridx6计算后,虚频反而比之前误用grid5更严重了,发现opt结束时不是全yes,是有一个no的,只 ...

分清楚gridx和grid,各司其职
北京科音自然科学研究中心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

98

帖子

0

威望

823

eV
积分
921

Level 4 (黑子)

62#
发表于 Post on 2021-3-29 10:14:16 | 只看该作者 Only view this author
在加入recalc=3之后,结构反而越来越难收敛了,基本只有1个yes或都是no,跑了快一个星期了。请问是继续等呢,还是尝试换个方法。(一开始,B3LYP-D3/6-311G**下opt四个yes,freq两个yes)

601

帖子

1

威望

4776

eV
积分
5397

Level 6 (一方通行)

63#
发表于 Post on 2021-3-29 16:52:36 | 只看该作者 Only view this author
li447fan 发表于 2021-3-29 10:14
在加入recalc=3之后,结构反而越来越难收敛了,基本只有1个yes或都是no,跑了快一个星期了。请问是继续等呢 ...

优化的啥,能贴一下输入输出文件不
现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

31

帖子

0

威望

362

eV
积分
393

Level 3 能力者

64#
发表于 Post on 2021-4-9 14:00:55 | 只看该作者 Only view this author
想请问老师,计算频率的时候虽然出现了两个NO,但无虚频,这样优化后的结构是否能用?谢谢您
         Item               Value     Threshold  Converged?
Maximum Force            0.000021     0.000450     YES
RMS     Force            0.000003     0.000300     YES
Maximum Displacement     0.008541     0.001800     NO
RMS     Displacement     0.002329     0.001200     NO

1657

帖子

5

威望

4554

eV
积分
6311

Level 6 (一方通行)

喵星人

65#
发表于 Post on 2021-4-9 19:47:42 | 只看该作者 Only view this author
HJC 发表于 2021-4-9 14:00
想请问老师,计算频率的时候虽然出现了两个NO,但无虚频,这样优化后的结构是否能用?谢谢您
         Ite ...

说明势能面比较平滑,力已经充分收敛了但是位移还没有,影响不大,如果介意也可以用这次振动分析得到的heseeian再优化一下,基本就OK了

10

帖子

0

威望

375

eV
积分
385

Level 3 能力者

66#
发表于 Post on 2022-10-19 01:25:29 | 只看该作者 Only view this author
请问老师,我在做激发态优化时,opt的4个都是yes,但是用opt后的结构做频率计算,4个都是no,均是在同一计算级别下做的,这是怎么回事?

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

67#
 楼主 Author| 发表于 Post on 2022-10-19 01:26:46 | 只看该作者 Only view this author
jiuyangli 发表于 2022-10-19 01:25
请问老师,我在做激发态优化时,opt的4个都是yes,但是用opt后的结构做频率计算,4个都是no,均是在同一计 ...

freq的时候guess=read读取opt最后的波函数再试。若还有这个问题,把opt和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

10

帖子

0

威望

375

eV
积分
385

Level 3 能力者

68#
发表于 Post on 2022-10-19 09:00:48 | 只看该作者 Only view this author
sobereva 发表于 2022-10-19 01:26
freq的时候guess=read读取opt最后的波函数再试。若还有这个问题,把opt和freq的输入、输出文件都压缩后上 ...

老师,这是guess=read重新计算的,还是这样。 DOMB-T1-OptFreq.zip (561.65 KB, 下载次数 Times of downloads: 2)


6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

69#
 楼主 Author| 发表于 Post on 2022-10-19 18:58:54 | 只看该作者 Only view this author
jiuyangli 发表于 2022-10-19 09:00
老师,这是guess=read重新计算的,还是这样。

opt的gjf也上传,否则都不知道你用modredundant读的是什么
北京科音自然科学研究中心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

10

帖子

0

威望

375

eV
积分
385

Level 3 能力者

70#
发表于 Post on 2022-10-22 09:34:22 | 只看该作者 Only view this author
sobereva 发表于 2022-10-19 18:58
opt的gjf也上传,否则都不知道你用modredundant读的是什么

老师,这是重新上传后的。

DOMB-T1-OptFreq.zip

563.06 KB, 阅读权限: 100, 下载次数 Times of downloads: 1

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

71#
 楼主 Author| 发表于 Post on 2022-10-23 04:17:36 | 只看该作者 Only view this author
jiuyangli 发表于 2022-10-22 09:34
老师,这是重新上传后的。

写modredundant是完全多余的

使用IOp自定义泛函的时候不要用opt freq,因为IOp没法传递到后续子任务,导致freq部分和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

81

帖子

0

威望

653

eV
积分
734

Level 4 (黑子)

72#
发表于 Post on 2024-3-15 09:46:47 | 只看该作者 Only view this author
老师您好,我在做一个由CHONP元素组成的体系,使用orca501优化结构,关键词是! b3lyp/G cc-pVTZ d3 opt freq,计算的结果总是有虚频,使用tihgtopt,最后的结果仍然有虚频,利用“Recalc_Hess 3”则无法收敛;使用高斯用相同的关键词最后得到的结构没有虚频,拿该结构再用orca仍然会有虚频;使用cc-pvdz的话则可以收敛且没有虚频,但是使用该结构再去优化,仍然会有虚频;且上述虚频都是4个,大小和振动模式十分相似,想问一下老师这种情况有什么好的解决办法吗。

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

73#
 楼主 Author| 发表于 Post on 2024-3-15 17:35:00 | 只看该作者 Only view this author
Sui1999 发表于 2024-3-15 09:46
老师您好,我在做一个由CHONP元素组成的体系,使用orca501优化结构,关键词是! b3lyp/G cc-pVTZ d3 opt fre ...

提高积分格点精度再试

我自己基本不用ORCA做opt freq任务

cc-pVTZ给B3LYP用不是好主意,下文说了
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html

北京科音自然科学研究中心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

81

帖子

0

威望

653

eV
积分
734

Level 4 (黑子)

74#
发表于 Post on 2024-3-15 18:42:40 | 只看该作者 Only view this author
sobereva 发表于 2024-3-15 17:35
提高积分格点精度再试

我自己基本不用ORCA做opt freq任务

感谢老师的回复,已经试过关键词“defgrid3”,已经是orca最高的格点精度了,还是会产生虚频,且虚频和之前的还是很类似。

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

75#
 楼主 Author| 发表于 Post on 2024-3-15 20:06:36 | 只看该作者 Only view this author
Sui1999 发表于 2024-3-15 18:42
感谢老师的回复,已经试过关键词“defgrid3”,已经是orca最高的格点精度了,还是会产生虚频,且虚频和之 ...

懒得折腾就用Gaussian吧
北京科音自然科学研究中心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 17:25 , Processed in 0.222353 second(s), 22 queries , Gzip On.

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