半导体技术是现代信息社会的基石,从微处理器、存储器到光伏电池和发光二极管,其核心物理性质都源于其独特的电子能带结构。能带结构清晰地描绘了电子在晶体动量空间中的能量分布,决定了材料是金属、半导体还是绝缘体,并直接关联其电学、光学、热学和磁学性质。因此,准确获取并深入理解半导体的能带结构,对于新材料设计和器件应用至关重要。
第一性原理计算,特别是基于密度泛函理论(Density Functional Theory, DFT)的方法,已经成为计算材料科学中研究电子结构最强大和最常用的工具。它仅基于原子种类和初始结构,通过求解量子力学基本方程来预测材料的各种性质。然而,DFT计算输出的能带结构数据包含了丰富的信息,如何正确地分析这些数据,并理解其与宏观物理性质的联系,是计算工作者必须具备的核心技能。
华算科技朱老师将系统性地阐述如何对半导体材料的DFT能带计算结果进行全面的分析。内容将涵盖从基本参数的读取、关键特性的判断,到深入物性的关联与预测,并讨论DFT方法的固有局限性及其解决方案。
在开始分析之前,我们首先需要理解能带结构是什么,以及DFT是如何计算它的。
在孤立原子中,电子处于分立的能级上。当大量原子周期性排列形成晶体时,由于原子轨道的重叠和电子波的干涉,原本分立的能级会展宽成一系列连续的能级,称为能带。允带(导带和价带)之间是禁止电子存在的能量区域,称为禁带。
能带结构图 通常以能量E为纵轴,以晶体动量 k(在倒易空间中)为横轴进行绘制。横轴上的特殊点(如Γ, X, L, K, W等)是高对称性点,它们之间的路径反映了不同晶体方向上的能量-动量关系。
对于半导体而言,其能带结构的核心特征是存在一个带隙,即价带顶(Valence Band Maximum, VBM)和导带底(Conduction Band Minimum, CBM)之间的能量差。这个带隙的存在是半导体一切特性(如可控导电性、光吸收阈值等)的物理根源。
结构优化:首先,使用DFT对晶胞的原子位置和晶格常数进行弛豫,找到能量最低的稳定结构。这是后续所有性质计算的基础。
自洽计算:在优化后的结构上,进行自洽场(SCF)计算,以确定体系的基态电荷密度和总能。
非自洽计算:基于SCF得到的电荷密度,沿着布里渊区预设的高对称性路径,计算不同k点的本征值(即能级),从而绘制出能带结构图。
在这个过程中,交换关联泛函的选择至关重要,它直接决定了计算结果的准确性。
拿到一张DFT计算出的能带图后,我们应系统地审视以下几个核心要素。
直接带隙:价带顶和导带底位于布里渊区同一个k点。例如,GaAs的VBM和CBM都位于Γ点。
物理意义:电子从价带顶跃迁到导带底时,动量变化为零。这种跃迁只需要吸收或发射一个光子,概率很大,因此直接带隙半导体通常具有很高的光吸收系数和发光效率,是光电子器件(如LED、激光器)的理想材料。
间接带隙:价带顶和导带底位于布里渊区不同的k点。例如,Si的VBM在Γ点,而CBM在Γ-X方向上的某个点。
物理意义:电子跃迁过程涉及巨大的动量变化。为了满足动量守恒,这种跃迁必须同时吸收或发射一个声子(晶格振动量子)。这是一个三体过程,发生的概率远低于直接跃迁。因此,间接带隙半导体通常发光效率很低,但广泛应用于微电子领域。
分析方法:在能带图上,直接观察VBM和CBM所处的k点位置。如果两者对齐,则为直接带隙;否则为间接带隙。
带隙值(Eg)是半导体最重要的参数之一。从能带图上,Eg = CBM的能量 – VBM的能量(通常将VBM设为能量零点)。
DFT的“带隙问题”:这是DFT分析中必须面对的核心挑战。使用传统的局域密度近似(LDA)或广义梯度近似(GGA)泛函时,DFT会系统性低估半导体的带隙,通常只有实验值的50%-80%。例如,计算Si的带隙可能只有0.6 eV左右,而实验值为1.12 eV。
原因:DFT理论本身是针对基态总能量的,而带隙是一个激发态性质。Kohn-Sham方程中的本征值本不应被严格解释为准粒子能级,但实践中我们仍将其近似作为能带结构。LDA/GGA泛函无法准确描述电子-空穴之间的相互作用,导致了带隙低估。
杂化泛函:如HSE06,混合了部分精确的Hartree-Fock交换能,能显著改善带隙计算精度,是目前平衡精度与计算成本的最佳选择之一。
GW近似:这是一类更高级的多体微扰理论方法,可以非常精确地计算准粒子能级和带隙,但计算成本极其高昂。
DFT+U:对于强关联体系(如过渡金属氧化物),通过在局域d或f轨道上加上一个有效库仑排斥势U来修正,但对普通半导体作用有限。
在分析报告中,必须明确指出所使用的泛函,并对计算带隙与实验值的差异进行讨论。
能带在k空间的弯曲程度,即能带色散,直接决定了载流子(电子和空穴)的有效质量。
物理意义:在晶体中,电子不像自由电子那样运动,它感受到周期势场的作用。其对外场的响应(加速度)可以用一个“有效质量”m*来描述。
“平”的能带(色散小,曲率小)意味着大的有效质量。载流子迁移困难,电导率低。
“陡”的能带(色散大,曲率大)意味着小的有效质量。载流子迁移容易,电导率高。
识别主导能带:在导带底附近,找出能量最低的那条能带(CBM所在能带)。在价带顶附近,找出能量最高的那条或几条能带(VBM所在能带)。
评估曲率:目测或通过计算拟合CBM和VBM附近的能带曲率。
对于电子(n型),我们关心CBM的曲率。例如,GaAs的导带底非常陡峭,预示着其电子迁移率很高。
对于空穴(p型),我们关心VBM的曲率。价带顶通常比较复杂,可能由多条能带简并或接近。需要分析重空穴带和轻空穴带(在Γ点通常简并),以及自旋轨道耦合分裂带。轻空穴带曲率大,有效质量小,对导电贡献大;重空穴带则相反。
仅凭能带图无法知道特定能带是由哪些原子的哪些轨道贡献的。这时需要结合投影态密度 或轨道投影能带图 进行分析。
通过计算,将总的态密度分解为不同原子种类、不同角动量轨道(s, p, d, f)的贡献,并用不同颜色或线条表示在能带图上。
判断能带性质:可以清晰地看到VBM主要由哪些轨道构成(如O-2p, Se-4p),CBM主要由哪些轨道构成(如Ti-3d, Si-3p)。这有助于理解材料的化学键合和光学跃迁选择定则。
设计能带工程:例如,如果发现CBM主要由阳离子的d轨道贡献,那么通过替换不同的阳离子,可以系统地调控导带的位置和色散。
理解掺杂行为:如果引入一个掺杂原子,其能级是位于价带、导带还是带隙中(形成缺陷能级),取决于其轨道与宿主材料能带的相互作用。态密度投影为此提供了微观图像。
🏅 华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。
🎯500+博士团队护航,累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果,计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!