计算化学公社

标题: 求助：Fe-Co-Ni-Cr-Cu 合金结构优化scf不收敛 [打印本页]

作者
Author: 37010405 时间: 2025-3-2 13:10
标题: 求助：Fe-Co-Ni-Cr-Cu 合金结构优化scf不收敛
本帖最后由 37010405 于 2025-3-2 13:17 编辑

用面心立方的镍的晶胞，建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个原子，合金化之后每种原子周围都是其他原子，会导致化学意义太差吗。这样的话需要建多大晶胞合适。

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

作者
Author: kaikai 时间: 2025-3-2 14:29
MULTIPLICITY 89 #Spin multiplicity
第160行里面这个89，这个数值？？？？？？

作者
Author: 37010405 时间: 2025-3-2 14:37

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)