24小时热门版块排行榜    

查看: 2896  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

空谷新苗

木虫 (小有名气)

[求助] 求助:vasp计算TiO2薄膜无带隙。已有8人参与

vasp计算锐钛矿相TiO2薄膜,具体步骤是:
  1.用优化好的体材料的晶格参数建模。因计算的是纯薄膜的性质,故只建了1*1的三层薄膜(后来又建了一个5层的薄膜,但结果还是无带隙,感觉应该不是薄膜层数的问题)。
  2.由建好的薄膜模型得到POSCAR。POSCAR采用Selective Dynamic,固定最下层电子,让上面两层原子只在z方向弛豫(设置为FFT).
  3.进行结构优化。INCAR文件为:
ENCUT = 500
PREC = Accurate
EDIFF = 1E-6
EDIFFG = -0.01
NPAR=6
IALOGO=48
ISTART = 0
ISIF=2
ICHARG = 2
ISMEAR = 1
SIGMA = 0.1
IBRION = 1
NSW =120
POTIM = 0.2

K点取3*7*1
POSCAR文件为:
anatase CIF
1.00000000000000   
       10.304135548894921       0.000000000000000      -0.000000000000000
       0.000000000000000       3.797100000000000       0.000000000000000
       0.000000000000000       0.000000000000000      23.440695973720747
O    Ti
24    12
Selective dynamics
Direct
0.4130356769206680   0.0000000000000001   0.0249373600250313   F   F   F

0.5488297901615880   0.0000000000000001   0.1755252828815010   F   F   T

0.6846239034025090   0.0000000000000001   0.3261132057379700   F   F   T

0.9130356769206680   0.5000000000000000   0.0249373600250313   F   F   F

0.0488297901615880   0.5000000000000000   0.1755252828815010   F   F   T

0.1846239034025090   0.5000000000000000   0.3261132057379700   F   F   T

0.4079937566986610   0.5000000000000000   0.1129409421423520   F   F   F

。。。。。
。。。。。
  4.结构优化后将CONTCAR文件复制到POSCAR中,进行静态自洽计算。INCAR文件为:
ENCUT = 500
PREC = Accurate
EDIFF = 1E-6
EDIFFG = -0.01
NELM=150
IALOGO=48
NPAR=6
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
IBRION=-1
NSW=0
ISPIN=2

K点取5*9*1

4.计算DOS。静态自洽计算后修改INCAR:ISTART = 1;ICHARG = 11;ISMEAR = -5,LORBIT=11.增加K点,然后计算。
最后用p4vasp得出的DOS图为:求助:vasp计算TiO2薄膜无带隙。

我想请教一下为什么DOS图无带隙,我哪个环节出错了,该如何修改,谢谢。
回复此楼

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

光催化专辑

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

读万卷书,行万里路。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

空谷新苗

木虫 (小有名气)

引用回帖:
2楼: Originally posted by xxxxxxxfc at 2015-05-22 14:41:09
我很好奇,你们要研究的薄膜是不是只有3层原子,实际薄膜有这么薄吗?我们这儿磁控溅射一般也有几十到几百纳米。如果你是研究表面那可能确实应该这样做。

谢谢您的回复!因为我现在刚开始计算薄膜,所以想先建一个小点的体系练习一下并检验一下计算方法的正确性,所以只建了3层。虽然体系小了点,但也不该由半导体变为导体呀。并且我看有些文献上只建了5层薄膜,结果仍是半导体,而我后来用5层薄膜计算的结果仍是无带隙,所以感觉不是层数的原因。您做过薄膜计算吗?一般取几层合适?

[ 发自小木虫客户端 ]
读万卷书,行万里路。
3楼2015-05-22 19:09:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

xxxxxxxfc

金虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
ljw4010: 金币+1, 谢谢交流! 2015-05-23 20:46:16
我很好奇,你们要研究的薄膜是不是只有3层原子,实际薄膜有这么薄吗?我们这儿磁控溅射一般也有几十到几百纳米。如果你是研究表面那可能确实应该这样做。
2楼2015-05-22 14:41:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

golddoushi

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
fzx2008: 金币+2, 谢谢回复 2015-05-23 08:57:43
主要是你没加U,TiO2即使块体模型按你这么算也才0.8eV gap 远小于实验3.6,变表面后不饱和了进一步缩小gap就成你这样了
+U,GW或者杂化泛函,具体参考说明书相关章节
4楼2015-05-22 22:36:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

空谷新苗

木虫 (小有名气)

引用回帖:
4楼: Originally posted by golddoushi at 2015-05-22 22:36:44
主要是你没加U,TiO2即使块体模型按你这么算也才0.8eV gap 远小于实验3.6,变表面后不饱和了进一步缩小gap就成你这样了
+U,GW或者杂化泛函,具体参考说明书相关章节

谢谢您的回复! 我明天+U计算一下,不过我在计算体材料时未加U计算的带隙也有2.2eV。希望明天能出现奇迹。

[ 发自小木虫客户端 ]
读万卷书,行万里路。
5楼2015-05-22 22:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见