计算化学公社

标题: 关于含碘的赝势计算出现虚频 [打印本页]

作者
Author: suosuosky 时间: 2019-3-31 14:16
标题: 关于含碘的赝势计算出现虚频
请问各位老师，在计算碘取代的乙炔分子时，输入文件：
#p opt freq mp2/genecp
Title Card Required
0 1
C                -1.85915491 -0.15492958 0.00000000
C                -0.65795491 -0.15492958 0.00000000
H                0.41204509 -0.15492958 0.00000000
I                -3.95915491 -0.15492958 0.00000000
I 0
SDD
****
H C 0
6-311++G**
****
I 0
SDD

但是结果分子出现两个虚频，请问是什么原因？应如何修改？谢谢

作者
Author: klklklzzd 时间: 2019-3-31 14:30
http://sobereva.com/278

作者
Author: suosuosky 时间: 2019-3-31 16:45
本帖最后由 suosuosky 于 2019-3-31 17:01 编辑

我已经微调了构型，还是有两个虚频，分子又变回直线。
不知道这个问题应该如何解决？

作者
Author: suosuosky 时间: 2019-3-31 17:04
当调整构型较大后，虚频消除了，但是分子明显偏离直线，不符合常理。附件是没有虚频的输出文件。请问应如何理解？

作者
Author: sobereva 时间: 2019-4-1 02:56
计算级别用得太烂，MP2又贵精度又烂，仔细看此文末尾
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

给H、C加弥散完全是在浪费时间。而且当弥散函数多而高角动量极化函数不够的时候，MP2甚至能把苯优化成弯的，这早就有过报道：J. AM. CHEM. SOC. 2006, 128, 9342-9343。你当前的情况再次印证了这一点。该用什么基组、弥散函数什么时候才该加仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈弥散函数和“月份”基组
http://sobereva.com/119

哪怕在你原先的级别上，把弥散函数给去了，都能得到没虚频的直线结构。

所以说你当前的计算级别属于超级典型的反面教材：乱用MP2、乱加弥散、基组和理论方法不搭

这种体系直接用个B3LYP/def2TZVP opt freq就完事了，直线，没虚频，结果足够精确
(, 下载次数 Times of downloads: 63)

作者
Author: suosuosky 时间: 2019-4-1 09:54

sobereva 发表于 2019-4-1 02:56
计算级别用得太烂，MP2又贵精度又烂，仔细看此文末尾
简谈量子化学计算中DFT泛函的选择
http://sobereva. ...

谢谢sob老师，确实是反面典型 ...

我选择这个方法和基组是因为想比较单体和复合物中C-H键长的变化，同时计算结合能。因为复合物的构型优化全部使用的都是MP2/6-311++g**，I使用SDD。为了有可比性，所以单体也选择这种方法和基组。

那按照您的方法，我优化的复合物是不是全部要修改计算方法和基组？复合物不同构型中有氢键、卤键，10多个原子，基组是不是应该要加上弥散呢？如果我使用MP2/aug-cc-pVTZ, 碘使用aug-cc-pVTZ-PP，是否可以？

作者
Author: sobereva 时间: 2019-4-1 11:46

suosuosky 发表于 2019-4-1 09:54
谢谢sob老师，确实是反面典型 ...

我选择这个方法和基组是因为想比较单体和复合物中C-H键长的变化，同 ...

复合物就明显不应该用那个级别优化，又贵又差，很容易遭到对计算级别比较讲究的审稿人的吐槽

opt freq时候，方法用B3LYP-D3(BJ)，基组用def-TZVP或者6-311G**，碘用Lanl08(d)，又好又快。弥散不是必需的。
算能量时候再进一步提升级别，诸如用双杂化泛函，结合jun-cc-pVTZ和jun-cc-pVTZ-PP。

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