计算化学公社

标题: 求助：对反铁磁体系加电场 [打印本页]

作者
Author: diffcwmol 时间: 2023-7-11 09:18
标题: 求助：对反铁磁体系加电场
各位同仁，向您请教在VASP计算反铁磁体系中出现的一个问题。

我通过设置初始磁矩得到了体系的反铁磁耦合状态，现在想加电场看体系的\alpha和\beta能带迁移。
但是加了电场之后发现体系的原子的磁矩都变成0了，不知出了啥问题。

加电场使用的是：
LDIPOL =.TRUE.
IDIPOL = 2
DIPOL = 0.5 0.5 0.5
EFIELD = 0.1

使用（或者不使用）未加电场时收敛的波函数和电荷密度，体系都收敛到原子磁矩约为0的状态。
试了加NUPDOWN = 0； MAGMOM = 1.0 -1.0 或者不加，对结果没有影响。
收敛后能带变了，但是没有规律。

作者
Author: 阿司匹林泡腾片 时间: 2023-7-11 16:41
本帖最后由阿司匹林泡腾片于 2023-7-11 17:32 编辑

我的经验：1.NUPDOWN = 0（这个反铁磁体系最好加上，虽然官网的反铁磁计算好像也没加，但是加的话你要分析一下总磁矩是不是为0，例如：如果出现了体系被拿走一个H+或者得到个OH-的就不是0了）
2.MAGMOM对于反铁磁我认为必须加（愚见），例如有四个原子两个自旋向上两个自旋向下，就得：MAGMOM = -1 1 -1 1 类似这种方式给他们进行排布，主要是指定磁矩方向了，尽量别1，官网上也说尽量设置的大一点，但是正负的设置这个要按照这个体系的磁矩分布情况设置，每个体系都不一样。
一点愚见仅供参考，有不对的地方敬请指点一二

作者
Author: diffcwmol 时间: 2023-7-11 18:48

阿司匹林泡腾片发表于 2023-7-11 16:41
我的经验：1.NUPDOWN = 0（这个反铁磁体系最好加上，虽然官网的反铁磁计算好像也没加，但是加的话你要分析 ...

非常感谢您的耐心回复和建议。
请问您对加电场后反铁磁耦合消失（理论上不应该消失）有没有好的建议？

作者
Author: 阿司匹林泡腾片 时间: 2023-7-11 21:13

diffcwmol 发表于 2023-7-11 18:48
非常感谢您的耐心回复和建议。
请问您对加电场后反铁磁耦合消失（理论上不应该消失）有没有好的建议？

我是做催化的，没有研究过这方面的内容

，期待论坛里大佬们的回复

作者
Author: diffcwmol 时间: 2023-7-14 08:35
问题已经解决，需要破坏体系的对称性，加关键词ISYM=0。

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