Daniel_Arndt 发表于 2022-10-11 00:15
稳妥一点确实该这样。我个人为了省时间，都是最后才stable=opt的。

我个人经验是稍微大一点的分子，用 ...

嗯，确实。guess=mix出现这种情况实际上是由于高斯只能mix HOMO和LUMO，但是有时最低的解需要mix HOMO-1和LUMO，或者初猜的HOMO本身不靠谱；guess=fragment其实也可能不稳定的，但通常出现在片段对称性比较高的时候（比如只有一个原子），大一点的分子通常不会。

hebrewsnabla 发表于 2022-10-10 19:32
第一个问题就是由三楼提到的原因导致的（opt和stable=opt不能同时做），虽然看起来不像。

是的，opt和stable=opt不能同时出现

hebrewsnabla 发表于 2022-10-10 16:22
第二步就应该加上stable=opt

稳妥一点确实该这样。我个人为了省时间，都是最后才stable=opt的。

我个人经验是稍微大一点的分子，用guess=mix时容易出现输出文件中的自旋布居看起来对称破缺，但实际上波函数不稳定的情况。而用guess=(fragment=2)的时候，很少出现类似情形。如果第一步用的是guess=mix的话，第二步确实应该加stable=opt。

fatpig 发表于 2022-10-10 19:15
没人回答第一个问题？
这明显是没读到chk里的结构...

第一个问题就是由三楼提到的原因导致的（opt和stable=opt不能同时做），虽然看起来不像。

没人回答第一个问题？
这明显是没读到chk里的结构...

Daniel_Arndt 发表于 2022-10-10 11:49
你如果真的想在自旋极化单重态的几何优化中加隐式溶剂模型的话，有点复杂。

建议第一步先用“guess=(fra ...

第二步就应该加上stable=opt

你如果真的想在自旋极化单重态的几何优化中加隐式溶剂模型的话，有点复杂。

建议第一步先用“guess=(fragment=2)”，不加溶剂模型。第二步用“guess=read”，加上溶剂模型。第三步用“guess=read”，加上溶剂模型，加上几何优化的关键词。几何优化完成后是第四步，用“guess=read”，加上溶剂模型，加上“stable=opt”。

zjxitcc 发表于 2022-10-10 09:45
先算单点，检查波函数稳定性，然后读取稳定波函数做结构优化/找过渡态。没有“在opt过程中保证每一个结构 ...

收到非常清楚

chuan437 发表于 2022-10-10 09:27
谢谢，但是我还想知道在opt时如何控制波函数的稳定性

先算单点，检查波函数稳定性，然后读取稳定波函数做结构优化/找过渡态。没有“在opt过程中保证每一个结构的波函数均稳定”的关键词，也没法用一行关键词搞定；即使有，其耗时甚巨，大多数人无法承受。保险起见，结构优化/找过渡态完，还可对最后一帧结构再检测一次波函数稳定性。
另外，用guess=mix stable=opt时最好加上nosymm。

liyuanhe211 发表于 2022-10-10 08:53
stable=opt单点生成一个正确发生自旋极化的波函数。然后在opt的时候读取该波函数作为初猜优化，opt时不加st ...

谢谢，但是我还想知道在opt时如何控制波函数的稳定性

wzkchem5 发表于 2022-10-10 04:27
能量升高是正常的。如果你的初始结构在垂直于反应坐标的方向上优化得很好，但在反应坐标方向上很差，那么过 ...

好的，非常感谢你解答了我的疑惑。

stable=opt单点生成一个正确发生自旋极化的波函数。然后在opt的时候读取该波函数作为初猜优化，opt时不加stable=opt

能量升高是正常的。如果你的初始结构在垂直于反应坐标的方向上优化得很好，但在反应坐标方向上很差，那么过渡态能量比你初始结构能量高，优化时能量自然会升高。之所以你以前没怎么见过，是因为垂直于反应坐标的自由度数（3N-7，N为原子数）远多于反应坐标的自由度数（1），所以如果是随手画的一个结构，那么垂直于反应坐标的方向上的总误差几乎一定大于反应坐标的方向上的总误差，导致初始结构能量高于过渡态能量。