计算化学公社

标题: CP2K变胞优化亚铁磁尖晶石Fe3S4如何解决磁矩引起的不收敛问题 [打印本页]

作者
Author: TCN 时间: 2026-6-6 21:42
标题: CP2K变胞优化亚铁磁尖晶石Fe3S4如何解决磁矩引起的不收敛问题
CP2K变胞优化亚铁磁尖晶石Fe3S4单胞，尝试MP给出的磁矩以及采用&BS字段设置Fe2+和Fe3+均不收敛，该类尖晶石结构应该如何有效去尝试不同磁矩去收敛，inp如下

作者
Author: sobereva 时间: 2026-6-7 08:15
几何优化任务不收敛的时候始终自觉说清楚是SCF不收敛还是几何优化不收敛

假设是前者，记得初猜原子自旋磁矩不是影响SCF能否收敛的唯一因素，该怎么解决下文该说的都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）
MP那数据库里很多原子自旋磁矩都是错的，不要盲目以那个为参考，应该找相关的靠谱的文献里的信息

作者
Author: 阿夸夸 时间: 2026-6-7 11:07
本帖最后由阿夸夸于 2026-6-7 11:27 编辑

我这里按照反尖晶石（inverse spinel）结构，
8个四面体Fe为3+（3d5）自旋向下，16个八面体为Fe为3+（3d5）自旋向上（但其实是一半Fe为3+（3d5）一半Fe为2+（3d6），这里只是初猜不影响）的情况设置了&BS字段，
自旋多重度为（8个四面体8*5=40，16个八面体一半Fe3+，一半Fe2+即8*5+8*4=72，亚铁磁四面体与八面体自旋相反，72-40+1=33）33，设置了输入文件。
昨天结构优化了5步卡住了

，但从输出文件上看前几步确实scf不收敛但之后scf收敛趋势已经很明显了，继续跑scf收敛应该不是问题。

作者
Author: TCN 时间: 2026-6-7 11:55

阿夸夸发表于 2026-6-7 11:07
我这里按照反尖晶石（inverse spinel）结构，
8个四面体Fe为3+（3d5）自旋向下，16个八面体为Fe为3+（3d5 ...

感谢老师百忙中的回复，学到了不少

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