应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 计算石墨烯能带和态密度有问题
181/2返回列表12下一页
查看: 5477  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

紫海心蝶

木虫 (正式写手)

[求助] 计算石墨烯能带和态密度有问题 已有3人参与

各位大师帮我看看是哪儿有问题，计算石墨烯能带，图看起来不对劲，还有态密度，计算能带时我的INCAR:
SYSTEM=GRAPHENE
ENCUT= 800
ISTART= 1
ICHARG= 11
ISMEAR= 1
SIGMA= 0.2
NBANDS= 5
RWIGS= 0.863
EDIFF= 1e-5
PREC= Accurate
我的POSCAR:
graphene
   1.00252525250000
     2.1304215583187109    1.2299994601826920    0.0000000000000000
     0.0000000000000010   -2.4599989203653809    0.0000000000000000
     0.0000000000000000    0.0000000000000000  -20.0000000000000036
   C
     2
Direct
  0.3333333333333357  0.6666666666666643  0.0000000000000000
  0.6666666666666643  0.3333333333333357  0.0000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
KPOINTS:
k-points
20
Line-mode
Rec
0.000000000000 0.000000000000 0.000000000000 !G
0.000000000000 0.500000000000 0.000000000000 !M

0.000000000000 0.500000000000 0.000000000000 !M
-0.333333333333 0.666666666667 0.000000000000 !K

-0.333333333333 0.666666666667 0.000000000000 !K
0.000000000000 0.000000000000 0.000000000000 !G

请大家帮我看看是参数设置不对还是结构建的有问题

计算石墨烯能带和态密度有问题
能带.jpg


计算石墨烯能带和态密度有问题-1
态密度.jpg
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-21 19:39:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

白小纸bql

新虫 (初入文坛)
引用回帖:
9楼: Originally posted by 白小纸bql at 2014-05-25 01:09:41
楼主，我想请问一下，你的DOS图的比0小一点的地方，是出现了带隙么？而且还不小
能带图和带隙不应该是对应的么？在DOS图中出现一段能量值范围内，DOS=0，但是我看到楼主在BAND图对应的地方并没有这么大的带隙呢？
...

楼主，我发现如果想让DOS图在Dirac点处没有一段区间为0的话，在计算DOS的时候k点需要取密集，比如45*45*1，这样就可以了~因为k点取密集了，单位体积单位能量范围内k点就更多，为0的几率就更小，所以不存在为0的点~我刚学的，请多多指教~
赞一下(1人) 回复此楼
10楼2014-05-25 12:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

mywai520

铁杆木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2014-05-25 15:05:44
你的参数有些小问题：
阶段能太高，500足够了
SIGMA = 0.01这个值必须要小

NBANDS=    RWIGS= 0.863  去掉

祝你好运！！！！
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

紫海心蝶
2楼2014-04-22 00:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

紫海心蝶

木虫 (正式写手)
补充一下计算态密度时，INCAR:
SYSTEM=GRAPHENE
PREC= Accurate
ENCUT= 800
ISTART= 1
ICHARG= 11
RWIGS= 0.863
ISMEAR= -5
EDIFF= 1e-5
IBRION= -1

KPOINTS:
Automatic generation
0
MonhkorstPack
11 11 1
0 0 0
赞一下 回复此楼
3楼2014-04-22 08:25:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

紫海心蝶

木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by mywai520 at 2014-04-22 00:39:18
你的参数有些小问题：
阶段能太高，500足够了
SIGMA = 0.01这个值必须要小

NBANDS=    RWIGS= 0.863  去掉

祝你好运！！！！

谢谢，我用你的建议重新算一下，但是算能带的时候NBANDS这个参数不是必须要设置的吗？我刚开始计算的时候没有加RWIGS这个参数，但是用小程序处理数据，得画图数据时出现一系列unknown，加上这个参数就没有错误了。很感谢你的回复，我试试。
赞一下 回复此楼
4楼2014-04-22 08:30:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
引用回帖:
4楼: Originally posted by 紫海心蝶 at 2014-04-22 08:30:17
谢谢，我用你的建议重新算一下，但是算能带的时候NBANDS这个参数不是必须要设置的吗？我刚开始计算的时候没有加RWIGS这个参数，但是用小程序处理数据，得画图数据时出现一系列unknown，加上这个参数就没有错误了。 ...

那是你小程序的问题，若是你非的用。可以设置LORBIT = 11，我不知道你的RWIGS是否正确。

NBANDS不是一定要，若是空带不够需要设置。看你的设置值也太小了，石墨烯是不要设置的。
赞一下 回复此楼
5楼2014-04-22 08:53:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

紫海心蝶

木虫 (正式写手)
引用回帖:
5楼: Originally posted by mywai520 at 2014-04-22 08:53:16
那是你小程序的问题，若是你非的用。可以设置LORBIT = 11，我不知道你的RWIGS是否正确。

NBANDS不是一定要，若是空带不够需要设置。看你的设置值也太小了，石墨烯是不要设置的。...

我的NBANDS是按照NBANDS=NIONS/2+NELECT/2计算得到的，然后那个RWIGS是POTCAR中以A为单位的那个值，意思是我的小程序有问题吗？
赞一下 回复此楼
6楼2014-04-22 09:42:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
引用回帖:
6楼: Originally posted by 紫海心蝶 at 2014-04-22 09:42:28
我的NBANDS是按照NBANDS=NIONS/2+NELECT/2计算得到的，然后那个RWIGS是POTCAR中以A为单位的那个值，意思是我的小程序有问题吗？...

没有问题，NBAND取的多，只是增加空带。你不要纠正这个问题了，问题是你要算出结果。我不知道你是什么小程序。做能带的多的去了，但是你要理解你的程序是怎么去取能带的。没有程序自己一样可以取，只是程序代替了你人工的操作而已。RWIGS你是在要取也是在算DOS，但是加了LORBIT = 11是不需要取了。这样方便一些，至于理由去看说明。我只是告诉你怎么做。
赞一下 回复此楼
7楼2014-04-22 10:47:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

紫海心蝶

木虫 (正式写手)
引用回帖:
7楼: Originally posted by mywai520 at 2014-04-22 10:47:38
没有问题，NBAND取的多，只是增加空带。你不要纠正这个问题了，问题是你要算出结果。我不知道你是什么小程序。做能带的多的去了，但是你要理解你的程序是怎么去取能带的。没有程序自己一样可以取，只是程序代替了你 ...

我修正了下ISMEAR和SIGMA的值，倒是算对了，再仔细钻研一下这两个参数，很感谢你
赞一下 回复此楼
8楼2014-04-22 14:55:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

白小纸bql

新虫 (初入文坛)
楼主，我想请问一下，你的DOS图的比0小一点的地方，是出现了带隙么？而且还不小
能带图和带隙不应该是对应的么？在DOS图中出现一段能量值范围内，DOS=0，但是我看到楼主在BAND图对应的地方并没有这么大的带隙呢？
求指教
赞一下 回复此楼
9楼2014-05-25 01:09:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 紫海心蝶 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 01:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 跳槽后在研项目怎么办？ +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 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 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭