24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3621  |  回复: 30

haowenping

木虫 (著名写手)

引用回帖:
6楼: Originally posted by lonewolf659 at 2011-08-10 16:23:45:
DOS图是有smearing展宽的,峰会变得胖一点。
你减小smearing的值试试。
但是VASP具体怎么设置我就不清楚了,你再等等高手吧,呵呵。

嗯,非常感谢了,只是我取用ISMEAR=-5的时候算dos,跟展宽没关系了。依然十分感激了~
人有多大胆,地有多大产!
11楼2011-08-11 15:11:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
8楼: Originally posted by 水座 at 2011-08-11 09:13:54:
楼主计算的是不是类似graphene的二位结构,试试把非周期方向的k点改为1,或是使用smearing的方法。

是的,我算的是类似于石墨烯的系统,我也试着用K点为11x11x1做过scf,在scf里面用的ISMEAR=0,SIGMA=0.16。在dos计算的时候,我取的K点是也是11x11x1,ISMEAR=-5。但是画出来的dos图仍然是一样的,就是费米能级对应的是价带顶,而不是带隙的正中间。
人有多大胆,地有多大产!
12楼2011-08-11 15:27:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
9楼: Originally posted by jinhua1983 at 2011-08-11 09:33:11:
感觉你这个能带在输出band.dat的时候没有选择新的费米能级作为能量的零点哦,重新输出一下band.dat,选用new fermi level看看。DOS是输出的时候直接用费米能级作为能量的零点的,所以你的DOS估计问题是不大的,关 ...

你说的新的费米能级是什么意思呢?照你这么说,是我的bands计算需要重新计算吗?
人有多大胆,地有多大产!
13楼2011-08-11 15:29:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
10楼: Originally posted by y1ding at 2011-08-11 09:39:07:
也可能是DOS取得点太少。
2D体系有时候,k点取到36~48个才可以。

但是我的体系也有30多个原子,截断能取的有点大,取到了850,下面是我优化截断能的时候的一些数据,你帮我看看,有必要取850吗?
800 -292.301962
750 -292.272347
700 -292.241387
650 -292.174936
600 -292.094963
这是一个结点上算的从800到下面的截断能的取值,另一个结点上的截断能的取值是:
1000 -292.323362
950 -292.320887
900 -292.318806
850 -292.313084
800 -292.301214
750 -292.272892
像我这样取那么大的截断能,是否可以减小,我取850合理吗?如果可以减小的话,那样我应该可以加大K点了,像你说的那样加大到36~48。期待着你的回帖。
人有多大胆,地有多大产!
14楼2011-08-11 15:39:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
haowenping(金币+2): 谢谢了~ 2011-08-12 08:41:03
youzhizhe(金币+2): 谢谢交流。 2011-08-17 15:10:38
引用回帖:
14楼: Originally posted by haowenping at 2011-08-11 15:39:05:
但是我的体系也有30多个原子,截断能取的有点大,取到了850,下面是我优化截断能的时候的一些数据,你帮我看看,有必要取850吗?
800 -292.301962
750 -292.272347
700 -292.241387
650 -292.174936
600 -2 ...

不清楚你的具体体系,但650与600的截断能算出来的总能,只差了0.02%,其实都差不多够了。
VASP计算,除非用于算声子,一般在600都够了。
15楼2011-08-11 17:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinhua1983

金虫 (小有名气)

【答案】应助回帖


haowenping(金币+2): 谢谢 2011-08-12 08:44:40
gzqdyouxia(金币+1): 谢谢提示 2011-08-13 10:39:48
引用回帖:
13楼: Originally posted by haowenping at 2011-08-11 15:29:32:
你说的新的费米能级是什么意思呢?照你这么说,是我的bands计算需要重新计算吗?

不用重新计算,只要在得到band.dat的时候选择新的费米能级作为能量零点就可以。或者你把你的能带统一向上平移,使得你的价带顶位于费米能级的位置。
16楼2011-08-12 08:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
16楼: Originally posted by jinhua1983 at 2011-08-12 08:43:00:
不用重新计算,只要在得到band.dat的时候选择新的费米能级作为能量零点就可以。或者你把你的能带统一向上平移,使得你的价带顶位于费米能级的位置。

我把能带统一上移就能很好的对应上,但是这对吗?我的scf和bands里面的费米能级根本就不是那么大,我在bands和scf里面的费米能级的基础之上,需要再加上0.1547左右的费米能级才能对的上,但是,我觉得这不对吧!
人有多大胆,地有多大产!
17楼2011-08-12 08:46:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)

【答案】应助回帖


uuv2010(金币+1): 多谢提示 2011-08-12 17:31:02
你计算能带的时候换一个ISMEAR试试看
paper我要好多好的的paper
18楼2011-08-12 10:45:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
18楼: Originally posted by liuhuimeiwoo at 2011-08-12 10:45:17:
你计算能带的时候换一个ISMEAR试试看

我换了,ISMEAR=-5,1,0都换过的,只是还是对应不上,用ISMEAR=0的时候,SIGMA=0.01,做的scf,bands,费米能级相差很大了,scf中的费米能级是-2.5,而bands中的费米能级是-2.1左右。纳闷了现在~
人有多大胆,地有多大产!
19楼2011-08-13 08:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinhua1983

金虫 (小有名气)

【答案】应助回帖


gzqdyouxia(金币+1): 谢谢提示 2011-08-13 10:40:03
引用回帖:
17楼: Originally posted by haowenping at 2011-08-12 08:46:21:
我把能带统一上移就能很好的对应上,但是这对吗?我的scf和bands里面的费米能级根本就不是那么大,我在bands和scf里面的费米能级的基础之上,需要再加上0.1547左右的费米能级才能对的上,但是,我觉得这不对吧!

你要加的费米能级能量应该是能带计算时候的那个值,不是构型优化的那个能量哦,我所提到的使用新的费米能级作为能量零点也就是将能带统一向上平移的,道理是一样的,做法上没有问题,可能具体的可能是你那个能量取错了,好好检查一下看看。
20楼2011-08-13 10:10:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 haowenping 的主题更新
信息提示
请填处理意见