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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 虫友帮忙！band和dos计算得到的gap值差别很大，无法对应！
51/1返回列表
查看: 3002  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

first_yg

银虫 (小有名气)

[求助] 虫友帮忙！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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-08-20 17:18:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)
引用回帖:
8楼: Originally posted by haowenping at 2011-08-26 16:10:03:
你计算bands时候，K点怎么选取的？DOS的时候呢？K点又是怎么选取的？如果你在算DOS的时候，选用ISMEAR=-5，K点选取为在计算bands的时候的3倍再试试。

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-08-21 00:20:55
first_yg(金币+2): 谢谢 2011-08-22 09:13:54
把两者的ismear调到一样的试试
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
2楼2011-08-20 21:48:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)
引用回帖:
2楼: Originally posted by 贺仪 at 2011-08-20 21:48:29:
把两者的ismear调到一样的试试

但是DOS一般不是说用-5比较准吗？
请问你一般band和dos都设置什么呢？
赞一下 回复此楼
3楼2011-08-22 09:14:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

【答案】应助回帖

★ ★
first_yg(金币+3): 多谢你这么详细的回帖，我试一试！ 2011-08-22 09:34:27
youzhizhe(金币+2): 谢谢交流。 2011-08-22 20:37:30
引用回帖:
3楼: Originally posted by first_yg at 2011-08-22 09:14:35:
但是DOS一般不是说用-5比较准吗？
请问你一般band和dos都设置什么呢？

你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.这样试一试~
赞一下(1人) 回复此楼
人有多大胆，地有多大产！
4楼2011-08-22 09:23:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[教师之家] 山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐 +4 qut2026 2026-04-11 8/400 2026-04-17 16:10 by 会飞的猪157
[考研] 一志愿沪9，326求生物学调剂 +11 刘墨墨 2026-04-13 11/550 2026-04-17 15:35 by wanganpin
[考研] 291求调剂 +9 关忆北. 2026-04-14 9/450 2026-04-16 22:49 by cfdbai
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 一志愿A区211，22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 通信工程求调剂！！！ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 药学求调剂 +11 RussHu 2026-04-12 13/650 2026-04-15 19:07 by zhuwenxu
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 材料工程281还有调剂机会吗 +43 xaw. 2026-04-11 44/2200 2026-04-15 12:46 by 西北望—风沙
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 考研求调剂 +6 ban班小七 2026-04-11 6/300 2026-04-14 14:06 by 哆啦A梦只是个梦
[考研] 食品与营养（0955）271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
信息提示
关闭
请填处理意见
关闭