VASP计算高熵合金磁性体系（DFT+U）计算结果磁矩偏大

lll000

VASP计算高熵合金磁性体系（DFT+U）计算中，静态自洽求的高熵合金各个原子贡献磁矩，计算值偏大49-68号是铁原子（我求出Fe原子结果为3波尔，别人论文里是1-2），该如何去解决，以下是INCAR,POSCAR

Tosykie

如果不+U计算呢？结果如何。

876449830

没有氧化物，氮化物的形式存在，也可以不加吧

lll000

Tosykie 发表于 2025-2-26 01:54
如果不+U计算呢？结果如何。

这个我还没试过，等我试试再跟您交流

lll000

876449830 发表于 2025-2-26 09:19
没有氧化物，氮化物的形式存在，也可以不加吧

相同方向的论文都加了，所以我也得加

lll000

Tosykie 发表于 2025-2-26 01:54
如果不+U计算呢？结果如何。

我觉得我MAGMOM设置也有问题，U值也有问题，您看看有什么不合适的地方吗

lll000

@卡开发发 大佬可以解答一下吗

卡开发发

lll000 发表于 2025-2-26 10:14
@卡开发发 大佬可以解答一下吗

1、金属体系的U值通常可能比氧化物的情况来得低，我不确定你U值是怎样获取和设置的。
2、金属体系不建议使用ALGO=A。
3、自洽场是否收敛在这没办法判断。
4、LORBIT=11得到的局域磁矩本身依赖WS半径，是否文献一定合理或者你的结果合理这也很难说。
高熵合金第一性原理计算从模型到计算方法都存在不尽人意的问题，这种计算应当慎重进行。

lll000

卡开发发 发表于 2025-2-26 11:43
1、金属体系的U值通常可能比氧化物的情况来得低，我不确定你U值是怎样获取和设置的。
2、金属体系不建议 ...

感谢大佬回答
1.我找不到相关文献的U值，因为高熵合金磁性计算领域好多都是实验+计算，计算过程很不详细，我都是看了其他类型U值，J值，直接COPY的
2.ALGO之前换了V、F、N均不收敛，最后实在没法才换的A
3.换ALGO=A以后，自洽收敛了
4.我不太懂
我觉得我这个体系MAGMOM，U,J值都有问题，后续又重新进行结构优化INCAR，请您再帮忙看一下

卡开发发

lll000 发表于 2025-2-26 12:55
感谢大佬回答
1.我找不到相关文献的U值，因为高熵合金磁性计算领域好多都是实验+计算，计算过程很不详细 ...

1、确定U的方法有很多种，或许你可以采用这些方法。
2、3、ALGO=A并不推荐用在金属上，或许你可以调整为ALGO=D试试，或者调整自洽混合相关的参数，理论上说ALGO=A并非唯一解决方案。
4、各个原子的局域磁矩不是严格实验可测的，电子密度和自旋密度都不是严格空间局域的，因此计算上得到局域磁矩需要通过划分方案才能够得到。

Tosykie

lll000 发表于 2025-2-26 09:50
我觉得我MAGMOM设置也有问题，U值也有问题，您看看有什么不合适的地方吗

可能是初始磁矩的大小和/或方向导致的。磁性材料的磁矩设置很有tricks。

lll000

卡开发发 发表于 2025-2-26 21:29
1、确定U的方法有很多种，或许你可以采用这些方法。
2、3、ALGO=A并不推荐用在金属上，或许你可以调整为 ...

嗯嗯，感谢卡卡大佬，我去试试

lll000

Tosykie 发表于 2025-3-2 23:56
可能是初始磁矩的大小和/或方向导致的。磁性材料的磁矩设置很有tricks。

我的铁钴是3，铬是+-2交替，V不知道有没有磁性，假设一下设置为2，Al直接不设置了，您看这样可以吗

lll000

Tosykie 发表于 2025-2-26 01:54
如果不+U计算呢？结果如何。

不+U计算结果正常，加了磁矩偏大

土土拨鼠

lll000 发表于 2025-2-26 09:42
相同方向的论文都加了，所以我也得加

您好，这个论文方便分享一下吗？