用orca进行对称性破缺计算时Level Shifting的设置为什么会导致铁磁耦合分析结果相反？

pikachuupup

老师您好，我使用orca在B97-3c下对Fe2S2进行对称性破缺的单点计算时，输出文件显示先尝试高自旋态，然后尝试对称性破缺态，但是输入文件中使用shift shift 1.2 erroff 0 end后会得到完全不同的两种结果。

请问老师：

（1）如果计算完显示铁磁耦合，最后输出文件是高自旋态的能量，还是对称性破缺态的能量呢？

（2）shift shift 1.2 erroff 0 end使用后，bs态的能量不变而只有hs态的能量发生了变化，为什么不是两个一起变呢？有没有什么技巧判断如何合理设置此语句呢？

附上计算文件，谢谢老师。

pikachuupup

不好意思各位老师，这是我刚入门orca时提的疑问，现在回过头来看看不是什么值得问的问题，解答如下供后来人参考：
（1）使用brokensym或flipspin给出的最后结果是反铁磁耦合的结果；
（2）shift shift 1.2 erroff 0 end一般不会影响最后结果，本例结果不同是因为高自旋态收敛到了不同的波函数。