态密度峰越高越好吗？统计含义与输运性能的关联性辨析

说明：本文华算科技主要介绍态密度峰的统计含义，以及峰高与电荷输运、热电响应和自旋极化之间的关系。

态密度的定义与峰形来源

在周期晶体中，电子态由能带编号 n、波矢 k 和自旋共同标记。态密度 D(E) 将这些态按能量统计：D(E)dE 是窄能量区间内的态数，常用单位为 states·eV⁻¹·cell⁻¹。第一性原理计算通常从 Kohn–Sham 本征值 ε_nk 出发，按布里渊区权重 w_k 求和，写成 D(E)=Σ_nkw_kδ(E−ε_nk)；这里把自旋计入能带编号。

若一片 k 空间中的能带色散很小，许多电子态集中在相近能量，就会产生窄而高的 DOS 峰。能带的局部平坦、简并及极值附近的形状都能改变峰形。GaS 单层的价带顶偏离 Γ 点，沿高对称线形成中央略低、两侧较高的轮廓，价带顶附近出现尖锐的态密度特征；增加层数后，色散随层间耦合改变。

图1把能量沿竖轴排列，DOS 沿横轴展开，峰的“高低”实际表现为横向伸出的幅度。能量零点是价带顶 E_VBM，并非费米能级。原图按每晶胞归一化，不同厚度模型的晶胞含有不同数量的原子；比较每个 GaS 化学式单元贡献时，应先换算态数单位。原子投影曲线还取决于投影区域与轨道定义，局部峰高无法直接替代全体系态数。

GaS 块体在 Γ 点附近接近抛物线色散，单层的顶部则明显展宽。图2用 k²、k⁴ 形式拟合两种轮廓；单层拟合线没有重现 Γ 点的凹陷。这一偏差涉及价带顶附近的态数分布，精细峰形须由实际能带求和获得，用一条平顶拟合曲线无法恢复全部细节。

绘图时，离散能级常用归一化高斯函数等展宽。对于同一组能级，减小展宽可使孤立峰变窄、峰顶升高，全能区积分的态数保持不变。稀疏 k 点也可能产生不收敛的尖峰。物理色散变化与数值展宽变化都能改变峰高，只有固定归一化方式，并确认峰形随采样加密已收敛，才能比较两次计算中的能量分布。

区间积分 ∫D(E)dE 给出态数，占据电子数则由 ∫D(E)f(E)dE 计算，f(E) 是费米分布。深价带中的高峰通常已经占满，导带中的高峰可能仍为空态。自旋分辨图把某个自旋分量画在负轴时，负号只是显示约定；计算电子数应按两种自旋各自的非负态密度积分。金属的低温电荷输运主要涉及费米能级附近可改变占据的电子态。

高峰中的电子态能贡献多少电流？

在能带描述和弛豫时间近似下，电导率沿 i 方向的分量为 σ_ii=e²∫Ξ_ii(E)(−∂f/∂E)dE。输运分布 Ξ_ii 将单位体积内的电子态按速度分量平方 v_i² 与弛豫时间 τ 加权。DOS 只统计态数；Ξ_ii 还区分每个态沿测量方向的运动和散射。

群速度满足 v_nk=∇_kε_nk/ℏ。能带变平时，DOS 可以增大，该方向的群速度却可能减小；完全平坦的理想孤立能带在普通带内输运模型中没有群速度贡献。真实材料可同时具有平坦和陡峭的色散区，也可由多条能带共同导电。LaPt₂B 的多个费米面就具有不同的速度方向，图3a的颜色标示速度大小。

图3中的空穴型 α 费米面主要贡献沿 c 轴的正热电响应，近圆柱形的电子型 γ 费米面主要贡献面内负响应。图3b统计的是Peltier 电导率 α_ii=σ_iiS_ii，其中 S 为 Seebeck 系数；柱形的正负表示热电响应符号。该计算采用300 K及统一的 τ=10⁻¹⁴ s，用于分辨各能带的贡献，未计算不同能带的实际散射寿命。

Seebeck 系数涉及费米能级两侧输运分布的不对称。即使 DOS 在 E_F 附近很高，若两侧对热电流的贡献接近抵消，S 仍可很小。同一块 LaPt₂B 晶体沿 a、c 轴测得的 S 符号相反；旋转测量方向没有改变总 DOS，却改变了各费米面的输运权重。

图4中，c 轴电阻率低于 a 轴，两个方向均随温度升高而增大。电导率、热电系数和功率因子 S²σ是不同的评价量。将高 DOS 称为热电优势时，必须计算 S 与 σ 的乘积效应；若评价热电优值，还涉及电子和晶格热导率。

对于缺陷或声子散射，若费米能级附近增加的电子态可作为允许的散射末态，且散射矩阵元不小，散射率可能随可用末态数增加，抵消态数增多带来的电流贡献。反过来，矩阵元很小或动量、自旋选择限制较强时，相同的 DOS 增幅也可能只有较小影响。DOS 曲线没有包含这些矩阵元，计算 τ 或测量迁移率才能检验这一项。

费米能级处的高峰为何可能引起磁性？

带有周期性边缘凸起的折线形石墨烯纳米带，在横向电场下形成靠近带隙的窄带。图5c的带宽随电场呈非单调变化：价带与导带的最小带宽出现在不同电场，分别约为0.78和0.82 V·Å⁻¹，对应带宽约0.5和1.4 meV。

图5e、f中，价带和导带的局域态密度分布在相反侧的边缘凸起，这里的空间图对相应能带作了能量积分。电场改变碳原子的局部势，电子态在边缘形成局域分布，带宽随之缩小。未掺杂时窄带可处于全满或全空状态；电荷掺杂改变填充程度后，费米能级才会经过相关高峰。

在 Stoner 模型中，自旋极化带来交换能降低，也带来重新占据能级的能量代价。Iρ(E_F)>1描述非磁态对微小自旋极化失稳的条件，I 为交换作用参数，ρ 来自非磁参照。两者必须采用配套的晶胞与自旋归一化约定；把单自旋 DOS 换成双自旋总 DOS 后，原来的系数表达也须调整。高峰只增大其中的态密度因子。

在0.82 V·Å⁻¹电场、每晶胞增加一个电子的模型中，图6a的非磁参照在 E_F 附近有尖峰；图6b中，有限磁矩态的总能低于零磁矩态。后一个结果检验了自旋极化的能量收益。形成磁态后，自旋分裂会重新分配能级和占据，非磁态的高峰可能分裂或偏离 E_F，磁态 DOS 应作为新的自洽结果分析。

该纳米带的磁矩随掺杂先增大再减小，在每晶胞约增加或移去一个电子时接近1 μ_B，到约两个电子时消失。部分填充与填满窄带产生不同磁矩，掺杂量增加并未维持单调的磁性增强。图7同时给出两类掺杂下的空间自旋密度和自旋分辨能带，分别检验磁矩分布与能级分裂。

图7b、e的自旋密度是两种自旋电子密度之差，主要分布在一侧边缘凸起；掺杂符号改变时，磁矩分布转向另一侧。这些结果属于给定外电场、填充和结构下的自旋极化计算。若讨论纳米带能否维持有限温度下的长程磁序，还涉及磁交换、磁各向异性能和一维热涨落；上述能带与磁矩图没有给出居里温度。