4、BJT能带图（平衡态）：发射区、基区、集电区的能带连接与费米能级对齐

4.4 发射区与基区的能带连接

发射区（N⁺）与基区（P）形成的发射结是BJT工作的关键。在平衡态下，该结的能带连接具有以下特征：

1. 能带弯曲方向：从N⁺发射区到P基区，E_C和E_V向上弯曲（能量增加），形成电子势垒。
2. 势垒高度：势垒高度等于内建电势qV_bi,EB，其值由发射区和基区的掺杂浓度决定：

   V_bi_EB = (kT/q) * ln(N_DE * N_AB / n_i^2)

   其中N_DE和N_AB分别为发射区和基区的掺杂浓度，n_i为本征载流子浓度。
3. 空间电荷区：由于发射区掺杂远高于基区，空间电荷区主要向基区一侧扩展，基区一侧的耗尽层宽度远大于发射区一侧。

4.5 基区与集电区的能带连接

基区（P）与集电区（N）形成的集电结在平衡态下同样形成PN结，但其特征与发射结有所不同：

4.6 平衡态能带图的参数提取

从平衡态能带图中，我们可以提取以下关键器件物理参数：

• 内建电势V_bi：通过测量结区E_C或E_V的弯曲量直接获得。
• 掺杂浓度：结合中性区E_C相对于E_F的位置，利用公式N_D = N_C exp[-(E_C-E_F)/kT] 或 N_A = N_V exp[-(E_F-E_V)/kT] 计算。
• 耗尽层宽度：通过能带弯曲区域的水平跨度估算。
• 基区中性区宽度：从发射结耗尽层边缘到集电结耗尽层边缘的距离，是决定BJT电流增益的关键几何参数。

4.7 平衡态能带图的数值计算示例

以下Python代码演示了如何计算NPN型BJT在平衡态下的能带图关键参数：

import numpy as np

# 物理常数
k = 8.617e-5  # eV/K
T = 300       # K
q = 1.602e-19 # C
epsilon_si = 11.7 * 8.854e-14 # F/cm

# 硅材料参数
Eg = 1.12     # eV
Nc = 2.8e19   # cm^-3
Nv = 1.04e19  # cm^-3
ni = 1.5e10   # cm^-3

# 掺杂浓度 (cm^-3)
N_DE = 1e19   # 发射区 N+
N_AB = 1e17   # 基区 P
N_DC = 1e16   # 集电区 N

# 计算费米能级相对于本征费米能级的位置
# 发射区 (N+)
phi_F_E = k * T * np.log(N_DE / ni)
# 基区 (P)
phi_F_B = -k * T * np.log(N_AB / ni)
# 集电区 (N)
phi_F_C = k * T * np.log(N_DC / ni)

print("=== 费米势 (相对于本征费米能级) ===")
print(f"发射区: {phi_F_E:.3f} V")
print(f"基区:   {phi_F_B:.3f} V")
print(f"集电区: {phi_F_C:.3f} V")

# 计算内建电势
V_bi_EB = phi_F_E - phi_F_B
V_bi_CB = phi_F_B - phi_F_C

print("\n=== 内建电势 ===")
print(f"发射结 V_bi_EB: {V_bi_EB:.3f} V")
print(f"集电结 V_bi_CB: {V_bi_CB:.3f} V")

# 计算导带底相对于费米能级的位置
Ec_Ef_E = k * T * np.log(Nc / N_DE)  # 发射区
Ec_Ef_B = Eg - k * T * np.log(Nv / N_AB)  # 基区
Ec_Ef_C = k * T * np.log(Nc / N_DC)  # 集电区

print("\n=== 导带底 E_C 相对于 E_F 的位置 (eV) ===")
print(f"发射区: {Ec_Ef_E:.3f} eV")
print(f"基区:   {Ec_Ef_B:.3f} eV")
print(f"集电区: {Ec_Ef_C:.3f} eV")

# 计算耗尽层宽度 (单边突变结近似)
# 发射结 - 耗尽层主要向基区扩展
W_dep_EB_base = np.sqrt(2 * epsilon_si * V_bi_EB / (q * N_AB))
# 集电结 - 耗尽层主要向集电区扩展
W_dep_CB_collector = np.sqrt(2 * epsilon_si * V_bi_CB / (q * N_DC))

print("\n=== 耗尽层宽度 (cm) ===")
print(f"发射结在基区侧的耗尽层宽度: {W_dep_EB_base:.2e} cm")
print(f"集电结在集电区侧的耗尽层宽度: {W_dep_CB_collector:.2e} cm")

# 输出结果示例
# === 费米势 (相对于本征费米能级) ===
# 发射区: 0.527 V
# 基区:   -0.407 V
# 集电区: 0.347 V
# 
# === 内建电势 ===
# 发射结 V_bi_EB: 0.934 V
# 集电结 V_bi_CB: 0.754 V
# 
# === 导带底 E_C 相对于 E_F 的位置 (eV) ===
# 发射区: 0.047 eV
# 基区:   0.953 eV
# 集电区: 0.147 eV
# 
# === 耗尽层宽度 (cm) ===
# 发射结在基区侧的耗尽层宽度: 1.10e-05 cm
# 集电结在集电区侧的耗尽层宽度: 3.13e-05 cm

4.8 平衡态能带图与器件物理参数的关系总结

平衡态能带图是理解BJT所有非平衡态工作模式的基础。下表总结了能带图关键特征与器件物理参数之间的对应关系：