计算化学公社

标题: CO2阴离子进行角度扫描出现能量突变 [打印本页]

作者
Author: chenweiguang 时间: 2024-2-19 11:33
标题: CO2阴离子进行角度扫描出现能量突变
大家好，我初接触量子化学，操作上感觉跟第一性原理还是有挺大差别。刚接触orca，orca手册都1000多页，相当于是便算边学。
最近算CO2阴离子的能量随OCO夹角的关系，发现出现能量突变的情况，试了很多方法和基组都有这个问题。比如最开始用的这样的输入文件
! B3LYP def2-SVP Opt

%scf
ConvForced true
end

%geom Scan
A 1 0 2 = 180, 100, 81 #A for Angle
end
end

*xyz -1 2
C       0.00000000    0.00000000    0.00000000
O       0.00000000    0.00000000    1.17865800
O       0.00000000    0.00000000    -1.17865800
*
但是得到的结果是这样的

作者
Author: wjc404 时间: 2024-2-19 12:03
本帖最后由 wjc404 于 2024-2-29 16:20 编辑

检查152-180度的优化是不是跑到了不稳定的波函数上，使用关键词%scf stabperform true end

另外这么小的体系用def2-SVP有些抠门，且阴离子体系应该带上s/p的弥散函数。

用CCSD(T)/jun-cc-pVTZ对CO2阴离子进行柔性扫描，结果见下图。
(, 下载次数 Times of downloads: 32)

作者
Author: chenweiguang 时间: 2024-2-19 12:43

wjc404 发表于 2024-2-19 12:03
检查152-180度的优化是不是跑到了不稳定的波函数上，使用关键词%scf stabperform true end

另外这么小的 ...

多谢！已经在做波函数检查了

请问有什么推荐的基组吗？

作者
Author: chenweiguang 时间: 2024-2-19 12:46
我还试了一些其他基组，结果如下

作者
Author: sobereva 时间: 2024-2-19 14:23

chenweiguang 发表于 2024-2-19 12:43
多谢！已经在做波函数检查了

请问有什么推荐的基组吗？

这么点体系，而且还是阴离子，直接用ma-def2-TZVP
仔细看下文了解常识
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈弥散函数和“月份”基组
http://sobereva.com/119

没事甭用BP86，GGA泛函SIE问题很大，算阴离子极其不适合
仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
用wB97X-D3或wB97M-V

作者
Author: chenweiguang 时间: 2024-2-19 15:07

sobereva 发表于 2024-2-19 14:23
这么点体系，而且还是阴离子，直接用ma-def2-TZVP
仔细看下文了解常识
谈谈量子化学中基组的选择

多谢大神，我试试！

作者
Author: wjc404 时间: 2024-2-25 14:26
这个体系，用B-CCD(T)/may-cc-pVQZ做柔性扫描，结果如下：
(, 下载次数 Times of downloads: 33)

作者
Author: chenweiguang 时间: 2024-7-30 11:57

wjc404 发表于 2024-2-25 14:26
这个体系，用B-CCD(T)/may-cc-pVQZ做柔性扫描，结果如下：

好久没来了，刚看到消息。这个实验上是135度是能量最小的位置，与你的结果有点差异。

作者
Author: wzkchem5 时间: 2024-7-30 15:10

chenweiguang 发表于 2024-7-30 04:57
好久没来了，刚看到消息。这个实验上是135度是能量最小的位置，与你的结果有点差异。

仔细阅读http://sobereva.com/569

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