24小时热门版块排行榜

返回列表

当前主题已经存档。

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

jugengfans

金虫 (正式写手)

应助: 3 (幼儿园)
金币: 773.9
散金: 151
红花: 3
帖子: 537
在线: 205.6小时
虫号: 291609
注册: 2006-10-29
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】TiO2 的能带图，我画出来的，大家看看对么？

Vasp刚学，拿来入手。金红石结构的TIO2.
画出来的能带图很怪异
大家看看我的设置有问题么？

INCAR:
SYSTEM = TIO2
ISTART = 1
ICHARG = 11
ENCUT = 520
ISPIN = 2
#MAGMOM = 8*0.0 8*1.0  24*0.0
EDIFF = 1E-4
EDFIIG = -1E-2
ISMEAR = 2
NPAR = 2
AMIX    = 0.2
BMIX    = 0.00001
AMIX_MAG = 0.8
BMIX_MAG = 0.00001
LREAL = .TRUE.
Dynamic:
#IBRION = 2
#NSW = 500
#POTIM = 0.2
LMAXMIX = 4       ! Important: mix paw occupancies up to L=4

POSCAR:
POSCAR FOR FILM
1.000000000000000
   4.5940000000000000 0.0000000000000000 0.0000000000000000
   0.0000000000000000 4.5940000000000000 0.0000000000000000
   0.0000000000000000 0.0000000000000000 2.9590000000000000
2 4
Selective dynamics
Direct
0.5000000000000000 0.5000000000000000 0.5000000000000000 T T T
0.0000000000000000 0.0000000000000000 0.0000000000000000 T T T
-0.3047997387897260  -0.3047997387897260 0.0000000000000000 T T T
0.3047997387897260 0.3047997387897260 0.0000000000000000 T T T
0.8047997387897259 0.1952002612102740 0.5000000000000000 T T T
0.1952002612102740 0.8047997387897259 0.5000000000000000 T T T

POTCAR：
GGA-PAW势

KPOINTS：
AUTO
0
Monkhorst
10 10 10
0  0  0

就这些了，但是画出EIGENVAL感觉能带震荡特别剧烈，不像一般在文献中看到的，是不是我上面的设置出了问题了？请教大家！

[ Last edited by jugengfans on 2009-8-16 at 01:31 ]

回复此楼

» 猜你喜欢

1楼 2009-08-15 23:35:39

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

jugengfans

金虫 (正式写手)

应助: 3 (幼儿园)
金币: 773.9
散金: 151
红花: 3
帖子: 537
在线: 205.6小时
虫号: 291609
注册: 2006-10-29
专业: 凝聚态物性 II ：电子结构

★ ★ ★ ★ ★ ★ ★ ★
aylayl08(金币+5):感谢讨论交流 2010-03-13 08:43
wuchenwf(金币+3):谢谢 2010-03-13 21:31

引用回帖:

Originally posted by xiaojie7783 at 2009-9-8 11:10:
你好我想请问你用VASP计算能带的时候的KPOINTS文件的写法现在弄对了吗我也遇到你上面提到的问题就是K点分割的问题，不知道要按照什么原则分割倒空间，望指教，谢谢。

已经对了。关于高对称点的选取有一些简单的规则，
第一是对称性。不同的晶体结构有不同的对称性，我们应该根据晶体结构选择合适的K点，这个在各种教材中已经给出了。不过手动选择K点只在能带计算的时候是必要的，在做优化和性质计算时我们通常选择自动撒点。
第二是K点分割，K点分割越密，计算的结果就越精确，而耗时也会更多，因此我们要寻找一个比较合适的分割。这个可以在你计算的时候做一个测试，就是从小到大选择不同的分割，然后比较静态计算时的总能变化，自己画一个曲线图，当总能变化范围小于0.2eV时，我们就认为不用再取更多的K点了。

一般来说，K点的选取跟晶胞的边长有关，边长越长，所需的K点就越少。如果是块体的单胞，K点的分割要取密一些，比如我们的TiO2单胞，K点分割12*12*12就非常精确了。而如果我们将这个晶胞延xyz方向各扩展一倍，这时只要6*6*6足够。越大的晶胞，所需的K点数就越小。如果是表面，那么在某个方向有大于10A的真空层，那么这个方向上K点最多只须取2足矣。