用cp2k计算ZnS能带结果有些问题

萧瑟 · 发表于 Post on 2024-8-1 14:02:08

各位老师，我计算了一个ZnS的结构，在MP上给出的band gap有2.02eV,但是我算出来是0，不知道为什么。昨天在群里问过sob老师了，但感觉群里不太方便说明太多，所以在此发帖问一下各位老师。
优化的是一个ZnS（111）的结构，放到seek-path里给出了k点路径，在inp文件（Multiwfn生成的计算单点的文件）了加上了k点的信息，我看到下面提示让我同时替换计算的晶胞和坐标信息。我用了它给的晶胞坐标和k点路径，用PBE和HLE17都算过，k点7 7 7 和9 9 9也都算过了，结果都是带隙为零。

烦请各位老师帮我看看具体是什么问题
感谢各位！

sobereva · 发表于 Post on 2024-8-2 04:09:20

二维体系在垂直于表面方向考虑k点毫无意义
这么小体系，而且还是二维的，首选RI-HFx方式结合HSE06算能带，显然不应该用PBE。另外CP2K还可以做GW的能带计算

萧瑟 · 发表于 Post on 2024-8-2 13:10:21

sobereva 发表于 2024-8-2 04:09
二维体系在垂直于表面方向考虑k点毫无意义
这么小体系，而且还是二维的，首选RI-HFx方式结合HSE06算能带， ...

感谢sob老师的回复，我这就调整一下算法参数

萧瑟 · 发表于 Post on 2024-8-3 12:45:53

sobereva 发表于 2024-8-2 04:09
二维体系在垂直于表面方向考虑k点毫无意义
这么小体系，而且还是二维的，首选RI-HFx方式结合HSE06算能带， ...

sob老师，我按您说的RI-HFx结合HSE06计算了能带，算出来是0.24（cubic111面）还是比MP上给出的2.0（cubic结构）小了一个数量级，按这个计算精度来说就算再提升精度也不能改变这个差距了。
请问老师，这个计算结果就是我的这个结构的带隙了吗？还是说我计算的有问题，所以这个值就还差这么多？
烦请老师指教。

sobereva · 发表于 Post on 2024-8-4 07:29:24

萧瑟发表于 2024-8-3 12:45
sob老师，我按您说的RI-HFx结合HSE06计算了能带，算出来是0.24（cubic111面）还是比MP上给出的2.0（cubic ...

要和MP对比，就用MP上完全相同的模型，尽可能相同的计算级别（泛函保持一致，并且确保自己的设置没有硬伤）

萧瑟 · 发表于 Post on 2024-8-5 09:14:00

sobereva 发表于 2024-8-4 07:29
要和MP对比，就用MP上完全相同的模型，尽可能相同的计算级别（泛函保持一致，并且确保自己的设置没有硬伤 ...

谢谢sob老师，相同的结构计算了一下结果非常接近了。感谢老师回复

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

[CP2K] 用cp2k计算ZnS能带结果有些问题