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

li447fan · 发表于 Post on 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, 下载次数 Times of downloads: 82)

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

447951397 · 发表于 Post on 2021-3-5 17:41:29

本帖最后由 447951397 于 2021-3-5 17:42 编辑

没有加opt优化到稳定结构，建议取每次优化后的结构多opt几次。且Gaussian的b3lyp对应于ORCA的b3lyp/G

li447fan · 发表于 Post on 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的算

li447fan · 发表于 Post on 2021-3-5 17:52:27

或者ORCA里有像gaussian一样，可以从.gbw文件中readfc精确hessian再优化的方法？

sobereva · 发表于 Post on 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再试

量化小菜鸡 · 发表于 Post on 2021-3-5 19:12:09

有高斯用高斯，我也遇到过这种。G16算出来啥事没有，ORCA 4.2 关掉RI升高格点还是个虚频消不掉。我不知道社长Gaussian和ORCA联用那种办法对这种管不管用。

greatqi · 发表于 Post on 2021-3-5 19:43:20

有微小的虚频，可以忽略不计的。

量化小菜鸡 · 发表于 Post on 2021-3-5 19:49:43

greatqi 发表于 2021-3-5 19:43
有微小的虚频，可以忽略不计的。

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

喵星大佬 · 发表于 Post on 2021-3-5 19:52:37

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

显然不管用

wzkchem5 · 发表于 Post on 2021-3-5 21:16:00

greatqi 发表于 2021-3-5 19:43
有微小的虚频，可以忽略不计的。

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

greatqi · 发表于 Post on 2021-3-6 08:31:45

那你试试看，对称性是否关联，去掉对称性。或者按照wzkchem5解释的那样，稍微调整一下分子的坐标，看能否管用。

西村W · 发表于 Post on 2021-4-13 10:27:41

请问你后来用高斯做出来的还有虚频吗？我最近用Gaussian做的有机物也一直有一个小虚频，不知道怎么解决了。

li447fan · 发表于 Post on 2021-4-20 15:54:09

西村W 发表于 2021-4-13 10:27
请问你后来用高斯做出来的还有虚频吗？我最近用Gaussian做的有机物也一直有一个小虚频，不知道怎么解决了。

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

shenchaoren · 发表于 Post on 2021-4-28 14:12:14

我也在使用ORCA优化结构的时候遇到了类似问题，Sob的OFakeG转换以后，在GaussView上显示有小虚频（红线划出），但是用orca的out文件用Avogardo打开读取freq的结果，没有显出虚频（红圈划出）。

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

li447fan · 发表于 Post on 2021-4-28 15:57:14

shenchaoren 发表于 2021-4-28 14:12
我也在使用ORCA优化结构的时候遇到了类似问题，Sob的OFakeG转换以后，在GaussView上显示有小虚频（红线划出 ...

这个估计要@ 一下社长了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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