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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 计算石墨烯能带和态密度有问题
181/2返回列表12下一页
查看: 5055  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 0703化学调剂 ，六级已过，有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
关闭
请填处理意见
关闭