24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Liz_Tesla

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1435.8
散金: 902
红花: 1
帖子: 285
在线: 424.7小时
虫号: 2445360
注册: 2013-05-02
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 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

回复此楼

1楼 2015-03-29 17:52:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Liz_Tesla

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1435.8
散金: 902
红花: 1
帖子: 285
在线: 424.7小时
虫号: 2445360
注册: 2013-05-02
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

5楼: Originally posted by skydark at 2015-03-30 14:05:25
你試試看照文獻的原子位置去設定POSCAR
去掉relax那一步
只做一次自洽計算直接開SOC
其他設定還是照我之前告訴你的那樣去改

對了你是不是POSCAR是不是設15顆原子的超胞?
我建議你用5顆原子的最小晶胞
...

嗯，我试试
POSCAR用的是含五个原子的最小晶胞

赞一下

回复此楼

6楼2015-03-30 14:21:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 12 个回答

skydark

木虫 (正式写手)

1ST强帖: 1
应助: 44 (小学生)
金币: 3652.2
散金: 30
红花: 6
帖子: 302
在线: 249.8小时
虫号: 1927082
注册: 2012-08-06
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+4, 谢谢回帖交流。。。。 2015-03-31 10:23:44

第2步就弄成 ICHARG = 11?
那我是不會這樣做
算Bi2Se3 第2步時
我都設ICHARG = 2
另外MAGMOM 我也是直接用預設值的
LMAXMIX=2 也是有問題
應該要等於4

赞一下

回复此楼

2楼2015-03-30 02:27:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Liz_Tesla

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1435.8
散金: 902
红花: 1
帖子: 285
在线: 424.7小时
虫号: 2445360
注册: 2013-05-02
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by skydark at 2015-03-30 02:27:03
第2步就弄成 ICHARG = 11?
那我是不會這樣做
算Bi2Se3 第2步時
我都設ICHARG = 2
另外MAGMOM 我也是直接用預設值的
LMAXMIX=2 也是有問題
應該要等於4

1）ICHARG=11是按照VASP教程设置的，用于读取第一步的WAVECAR和CHGCAR
2）LMAXMIX=4不是对于d element吗？
3）你用这些参数能看出SO打开能隙？
我用你的参数试试