24小时热门版块排行榜

返回列表

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

utimasa

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3
帖子: 10
在线: 16.4小时
虫号: 1243869
注册: 2011-03-24
专业: 凝聚态物性 II ：电子结构

[求助] 求助VASP 处理带电缺陷的问题！能量变化太大

本人正在算一个半导体氧化物中引入N杂质，N替换O，计算其形成能。
首先计算了中性缺陷，体系的总能： ETOT= -648.481449 eV
然后计算N（-1）,体系的总能，将NELECT设为NELECT(原来的中性体系）+1
得到ETOT=-643.658924 eV

竟然相差4.82eV。
也就是说如果N得到一个电子，体系的能量竟然增加了4.82eV.
其他参数都是一样，没有变化。
根据实验，N是很好的受主，应该是容易得到电子。
百思不得其解。
望有高手指点一二！不胜感激！

以下是输入文件INCAR：

NWRITE = 2 # verbosity write-tag(how much is written)
PREC = Normal
ISTART = 0 # startjob 0=new 1=cont 2=samecut
ISPIN = 2 # spin polarized calculation (2=yes,1=no)
INIWAV = 1 # initial electr wf 0=lowe 1=rand
ENCUT = 450.0 eV
EDIFF = 1.0E-06
NSW = 0 # number of steps for ionic upd
IBRION = -1 # ionic relaxation 0=MD 1=quasi-New 2=CG
ISIF = 2
ISMEAR = 0 # partial occupancies -5=Blochl(energy) -4=tet -1=fermi 0=gaussian >0 MP(force)
SIGMA = 0.05
LORBIT = 11 # 0 simple, 1 ext, 2 COOP (PROOUT)
# NELECT = 702.0000
#plus U

LDAU=.TRUE.
LDAUTYPE=2
LDAUL= -1 -1 -1 2
LDAUU=0.00 0.00 0.00 5.2
LDAUPRINT=2
LMAXMIX=4

[ Last edited by utimasa on 2012-2-12 at 21:31 ]

回复此楼

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

» 猜你喜欢

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

1楼 2012-02-12 21:29:50

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

WDD880227

木虫 (正式写手)

1ST强帖: 1
应助: 25 (小学生)
贵宾: 0.098
金币: 1348.5
散金: 194
红花: 16
帖子: 761
在线: 371.6小时
虫号: 1129580
注册: 2010-10-22
性别: MM
专业: 凝聚态物性 II ：电子结构

引用回帖:

6楼: Originally posted by acridine at 2012-02-15 17:46:46:
两个体系电子数目不等，直接比较总能量是没有意义的。
如果真的要比带电前后的能量差别，那也应该是比较：（"中性体系+一个额外电子"的总能量)与（"带一个电子体系"的总能量)
后一项自然是 ...

请问band align您是怎么做的呢？