24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1303  |  回复: 2
【悬赏金币】回答本帖问题,作者happyzyf将赠送您 20 个金币

happyzyf

铜虫 (小有名气)

[求助] vasp手册GW计算Si带隙有问题已有1人参与

各位高手,

我最近根据vasp手册里面的GW计算Si的带隙的例子(http://cms.mpi.univie.ac.at/wiki/index.php/Bandgap_of_Si_in_GW)做计算,结果发现我计算出来的Si的HOMO=156.8865 eV,而LUMO=142.2128 eV,OUTCAR中费米能级的位置在144.9269 eV。

[iklhrd]$ ./gap_GW.sh OUTCAR
HOMO: band: 4  E= 156.8865
LUMO: band: 5  E= 142.2128


这里HOMO比LUMO要高,这样的话,说明是金属性质啊 我是用的例子里面的输入文件,到底是哪里的问题呢?

三步的INCAR分别是

1groundstate.dft

ISMEAR =  0
SIGMA  =  0.05
EDIFF  = 1E-8
ALGO   = Very_Fast

2DFT.virtual.orbitals

# Frequency dependent dielectric tensor without
# local field effects
ALGO = Exact
NBANDS  = 64
LOPTICS = .TRUE. ; CSHIFT = 0.1
NEDOS = 2000
# you might try
#LPEAD = .TRUE.

ISMEAR =  0
SIGMA  =  0.05
EDIFF  = 1E-8

3GW

# Frequency dependent dielectric tensor including
# local field effects within the RPA (default) or
# including changes in the DFT xc-potential (LRPA=.FALSE.).
# N.B.: beware one first has to have done a
# calculation with ALGO=Exact, LOPTICS=.TRUE.
# and a reasonable number of virtual states (see above)
ALGO = GW0 ; LSPECTRAL = .TRUE. ; NOMEGA = 50

# be sure to take the same number of bands as for
# the LOPTICS=.TRUE. calculation, otherwise the
# WAVEDER file is not read correctly
NBANDS = 64

# Add this to update the quasiparticle energies
# in the Green's function (GW0)
#NELM = 4

ISMEAR =  0
SIGMA  =  0.05
EDIFF  = 1E-8
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taichizsf

新虫 (初入文坛)

算错了吧,可以先看看DOS

发自小木虫Android客户端
2楼2017-10-22 19:57:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

IanTsian

新虫 (初入文坛)

【答案】应助回帖

应该是你用了官网上算例自带的gap.sh和gap_GW.sh吧,那两个程序求LUMO和HOMO本身就是有问题的,还是自己写个程序抽取LUMO和HOMO吧。
3楼2018-12-04 21:37:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 happyzyf 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 24计算机申博 +4 学无止境er 2024-04-13 6/300 2024-04-16 19:15 by 学无止境er
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 13/650 2024-04-16 18:02 by 鱼翔浅底1
[教师之家] 北邮女博导并没有强制学生必须参与项目,说了学生可以不再参与项目 +38 苏东坡二世 2024-04-10 59/2950 2024-04-16 17:05 by dujianchao
[考研] 345 求调剂 +4 南风轻语5 2024-04-12 6/300 2024-04-16 16:36 by chenweiwade
[找工作] 合肥新成立市属本 VS 安徽老牌省属本 50+3 masterofog 2024-04-11 14/700 2024-04-16 15:48 by 塌棵菜
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[考研] 一志愿南京理工大学材料工程282求调剂 +13 abc16785 2024-04-11 13/650 2024-04-15 11:26 by coolfishwll
[晶体] 单晶仪在50kv 30mA ,电压电流没降,能开门上样吗? +3 15665390996 2024-04-10 3/150 2024-04-15 07:36 by gyliu
[考研] 367求调剂 +8 gjf峰ab 2024-04-10 9/450 2024-04-14 22:42 by 永字号
[考研] 化工专硕求调剂 +9 青梅duoduo 2024-04-12 11/550 2024-04-13 09:33 by lincunhui
[考研] 生物与医药调剂 +8 祁佳慧 2024-04-10 8/400 2024-04-12 22:11 by achao9525
[考研] 北信科第二批调剂,工科相关 +4 biser 科研人 2024-04-11 6/300 2024-04-12 21:26 by biser 科研人
[考研] 284求调剂 +5 Luvkyggggg 2024-04-11 6/300 2024-04-12 21:11 by coco1981
[考研] 312求调剂 +16 dhdhjdhd 2024-04-10 20/1000 2024-04-11 21:34 by kaiyanjie
[考研] 考研调剂 +5 夏末的晨曦12 2024-04-10 6/300 2024-04-11 18:51 by 17814699553
[考研] 327求调剂 +5 某十七 2024-04-11 5/250 2024-04-11 18:48 by shangxh
[考研] 321求调剂 +7 MengdieW 2024-04-10 8/400 2024-04-11 00:18 by xuxiaodai
[考研] 312求调剂 +9 以前老狠啦 2024-04-10 10/500 2024-04-10 22:42 by 15081537629
[考研] 求调剂 +3 LED绿光灯 2024-04-10 6/300 2024-04-10 17:22 by LED绿光灯
[考研] 求调剂 +4 吃了一口德芙 2024-04-10 4/200 2024-04-10 13:48 by 缘子你好
信息提示
请填处理意见