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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 计算石墨烯能带和态密度有问题
181/2返回列表12下一页
查看: 5050  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +8 不知道叫什么！ 2026-03-15 9/450 2026-03-16 07:45 by L135790
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 求调剂！ +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 328，0703考生求调剂，一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
信息提示
关闭
请填处理意见
关闭