态密度（DOS）峰位、高度、宽度代表什么？从计算来源、相互约束到物理成因与判读条件解析

说明：本文华算科技主要介绍态密度峰位、高度和宽度的计算来源、相互约束、物理成因与判读条件。

DOS峰怎样从本征值分布形成？

周期性固体计算会在每个 k 点得到一组能带本征值 ε_nk。总态密度可写成 D(E)=Σ_n,kw_kδ(E−ε_nk)，其中 w_k 是布里渊区采样权重。DOS 统计单位能量区间内有多少电子态，它把带结构中的动量信息积分掉，只保留能量分布。

图1. 抛物线能带、Kane 能带及多条能带的色散关系与对应态密度。DOI：10.1016/j.mtelec.2022.100002

图1清楚展示了这种映射。抛物线带边在三维固体中产生近似 D(E)∝√|E−E_b| 的起伏；Kane 型非抛物线能带离开带边后呈现另一种增长率；第二条能带进入同一能量区间时，总 DOS 是各条能带贡献之和。峰并不等同于一条孤立能级，它常由许多 k 点、多个能带或若干轨道投影叠加而成。

理想 δ 函数无法直接画成连续曲线，程序会采用四面体积分，或用高斯、Methfessel–Paxton 等核函数替代 δ 函数。核函数带有有限能宽，离散本征值便成为可见曲线。普通 Kohn–Sham DOS 的有限峰宽含有作图与积分参数，不宜直接当成电子寿命。

带结构与 DOS 的对应还能解释尖峰来源。能带在某个能量附近较平坦时，大量 k 点落入很窄的能量区间，DOS 会快速升高；能带色散较强时，同一批态铺在更宽的能量范围内，曲线趋于低而宽。鞍点、极值点与近乎平坦的色散可产生 Van Hove 奇点，尖峰记录的是本征值在能量轴上的聚集。

图2. SrTiO₃、PbS 和 CoSb₃ 的能带色散与态密度特征。DOI：10.1016/j.mtelec.2022.100002

图2中 PbS 价带内的尖锐结构与能带斜率接近零的位置相对应，CoSb₃ 则显示多条能带叠加后 DOS 增长率发生改变。看峰之前应核对能量轴方向、DOS 轴方向与零点定义；坐标方向改变外观，却不改变单位能量内的态数含义，有的图把能量画在横轴，有的画在纵轴。

峰位、高度和宽度受同一谱权重约束

峰位先要写明采用哪一种定义。孤立且近似对称的峰可以用最大值所在能量 E_max；不对称峰、肩峰或重叠峰更适合报告能量质心 E_c=∫ED(E)dE/∫D(E)dE。两种数值可能相差明显。直接说“峰向低能移动”，却不写能量参考与取值方法，无法支持轨道能级下降之类的判断。

跨结构比较时，把每个体系各自的 E_F 设为 0 eV 会消去绝对能量差。若关注掺杂、界面或表面引起的能级移动，应采用一致参考，例如真空能级、深层芯能级或体相平均静电势，并说明对齐步骤。峰位差只有在共同能量标尺下才有可比含义。

峰高是 D(E) 在某个能量点的函数值，单位常见为 states·eV⁻¹·cell⁻¹、states·eV⁻¹·atom⁻¹ 或归一化后的任意单位。它受晶胞原子数、自旋求和、投影球半径、能量步长和展宽参数影响。纵轴归一化不同，峰高的数值比较随即失效；同图、同设置下的相对变化才适合讨论。

峰宽常用半高全宽 FWHM 表示，仅适合基线清楚、近似单峰的情形。Li₂GeO₃ 的原子与轨道投影 DOS 含有不对称峰、肩部以及多个轨道的重合，硬套一个 FWHM 会把多种成分压成一个数。肩峰应先辨认成分，再决定是否拟合；拟合函数与能量窗口会改变分峰结果。

图3. Li₂GeO₃ 的总态密度、原子投影与轨道投影态密度。DOI：10.1038/s41598-021-84506-0

选定能量窗口后，曲线面积 A=∫D(E)dE 给出该窗口内的态数；对 PDOS 而言，它是指定原子或轨道的投影权重，数值还受投影方案限制。对于面积近似不变且线形相似的一组峰，宽度增大通常伴随高度降低，可粗略理解为 H∼A/W。面积约束把峰高与峰宽连在一起，高峰不自动对应更多电子，也不自动对应更强成键。

位移/分裂/展宽与权重转移怎样识别？

若一组对应峰在共同能量标尺下整体平移，线形与积分面积变化很小，局域静电势、晶体场平均能量或化学势改变是可检验的来源。若只有某个轨道分量移动，总 DOS 变化较弱，应回到 PDOS 与能带查看该轨道。整体平移与局部权重转移是两类现象，两者对电荷转移的证据强度不同。

一个峰分成两个或多个峰时，常见来源包括配位环境改变引起的晶体场分裂、磁性体系中的交换分裂、结构畸变解除简并，以及两个轨道杂化后形成成键态与反键态。此时峰顶会下降、表观宽度会增加，原峰的权重分散到多个能量区间。分裂要在轨道、原子或自旋分量中找到各自归属，仅凭总 DOS 的双峰形状无法指定机理。

图4. 有序、随机及两类隐序模型的电子态密度、声子态密度与总能。DOI：10.1038/s41467-023-40063-w

图4中组分相近而局域畸变关联方式不同的四种模型，在电子 DOS 中分别保留、填充或打开不同能隙；部分能量区的谱权重降低，别处则出现尖峰或稀疏态。结构关联改变的是整段电子态分布，用单个峰的高低概括这种重排会遗漏能隙与邻近峰的联动。

轨道杂化增强常使原先较集中的态分散到更宽能区，也可能产生新的成键—反键分裂；无序会让不同局域环境贡献稍有偏移的能级，构型平均后表现为展宽。宽峰可能包含强色散、杂化、简并解除或环境分布。判断哪一种占主导，要查看能带、分轨道 DOS、实空间电荷以及不同构型的统计结果。

缺陷共振态能把三种几何变化集中在一组曲线上。Bi₂Se_3−xTe_x 表面模型中，点缺陷势 U 改变时，DOS 与缺陷位置 LDOS 的峰位、峰高和线形共同变化；动量分辨谱函数还显示共振态如何扰动 Dirac 锥。局域峰应与总 DOS 和动量分辨谱互相核对，否则局域强峰可能被误写成整个体系新增一条平坦能带。

图5. 不同缺陷势下拓扑绝缘体表面的 DOS、LDOS 与动量分辨谱函数。DOI：10.1038/ncomms14081

物理展宽和数值展宽怎样分开？

计算曲线变宽时，可先做控制变量：保持展宽核、能量网格、k 点密度、晶胞大小、交换关联泛函和投影设置一致，只改变待研究的结构、温度模型或缺陷浓度。若峰宽仍随物理变量有规律地变化，才有理由追查能带色散、构型分布或相互作用。同一展宽参数是比较峰宽的最低条件。

图6的低能 DOS 在计算中采用恒定展宽参数，并排除了费米–狄拉克分布造成的温度展宽。随模型温度降低，零能峰的 FWHM 仍发生变化，随后低温相出现完整能隙；论文把前者联系到低能态的 Doppler 位移分布。固定数值核后保留下来的系统变化才进入物理解释，模型 DOS 的 FWHM 仍须与实验仪器线宽分开表述。

图6. 洁净与无序模型的低能态密度及零能峰半高全宽。DOI：10.1038/s41535-022-00450-w

电子–声子散射、电子–电子相互作用和有限寿命对应更严格的谱函数 A(k,E) 线宽，需要包含自能的频率依赖；普通基态 DFT 本征值 DOS 没有自动包含这些寿命效应。静态无序超胞产生的峰展宽，多半反映不同局域环境的能级分布。“物理展宽”要说明所用理论量与散射来源，不能从一条平滑 DOS 曲线反推寿命。

数值设置也能同时改写峰位、高度和宽度。TiNiSn 的对照计算显示，稀疏 3×3×3 k 网格配合高斯展宽会产生多处伪峰；即使增至 21×21×21 网格，高斯核仍会抹平带边，而同一密网格下的四面体法保留更尖锐的谱结构。增加 k 点与减小展宽应分别做收敛检查，两项不能互相替代。

图7. TiNiSn 采用四面体法及不同 k 点网格高斯展宽得到的态密度。DOI：10.1016/j.mtelec.2022.100002

实际判读可报告一组能复核的量：共同能量参考下的 E_max 或 E_c，明确能量窗口内的积分面积，适用时报告 FWHM，并附上归一化方式、k 网格与展宽核。峰位给出态集中在哪段能量，峰高给出该点的态密度，峰宽给出谱权重铺开的能量范围；材料机理由三者的联动及其配套电子结构证据限定。