| 查看: 3755 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
Liz_Tesla木虫 (小有名气)
|
[求助]
VASP 自旋轨道耦合 问题 已有1人参与
|
||
|
最近在尝试用VASP计算自旋轨道耦合，用拓扑绝缘体Bi2Se3来检验，发现加入SO后能隙并没有打开而是减少。下面是我的步骤: 1）正常优化结构，进行自洽静态计算，得到WAVECAR和CHGCAR（这一步没有加SO） 2）加入SO，其中将第一步得到的IBZKPT拷贝成KPOINTS，并将INCAR改写成 SYSTEM = Bi2Se3 ENCUT = 400 ISTART = 1 ; ICHARG = 11 #ISMEAR = -5 PREC = Normal LORBIT=11 #LELF=.TRUE LSORBIT= .TRUE LMAXMIX=2 MAGMOM=15*0 SAXIS= 0 0 1 NBANDS=64 ISYM=0 GGA_COMPAT=.FALSE 其中LMAXMIN=2，是因为Bi和Se都是 s element 3）计算能带。将KPOINTS修改为line-mode取点。INCAR改为 SYSTEM = Bi2Se3 ENCUT = 400 ISTART = 1 ; ICHARG = 11 #ISMEAR = -5 PREC = Normal LORBIT=11 #LELF=.TRUE LSORBIT= .TRUE LMAXMIX=2 #MAGMOM=15*0 #SAXIS= 0 0 1 #NBANDS=64 ISYM=0 GGA_COMPAT=.FALSE 下面两幅图是我得到的能带图。不加SO时能带图与文献符合。加入SO后能隙不是打开而是有点闭合（费米能级为3.65） 另外我的步骤有没有问题？谢谢  no SO.png  with SO.png |
回复此楼» 收录本帖的淘帖专辑推荐
 自旋轨道耦合 |
» 猜你喜欢
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有7人回复
-
26申博
已经有3人回复
-
求调剂，总分315，考的生物医药，一志愿湖南师范大学。调剂到任何专业都可以
已经有3人回复
-
材料专硕306英一数二
已经有7人回复
-
材料工程专硕274一志愿211求调剂
已经有6人回复
-
一志愿南京大学，080500材料科学与工程，调剂
已经有4人回复
-
332求调剂
已经有5人回复
-
271求调剂
已经有14人回复
-
275求调剂
已经有4人回复
-
302求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP 铁的磁性设置
已经有20人回复
- 用VASP计算了一个结构，怎么看磁矩呢？
已经有14人回复
- 关于自旋轨道耦合能量最低方向的判断。
已经有5人回复
- vasp中有关自旋
已经有31人回复
- 关于VASP中MAGMON的磁矩设置问题
已经有6人回复
- 自旋轨道耦合后怎么算态密度
已经有4人回复
- VASP自旋轨道耦合
已经有26人回复
- NBANDS 的设置问题
已经有6人回复
- 关于 VASP 的编译经验分享
已经有16人回复
- 自旋轨道耦合第一性原理计算
已经有11人回复
- vasp自旋极化计算时，初始设置自旋能否不被优化
已经有6人回复
- VASP 中NELECT 的设置
已经有8人回复
- VASP如何读取分子轨道？
已经有18人回复
- 使用vasp计算cineb一定要编译vtstcode吗？
已经有8人回复
- 求助vasp编译好后无法运行问题？？
已经有24人回复
- ＶＡＳＰ计算出错，请高手指点
已经有19人回复
- molpro试做旋轨耦合,怎么在输出中找omega及其对应的能量
已经有5人回复
- vasp为何每次遇见诸如wrap around error must be expected set NGX to 这类问题
已经有12人回复
- 求助vasp中的自旋轨道耦合
已经有3人回复
- vasp计算磁矩的困惑
已经有10人回复
- 【求助】VASP考虑自旋轨道耦合的话，如何编译
已经有17人回复
Liz_Tesla
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1435.8
- 散金: 902
- 红花: 1
- 帖子: 285
- 在线: 424.7小时
- 虫号: 2445360
- 注册: 2013-05-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2015-03-30 09:12:55
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2015-03-30 02:27:03
Liz_Tesla
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1435.8
- 散金: 902
- 红花: 1
- 帖子: 285
- 在线: 424.7小时
- 虫号: 2445360
- 注册: 2013-05-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-03-30 09:49:23
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。。 2015-03-31 10:23:52
Liz_Tesla: 金币+30, ★★★★★最佳答案 2015-03-31 12:24:56
franch: 金币+3, 谢谢回帖交流。。。。 2015-03-31 10:23:52
Liz_Tesla: 金币+30, ★★★★★最佳答案 2015-03-31 12:24:56
|
你試試看照文獻的原子位置去設定POSCAR 去掉relax那一步 只做一次自洽計算 直接開SOC 其他設定還是照我之前告訴你的那樣去改 對了 你是不是POSCAR是不是設15顆原子的超胞? 我建議你用5顆原子的 最小晶胞 因為?埵钻苫蚱渌鸅i2Se3的文獻大部分在算 bulk band structure時 都是用最小晶胞的 你要跟他們的bulk band structure比較 你一定要用最小晶胞 這樣KPOINTS路?讲拍茉O一樣的 除非你是要算surface state 不然不該用超胞 |
5楼2015-03-30 14:05:25
