求助：a-Fe2O3（001）面的结构优化问题

fanzhichao001

求助：如图是对a-Fe2O3的（001）面进行切面、扩胞、建真空层后的结构，对此结构用以下INCAR参数进行结构优化。优化后OSZICAR的最后一步mag=62.7321正常吗？a-Fe2O3是反铁磁体系，总磁矩不应该为0嘛？

梁雪彤

WilliamH 发表于 2025-3-12 11:17
如果是晶胞中同种原子存在不同的自旋朝向，在VESTA中能够直观的看到原子的配位情况，很容易区分自旋向上/ ...

好的好的，感谢感谢

WilliamH

梁雪彤 发表于 2025-3-12 09:43
我没太看懂，为什么要将磁性向下的进行替换成别的原子，可以麻烦您解释一下吗，我算的也是a-Fe2O3,算出来的 ...

如果是晶胞中同种原子存在不同的自旋朝向，在VESTA中能够直观的看到原子的配位情况，很容易区分自旋向上/向下的原子，但是在POSCAR中就很难根据坐标来区分。在VESTA中建模的时候就替换掉其中一种自旋朝向的原子，然后在导出的POSCAR中直接把元素符合改一下就好了。   α-Fe2O3在论坛中应该有过很多讨论了，可以先搜索之前的相关帖子看看。

梁雪彤

我没太看懂，为什么要将磁性向下的进行替换成别的原子，可以麻烦您解释一下吗，我算的也是a-Fe2O3,算出来的总磁矩也不对。晶胞优化完是0，用优化完得切面扩胞后再进行结构优化就不对了，而且感觉很难收敛，请问有什么解决办法吗

123wo

卡开发发 发表于 2024-4-13 02:14
会，随着接近收敛应该会越算越快。

感谢O(∩_∩)O哈哈~

卡开发发

123wo 发表于 2024-4-12 10:39
卡老师，请问如果不能正确的设置初始磁矩，是否会导致计算速度显著减慢呢

会，随着接近收敛应该会越算越快。

123wo

卡开发发 发表于 2024-2-28 18:14
不正常。属于不设置MAGMOM的一个后果。

卡老师，请问如果不能正确的设置初始磁矩，是否会导致计算速度显著减慢呢

卡开发发

woshixuebaxixi 发表于 2024-4-3 21:47
好滴发发老师 还有一个问题想请教一下 vasp用ISIF=2理论上可算到自旋相反的电子 我该怎么确定某一个原子的 ...

严格来说是没办法真正确定相互作用中多原子体系某个原子的电子结构，当然通过一些布居分析（也就是划分/投影到局域空间上）粗略看看可以。

woshixuebaxixi

好滴发发老师 还有一个问题想请教一下 vasp用ISIF=2理论上可算到自旋相反的电子 我该怎么确定某一个原子的电子组态呢？您有相关文章或教程推荐吗？

卡开发发

woshixuebaxixi 发表于 2024-3-7 15:12
发发老师，我是直接用vaspkit103生成的POTCAR文件 想问一下您说的这个PAW-PBE应该怎么修改呢？

我不太清楚，我一般不太用vaspkit，我自己一般是自己写一套ase的工作流。

woshixuebaxixi

发发老师，我是直接用vaspkit103生成的POTCAR文件 想问一下您说的这个PAW-PBE应该怎么修改呢？

卡开发发

woshixuebaxixi 发表于 2024-3-7 14:55
发发老师，“你的POTCAR用的是LDA（CA-PZ）的，但实际上指定的交换-关联泛函是GGA=PE” 是不是直接在INCAR ...

把你的POTCAR换成PAW-PBE的，这样不写GGA=PE默认用的也是PBE，别用LDA/CA-PZ算。

woshixuebaxixi

发发老师，“你的POTCAR用的是LDA（CA-PZ）的，但实际上指定的交换-关联泛函是GGA=PE” 是不是直接在INCAR里面注释掉GGA=PE这句话即可？

卡开发发

woshixuebaxixi 发表于 2024-3-7 10:22
发发老师 我尝试关闭Umag变成了0.0017

应该可以忽略，能量差别和之前的不太大。另外Fe2O3需要DFT+U，建议进行设置。
另外，你的POTCAR用的是LDA（CA-PZ）的，但实际上指定的交换-关联泛函是GGA=PE，建议不要这么做。

woshixuebaxixi

卡开发发 发表于 2024-3-7 10:17
数值误差是会有的，不过我建议你先检查下你那个DFT+U的设置。

发发老师 我尝试关闭Umag变成了0.0017