Administrator
发布于 2026-07-31 / 0 阅读
0

陷阱辅助复合理论与扩展界面模型

04 · 陷阱辅助复合理论与扩展界面模型

系列文章:太阳电池晶体硅界面钝化研究
本文定位:从"陷阱能级"这个最基础概念开始,用统计力学严格推导陷阱辅助复合公式,然后推导论文第 3 章的扩展 SRH 模型
前置知识:建议先读 01 · 半导体物理基础02 · 表面与界面化学
写作原则:每个公式都从最基本的物理概念开始推导,不跳步;每个符号都明确物理意义


目录

  1. 从单个陷阱开始:陷阱能级的物理图像
  2. 陷阱的占据率:费米-狄拉克统计
  3. 陷阱辅助复合的 4 个基本过程
  4. 稳态条件与 R_SRH 公式推导
  5. 本征费米能级 E_i 与特征浓度 n₁, p₁
  6. 电子/空穴寿命 τ_n, τ_p 的物理含义
  7. 从体复合到表面复合:边界条件
  8. 论文扩展模型:5 个新公式(3.27-3.31)
  9. 反型层条件:Q_f 临界值计算
  10. D_it 谱的物理意义:为什么禁带中央最致命
  11. 从 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

复合速率 = 单位时间、单位体积内,一个电子-空穴对通过陷阱消失

每个复合事件需要:

  1. 过程 ①(电子被捕获)或过程 ③(空穴被中和)
  2. 两个过程都发生,才完成一次 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 表面复合的特殊性

表面 = 半无限半导体的边界

在表面,有几个变化:

  1. 载流子分布不均:表面有电场(场效应)时,载流子在表面聚集
  2. 陷阱密度极高:表面态 D_it 通常比体内 N_T 高很多
  3. 没有"另一侧":表面只能向体内复合(单向)

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——钝化的边际收益递减但仍然显著。



本篇完