| 查看: 1745 | 回复: 2 | ||
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 |
回复此楼» 猜你喜欢
基金未中，这种答复是模板吗？
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
-
面上合作单位盖章
已经有3人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有5人回复
-
基金不中，共勉
已经有7人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
-
让我中一个面上吧！
已经有17人回复
taichizsf
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 6.8
- 帖子: 43
- 在线: 64.7小时
- 虫号: 3450928
- 注册: 2014-09-30
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2017-10-22 19:57:38
3楼2018-12-04 21:37:48
投票:
10