计算化学公社

 找回密码 Forget password
 注册 Register

求助:a-Fe2O3(001)面的结构优化问题

查看数: 1133 | 评论数: 28 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-2-28 16:23

正文摘要:

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

回复 Reply

123wo 发表于 Post on 2024-4-13 10:11:10
卡开发发 发表于 2024-4-13 02:14
会,随着接近收敛应该会越算越快。

感谢O(∩_∩)O哈哈~
卡开发发 发表于 Post on 2024-4-13 02:14:04
123wo 发表于 2024-4-12 10:39
卡老师,请问如果不能正确的设置初始磁矩,是否会导致计算速度显著减慢呢

会,随着接近收敛应该会越算越快。
123wo 发表于 Post on 2024-4-12 10:39:44
卡开发发 发表于 2024-2-28 18:14
不正常。属于不设置MAGMOM的一个后果。

卡老师,请问如果不能正确的设置初始磁矩,是否会导致计算速度显著减慢呢
卡开发发 发表于 Post on 2024-4-4 11:02:17
woshixuebaxixi 发表于 2024-4-3 21:47
好滴发发老师 还有一个问题想请教一下 vasp用ISIF=2理论上可算到自旋相反的电子 我该怎么确定某一个原子的 ...

严格来说是没办法真正确定相互作用中多原子体系某个原子的电子结构,当然通过一些布居分析(也就是划分/投影到局域空间上)粗略看看可以。
woshixuebaxixi 发表于 Post on 2024-4-3 21:47:01
好滴发发老师 还有一个问题想请教一下 vasp用ISIF=2理论上可算到自旋相反的电子 我该怎么确定某一个原子的电子组态呢?您有相关文章或教程推荐吗?
卡开发发 发表于 Post on 2024-4-3 15:36:28
woshixuebaxixi 发表于 2024-3-7 15:12
发发老师,我是直接用vaspkit103生成的POTCAR文件 想问一下您说的这个PAW-PBE应该怎么修改呢?

我不太清楚,我一般不太用vaspkit,我自己一般是自己写一套ase的工作流。
woshixuebaxixi 发表于 Post on 2024-3-7 15:12:56
发发老师,我是直接用vaspkit103生成的POTCAR文件 想问一下您说的这个PAW-PBE应该怎么修改呢?
卡开发发 发表于 Post on 2024-3-7 15:01:54
woshixuebaxixi 发表于 2024-3-7 14:55
发发老师,“你的POTCAR用的是LDA(CA-PZ)的,但实际上指定的交换-关联泛函是GGA=PE” 是不是直接在INCAR ...

把你的POTCAR换成PAW-PBE的,这样不写GGA=PE默认用的也是PBE,别用LDA/CA-PZ算。
woshixuebaxixi 发表于 Post on 2024-3-7 14:55:23
发发老师,“你的POTCAR用的是LDA(CA-PZ)的,但实际上指定的交换-关联泛函是GGA=PE” 是不是直接在INCAR里面注释掉GGA=PE这句话即可?
卡开发发 发表于 Post on 2024-3-7 12:13:20
woshixuebaxixi 发表于 2024-3-7 10:22
发发老师 我尝试关闭Umag变成了0.0017

应该可以忽略,能量差别和之前的不太大。另外Fe2O3需要DFT+U,建议进行设置。
另外,你的POTCAR用的是LDA(CA-PZ)的,但实际上指定的交换-关联泛函是GGA=PE,建议不要这么做。
woshixuebaxixi 发表于 Post on 2024-3-7 10:22:14
卡开发发 发表于 2024-3-7 10:17
数值误差是会有的,不过我建议你先检查下你那个DFT+U的设置。

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

093f635bacd6f92abe9b05521d34193.png (56.85 KB, 下载次数 Times of downloads: 14)

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

数值误差是会有的,不过我建议你先检查下你那个DFT+U的设置。
fanzhichao001 发表于 Post on 2024-3-7 09:51:49
卡开发发 发表于 2024-3-5 19:17
要。不过你最好看看跑不起来的原因是什么。

好的,谢谢老师。
woshixuebaxixi 发表于 Post on 2024-3-6 14:57:13
本帖最后由 woshixuebaxixi 于 2024-3-6 14:58 编辑

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

INCAR

2.93 KB, 下载次数 Times of downloads: 4

OSZICAR

60.21 KB, 下载次数 Times of downloads: 5

CONTCAR

6.74 KB, 下载次数 Times of downloads: 3

KPOINTS

86 Bytes, 下载次数 Times of downloads: 2

OUTCAR

2.27 MB, 下载次数 Times of downloads: 3

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 22:36 , Processed in 0.245866 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list