24小时热门版块排行榜    

查看: 6430  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见