计算化学公社

标题: 结构优化收敛后出现NO和虚频 [打印本页]

作者
Author: 吃饭睡觉打豆豆 时间: 2023-4-2 15:40
标题: 结构优化收敛后出现NO和虚频
大家好，我最近做计算时碰到结构优化后，计算频率时出现了NO和虚频，我也看过sob老师相关的博文，通过我做类似结构优化时的经验来看，分子在构象上应该是有一个较大的二面角，我在调整初始结构的时候也尝试扭转一个二面角，但是每次优化到最后都会优化到一个近乎平面的结构，这应该是出现虚频的原因，是不是因为步长太大，导致走过了，请问这种情况可以通过减小步长上限的方法解决吗？

作者
Author: wzkchem5 时间: 2023-4-2 16:09
这个分子本来应该是平面吗？扭转之前是严格平面、接近平面还是非平面？二面角扭转了多少？
把这些细节都交代清楚，否则无法回答

作者
Author: sobereva 时间: 2023-4-2 16:15
上传压缩后的输出文件。当前描述过于抽象

下文里已经罗列了所有可能有效的解决方法，没有什么可补充的。没有更具体细节没法针对性地回答
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

作者
Author: 吃饭睡觉打豆豆 时间: 2023-4-2 16:24

wzkchem5 发表于 2023-4-2 16:09
这个分子本来应该是平面吗？扭转之前是严格平面、接近平面还是非平面？二面角扭转了多少？
把这些细节都交 ...

由计算类似结构的分子来看，这个分子应该有一个较大的二面角。计算类似结构的分子也出现过虚频的情况，结构优化得到一个近乎平面，二面角位0.007°的结构，我通过扭转二面角最终得到一个二面角约为20°的结构，成功消掉了虚频。目前正在计算的分子结构优化后得到一个二面角为0.02°的近平面结构，出现虚频，我将初始结构扭转角定位15°，想通过同样的方法消掉虚频，但是每次都会优化为近平面结构，出现虚频。

作者
Author: wzkchem5 时间: 2023-4-2 16:31

吃饭睡觉打豆豆发表于 2023-4-2 09:24
由计算类似结构的分子来看，这个分子应该有一个较大的二面角。计算类似结构的分子也出现过虚频的情况，结 ...

试过计算初始Hessian了吗？

作者
Author: 吃饭睡觉打豆豆 时间: 2023-4-2 16:33

sobereva 发表于 2023-4-2 16:15
上传压缩后的输出文件。当前描述过于抽象

下文里已经罗列了所有可能有效的解决方法，没有什么可补充的。 ...

谢谢sob老师，结构优化和频率计算的输出文件已上传

作者
Author: 吃饭睡觉打豆豆 时间: 2023-4-2 16:35

wzkchem5 发表于 2023-4-2 16:31
试过计算初始Hessian了吗？

关键词写了opt=calcfc，因为用的办公电脑，而且软件为gaussian09，暂时只能这样了。

作者
Author: 吃饭睡觉打豆豆 时间: 2023-4-2 16:52

wzkchem5 发表于 2023-4-2 16:31
试过计算初始Hessian了吗？

不好意思，刚刚发现扭转二面角的初始结构里面关键词没有calcfc 只写了int=ultrafine 我再试一下

作者
Author: sobereva 时间: 2023-4-3 07:07

吃饭睡觉打豆豆发表于 2023-4-2 16:52
不好意思，刚刚发现扭转二面角的初始结构里面关键词没有calcfc 只写了int=ultrafine 我再试一下

按虚频微调结构，opt(calcfc,tight,maxstep=5,notrust) freq int=ultrafine再试

作者
Author: 吃饭睡觉打豆豆 时间: 2023-4-3 08:49

sobereva 发表于 2023-4-3 07:07
按虚频微调结构，opt(calcfc,tight,maxstep=5,notrust) freq int=ultrafine再试

感谢sob老师，我再试一下

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