24小时热门版块排行榜    

查看: 4502  |  回复: 9

lvyan307

金虫 (小有名气)


[交流] vasp算磁性时关于MAGMOM设置的一个奇怪问题

最近在学材料的磁性计算,假设元胞内有6个原子,按照手册设置ISPIN = 2,MAGMOM = 6*1,那么此时算的是体系的铁磁性。查文献得知这种材料反铁磁时能量更低,于是做反铁磁计算,我重新调整了POSCAR和POTCAR中原子顺序,把自旋向上的3个原子放在前面,自旋向下的原子放后面,相应的设置ISPIN = 2, MAGMOM = 3*1 3*-1。计算结果很奇怪,此时算得材料的总磁矩和能量与前面算铁磁性时得到的是一样的。
仔细分析OUTCAR, 发现MAGMOM如上面两种设置后体系的初始磁矩magnetization都是等于元胞内原子数6,也就是说计算过程中没有读出MAGMOM中的设置。
我又算了一次反铁磁性,这时设置ISPIN = 2, MAGMOM = 1 1 1 -1 -1 -1,奇怪的是这个时候得到了正确的结果,即OUTCAR中初始磁矩为0,弛豫后体系的总磁矩也为0,能量也比铁磁性时更低。
不知有谁能解释下为什么会这样?
谢谢!
回复此楼

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

strive123123

禁虫 (文坛精英)

本帖内容被屏蔽

2楼2012-07-18 08:40:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
fzx2008: 金币+1, 谢谢指教! 2012-07-18 10:04:18
lvyan307: 金币+2 2012-07-19 14:40:16
改成1*3 -1*3
3楼2012-07-18 09:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
sunyang1988: 金币+1, 不知道为什么了,我没出过问题啊 2012-07-18 15:50:11
lvyan307: 金币+2 2012-07-19 14:40:23
奇怪,我的这样写没有问题
MAGMOM = 48*0  4*1.0  4*-1.0
4楼2012-07-18 10:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvyan307

金虫 (小有名气)


引用回帖:
3楼: Originally posted by sunyang1988 at 2012-07-18 09:35:43
改成1*3 -1*3

我这样改过,不行啊!
5楼2012-07-18 11:38:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvyan307

金虫 (小有名气)


引用回帖:
4楼: Originally posted by fzx2008 at 2012-07-18 10:15:03
奇怪,我的这样写没有问题
MAGMOM = 48*0  4*1.0  4*-1.0

我觉得理论上是可以这样写的,但可能不同的系统装vasp编译时有些没编译好。当然,我没编译过vasp, 也有可能是其它原因。
我发现本论坛其他几个帖子说做反铁磁性不成功应该与我的问题是一样的。
6楼2012-07-18 11:42:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ip0371

木虫 (正式写手)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by lvyan307 at 2012-07-18 11:42:50
我觉得理论上是可以这样写的,但可能不同的系统装vasp编译时有些没编译好。当然,我没编译过vasp, 也有可能是其它原因。
我发现本论坛其他几个帖子说做反铁磁性不成功应该与我的问题是一样的。...

问题解决了吗?楼主
7楼2014-06-23 22:14:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lvyan307

金虫 (小有名气)


引用回帖:
7楼: Originally posted by ip0371 at 2014-06-23 22:14:05
问题解决了吗?楼主...

问题解决了,不过不知道为什么,现在我用一楼提到的那两种MAGMOM的设置方式都可以得到正确的计算结果。
8楼2014-06-24 12:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1499526749

新虫 (初入文坛)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主,请问怎么知道哪个原子是自旋向上的,哪个原子是自旋向下的?急求!谢谢
9楼2016-09-09 11:06:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1499526749

新虫 (初入文坛)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by 1499526749 at 2016-09-09 11:06:51
楼主,请问怎么知道哪个原子是自旋向上的,哪个原子是自旋向下的?急求!谢谢

我的结构都是12个同种原子,
10楼2016-09-09 11:10:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lvyan307 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 PYZmcHr5kSaq 2026-11-05 8/400 2026-11-06 06:16 by qL1cU46YIdK0
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 05:54 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 05:44 by qL1cU46YIdK0
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[基金申请] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 05:33 by qL1cU46YIdK0
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 05:21 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 05:15 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-06 04:58 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 01:29 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 01:20 by qL1cU46YIdK0
信息提示
请填处理意见