24小时热门版块排行榜    

查看: 4258  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 我感觉当老师好累啊。。 +16 fairy1122 2024-11-15 17/850 2024-11-16 11:59 by ma0526
[硕博家园] 硕博人来此 +3 ColonMo 2024-11-16 4/200 2024-11-16 11:26 by yanghu01
[硕博家园] 和导师撕破脸了 +18 圆满123456 2024-11-13 22/1100 2024-11-16 09:49 by zyhou
[第一性原理] 有没有要一起买VASP版权的,可以6个人 5+3 变成小神 2024-11-14 8/400 2024-11-16 09:31 by 变成小神
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+4 猪小耍 2024-11-13 15/750 2024-11-16 09:24 by guxue
[论文投稿] 国内期刊审稿人数量 +4 新时代核动力驴 2024-11-13 6/300 2024-11-16 08:24 by pityzhou
[教师之家] 这种人痛苦吗 +12 2671 2024-11-12 14/700 2024-11-15 22:15 by 仇恨毁灭财富
[论文投稿] chemical science和advanced science哪个好一点 +6 yly150 2024-11-12 6/300 2024-11-15 21:06 by hcj9909053
[论文投稿] 找不到审稿人 +6 吉祥如意141913 2024-11-13 6/300 2024-11-15 16:48 by LEyblyt
[基金申请] 博后面上太难了 +14 不安定因素 2024-11-14 17/850 2024-11-15 16:05 by lupbao
[基金申请] 去年七月底入站的还能申请下一批吗? +4 brightwo 2024-11-14 4/200 2024-11-15 14:18 by gazi1111
[基金申请] 广东省自然科学基金-面上项目消息有了么? +4 pachang 2024-11-14 6/300 2024-11-15 09:39 by zhoulei7170
[基金申请] 博后基金分组排名 +7 攻城2025 2024-11-14 7/350 2024-11-14 21:19 by 实验小白ha
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[基金申请] 浙江省青基怎么看A还是B? +8 ToallGoof 2024-11-11 14/700 2024-11-13 09:11 by Praybest
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[论文投稿] 投稿求助 40+3 liukun11 2024-11-12 4/200 2024-11-12 19:13 by liukun11
信息提示
请填处理意见