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 ]

回复此楼

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

VASP资料

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

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

whzhang81

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 217.1
帖子: 32
在线: 9.8小时
虫号: 815235
注册: 2009-07-24
专业: 理论和计算化学

★ ★ ★
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