24小时热门版块排行榜    

查看: 2892  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 321求调剂 +4 大米饭! 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +4 Liwangman 2026-03-15 4/200 2026-03-16 08:01 by wang_dand
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
请填处理意见