| 查看: 1075 | 回复: 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 这些从哪个文件得到啊？ |
回复此楼» 猜你喜欢
中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培
已经有4人回复
-
调剂推荐
已经有4人回复
-
总分322求生物学/生化与分子/生物信息学相关调剂
已经有5人回复
-
0703化学338求调剂！
已经有3人回复
-
284求调剂
已经有7人回复
-
271求调剂
已经有6人回复
-
312求调剂
已经有9人回复
-
352求调剂
已经有4人回复
-
085602化学工程求调剂。
已经有3人回复
-
085600材料与化工306
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP+SOC后的能带计算
已经有21人回复
- 最新计算能带的方法 得出的带隙和实验非常的吻合。
已经有9人回复
- HSE（GW类似）计算能带总结
已经有51人回复
- 求石墨烯能带计算出现狄拉克点的K点路径
已经有13人回复
- 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
已经有13人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 【其他】关于GW（准粒子计算能带结构等）方法的工作
已经有6人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【求助】初学者求助： 紧束缚方法计算能带结构
已经有5人回复
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
4楼2012-03-21 22:23:00
thinkllh
铁杆木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 6569.1
- 散金: 50
- 红花: 1
- 帖子: 538
- 在线: 157小时
- 虫号: 715800
- 注册: 2009-03-05
- 性别: GG
- 专业: 极端条件下使用的金属材料
5楼2012-03-22 09:46:37
