加了calcall还会发生震荡，请问还有什么方法可以用吗，谢谢老师指教 - 计算化学公社

Gaussian 常见报错及解决方法；新手求助报错时的注意事项 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)
L9999

Error termination via Lnk1e in l9999.exe
复制代码
此报错需查看输出文件中的额外信息，自结尾处向上检索“Optimization stop”可见如下几种情况

Optimization stopped.
-- Wrong number of Negative eigenvalues: Desired= 1 Actual= 4
-- Flag reset to prevent archiving.
复制代码
成因：Gaussian会在过渡态优化时自动检查虚频的数量，若不为1即终止。这一检查没有必要
解决：将NoEigenTest关键词加入到opt的选项中，如 opt=(TS, CalcFC, NoEigenTest)

Optimization stopped.
-- Number of steps exceeded, NStep= 100
-- Flag reset to prevent archiving.
复制代码
成因：几何优化未能在指定步数内收敛
解决：打开GaussView，File-Open，按照如下设置：

然后点击Result-Optimization，可以看到优化时的能量变化曲线。观察优化曲线应纵向放大，方法为用鼠标反复选择纵轴方向仅包含最末几个点、但横轴方向包含所有点的“宽而矮”的矩形范围，直到能看清最末几个点的变化情况为止。根据图形判断是否震荡。如下图，需纵向放大至右边才能看清确实震荡了。（注意图中偶尔有突然的、无“周期性的”某1~2个点的spike是正常现象，不是震荡，下方小幅度但有“周期性”的才是震荡。）

左侧看不出是否震荡。需纵向放大至右侧才能看出。

如果没有震荡：opt=maxcycle=当前步数的2~3倍。可以读取当前最末/最优的一个构象做为新的初猜。注意检查构象、轨迹是否合理，如果结构已经不合理了，则应适当调整初猜后再做（尤其是找过渡态的时候）。

如果震荡了：http://sobereva.com/164。

文中的众多解决方法中，我个人倾向于找到能量较低、4个收敛标准离收敛限较近，且结构合理的某点为新的初猜，并首先尝试在 opt 选项中加上(MaxStep=5, NoTrustUpdate, GDIIS) 选项；如果有计算频率的资源，可再加上opt=calcfc选项。特别仅对于对称性较高的结构，如T、C3 等群的复杂分子，若反复调整关键词后仍然震荡，可以尝试构建或破坏具有对称性的初猜结构，分别结合opt=Cartesian 优化，如 opt=(MaxStep=5,NoTrustUpdate,Cartesian), 有时有效果。
注意上一段仅是我个人首先尝试的建议，不要仅仅试了这一个不行就又把问题拿出来问。http://sobereva.com/164 有的是其他办法。

乐平发表于 2024-7-22 20:47
先把基组里的弥散函数给去掉，也就是改成 6-31G(d,p), 甚至对 H 的极化也可以先不用。也就是直接 6-31G(d ...

藤井树发表于 2024-7-22 16:53

zlisa 发表于 2024-7-22 20:57
谢谢老师！请问老师：同一种体系的其他构型用这个方法优化完成了的还需要这样改吗，如果不改会有影响吗？ ...

乐平发表于 2024-7-22 21:19
我前面说了啊，先用小基组优化完，然后再读优化后的 .chk 文件

zlisa 发表于 2024-7-23 13:32
谢谢老师！已经按照您的方法在重新优化啦！
但是我之前有一些用之前的方法优化完成的构型，它们也需要重 ...

藤井树发表于 2024-7-22 16:53

sobereva 发表于 2024-7-26 05:44
又不是算热力学量，写scale毫无意义

M06-2X比其它大多数泛函更容易震荡不收敛，遇到难收敛时避免用M06-2 ...