激子结合能计算

我做过一个二维光电材料的激子结合能计算项目，客户想搞清楚他的单层半导体在室温下能不能稳定存在自由激子，还是说激子结合能太大导致光电响应被束缚态拖慢。这个案例让我把里德伯系列能级和介电屏蔽的关系重新梳理了一遍。

我拿到的体系是一种类过渡金属硫族化合物的单层材料，带隙在1.6 eV附近。先明确一下概念：激子就是激发态电子-空穴对（Wannier激子），电子被空穴通过库仑作用束缚，形成一系列类氢里德伯态。激子结合能 Eb 通常定义为第一激子态（n=1）到导带底（CBM）的能量差。图中左侧给出的就是这套里德伯系列：VBM 在下方，CBM 在上方，n=1、n=2、n=3 三条绿色能级逐级靠近 CBM，n=2→CBM 的束缚能标注为 Eb = 0.42 eV。这个 0.42 eV 在二维材料里算中等偏大，说明激子比较稳定，室温 kT 约 26 meV 远小于它，热离化不显著，光电导可能受激子主导。

计算上，我没有只用简单有效质量近似就交差，而是走了 GW+BSE 路线。先用 VASP 做 G0W0 修正准粒子能隙，这一步把 PBE 下偏小的带隙从约 1.1 eV 拉到 1.8 eV；再用 VASP 的 BSE 计算激发态，得到激子能级。BSE 的核心是把电子-空穴相互作用写成 Bethe-Salpeter 哈密顿量，解出来的本征值就是激子束缚能。实际跑的时候，NBANDS 要足够大（我设到 400），否则高阶里德伯态（n=3 以上）根本出不来；介电函数截断也要兼顾效率，我用的能量截断 300 eV，K 点 12×12×1。

右侧图讲的是介电屏蔽依赖。横轴是介电常数 ε，纵轴是激子结合能 Eb，取了对数坐标。三条曲线分别对应约化质量 μ = 0.10、0.25、0.50 m0，虚线是二维材料由于介电受限带来的 ×4 增强。这个图的物理非常清楚：Eb ∝ μ / ε²。有效质量越大、介电常数越小，激子结合能就越高。二维材料因为外部真空环境屏蔽弱，Eb 会比三维体材料大 3-4 倍。 client’s 材料 μ≈0.25 m0、ε≈5 时，Eb 落在 0.3-0.5 eV 区间，和左侧 BSE 结果一致。

踩坑经验有两个。第一个是 K 点收敛：一开始我用 6×6×1，n=1 和 n=2 能级差只有 0.31 eV，放到 12×12×1 后 Eb 才稳定在 0.42 eV，低频介电函数的虚部尾巴对 Eb 影响很大。第二个是真空层厚度：单层/slab 模型真空层至少 20 Å，我试过 15 Å，镜像 slab 的激子波函数有重叠，导致 n=1 能级被人为压低近 50 meV。第三个是轨道投影：客户后来问激子是 s 型还是 p 型，我回头做 BSE 波函数投影，发现第一激子态主要来自 d 轨道电子和 p 轨道空穴，这在写结论时补了进去。

从我的工程经验看，激子结合能不能只看一个数值，得同时给带隙、里德伯系列和 ε-μ 依赖关系。单层材料如果 Eb > 150 meV，室温激子效应就不可忽略；如果 Eb > 400 meV，发光峰位会显著偏离带隙，光伏器件的开路电压损失也要重新评估。最后我也坦诚地说了局限：GW+BSE 计算量大，对几百原子的体系不太现实；更大数据集还是得靠紧束缚或机器学习的介电模型近似。

关于这类第一性原理光电器件计算，我在 [激子结合能计算](https://www.keyanxueshu.com/category/dft/) 栏目里整理过更多带隙、介电函数和激子性质的实战经验，感兴趣可以先看看。[激子结合能计算](https://www.keyanxueshu.com/)

补充一点：看左侧里德伯系列图时，我常教学生把握两个关键比例。第一是 n=1 到 n=2 的能级间距。如果体系类似氢原子，能级应按 1/n² 缩放，n=1 与 n=2 差应约为 Eb 的 3/4。我们这张图里 n=2 距 CBM 还有一段距离，说明它不是纯氢模型，电子-空穴相互作用有短程修正。第二是右侧对数坐标的斜率。μ=0.25 m0 的曲线在 ε=5 时 Eb≈0.4 eV，到 ε=20 时降到约 0.02 eV，斜率接近 −2，正好对应 Eb∝1/ε²。做报告时把这两个比例摆出来，审稿人会觉得你不仅算了数，还懂背后的标度律。如果材料是二维，别忘记在图中把 2D 增强因子 ×4 标出来，否则和三维体材料对比会吃亏。