24小时热门版块排行榜    

查看: 5621  |  回复: 14

Prettyswan

木虫 (小有名气)

[交流] 【求助成功】如何用VASP计算反铁磁 已有7人参与

第一次使用VASP计算反铁磁,不知道如何入手。请教。
主要想知道:
反铁磁计算的一般步骤;
反铁磁如何建模,POSCAR如何得到;
INCAR文件中关于磁矩如何得到、如何设置;
等一系列问题。
请做过反铁磁计算的踊跃发言,越详细越好。如果有可能给出点参考资料更好。
谢谢。

[ Last edited by Prettyswan on 2010-7-5 at 09:19 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 第一性原理相关文档 vasp计算资料 VASP软件从入门到高深
vasp 第一性原理杂 VASP

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

valenhou001

至尊木虫 (职业作家)

★ ★
lvjian8596(金币+2):感谢侯老师,希望以后多多交流! 2010-06-22 08:43:06
反铁磁:需要根据你所要计算的磁构型构造超原胞。并在INCAR中将该自旋是相反的两个原子的初始磁矩设置为大小相等但符号相反的数。

比如磁性原子为A,非磁性态时的结构为简单立方,里面只含一个原子A。那么考虑它的反铁磁,需将它扩展成1x1x2, 1x1x2, 1x2x1或者改为FCC基矢样子:
原来的晶格常数
  0.00000000000000    1.00000000000000    1.00000000000000
    1.00000000000000    0.00000000000000    1.00000000000000
    1.00000000000000    1.00000000000000    0.00000000000000
原子坐标是(0,0,0)和(0.5, 0.5, 0.5)。
ISPIN = 2
MAGMOM =  2.0   -2.0
3楼2010-06-22 08:39:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

congmingbing

木虫 (正式写手)


gavinliu7390(金币+1):谢谢交流! 2010-06-29 15:43:26
引用回帖:
Originally posted by xiaojie7783 at 2010-06-22 16:49:26:

我把体系中可能出现的磁构型就是磁矩的政府情况都计算了一遍,结果计算出来的总能却是一样的,这样的情况要怎么处理呢,是我设置的问题吗?

我也遇到了这样的问题。所有的磁构型能量都是一样的,而且和顺磁,也就是非磁计算的总能也是相同的。加上lorbit=11 参数后,在outcar中看到的磁构型不是我在INCAR中所设置的磁构型,比如我在INCAR中设置有磁性的两个原子的磁矩是1和-1,最后给出的总磁矩反而不是0了,这是什么原因导致的呢?

[ Last edited by congmingbing on 2010-6-25 at 10:19 ]
9楼2010-06-25 10:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

锐利的碎片

木虫 (正式写手)

star watcher


hedaors(金币+1):谢谢回复 2010-06-22 07:58:53
机器好wien2k或者elk这些算磁性比较好,后者的例子里有几个很详细的例子。
2楼2010-06-21 23:21:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-06-22 08:39:08:
反铁磁:需要根据你所要计算的磁构型构造超原胞。并在INCAR中将该自旋是相反的两个原子的初始磁矩设置为大小相等但符号相反的数。

比如磁性原子为A,非磁性态时的结构为简单立方,里面只含一个原子A。那么考虑 ...

想请问一下,我现在有一个晶胞但是不知道怎么确定其中原子的自旋方向,这个问题要怎么解决呢?非常谢谢您
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
4楼2010-06-22 10:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


zzy870720z(金币+1):感谢侯老师指点 2010-06-22 19:40:21
如果都是共线的计算,很多程序都是将磁轴的方向取为z轴的。那给定原子的磁矩时,也就是默认的沿着z轴的。磁矩正的值,表示沿着z,负值就是沿着-z的方向。至于就是哪种磁构型稳定,要尽可能的考虑多的不同构型,然后比较总能。常见的磁构型,在一些将铁磁学的书都可以找到的。比如郭贻诚先生的<铁磁学>,姜寿亭,李卫 的<凝聚态磁性物理>.

[ Last edited by valenhou001 on 2010-6-22 at 13:29 ]
5楼2010-06-22 13:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-06-22 13:28:12:
如果都是共线的计算,很多程序都是将磁轴的方向取为z轴的。那给定原子的磁矩时,也就是默认的沿着z轴的。磁矩正的值,表示沿着z,负值就是沿着-z的方向。至于就是哪种磁构型稳定,要尽可能的考虑多的不同构型,然 ...

我把体系中可能出现的磁构型就是磁矩的政府情况都计算了一遍,结果计算出来的总能却是一样的,这样的情况要怎么处理呢,是我设置的问题吗?
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
6楼2010-06-22 16:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


zzy870720z(金币+1):辛苦了,谢谢交流 2010-06-22 19:40:32
有可能是设置的问题。初始磁矩的绝对值不要太小。对于非磁性原子的初始磁矩可以设置为0.

将每个构型的计算结果仔细检查一下,看每个原子的局域磁矩的分布情况,也就是所得到磁构型是不是你所设置的磁构型。

可以的话,将你的INCAR, OUTCAR, POSCAR发送我签名档里面留的信箱中,帮你瞧瞧。

[ Last edited by valenhou001 on 2010-6-22 at 19:16 ]
7楼2010-06-22 19:14:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-06-22 08:39:08:
反铁磁:需要根据你所要计算的磁构型构造超原胞。并在INCAR中将该自旋是相反的两个原子的初始磁矩设置为大小相等但符号相反的数。

比如磁性原子为A,非磁性态时的结构为简单立方,里面只含一个原子A。那么考虑 ...

1、大小为a×b×c超胞的建立,如何选择a,b,c的值?遵循一个怎么样的原则?
2、初始磁矩MAGMOM的设置中,对于反铁磁来说自旋反平行,MAGMON数值的选择是怎么定的(就是说,MAGMOM是选1还是选2还是选其他的数有没有什么规律)?
3、你说到“如果都是共线的计算,很多程序都是将磁轴的方向取为z轴的。那给定原子的磁矩时,也就是默认的沿着z轴的。”;而我计算的CaFe2As2中“Fe moments are oriented along the orthorhombic a axis ”(a轴沿着x方向)。那么,如何解决这个问题呢?

[ Last edited by Prettyswan on 2010-6-24 at 22:36 ]
8楼2010-06-24 19:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
zzy870720z(金币+2):谢谢侯老师指点 2010-06-26 08:34:43
Prettyswan(金币+15):谢谢耐心讲解。 2010-06-26 21:33:23
引用回帖:
Originally posted by Prettyswan at 2010-06-24 19:01:17:

1、大小为a×b×c超胞的建立,如何选择a,b,c的值?遵循一个怎么样的原则?
2、初始磁矩MAGMOM的设置中,对于反铁磁来说自旋反平行,MAGMON数值的选择是怎么定的(就是说,MAGMOM是选1还是选2还是选其他的数有 ...

1. 主要是考虑磁性原子之间的磁相互作用的情况:比如参考这里
http://folk.uio.no/ravi/activity/ordering/spinordering.html
http://folk.uio.no/ravi/activity/ordering/spin.gif
考虑层间反铁磁而面内铁磁、层间铁磁而面内反铁磁、层间反铁磁而面内也反铁磁等这些来构造超原胞。

2. MAGMOM设置的是的初始磁距,一般设置的大些,默认的是1.0 uB,往往是不够的,需达到2-3倍,或者根据原子的价电子组态估计出可能的最大磁极化(也就是价电子排布出现的孤对电子最多的情况)。计算会自洽迭代收敛体系处于某个稳定态时的磁矩。
10楼2010-06-25 17:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
信息提示
请填处理意见