负载过渡金属计算磁矩问题

Deepast

计算一个二维材料负载单个3d金属原子，超胞大概12Ax10Ax20A，ISPIN=2。目前想研究小分子吸附和进行电子结构分析。因为资源有限，做了下测试
收敛性测试（Ni原子为例）发现Gamma 3x3x1 的时候相对能量变化已经小于10meV/atom了。但K点取到551，计算的磁矩与高精度计算（991）相比，也不见收敛。
也测试过展宽和赝势，对磁矩影响比较小，主要问题在于k点。
让我不解的是我不知道这个磁矩对电子结构有什么影响，收不收敛是否重要。
也看过文献，有趣的是明明是同一种材料，文献上dmol3和vasp计算的磁矩完全不同，甚至定性（有无磁性）都不一样。

876449830

INCAR、POSCAR贴出来

Deepast

INCAR POSCAR在这
ISMEAR=0 也试过了@876449830

876449830

Deepast 发表于 2022-6-14 21:58
INCAR POSCAR在这
ISMEAR=0 也试过了@876449830

MAGMOM注释掉，ISMEAR=0，IBRION=-1，EDIFFG=-0.05，ISIF=2，试试

Deepast

876449830 发表于 2022-6-15 08:45
MAGMOM注释掉，ISMEAR=0，IBRION=-1，EDIFFG=-0.05，ISIF=2，试试

单点计算
441 mag=0.3089
551 mag=0.4007
991 mag=0.4975
我不太确定这些磁矩对电子态影响有多大，尝试别的体系磁矩对K点似乎没这么敏感
另外我最开始试过不设ALGO=Fast和MAGMOM计算出无磁矩的，但没有输出相关的文件（包括LORBIT=11）。复现的时候发现一定要设置ALGO和MAGMOM之一才不报错，然后试过各种初猜也没找到无磁矩的。

876449830

Deepast 发表于 2022-6-15 14:51
单点计算
441 mag=0.3089
551 mag=0.4007

不应该呐，设置ISPIN=2、LORBIT=10or11的话OUTCAR里肯定有magnetization(x)这部分输出的，K点密度看上去对磁矩有一定影响，但是测试一下K点和截断能即可。

Deepast

876449830 发表于 2022-6-15 15:07
不应该呐，设置ISPIN=2、LORBIT=10or11的话OUTCAR里肯定有magnetization(x)这部分输出的，K点密度看上去 ...

不好意思，是我没表达清楚。刚开始的时候我还没有把LORBIT加上去，发现有问题后才加的（所以最开始那个无磁矩只有结构文件）。
我再测试下ENCUT试试，可能是其它因素没收敛

876449830

Deepast 发表于 2022-6-15 16:07
不好意思，是我没表达清楚。刚开始的时候我还没有把LORBIT加上去，发现有问题后才加的（所以最开始那个无 ...

截断能和K点密度要先测试好再去进行后续工作的。