计算化学公社

标题: oniom计算不收敛想请教一下 [打印本页]

作者
Author: a5667499 时间: 2024-2-27 19:39
标题: oniom计算不收敛想请教一下
用Gaussian计算oniom模型进行结构优化时，低层原子固定，用GV看optmization基本都趋于平稳，为什么会出现这种情况，差距这么大，真心向各位老师请教一下，感谢各位老师。

log文件显示
      Item             Value    Threshold  Converged?
Maximum Force          0.250936    0.000450    NO
RMS    Force          0.051033    0.000300    NO
Maximum Displacement    0.000362    0.001800    YES
RMS    Displacement    0.000116    0.001200    YES
Maximum MM Force       0.000000    0.000450    YES
RMS    MM Force       0.000000    0.000300    YES
Predicted change in Energy=-5.881492D-05

Rotational constants (GHZ):          0.0018502          0.0014548          0.0011145
Leave Link  202 at Tue Feb 27 17:48:47 2024, MaxMem= 104857600 cpu:             0.1 elap:             0.0
(Enter /opt/apps/G16-C01-AVX/g16/l9999.exe)

I am not a vegetarian because I love animals;
I am a vegetarian because I hate plants.
                           -- A. Whitney Brown
Error termination request processed by link 9999.
Error termination via Lnk1e in /opt/apps/G16-C01-AVX/g16/l9999.exe at Tue Feb 27 17:48:48 2024.

作者
Author: sobereva 时间: 2024-2-28 06:12
观察收敛情况窗口把末期拉大了看，仔细观察优化轨迹弄清楚情况，并且结合下文说的尝试帮助收敛
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

这种体系用ONIOM也没多大意义，直接用CP2K当周期性体系优化省心得多。就算是要用Gaussian算高层区域的反应之类，把当前体系截1/3或一半，冻结边界原子并饱和，不用ONIOM直接做DFT计算，也一点不费劲

作者
Author: a5667499 时间: 2024-2-28 18:56

sobereva 发表于 2024-2-28 06:12
观察收敛情况窗口把末期拉大了看，仔细观察优化轨迹弄清楚情况，并且结合下文说的尝试帮助收敛
量子化学计 ...

好的谢谢老师，我试一下

作者
Author: 287192454 时间: 2024-2-29 08:55
大概率能量震荡，经常碰到这种事

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