ORCA6.0中使用wB97M(2)计算单重态开壳层氧气分子单点能出错求助

MCF

       我用wB97M(2)计算三重态和单重态开壳层氧气分子的单点能的时候发现一个问题，就是算单重态开壳层氧气发现进行第二步wB97M(2)计算的时候ORCA6.0好像忽略了一半的电子数，体系只有8个电子被填充在轨道内，而且能量也差很多，第一步的wB97M-V计算出来是-150多Hartree，但是第二步wB97M(2)算出来是-80多Hartree。三重态计算一切正常，还请各位大佬指点一下是哪里出错了。
       下面分别是输入文件和部分计算结果，在算单点能之前两种电子态的氧气分子已经经过B3LYP-D4/def2-TZVP几何优化过了。
单重态开壳层氧气分子单点能输入文件

单重态wB97M-V计算能量结果

单重态wB97M(2)计算能量结果

单重态电子占据轨道情况

三重态输入文件

三重态wB97M-V计算能量结果

三重态wB97M(2)计算能量结果

zjxitcc

你自己说“单重态开壳层氧气分子”，但帖子内容似乎没看到你构建了对称破缺初猜，也没看到检验单重态UKS波函数稳定性，也没看到读取对称破缺UKS轨道，如何证明你算的是开壳层单重态？别用双杂化算了半天，发现概念搞错了。
后记：我看到guessmix了，你是想通过它来构建对称破缺初猜？

MCF

zjxitcc 发表于 2024-12-24 20:34
你自己说“单重态开壳层氧气分子”，但帖子内容似乎没看到你构建了对称破缺初猜，也没看到检验单重态UKS波 ...

图中单重态输入文件里面关键词加了guessmix，而且wB97M-V计算S平方算符的期望值接近于1，应该能说明算的是对称性破缺的电子态吧。

zjxitcc

MCF 发表于 2024-12-24 20:51
图中单重态输入文件里面关键词加了guessmix，而且wB97M-V计算S平方算符的期望值接近于1，应该能说明算的 ...

嗯，这看起来是合理的。建议在双杂化的任务里也加上UKS关键词，看是否还有你帖子里描述的问题。
另外，为何要写MaxIter 1和IgnoreConv 1？

MCF

zjxitcc 发表于 2024-12-24 21:03
嗯，这看起来是合理的。建议在双杂化的任务里也加上UKS关键词，看是否还有你帖子里描述的问题。
另外 ...

MaxIter 1和IgnoreConv 1这俩关键词是手册原文里提到的，应该有特殊作用。

MCF

zjxitcc 发表于 2024-12-24 21:03
嗯，这看起来是合理的。建议在双杂化的任务里也加上UKS关键词，看是否还有你帖子里描述的问题。
另外 ...

加了uks以后结果正常了，感谢。