24小时热门版块排行榜    

查看: 4906  |  回复: 17

紫海心蝶

木虫 (正式写手)

[求助] 计算石墨烯能带和态密度有问题 已有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
回复此楼

» 猜你喜欢

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

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

已阅   回复此楼   关注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的点~我刚学的,请多多指教~
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的回帖
相关版块跳转 我要订阅楼主 紫海心蝶 的主题更新
信息提示
请填处理意见