请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 603|回复: 6

[新手求助] 关于含碘的赝势计算出现虚频

[复制链接]

30

帖子

0

威望

207

eV
积分
237

Level 3 能力者

发表于 2019-3-31 14:16:11 | 显示全部楼层 |阅读模式
请问各位老师,在计算碘取代的乙炔分子时,输入文件:
#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

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

12

帖子

0

威望

386

eV
积分
398

Level 3 能力者

发表于 2019-3-31 14:30:30 | 显示全部楼层

30

帖子

0

威望

207

eV
积分
237

Level 3 能力者

 楼主| 发表于 2019-3-31 16:45:05 | 显示全部楼层
本帖最后由 suosuosky 于 2019-3-31 17:01 编辑

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

C2HI.log

141.15 KB, 下载次数: 4

30

帖子

0

威望

207

eV
积分
237

Level 3 能力者

 楼主| 发表于 2019-3-31 17:04:15 | 显示全部楼层
当调整构型较大后,虚频消除了,但是分子明显偏离直线,不符合常理。 附件是 没有虚频的输出文件。请问应如何理解?

C2HI-no.log

325.48 KB, 下载次数: 4

2万

帖子

25

威望

2万

eV
积分
45935

管理员

公社社长

发表于 2019-4-1 02:56:50 | 显示全部楼层
计算级别用得太烂,MP2又贵精度又烂,仔细看此文末尾
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html

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

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

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

这种体系直接用个B3LYP/def2TZVP opt freq就完事了,直线,没虚频,结果足够精确
1.out (81.67 KB, 下载次数: 8)

评分

参与人数 2eV +4 收起 理由
王二葛 + 3 谢谢
Novice + 1 我很赞同

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

30

帖子

0

威望

207

eV
积分
237

Level 3 能力者

 楼主| 发表于 2019-4-1 09:54:22 | 显示全部楼层
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,是否可以?

2万

帖子

25

威望

2万

eV
积分
45935

管理员

公社社长

发表于 2019-4-1 11:46:53 | 显示全部楼层
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://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2019-9-22 08:24 , Processed in 0.171323 second(s), 27 queries .

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