castep算窄带隙材料能带时，选hse06的结果却比实验大很多，应该如何调整

zhuyl · 发表于 Post on 2023-4-19 01:55:01

castep计算窄带隙半导体InSb材料能带，没有找到文献参数，但是有很多理论公式和实物实验显示带隙为0.24eV左右。选用HSE06和模守恒，结果出入很大如下图。也计算过GGAPBE的，结果为0.001eV。对参数还是有点一知半解，希望求助该如何调整。

卡开发发 · 发表于 Post on 2023-4-19 06:29:16

模守恒赝势的动能截断大概率400eV是不够的，你可以测试看看。

zhuyl · 发表于 Post on 2023-4-19 08:19:06

卡开发发发表于 2023-4-19 06:29
模守恒赝势的动能截断大概率400eV是不够的，你可以测试看看。

实际上为了避免这个我试过设置1000eV，但是由于HSE06计算时间太长了，用的是GGAPBE。关于这里我再验证一下。那么假如我得到个合理的GGAPBE截断能，在都使用同样的赝势比如模守恒，是否可以直接套用到HSE006呢？

zhuyl · 发表于 Post on 2023-4-19 10:24:07

卡开发发发表于 2023-4-19 06:29
模守恒赝势的动能截断大概率400eV是不够的，你可以测试看看。

我进行了新一轮测试，GGAPBE 模守恒，截断能和k点精度在0.01eV，450eV和8*8*8，得到的能带1.656eV，和理论的0.24eV还是区别蛮大的吧

卡开发发 · 发表于 Post on 2023-4-19 11:57:20

zhuyl 发表于 2023-4-19 10:24
我进行了新一轮测试，GGAPBE 模守恒，截断能和k点精度在0.01eV，450eV和8*8*8，得到的能带1.656eV，和理 ...

我测试了一下，应该是赝势问题，你换OTFG的模守恒。另外，HSE06计算是比较昂贵，除了算能带结果还行，结构优化犯不上用HSE06做。

zhuyl · 发表于 Post on 2023-4-19 13:03:56

卡开发发发表于 2023-4-19 11:57
我测试了一下，应该是赝势问题，你换OTFG的模守恒。另外，HSE06计算是比较昂贵，除了算能带结果还行，结 ...

还是有点不太理想，OTFG的模守恒下算出600eV，9x9x9，然后GGAPBE结构优化的能带看不出带隙。我再用这个结构优化后的文件选HSE06试试能带。

卡开发发 · 发表于 Post on 2023-4-19 13:23:39

zhuyl 发表于 2023-4-19 13:03
还是有点不太理想，OTFG的模守恒下算出600eV，9x9x9，然后GGAPBE结构优化的能带看不出带隙。我再用这个结 ...

GGA-PBE正常结果gap应该是0 eV（Materials Project结果），如果用你用那个模守恒赝势结果PBE算出来的gap是1.2 eV，动能截断和k点我也调整过，包括In换成01的recpot也不行，那么显然可能那个赝势的可转移特性不行并不适合你要研究这个体系。

卡开发发 · 发表于 Post on 2023-4-19 14:06:56

zhuyl 发表于 2023-4-19 13:03
还是有点不太理想，OTFG的模守恒下算出600eV，9x9x9，然后GGAPBE结构优化的能带看不出带隙。我再用这个结 ...

我测试了一下，按你那个路径做出来HSE06结果是0.126 eV

晶胞我没优化，参数也没测试，MS的路径不是好的，但应该包括了该有的k点。

zhuyl · 发表于 Post on 2023-4-20 13:47:18

本帖最后由 zhuyl 于 2023-4-20 13:49 编辑

卡开发发发表于 2023-4-19 14:06
我测试了一下，按你那个路径做出来HSE06结果是0.126 eV

晶胞我没优化，参数也没测试，MS的路径不是 ...

我还是没有跑出来和您这个类似的，至少能说明是窄带隙的结果。我的理解是，从定量上确实还有很多慢慢调整的参数，但是定性上我没得到有打开带隙或者看起来偏差不大的带隙，就像您这个感觉就挺合理的。您这个大概的参数是咋样的您还记得吗，就HSE06 OTFG模守恒，还有什么特别的参数吗

卡开发发 · 发表于 Post on 2023-4-20 14:30:09

zhuyl 发表于 2023-4-20 13:47
我还是没有跑出来和您这个类似的，至少能说明是窄带隙的结果。我的理解是，从定量上确实还有很多慢慢调整 ...

复制代码

zhuyl · 发表于 Post on 2023-4-20 15:33:15

卡开发发发表于 2023-4-20 14:30

非常谢谢您，学习到了，我再调整调整。

zhuyl · 发表于 Post on 2023-4-20 16:10:54

卡开发发发表于 2023-4-20 14:30

确实是赝势有很大问题，我重新分别测试一下了In和Zn的赝势文件，有所改善，只是精度我打算再提高一点，谢谢卡哥

leeru · 发表于 Post on 2023-4-20 18:22:50

HSE06不保证可以获得准确的带隙，试试koopmans这个qe的软件包，可以算出与GW一样精度的带隙，而且对IP和EA都算的比较准

zhuyl · 发表于 Post on 2023-4-20 21:38:51

卡开发发发表于 2023-4-19 13:23
GGA-PBE正常结果gap应该是0 eV（Materials Project结果），如果用你用那个模守恒赝势结果PBE算出来的gap ...

我做了基本的几个赝势测试，在ggapbe下的，像图片里的组合，标黄的组合得到了0.38eV的带隙，但是很纠结总能量很低，包括HSE的两种赝势能量也在-1700eV左右。是不是我应该认为这个黄标的测试是不可靠的呢

zhuyl · 发表于 Post on 2023-4-20 21:49:16

leeru 发表于 2023-4-20 18:22
HSE06不保证可以获得准确的带隙，试试koopmans这个qe的软件包，可以算出与GW一样精度的带隙，而且对IP和EA ...

好的谢谢您，我查查这个方面

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CASTEP/Dmol3/MS] castep算窄带隙材料能带时，选hse06的结果却比实验大很多，应该如何调整