铁是铁磁性的，所以单电子自旋平行。如果强行当自旋反平行优化，等于把铁强行搞成顺磁或者反铁磁相来算，结果当然不能用

sobereva 发表于 2023-11-3 08:50
别用量化程序结合簇模型试图代替第一性原理程序计算d族（不含ds族）金属单质，自找罪受

谢谢解答。我同时也在用CP2K跑这个模型，等我再折腾一段时间，再来请教。

别用量化程序结合簇模型试图代替第一性原理程序计算d族（不含ds族）金属单质，自找罪受

建议放弃。这种团簇使用孤立体系的程序目前算不了，即使侥幸有哪次计算normal termination，结果也丝毫不可靠。

（1）“无脑就把金属团簇的原子总数量弄成偶数个，那么该团簇就必然是闭壳层”完全不是这样的。偶数个电子的体系不一定是单重态，也可以是三重态、五重态等等，因此未必是闭壳层。至于到底什么自旋还得看切的方式、边缘饱和方式、体系特点，所以十分麻烦。即使侥幸是单重态（或者你就只想算单重态），则它必然不是闭壳层单重态，而是拥有巨量自由基特征的开壳层单重态，因为这里面每个Fe是开壳层的，要自旋反平行耦合成整体单重态，会有大量未成对电子。DFT无法描述这样的孤立体系，无论哪种泛函，所以不用特意问PBE0。

（2）边缘没有处理或饱和，边界效应过强，离真实周期性体系差距过大。当然，这种体系边缘怎么饱和才合适是个问题。

（3）与O2复合的话，得看你想研究三重态O2，还是研究单重态O2。三重态O2与一个三重态体系复合，整体可以是五重态，也可以是单重态，算了才知道哪个能量低。注意单重态O2是完美双自由基，不是闭壳层的。

（4）这样的孤立体系有大量的SCF解，其中有多个SCF稳定解，不知道该取哪个解进行计算。（这里要黑一下一些人算多核过渡金属连波函数稳定性检验都不做，那这种计算建议彻底不做，节省碳排放）

（5）那么该如何十分严谨地研究这种体系？在遥远的将来使用周期性多参考方法，并考虑多态反应性，考虑多态自旋混合模型。