自旋极化单重态的自由能计算流程应该是什么样的

ra2123456

计算某物种自旋极化单重态的吉布斯自由能时,是否一定需要遵循这样的流程,即:
从某个结构出发用guess=mix nosymm stable=opt这组关键词得到极化单重态波函数后再进行结构优化

能否直接将(电荷,自旋多重度)设置为0,1对该物种执行 opt freq +(U)B3LYP呢?
我采取这种方法对一个非过渡态物种最终得到了S**2等于1的结果,对过渡态物种得到了S**2= 0.4851的结果

wzkchem5

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

ra2123456

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

好的,谢谢您

zjxitcc

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

sobereva

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

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

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