为何opt(calcall)+freq的输出结果中，优化完成一步收敛了但频率计算时未收敛？

镜天心 · 发表于 Post on 2024-8-31 20:31:29

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

Huschein · 发表于 Post on 2024-9-2 05:17:45

你这force不都是一样的吗，位移是和上一步比较吧？freq没有上一步，我也不知道位移咋来的，但force一致的话应该就是收敛了

Loading0760 · 发表于 Post on 2024-9-2 14:06:01

用的是什么版本的

mfdsrax2 · 发表于 Post on 2024-9-2 15:17:27

看这里：
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

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

镜天心 · 发表于 Post on 2024-9-2 15:29:46

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

Gaussian 16, Revision A.03

镜天心 · 发表于 Post on 2024-9-2 15:31:38

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

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

mfdsrax2 · 发表于 Post on 2024-9-2 16:05:00

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

是的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 为何opt(calcall)+freq的输出结果中，优化完成一步收敛了但频率计算时未收敛？