24小时热门版块排行榜    

查看: 1879  |  回复: 3

李蕊199387

木虫 (小有名气)

[求助] DOS和能带图为什么和文献上不一样呢 已有1人参与

第一幅和第二幅图是我做的能带图和dos图,第三张图是一篇文献上的能带和dos图,为什么会差那么多呢?
DOS的INCAR:
SYSTEM= Al
ISTART=0;ICHARG=11
ISMEAR=-5;SIGMA=0.2
ENCUT=320
#NSW=300
ISIF=2
#IBRION=2
#POTIM=0.05
LWAVE=.FALSE.
LREAL=.FALSE.
#VOSKOWN=1
LCHARG=.F.
EDIFFG=-0.02
NCORE=4
LORBIT=11
能带的INCAR:
SYSTEM= Al
ISTART=0;ICHARG=11
ISMEAR=0;SIGMA=0.2
ENCUT=320
#NSW=300
ISIF=2
#IBRION=2
#POTIM=0.05
LWAVE=.FALSE.
LREAL=.FALSE.
#VOSKOWN=1
LCHARG=.F.
EDIFFG=-0.02
NCORE=4
POSCAR都是:
Al                                      
   4.05000000000000     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
   Al
     4
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.5000000000000000  0.5000000000000000
  0.5000000000000000  0.0000000000000000  0.5000000000000000
  0.5000000000000000  0.5000000000000000  0.0000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
高对称点是通过MS直接看的:
K-points
20
Line-mode
rec
0.5  0.25  0.75  !W
0.5  0.5   0.5   !L

0.5  0.5   0.5   !L
0.0  0.0   0.0   !G

0.0  0.0   0.0   !G
0.5  0.0   0.5   !X

0.5  0.0   0.5   !X
0.5  0.25  0.75  !W

0.5  0.25  0.75  !W
0.375  0.375  0.75  !K
我在自洽和算dos时都是用的惯用晶胞,而在算能带时,在MS里把晶胞转化为了原胞,然后取得高对称点(为了和文献里的高对称点一样),并且我也算了惯用晶胞的能带,出来的和上面的差不多。是我哪里算的不对吗?怎么能判断能带和dos计算的到底对不对呢?为什么文献上的能带和dos能量那么大,还有他们画能带时又是怎么设置让其横坐标高对称点不用均匀分布,我的每两个点之间距离是一样的。

DOS和能带图为什么和文献上不一样呢
dosme.png


DOS和能带图为什么和文献上不一样呢-1
dosmee.png


DOS和能带图为什么和文献上不一样呢-2
wenxian.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
用原胞计算,原胞静态计算之后算态密度,K点数目增大一倍,ISTART=1。
相互学习,共同进步。
2楼2017-03-06 15:15:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
3楼2017-03-06 15:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李蕊199387

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 对马拉二胡 at 2017-03-06 15:15:06
用原胞计算,原胞静态计算之后算态密度,K点数目增大一倍,ISTART=1。

ISTART=1不是从WAVECAR中读取吗,可是我之前都设置的不输出WAVECAR。在进行DOS和能带时都分别将静态自洽时的哪些文件粘贴过去呢?我是只将静态计算时的CHG*文件复制过去,然后再修改INCAR中的参数,这样可以吗?
4楼2017-03-07 20:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 李蕊199387 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 03:26 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:49 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见