24小时热门版块排行榜    

查看: 2768  |  回复: 12

wmh19889

铜虫 (初入文坛)

[交流] 【求助】vasp noncollinear计算能带的时候只能算spin-unpolarized吗 已有4人参与

vasp  noncollinear计算能带的时候只能算spin-unpolarized,即使设成spin-polarized
也会被改回spin-unpolarized,觉得很奇怪,
我是按VASP推荐的做法,先算non-magnetic groundstate to get a WAVECAR and CHGCAR file,然后再加磁,MAGMOM,ISPIN,LSORBIT,SAXIS等设置
但结果还只有spin-unpolarized,
大多磁性体系都是spin-polarized才对,有办法算spinplarized bandstructure吗?

[ Last edited by wmh19889 on 2009-9-18 at 21:20 ]
回复此楼

» 猜你喜欢

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

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

wei.ji

木虫 (正式写手)

★ ★
wmh19889(金币+1):谢谢参与
zxzj05(金币+1,VIP+0):鼓励讨论! 9-28 17:18
只需要电荷密度,不需要波函数,另外,你肯定是什么地方没搞对,呵呵~
http://sim.phys.ruc.edu.cn
2楼2009-09-23 21:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apzhough

铁虫 (小有名气)


wmh19889(金币+1):谢谢参与
路过,踩踩
3楼2009-09-28 15:33:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emily6526

铁虫 (初入文坛)

★ ★
wmh19889(金币+1):谢谢参与
zxzj05(金币+1,VIP+0):3Q! 鼓励讨论! 10-13 14:30
我在做非共线磁性,好像得到的DOSCAR里面没有spin polarized ,即使INCAR里ISPIN=2,OUTCAR中SPIN会是1。不知是怎么回事?期待高手~
4楼2009-10-12 19:42:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mingyan988

金虫 (小有名气)


wmh19889(金币+1):谢谢参与
我也遇到这样的问题,不知道怎么处理了。我认为应该是INCAR文件里的参数设置存在冲突了,寻求高手指点!
岁月轻狂
5楼2011-01-26 13:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gloryshine

铜虫 (正式写手)


wmh19889(金币+1): 谢谢参与
引用回帖:
5楼: Originally posted by mingyan988 at 2011-01-26 13:53:18
我也遇到这样的问题,不知道怎么处理了。我认为应该是INCAR文件里的参数设置存在冲突了,寻求高手指点!

我算的能带结果的OUTCAR中显示实际计算中ISPIN=1,尽管我在INCAR中明明输入的是ISPIN=2。不知道怎么解决,现在画能带图就分离不了上下自旋结果。请高手解答。
6楼2015-04-03 22:50:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by gloryshine at 2015-04-03 22:50:08
我算的能带结果的OUTCAR中显示实际计算中ISPIN=1,尽管我在INCAR中明明输入的是ISPIN=2。不知道怎么解决,现在画能带图就分离不了上下自旋结果。请高手解答。...

贴出你的INCAR啊
凝聚态理论与计算物理
7楼2015-04-03 23:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
ljw4010: 金币+3, 谢谢指导解答! 2015-04-04 21:32:13
考虑非共线极化时,哈密顿量不再与Sz对易,因此Sz不再是守恒量。也就是说,如果某个态有确定的能量,则一定没有确定的自旋Z分量。能带图中的态都是哈密顿算符本征态,因此都没有确定的自旋Z分量。

自旋非极化时电荷密度表示为rho(r),自旋极化时变成一个对角矩阵[rho_aa(r), 0; 0, rho_bb(r)],非共线磁化时变成一般矩阵[rho_aa(r), rho_ab(r); rho_ba(r), rho_bb(r)]。

这个可以用泡利算符验证。SigmaZ本征态为[ 1; 0 ]和[0; 1],任取一个对角矩阵[ 2, 0; 0, 3 ],可以验证[ 1; 0 ]和[0; 1]仍是这个对角矩阵的本征态,本征值分别为2和3.而任取一个非对角矩阵,[ 1; 0 ]和[0; 1]就不一定是它的本征态。
8楼2015-04-04 10:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
另一个涉及自旋的磁化现象是自选轨道耦合。这时自旋角动量不是守恒量,但总角动量仍是守恒量,本征值是j +- 1/2。因此考虑自旋轨道仍能区分出两组不同的能带。

非共线磁化时轨道角动量没被考虑进来,自旋角动量又不守恒,因此总角动量不守恒,分不出自旋向上和自旋向下。
9楼2015-04-04 10:39:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
ljw4010: 金币+3, 谢谢指导解答! 2015-04-04 21:32:28
vasp  noncollinear计算能带的时候只能算spin-unpolarized,即使设成spin-polarized
也会被改回spin-unpolarized,觉得很奇怪,
>   簡單說就是
     OUTCAR 內看到的 ISPIN=1 ?K不是指 spin-unpolarized calculation 的意思
     在 noncollinear 計算下,INCAR 內 ISPIN 不管給多少~最後在 OUTCAR 內看到的一定是 ISPIN=1
     原因很簡單
    如同 liqizuiyang 所言
    在 noncollinear 計算下,是沒有 up/down 的概念存在
    所以 VASP 很自然的就以 ISPIN=1 呈現
    你只需要確定OUTCAR內的 noncollinear 相關參數 (如 LSORBIT, or LNONCOLLINEAR) 是否為 .TRUE. 即可



大多磁性体系都是spin-polarized才对,有办法算spinplarized bandstructure吗?
>  可以
    就跟一般算 non-spin polarized 計算的 band structure 步驟是一樣的
10楼2015-04-04 21:19:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wmh19889 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见