计算化学公社

标题: 请教一下高斯的opt+freq时，力已经收敛到接近0，但freq仍然有虚频应当怎样处理？ [打印本页]

作者
Author: rugals 时间: 2022-11-19 11:27
标题: 请教一下高斯的opt+freq时，力已经收敛到接近0，但freq仍然有虚频应当怎样处理？
如题，体系为无共轭的柔性聚合物分子链。大约20个链节。b3lyp/6-311g** em=gd3bj
第一次优化时使用calcfc，最大力和rms力都已经接近0，位移没收敛，软件认定收敛。但freq有一个绝对值小于-1的虚频。
第二次优化在第一次基础上加了recalc=5,maxstep=5,notrust，结果同上。
第三次优化在第二次基础上改成了calcall,gdiis,maxstep=3,notrust, 结果仍然同上。opt只跑了一步就因为力接近0而认定收敛，故freq仍然有虚频。
第四次优化在第三次基础上又加了tight或very tight。但用时甚巨，一天只能跑10步左右。看优化曲线，四个指标反而跑得越来越远了。如下图。
(, 下载次数 Times of downloads: 10)
也试过了在虚频振动方向上调节结构，依然不起作用。(优化的结构是螺旋形链，虚频方向是螺旋截面的顺/逆时针扭转方向)

请问这种势能面太缓的结构应该如何消掉虚频？另外calcall搭配maxstep+notrust是否会起到反作用？因为我之前单用recalc时也不收敛，得recalc+maxstep+notrust才会收敛。

作者
Author: 七尺贱 时间: 2022-11-19 12:04
你把em=gd3bj改成em=gd3去算freq有奇效

作者
Author: rugals 时间: 2022-11-19 12:26

七尺贱发表于 2022-11-19 12:04
你把em=gd3bj改成em=gd3去算freq有奇效

谢谢老师，但可以这么做吗？
按理说opt和freq的理论方法得保持一致才可以吧？去了阻尼势能面不就变了？

作者
Author: sobereva 时间: 2022-11-19 12:58
对于G09 D.01之后的版本，用em=gd3bj完全没问题，而且推荐用em=gd3bj。对于G09 D.01，em=gd3bj算的Hessian有问题，此时应当opt freq全用em=gd3

不上传输出文件没法判断

对于震荡导致的难收敛，calcall和maxstep+notrust结合没有反效果。

作者
Author: rugals 时间: 2022-11-19 13:22

sobereva 发表于 2022-11-19 12:58
对于G09 D.01之后的版本，用em=gd3bj完全没问题，而且推荐用em=gd3bj。对于G09 D.01，em=gd3bj算的Hessian ...

谢谢老师，我用的g16；
第2-4次优化的log文件附上，麻烦sob老师看下。(第4次还没跑完，目前21步)

作者
Author: 七尺贱 时间: 2022-11-19 16:21

sobereva 发表于 2022-11-19 12:58
对于G09 D.01之后的版本，用em=gd3bj完全没问题，而且推荐用em=gd3bj。对于G09 D.01，em=gd3bj算的Hessian ...

我的版本就是G09 D.01,发现有小虚频消不掉，发em=gd3bj换成em=gd3就没虚频了

作者
Author: sobereva 时间: 2022-11-19 22:53

rugals 发表于 2022-11-19 13:22
谢谢老师，我用的g16；
第2-4次优化的log文件附上，麻烦sob老师看下。(第4次还没跑完，目前21步)

加上int=superfine

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