计算化学公社

标题: 求助：a-Fe2O3（001）面的结构优化问题 [打印本页]

作者
Author: fanzhichao001 时间: 2024-2-28 16:23
标题: 求助：a-Fe2O3（001）面的结构优化问题
求助：如图是对a-Fe2O3的（001）面进行切面、扩胞、建真空层后的结构，对此结构用以下INCAR参数进行结构优化。优化后OSZICAR的最后一步mag=62.7321正常吗？a-Fe2O3是反铁磁体系，总磁矩不应该为0嘛？
(, 下载次数 Times of downloads: 13)

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

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

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

作者
Author: fanzhichao001 时间: 2024-2-29 09:31

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

那老师像我这样的结构是要一个一个原子设置MAGMOM嘛，具体应该如何设置此参数呢？

作者
Author: 卡开发发 时间: 2024-2-29 09:33

fanzhichao001 发表于 2024-2-29 09:31
那老师像我这样的结构是要一个一个原子设置MAGMOM嘛，具体应该如何设置此参数呢？

你看看我答复你的另外一个帖子呢[doge]？

作者
Author: fanzhichao001 时间: 2024-2-29 09:40

卡开发发发表于 2024-2-29 09:33
你看看我答复你的另外一个帖子呢[doge]？

老师是要把自旋向下的换一种原子来跟Fe区分，然后再去设置MAGMOM嘛？

作者
Author: 卡开发发 时间: 2024-2-29 11:16

fanzhichao001 发表于 2024-2-29 09:40
老师是要把自旋向下的换一种原子来跟Fe区分，然后再去设置MAGMOM嘛？

对，建模的时候临时换一下。

作者
Author: fanzhichao001 时间: 2024-2-29 16:59

卡开发发发表于 2024-2-29 11:16
对，建模的时候临时换一下。

好的谢谢老师解答，那老师我去设置MAGMOM的时候一定要一个一个原子对应POSCAR里的坐标去设置嘛？还是说有其他简便一点的设置方式呢？

作者
Author: 卡开发发 时间: 2024-2-29 17:55

fanzhichao001 发表于 2024-2-29 16:59
好的谢谢老师解答，那老师我去设置MAGMOM的时候一定要一个一个原子对应POSCAR里的坐标去设置嘛？还是说有 ...

不需要，前面提到更换原子过后vesta所产生的*.vasp（也就是POSCAR）格式会按照相同的元素排到一起，假定有10个磁矩为+2的A原子、1个没磁性的B原子和14个磁矩为-4的C原子，这样MAGMOM可以写成10*2 0 14*-4。

作者
Author: woshixuebaxixi 时间: 2024-3-1 16:06
请问卡总 Fe是不是应该写成一般自旋向上用+5 一半自旋向下用-5呀？

作者
Author: fanzhichao001 时间: 2024-3-1 17:27

卡开发发发表于 2024-2-29 17:55
不需要，前面提到更换原子过后vesta所产生的*.vasp（也就是POSCAR）格式会按照相同的元素排到一起，假定 ...

好的，谢谢老师耐心解答，我去试一下

作者
Author: 卡开发发 时间: 2024-3-1 20:13

woshixuebaxixi 发表于 2024-3-1 16:06
请问卡总 Fe是不是应该写成一般自旋向上用+5 一半自旋向下用-5呀？

对，并且需要注意顺序。前面提到替换的原因是因为这样某些转换程序在产生POSCAR格式的时候会将同一类原子划分到一组，这样写MAGMOM会简单很多。

作者
Author: woshixuebaxixi 时间: 2024-3-1 20:58
感谢卡总！我也是算了α-Fe2O3（110）晶面结果mag= -0.0274 这样合理吗？

作者
Author: 卡开发发 时间: 2024-3-2 09:16

woshixuebaxixi 发表于 2024-3-1 20:58
感谢卡总！我也是算了α-Fe2O3（110）晶面结果mag= -0.0274 这样合理吗？

一般应该是0，我不确定你的体系是否符合计量比（包括缺陷和其他物质吸附等问题），以及自洽场是否真的达到收敛标准。

作者
Author: fanzhichao001 时间: 2024-3-5 16:30
本帖最后由 fanzhichao001 于 2024-3-5 16:31 编辑

卡开发发发表于 2024-3-2 09:16
一般应该是0，我不确定你的体系是否符合计量比（包括缺陷和其他物质吸附等问题），以及自洽场是否真的达 ...

老师好，我设置MAGMOM优化后mag=0了，但是我进行下一步引入C5H8在表面吸附计算时，用同样的INCAR却跑不起来是什么原因呢？是还需要设置C和O的磁矩=0嘛
(, 下载次数 Times of downloads: 9) (, 下载次数 Times of downloads: 13)

作者
Author: 卡开发发 时间: 2024-3-5 19:17

fanzhichao001 发表于 2024-3-5 16:30
老师好，我设置MAGMOM优化后mag=0了，但是我进行下一步引入C5H8在表面吸附计算时，用同样的INCAR却跑不 ...

要。不过你最好看看跑不起来的原因是什么。

作者
Author: woshixuebaxixi 时间: 2024-3-6 14:57
本帖最后由 woshixuebaxixi 于 2024-3-6 14:58 编辑

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

作者
Author: fanzhichao001 时间: 2024-3-7 09:51

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

好的，谢谢老师。

作者
Author: 卡开发发 时间: 2024-3-7 10:17

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

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

作者
Author: woshixuebaxixi 时间: 2024-3-7 10:22

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

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

作者
Author: 卡开发发 时间: 2024-3-7 12:13

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

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

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

作者
Author: 卡开发发 时间: 2024-3-7 15:01

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

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

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

作者
Author: 卡开发发 时间: 2024-4-3 15:36

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

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

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

作者
Author: 卡开发发 时间: 2024-4-4 11:02

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

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

作者
Author: 123wo 时间: 2024-4-12 10:39

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

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

作者
Author: 卡开发发 时间: 2024-4-13 02:14

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

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

作者
Author: 123wo 时间: 2024-4-13 10:11

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

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

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