24小时热门版块排行榜    

查看: 3584  |  回复: 16

韩D小希

铁虫 (正式写手)

[求助] 杂化泛函求大家帮助

刚开始接触杂化泛函计算带隙,自己也查了一些东西,有点疑问请教大家,先感谢。
我先列一下计算硅能带的步骤,希望大家给看一下:
1先正常结构优化,GGA-PBE;
2.将CONTCAR 转为POSCAR,INCAR中去掉离子步做静态计算,得CHGCAR
3.将CHGCAR 文件拷过来,设置ICHGCAR=11,LWAVE=.TURE.,IBRION=-1,设置杂化泛函参数:
LHFCALC = .TRUE.
HFSCREEN = 0.2  (To conform with the HSE06 functional you need to select (HFSCREEN=0.2))
ALGO = Damped
TIME = 0.4
ENCUTFOCK = 0
AEXX=?
4.BAND计算;
拷贝上一步的WAVECAR,然后设置新的KPOINTS,计算即可。

关于以上步骤,大侠们帮忙看下有没有不妥的地方,另外有个小白问题,就是HSE00.HSE03,HSE06有什么区别么,是不是指AEXX的比列不同啊?
希望大家踊跃讨论下杂化泛函,你理解的杂化泛函,谢谢各位。
回复此楼

» 收录本帖的淘帖专辑推荐

VASP相关 杂化泛函方面的研究

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liliangfang

荣誉版主 (著名写手)

引用回帖:
9楼: Originally posted by wkfan8988 at 2013-04-28 12:12:48
将高对称点加在IBZKPT后面,但是高对称点很有限,画出来的能带图很粗略,

http://cms.mpi.univie.ac.at/wiki/index.php/Si_bandstructure

请问这种情况怎么解决?
另外,你的方法算得带隙准吗?http://emuc ...

在你取出的高对称点之间可以在插入一些点,这样就好了,小木虫以前有人发过可以取高对称点的小程序,可以取多一点,自己写也是一样的
10楼2013-04-28 14:40:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
韩D小希: 金币+10, ★★★★★最佳答案 2013-03-23 10:22:07
http://muchong.com/bbs/viewthread.php?tid=3613684
仅供参考,至于参数可参照vasp说明书,03 06 有一点不同,说明书上面有介绍
2楼2013-03-22 15:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩D小希

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by liliangfang at 2013-03-22 15:42:16
http://muchong.com/bbs/viewthread.php?tid=3613684
仅供参考,至于参数可参照vasp说明书,03 06 有一点不同,说明书上面有介绍

谢谢您的回答,我想问下为什么我在做HSE这一步的时候,输出文件没有IIBZKPT这一文件呢,您是否有遇到这样的情况,期待您的回答
3楼2013-03-22 20:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 韩D小希 at 2013-03-22 20:32:02
谢谢您的回答,我想问下为什么我在做HSE这一步的时候,输出文件没有IIBZKPT这一文件呢,您是否有遇到这样的情况,期待您的回答...

不能啊,有的,你是不是没计算完啊
4楼2013-03-23 09:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩D小希

铁虫 (正式写手)

引用回帖:
4楼: Originally posted by liliangfang at 2013-03-23 09:26:40
不能啊,有的,你是不是没计算完啊...

恩  是我参数的失误,把ENCUTFOCK=o  换成PRECFOCK=F 就好了  VASP版本的问题
另外在麻烦您帮看下我按照您的步骤得到的硅的能带图,
您可以看到,十分的尖锐,这是否是由于ISMEAR=0导致没有展宽造成的呢。这样子图片未免有点丑,请问您是怎么处理的呢。
另外我们都知道能带横坐标是K点,在用DFT计算的时候很容易标定,但是现在HSE后,我的IBZKPT内有36个kpoints,高对称点为81个(就是之前不用杂化泛函的时候选择的k点路径)加起来一共117个,但是现在从EIGENVAL中的数据看,每一条带分成了235个K点,那此时怎么对横坐标进行标定呢,谢谢您

未命名.jpg

5楼2013-03-23 10:21:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

引用回帖:
5楼: Originally posted by 韩D小希 at 2013-03-23 10:21:53
恩  是我参数的失误,把ENCUTFOCK=o  换成PRECFOCK=F 就好了  VASP版本的问题
另外在麻烦您帮看下我按照您的步骤得到的硅的能带图,
您可以看到,十分的尖锐,这是否是由于ISMEAR=0导致没有展宽造成的呢。这样子 ...

ISMEAR=0 默认的SIGMA是0.1,你可以设置为0.05或者你可以测试一下选小一点。另外高对称点81个不算少了,算完能带只取高对称点的部分画图。另外请检查一下没有办法的参数设置
6楼2013-03-23 10:59:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩D小希

铁虫 (正式写手)

引用回帖:
6楼: Originally posted by liliangfang at 2013-03-23 10:59:42
ISMEAR=0 默认的SIGMA是0.1,你可以设置为0.05或者你可以测试一下选小一点。另外高对称点81个不算少了,算完能带只取高对称点的部分画图。另外请检查一下没有办法的参数设置...

好的,我试一下,谢谢您,另外还有个小问题,就是K点后面的权重是什么意思呢,为什么HSE的时候一定要用到IBZKPT里的k点呢,只是用高对称点会有什么区别么?再次打扰,期待指教。
7楼2013-03-23 12:32:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

引用回帖:
7楼: Originally posted by 韩D小希 at 2013-03-23 12:32:59
好的,我试一下,谢谢您,另外还有个小问题,就是K点后面的权重是什么意思呢,为什么HSE的时候一定要用到IBZKPT里的k点呢,只是用高对称点会有什么区别么?再次打扰,期待指教。...

为了能计算所以用IBZKPT里面的,后面加高对称点,后面的权重原来时1,改为0
8楼2013-03-23 15:25:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wkfan8988

禁虫 (小有名气)

本帖内容被屏蔽

9楼2013-04-28 12:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 韩D小希 的主题更新
信息提示
请填处理意见