24小时热门版块排行榜    

查看: 3767  |  回复: 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小希 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 本人考085602 化学工程 专硕 +19 不知道叫什么! 2026-03-15 21/1050 2026-03-20 20:48 by zhukairuo
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
信息提示
请填处理意见