oniom计算不收敛想请教一下

a5667499

用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.

sobereva

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

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

a5667499

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

好的谢谢老师，我试一下

287192454

大概率能量震荡，经常碰到这种事