计算化学公社

标题: 为何opt(calcall)+freq的输出结果中，优化完成一步收敛了但频率计算时未收敛？ [打印本页]

作者
Author: 镜天心 时间: 2024-8-31 20:31
标题: 为何opt(calcall)+freq的输出结果中，优化完成一步收敛了但频率计算时未收敛？
各位老师好，我在使用opt(calcall)+freq优化时，由于使用的是calcall，因此在优化完毕的同时进行了频率运算，此时四个收敛限都是Yes。而额外进行进一步的freq时，最大位移和均方根位移有了显著的偏移，高出收敛限很多。calcall是在每步都精确计算Hessian矩阵的，按理说优化完成后的结果是严格收敛的。但这里freq的结果与优化后的结果明显不同，这是什么原因，哪个结果是正确的，还请老师们指点。

作者
Author: Huschein 时间: 2024-9-2 05:17
你这force不都是一样的吗，位移是和上一步比较吧？freq没有上一步，我也不知道位移咋来的，但force一致的话应该就是收敛了

作者
Author: Loading0760 时间: 2024-9-2 14:06
用的是什么版本的

作者
Author: mfdsrax2 时间: 2024-9-2 15:17
看这里：
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

PS2：有人问，为什么有的时候优化时明明用了calcall，但是基于其结构再做freq任务的时候显示的判断标准却依然和优化最后一步不同，甚至又出现了NO。可能原因有三 (1)优化时用的级别、数值设定（如DFT积分格点）等因素和几何优化时不完全一致 (2)用GaussView读取优化的输出文件，然后又保存成freq任务的输入文件的时候，由于小数位数有限，因此造成了一点数值误差 (3)几何优化的时候用了GDIIS（默认的GEDIIS也有一定GDIIS的成份），由于这种方法预测下一步位移的时候还会参考之前步的信息，因此会和freq时基于牛顿法判断出的位移有所不同。

作者
Author: 镜天心 时间: 2024-9-2 15:29

Loading0760 发表于 2024-9-2 14:06
用的是什么版本的

Gaussian 16, Revision A.03

作者
Author: 镜天心 时间: 2024-9-2 15:31

mfdsrax2 发表于 2024-9-2 15:17
看这里：
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

好的，非常感谢。所以说位移的区别是源于GEDIIS内的GDIIS与牛顿法不一致而导致的是吗？

作者
Author: mfdsrax2 时间: 2024-9-2 16:05

镜天心发表于 2024-9-2 15:31
好的，非常感谢。所以说位移的区别是源于GEDIIS内的GDIIS与牛顿法不一致而导致的是吗？

是的

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