计算化学公社

标题: VASP计算反铁磁体系问题 [打印本页]

作者
Author:
linuxprobe    时间: 2017-4-7 15:42
标题: VASP计算反铁磁体系问题
在用VASP 计算反铁磁体系(含8个铁原子)时,对于MAGMOM该怎么设置磁矩的方向呢?还有对于这样的体系优化时需要加上SPIN和MAGMOM参数吗?谢谢

作者
Author:
jiangning198511    时间: 2017-4-7 16:26
本帖最后由 jiangning198511 于 2017-4-7 16:28 编辑

NiO的例子:

SYSTEM = NiO
   
ISTART = 0
   
ISPIN = 2
MAGMOM = 2.0 -2.0 2*0
     
ENMAX = 250.0
EDIFF = 1E-3
   
ISMEAR = -5
   
AMIX = 0.2
BMIX = 0.00001
AMIX_MAG = 0.8
BMIX_MAG = 0.00001
   
LORBIT = 11

POSCAR:
AFM  NiO 4.17
1.0 0.5 0.5
0.5 1.0 0.5
0.5 0.5 1.0
2 2
Cartesian
0.0 0.0 0.0
1.0 1.0 1.0
0.5 0.5 0.5
1.5 1.5 1.5


作者
Author:
linuxprobe    时间: 2017-4-7 21:33
jiangning198511 发表于 2017-4-7 16:26
NiO的例子:

SYSTEM = NiO

谢谢回复!对于二个原子的MAGMOM可以参考这样的设置,但是对于多个原子体系(内含8个铁),这个MAGMOM中的磁矩方向该如何定呢?
作者
Author:
jiangning198511    时间: 2017-4-8 16:06
linuxprobe 发表于 2017-4-7 21:33
谢谢回复!对于二个原子的MAGMOM可以参考这样的设置,但是对于多个原子体系(内含8个铁),这个MAGMOM中 ...

和例子一样,根据你的POSCAR看看哪个向上,哪个向下,例如
5 5 5 5 -5 -5 -5 -5
前四个向上,后四个向下,总自旋为0
这个值你可以设置大一些,程序会自动调节的
作者
Author:
linuxprobe    时间: 2017-4-10 17:28
jiangning198511 发表于 2017-4-8 16:06
和例子一样,根据你的POSCAR看看哪个向上,哪个向下,例如
5 5 5 5 -5 -5 -5 -5
前四个向上,后四个向 ...

非常感谢老师!那如何看POSCAR来确定方向是向上还是下呢?比如对下面的POSCAR来说:
POSCAR
1.0
        8.7869997025         0.0000000000         0.0000000000
        0.0000000000        11.2250003815         0.0000000000
        0.0000000000         0.0000000000         5.2880001068
   Ca    O   Fe
    4   12    8
Direct
     0.500000000         0.184499994         0.750000000
     0.500000000         0.815500021         0.250000000
     0.000000000         0.684499979         0.750000000
     0.000000000         0.315500021         0.250000000
     0.500000000         0.615599990         0.750000000
     0.500000000         0.384400010         0.250000000
     0.000000000         0.115599990         0.750000000
     0.000000000         0.884400010         0.250000000
     0.207399994         0.382400006         0.750000000
     0.792600036         0.617599964         0.250000000
     0.792600036         0.382400006         0.750000000
     0.207399994         0.617599964         0.250000000
     0.707399964         0.882400036         0.750000000
     0.292600006         0.117599994         0.250000000
     0.292600006         0.882400036         0.750000000
     0.707399964         0.117599994         0.250000000
     0.346100003         0.500000000         0.500000000
     0.653900027         0.500000000         0.500000000
     0.653900027         0.500000000         0.000000000
     0.346100003         0.500000000         0.000000000
     0.846099973         0.000000000         0.500000000
     0.153899997         0.000000000         0.500000000
     0.153899997         0.000000000         0.000000000
     0.846099973         0.000000000         0.000000000
8个Fe的方向该如何选取呢?谢谢您
作者
Author:
jiangning198511    时间: 2017-4-10 19:36
linuxprobe 发表于 2017-4-10 17:28
非常感谢老师!那如何看POSCAR来确定方向是向上还是下呢?比如对下面的POSCAR来说:
POSCAR
1.0

4*0.6 12*0.6 5 5 5 5 -5 -5 -5 -5

这是一种,你可以多试试几种排列,看看哪种能量最低

作者
Author:
linuxprobe    时间: 2017-4-10 20:53
jiangning198511 发表于 2017-4-10 19:36
4*0.6 12*0.6 5 5 5 5 -5 -5 -5 -5

这是一种,你可以多试试几种排列,看看哪种能量最低

好的,谢谢老师解答!
作者
Author:
guofengxiao    时间: 2017-4-11 10:05
linuxprobe 发表于 2017-4-10 17:28
非常感谢老师!那如何看POSCAR来确定方向是向上还是下呢?比如对下面的POSCAR来说:
POSCAR
1.0

......你就不能找个可视化软件比如MS或者VESTA什么的,根据可视化出来的图和坐标一一对应,来判断吗。。。
作者
Author:
error_323    时间: 2022-3-8 14:11
jiangning198511 发表于 2017-4-7 16:26
NiO的例子:

SYSTEM = NiO

我对类似的FeO体系进行结构优化,优化后晶格常数不一致,是为什么呢
作者
Author:
jiangning198511    时间: 2022-3-8 14:44
error_323 发表于 2022-3-8 14:11
我对类似的FeO体系进行结构优化,优化后晶格常数不一致,是为什么呢

不同自旋态导致不同的结构 这个可以理解

作者
Author:
error_323    时间: 2022-3-15 16:22
jiangning198511 发表于 2022-3-8 14:44
不同自旋态导致不同的结构 这个可以理解

您好!我想请教一下,对FeO反铁磁体系结构优化,应该优化原胞还是单胞呢,我优化之后空间群就变了。正常来说优化晶体结构,空间群是不会变化的吧
作者
Author:
jiangning198511    时间: 2022-3-17 11:06
error_323 发表于 2022-3-15 16:22
您好!我想请教一下,对FeO反铁磁体系结构优化,应该优化原胞还是单胞呢,我优化之后空间群就变了。正常 ...

根据你的计算而定,比如你要计算反铁磁结果,如果结果中的Fe原子能够满足反铁磁的排布,就直接用这个结构,如果不能满足 就扩胞,关于对称性 一般计算会有偏移,可能会破坏之前的对称性
作者
Author:
error_323    时间: 2022-3-17 14:56
jiangning198511 发表于 2022-3-17 11:06
根据你的计算而定,比如你要计算反铁磁结果,如果结果中的Fe原子能够满足反铁磁的排布,就直接用这个结构 ...

不同的磁矩排列会计算出不同的能量,是不是要取能量最低的那种磁矩排列优化结构呢?因为不同的磁矩排列,优化后的结构不一样,差别还挺大(如:原始结构晶格常数4.33 4.33 4.33,一种磁矩优化结果4.44 4.44 4.45 能量-50;另一种4.45  4.45  4.20 能量-53)
作者
Author:
jiangning198511    时间: 2022-3-17 16:33
error_323 发表于 2022-3-17 14:56
不同的磁矩排列会计算出不同的能量,是不是要取能量最低的那种磁矩排列优化结构呢?因为不同的磁矩排列, ...

取能量最低
作者
Author:
wrizad    时间: 2023-5-9 09:44
我想请教一下 超孢一般是往2*2*2超孢后 然后进行结构优化么

作者
Author:
wangyuxiao    时间: 2023-5-9 14:18
guofengxiao 发表于 2017-4-11 10:05
**** 作者被禁止或删除 内容自动屏蔽 ****

你好,用MS或VESTA怎么显示出是向上还是向下




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3