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

woshixuebaxixi

发发老师，我尝试更改参数(磁矩试过5)，Fe2O3的mag还是停留在0.0024 。pdf卡片里面确实存在(110)这个晶面，您能不能帮我看看哪里出问题了

fanzhichao001

卡开发发 发表于 2024-3-5 19:17
要。不过你最好看看跑不起来的原因是什么。

好的，谢谢老师。

卡开发发

woshixuebaxixi 发表于 2024-3-6 14:57
发发老师，我尝试更改参数(磁矩试过5)，Fe2O3的mag还是停留在0.0024 。pdf卡片里面确实存在(110)这个晶面， ...

数值误差是会有的，不过我建议你先检查下你那个DFT+U的设置。

woshixuebaxixi

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

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

卡开发发

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

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

woshixuebaxixi

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

卡开发发

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

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

woshixuebaxixi

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

卡开发发

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

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

woshixuebaxixi

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

卡开发发

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

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

123wo

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

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

卡开发发

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

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

123wo

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

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