| 查看: 3831 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 自旋轨道耦合 |
» 猜你喜欢
297求调剂
已经有16人回复
-
326分，一志愿沪9，求生物学调剂
已经有4人回复
-
288求调剂，一志愿华南理工大学071005
已经有3人回复
-
304求调剂
已经有13人回复
-
调剂
已经有18人回复
-
293调剂
已经有4人回复
-
338求调剂
已经有7人回复
-
263求调剂
已经有3人回复
-
求调剂
已经有15人回复
-
机械工程264学硕求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
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
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 ：电子结构
3楼2015-03-30 09:12:55
Liz_Tesla
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1435.8
- 散金: 902
- 红花: 1
- 帖子: 285
- 在线: 424.7小时
- 虫号: 2445360
- 注册: 2013-05-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-03-30 09:49:23
Liz_Tesla
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1435.8
- 散金: 902
- 红花: 1
- 帖子: 285
- 在线: 424.7小时
- 虫号: 2445360
- 注册: 2013-05-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-03-30 14:21:49
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 26041.2
- 散金: 244
- 红花: 13
- 帖子: 4087
- 在线: 939.1小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
7楼2015-03-30 21:12:53
Liz_Tesla
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1435.8
- 散金: 902
- 红花: 1
- 帖子: 285
- 在线: 424.7小时
- 虫号: 2445360
- 注册: 2013-05-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2015-03-31 09:27:42
Liz_Tesla
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1435.8
- 散金: 902
- 红花: 1
- 帖子: 285
- 在线: 424.7小时
- 虫号: 2445360
- 注册: 2013-05-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-03-31 12:24:51
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2015-03-31 12:39:47
