24小时热门版块排行榜    

查看: 5065  |  回复: 13

congmingbing

木虫 (正式写手)

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

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

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

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的回帖

bingmou

金虫 (著名写手)


gavinliu7390(金币+1):谢谢交流! 2010-06-29 15:42:42
引用回帖:
Originally posted by congmingbing at 2010-06-29 15:23:23:


你说的这个问题我也考虑过,因为不确定vasp默认的磁矩方向(不过好像记得在什么地方看到过,说是沿z方向的),所以我对晶格做了不同的定义,让磁矩沿x,y,z不同的方向,结果是一样的,能量也基本完全相同。而实 ...

如果你的晶胞具有O群对称性,那么沿xyz能量肯定都一样的。
实验缺的方向在VASP中不一定能完全重现的
5楼2010-06-29 15:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

引用回帖:
Originally posted by congmingbing at 2010-06-29 15:29:44:
其他输入都不变,只是在INCAR中加入LORBIT参数,得到的结果怎么会和不加这个差别很大呢?

这一个我觉得很奇怪。LORBIT只是在输出DOS的时候乎变化,不应该影响体系的能量的阿。
是否可能是RWIGS的影响(这个不太可能)
6楼2010-06-29 15:43:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

congmingbing

木虫 (正式写手)

[quote]Originally posted by bingmou at 2010-06-29 15:40:17:

加上LORBIT后总能是没有变化的,不过OSZICAR中的总磁矩完全变了。
7楼2010-06-29 15:47:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

congmingbing(金币+3):非常感谢! 2010-06-29 16:09:39
引用回帖:
Originally posted by congmingbing at 2010-06-29 15:47:22:
[quote]Originally posted by bingmou at 2010-06-29 15:40:17:

加上LORBIT后总能是没有变化的,不过OSZICAR中的总磁矩完全变了。

这个有可能是RWIGS的问题
8楼2010-06-29 16:05:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

congmingbing

木虫 (正式写手)

引用回帖:
Originally posted by congmingbing at 2010-06-29 15:47:22:
[quote]Originally posted by bingmou at 2010-06-29 15:40:17:

加上LORBIT后总能是没有变化的,不过OSZICAR中的总磁矩完全变了。

这个一般要怎么设呢,我都是用的默认值,没有设置过这个参数?是在自恰时同时加上lorbit和RWIGS吗?
另外,默认的磁矩方向是沿x吗,手册上没有说呢?我一直都认为沿z方向的。
9楼2010-06-29 16:11:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)


xiaohunhun(金币+1):xiexie 2010-06-29 23:44:49
引用回帖:
Originally posted by congmingbing at 2010-06-29 16:11:38:


这个一般要怎么设呢,我都是用的默认值,没有设置过这个参数?是在自恰时同时加上lorbit和RWIGS吗?
另外,默认的磁矩方向是沿x吗,手册上没有说呢?我一直都认为沿z方向的。

我问过一个老师,据说是朝X。但是我觉得沿Z方向也可以,无所谓。你可以问问开发组。
RWIGS是可以不设的,会采用势里面的默认值。但是我还是不能想出RWIGS对磁矩的影响。你可以看看OUTCAR里面的RWIGS,看看是不是有可能出现重叠了。而且在OUTCAR里面最后的部分有给出每个原子的local magnetization,看看
10楼2010-06-29 16:23:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 congmingbing 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +53 医学老男孩 2026-08-13 118/5900 2026-08-19 18:14 by 医学老男孩
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 3/150 2026-08-19 18:04 by oCMhZAsgNFoO
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 17:42 by oCMhZAsgNFoO
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +4 瞬息宇宙 2026-08-15 4/200 2026-08-19 16:23 by allmydream
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见