04 · 陷阱辅助复合理论与扩展界面模型
系列文章:太阳电池晶体硅界面钝化研究
本文定位:从"陷阱能级"这个最基础概念开始,用统计力学严格推导陷阱辅助复合公式,然后推导论文第 3 章的扩展 SRH 模型
前置知识:建议先读 01 · 半导体物理基础 和 02 · 表面与界面化学
写作原则:每个公式都从最基本的物理概念开始推导,不跳步;每个符号都明确物理意义
目录
- 从单个陷阱开始:陷阱能级的物理图像
- 陷阱的占据率:费米-狄拉克统计
- 陷阱辅助复合的 4 个基本过程
- 稳态条件与 R_SRH 公式推导
- 本征费米能级 E_i 与特征浓度 n₁, p₁
- 电子/空穴寿命 τ_n, τ_p 的物理含义
- 从体复合到表面复合:边界条件
- 论文扩展模型:5 个新公式(3.27-3.31)
- 反型层条件:Q_f 临界值计算
- D_it 谱的物理意义:为什么禁带中央最致命
- 从 D_it 到 S_eff:钝化设计的核心公式
1 · 从单个陷阱开始:陷阱能级的物理图像
1.1 什么是"陷阱"?
陷阱(trap) = 禁带中的孤立能级,可以临时捕获载流子。
物理起源:
- 晶格缺陷:空位(vacancy)、间隙原子(interstitial)、位错(dislocation)
- 杂质原子:金属杂质(Fe、Cu、Au)、掺杂剂(P、B)
- 界面态:Si/SiO₂、Si/Al₂O₃ 界面的悬挂键
- 复合中心:本论文关心的 P_b0、P_b1
1.2 陷阱的能级位置
每个陷阱在能带图中用一个水平短线表示,位于禁带中的某个能量 E_T:
Ec(导带底)
╱‾‾‾‾‾‾‾‾‾‾‾‾╲
╱ ╲
╱ ╲
│ ET │ ← 陷阱能级
│ │
│ │
╲ ╱
╲ ╱
╲____________╱
Ev(价带顶)
E_T 在禁带中的位置至关重要:
- E_T 靠近 E_C:陷阱主要捕获电子(类受主)
- E_T 靠近 E_V:陷阱主要捕获空穴(类施主)
- E_T 在禁带中央附近:对复合最致命(下面会详细解释)
1.3 陷阱的物理参数
每个陷阱有几个关键参数:
| 参数 | 物理含义 | 单位 | 典型值 |
|---|---|---|---|
| E_T | 陷阱能级 | eV | 0.0-0.5 eV(Si) |
| N_T | 陷阱密度 | cm⁻³(体)或 cm⁻²(表面) | 10¹⁰-10¹⁵ |
| σ_n | 电子捕获截面 | cm² | 10⁻¹⁵-10⁻¹³ |
| σ_p | 空穴捕获截面 | cm² | 10⁻¹⁵-10⁻¹³ |
| v_th | 载流子热速度 | cm/s | ~ 10⁷(300K) |
捕获截面 = 陷阱"捕捉"载流子的"靶面积"。如果一个陷阱"看上去"面积 σ,载流子撞上去就被捕获。
类比:陷阱像一张蜘蛛网,σ 是网的面积,载流子像飞虫,撞上网就被捕获。
1.4 单个陷阱的 4 个基本过程
对于一个被占据/未被占据的陷阱,有 4 个可能的过程:
① 电子捕获 ② 电子发射
e⁻ ↓ ↑ e⁻
↓ ↓ ↑ ↑
════╪══╪════════ ════╪══╪════════
│ │ │ │
│ │ │ │
═══════════════ ═══════════════
(陷阱空着, (陷阱被电子占据,
等着捕获电子) 准备释放电子)
③ 空穴捕获 ④ 空穴发射
h⁺ ↓ ↑ h⁺
↓ ↓ ↑ ↑
════╪══╪════════ ════╪══╪════════
│ │ │ │
│ │ │ │
═══════════════ ═══════════════
(陷阱被电子占据, (陷阱空着,
等着捕获空穴) 准备释放空穴)
| 过程 | 含义 | 速率 |
|---|---|---|
| ① 电子捕获 | 自由电子被陷阱捕获 | R_capture,n = σ_n × v_th × n × N_T × (1 - f_T) |
| ② 电子发射 | 陷阱释放电子到导带 | R_emission,n = σ_n × v_th × n₁ × N_T × f_T |
| ③ 空穴捕获 | 自由空穴被陷阱捕获(陷阱中的电子落回价带) | R_capture,p = σ_p × v_th × p × N_T × f_T |
| ④ 空穴发射 | 陷阱从价带捕获电子(等价于释放空穴) | R_emission,p = σ_p × v_th × p₁ × N_T × (1 - f_T) |
注意:
- f_T = 陷阱被电子占据的概率(0 ~ 1)
- (1 - f_T) = 陷阱空着的概率
- n = 自由电子浓度
- p = 自由空穴浓度
2 · 陷阱的占据率:费米-狄拉克统计
2.1 占据率公式
陷阱被电子占据的概率 f_T 由类似费米-狄拉克分布给出:
$$f_T = \frac{1}{1 + \frac{g_0}{g_1} \exp\left(\frac{E_T - E_F}{kT}\right)}$$
其中:
- g_0, g_1:陷阱空着/被占据时的简并度(通常 g_0 = 1, g_1 = 2 for 自旋)
- E_T:陷阱能级
- E_F:费米能级
特殊情况(本论文用):
如果只考虑自旋简并,类受主型陷阱的占据率:
$$f_T = \frac{1}{1 + \frac{1}{2} \exp\left(\frac{E_T - E_F}{kT}\right)}$$
2.2 占据率的物理含义
| 情况 | f_T | 含义 |
|---|---|---|
| E_F >> E_T | f_T ≈ 1 | 陷阱被电子填满(陷阱能级远低于费米能级) |
| E_F = E_T | f_T ≈ 0.67 | 半填充(具体取决于简并度) |
| E_F << E_T | f_T ≈ 0 | 陷阱空着(陷阱能级远高于费米能级) |
2.3 占据率随 E_F 变化
f_T(占据率)
1.0 ┤●●●●●●●●●
│ ╱
│ ╱
0.5 ┤───────●──────
│ ╱
│ ╱
0.0 ┤●●●╱
└─────────────────→ E_F - E_T
负 0 正
↑
此处占据率
急剧变化
关键观察:当 E_F ≈ E_T 时,f_T 急剧变化——这是费米能级"扫过"陷阱能级时的行为。
3 · 陷阱辅助复合的 4 个基本过程(详细推导)
3.1 过程 ①:电子捕获(自由电子 → 陷阱)
物理图像:一个导带自由电子在陷阱附近"撞上"陷阱,被陷阱捕获。
前提条件:
- 陷阱空着(否则电子去不了)
- 自由电子有足够热运动能接近陷阱
速率公式推导:
单位时间、单位体积内,有多少电子被陷阱捕获?
- 单个陷阱的捕获率:σ_n × v_th(单位时间,陷阱"扫过"面积为 σ_n × v_th 的空间)
- 自由电子密度:n(单位体积内的电子数)
- 陷阱密度:N_T(单位体积内的陷阱数)
- 陷阱空着的概率:(1 - f_T)
所以单位体积的电子捕获速率:
$$R_{capture,n} = \sigma_n \cdot v_{th} \cdot n \cdot N_T \cdot (1 - f_T)$$
量纲检查:
- σ_n:cm²
- v_th:cm/s
- n:cm⁻³
- N_T:cm⁻³
- (1 - f_T):无量纲
- 总体:cm² × cm/s × cm⁻³ × cm⁻³ = cm⁻³/s ✓ (单位体积内,单位时间的捕获次数)
3.2 过程 ②:电子发射(陷阱 → 导带)
物理图像:陷阱中的电子获得热能,逃出陷阱回到导带。
前提条件:
- 陷阱有电子占据
- 电子能获得足够能量克服陷阱束缚
速率公式:
陷阱中电子的发射率与"被电子占据的陷阱数"和"电子能逃出来的概率"成正比。
电子能逃出来的概率(类比"反应速率常数")可以用细致平衡原理推导——平衡时,过程 ① = 过程 ②(否则陷阱会无限填充或清空)。
平衡时:
- n = n_0(平衡电子浓度)
- f_T 由费米-狄拉克分布决定
- 平衡条件:过程 ① 速率 = 过程 ② 速率
平衡时 n_0 = N_C × exp(-(E_C - E_F)/kT),f_T = 1/(1 + 0.5 exp((E_T - E_F)/kT))。
代入细致平衡条件,可以推出:
$$R_{emission,n} = \sigma_n \cdot v_{th} \cdot n_1 \cdot N_T \cdot f_T$$
其中 n_1 是"特征电子浓度":
$$n_1 = N_C \cdot \exp\left(-\frac{E_C - E_T}{kT}\right)$$
n_1 的物理意义:当 E_F = E_T 时,电子浓度等于 n_1。n_1 反映"陷阱束缚电子的强度"——E_T 越深(远离 E_C),n_1 越小,陷阱越难释放电子。
3.3 过程 ③:空穴捕获(陷阱 → 价带)
物理图像:陷阱中已占据的电子"掉"到价带,与价带空穴复合。这就是空穴被"中和"。
物理上,等价于:陷阱中已占据的电子"逃"到价带(同时价带空穴被"消除")。
速率公式:
$$R_{capture,p} = \sigma_p \cdot v_{th} \cdot p \cdot N_T \cdot f_T$$
注意:
- p 是自由空穴浓度
- f_T 是陷阱被电子占据的概率(必须先有电子才能"释放到价带")
- R 是"电子落回价带"的速率,等价于"空穴被中和"的速率
3.4 过程 ④:空穴发射(价带 → 陷阱)
物理图像:价带电子获得热能跳到陷阱(等价于"陷阱从价带捕获电子"),在价带留下空穴。这就是空穴被"发射"。
$$R_{emission,p} = \sigma_p \cdot v_{th} \cdot p_1 \cdot N_T \cdot (1 - f_T)$$
其中 p_1 是"特征空穴浓度":
$$p_1 = N_V \cdot \exp\left(-\frac{E_T - E_V}{kT}\right)$$
3.5 4 个速率的汇总
| 过程 | 公式 | 含义 |
|---|---|---|
| ① 电子捕获 | σ_n v_th n N_T (1 - f_T) | 自由电子 → 陷阱 |
| ② 电子发射 | σ_n v_th n_1 N_T f_T | 陷阱 → 自由电子 |
| ③ 空穴捕获 | σ_p v_th p N_T f_T | 陷阱(电子) → 价带空穴中和 |
| ④ 空穴发射 | σ_p v_th p_1 N_T (1 - f_T) | 价带 → 陷阱(陷阱从价带"拉"电子) |
4 · 稳态条件与 R_SRH 公式推导
4.1 稳态条件:陷阱占据率不变
稳态(steady state) = 系统中陷阱的占据率不随时间变化(陷阱既不无限填充也不无限清空)。
数学表达:
$$\frac{df_T}{dt} = 0$$
这意味着:单位时间内,进入陷阱的电子数 = 离开陷阱的电子数。
进入:过程 ① 电子捕获 + 过程 ④ 空穴发射
离开:过程 ② 电子发射 + 过程 ③ 空穴捕获
稳态方程:
$$R_{capture,n} + R_{emission,p} = R_{emission,n} + R_{capture,p}$$
代入:
$$\sigma_n v_{th} n N_T (1-f_T) + \sigma_p v_{th} p_1 N_T (1-f_T) = \sigma_n v_{th} n_1 N_T f_T + \sigma_p v_{th} p N_T f_T$$
两边除以 N_T:
$$\sigma_n v_{th} n (1-f_T) + \sigma_p v_{th} p_1 (1-f_T) = \sigma_n v_{th} n_1 f_T + \sigma_p v_{th} p f_T$$
两边除以 v_th:
$$\sigma_n n (1-f_T) + \sigma_p p_1 (1-f_T) = \sigma_n n_1 f_T + \sigma_p p f_T$$
解出 f_T:
$$f_T = \frac{\sigma_n n + \sigma_p p_1}{\sigma_n (n + n_1) + \sigma_p (p + p_1)}$$
4.2 SRH 复合速率 R_SRH
复合速率 = 单位时间、单位体积内,一个电子-空穴对通过陷阱消失。
每个复合事件需要:
- 过程 ①(电子被捕获)或过程 ③(空穴被中和)
- 两个过程都发生,才完成一次 e-h 复合
净复合速率:
$$R_{SRH} = R_{capture,n} - R_{emission,n} = R_{capture,p} - R_{emission,p}$$
(两个表达式应该相等,这是细致平衡的另一个表述)
代入过程 ①:
$$R_{SRH} = R_{capture,n} - R_{emission,n} = \sigma_n v_{th} N_T [\sigma_n n (1-f_T) - \sigma_n n_1 f_T]$$
化简:
$$R_{SRH} = \sigma_n v_{th} N_T [n (1-f_T) - n_1 f_T]$$
代入 f_T 的表达式,经过代数化简:
$$\boxed{R_{SRH} = \frac{N_T \sigma_n \sigma_p v_{th} (np - n_i^2)}{\sigma_n (n + n_1) + \sigma_p (p + p_1)}}$$
这就是著名的陷阱辅助复合速率公式。
4.3 简化为常用形式
定义载流子寿命:
$$\tau_p = \frac{1}{\sigma_p v_{th} N_T}, \quad \tau_n = \frac{1}{\sigma_n v_{th} N_T}$$
代入 R_SRH 公式:
$$R_{SRH} = \frac{np - n_i^2}{\tau_p (n + n_1) + \tau_n (p + p_1)}$$
这就是最常见的陷阱辅助复合公式。
4.4 特殊情况下的简化
情况 1:本征载流子浓度附近(低注入)
n ≈ p ≈ n_i,且 E_T 在禁带中央 → n_1 ≈ p_1 ≈ n_i
$$R_{SRH} \approx \frac{n_i^2 - n_i^2}{\tau_p (n_i + n_i) + \tau_n (n_i + n_i)} = 0$$
但加上非平衡载流子 Δn 后:
$$R_{SRH} \approx \frac{\Delta n}{\tau_p + \tau_n}$$
如果 σ_n ≈ σ_p(对称陷阱),则:
$$R_{SRH} \approx \frac{\Delta n}{2\tau_0} \quad \text{其中 } \tau_0 = \frac{1}{\sigma v_{th} N_T}$$
情况 2:陷阱在禁带中央(E_T = E_i)
此时 n_1 = p_1 = n_i,禁带中央的陷阱最致命——后面 10 节会详细解释。
5 · 本征费米能级 E_i 与特征浓度 n₁, p₁
5.1 什么是 E_i?
E_i = 本征费米能级 = 没有任何掺杂时硅的费米能级。
对于本征硅:
- n = p = n_i
- E_F = E_i = 禁带中央偏下一点(因为 N_V ≠ N_C,导致 E_i 不在正中央)
对硅在 300K:
- E_g = 1.12 eV
- E_i ≈ E_V + 0.56 eV(近似禁带中央)
5.2 特征浓度 n_1 的物理含义
n_1 = N_C × exp(-(E_C - E_T)/kT)
解释:
- 如果 E_F = E_T,则 n = n_1
- n_1 反映"陷阱与导带的距离"
- E_T 越靠近 E_C,n_1 越大,陷阱越"像导带"(电子容易上去也容易下来)
5.3 特征浓度 p_1 的物理含义
p_1 = N_V × exp(-(E_T - E_V)/kT)
解释:
- 如果 E_F = E_T,则 p = p_1
- p_1 反映"陷阱与价带的距离"
- E_T 越靠近 E_V,p_1 越大,陷阱越"像价带"
5.4 n_1 × p_1 = n_i²
重要关系:
$$n_1 \cdot p_1 = N_C N_V \exp\left(-\frac{E_C - E_V}{kT}\right) = n_i^2$$
意义:n_1 和 p_1 不是独立的——它们相乘总是等于 n_i²。
特例:E_T = E_i → n_1 = p_1 = n_i(陷阱在禁带中央时,两个特征浓度都等于本征浓度)
6 · 电子/空穴寿命 τ_n, τ_p 的物理含义
6.1 寿命的定义
$$\tau_n = \frac{1}{\sigma_n v_{th} N_T}, \quad \tau_p = \frac{1}{\sigma_p v_{th} N_T}$$
τ_n 的物理含义:如果只有电子-陷阱相互作用(没有空穴),电子的平均寿命。
推导:
如果没有空穴,只有电子-陷阱相互作用,电子被陷阱捕获的速率:
$$R_{capture,n} = \sigma_n v_{th} n N_T$$
每个电子的"死亡率":
$$\frac{R_{capture,n}}{n} = \sigma_n v_{th} N_T = \frac{1}{\tau_n}$$
所以 τ_n 就是电子被陷阱"抓住"的平均时间(在没有空穴的情况下)。
6.2 实际复合寿命与 τ_n, τ_p 的关系
在实际陷阱辅助复合中,电子和空穴都被陷阱捕获:
$$R_{SRH} = \frac{np - n_i^2}{\tau_p (n + n_1) + \tau_n (p + p_1)}$$
有效寿命:
$$\tau_{eff} = \frac{\Delta n}{R_{SRH}} = \frac{\Delta n (\tau_p (n + n_1) + \tau_n (p + p_1))}{np - n_i^2}$$
当陷阱在禁带中央(n_1 = p_1 = n_i)且在低注入下(n ≈ p >> n_i):
$$\tau_{eff} \approx \tau_n + \tau_p$$
当陷阱在带边附近(如 n_1 >> n_i)且在低注入下(n < n_1):
$$\tau_{eff} \approx \frac{\tau_p n_1}{n}$$
→ 高浓度 n 下寿命反而短(因为陷阱更容易捕获电子)
7 · 从体复合到表面复合:边界条件
7.1 表面复合的特殊性
表面 = 半无限半导体的边界。
在表面,有几个变化:
- 载流子分布不均:表面有电场(场效应)时,载流子在表面聚集
- 陷阱密度极高:表面态 D_it 通常比体内 N_T 高很多
- 没有"另一侧":表面只能向体内复合(单向)
7.2 表面复合速率 U_s
定义:单位面积、单位时间内,在表面复合的电子-空穴对数。
类比体复合:
$$U_s = \frac{n_s p_s - n_i^2}{\frac{n_s + n_1}{s_p} + \frac{p_s + p_1}{s_n}}$$
其中:
- n_s, p_s:表面电子/空穴浓度(不是体内的!)
- s_n, s_p:表面复合速度(单位 cm/s,类比 v_th × σ × N_T 但单位是 cm/s)
7.3 表面复合速度 S
定义:单位时间内,流过表面单位面积的载流子数。
$$S = \frac{U_s}{\Delta n_s}$$
其中 Δn_s 是表面过剩载流子浓度。
简化情况(陷阱在禁带中央,场钝化强):
$$S \approx \frac{\pi kT D_{it}}{q^2 N_A L_D}$$
其中:
- D_it:界面态密度(cm⁻²·eV⁻¹)
- N_A:掺杂浓度
- L_D:德拜长度
7.4 有效寿命与 S_eff
对于厚度 W 的双面钝化硅片:
$$\frac{1}{\tau_{eff}} = \frac{1}{\tau_{bulk}} + \frac{2S_{eff}}{W}$$
关键洞见:S_eff 越小,τ_eff 越大——钝化的目标就是把 S_eff 压低。
本论文实测:
- 衬底:500 μm 厚 n 型硅
- PE-ALD Al₂O₃:S_eff ≈ 13.8 cm/s
- T-ALD Al₂O₃:S_eff ≈ 23.5 cm/s
代入:
- 体寿命 1 ms(高质量硅)
- 1/τ_eff = 1/1ms + 2×13.8/0.05 cm = 1 ms + 552 cm/s/cm
如果体寿命足够长,τ_eff ≈ W/(2S_eff):
- PE-ALD:τ_eff ≈ 500/(2×13.8) = 1800 μs ✓ 与论文 1802 μs 吻合
- T-ALD:τ_eff ≈ 500/(2×23.5) = 1064 μs ✓ 与论文 1063.4 μs 吻合
完美!公式和实验完全一致。
8 · 论文扩展模型:5 个新公式(3.27-3.31)
8.1 论文 3.4 节的研究目标
论文 3.4 节的研究目标:推导一个能描述场钝化效应的扩展陷阱辅助模型。
标准 SRH 公式没有考虑场钝化(假设表面是平的、没有电荷)。但实际上,Al₂O₃ 钝化层有大量负电荷,会在表面建立电场,强烈影响载流子分布。
8.2 公式 3.27:表面电荷与表面势的关系
论文建立电荷守恒(Gauss 定律):
$$Q_G + Q_{it}(\psi_s) + Q_f + Q_{ox} = 0$$
其中:
- Q_G:硅表面空间电荷区电荷
- Q_it(ψ_s):界面态电荷(依赖于表面势)
- Q_f:钝化层固定电荷(Al₂O₃ 中的负电荷)
- Q_ox:氧化层陷阱电荷(可移动)
当 Q_ox 可忽略(质量好的钝化层),简化:
$$Q_G = -Q_{it}(\psi_s) - Q_f$$
8.3 公式 3.28:表面势 ψ_s 的数值求解
Q_G 是 ψ_s 的函数(通过泊松方程):
$$Q_G(\psi_s) = \pm \frac{\sqrt{2}\epsilon_s kT}{q L_D} \cdot F(\psi_s)$$
其中 F(ψ_s) 是表面势的复杂函数,没有解析解。
论文 3.4 节用数值方法求解(对不同 ψ_s 算 Q_G,再与 -Q_it - Q_f 比较)。
8.4 公式 3.29-3.30:扩展表面复合速率
在强场效应下,n_s 和 p_s 不再等于体内的 n 和 p。
如果表面是反型层(P 型 Si 表面被负电荷吸引电子):
- 表面电子浓度:n_s >> n(电子积累)
- 表面空穴浓度:p_s << p(空穴被推走)
代入陷阱辅助复合公式:
$$U_s = \frac{n_s p_s - n_i^2}{\frac{n_s + n_1}{s_p} + \frac{p_s + p_1}{s_n}}$$
由于 n_s >> p_s,复合速率主要由 p_s 决定:
$$U_s \approx \frac{n_s p_s}{\frac{n_s}{s_p} + \frac{p_s + p_1}{s_n}}$$
如果 s_p >> s_n(空穴比电子难被捕获),化简:
$$U_s \approx s_n p_s$$
关键结论:在强反型条件下,U_s 与 p_s 成正比——把 p_s 压低,就能压低 U_s。
8.5 公式 3.31:总表面复合
论文将表面复合分解为3 个分量:
$$U_{s,total} = U_{s,field} + U_{s,defect} + U_{s,auger}$$
| 分量 | 物理起源 | 主导条件 |
|---|---|---|
| U_s,field | 场效应相关(电荷分离不足) | 弱场钝化 |
| U_s,defect | 界面态 D_it 直接复合 | 强场 + 弱化学钝化 |
| U_s,auger | 表面高载流子浓度的俄歇复合 | 强反型 + 高载流子 |
论文关键洞见:最佳钝化是"低 D_it + 高 |Q_f|"同时实现——单纯场钝化或单纯化学钝化都不够。
8.6 公式 3.35:表面势对表面电荷的影响(论文图 3.11)
论文 3.4 节末尾的图 3.11 用公式 3.35 计算:
$$\psi_s(Q_G) = \text{泊松方程数值解}$$
论文给出三种表面状态:
| 表面状态 | ψ_s | 载流子 | 复合 |
|---|---|---|---|
| 积累 | 强负偏置 | 一种载流子(空穴)主导 | 极低 |
| 耗尽 | ψ_s ≈ 0 | 电子空穴接近 | 最高 |
| 反型 | 强正偏置 | 一种载流子(电子)主导 | 极低 |
反直觉洞见:耗尽状态(平带附近)是复合最严重的——因为电子和空穴密度相近,陷阱辅助复合最活跃。
这就是为什么场钝化要"推到底"——要么强积累,要么强反型,绝对不能停在平带附近。
9 · 反型层条件:Q_f 临界值计算
9.1 反型层的物理条件
P 型 Si 表面形成反型层 = 表面电子浓度 > 表面空穴浓度。
数学条件:
$$n_s = p_s$$
$$n_i \exp(q\psi_s/kT) = N_A \exp(-q\psi_s/kT)$$
(假设 P 型 Si 掺杂 N_A,表面势 ψ_s 为正时吸引电子)
求解:
$$\exp(2q\psi_s/kT) = N_A / n_i$$
$$\psi_{s,inv} = \frac{kT}{q} \ln(N_A / n_i) = 2\phi_B$$
其中 φ_B = (kT/q) ln(N_A/n_i) 是费米势。
数值示例(N_A = 10¹⁵ cm⁻³):
$$\phi_B = 0.026 \times \ln(10^{15}/10^{10}) = 0.026 \times 11.5 = 0.30 \text{ V}$$
$$\psi_{s,inv} = 2\phi_B = 0.60 \text{ V}$$
9.2 临界 Q_f 的推导
要在 P 型 Si 上建立 ψ_s = 0.60 V 的反型层,需要多少 Q_f?
用泊松方程(简化):
空间电荷区的电荷密度 ρ(x) 在 P 型 Si 中:
- 当 ψ_s < φ_B:ρ = -qN_A(电离受主,负)
- 当 ψ_s > φ_B:ρ = -qN_A + qn(电离受主 + 自由电子,但仍约 -qN_A)
空间电荷区宽度 W:
$$W = \sqrt{\frac{2\epsilon_s \psi_s}{qN_A}}$$
空间电荷区总电荷(单位面积):
$$Q_{sc} = -qN_A W = -\sqrt{2\epsilon_s q N_A \psi_s}$$
Q_f 临界值 = |Q_sc| = √(2ε_s q N_A ψ_s)
数值代入(N_A = 10¹⁵ cm⁻³, ψ_s = 0.6 V, ε_s = 1.04 × 10⁻¹² F/cm):
$$Q_{f,crit} = \sqrt{2 \times 1.04 \times 10^{-12} \times 1.6 \times 10^{-19} \times 10^{15} \times 0.6}$$
$$= \sqrt{2 \times 10^{-16}} = \sqrt{2} \times 10^{-8} \approx 1.4 \times 10^{-8} \text{ C/cm}^2$$
换算为 cm⁻²:
$$Q_{f,crit} = \frac{1.4 \times 10^{-8}}{1.6 \times 10^{-19}} = 8.9 \times 10^{10} \text{ cm}^{-2}$$
所以 Q_f 临界值约 10¹¹ cm⁻²——任何 |Q_f| > 10¹¹ cm⁻² 的钝化层,理论上都可在 P 型 Si 上形成反型层。
9.3 论文实测与临界值对比
| 钝化层 | Q_f(cm⁻²) | vs Q_f,crit(10¹¹) | 反型能力 |
|---|---|---|---|
| SiO₂ | +5×10¹⁰ | < (正电荷方向反) | 不能 |
| SiNₓ | +5×10¹¹ | 足够但方向反 | 不能 |
| T-ALD Al₂O₃ | -5.3×10¹² | 50× | 强 |
| PE-ALD Al₂O₃ | -1×10¹³ | 100× | 极强 |
| a-Si:H(i) | 弱 | — | 弱 |
| Ga₂O₃(< 650°C) | -10¹⁰ | < 临界 | 弱 |
| Ga₂O₃(≥ 650°C) | +10¹⁰ | 反转 | 不能 |
关键洞见:Al₂O₃ 的 Q_f 比 Q_f,crit 大 50-100 倍——所以 Al₂O₃ 不仅能形成反型层,而且形成的反型层非常强(表面电子浓度比体内高 50-100 倍)。
10 · D_it 谱的物理意义:为什么禁带中央最致命
10.1 几个特殊位置的 D_it
禁带中不同位置的陷阱,复合能力差别很大。考虑对称陷阱(σ_n = σ_p = σ)和不同 E_T:
$$R_{SRH}(\Delta n \text{ small}) = \frac{\sigma v_{th} N_T \Delta n}{1 + \cosh[(E_T - E_i)/kT]}$$
当 E_T = E_i(禁带中央):
- cosh(0) = 1
- R_SRH 最大
当 E_T = E_i + 0.2 eV:
- cosh(0.2/0.026) = cosh(7.7) ≈ 1096
- R_SRH 下降 ~ 1000 倍
10.2 数值对比:不同位置陷阱的相对复合能力
| 陷阱位置(E - E_i) | cosh 值 | 相对 R_SRH |
|---|---|---|
| 0 eV(禁带中央) | 1.00 | 1.0(基准) |
| ±0.1 eV | 46 | 0.022 |
| ±0.2 eV | 1096 | 0.0009 |
| ±0.3 eV | 26000 | 4 × 10⁻⁵ |
| ±0.5 eV | 10⁷ | 10⁻⁷ |
关键洞见:禁带中央 ± 0.1 eV 内的陷阱,占复合率的 99%。远离禁带中央的陷阱几乎不参与复合。
10.3 P_b0 中心的双 U 型
P_b0 中心在禁带中有两个能级(类受主 + 类施主),间距 0.7 eV。
如果类受主能级在 E_i + 0.35 eV,类施主能级在 E_i - 0.35 eV:
- 两个能级都在"危险区"内(虽然不在 E_i,但 cosh(0.35/0.026) ≈ 6.6 × 10⁵,相对 R_SRH 仍 ~ 10⁻⁵,不致命)
但如果 P_b0 中心能级有部分分布靠近 E_i(由于重构应变、界面缺陷):
- 这些"中间能级"的 P_b0 中心可以非常致命
- 论文 3.4 节将 P_b0 模拟为两个高斯分布(不是离散能级),考虑到了这种能量分散
10.4 论文的 D_it 拟合(3.4 节,图 3.14)
论文将 Al₂O₃/Si 界面的 D_it 谱拟合为 3 个高斯函数:
| 高斯函数 | 物理起源 | 中心位置 | 用途 |
|---|---|---|---|
| G1 | 类 UT(导带尾) | E_i + 0.4-0.5 eV | 拟合带边 D_it 上升 |
| G2 | 类 Pb0(类受主) | E_i + 0.3 eV | 拟合禁带中央 D_it |
| G3 | 类 Pb0(类施主) | E_i - 0.3 eV | 拟合禁带中央 D_it |
加上两个带尾态 G1 和 G4,总共 3-4 个高斯。
拟合优度:论文实际测量 vs 拟合曲线在 ±0.5 eV 范围内吻合。
10.5 D_it 对 V_oc 的影响
经验关系(本论文体系):
$$V_{oc} \approx 700 \text{ mV} - 60 \text{ mV} \times \log_{10}\left(\frac{D_{it}}{10^{10}}\right)$$
| D_it(cm⁻²·eV⁻¹) | 预期 V_oc | 实际(论文) |
|---|---|---|
| 10¹⁰ | 700 mV | — |
| 10¹¹ | 640 mV | — |
| 10¹² | 580 mV | ~ 600 mV |
| 10¹³ | 520 mV | — |
关键观察:D_it 每升高 10 倍,V_oc 降低 ~ 60 mV——这是为什么 SiO₂(D_it = 10¹⁰)的 V_oc 高达 720 mV,而未钝化硅(D_it = 10¹³)V_oc 不到 500 mV。
11 · 从 D_it 到 S_eff:钝化设计的核心公式
11.1 S_eff 的完整公式
本论文 4.6 节 + 5 章使用的核心公式(综合陷阱辅助复合 + 场效应):
$$S_{eff} = \frac{1}{\frac{1}{S_{it}} + \frac{1}{S_{field}} + \frac{1}{S_{auger}}}$$
其中:
- S_it:界面态相关复合速度
- S_field:场效应相关复合速度
- S_auger:表面俄歇复合速度
11.2 三个分量的具体表达
S_it(界面态主导)
$$S_{it} = \frac{\pi kT D_{it}}{q^2 N_A L_D} \cdot f(Q_f)$$
其中 f(Q_f) 是场效应修正函数(Q_f 越大,f 越小,S_it 越小)。
S_field(场效应主导)
$$S_{field} = \frac{q \mu_n E_s}{2}$$
其中 E_s 是表面电场,由 Q_f 决定。
S_auger(表面俄歇)
$$S_{auger} = \gamma \cdot n_s^2$$
11.3 钝化设计的核心思想
最小化 S_eff = 同时降低 D_it + 增大 |Q_f| + 控制表面载流子浓度。
| 优化目标 | 物理手段 | 论文做法 |
|---|---|---|
| 降低 D_it | 化学钝化(饱和悬挂键) | HF 漂洗 + Al₂O₃ 沉积 |
| 增大 |Q_f| | 场钝化(优化 Al₂O₃ 配位) | 退火 450°C 激活负电荷 |
| 控制 n_s | 场钝化(避免过度反型) | Q_f < 1.5×10¹³ cm⁻²(避免表面俄歇) |
11.4 论文 Al₂O₃ 的实测参数
| 试样 | D_it | Q_f | S_eff(实测) | S_eff(理论) |
|---|---|---|---|---|
| PE-ALD | 7.7×10¹¹ | -1×10¹³ | 13.8 cm/s | 14.5 cm/s ✓ |
| T-ALD | 4.0×10¹¹ | -5.3×10¹² | 23.5 cm/s | 22.8 cm/s ✓ |
理论值和实测值高度吻合——证明论文的扩展陷阱辅助模型(公式 3.27-3.31)是准确的。
11.5 S_eff 与 V_oc 的最终关系
本论文体系下(衬底电阻率 3 Ω·cm, 厚度 500 μm):
| S_eff(cm/s) | τ_eff(μs) | V_oc(mV) | 转换效率(%) |
|---|---|---|---|
| 1000 | 25 | 580 | 18% |
| 100 | 250 | 640 | 22% |
| 20 | 1200 | 680 | 25% |
| 10 | 2500 | 700 | 26% |
| 1 | 25000 | 720 | 26.5% |
每降低 1 个数量级的 S_eff,V_oc 提升 ~ 20-30 mV——钝化的边际收益递减但仍然显著。
本篇完