计算化学公社

标题: Orca5在B3LYP/G D3 def2-TZVP计算级别下做结构优化收敛不到最优结构的疑问 [打印本页]

作者
Author: we000000 时间: 2023-10-11 16:42
标题: Orca5在B3LYP/G D3 def2-TZVP计算级别下做结构优化收敛不到最优结构的疑问
本帖最后由 we000000 于 2023-10-11 16:43 编辑

我最近尝试在Orca5.0.3上用B3LYP/G D3 def2-TZVP这个计算级别对有机分子1,1,2,2-四氟乙基-2,2,3,3-四氟丙醚（C5H4F8O）做结构优化。这个分子在PubChem数据库中查询到的、以及用VASP跑AIMD稳定后的结构如下所示：红色是O，褐色是C，白色是H，蓝灰色是F

(, 下载次数 Times of downloads: 8)
另外，我用Orca在b97-3c下也运行了一个3000fs的AIMD计算任务，不知道为什么1400fs左右就自动正常中断了，最后的结构大致也和上面的接近。

奇怪的是，因为一开始我没有查数据库中的结构也没有跑AIMD，不确定这个分子两端的-CF2H是怎样旋转着的状态，所以构建了两种-CF2H基团状态的初始结构（都是锯齿状的直线结构）。随后用Orca在B3LYP/G D3 def2-TZVP这个计算级别下分别做了结构优化。但结构优化完毕之后仍然是锯齿状直线型分子，两种不同初始结构在结构优化后得到的结果分别如下所示：
(, 下载次数 Times of downloads: 7) (, 下载次数 Times of downloads: 7)
这和VASP以及ORCA用AIMD跑出来的结构明显不一样，起码不应该是直线型的。

根据这个帖子http://bbs.keinsci.com/thread-16705-2-1.html中楼主的说法，他在关闭RICOSX后结构优化就正常了，我也尝试在inp文件中加了!NoCOSX ，结构优化的结果依然没有明显区别。使用DEFGRID3关键字也尝试过，也是没什么区别。

想请教一下大家是因为我泛函+基组选择不合适的原因吗还是什么原因

作者
Author: we000000 时间: 2023-10-11 16:50
使用上述B3LYP/G计算出的结果作为初始结构，在B97-3C下做结构优化我也尝试过，结果还是原样子。
当然，用正确的初始结构作为输入，然后在上述B3LYP/G D3 def2-TZVP下做结构优化，也不会变成错误的直线型，可以输出正确的结构。就是不能理解为什么初始结构不对的时候做结构优化无法找到最稳定的结构

作者
Author: pal 时间: 2023-10-11 16:56

we000000 发表于 2023-10-11 16:50
使用上述B3LYP/G计算出的结果作为初始结构，在B97-3C下做结构优化我也尝试过，结果还是原样子。
当然，用 ...

结构优化只会根据梯度找最接近的稳定结构，并不保证是最稳定结构没啥毛病

作者
Author: wzkchem5 时间: 2023-10-11 17:04
仔细阅读http://bbs.keinsci.com/thread-577-1-1.html，了解结构优化和构象搜索的区别
如果想从一个错误的构象出发，不经人为调节而自动得到正确的构象，一般必须做构象搜索，只做结构优化是不够的

作者
Author: we000000 时间: 2023-10-11 17:13

pal 发表于 2023-10-11 16:56
结构优化只会根据梯度找最接近的稳定结构，并不保证是最稳定结构没啥毛病

如果不十分确定初始结构的情况似乎只能按照楼下老师说的做构像搜索了，我去了解下谢谢

作者
Author: we000000 时间: 2023-10-11 17:14

wzkchem5 发表于 2023-10-11 17:04
仔细阅读http://bbs.keinsci.com/thread-577-1-1.html，了解结构优化和构象搜索的区别
如果想从一个错误的 ...

我去了解一下，感谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)