| 查看: 3431 | 回复: 11 | |||
[求助]
虫友帮忙！band和dos计算得到的gap值差别很大，无法对应！
|
|
计算过很多体系的Band和dos，有时能很好的对应，但有的体系怎么算都对不上，非常着急，困扰我很久了。请虫友们帮忙！ 对于这个对应不上的体系，band计算时ismear =0 和默认值都测试过，但是无法用-5. dos计算一直采用的是-5.另外band计算时高对称点的路径也考查了多种，都不行。 下面是我的INCAR : 能带计算的： NWRITE = 2; #LPETIM=.TRUE. ISTART = 1 job : 0-new 1-cont 2-samecut LCHARG = F LWAVE = F ISPIN = 1 ISYM = 2 PREC = Accurate ENCUT = 400 ICHARG = 11 Electronic Relaxation NELM = 60; EDIFF = 1E-05 stopping-criterion for ELM LREAL = A Ionic Relaxation # IBRION = -1 ionic relax: 0-MD 1-quasi-New 2-CG ISIF = 2 relax z axiz NSW = 0 EDIFFG = -0.01 POTIM = 0.50 DOS related values: #ISMEAR = 0 #SIGMA = 0.1 Algorithm: #IALGO = 48 这个是DOS的： LORBIT=11 NEDOS=1001 ICHARG =11 NSW = 1 IBRION = -1 POTIM = 0.05 ENCUT = 400.000000 #SIGMA = 0.100000 PREC = Accurate ISMEAR = -5 EMIN = -10 EMAX = 10 ADDGRID = .TRUE. LPLANE = .TRUE. LCHARG = .FALSE. LWAVE =.FALSE. ISPIN = 2  [ Last edited by first_yg on 2011-8-20 at 17:19 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 lywiailyw |
» 猜你喜欢
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算能带，刚开始就服务器就卡了
已经有20人回复
- 蛋白酶解后 得到的肽 比蛋白还多 求高人解释一下
已经有12人回复
- 金簇浓度的计算
已经有6人回复
- 第一性原理形成能计算
已经有6人回复
- 纯理论计算投稿的问题
已经有21人回复
- 请问发电机中的sideband current如何翻译
已经有1人回复
- vasp 计算表面DOS与文献对应不上 请求大家帮助.......
已经有20人回复
- 【求助】理论计算新手，求帮忙解决带隙及DOS的问题
已经有21人回复
- 【求助成功】supercell中DOS计算如何得到合理的电荷密度和波函数？
已经有3人回复
- 【求助】计算dos时很慢，而且WAVECAR很大
已经有21人回复
2楼2011-08-20 21:48:29
3楼2011-08-22 09:14:35
haowenping
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 2330.2
- 散金: 384
- 红花: 5
- 帖子: 1104
- 在线: 244.5小时
- 虫号: 1124850
- 注册: 2010-10-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★
first_yg(金币+3): 多谢你这么详细的回帖，我试一试！ 2011-08-22 09:34:27
youzhizhe(金币+2): 谢谢交流。 2011-08-22 20:37:30
first_yg(金币+3): 多谢你这么详细的回帖，我试一试！ 2011-08-22 09:34:27
youzhizhe(金币+2): 谢谢交流。 2011-08-22 20:37:30
|
你bands里面的ISMEAR应该是取的0吧，还有SIGMA是默认的0.2，这个0.2可能有点大，你可以这样试一试，就是重新做一次计算，在弛豫计算的时候，设置，ISMEAR=0，SIGMA=0.01，SCF计算的时候，设置ISMEAR=0，SIGMA=0.01，bands也是ISMEAR=0，SIGMA=0.01，dos计算分别用（1）、ISMEAR=0，SIGMA=0.01，（2）、ISMEAR=-5.这样试一试~ |
4楼2011-08-22 09:23:22
tider
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1950.4
- 散金: 40
- 红花: 5
- 帖子: 313
- 在线: 87.2小时
- 虫号: 1190614
- 注册: 2011-01-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2011-08-22 23:43:50
6楼2011-08-23 10:48:06
7楼2011-08-26 14:59:36
haowenping
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 2330.2
- 散金: 384
- 红花: 5
- 帖子: 1104
- 在线: 244.5小时
- 虫号: 1124850
- 注册: 2010-10-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2011-08-26 16:10:03
|
band的KPOINT kpoints for bandstructure G-M2-K2-M1-G-K1-M1 41 line reciprocal 0.000000000000000 0.000000000000000 0.000000000000000 1 ! G 0.500000000000000 0.500000000000000 0.000000000000000 1 ! M2 0.500000000000000 0.500000000000000 0.000000000000000 1 ! M2 0.333333333333333 0.666666666666666 0.000000000000000 1 ! K2 0.333333333333333 0.666666666666666 0.000000000000000 1 ! K2 0.500000000000000 0.000000000000000 0.000000000000000 1 ! M1 0.500000000000000 0.000000000000000 0.000000000000000 1 ! M1 0.000000000000000 0.000000000000000 0.000000000000000 1 ! G 0.000000000000000 0.000000000000000 0.000000000000000 1 ! G 0.666666666666666 0.333333333333333 0.000000000000000 1 ! K1 0.666666666666666 0.333333333333333 0.000000000000000 1 ! K1 0.500000000000000 0.000000000000000 0.000000000000000 1 ! M1 dos时候的： 0 Monkhorst-Pack 21 21 1 0 0 0 |
9楼2011-08-26 16:50:11
10楼2011-09-09 07:52:49
20