24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3856  |  回复: 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小希 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 6/300 2026-04-19 13:00 by Aaron_zyn
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考博] 申博/考博 +3 啃面包的小书虫 2026-04-17 4/200 2026-04-17 23:54 by 阳阳阳^_^
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
信息提示
请填处理意见