24小时热门版块排行榜     意得辑论文翻译润色立减200,排版五折

查看: 4750  |  回复: 95

WDD880227

木虫 (正式写手)

[交流] HSE计算能带经验小总结已有42人参与

http://muchong.com/bbs/viewthread.php?tid=4016318首先给出以前总结的HSE计算的链接,这段时间自己用HSE方法计算了能带,72个原子(掺杂体系),42个K点,现在给出自己的一点小经验,希望对大家有帮助。
1.hse band calculation之前不需要进行hse scf calculation,只需要进行pbe scf calculation,并且注意pbe scf calculation时就需要将KPOINTS设置为IBZKPT+inp.kpt,这样hse band calculation时读取的上一步的波函数才是正确的。
2.hse band calculation时,ISTART=1,ICHARG不能取11而是默认值。因为hse时哈密顿量变了,也就是CHARG变了,不能从上一步pbe自洽计算读。
3.能带计算完成后数据处理时,遇到的问题是能级出现交叉,如图中的红线标记,而实际上是不应该交叉的,这里图中交叉是由于显示的问题。本人解决方式是:在origin里每一列Sort Column--Ascending,这样重新排序之后再做图即可。注意:线性产生的点也就是权重不为零的点不能用来做图,只用那些权重为0的高对称点来画能带图。
回复此楼

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

精华网帖收集 VASP and MS 第一性原理 材料计算模拟实用技巧
VASP VASP 第一性原理相关文档 第一性原理计算-tg
第一性原理计算经验 计算-vasp 第一性原理 VASP算例
first principle vasp光学 HSE & 杂化函数 有效质量
杂化泛函方面的研究 vasp 模拟,理论 VASP学习
第一性原理计算 方法、理论、参数设置 凝聚态物理 第一性原理 科研
VASP 计算 VASP报错建模问题集 VASP自学成才 VASP 新计算方法
留学海外

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

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

WDD880227

木虫 (正式写手)

引用回帖:
27楼: Originally posted by repent03 at 2012-03-15 00:22:49:
请问这个inp.k是什么程序?该怎么去取得呢?感谢

inp.kpt是由gk_gw.x产生的,需要输入文件syml

» 本帖已获得的红花(最新10朵)

smiling&fighting...
28楼2012-03-15 08:41:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

引用回帖:
13楼: Originally posted by franch at 2012-03-12 15:20:26:
12×2×9=216个核?

这资源很给力啊。。  算了多久?

O(∩_∩)O~算了80个小时。。。
smiling&fighting...
14楼2012-03-12 15:23:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lightgjx

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主,你好,谢谢你的分享,你的分享帮助了很多人,包括我这个初学者。有个问题想问一下:
1,我还看到一个帖子,http://muchong.com/bbs/viewthread.php?tid=3723568&target=self&page=1   讲的是步骤是先算标准DFT计算,然后HSE计算(前边两步都不改变kpoints),第三步算band。这和楼主的方法不太一样,我不知道为什么?算个DFT然后再算HSE这个有意义吗?
2、楼主的方法是先进行PBE算法(这个时候改变kpoints),如果用HSE的方法的话,kpoints用改变吗?
希望楼主能帮忙告诉我一下,真的不太明白。十分感谢!!!
79楼2015-10-02 19:39:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lxk123572

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
28楼: Originally posted by WDD880227 at 2012-03-15 08:41:56
inp.kpt是由gk_gw.x产生的,需要输入文件syml...

前辈您好!inp.kpt是由gk_gw.x产生的,需要输入文件syml,怎样具体得到inp.kpt的,计算能带的时候除了4个基本文件时,还需要gk_gw.x,和syml吗,请求指点
82楼2016-09-28 11:06:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mika

金虫 (正式写手)

纳米宝贝


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问一下楼主,用HSE算出能带结构之后,用哪一步的E-fermi值来处理能带数据呢?在普通DFT计算是,一般都是取电子自洽的那一步的E-fermi。谢谢。
86楼2017-12-05 20:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

引用回帖:
2楼: Originally posted by uuv2010 at 2012-03-12 10:11:56:
第3条是个挺好的做法
有一个疑问:pbe scf calculation和利用pbe产生的波函数计算出hse的波函数然后计算能带这两种方法给出的结果有差异吗?如果没有差异的话,是什么原因?

没有做利用pbe产生的波函数计算出hse的波函数然后计算能带,接下来准备做测试。区长是怎么认为的?

» 本帖已获得的红花(最新10朵)

smiling&fighting...
3楼2012-03-12 10:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caesarly

金虫 (小有名气)

★ ★
小木虫(金币+0.5): 给个红包,谢谢回帖
franch(金币+1): 鼓励交流。。。。。 2012-03-12 12:41:29
问一下,第一条方法计算结果,是否有实验值可进行参考?有的话,误差有多大?
4楼2012-03-12 11:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)


小木虫(金币+0.5): 给个红包,谢谢回帖
送鲜花一朵
顶楼主,太感谢你的分享,多发此类帖子啊
5楼2012-03-12 11:42:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)


liliangfang(金币+1): 谢谢交流 2012-03-12 14:32:41
引用回帖:
4楼: Originally posted by caesarly at 2012-03-12 11:15:26:
问一下,第一条方法计算结果,是否有实验值可进行参考?有的话,误差有多大?

未掺杂体系的计算结果和实验值几乎吻合,掺杂体系没有实验值~~~
smiling&fighting...
6楼2012-03-12 14:06:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沙风

铁杆木虫 (正式写手)

有为青年


小木虫(金币+0.5): 给个红包,谢谢回帖
引用回帖:
11楼: Originally posted by WDD880227 at 2012-03-12 14:32:48:
12个节点 一个节点2个CPU,一个CPU用了9个核。。。

这计算能力太牛了
15楼2012-03-12 15:25:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

★ ★
franch: 金币+2, 鼓励交流。。。。。呵呵 2012-03-19 10:31:49
引用回帖:
32楼: Originally posted by chuanghua304 at 2012-03-17 15:43:40:
"1.hse band calculation之前不需要进行hse scf calculation,只需要进行pbe scf calculation,并且注意pbe scf calculation时就需要将KPOINTS设置为IBZKPT+inp.kpt,这样hse band calculation时读取的上一步的 ...

IBZKPT在第一步优化计算时就有生成,我在有产生高对称点时,输入的费米能级的位置一直设置为0,只在最后作图的时候把能量零点调到价带顶的位置
smiling&fighting...
35楼2012-03-19 10:15:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

引用回帖:
38楼: Originally posted by 铁昂笑 at 2012-04-08 13:45:44:
我有个疑问:为什么不用HSE-SCF的CHG结果来计算能带?

混合泛函,哈密顿量变了,也就是CHARG变了,导致charge density不准,所以不能用HSE-SCF的CHG结果来计算能带~~
smiling&fighting...
39楼2012-04-09 09:04:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjh9160

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
要把高对称K点的权重设为0吗?
54楼2014-12-19 20:59:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

★ ★ ★
小木虫(金币+0.5): 给个红包,谢谢回帖
franch(金币+2): 鼓励交流。。。。。 2012-03-12 12:41:16
第3条是个挺好的做法
有一个疑问:pbe scf calculation和利用pbe产生的波函数计算出hse的波函数然后计算能带这两种方法给出的结果有差异吗?如果没有差异的话,是什么原因?
2楼2012-03-12 10:11:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主


小木虫(金币+0.5): 给个红包,谢谢回帖
送鲜花一朵
引用回帖:
3楼: Originally posted by WDD880227 at 2012-03-12 10:46:58:
没有做利用pbe产生的波函数计算出hse的波函数然后计算能带,接下来准备做测试。区长是怎么认为的?

没有认真做过,只是测试过几个小例子,所以没搞清楚这些事情
等您昨晚测试,记得来继续分享哦~
7楼2012-03-12 14:08:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

uuv2010: 编辑内容 2012-03-12 14:11
引用回帖:
7楼: Originally posted by uuv2010 at 2012-03-12 14:08:48:
没有认真做过,只是测试过几个小例子,所以没搞清楚这些事情
等您做完测试,记得来继续分享哦~



[ Last edited by uuv2010 on 2012-3-12 at 14:11 ]
8楼2012-03-12 14:10:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caesarly

金虫 (小有名气)


小木虫(金币+0.5): 给个红包,谢谢回帖
由于ISTART=1,第二条中写的ICHARG默认是=0吧?之前参考论坛写的HSE能带计算ICHARG=2.。。。。。问一下这区别大吗?
9楼2012-03-12 14:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沙风

铁杆木虫 (正式写手)

有为青年


小木虫(金币+0.5): 给个红包,谢谢回帖
楼主用HSE计算72个原子体系,你的机器是多少cpu的?
10楼2012-03-12 14:21:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 WDD880227 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 大学老师招聘必须是应届博士吗 +3 魅力女王 2021-04-20 4/200 2021-04-20 23:44 by 石油抄人
[硕博家园] 快毕业了,投稿要不要加导师名字 +14 调剂生2018 2021-04-20 21/1050 2021-04-20 23:39 by 123春
[论文投稿] 论文一审进入Decision了 +6 Felix000 2021-04-19 6/300 2021-04-20 22:51 by quan3145
[考博] 博士调剂 +7 思-勤 2021-04-20 13/650 2021-04-20 22:09 by dxj900720
[考博] 读博的话就要异地了 +13 云朵ing 2021-04-20 19/950 2021-04-20 21:43 by jqy2019
[考博] 一篇ccf c期刊 +4 charon_s 2021-04-20 5/250 2021-04-20 14:24 by charon_s
[考博] 深圳的研究生院哪个好一些? +5 小白zhen菜 2021-04-19 5/250 2021-04-19 21:27 by yze123
[考博] 博导回复欢迎报考 还不确定有没有名额 是婉拒吗 +8 浊上归去 2021-04-19 9/450 2021-04-19 19:37 by 偏振片
[论文投稿] 询问期刊 10+3 雨启青晨 2021-04-17 6/300 2021-04-19 13:41 by 咿呀学语
[考研] 07调剂最难能有多难 +23 gmst8853 2021-04-16 46/2300 2021-04-19 09:17 by 上岸188
[有机交流] 有机化学就业怎样 +10 角动量神秘 2021-04-18 12/600 2021-04-18 20:16 by upc加油中国
[考研] 391求调剂 +3 单佳缘 2021-04-16 4/200 2021-04-18 18:18 by 单佳缘
[硕博家园] ..... +10 xielan40 2021-04-17 11/550 2021-04-17 22:44 by 河大-芒果
[硕博家园] 有些迷茫 +14 努力变优秀吧 2021-04-14 18/900 2021-04-17 14:18 by 行者武松
[考研] 初试357分没学上了,真难!谁能给孩子一个机会? +3 阿雪1393 2021-04-16 4/200 2021-04-16 20:53 by 背多分0423
[硕博家园] DFT有偿计算 +3 001陌上花开 2021-04-15 3/150 2021-04-16 10:08 by fanqinghai
[考研] 很迷茫 +5 嗯嗯嗯啊! 2021-04-15 6/300 2021-04-16 02:25 by fnnfcw
[考研] 370求调剂 +7 wsy41 2021-04-14 10/500 2021-04-15 13:12 by ccf水叮当
[考研] 299求调剂 +3 5791 2021-04-14 3/150 2021-04-14 17:22 by zhujunfly
[考研] 308求调剂 +3 mkqj5336 2021-04-14 3/150 2021-04-14 10:28 by 国国1112
信息提示
请填处理意见