版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zhoudayan2

木虫 (小有名气)

sherry

应助: 1 (幼儿园)
金币: 1550.3
帖子: 109
在线: 67小时
虫号: 1098641
注册: 2010-09-14
性别: MM
专业: 能源化工

[求助] 能带计算过程

第一次算能带，过程如下，欢迎指正！！首先对体系构型优化，INCAR如下：
SYSTEM = default # Selfconsistent calculation is firstly performed.

Start parameter for this Run:
  ISTART = 0
  ICHARG = 2
  INIWAV = 1

Electronic Relaxation:
  ENCUT = 400
  ENAUG = 650
  PREC = Accurate
  IALGO = 38
  NELM = 600
  NELMIN = 8 # surface = 8, bulk = 4
#  NELMDL = -12 # When IALGO = 38, 0
  EDIFF = 1E-04 # If high accuracy is required, set 1E-06.
#  NBANDS = 7
  GGA = PE
#  VOSKOWN = 1
  LREAL = Auto
  WEIMIN = 0

Ionic Relaxation:
  EDIFFG = -0.03
  NSW = 1000
  IBRION = 3
  ISIF = 2
  POTIM = 0.1
#  IOPT = 2
#  ISYM = 0

DOS related values:
  SIGMA = 0.05
  ISMEAR = 0

Spin polarized:
ISPIN = 2
  MAGMOM =13*4

File writing
  LWAVE = .FALSE.
  LCHARG = .FALSE.

Calculation of frequency:
#  POTIM = 0.02
#  NFREE = 2
NPAR=1
优化好后，将CONTCAR copy 成POSCAR进行静态计算，INCAR如下：
SYSTEM = default # Selfconsistent calculation is firstly performed.

Start parameter for this Run:
  ISTART = 0
  ICHARG = 2
  INIWAV = 1

Electronic Relaxation:
  ENCUT = 400
  ENAUG = 650
  PREC = Accurate
  IALGO = 38
  NELM = 600
  NELMIN = 8 # surface = 8, bulk = 4
#  NELMDL = -12 # When IALGO = 38, 0
  EDIFF = 1E-07 # If high accuracy is required, set 1E-06.
#  NBANDS = 7
  GGA = PE
#  VOSKOWN = 1
  LREAL = Auto
  WEIMIN = 0

Ionic Relaxation:
#  EDIFFG = -0.03
  NSW = 1
  IBRION = -1
#  ISIF = 2
#  POTIM = 0.1
#  IOPT = 2
#  ISYM = 0

DOS related values:
  SIGMA = 0.05
  ISMEAR = 0

Spin polarized:
ISPIN = 2
  MAGMOM =13*4

File writing
  LWAVE = .FALSE.
  LCHARG = .TRUE.

Calculation of frequency:
#  POTIM = 0.02
#  NFREE = 2
NPAR=1
得到的CHG,与CHGCAR 在copy到算能带中，同时静态计算中的4个输入文件，INCAR 如下：
SYSTEM = default # Selfconsistent calculation is firstly performed.

Start parameter for this Run:
  ISTART = 1
  ICHARG = 11
  INIWAV = 1

Electronic Relaxation:
  ENCUT = 400
  ENAUG = 650
  PREC = Accurate
  IALGO = 38
  NELM = 600
  NELMIN = 8 # surface = 8, bulk = 4
#  NELMDL = -12 # When IALGO = 38, 0
  EDIFF = 1E-07 # If high accuracy is required, set 1E-06.
NBANDS = 100
  GGA = PE
#  VOSKOWN = 1
  LREAL = Auto
  WEIMIN = 0

Ionic Relaxation:
#  EDIFFG = -0.03
  NSW = 1
  IBRION = -1
#  ISIF = 2
#  POTIM = 0.1
#  IOPT = 2
#  ISYM = 0

DOS related values:
  SIGMA = 0.05
  ISMEAR = 0

Spin polarized:
ISPIN = 2
  MAGMOM =13*4

File writing
  LWAVE = .FALSE.
  LCHARG = .TRUE.

Calculation of frequency:
#  POTIM = 0.02
#  NFREE = 2
NPAR=1
LORBIT=11，
这样算对吗，步骤我是按照侯老师手册上的。
这三步，的KPOINTS 都没变，是：
K-Points
0
Gamma
1  1  1
0 0 0.
欢迎指正，然后怎么画能带图啊？HOMO-LUMO gap 这些从哪个文件得到啊？

回复此楼

» 猜你喜欢

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

VASP+SOC后的能带计算已经有21人回复
最新计算能带的方法得出的带隙和实验非常的吻合。已经有9人回复
HSE（GW类似）计算能带总结已经有51人回复
求石墨烯能带计算出现狄拉克点的K点路径已经有13人回复
【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？已经有13人回复
【求助】求VASP计算ZnO能带和态密度的详细输入设置已经有15人回复
【其他】关于GW（准粒子计算能带结构等）方法的工作已经有6人回复
【求助】PBE0 计算Si、Ge能带结构已经有9人回复
【求助】初学者求助：紧束缚方法计算能带结构已经有5人回复

人生就像河流，无法回头...

1楼 2012-03-21 15:15:16

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

贺仪

铁杆木虫 (著名写手)

1ST强帖: 5
应助: 93 (初中生)
贵宾: 0.636
金币: 5937
散金: 1063
红花: 38
沙发: 2
帖子: 2499
在线: 588.7小时
虫号: 995125
注册: 2010-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★
zzy870720z: 金币+2, 谢谢指教 2012-03-21 16:20:27

这三步，的KPOINTS 都没变，

不变不行
算能带这个KPOINTS是要根据晶体类型来写的

赞一下

回复此楼

嘟啦啦嘟啦啦嘟

2楼2012-03-21 16:16:35

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

zhoudayan2

木虫 (小有名气)

sherry

应助: 1 (幼儿园)
金币: 1550.3
帖子: 109
在线: 67小时
虫号: 1098641
注册: 2010-09-14
性别: MM
专业: 能源化工

引用回帖:

2楼: Originally posted by 贺仪 at 2012-03-21 16:16:35:
这三步，的KPOINTS 都没变，

不变不行
算能带这个KPOINTS是要根据晶体类型来写的

嗯，在第三步算能带时需要改变K点，但如何在kpoints文件中设置啊，侯老师写的没太看明白，能否直接用vasp手册中的：k-points along high symmetry lines
10  ! 10 intersections
Line-mode
cart
  0 0 0 ! gamma
  0 0 1 ! X

  0 0 1 ! X
  0.5 0 1 ! W

  0.5 0 1 ! W
  0 0 1 ! gamma

谢谢啊~实在是没弄懂~

赞一下

回复此楼

人生就像河流，无法回头...

3楼2012-03-21 20:44:13

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

贺仪

铁杆木虫 (著名写手)

1ST强帖: 5
应助: 93 (初中生)
贵宾: 0.636
金币: 5937
散金: 1063
红花: 38
沙发: 2
帖子: 2499
在线: 588.7小时
虫号: 995125
注册: 2010-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
感谢参与，应助指数 +1
WDD880227: 金币+1, 感谢提示 2012-03-22 10:37:02

引用回帖:

3楼: Originally posted by zhoudayan2 at 2012-03-21 20:44:13:
嗯，在第三步算能带时需要改变K点，但如何在kpoints文件中设置啊，侯老师写的没太看明白，能否直接用vasp手册中的：k-points along high symmetry lines
10 ! 10 intersections
Line-mode
cart
0 0 ...

这个根据你算体系的具体对称性来决定的

赞一下

回复此楼

嘟啦啦嘟啦啦嘟

4楼2012-03-21 22:23:00

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