请选择 进入手机版 | 继续访问电脑版
第9届北京科音分子动力学与GROMACS培训班将于4月17~20日于北京举办,请点击此链接查看培训详情,欢迎参加和相互转告!

计算化学公社

 找回密码
 现在注册!
查看: 319|回复: 10

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

[复制链接]

62

帖子

0

威望

513

eV
积分
575

Level 4 (黑子)

发表于 2021-3-5 16:54:54 | 显示全部楼层 |阅读模式
分别利用Orca在b3lyp-d3bj/def2-SVP,b3lyp-d3bj/def-TZVP下利用高格点精度(如Grid6,GridX6)以及高收敛标准(tightopt和verytightopt)优化了分子结构,并做freq,但始终会有一处有微小虚频,不知道该怎么处理了,希望有人能帮助我一下。 freq.zip (4.02 MB, 下载次数: 7)

46

帖子

0

威望

764

eV
积分
810

Level 4 (黑子)

发表于 2021-3-5 17:41:29 | 显示全部楼层
本帖最后由 447951397 于 2021-3-5 17:42 编辑

没有加opt优化到稳定结构,建议取每次优化后的结构多opt几次。且Gaussian的b3lyp对应于ORCA的b3lyp/G
专注于面向个人的量化问题解决,LAMMPS ReaxFF力场。一对一提供服务,本人个人电脑提供计算资源,双路8168,128G内存。

62

帖子

0

威望

513

eV
积分
575

Level 4 (黑子)

 楼主| 发表于 2021-3-5 17:50:00 | 显示全部楼层
447951397 发表于 2021-3-5 17:41
没有加opt优化到稳定结构,建议取每次优化后的结构多opt几次。且Gaussian的b3lyp对应于ORCA的b3lyp/G

已经opt了挺多次了,各种初始结构和收敛限标准。
orca的是只在orca的算,的确没有用gaussian的b3lyp定义形式。gaussian的在一直用gaussian的算

62

帖子

0

威望

513

eV
积分
575

Level 4 (黑子)

 楼主| 发表于 2021-3-5 17:52:27 | 显示全部楼层
或者ORCA里有像gaussian一样,可以从.gbw文件中readfc精确hessian再优化的方法?

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

发表于 2021-3-5 18:03:20 | 显示全部楼层
既然有Gaussian就直接用Gaussian算就完了,本来ORCA的freq速度就慢,为虚频折腾还耽误时间。之前论坛里有说ORCA 4.2.1的RIJCOSX的解析Hessian还有bug(https://orcaforum.kofo.mpg.de/viewtopic.php?f=26&t=5566),非要用的话加上%scf thresh 0 end再试
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

109

帖子

0

威望

1515

eV
积分
1624

Level 5 (御坂)

发表于 2021-3-5 19:12:09 | 显示全部楼层
有高斯用高斯,我也遇到过这种。G16算出来啥事没有,ORCA 4.2 关掉RI升高格点还是个虚频消不掉。我不知道社长Gaussian和ORCA联用那种办法对这种管不管用。

14

帖子

0

威望

610

eV
积分
624

Level 4 (黑子)

发表于 2021-3-5 19:43:20 | 显示全部楼层
有微小的虚频,可以忽略不计的。

109

帖子

0

威望

1515

eV
积分
1624

Level 5 (御坂)

发表于 2021-3-5 19:49:43 | 显示全部楼层
greatqi 发表于 2021-3-5 19:43
有微小的虚频,可以忽略不计的。

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

703

帖子

1

威望

839

eV
积分
1562

Level 5 (御坂)

喵星人

发表于 2021-3-5 19:52:37 | 显示全部楼层
量化小菜鸡 发表于 2021-3-5 19:12
有高斯用高斯,我也遇到过这种。G16算出来啥事没有,ORCA 4.2 关掉RI升高格点还是个虚频消不掉。我不知道社 ...

显然不管用

1424

帖子

0

威望

854

eV
积分
2278

Level 5 (御坂)

发表于 2021-3-5 21:16:00 | 显示全部楼层
greatqi 发表于 2021-3-5 19:43
有微小的虚频,可以忽略不计的。

这个不一定。首先得证明这个虚频确实是假虚频,方法是把分子往虚频的方向扰动一下,看梯度是指向平衡结构还是远离平衡结构。其次只有计算零点能、焓之类的可以忽略小虚频,对于自由能,即使已知虚频是假虚频也不能简单忽略,因为把一个很低的实频算成虚频对自由能影响非常大。

14

帖子

0

威望

610

eV
积分
624

Level 4 (黑子)

发表于 2021-3-6 08:31:45 | 显示全部楼层
那你试试看,对称性是否关联,去掉对称性。或者按照wzkchem5解释的那样,稍微调整一下分子的坐标,看能否管用。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )

GMT+8, 2021-4-11 17:37 , Processed in 0.289233 second(s), 28 queries .

快速回复 返回顶部 返回列表