24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2940  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物与医药273求调剂 +6 荔题南墙 2026-04-05 6/300 2026-04-06 08:35 by zhangdingwa
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 338求调剂 +3 我想上岸ii 2026-04-05 3/150 2026-04-05 19:59 by nepu_uu
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 288求调剂 +7 没有答案_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[考研] 材料调剂 +12 一样YWY 2026-04-04 12/600 2026-04-05 08:24 by 544594351
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 数二英二348求调剂 +4 hxdzj1 2026-04-03 5/250 2026-04-03 21:25 by zhq0425
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考研] 一志愿厦门大学材料工程专硕354找调剂!!! +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 0805求调剂 +8 是水分 2026-03-31 8/400 2026-04-02 10:46 by guanxin1001
[考研] 302求调剂一志愿北航070300,本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
信息提示
请填处理意见