结合案例图分析,第一张案例图横轴是高对称点Γ→M→K→Γ,纵轴是能量。红色曲线是导带底(CBM),蓝色曲线是价带顶(VBM)。可以看到VBM和CBM都位于K点,带隙E_g≈1.45 eV(图中绿色双向箭头标注),这是一个典型的直接带隙半导体特征。直接带隙对光电器件很重要,因为光跃迁不需要声子参与,辐射复合效率高。Γ→M方向的能带色散比较平缓,说明有效质量较大;而K点附近色散较陡,意味着载流子迁移率可能较高。

第二张案例图横轴是DOS(states·eV^{-1}),纵轴是能量。蓝色曲线代表价带,红色曲线代表导带。价带顶在0 eV附近有一个强峰,主要来自p_z轨道;导带底在约1.45 eV处也有一个峰,主要由d轨道贡献。费米能级E_F位于0 eV。DOS在带隙区间几乎为零,说明材料本征情况下绝缘性较好,适合做场效应晶体管沟道。
从能带结构和态密度的组合来看,这是一个二维材料的DFT计算项目——客户想确认他合成的单层半导体是直接带隙还是间接带隙,以及态密度分布是否适合做光电器件。能带图回答了带隙类型和载流子有效质量的问题,DOS图补充了轨道成分和态密度峰位的信息。
项目计算模型我用的是单层2H相,超胞约20×20×1 Å,真空层18 Å,避免层间相互作用。泛函用PBE,因为客户主要关心带隙类型而非精确数值;如果需要更准的带隙,我会加GW或HSE06修正。K点路径取Γ-M-K-Γ,每个方向采20个点。平面波截断能500 eV,总能收敛到10^{-6} eV。
数据处理时我重点关注了三个量。第一个是带隙类型:通过比较CBM/VBM是否在同一k点判断,这里clearly直接带隙。第二个是有效质量:我在K点附近做了抛物线拟合,电子有效质量约0.35 m_0,空穴约0.42 m_0,都属于中等。第三个是激子结合能:虽然这张图没给出,但我用Wannier-Mott模型估算,由于二维介电屏蔽弱,Eb可能在0.3-0.5 eV量级,后续如果做发光器件需要单独算BSE。
实际计算中有两个坑。一个是真空层厚度:我用12 Å时带隙被镜像层影响偏小约0.1 eV,换到18 Å才稳定。另一个是范德华修正:层内是共价键,层间没有相互作用,所以不需要DFT-D;但如果算的是多层堆叠,必须加DFT-D3。另外单层模型如果有电荷转移,真空层内电势会有斜率,建议加LDIPOL和IDIPOL=3做偶极修正。
这张图对客户的核心价值是:材料是直接带隙、带隙约1.45 eV,落在可见光到近红外区间,适合做光探测器或LED。关于二维材料DFT计算的更多经验,我在二维材料dft计算栏目里做了汇总
cp2k第一性原理计算
CP2K模拟计算
CP2K模拟计算
castep计算吸附能:Materials Studio里的表面建模与能量收敛复盘
CP2K计算能带:大体系电子结构模拟的混合基组方案
能带理论计算:固体能带结构的DFT模拟方法与工程应用
CP2K计算能带:混合基组DFT方法在周期性体系中的实战应用
CP2K分子动力学模拟详解:大体系加速策略与GPW方法实战