计算化学公社

标题: VASP计算高熵合金磁性体系（DFT+U）计算结果磁矩偏大 [打印本页]

作者
Author: lll000 时间: 2025-2-25 21:24
标题: VASP计算高熵合金磁性体系（DFT+U）计算结果磁矩偏大
VASP计算高熵合金磁性体系（DFT+U）计算中，静态自洽求的高熵合金各个原子贡献磁矩，计算值偏大49-68号是铁原子（我求出Fe原子结果为3波尔，别人论文里是1-2），该如何去解决，以下是INCAR,POSCAR
(, 下载次数 Times of downloads: 4)

(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 4)

作者
Author: Tosykie 时间: 2025-2-26 01:54
如果不+U计算呢？结果如何。

作者
Author: 876449830 时间: 2025-2-26 09:19
没有氧化物，氮化物的形式存在，也可以不加吧

作者
Author: lll000 时间: 2025-2-26 09:42

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

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

作者
Author: lll000 时间: 2025-2-26 09:42

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

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

作者
Author: lll000 时间: 2025-2-26 09:50

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

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

作者
Author: lll000 时间: 2025-2-26 10:14
@卡开发发大佬可以解答一下吗

作者
Author: 卡开发发 时间: 2025-2-26 11:43
本帖最后由卡开发发于 2025-2-26 11:44 编辑

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

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

作者
Author: lll000 时间: 2025-2-26 12:55

卡开发发发表于 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，请您再帮忙看一下

作者
Author: 卡开发发 时间: 2025-2-26 21:29

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

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

作者
Author: Tosykie 时间: 2025-3-2 23:56

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

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

作者
Author: lll000 时间: 2025-3-5 11:43

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

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

作者
Author: lll000 时间: 2025-3-5 11:45

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

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

作者
Author: lll000 时间: 2025-3-20 17:18

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

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

作者
Author: 土土拨鼠 时间: 2025-5-23 15:32

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

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

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