计算化学公社

标题: 优化分子为啥总有一个甲基转动的虚频 [打印本页]

作者
Author: 剑与玫瑰 时间: 2023-12-16 12:10
标题: 优化分子为啥总有一个甲基转动的虚频
我优化一个分子时后算出来总是出现一个甲基转动的虚频，不知道怎么回事，大家看一下这是我的输入和输出文件 (, 下载次数 Times of downloads: 8)
(, 下载次数 Times of downloads: 5)

作者
Author: GoldenBaby 时间: 2023-12-16 13:50
我也发现过这个问题，可能是Gaussian的一个bug吧。我倒是能提供几个解决方案，不保证一定能成功。
1.改变几何优化步长，可以搭配使用更严格的收敛限制，但是改变收敛限这个方法不一定能起作用，这个得看运气
2.尝试使用opt=calcall，我试过，只计算一次calcfc或者使用recalc都不管用，最有效的还是calcall，但是也不能保证每次都有效，所以这个并不作为第一推荐方案
3.改用orca，当时我的那个结构用Gaussian死活都消不掉虚频，但是orca轻松解决，所以我怀疑是Gaussian的一个bug

作者
Author: Loading0760 时间: 2023-12-16 17:47
用g09，加了int=ultrafine，收敛了，无虚频。

作者
Author: sobereva 时间: 2023-12-16 18:16
照常解决，没有特殊性：
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

opt freq不很建议用SMD，否则容易出虚假的虚频，上文和下文都说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

作者
Author: sobereva 时间: 2023-12-16 18:16

GoldenBaby 发表于 2023-12-16 13:50
我也发现过这个问题，可能是Gaussian的一个bug吧。我倒是能提供几个解决方案，不保证一定能成功。
1.改变 ...

不是bug。按照下文处理，Gaussian不存在消不掉的虚频，我遇见过的众多无论多讨厌的虚频最终用合适的做法总能合理地解决掉。唯一能说得上是bug的也就是G09 D.01的DFT-D3(BJ)的解析Hessian不准确造成的虚频
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

作者
Author: 剑与玫瑰 时间: 2023-12-17 16:37

GoldenBaby 发表于 2023-12-16 13:50
我也发现过这个问题，可能是Gaussian的一个bug吧。我倒是能提供几个解决方案，不保证一定能成功。
1.改变 ...

感谢回答

作者
Author: 剑与玫瑰 时间: 2023-12-17 16:37

sobereva 发表于 2023-12-16 18:16
照常解决，没有特殊性：
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva. ...

谢谢老师

作者
Author: 剑与玫瑰 时间: 2023-12-18 09:02

Loading0760 发表于 2023-12-16 17:47
用g09，加了int=ultrafine，收敛了，无虚频。

感谢，我加上去用G16试一下看看是不是也可行

作者
Author: Loading0760 时间: 2023-12-18 14:02

剑与玫瑰发表于 2023-12-18 09:02
感谢，我加上去用G16试一下看看是不是也可行

Gaussian 16里面int=ultrafine应该是默认的

作者
Author: 剑与玫瑰 时间: 2023-12-18 14:49

Loading0760 发表于 2023-12-18 14:02
Gaussian 16里面int=ultrafine应该是默认的

我试过了，还是有甲基转动的虚频

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