24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 3163  |  回复: 16

callen123

金虫 (小有名气)

[求助] 关于自旋轨道耦合的问题 已有1人参与

vasp已重新编译,我的算SOC的过程是开始做结构优化(ISPIN=2),然后再做一步共线磁性计算(我的理解是spin打开,加入LMAXMIX参数,做一步静态计算,得到CHGCAR和WAVECAR),之后在进行自旋轨道耦合计算(这一步也是静态的吗?还有如果算后面再算态密度DOS的话,是算完自旋轨道耦合在算DOS还是可以在算SOC时加入DOS参数一起算?)。但是我算完SOC后,print-out里面提示WAVECAR未读取,如图片所示。还有我算SOC时,KPOINTS直接copy的是共线磁性计算那一步的KPOINTS,不是copy那个IBZKPT, 之前也有试过直接copy KPOINTS到IBZKPT,但是好像不行,难道是copy完还要往KPOINTS里面再加什么吗? 求指导 十分迷茫 555

关于自旋轨道耦合的问题
QQ Photo20141121231509.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

VASP and MS 自旋轨道耦合

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dx0620

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
callen123: 金币+6, ★★★很有帮助 2014-11-24 14:46:30
ljw4010: 金币+2, 谢谢交流! 2014-11-26 14:24:54
我计算时加SOC的处理过程(亲测可行):
1、正常优化
2、加SOC静态计算,得到电荷密度,不用管WAVECAR
3、加SOC,用上一步得到的电荷密度计算能带、态密度(侯柱峰手册方法即可)

需要注意的是,我不计算磁性性质,只是需要考虑SOC对能带的影响,INCAR中只加入了下列参数:

LSORBIT = .TRUE.
GGA_COMPAT = .FALSE.
LMAXMIX = 4
ISYM = 0
NBANDS = 48

与文献对比,能带计算的没有问题,态密度有点问题(这个还在检验中)
2楼2014-11-23 20:34:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

callen123

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dx0620 at 2014-11-23 20:34:35
我计算时加SOC的处理过程(亲测可行):
1、正常优化
2、加SOC静态计算,得到电荷密度,不用管WAVECAR
3、加SOC,用上一步得到的电荷密度计算能带、态密度(侯柱峰手册方法即可)

需要注意的是,我不计算磁性 ...

我也不考虑磁性 但是我算完SOC之后,在算DOS时,没有他们所谓的DOSCAR里面有37列,而且数据主要集中在S轨道,感觉不对啊。。 而且SOC算完的OUTCAR里面 ISPIN变成1了,INCAR我设的是2  求指教
4楼2014-11-24 14:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

callen123

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dx0620 at 2014-11-23 20:34:35
我计算时加SOC的处理过程(亲测可行):
1、正常优化
2、加SOC静态计算,得到电荷密度,不用管WAVECAR
3、加SOC,用上一步得到的电荷密度计算能带、态密度(侯柱峰手册方法即可)

需要注意的是,我不计算磁性 ...

请问你算完DOSCAR里面有他们说的37列吗? 还有你的DOS是算完SOC之后在进行的  还是加到SOC计算的INCAR里面的 我试了下发现ISMEAR=-5情况下 SOC运行会出错  如果是跑完SOC在计算DOS的话 请问你应该不是用的split_dos处理的吧 那该用什么脚本处理的  不好意思 刚刚接触问题有点多 期待您的指点
3楼2014-11-23 22:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
4楼: Originally posted by callen123 at 2014-11-24 14:36:50
我也不考虑磁性 但是我算完SOC之后,在算DOS时,没有他们所谓的DOSCAR里面有37列,而且数据主要集中在S轨道,感觉不对啊。。 而且SOC算完的OUTCAR里面 ISPIN变成1了,INCAR我设的是2  求指教...

对,我计算DOS的错误也差不多,集中在一个轨道,这个问题还在解决中。你最好确认下这个37列是怎么回事,我看到过这个帖子,后面貌似有说明,不一定会有37列数据。
加SOC的话,ISPIN这个参数好像没用吧
5楼2014-11-24 15:04:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

callen123

金虫 (小有名气)

引用回帖:
5楼: Originally posted by dx0620 at 2014-11-24 15:04:12
对,我计算DOS的错误也差不多,集中在一个轨道,这个问题还在解决中。你最好确认下这个37列是怎么回事,我看到过这个帖子,后面貌似有说明,不一定会有37列数据。
加SOC的话,ISPIN这个参数好像没用吧...

SOC做完分割DOSCAR 用的也是split_dos吗 会不会不一样的 还有手册上说是先是线性静态计算 然后才是SOC非共线计算(读取CHGCAR和WAVECAR)  怎么大家都在静态时加了SOC参数 这是对还是不对啊  真纠结
6楼2014-11-24 15:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

callen123

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dx0620 at 2014-11-23 20:34:35
我计算时加SOC的处理过程(亲测可行):
1、正常优化
2、加SOC静态计算,得到电荷密度,不用管WAVECAR
3、加SOC,用上一步得到的电荷密度计算能带、态密度(侯柱峰手册方法即可)

需要注意的是,我不计算磁性 ...

你计算能带时,是直接用的SOC算完的数据,还是算完SOC再算一次能带计算的?
7楼2014-11-25 08:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
7楼: Originally posted by callen123 at 2014-11-25 08:49:08
你计算能带时,是直接用的SOC算完的数据,还是算完SOC再算一次能带计算的?...

SOC不是特别算的,2L说的很明白了,静态计算的时候加SOC,能带也加
8楼2014-11-25 09:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

callen123

金虫 (小有名气)

引用回帖:
8楼: Originally posted by dx0620 at 2014-11-25 09:03:17
SOC不是特别算的,2L说的很明白了,静态计算的时候加SOC,能带也加...

我明白了 我所谓的SOC计算(静态算完之后那一步)就是你说的能带计算 谢谢了
9楼2014-11-25 09:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

callen123

金虫 (小有名气)

引用回帖:
5楼: Originally posted by dx0620 at 2014-11-24 15:04:12
对,我计算DOS的错误也差不多,集中在一个轨道,这个问题还在解决中。你最好确认下这个37列是怎么回事,我看到过这个帖子,后面貌似有说明,不一定会有37列数据。
加SOC的话,ISPIN这个参数好像没用吧...

LORBIT=11  DOSCAR就有37列 =10就没有那么多
10楼2014-11-25 10:44:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 callen123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 347材料专硕求调剂 +18 zj8215216 2026-04-06 18/900 2026-04-08 16:27 by luoyongfeng
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 293调剂 +18 yj1221 2026-04-08 19/950 2026-04-08 14:33 by 418490947
[考研] 085404,285分求调剂 +4 薇薇考研 2026-04-07 4/200 2026-04-08 14:32 by zhq0425
[考研] 求调剂 +13 熊二想上岸 2026-04-04 13/650 2026-04-08 10:33 by 不吃魚的貓
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 一志愿华南师范大学0702物理学305调剂 +3 念常安 2026-04-07 5/250 2026-04-08 06:45 by lijunpoly
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 材料与化工371求调剂 +14 陪琳看海 2026-04-04 15/750 2026-04-06 06:59 by houyaoxu
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 324求调剂 +14 想上学求调 2026-04-02 15/750 2026-04-04 20:31 by 无际的草原
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 一志愿华中农业071010,总分320求调剂 +7 困困困困坤坤 2026-04-02 7/350 2026-04-03 17:26 by Yuena_Wang
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
信息提示
请填处理意见