计算化学公社

标题: 结构优化出现虚频，采用calcall消除虚频后，结构不再是对称的该如何解决 [打印本页]

作者
Author: wangwangshuidi 时间: 2018-12-26 09:28
标题: 结构优化出现虚频，采用calcall消除虚频后，结构不再是对称的该如何解决
如题目描述
当我用关键词“opt=tight freq b3lyp/6-31g(d) int=ultrafine”时，计算得到的碳结构是对称的（图1），但是一直有一个很大的虚频，我用了网上的一些消虚频的方法以后发现虚频没有了，但是优化出来的结构又不是对称的了。最后，万不得已，我用关键词“ opt=(tight,calcall) freq b3lyp/6-31g(d) int=ultrafine”时，直接就得到了没有虚频的结构，但是最终结构也不是对称的（图2）。但是我看了好多相关文献，发现他们最后优化的结果都是对称的。
所以想问问各位大神遇没遇到过这种问题，对此有没有什么建议？

作者
Author: wangwangshuidi 时间: 2018-12-26 09:31
上图中的原子编号不一样是由于我画图画了两次，原子添加顺序不一样导致的，当时初始的键长键角都是一样的，而且是平面结构，二面角都为0°或者180°。

作者
Author: loovfnd 时间: 2018-12-26 14:59
俺也遇到过类似的问题，加上calcall算的很慢，还在摸索呢

作者
Author: loovfnd 时间: 2018-12-26 15:08
考虑换M062x试试，不知道效果如何。

作者
Author: wangwangshuidi 时间: 2018-12-26 16:02

loovfnd 发表于 2018-12-26 15:08
考虑换M062x试试，不知道效果如何。

我参考的文献里面没有人用这个泛函的，怕算出来也不好进行对比验证。

作者
Author: sobereva 时间: 2018-12-27 07:11
M06-2X算这种体系毫无问题。有点理论计算经验和常识的审稿人不可能会质疑这点

不对称就不对称，本来也未必是对称的
B3LYP算这种体系毫无问题，如果怀疑合理性，用6-311G*再试。
另外对最终的波函数做波函数稳定性测试，确认没有收敛到非基态波函数上

作者
Author: wangwangshuidi 时间: 2018-12-27 08:48

sobereva 发表于 2018-12-27 07:11
M06-2X算这种体系毫无问题。有点理论计算经验和常识的审稿人不可能会质疑这点

不对称就不对称，本来也未 ...

好的，谢谢社长。我按您说的方法再试一下看看。

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