24小时热门版块排行榜    

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

lu_wei07

银虫 (小有名气)

[求助] VASP 磁矩 mag跟OUTCAR最后total 不相等??

算单个Pt原子
System = Pt
ENCU = 300
MAGMOM = 3
ISPIN =2
ISMEAR = 0
SIGMA = 0.1
PREC = Normal
LWAVE= F
LREAL= F
LORBIT=12


OUTCAR的输出文件出现这么两行
(1), number of electron       9.9999960 magnetization       2.0000512

(2) ----------------------------------------
  1        0.442   0.000   7.905   8.347
magnetization (x)

# of ion     s       p       d       tot
----------------------------------------
  1        0.442   0.000   0.888   1.330
为什么上面的2 跟下面的1.33不想等呢??按照洪特定则Pt的磁矩是2没错,
那1.33又是什么意思呢??
我试了下其他的原子也是这样不等,就解??
回复此楼

» 猜你喜欢

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

这只是一场游戏
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
lu_wei07: 金币+10, ★★★很有帮助 2013-05-06 13:01:50
franch: 金币+2, 谢谢回帖交流,, 2013-05-06 22:07:36
这个问题是和RWIGS有关,你设置了LORBIT=12分解每个原子spd轨道的磁矩,这必然涉及原子半径问题,或重叠或有间隙导致两个磁矩不同。
3楼2013-05-06 10:36:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-05-06 22:07:20
引用回帖:
5楼: Originally posted by lu_wei07 at 2013-05-06 13:02:55
也就是说2 是在整个空间的积分,而1.33是在RWIGS 中的积分了?? 我没设置所以就默认贋势中的RWIGS???...

我的理解是,只要包含每个原子的部分,就必须有个截止,这个截止就是RWIGS,不然没法判断这个原子的磁矩是多少。我之前用bader分析,得到了不同原子的原子磁矩,那个结果也差不多,也存在误差。
6楼2013-05-06 17:33:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lu_wei07 的主题更新
信息提示
请填处理意见