24小时热门版块排行榜    

查看: 1725  |  回复: 8

梦里花

新虫 (初入文坛)

[求助] 计算空位石墨烯得到的DOS与文献不符

最近在计算空位石墨烯的态密度,
K点设置:
0 0 0 !G
0 0.5 0!F
0 0.5 0!F
0.333 0.667 0!Q
0.333 0.667 0!G
0 0 0 !G
结果的到的态密度图与文献差别很大,不知道是什么原因造成的,有没有人遇到类似情况,十分需要帮助,在此先行谢过!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

xiaoluoj

金虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
WDD880227: 金币+1, 多谢提示 2012-07-02 10:12:42
这个K点格式是计算band时用的吧,计算dos时好像不是这种格式
2楼2012-07-01 14:31:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
WDD880227: 金币+1, 谢谢提示,鼓励参与 2012-07-02 10:12:59
K点应该用Monkhorst-Pack
但这样的K点dos也不会不会差别“巨大”,请贴INCAR。
3楼2012-07-01 21:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦里花

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by xiaoluoj at 2012-07-01 14:31:01
这个K点格式是计算band时用的吧,计算dos时好像不是这种格式

因为要结合能带一起看,就把两者的K点设成一样的了,不知道也影响大不大
4楼2012-07-02 10:58:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦里花

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by mousekingadv at 2012-07-01 21:30:01
K点应该用Monkhorst-Pack
但这样的K点dos也不会不会差别“巨大”,请贴INCAR。

SYSTEM    = Cu/graphene

     PREC    = Medium
     ENCUT = 400
     ISTART  = 0  
     EDIFF   = 1e-5
     EDIFFG  = -2e-2     
    IALGO   = 48         residue pre-conditioned conjugate gradient
    NPAR   =  1
     LPLANE = .TRUE.
     LCHARG =.TRUE

#     NELMIN  = 3         
     LREAL   = On   optimized real space projection
#    ROPT    = 1E-3
#    ISYM    = 0         symmetry off
#    NWRITE  = 1
    ICHARG=11
RWIGS=0.863 1.312
     DIPOLE CORRECTION
#    IDIPOL = 3
  SPIN POLARIZATION SPECIFICATION
  ISPIN   = 2
# NUPDOWN = 2

  IONIC RELAXATION
     NSW     = 0        number of steps for IOM
     NBLOCK  = 1          we don't need to write a lot during geo opt
     IBRION  = -1         (set to 1-2 for geo opt)
     POTIM   = 0.5      1fs for time step = around 40 a.u.
     ISIF    = 2          calculate force, but not stress
     LWAVE =.F.
  LCHG  = .T.

  ELECTRON OCCUPATION
     ISMEAR = 0  occupation forced by INCAR
#    NBANDS = 60

  DOS Stuff
     SIGMA   = 0.1       check T*S to make sure

这是算DOS的INCAR,加自旋算的,文献上得到的结果自旋向上和向下的态密度是对称的,而我偶的到的不对称,不知道是什么原因造成的。感谢帮助
5楼2012-07-02 11:02:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦里花

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by mousekingadv at 2012-07-01 21:30:01
K点应该用Monkhorst-Pack
但这样的K点dos也不会不会差别“巨大”,请贴INCAR。

额,还有个问题,我同学赝势用的PBE,结果与文献上符合的挺好,而我用的是PW91,不知道赝势对结果影响大不大?谢谢
6楼2012-07-02 11:22:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2012-07-02 18:46:40
引用回帖:
6楼: Originally posted by 梦里花 at 2012-07-02 11:22:59
额,还有个问题,我同学赝势用的PBE,结果与文献上符合的挺好,而我用的是PW91,不知道赝势对结果影响大不大?谢谢...

计算dos还是不能用高对称K点的...实在想dos和band一锅端,请参考http://muchong.com/bbs/viewthread.php?tid=3613684&page=1的K点设置方法。个人无责任无压力认为同样适用GGA。
水平有限,没看出INCAR有什么问题...
我自己的体系上测试过两种赝势没什么区别,但现在大家常用PAW_PBE。
7楼2012-07-02 18:35:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦里花

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by mousekingadv at 2012-07-02 18:35:21
计算dos还是不能用高对称K点的...实在想dos和band一锅端,请参考http://muchong.com/bbs/viewthread.php?tid=3613684&page=1的K点设置方法。个人无责任无压力认为同样适用GGA。
水平有限,没看出INCAR有什么问题 ...

Thank you!
8楼2012-07-03 10:28:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦里花

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by mousekingadv at 2012-07-02 18:35:21
计算dos还是不能用高对称K点的...实在想dos和band一锅端,请参考http://muchong.com/bbs/viewthread.php?tid=3613684&page=1的K点设置方法。个人无责任无压力认为同样适用GGA。
水平有限,没看出INCAR有什么问题 ...

能否留下qq,问题比较多,呵呵
9楼2012-07-03 11:28:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 梦里花 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见