24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1426  |  回复: 13
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 jugengfans 的 3 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jugengfans

金虫 (正式写手)

[交流] 【求助】ISMEAR对费米能级的影响很大!求助

我用同样的条件下改变ISMEAR的值画了三幅TiO2的 DOS图。
图中的能量均没有减去费米能时可以重合
但是!
我发现ISMEAR=0时的费米能为2.6,而其余两个的费米能为1.8
vasp手册上不是说半导体应该用ISMEAR=0么?
然而如果ISMEAR=0的时候,我减去这个费米能,则费米能级跑到价带和导带中间去了!而其余两个跟文献符合的比较好!为什么会有这么大的差异?请教大家!
到底应该用哪个ISMEAR?

费米能级在取不同的ISMEAR时有明显差异,大概差了1eV。这也太离谱了吧

[ Last edited by jugengfans on 2009-8-19 at 01:32 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
Originally posted by fanchen021 at 2009-8-18 14:53:
应该用后面的,DOS计算和 Total Energy 的计算时ISMEAR 的取值不一样

我算总能的时候用不同的ISMEAR也不一样,奇怪了!
3楼2009-08-18 14:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
Originally posted by cationly at 2009-8-18 15:17:
半导体算能带应该用-5吧

谢谢,我的问题是计算总能的时候用ismear=0时费米能级与ismear=-5就会有很大差异,到底以哪个为准呢?请教了!

[ Last edited by jugengfans on 2009-8-18 at 16:19 ]
6楼2009-08-18 15:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
Originally posted by netx_ray at 2009-8-18 16:57:
没怎么用vasp

但是从smearing的作用来说,不同的smearing方法对于fermi能级和总能量是有一定影响的。但是影响不是很大,我觉得应该不会有0.1eV这个量级。

前提条件是 其他参数都收敛很好 例如截断能 K点都足够

我的截断能已经取到700了,K点我各个方向20个,应该已经足够。
问题就是在做静态计算时,用ismear=2,ismear=-5时都能给出比较不错的结果,而在ismear=0(vasp推荐计算半导体用这个值)时竟然会有如此大的差异!

这个费米能级在ismear=0时都跑到带隙最中间去了!这跟文献,跟dmol做出的结果极大的不符啊!

然而ismear=-5 或ismear=0都能给出非常好的结果,让我觉得很郁闷,

[ Last edited by jugengfans on 2009-8-19 at 01:01 ]
10楼2009-08-19 00:55:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
Originally posted by spur at 2009-8-18 19:10:
vasp manual: For the calculations of the DOS and very accurate total energy calculations (no relaxation in metals) use the tetrahedron method (ISMEAR= -5)

谢谢,我是在静态计算的时候就发现总能不一样了,这时还没做DOS呢。
ismear=-5是不是在给半导体做静态计算时也能给出一个准确的值呢?
11楼2009-08-19 00:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

我把ismear=0,取不同sigma时的DOS文件贴出来给大家看看

我取了一个周期的TiO2 ,6个原子,难道是因为周期不够大的原因么?

另外,我静态计算的时候取ismear=0,DOS时取ismear=-5 ,这个时候的费米能级究竟是按照先前静态计算的为准,还是从DOS计算的EIGENVAL文件里读呢?

[ Last edited by jugengfans on 2009-8-19 at 01:25 ]
12楼2009-08-19 01:20:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

看了候柱峰老师的教程。发现,对于半导体来说,本来在带隙之间就没有态,DOS的费米能级变得没有任何意义,因此只要费米能级在在bandgap 之间,计算就算不上有问题。
这是我的理解,不知道对不对。
13楼2009-08-19 01:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jugengfans 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见