计算化学公社

标题: 苯硅烷的虚频如何消除？ [打印本页]

作者
Author: 得浮一大白 时间: 2019-6-20 16:52
标题: 苯硅烷的虚频如何消除？
各位朋友们，你们好，想请教一个关乎消除虚频的问题。本人利用b3lyp/def2svp em=gd3bj来优化苯硅烷的构型，但是无论如何都会出现虚频。参照sob老师之前的那篇文章分享，也曾尝试过添加"tight"、“int=ultrafine”等帮助收敛的关键词，也试过在虚频振动方向上调节，甚至加上calcall来计算，都没有把虚频给除掉。想问各位有什么好建议吗？计算文件放在了附件中，请各位不吝赐教，谢谢！

作者
Author: h840473807 时间: 2019-6-20 19:11
同时添加opt=tight int=ultrafine 没有虚频。我用的是G09 E01
G09 D01 使用D3(bj)可能会有莫名其妙的虚频，社长博文里也提到过

作者
Author: sobereva 时间: 2019-6-21 00:08
此文明确说了
DFT-D色散校正的使用
http://sobereva.com/210

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

作者
Author: 得浮一大白 时间: 2019-6-21 09:07

h840473807 发表于 2019-6-20 19:11
同时添加opt=tight int=ultrafine 没有虚频。我用的是G09 E01
G09 D01 使用D3(bj)可能会有莫名其妙的虚频 ...

原来如此，感谢！

作者
Author: 得浮一大白 时间: 2019-6-21 09:14

sobereva 发表于 2019-6-21 00:08
此文明确说了
DFT-D色散校正的使用
http://sobereva.com/210

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

作者
Author: sobereva 时间: 2019-6-21 10:49

得浮一大白发表于 2019-6-21 09:14
也就是说如果使用零阻尼形式的D3就可能不出现虚频，但这两种结果是否会有很大的不同吗？

对结构影响很小

作者
Author: granvia 时间: 2019-6-28 19:39
顺便扯一篇老JACS文章，发现MP2, MP3, CISD和CCSD等post HF方法在某些常用基组（主要是Pople基组）下计算苯分子，会得到诡异的虚频。Popular Theoretical Methods Predict Benzene and Arenes To Be Nonplanar （J. Am. Chem. Soc. 128, 9342）

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