24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 5832  |  回复: 20

first_yg

银虫 (小有名气)

[交流] 【求助成功】VASP算出来的能带和DOS图对应不上,哪里出了问题? 已有13人参与

我画Band的INCAR:
NWRITE =      2;
  LPETIM=.TRUE.                                          
  ISTART =      0    job   : 0-new  1-cont  2-samecut                                             
  LCHARG =      F
  LWAVE  =      F                                                                                 
  ISPIN  =      1                                                                                 
  ISYM   =      2                                                                                 
  PREC   =      Normal                                                                           
  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.05                                                                              
                                                                                                  
Algorithm:                                                                                       
  #IALGO  =     48

画DOS的脚本:
NWRITE =      2
LORBIT=11
NEDOS=1001
EMIN = -20
EMAX = 20
ICHARG = 11
NSW     = 1      
IBRION  = -1         
POTIM   = 0.05
ENCUT = 850.000000
SIGMA = 0.100000
PREC = Accurate
ISMEAR = -5
ADDGRID = .TRUE.
LPLANE = .TRUE.
ISPIN = 1

最后图形,差了很多



[ Last edited by first_yg on 2010-12-8 at 14:35 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 模拟,理论

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

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

first_yg

银虫 (小有名气)

哪位高手能帮忙解答一下啊
2楼2010-12-03 12:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondstar

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-12-03 14:08:08
能带是沿着高对称点走的,DOS是全空间的,应该以DOS的为准
3楼2010-12-03 13:34:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)

引用回帖:
Originally posted by beyondstar at 2010-12-03 13:34:35:
能带是沿着高对称点走的,DOS是全空间的,应该以DOS的为准

可是在我所看过的文献里,别人的能带结构和dos都是完全吻合的啊
4楼2010-12-03 14:55:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢指教 2010-12-03 22:33:55
很明显, 这是由于你在做能带图的时候,减去的费米能不对造成的。
真理是一点点接近的!
5楼2010-12-03 21:52:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaogh

银虫 (初入文坛)


youzhizhe(金币+1): 谢谢交流。 2011-03-21 10:05:00
引用回帖:
Originally posted by gavinliu7390 at 2010-12-03 21:52:41:
很明显, 这是由于你在做能带图的时候,减去的费米能不对造成的。

正解
6楼2010-12-03 23:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangguiyang

新虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang1988(金币+1):平移要有根据,不能为对齐而平移吧 2010-12-04 00:25:43
自己平移一下,不就能对齐了吗?
7楼2010-12-04 00:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)

楼上各位,是这样:
其实左图和右图为了归零减去的fermi能级值是不同的。左图计算能带的时候,我用的是ISMEAR=0,而右图计算dos的时候用的是-5。因为考虑到我的体系结构优化用0比较好,但是DOS又必须用-5.这样算出来的fermi能级数值是不同的。
我现在不知道fermi能该取哪个了

[ Last edited by first_yg on 2010-12-6 at 13:10 ]
8楼2010-12-06 10:53:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyg3395

金虫 (小有名气)

★ ★
zzy870720z(金币+2):谢谢指教 2010-12-06 12:02:42
引用回帖:
Originally posted by first_yg at 2010-12-06 10:53:44:
楼上各位,是这样:
其实左图和右图为了归零减去的fermi能级值是不同的。左图计算能带的时候,我用的是ISMEAR=0,而右图计算dos的时候用的是-5。因为考虑到我的体系结构优化用0比较好,但是DOS又必须用-5.这样算 ...

ISMEAR=0是K点数目较少的时候用的,一般情况下,计算band时用系统默认的就可以了。静态计算ISMEAR是要设置为-5的
9楼2010-12-06 11:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

first_yg

银虫 (小有名气)

★ ★
youzhizhe(金币+2):谢谢交流。 2010-12-06 19:53:14
对于DOS,因为ISMEAR值只能取-5,所以更改能带计算中的ISMEAR值为默认值(1),做了计算。发现能带形状不会受到ISMEAR取值的影响,但是fermi能级位置会有很大变化。如果考虑到fermi能级的平移,确实能带和DOS就能对应上了。
我现在的问题是fermi能级到底以ISMEAR取何值时的为准呢?ISMEAR= -5(DOS) or 1(能带)?
另:能带计算中ISMEAR取不了-5,所以二者无法统一。
不知我的问题说清楚没有。

[ Last edited by first_yg on 2010-12-6 at 13:10 ]
10楼2010-12-06 13:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 first_yg 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703 化学 求调剂,一志愿山东大学 342 分 +5 Shern—- 2026-03-28 5/250 2026-03-29 23:32 by 飞行日记西
[考博] 26申博自荐 +6 whh869393 2026-03-24 6/300 2026-03-29 21:05 by nxgogo
[考研] 材料与化工328分调剂 +7 。,。,。,。i 2026-03-23 7/350 2026-03-29 20:58 by 唐沐儿
[考研] 化学0703 调剂 306分 一志愿211 +7 26要上岸 2026-03-28 7/350 2026-03-29 20:04 by 无际的草原
[考研] 化学工程085602 305分求调剂 +22 RichLi_ 2026-03-25 22/1100 2026-03-29 19:51 by 无际的草原
[考研] 考研调剂 +7 小蜡新笔 2026-03-29 7/350 2026-03-29 19:00 by 学员8dgXkO
[考研] 一志愿郑州大学,080500学硕,总分317分求调剂 +8 举个栗子oi 2026-03-24 9/450 2026-03-29 13:08 by peike
[考研] 375求调剂 +4 雨夏整夜 2026-03-29 4/200 2026-03-29 12:09 by 无际的草原
[考研] 086000生物与医药调剂 +5 Feisty。 2026-03-28 9/450 2026-03-29 12:02 by longlotian
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 调剂求院校招收 +6 鹤鲸鸽 2026-03-28 6/300 2026-03-29 08:15 by fmesaito
[考研] 调剂考研 +3 王杰一 2026-03-29 3/150 2026-03-29 08:09 by fmesaito
[考研] 329求调剂 +7 星野? 2026-03-26 7/350 2026-03-29 06:43 by 544594351
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 食品工程专硕一志愿中海洋309求调剂 +4 小张zxy张 2026-03-26 8/400 2026-03-28 19:25 by lbsjt
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考博] 26申博 +3 加油冲啊! 2026-03-26 3/150 2026-03-27 15:38 by cls512
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 085601求调剂总分293英一数二 +4 钢铁大炮 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
信息提示
请填处理意见