| 查看: 2441 | 回复: 18 | ||
| 当前主题已经存档。 | ||
jugengfans金虫 (正式写手)
|
[交流]
【求助】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 ] |
回复此楼
2楼2009-08-15 23:54:33
3楼2009-08-15 23:57:15
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
4楼2009-08-16 00:10:23
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
★ ★
fegg7502(金币+2,VIP+0):鼓励交流 8-16 17:43
fegg7502(金币+2,VIP+0):鼓励交流 8-16 17:43
|
你说的对，但具体怎么设置？ 我手动生成K点用的是如下的KPOINTS band 10 Line-mode Cartesian 0 0 0 !GAMMA 0 1 0 !F 0 1 0 !F 0 1 1 !Q 0 1 1 !Q 0 0 1 !Z 0 0 1 !Z 0 0 0 !GAMMA 0 0 0 !GAMMA 1 0 0 !B 是按照CASTEP给出的K点顺序生成的 图很怪异！震荡非常古怪！请教大家了！ [ Last edited by jugengfans on 2009-8-16 at 01:31 ] |
5楼2009-08-16 00:20:57
chuzhaonan
铁杆木虫 (著名写手)
第一性原理计算大师
- 应助: 3 (幼儿园)
- 金币: 7049.8
- 散金: 120
- 红花: 3
- 帖子: 2282
- 在线: 129小时
- 虫号: 515926
- 注册: 2008-03-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2009-08-16 10:49:38
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
7楼2009-08-16 21:34:22
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
8楼2009-08-16 23:32:22
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
★ ★
fegg7502(金币+2,VIP+0):鼓励交流2 9-6 18:37
fegg7502(金币+2,VIP+0):鼓励交流2 9-6 18:37
|
上文中所画的 TiO2能带图我们采用的是GGA-PAW势， Ti元素有三种POTENTIAL可选，分别是Ti Ti_sv Ti_pv 其中_sv和_pv表示把s或p电子作为半芯态来处理，（这句话是什么意思还未查明） 我们选用的是Ti 接下来我会用Ti_pv的赝势进行对比计算。 最下面的图是Ti_Pv的浅能级图，感觉和我第一次回复时用Ti势画出的图波动很相似，看来没有太大的问题了。 最后一个问题，我的KPOINTS只能在每两个K点间统一插入10个值。那么如何在计算band时让vasp在不同的高对称点之间插入不同的分割值？请大家不吝赐教！ 还有一个图是用castep计算得到的，也是论坛里的朋友发的。 不知道他用的是什么赝势，不过能带分布和gap看起来差不多的。 [ Last edited by jugengfans on 2009-8-17 at 00:26 ] |
9楼2009-08-16 23:49:16
★ ★ ★
yjr(金币+3,VIP+0):谢谢！！ 8-17 21:38
yjr(金币+3,VIP+0):谢谢！！ 8-17 21:38
|
你这个，不知道怎么算的，你的INCAR里边首先有个问题，对于Rutile TiO2非金属计算的时候ISMEAR 应该用0，不能用2，不过这个应该不是你能带成这样的原因。 我用DMol3低精度跑了一个能带图 对于高对称点的坐标和命名，以前有人贴过，对于你这个晶体，如下图所示，你可以对照一下看看。 k-vector label Wyckoff position Parameters CDML ITA ITA GM 0,0,0 1 a 4/mmm 0,0,0 Z 0,0,1/2 1 b 4/mmm 0,0,1/2 M 1/2,1/2,0 1 c 4/mmm 1/2,1/2,0 A 1/2,1/2,1/2 1 d 4/mmm 1/2,1/2,1/2 R 0,1/2,1/2 2 e mmm. 0,1/2,1/2 X 0,1/2,0 2 f mmm. 0,1/2,0 |
10楼2009-08-17 19:53:34
5