计算化学公社

标题: 用cp2k计算ZnS能带结果有些问题 [打印本页]

作者
Author: 萧瑟 时间: 2024-8-1 14:02
标题: 用cp2k计算ZnS能带结果有些问题
各位老师，我计算了一个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也都算过了，结果都是带隙为零。

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

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

作者
Author: 萧瑟 时间: 2024-8-2 13:10

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

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

作者
Author: 萧瑟 时间: 2024-8-3 12:45

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

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

作者
Author: sobereva 时间: 2024-8-4 07:29

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

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

作者
Author: 萧瑟 时间: 2024-8-5 09:14

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)