计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8069|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 苯硅烷的虚频如何消除?

[复制链接 Copy URL]

3

帖子

0

威望

23

eV
积分
26

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位朋友们,你们好,想请教一个关乎消除虚频的问题。本人利用b3lyp/def2svp em=gd3bj来优化苯硅烷的构型,但是无论如何都会出现虚频。参照sob老师之前的那篇文章分享,也曾尝试过添加"tight"、“int=ultrafine”等帮助收敛的关键词,也试过在虚频振动方向上调节,甚至加上calcall来计算,都没有把虚频给除掉。想问各位有什么好建议吗?计算文件放在了附件中,请各位不吝赐教,谢谢!

PhSiH3.zip

42.46 KB, 下载次数 Times of downloads: 6

输入和输出文件

124

帖子

0

威望

2484

eV
积分
2608

Level 5 (御坂)

2#
发表于 Post on 2019-6-20 19:11:42 | 只看该作者 Only view this author
同时添加opt=tight int=ultrafine 没有虚频。我用的是G09 E01
G09 D01 使用D3(bj)可能会有莫名其妙的虚频,社长博文里也提到过

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

3#
发表于 Post on 2019-6-21 00:08:39 | 只看该作者 Only view this author
此文明确说了
DFT-D色散校正的使用
http://sobereva.com/210

注意G09 D.01版D3(BJ)计算频率时色散校正的贡献是不完全准确的,因为忽略掉了依赖于r^-8的那项。对于力常数较大(较硬)的振动模式这个问题的影响可忽略,但对于力常数较小(较软)的振动模式就会造成一定误差,往往会导致虚频的出现(即便优化后已经准确处在极小点位置了)。但是零阻尼没有这个问题,所以G09 D.01的用户如果要算频率的话,还是在零阻尼下优化为宜。不过D3(BJ)的这个问题在G09 E.01版已经得到了修正。


北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

3

帖子

0

威望

23

eV
积分
26

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2019-6-21 09:07:05 | 只看该作者 Only view this author
h840473807 发表于 2019-6-20 19:11
同时添加opt=tight int=ultrafine 没有虚频。我用的是G09 E01
G09 D01 使用D3(bj)可能会有莫名其妙的虚频 ...

原来如此,感谢!

3

帖子

0

威望

23

eV
积分
26

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2019-6-21 09:14:02 | 只看该作者 Only view this author
sobereva 发表于 2019-6-21 00:08
此文明确说了
DFT-D色散校正的使用
http://sobereva.com/210

也就是说如果使用零阻尼形式的D3就可能不出现虚频,但这两种结果是否会有很大的不同吗?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

6#
发表于 Post on 2019-6-21 10:49:59 | 只看该作者 Only view this author
得浮一大白 发表于 2019-6-21 09:14
也就是说如果使用零阻尼形式的D3就可能不出现虚频,但这两种结果是否会有很大的不同吗?

对结构影响很小
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

7#
发表于 Post on 2019-6-28 19:39:51 | 只看该作者 Only view this author
顺便扯一篇老JACS文章,发现MP2, MP3, CISD和CCSD等post HF方法在某些常用基组(主要是Pople基组)下计算苯分子,会得到诡异的虚频。Popular Theoretical Methods Predict Benzene and Arenes To Be Nonplanar (J. Am. Chem. Soc. 128, 9342)

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-25 03:31 , Processed in 0.210329 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list