求助：Fe-Co-Ni-Cr-Cu 合金结构优化scf不收敛

37010405

用面心立方的镍的晶胞，建2x2x2的超胞 然后随机替换原子，差不多每个原子占1/5. 然后用multiwfn生成inp文件做晶胞的结构优化，但是scf不收敛。按照CP2K中遇到SCF难收敛时的解决方法 - 第一性原理 (First Principle) - 计算化学公社贴子说的把k-point 加到5x5x5，cutoff加到800，smear alpha 值和温度调了都没有收敛迹象。
1.考虑是不是磁矩没设置好。没查到这个合金的磁性，猜想应该是铁磁性。设原子的初始磁矩时候应该怎么设置 是按照ms里high spin的值设Fe 4  ；Co 3  ；Ni 2； Cr 5 ；还是设成体相的磁矩 fe2.2 ；Co1.7； ni 0.6 ；Cr 2；或者 把后面这个值乘1.5？
2.  是不是本身结构太小导致的。看见别人的高熵合金晶胞都比较大，我建的这个2x2x2的晶胞只有32个原子，合金化之后每种原子周围都是其他原子，会导致化学意义太差吗。这样的话需要建多大晶胞合适。

不知道 这两个方向对不对？还有没有其他的地方我没设置好的也烦请各位大神告之，不胜感谢。

kaikai

    MULTIPLICITY   89 #Spin multiplicity
第160行里面这个89，这个数值？？？？？？

37010405

kaikai 发表于 2025-3-2 14:29
MULTIPLICITY   89 #Spin multiplicity
第160行里面这个89，这个数值？？？？？？

这个值是multiwfn里设置了原子的 magnetization之后自动生成的值。这个值不合理吗？
我也试过把multiplicity改成成relax_multiplicity  但还是不行。
我记得社长有说过开了smear之后  这个值的设成什么都不重要，不知道我的理解有问题吗