24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3806  |  回复: 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小希 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +12 李rien 2026-04-04 12/600 2026-04-10 21:19 by zhouxiaoyu
[考研] 一志愿西北工业大学289 085602 +30 yang婷 2026-04-10 30/1500 2026-04-10 20:02 by jefferylhl
[考研] 362求调剂 +10 我要考大 2026-04-06 14/700 2026-04-10 17:00 by luoyongfeng
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[基金申请] 有爆料,一个青年教师卖房得400万,然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +4 m永o不v言o弃m 2026-04-09 4/200 2026-04-10 11:30 by asy1wn
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 生物学求调剂 一志愿沪9,326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +16 路痴小琪 2026-04-05 16/800 2026-04-10 08:08 by kangsm
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +18 努力奋斗112 2026-04-08 18/900 2026-04-09 23:28 by wolf97
[考研] 269电子信息求调剂,可转专业 +9 独酌wl 2026-04-06 9/450 2026-04-09 20:55 by laoshidan
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 求调剂 +11 wwwwabcde 2026-04-07 11/550 2026-04-07 23:16 by JourneyLucky
[考研] 专硕304找调剂,一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
信息提示
请填处理意见