4. 能带结构模型选择

在 Sentaurus Device 中,能带结构模型决定了载流子在电场与浓度梯度下的输运参数如何被计算。它并不直接求解薛定谔方程,而是通过态密度、群速度与有效质量三个中间量,把能带形状的物理信息传递给漂移-扩散或流体动力学方程。因此,模型选择的核心问题是:在当前材料、尺寸与偏置条件下,抛物带近似是否仍然成立?若不成立,误差出现在哪里,又该用哪种修正来补偿。

一、三种模型的物理含义与适用边界

最基础的是抛物带有效质量模型。它假设 E-k 关系为 E = ħ²k²/2m*,态密度在三维下正比于 √(E-Ec),有效质量 m* 为常数。对硅、锗在低场、低载流子能量(约低于 0.3 eV)时,这一近似误差通常在几个百分点以内,计算量最小,是默认起点。

当载流子被高电场加热,或材料本身导带底附近曲率变化剧烈时,E-k 曲线偏离抛物线,态密度不再按 √E 增长,群速度也不再线性于 k。此时需要非抛物带模型。Sentaurus 中常用 Kane 型色散关系:

E(1 + αE) = ħ²k² / (2m*)
α 为非抛物因子,Si 电子约 0.5 eV⁻¹,GaAs 电子约 0.64 eV⁻¹

非抛物带修正会压低高能态密度、改变饱和速度与碰撞电离系数,对 0.5 eV 以上能量的载流子尤为关键。

第三类是量子修正模型。当器件尺寸缩小到沟道或势垒宽度与载流子德布罗意波长可比(硅中电子约 5~10 nm)时,即使能带本身接近抛物线,载流子在势阱中的量子限制也会使电荷重心偏离界面、等效开启电压升高、有效势垒变厚。Sentaurus 提供密度梯度(Density Gradient)与改进的局域密度近似(eQuantumPotential)两种主流修正。它们不改变 E-k 关系,而是在漂移-扩散方程中引入量子势 Λ:

Λ = - (ħ² / 6m*q) · (∇²√n) / √n
n 为载流子浓度,q 为电子电荷

量子势把载流子从势垒顶部“推开”,从而在经典框架下复现量子限制效应。

二、不同材料与尺寸下的选择依据

选择不是凭经验拍板,而是由三个可量化判据共同决定:载流子平均能量、器件特征尺寸、材料非抛物因子。

  • 硅基器件,特征尺寸 > 50 nm,低场:抛物带有效质量模型足够,态密度质量取 1.08 m₀(电子)、0.81 m₀(空穴)。
  • 硅基器件,特征尺寸 10~50 nm,或漏端电场 > 10⁵ V/cm:需开启非抛物带修正,α 取 0.5 eV⁻¹,同时保留有效质量模型作为基础。
  • 特征尺寸 < 10 nm 的 FinFET、纳米片或超薄体:必须叠加量子修正。此时即使非抛物带已开启,仍会因忽略量子限制而高估开态电流、低估阈值电压。
  • GaAs、InGaAs 等高迁移率材料:非抛物因子普遍大于硅,α 常取 0.6~0.7 eV⁻¹,且 Γ 谷与 L 谷能量间隔小,高场下需考虑谷间转移,非抛物带模型应作为默认选项。
  • SiC、GaN 等宽禁带材料:击穿场强高,载流子能量可超过 1 eV,非抛物带修正是碰撞电离系数准确性的前提;量子修正则在栅极附近势垒极薄时启用。

三、Sentaurus 中的实际配置与判据

在 Physics 段中,能带模型通过以下语句组合启用:

Physics {
  EffectiveIntrinsicDensity( OldSlotboom )   # 或 NewSlotboom
  eQuantumPotential                          # 量子修正
  hQuantumPotential
  Bandgap narrowing
  Recombination( SRH Auger )
}
Mobility( DopingDependence HighFieldSaturation )
EffectiveMass( NonParabolic( alpha = 0.5 ) )

判断是否需要量子修正,可在同一结构上做两次仿真:一次关闭量子势,一次开启,比较阈值电压漂移量 ΔVth。若 ΔVth > 30 mV,或反型层电荷重心与界面距离偏差超过 1 nm,则量子修正不可省略。判断非抛物带是否必要,则观察高场区载流子平均能量:若超过 0.4 eV,关闭非抛物带会使饱和速度偏高 10% 以上。

实用判据:先以抛物带 + 有效质量跑基线,提取漏端平均能量与沟道特征长度;若能量 > 0.4 eV 则加非抛物带,若特征长度 < 10 nm 则加量子修正。两者可同时启用,互不冲突。

四、模型组合与结果判读

三种模型并非互斥,而是分层叠加:有效质量提供基础态密度,非抛物带修正高能尾部,量子修正调整空间电荷分布。实际校准中,建议按以下顺序固定变量:先固定有效质量与材料参数,再扫描 α 观察饱和速度与电离系数,最后开启量子势并对比 C-V 曲线与阈值电压。若开启量子修正后亚阈值摆幅反而变差,通常说明量子势系数过大或网格在界面处过粗,应检查界面附近网格是否小于 0.5 nm。

能带结构模型选择决策图 依据载流子平均能量、特征尺寸与材料非抛物因子逐层判断 开始:确定材料与结构 判据 1:载流子平均能量 是否 > 0.4 eV? 否 抛物带有效质量 m* 常数,√E 态密度 是 启用非抛物带 Kane 型,α 修正 判据 2:器件特征尺寸 是否 < 10 nm? 否 经典漂移-扩散 无需量子势 是 叠加量子修正 Density Gradient 判据 3:材料非抛物因子 α 是否 > 0.5 eV⁻¹? 否 Si / Ge 基础参数 α ≈ 0.5 eV⁻¹ 是 GaAs / InGaAs α ≈ 0.6~0.7 eV⁻¹ 最终模型组合 有效质量(基础) + 非抛物带(高能修正) + 量子势(尺寸修正) 三者分层叠加,互不冲突,按判据逐层启用

图中三条判据分别对应能量、尺寸与材料三个维度,实际使用时不必严格按顺序,但建议先判断尺寸,因为量子修正对阈值电压的影响往往比非抛物带更直接。若器件同时满足高场与纳米尺寸,则三种模型全部启用,此时应重点检查网格密度与量子势系数是否匹配,避免因数值误差掩盖物理修正效果。

常见误区:把非抛物带模型当作量子修正使用。前者修正的是能带色散关系,后者修正的是载流子空间分布。在 5 nm 沟道中只开非抛物带而不开量子势,阈值电压仍会偏低约 50~80 mV。

最后需要强调,模型选择不是越复杂越好。每增加一层修正,都会引入新的参数(α、量子势系数、有效质量张量),这些参数需要与实验数据或更高精度计算校准。若缺乏校准依据,优先保证基础有效质量与材料参数准确,再逐步叠加修正,并通过 C-V、I-V 与饱和速度三组曲线交叉验证。