24小时热门版块排行榜

返回列表

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

navyzhang

木虫 (正式写手)

应助: 20 (小学生)
金币: 1411.8
散金: 1197
红花: 5
帖子: 640
在线: 187.1小时
虫号: 883307
注册: 2009-10-25
性别: GG
专业: 理论和计算化学

[求助] GW计算，在哪一步找费米能级？已有1人参与

各位专家：
我最近在测试，wennier插值法进行GW计算，已经按照网上的教程跑完全部步骤
http://blog.sciencenet.cn/blog-567091-802483.html

最后画出来的能带跟网上晒出来的结果也是非常相近的（可能因为我的计算精度低些，能带有些许不一样），具体两个能带见下面附图。

我想，两个能带的区别在于费米能级有没有放在 0 eV位置吧？

那么这个网上GW方法计算Si的例子里面，费米能级到底是在哪一步？哪个文件里面找到的？我是在第三步（ALGO=GW0）的OUTCAR文件中找不到Fermi level.

还望各位牛人不吝赐教！在下感激不尽！！

前几步的VASP计算对应INCAR一并贴在这里吧

step1
ISTART = 0
ICHARG = 2

IBRION = -1
NSW = 0

   GGA = PE
   ENCUT = 600

PREC = Accurate
NELM = 100
NELMIN = 2
EDIFF  = 1.0e-5
ISMEAR =  0
SIGMA  = 0.05
LWAVE = .TRUE.
  LCHARG = .TRUE.
ALGO  = NORMAL

LREAL = .False.

step2
ISTART = 1

IBRION = -1
NSW = 0

   GGA = PE
   ENCUT = 600

PREC = Accurate
NELM = 100
NELMIN = 2
EDIFF  = 1.0e-5
ISMEAR =  0
SIGMA  = 0.05
LWAVE = .TRUE.
  LCHARG = .TRUE.

ALGO  = Exact
NBANDS = 64
LOPTICS = .TRUE.
NEDOS = 2000

LREAL = .False.

step3

ISTART = 1
IBRION = -1
NSW = 0

   GGA = PE
   ENCUT = 600

PREC = Accurate
NELM = 100
NELMIN = 2
EDIFF  = 1.0e-5
ISMEAR =  0
SIGMA  = 0.05
LWAVE = .TRUE.
  LCHARG = .TRUE.

ALGO = GW0
LSPECTRAL = .TRUE.
NOMEGA = 50
NBANDS = 64
LRPA = .FALSE.
LWANNIER90 = .TRUE.

# ALGO  = Exact
# NBANDS = 64
# LOPTICS = .TRUE.
# NEDOS = 2000

LREAL = .False.

未标题-1.png

回复此楼

» 猜你喜欢

1楼 2016-01-18 06:00:56

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

对马拉二胡

至尊木虫 (职业作家)

应助: 27 (小学生)
金币: 25552.8
散金: 244
红花: 13
帖子: 4006
在线: 937.9小时
虫号: 2471734
注册: 2013-05-19
性别: GG
专业: 无机非金属材料

楼主，你好。我按着官网例子计算Si时，在运行wannier90.x wannier90.win 这步出现了
Wannier90: Execution started on 21Jan2016 at 14:42:15
Exiting.......
param_get_projections: Fewer projections defined than the number of Wannier functions requested
我想请教一下这是什么原因呢？谢谢

赞一下

回复此楼

相互学习，共同进步。

4楼2016-01-21 14:43:31

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

查看全部 8 个回答

navyzhang

木虫 (正式写手)

应助: 20 (小学生)
金币: 1411.8
散金: 1197
红花: 5
帖子: 640
在线: 187.1小时
虫号: 883307
注册: 2009-10-25
性别: GG
专业: 理论和计算化学

不能让帖子沉下去

赞一下

回复此楼

2楼2016-01-18 20:53:56

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

navyzhang

木虫 (正式写手)

应助: 20 (小学生)
金币: 1411.8
散金: 1197
红花: 5
帖子: 640
在线: 187.1小时
虫号: 883307
注册: 2009-10-25
性别: GG
专业: 理论和计算化学

引用回帖:

4楼: Originally posted by 对马拉二胡 at 2016-01-21 14:43:31
楼主，你好。我按着官网例子计算Si时，在运行wannier90.x wannier90.win 这步出现了
Wannier90: Execution started on 21Jan2016 at 14:42:15
Exiting.......
param_get_projections: Fewer projections defi ...

wannier90.win文件能贴出来看一下吗？

赞一下

回复此楼

5楼2016-01-22 01:47:43

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