24小时热门版块排行榜    

查看: 2891  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
信息提示
请填处理意见