24小时热门版块排行榜    

查看: 2308  |  回复: 11

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

» 猜你喜欢

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

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

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


uuv2010(金币+1): 多谢提示 2011-08-21 00:20:55
first_yg(金币+2): 谢谢 2011-08-22 09:13:54
把两者的ismear调到一样的试试
嘟啦啦嘟啦啦嘟
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.这样试一试~
人有多大胆,地有多大产!
4楼2011-08-22 09:23:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tider

木虫 (正式写手)

这个问题我也遇到过呢,LZ算出结构来可否分享一下
5楼2011-08-22 23:43:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)

引用回帖:
5楼: Originally posted by tider at 2011-08-22 23:43:50:
这个问题我也遇到过呢,LZ算出结构来可否分享一下

好的,等我算对了把经验贴上来。这两天空调坏了,机器关了,唉。
6楼2011-08-23 10:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)

引用回帖:
4楼: Originally posted by haowenping at 2011-08-22 09:23:22:
你bands里面的ISMEAR应该是取的0吧,还有SIGMA是默认的0.2,这个0.2可能有点大,你可以这样试一试,就是重新做一次计算,在弛豫计算的时候,设置,ISMEAR=0,SIGMA=0.01,SCF计算的时候,设置ISMEAR=0,SIGMA=0. ...

按照这种方法算了一下,还是不行,band 算出来gap依然远大于dos
7楼2011-08-26 14:59:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

★ ★
gzqdyouxia(金币+2): 谢谢提示 2011-09-09 15:47:51
引用回帖:
7楼: Originally posted by first_yg at 2011-08-26 14:59:36:
按照这种方法算了一下,还是不行,band 算出来gap依然远大于dos

你计算bands时候,K点怎么选取的?DOS的时候呢?K点又是怎么选取的?如果你在算DOS的时候,选用ISMEAR=-5,K点选取为在计算bands的时候的3倍再试试。
人有多大胆,地有多大产!
8楼2011-08-26 16:10:03
已阅   回复此楼   关注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的回帖

utdwang

银虫 (初入文坛)

★ ★
gzqdyouxia(金币+1): 谢谢指教 2011-09-09 15:46:29
fzx2008(金币+1): 谢谢指教 2011-12-11 19:54:13
It could also be possible that your K-point in band calculation is not complete. For instance, there might be one k-point missing which produce the smaller band gap. FYR.
Thekeytosuccessinnearlyeverythingisdailyaction
10楼2011-09-09 07:52:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 first_yg 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 怎么带研究生? +22 豆豆小小2 2024-06-20 24/1200 2024-06-26 09:10 by Xiaolin81
[考研] 刚当完兵回来想考研 +5 五官在线 2024-06-23 16/800 2024-06-26 09:01 by Xiaolin81
[硕博家园] 申博雅思求助 +6 刘嘉研究究 2024-06-24 7/350 2024-06-26 07:08 by 小李醒yy
[基金申请] 基金申请书名称有变化 +5 xuel2011 2024-06-25 7/350 2024-06-26 00:07 by 老虎当猫养
[有机交流] 对苯乙烯磺酰氯的合成机理 25+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[考博] 一名额难求啊 +14 pinbo拼搏 2024-06-19 14/700 2024-06-25 16:55 by wjykycg
[基金申请] 今天能不能出来名单 +7 地理学1995 2024-06-25 9/450 2024-06-25 15:20 by lyfbangong
[教师之家] 复旦夏同学退学理由说明,哪儿可以下载? +7 苏东坡二世 2024-06-21 8/400 2024-06-25 12:42 by 药嘿1233
[硕博家园] 数据不好 +5 Hetai 2024-06-23 7/350 2024-06-25 12:37 by 1591099
[基金申请] 今年什么时候会评啊 +7 lancet0903 2024-06-24 7/350 2024-06-24 18:48 by Pickfoot
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 国自然资助比率是不是要下降了?? +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[基金申请] 青年和面上,哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[有机交流] 三颈瓶加热 5+4 wzjwx888 2024-06-20 4/200 2024-06-24 10:11 by w87437985
[公派出国] 博士csc联培会看重第一学历学校层次吗 +4 也就这样 2024-06-23 4/200 2024-06-24 08:18 by 晓目崇
[基金申请] 江苏省333人才工程 出校后被刷的概率大不大? +8 maxbirdzhang 2024-06-19 12/600 2024-06-23 20:39 by 查拉图斯特拉
[基金申请] 工材口青年基金大概什么样能上会? +15 今晚推荐22 2024-06-20 21/1050 2024-06-22 23:04 by qbn0326
[基金申请] 教育部基金 +5 m1393 2024-06-21 5/250 2024-06-21 21:13 by odes
[基金申请] 面上博士后 +11 jsqy 2024-06-19 13/650 2024-06-21 17:12 by sizhouyi
[论文投稿] ACS AMI 返回审稿意见,一个大修,两个据稿,编辑给的修改重投 +7 智商已更新 2024-06-19 7/350 2024-06-20 19:54 by kinlin13
信息提示
请填处理意见