24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2980  |  回复: 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 ]
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料334求调剂 +21 Eecho# 2026-04-03 21/1050 2026-04-08 22:55 by 猪会飞
[考研] 285求调剂 +7 AZMK 2026-04-07 7/350 2026-04-08 11:10 by 逆水乘风
[考研] 274求调剂 +4 山阿蔓 2026-04-07 4/200 2026-04-08 09:49 by eyeliyang
[考研] 288求调剂 +12 没有答案_ 2026-04-05 12/600 2026-04-08 00:17 by T可可西里T
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] (调剂)一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 8/400 2026-04-06 16:52 by momo皓
[考研] 求调剂 +11 xzghyuj 2026-04-04 11/550 2026-04-06 11:49 by lijunpoly
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 26考研调剂0710 0860 +9 补补不补 2026-04-03 14/700 2026-04-04 23:32 by 果冻大王
[考研] 298求调剂 +5 zzz,,r 2026-04-02 8/400 2026-04-04 19:55 by 蓝云思雨
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
信息提示
请填处理意见