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

fanzhichao001

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

卡开发发

不正常。属于不设置MAGMOM的一个后果。

fanzhichao001

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

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

卡开发发

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

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

fanzhichao001

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

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

卡开发发

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

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

fanzhichao001

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

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

卡开发发

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。

woshixuebaxixi

请问卡总 Fe是不是应该写成一般自旋向上用+5 一半自旋向下用-5呀？

fanzhichao001

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

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

卡开发发

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

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

woshixuebaxixi

感谢卡总！我也是算了α-Fe2O3（110）晶面   结果mag=    -0.0274 这样合理吗？

卡开发发

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

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

fanzhichao001

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

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

卡开发发

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

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