24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4241  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

congmingbing

木虫 (正式写手)

[交流] 【求助】用vasp计算磁矩的问题已有2人参与

用vasp计算磁性材料时,自恰计算的INCAR中加上ISPIN和MAGMOM参数,总磁矩除了在OSZICAR中,在OUTCAR中能给出吗,我没有找到呢。
对反铁磁的计算,总磁矩在OSZICAR中是0,但是如果在自恰计算时INCAR中加上LORBIT=11后,OSZICAR中的总磁矩就不是0了,不知道这是什么原因,同时OUTCAR中给出两个磁性原子的磁矩都是正值,这和实际情况相悖,这个问题请高手们指点。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyniu

金虫 (初入文坛)


zzy870720z(金币+1):鼓励交流 2010-12-28 23:52:31
VASP 说明书第55页有这么一段话:With the MAGMOM line specified above, VASP should converge to the proper groundstate. In this example, the total net magnetisation is matter of factly zero, but it is possible to determine the local magnetic moments by using the
RWIGS or LORBIT tags (see sections 6.33 6.32). 加上LORBIT=11可能算的就是局域了吧!
11楼2010-12-28 21:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

bingmou

金虫 (著名写手)

congmingbing(金币+1):感谢您的回复 2010-06-29 15:19:41
我只能说,磁性计算是很复杂,VASP是很有限的。
我就最后一个问题说一下吧。你在给MAGMOM有没有想过这是什么意思呢?MAGMON=71*0 1,这样的式子表示:那个磁性原子的初始磁矩方向是向着x方向(内建坐标),并且为1。但是实际上,可能这个方向并不是原子的易磁化方向,这个原子的磁矩并不向这个方向。
同样的,MAGMOM=70*0 1 -1,这样强行要求这两个原子的磁矩朝着x方向反向排列,体系的能量很高。所以VASP会让他们朝别的方向,例如一个(111),一个(1-1-1),而不是你想要的反向排列。除非采用复杂的方向磁化计算,进行精心的“人工“设置,不然你对他没有办法。
在磁性面前,我们只能杯具
2楼2010-06-29 15:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

congmingbing

木虫 (正式写手)

引用回帖:
Originally posted by congmingbing at 2010-06-28 20:49:48:
用vasp计算磁性材料时,自恰计算的INCAR中加上ISPIN和MAGMOM参数,总磁矩除了在OSZICAR中,在OUTCAR中能给出吗,我没有找到呢。
对反铁磁的计算,总磁矩在OSZICAR中是0,但是如果在自恰计算时INCAR中加上LORBIT= ...

你说的这个问题我也考虑过,因为不确定vasp默认的磁矩方向(不过好像记得在什么地方看到过,说是沿z方向的),所以我对晶格做了不同的定义,让磁矩沿x,y,z不同的方向,结果是一样的,能量也基本完全相同。而实验确定的磁矩确实是沿某一晶向的,如(1 0 0),(0 1 0)或(0 0 1)。
3楼2010-06-29 15:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

congmingbing

木虫 (正式写手)

其他输入都不变,只是在INCAR中加入LORBIT参数,得到的结果怎么会和不加这个差别很大呢?
4楼2010-06-29 15:29:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见