24小时热门版块排行榜    

查看: 2327  |  回复: 18
当前主题已经存档。

whzhang81

铜虫 (初入文坛)

★ ★ ★
jugengfans(金币+1,VIP+0):谢谢你的热心!非常感谢,只剩一个金币了,不足以表示我的感激 8-17 23:49
yjr(金币+2,VIP+0):在给个金币!! 8-18 15:31
你这个,不知道怎么算的,你的INCAR里边首先有个问题,对于Rutile TiO2非金属计算的时候ISMEAR 应该用0,不能用2,不过这个应该不是你能带成这样的原因。

我用DMol3低精度跑了一个能带图,很正常,不知道能不能上传图片。。。。

对于高对称点的坐标和命名,以前有人贴过,对于你这个晶体,打开这个链接,
http://www.cryst.ehu.es/cgi-bin/ ... 136&fig=f4ommmp

你可以对照一下看看。


k-vector label           Wyckoff position           Parameters
CDML         ITA         ITA
GM         0,0,0         1         a         4/mmm         0,0,0
Z         0,0,1/2         1         b         4/mmm         0,0,1/2
M         1/2,1/2,0         1         c         4/mmm         1/2,1/2,0
A         1/2,1/2,1/2         1         d         4/mmm         1/2,1/2,1/2
R         0,1/2,1/2         2         e         mmm.         0,1/2,1/2
X         0,1/2,0         2         f         mmm.         0,1/2,0
11楼2009-08-17 20:00:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlshi

金虫 (正式写手)


yjr(金币+1,VIP+0):谢谢!! 8-18 15:31
增加K 点 和 截断能, 增加SIGMA的值
能使能带图变得光滑
12楼2009-08-17 20:20:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

引用回帖:
Originally posted by whzhang81 at 2009-8-17 20:00:
你这个,不知道怎么算的,你的INCAR里边首先有个问题,对于Rutile TiO2非金属计算的时候ISMEAR 应该用0,不能用2,不过这个应该不是你能带成这样的原因。

我用DMol3低精度跑了一个能带图,很正常,不知道能不能 ...

你好,我先调出这次计算的DOS,红色曲线的是我用GGA - PAW势计算出的结果,在费米面下看来差异不大,费米面上的差异估计是DMOL计算的空带更多导致的,而且带隙也差不多。证明我的计算应该没太大问题,可能能带这方面单独计算的设置有些问题,估计是手动设置K点导致的。让我用你提供的K点再试一试,谢谢·

[ Last edited by jugengfans on 2009-8-18 at 00:41 ]
13楼2009-08-18 00:39:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

用不同的ismear值画出的DOS对比

★ ★
fegg7502(金币+2,VIP+0):鼓励交流2 9-6 18:37
红绿蓝分别对应 ismear=0 ismear=-5 ismear=2 的值
发现ismear只是改变了曲线的平滑度而已。看来只是对本征值的展宽方式有影响,并不影响到本征值。
问题来了:三幅图的设置只有ismear不同,但是ismear=-5 和2时,费米能都在1.86左右,而ismear=0时费米能跑到了2.6!导致费米能级不在价带顶部!按理来说绝缘体或半导体应该用ismear=0的,然而结果差异这么大,到底费米能应该用哪个ismear ?

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

jugengfans

金虫 (正式写手)

顶,问题仍未解决,不知道到底是哪里出了错,能带就是很怪异,
还请大家帮忙看看啊!
15楼2009-09-05 21:09:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

★ ★ ★
wuchenwf(金币+3):谢谢 2010-03-13 21:31
找出错误了,把KPOINTS里的Cartesian改成REC就可以了。原来是写错了空间。
我以为直角坐标是倒空间里的直角坐标系呢,而REC是相对倒格空间的基失的坐标,而其实完全不是这样。Cartesian是真实空间?

[ Last edited by jugengfans on 2009-9-5 at 21:51 ]
16楼2009-09-05 21:49:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaojie7783

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by jugengfans at 2009-9-5 21:49:
找出错误了,把KPOINTS里的Cartesian改成REC就可以了。原来是写错了空间。
我以为直角坐标是倒空间里的直角坐标系呢,而REC是相对倒格空间的基失的坐标,而其实完全不是这样。Cartesian是真实空间?

[ Last  ...

你好 我想请问你用VASP计算能带的时候的KPOINTS文件的写法现在弄对了吗 我也遇到你上面提到的问题 就是K点分割的问题,不知道要按照什么原则分割倒空间,望指教,谢谢。
成功其实很简单,就是当你坚持不住的时候,再坚持一下!
17楼2009-09-08 11:10:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★
aylayl08(金币+5):感谢讨论交流 2010-03-13 08:43
wuchenwf(金币+3):谢谢 2010-03-13 21:31
引用回帖:
Originally posted by xiaojie7783 at 2009-9-8 11:10:
你好 我想请问你用VASP计算能带的时候的KPOINTS文件的写法现在弄对了吗 我也遇到你上面提到的问题 就是K点分割的问题,不知道要按照什么原则分割倒空间,望指教,谢谢。

已经对了。关于高对称点的选取有一些简单的规则,
第一是对称性。不同的晶体结构有不同的对称性,我们应该根据晶体结构选择合适的K点,这个在各种教材中已经给出了。不过手动选择K点只在能带计算的时候是必要的,在做优化和性质计算时我们通常选择自动撒点。
第二是K点分割,K点分割越密,计算的结果就越精确,而耗时也会更多,因此我们要寻找一个比较合适的分割。这个可以在你计算的时候做一个测试,就是从小到大选择不同的分割,然后比较静态计算时的总能变化,自己画一个曲线图,当总能变化范围小于0.2eV时,我们就认为不用再取更多的K点了。

一般来说,K点的选取跟晶胞的边长有关,边长越长,所需的K点就越少。如果是块体的单胞,K点的分割要取密一些,比如我们的TiO2单胞,K点分割12*12*12就非常精确了。而如果我们将这个晶胞延xyz方向各扩展一倍,这时只要6*6*6足够。越大的晶胞,所需的K点数就越小。如果是表面,那么在某个方向有大于10A的真空层,那么这个方向上K点最多只须取2足矣。
18楼2009-09-08 16:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

幽雅百合

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
侯柱峰的教程上说在进行能带结构计算时,ISMEAR和SIGMA用默认值就好,它们的默认值分别是1和0.2啊
19楼2010-03-12 23:20:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jugengfans 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见