24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1007  |  回复: 6

yangbs2009

木虫 (小有名气)

[求助] vasp间隙磁矩已有1人参与

我在用vasp计算一个层状的体系,当计算一个金属空位时有如下的问题:
1、在outcar中显示磁矩是1,但OSZICAR中显示的是3;两者相差的原因可能是间隙磁矩;而且间隙磁矩是原子总磁矩的两倍
2、在计算这个体系的金属原子被替代时也遇到了间隙磁矩是原子总磁矩的两倍这种情况
3、计算的体系共价键更强一些,即两个相邻原子之间的电子应该是反平行排列,既然是反平行排列,间隙处怎么又会有磁矩呢

以上的间隙磁矩合理么,这个结果又如何来解释呢,
vasp中有magmom参数对原子的初始磁矩进行设置,有没有对间隙磁矩设置的参数呢?
回复此楼

» 猜你喜欢

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

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

golddoushi

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
yangbs2009: 金币+20, ★★★很有帮助 2014-02-18 14:07:05
无缺陷无替换时是共价键,不代表你弄过之后还是共价键。这个要看具体情况。
OUTCAR里只有wigner半径内的磁矩,通常都是不准的。一般bader方法算的极化电荷就能跟OSZICAR对上

间隙诱发/局域/磁矩,虽然大家都这么说,个人觉得还是容易误导的。间隙本没电子,何来磁矩?计算出极化电荷在间隙里,只能说明是原子互相作用,在间隙产生了一个新轨道。所以间隙磁矩说到底还是多个原子局域磁矩的整体表现,由轨道所有者共同控制。要调整间隙里的电子态,调整周围原子的初始磁矩就行了
2楼2014-02-17 17:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
2楼: Originally posted by golddoushi at 2014-02-17 17:47:59
无缺陷无替换时是共价键,不代表你弄过之后还是共价键。这个要看具体情况。
OUTCAR里只有wigner半径内的磁矩,通常都是不准的。一般bader方法算的极化电荷就能跟OSZICAR对上

间隙诱发/局域/磁矩,虽然大家都这么 ...

bader电荷分析好像只能显示电荷的分布情况,但是不能显示磁性,我最终想要的还是磁性方面的结果,就没有别的办法么
而且初始磁矩也调整过几种情况,但是结果大同小异
3楼2014-02-18 08:52:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

golddoushi

木虫 (正式写手)

引用回帖:
3楼: Originally posted by yangbs2009 at 2014-02-18 08:52:19
bader电荷分析好像只能显示电荷的分布情况,但是不能显示磁性,我最终想要的还是磁性方面的结果,就没有别的办法么
而且初始磁矩也调整过几种情况,但是结果大同小异...

磁矩分布不就是电荷分布么,上下旋电荷差分就行。你调初始磁矩要先弄清楚间隙轨道是哪几个原子公有
4楼2014-02-18 13:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
4楼: Originally posted by golddoushi at 2014-02-18 13:03:55
磁矩分布不就是电荷分布么,上下旋电荷差分就行。你调初始磁矩要先弄清楚间隙轨道是哪几个原子公有...

初始磁矩要先弄清楚间隙轨道是哪几个原子共有,这个知道
bader电荷分析结果中可以分别显示上下自旋的电荷么?
5楼2014-02-18 13:54:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

golddoushi

木虫 (正式写手)

引用回帖:
5楼: Originally posted by yangbs2009 at 2014-02-18 13:54:53
初始磁矩要先弄清楚间隙轨道是哪几个原子共有,这个知道
bader电荷分析结果中可以分别显示上下自旋的电荷么?...

嗯可以。自旋计算出来的CHGCAR里是分自旋态记录的,写个小程序(或网上搜个对应版本的chgsplit)分开,比如chg1, chg2, 然后用bader减一下, bader chg1 -ref chg2 ,这样给出来的就是净上/下旋电荷,也就是磁矩啦
6楼2014-02-18 14:22:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangbs2009

木虫 (小有名气)

引用回帖:
6楼: Originally posted by golddoushi at 2014-02-18 14:22:53
嗯可以。自旋计算出来的CHGCAR里是分自旋态记录的,写个小程序(或网上搜个对应版本的chgsplit)分开,比如chg1, chg2, 然后用bader减一下, bader chg1 -ref chg2 ,这样给出来的就是净上/下旋电荷,也就是磁矩啦...

ok,谢谢啦,我先试试看
7楼2014-02-18 14:25:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangbs2009 的主题更新
信息提示
请填处理意见