24小时热门版块排行榜    

查看: 2950  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 一志愿武汉理工材料工程专硕调剂 +5 Doleres 2026-03-19 5/250 2026-03-19 20:14 by 制度的
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
信息提示
请填处理意见