计算化学公社

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

[ORCA] 无论怎样优化的有机体系始终有微小虚频

[复制链接 Copy URL]

98

帖子

0

威望

823

eV
积分
921

Level 4 (黑子)

分别利用Orca在b3lyp-d3bj/def2-SVP,b3lyp-d3bj/def-TZVP下利用高格点精度(如Grid6,GridX6)以及高收敛标准(tightopt和verytightopt)优化了分子结构,并做freq,但始终会有一处有微小虚频,不知道该怎么处理了,希望有人能帮助我一下。 freq.zip (4.02 MB, 下载次数 Times of downloads: 82)



PS:用Gaussian在b3lyp-d3bj/6-311G*下的计算还在算,不知道结果。

312

帖子

0

威望

2274

eV
积分
2586

Level 5 (御坂)

19#
发表于 Post on 2021-11-6 09:51:07 | 只看该作者 Only view this author
刚才用Orca 5.0.1 优化邻甲基乙苯也出现了虚频,加上TightOPT也不管用。后来调了一下乙基末端的甲基位置,重新优化,才消除了虚频。

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

18#
发表于 Post on 2021-5-3 02:47:38 | 只看该作者 Only view this author
shenchaoren 发表于 2021-4-28 14:12
我也在使用ORCA优化结构的时候遇到了类似问题,Sob的OFakeG转换以后,在GaussView上显示有小虚频(红线划出 ...

avogadro该显示虚频没显示,是avogadro的问题
北京科音自然科学研究中心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

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

17#
发表于 Post on 2021-4-28 18:50:52 | 只看该作者 Only view this author
shenchaoren 发表于 2021-4-28 14:12
我也在使用ORCA优化结构的时候遇到了类似问题,Sob的OFakeG转换以后,在GaussView上显示有小虚频(红线划出 ...

你在Avogardo里数数看振动模式数量对不对。。。。

334

帖子

1

威望

6159

eV
积分
6513

Level 6 (一方通行)

16#
发表于 Post on 2021-4-28 17:02:53 | 只看该作者 Only view this author
本帖最后由 paramecium86 于 2021-4-28 17:07 编辑
shenchaoren 发表于 2021-4-28 14:12
我也在使用ORCA优化结构的时候遇到了类似问题,Sob的OFakeG转换以后,在GaussView上显示有小虚频(红线划出 ...

我以前也试过,ORCA计算完TS之后 在末尾会输出一个IR 光谱的振动模式。这里面会把虚频的振动模式给去掉。也就是振动模式会少显示一个,而Avogardo读的就是末尾的这个振动输出。数一下个数就能发现。不是3N-6个振动模式。
而原始的频率计算结果(包含虚频的)在频率计算的输出文件靠前的位置,标着 “VIBRATIONAL FREQUENCIES”。ORCA会输出3N个模式,其中0-5号是平动转动忽略掉, 从6号开始是振动频率。

98

帖子

0

威望

823

eV
积分
921

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2021-4-28 15:57:14 | 只看该作者 Only view this author
shenchaoren 发表于 2021-4-28 14:12
我也在使用ORCA优化结构的时候遇到了类似问题,Sob的OFakeG转换以后,在GaussView上显示有小虚频(红线划出 ...

这个估计要@ 一下社长了

19

帖子

2

威望

1672

eV
积分
1731

Level 5 (御坂)

14#
发表于 Post on 2021-4-28 14:12:14 | 只看该作者 Only view this author
我也在使用ORCA优化结构的时候遇到了类似问题,Sob的OFakeG转换以后,在GaussView上显示有小虚频(红线划出),但是用orca的out文件用Avogardo打开读取freq的结果,没有显出虚频(红圈划出)。

ketone-2_orca_fake.out (38.39 KB, 下载次数 Times of downloads: 8)

屏幕截图 2021-04-28 140749.png (124.51 KB, 下载次数 Times of downloads: 49)

屏幕截图 2021-04-28 140749.png

98

帖子

0

威望

823

eV
积分
921

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2021-4-20 15:54:09 | 只看该作者 Only view this author
西村W 发表于 2021-4-13 10:27
请问你后来用高斯做出来的还有虚频吗?我最近用Gaussian做的有机物也一直有一个小虚频,不知道怎么解决了。

就recalc了很久,估计用精确hessian也可以解决,就是会更慢

17

帖子

0

威望

124

eV
积分
141

Level 2 能力者

12#
发表于 Post on 2021-4-13 10:27:41 | 只看该作者 Only view this author
请问你后来用高斯做出来的还有虚频吗?我最近用Gaussian做的有机物也一直有一个小虚频,不知道怎么解决了。

14

帖子

0

威望

1586

eV
积分
1600

Level 5 (御坂)

11#
发表于 Post on 2021-3-6 08:31:45 | 只看该作者 Only view this author
那你试试看,对称性是否关联,去掉对称性。或者按照wzkchem5解释的那样,稍微调整一下分子的坐标,看能否管用。

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

10#
发表于 Post on 2021-3-5 21:16:00 | 只看该作者 Only view this author
greatqi 发表于 2021-3-5 19:43
有微小的虚频,可以忽略不计的。

这个不一定。首先得证明这个虚频确实是假虚频,方法是把分子往虚频的方向扰动一下,看梯度是指向平衡结构还是远离平衡结构。其次只有计算零点能、焓之类的可以忽略小虚频,对于自由能,即使已知虚频是假虚频也不能简单忽略,因为把一个很低的实频算成虚频对自由能影响非常大。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

9#
发表于 Post on 2021-3-5 19:52:37 | 只看该作者 Only view this author
量化小菜鸡 发表于 2021-3-5 19:12
有高斯用高斯,我也遇到过这种。G16算出来啥事没有,ORCA 4.2 关掉RI升高格点还是个虚频消不掉。我不知道社 ...

显然不管用

199

帖子

2

威望

3559

eV
积分
3798

Level 5 (御坂)

8#
发表于 Post on 2021-3-5 19:49:43 | 只看该作者 Only view this author
greatqi 发表于 2021-3-5 19:43
有微小的虚频,可以忽略不计的。

这种发文章的时候怎么说?我们老板要求必须全去掉才行。

14

帖子

0

威望

1586

eV
积分
1600

Level 5 (御坂)

7#
发表于 Post on 2021-3-5 19:43:20 | 只看该作者 Only view this author
有微小的虚频,可以忽略不计的。

199

帖子

2

威望

3559

eV
积分
3798

Level 5 (御坂)

6#
发表于 Post on 2021-3-5 19:12:09 | 只看该作者 Only view this author
有高斯用高斯,我也遇到过这种。G16算出来啥事没有,ORCA 4.2 关掉RI升高格点还是个虚频消不掉。我不知道社长Gaussian和ORCA联用那种办法对这种管不管用。

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

GMT+8, 2026-2-26 01:20 , Processed in 0.195073 second(s), 24 queries , Gzip On.

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