计算化学公社

标题: 求助：k点采样到晶胞真空区得到了band gap是否正常 [打印本页]

作者
Author: 柠夏七开 时间: 2026-1-27 11:44
标题: 求助：k点采样到晶胞真空区得到了band gap是否正常
如题，我使用cp2k计算如下体系的band时，利用seekpath得到了建议的k点采样路径，其中明显band gap主要来自于A-E与C2-Y2，但是我对应到晶胞发现这是两个真空区。因为得到的gap太小，虽然没做TDDFT，但是实验测得的导电率（实验中相当低）、光学gap（比计算结果大1 eV以上）都与之相差甚远，我认为应该是可以说明一些问题的。
gap对应的HOCO和LUCO主要由真空区贡献，这真的对吗？
泛函选用PBE/HLE17结果相近，基组用的DZVP-MOLOPT-SR-GTH（用TZVP级别也是相近的），k点为1x3x3

作者
Author: 柠夏七开 时间: 2026-1-27 11:45
补充一下各个采样点的坐标： &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT GAMMA    0.0000000000    0.0000000000    0.0000000000
            SPECIAL_POINT    Z    0.0000000000    0.5000000000    0.0000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT    Z    0.0000000000    0.5000000000    0.0000000000
            SPECIAL_POINT    D    0.0000000000    0.5000000000    0.5000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT    D    0.0000000000    0.5000000000    0.5000000000
            SPECIAL_POINT    B    0.0000000000    0.0000000000    0.5000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT    B    0.0000000000    0.0000000000    0.5000000000
            SPECIAL_POINT GAMMA    0.0000000000    0.0000000000    0.0000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT GAMMA    0.0000000000    0.0000000000    0.0000000000
            SPECIAL_POINT    A -0.5000000000    0.0000000000    0.5000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT    A -0.5000000000    0.0000000000    0.5000000000
            SPECIAL_POINT    E -0.5000000000    0.5000000000    0.5000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT    E -0.5000000000    0.5000000000    0.5000000000
            SPECIAL_POINT    Z    0.0000000000    0.5000000000    0.0000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT    Z    0.0000000000    0.5000000000    0.0000000000
            SPECIAL_POINT C_2 -0.5000000000    0.5000000000    0.0000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT C_2 -0.5000000000    0.5000000000    0.0000000000
            SPECIAL_POINT Y_2 -0.5000000000    0.0000000000    0.0000000000
            NPOINTS 10
         &END KPOINT_SET
         &KPOINT_SET
            UNITS B_VECTOR
            SPECIAL_POINT Y_2 -0.5000000000    0.0000000000    0.0000000000
            SPECIAL_POINT GAMMA    0.0000000000    0.0000000000    0.0000000000
            NPOINTS 10

作者
Author: UW_0728. 时间: 2026-1-27 14:46
1,3,3的k点采样真的合理吗，从晶胞结构看a方向边长明显短于bc方向啊，你这应该也不是二维周期性计算

作者
Author: 柠夏七开 时间: 2026-1-27 15:31

UW_0728. 发表于 2026-1-27 14:46
1,3,3的k点采样真的合理吗，从晶胞结构看a方向边长明显短于bc方向啊，你这应该也不是二维周期性计算

因为实际上是二维层状结构，我之前也试过2x3x3，和3x2x2，结果没有显著的变化

作者
Author: sobereva 时间: 2026-1-28 16:03
你对k点采样的理解有问题。不需要顾虑你说的这个问题

作者
Author: 柠夏七开 时间: 2026-1-29 14:53

sobereva 发表于 2026-1-28 16:03
你对k点采样的理解有问题。不需要顾虑你说的这个问题

好的，谢谢老师，我后来发现了是自己的理解问题，把k点采样误认为为在实空间采样了，应该是在布里渊区中

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