| 查看: 1079 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
zhoudayan2木虫 (小有名气)
sherry
|
[求助]
能带计算过程
|
|
|
第一次算能带，过程如下，欢迎指正！！首先对体系构型优化，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 这些从哪个文件得到啊？ |
回复此楼» 猜你喜欢
调剂求收留
已经有7人回复
-
317求调剂
已经有6人回复
-
336材料求调剂
已经有9人回复
-
341求调剂
已经有7人回复
-
286求调剂
已经有5人回复
-
化学308分求调剂
已经有5人回复
-
0703化学338求调剂！
已经有6人回复
-
求调剂推荐 材料 304
已经有3人回复
-
材料学硕333求调剂
已经有6人回复
-
一志愿郑州大学，080500学硕，总分317分求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP+SOC后的能带计算
已经有21人回复
- 最新计算能带的方法 得出的带隙和实验非常的吻合。
已经有9人回复
- HSE（GW类似）计算能带总结
已经有51人回复
- 求石墨烯能带计算出现狄拉克点的K点路径
已经有13人回复
- 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
已经有13人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 【其他】关于GW（准粒子计算能带结构等）方法的工作
已经有6人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【求助】初学者求助： 紧束缚方法计算能带结构
已经有5人回复
4楼2012-03-21 22:23:00
2楼2012-03-21 16:16:35
zhoudayan2
木虫 (小有名气)
sherry
- 应助: 1 (幼儿园)
- 金币: 1550.3
- 帖子: 109
- 在线: 67小时
- 虫号: 1098641
- 注册: 2010-09-14
- 性别: MM
- 专业: 能源化工
3楼2012-03-21 20:44:13
