先仔细把此文看了
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

通常建议先参考上文，令初猜结构收敛到自旋极化单重态波函数，然后opt freq任务时用guess=read从chk里读取其波函数当初猜，并检查最终优化完的结构合理性，以及最后一步的波函数是否还是自旋极化波函数。如果是，且此结构也没虚频，就可以用了

opt freq guess=mix一次性做也不是不可以，但光靠guess=mix往往在初始结构下并不能收敛到应有的自旋极化单重态波函数，使得之后的优化和振动分析大概率没意义。

我在论坛上写过非常详细的计算步骤，可以参考http://bbs.keinsci.com/forum.php ... 57&fromuid=2632
（虽说计算步骤不唯一，但有一套写好的详细步骤，方便参考和讨论）

wzkchem5 发表于 2026-4-22 17:03
只要最后的结构既满足结构优化的收敛标准（且虚频数目正确），波函数又是自旋极化单重态，就行。达到这个目 ...

好的,谢谢您

只要最后的结构既满足结构优化的收敛标准（且虚频数目正确），波函数又是自旋极化单重态，就行。达到这个目的的方法不重要。也就是说，可以先stable=opt再结构优化到收敛，但要检查收敛的结构的波函数仍然是自旋极化单重态；可以先结构优化再stable=opt，但要检查stable=opt找到的波函数是否和stable=opt之前不同，如果不同，需要继续用新的波函数做结构优化至收敛，再检查波函数，等等。只要最后该检查的都检查了，不同方法只有需要的人力和机时不同，或者说返工次数不同的区别