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

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

收到，谢谢老师！！

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

收到，谢谢老师！

又不是算热力学量，写scale毫无意义

M06-2X比其它大多数泛函更容易震荡不收敛，遇到难收敛时避免用M06-2X
calcall不解决所有问题，该用小步长上限用小步长上限，下文都说了。文中又不是只说了用calcall，也没说calcall能保证100%收敛
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
没事甭给氢乱加弥散，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

本来就难达到收敛限居然还用tight，雪上加霜

别光看这种折线图，始终结合优化过程的轨迹的动画判断情况。

l9999 不是报错，只是 Gaussian 的一个模块，提示你任务完成了

看不到结构，什么关键词都没用，要以结构为主！！!, 发文件！

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

前面的方法不需要再同样小基组优化，然后都chk大基组再算一遍。能用大基组一遍过一样的。

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

方法和基组要一致才能够比较

我说的意思是，你这个结构先用小基组 6-31G(d) 优化，得到结果之后，再读取优化后的 .chk 文件，用你现在的大基组再优化。
其余的不用重算

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

谢谢老师！已经按照您的方法在重新优化啦！
但是我之前有一些用之前的方法优化完成的构型，它们也需要重新用小基组优化再读chk用大基组优化吗？
如果不改的话，这两种方法优化（一种是我帖子里提到的方法，是我之前一直用的，另一个就是您说的小基组先优化-再读chk文件的方法）得到的构型能一起做比较吗？不同的方法对优化的结果会有影响吗？
（如果有说的不对的地方还请老师赐教，原谅我是个刚开始摸索的菜鸟，真的麻烦老师了）

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

我前面说了啊，先用小基组优化完，然后再读优化后的 .chk 文件

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

谢谢老师！

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

谢谢老师！请问老师：同一种体系的其他构型用这个方法优化完成了的还需要这样改吗，如果不改会有影响吗？谢谢老师解答~~~

先把基组里的弥散函数给去掉，也就是改成  6-31G(d,p),  甚至对 H 的极化也可以先不用。也就是直接 6-31G(d)，先用小基组优化好之后，再读优化后的 .chk 文件。

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 有的是其他办法。