24小时热门版块排行榜    

查看: 3135  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +9 z1z2z3879 2026-03-16 12/600 2026-03-18 08:53 by zhukairuo
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 070300化学319求调剂 +4 锦鲤0909 2026-03-17 4/200 2026-03-17 18:21 by 重科小霸王
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 085600材料与化工 +4 安全上岸! 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
信息提示
请填处理意见