应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】用vasp计算磁矩的问题
141/2返回列表12下一页
查看: 4920  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

congmingbing

木虫 (正式写手)

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

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-06-28 20:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bingmou

金虫 (著名写手)
congmingbing(金币+1):感谢您的回复 2010-06-29 15:19:41
我只能说，磁性计算是很复杂，VASP是很有限的。
我就最后一个问题说一下吧。你在给ＭＡＧＭＯＭ有没有想过这是什么意思呢？MAGMON=71*0 1,这样的式子表示：那个磁性原子的初始磁矩方向是向着x方向（内建坐标），并且为1。但是实际上，可能这个方向并不是原子的易磁化方向，这个原子的磁矩并不向这个方向。
同样的，MAGMOM=70*0 1 -1，这样强行要求这两个原子的磁矩朝着x方向反向排列，体系的能量很高。所以VASP会让他们朝别的方向，例如一个（111），一个（1-1-1),而不是你想要的反向排列。除非采用复杂的方向磁化计算，进行精心的“人工“设置，不然你对他没有办法。
在磁性面前，我们只能杯具
赞一下(2人) 回复此楼
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中不一定能完全重现的
赞一下(1人) 回复此楼
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，看看
赞一下(1人) 回复此楼
10楼2010-06-29 16:23:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 congmingbing 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:53 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:33 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ph7bxbho8f 2026-07-09 3/150 2026-07-10 05:53 by j4frgf6zcw
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[基金申请] 祈祷青基必中 50+3 hadengry 2026-07-09 9/450 2026-07-09 17:37 by 文杰猪猪
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 关于如何从代码看上不上会 +11 且听虎啸 2026-07-09 14/700 2026-07-09 15:27 by libeny
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评？应用基础研究和基础研究，希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[硕博家园] 大龄残疾硕士的一点执念 +14 3776737374 2026-07-06 22/1100 2026-07-09 06:07 by lovechizi
[基金申请] 生命口会评 +4 wuke100666 2026-07-07 4/200 2026-07-08 16:38 by 逍遥不够自在
[教师之家] 网络举报学术不端成为了男女之间情杀的手段了，真荒唐。。。 +3 zju2000 2026-07-04 3/150 2026-07-07 21:00 by hyx552200
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联，生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
[论文投稿] 刑法学论文投稿求助 +3 高源678 2026-07-04 5/250 2026-07-06 12:20 by 平天尺
[论文投稿] Journal of Environmental Chemical Engineering +3 雪花天涯 2026-07-05 3/150 2026-07-06 09:00 by bobvan
[基金申请] 最后一年，祈求好运 +3 pyw1992 2026-07-06 3/150 2026-07-06 08:53 by 苏_frp
信息提示
关闭
请填处理意见
关闭