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

计算化学公社

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

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

[复制链接]

55

帖子

0

威望

405

eV
积分
460

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

威望

416

eV
积分
428

Level 3 能力者

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

55

帖子

0

威望

405

eV
积分
460

Level 3 能力者

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

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

C2HI.log

141.15 KB, 下载次数: 4

55

帖子

0

威望

405

eV
积分
460

Level 3 能力者

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

C2HI-no.log

325.48 KB, 下载次数: 8

2万

帖子

25

威望

2万

eV
积分
54438

管理员

公社社长

发表于 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, 下载次数: 13)

评分

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

查看全部评分

第12届北京科音初级量子化学培训班将于六月举办!是零基础和初学者速成量子化学计算的非常难得和宝贵的机会,欢迎参加!详情和报名方式见http://www.keinsci.com/workshop/12KEQC
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计约6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

Money and papers are rubbish, get a real life!

55

帖子

0

威望

405

eV
积分
460

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
积分
54438

管理员

公社社长

发表于 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。

第12届北京科音初级量子化学培训班将于六月举办!是零基础和初学者速成量子化学计算的非常难得和宝贵的机会,欢迎参加!详情和报名方式见http://www.keinsci.com/workshop/12KEQC
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计约6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2020-5-31 12:06 , Processed in 0.154781 second(s), 27 queries .

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