TOPCon 硼扩散温度与方阻的定量关系

硼扩散是 TOPCon 电池形成 PN 结的核心工艺，扩散温度通过 Arrhenius 关系决定硼原子的扩散系数，进而影响掺杂浓度分布、结深和方块电阻。本文面向光伏扩散工艺工程师和工艺整合工程师，从菲克扩散定律出发，推导方阻与扩散温度、时间的定量关系，给出硼在硅中的扩散激活能等关键参数，并讨论产线中方阻均匀性控制的工程约束，假定读者已了解石英管式扩散炉的基本操作和掺杂工艺概念。

1 理论基础：Arrhenius 扩散关系

1.1 扩散系数与温度

扩散系数 D 描述杂质原子在晶格中的迁移速率。硼在硅中以替位式扩散为主，通过占据空位进行迁移。D 与温度 T 的关系遵循 Arrhenius 公式：

D = D₀ · exp(−Ea / kT)

其中：

实验研究表明，硼在硅中的本征扩散激活能为 3.12 ± 0.04 eV，磷为 2.74 ± 0.07 eV。尽管两者均以硅自间隙原子介导扩散，磷因与 Si 自间隙原子结合更强而具有更低的激活能。这也是硼扩散温度（900 ~ 1100°C）高于磷扩散温度（800 ~ 900°C）的根本原因。

激活能的温度依赖性

Arrhenius 关系中 Ea 并非严格常数——硼的扩散机制在 900 ~ 1100°C 范围内以硅自间隙介导的"踢出"（kick-out）机制为主，Ea 约为 3.12 eV。但在低温段（< 850°C），部分扩散通过空位辅助机制进行，有效激活能偏高（~3.50 eV），D 的温度敏感性更强。这一"激活能分段"特性在快速热退火（RTA）工艺中需特别考虑，因为 RTA 的升降温速率（50 ~ 150°C/s）远快于常规管式炉（5 ~ 10°C/min），低温段（< 850°C）的累计扩散贡献通常可忽略，等效扩散时间可仅按 > 850°C 的高温段积分计算。

1.2 温度敏感性

Arrhenius 关系的指数特性决定了扩散系数对温度极度敏感。温度每升高 10°C，扩散系数 D 增加约 35% ~ 50%。这一特性具有双重意义：

• 工艺控制：温度漂移必须严格控制在 ±5°C 以内。某产线因热电偶老化导致实际温度偏高 22°C，结深从设计值 0.38 μm 漂移至 0.65 μm，方阻超标 2.4 倍
• 工艺设计：可利用温度梯度精确调控掺杂深度——不同温度档位的扩散行为差异巨大

温区均匀性的工程来源

管式扩散炉的温区分为恒温区（center zone）和过渡区（ramp zone），恒温区长度通常为 400 ~ 800 mm。即使恒温区设定温度一致，实际温度分布在轴向上可存在 ±3 ~ 5°C 的波动，源自加热丝布局密度、石英管的径向热辐射损失和气体对流冷却。径向温差（中心 vs 边缘）可达 2 ~ 4°C——管壁附近的硅片因与石英管壁的辐射换热而温度偏低。对于 D 增加 35%/10°C 的温度敏感性，径向 2°C 的温差对应 D 约 6% ~ 7% 的不均匀，直接反映为同一批次内不同位置硅片的方向阻差异。

2 方阻与扩散参数的关系

2.1 方阻的定义

方块电阻（Sheet Resistance, Rs 或 R□）定义为正方形薄膜两端之间的电阻：

Rs = ρ / xj = 1 / (q · μ · N_eff · xj)

其中 ρ 为电阻率，xj 为结深，q 为电子电荷，μ 为迁移率，N_eff 为有效掺杂浓度。

方阻是掺杂浓度和结深的综合表征——方阻越低，掺杂越重或结越深；方阻越高，掺杂越轻或结越浅。

四探针法测量原理

产线中方阻的标准测量工具为四探针测试仪（Four-Point Probe, 4PP）。四根等间距（s ≈ 1 mm）的钨探针以直线排列压在硅片表面，外侧两针通以恒定电流 I，内侧两针测量电压 V。对于无限大半无限厚样品，Rs = (π/ln2) · V/I = 4.532 · V/I。对于有限厚度样品（扩散层深度与探针间距可比拟），需引入几何修正因子 F(w/s, d/s)，其中 w 为样品边长，d 为扩散层深度。实际产线中使用校准片（已知 Rs 的标准样品）定期校验探头接触电阻和修正因子，测量精度可达 ±1%。

四探针测量的是扩散层内掺杂分布的加权平均电导——电流路径并非仅局限于表面，而是渗透至结深 xj 的整个扩散层，因此 Rs 反映的是从表面到 xj 整个剖面掺杂浓度的积分效应：

1/Rs = q · ∫₀^xj μ(N(x)) · N(x) dx

其中 N(x) 为掺杂浓度剖面。迁移率 μ 本身是掺杂浓度的函数（高掺杂时电离杂质散射增强，μ 下降），因此 Rs 与总掺杂剂量 Q_dose = ∫N(x)dx 并非严格线性关系，高剂量区域 μ 的下降导致 Rs 的降低幅度小于剂量的增加幅度。

2.2 结深的扩散模型

在有限源扩散（推进阶段）条件下，掺杂浓度分布为高斯函数：

C(x, t) = (Q / √(πDt)) · exp(−x² / 4Dt)

结深 xj 定义为掺杂浓度等于衬底本底浓度（NBG ≈ 5×10¹⁷ cm⁻³）处的深度：

xj ≈ 2 · √(D · t)

其中 t 为扩散时间。可见结深由扩散系数 D 和时间 t 共同决定。

高斯分布与误差函数分布的区分

扩散工艺中掺杂浓度剖面取决于边界条件：

• 有限源（limited source）——预沉积后在 BSG 层中固定剂量的硼原子进入推进阶段，无额外硼源补充，剖面为高斯分布：C(x) ∝ exp(−x²/4Dt)。峰值浓度随时间下降（C(0) = Q_dose/√(πDt)），总剂量守恒
• 恒定表面浓度（constant surface concentration）——扩散过程中表面硼浓度始终保持为固溶度极限 Csol，剖面为余误差函数：C(x) = Csol · erfc(x/2√(Dt))。积分剂量随时间增加，峰值浓度不变

TOPCon 硼扩散中推进阶段关闭硼源，属于有限源条件，浓度剖面以高斯分布为主。但若推进阶段 BSG 层中残留的 B₂O₃ 持续向硅体内供硼（BSG 作为准无限源），则剖面向余误差函数过渡，实际剖面介于两者之间。这一过渡行为在高温（> 1050°C）长推进时间条件下更为显著。

2.3 方阻与温度的关系

将结深表达式代入方阻定义，得到：

Rs ∝ 1 / √(D · t) ∝ 1 / √(exp(−Ea/kT) · t)

完整形式：

Rs = A · exp(Ea / 2kT)

即方阻与温度呈指数反比关系。温度升高 → D 指数增大 → 结深增加 → 方阻指数下降。

对于扩散时间的影响：Rs ∝ 1/√t，时间加倍，方阻降至原来的 1/√2 ≈ 71%。

3 定量仿真

3.1 典型工艺窗口

数据表明：温度从 950°C 升至 1100°C，结深增加约 3.8 倍，方阻降至约 30%。

仿真数据的工程验证

以上仿真基于高斯有限源模型和均匀掺杂假设。实际工艺中需要四探针实测验证，并根据实测结果反向校准扩散系数和激活能参数。典型校准流程：（1）选取 3 ~ 5 个不同温度/时间条件运行扩散实验；（2）四探针测量方阻；（3）SRP（Spreading Resistance Profiling，扩展电阻分析）或 ECV（Electrochemical Capacitance-Voltage，电化学 C-V）测量掺杂浓度剖面，提取实际结深和表面浓度；（4）将 Arrhenius 公式与实测数据拟合，获得该炉管在具体工艺条件下的有效激活能 Ea_eff 和有效指前因子 D₀_eff。

同一型号扩散炉的不同炉管之间因热电偶精度、石英管老化和气流分布的差异，有效扩散参数可偏离标称值 5% ~ 15%，校准是工艺转移和量产扩线的必要步骤。

3.2 TOPCon 硼扩的工程约束

TOPCon 电池追求高方阻 + 浅结以降低表面俄歇复合和死层效应。SE（Selective Emitter）工艺通过在栅线下方局部重掺杂解决高方阻带来的接触电阻问题，实现"高方阻（高 Voc）+ 低接触电阻（高 FF）"的最优组合。

氧化增强扩散效应（OED）

硼扩散中一个不可忽视的次级效应是氧化增强扩散（Oxidation-Enhanced Diffusion, OED）。硅表面在 O₂ 氛围中热氧化时，Si/SiO₂ 界面处注入大量硅自间隙原子（1 cm³ SiO₂ 生长约产生 8×10²² 个间隙原子）。这些过剩间隙原子增强替位硼的扩散速率——D_eff = D_intrinsic + ΔD_OED，ΔD_OED ∝ (dX_ox/dt)^n（氧化速率），n ≈ 0.3 ~ 0.5。

推进阶段若同时通入 O₂（原位氧化），OED 效应将导致结深略大于纯 N₂ 推进的预测值——典型偏差为 5% ~ 15%，取决于氧化温度和 O₂ 分压。工艺设计和仿真中若忽略 OED 效应，方阻的预测值将偏高（即实际扩散深度大于预期），导致结深超标。补救策略：（1）推进阶段使用 N₂ 为主、极少量 O₂（< 1%）的气氛——O₂ 仅用于维持 SiO₂/Si 界面质量，氧化速率和 OED 效应最小化；（2）在工艺模型中显式引入 OED 修正项。

4 多步扩散的方阻控制

传统单步扩散的温度-方阻关系存在固有矛盾：温度过高 → 方阻低、结深过大、死层严重；温度过低 → 方阻高、掺杂不足、发射极串联电阻大。

多步扩散通过分解沉积和推进阶段解耦表面浓度与结深的控制：

多步扩散 + 梯度降温氧化可有效降低表面非激活硼含量，减少死层厚度，在相同方阻条件下获得更浅的结深和更高的 Voc（+3 ~ 5 mV）。

4.1 两次沉积的剂量分配策略

第一次沉积（低流量）用于建立初始掺杂层，剂量约占总剂量的 30% ~ 50%。该层中的硼原子在推进阶段作为"前沿扩散源"率先向深处迁移，主导结深的形成。第二次沉积（略高流量）补充剩余剂量，硼原子主要分布在与第一次沉积层的过渡区内。

分配策略的优化目标：总剂量相同的前提下，两次沉积的剂量比和温度窗口共同决定最终掺杂剖面的形状——第一次剂量占比越高，剖面越"平顶"（表面至一定深度内浓度均匀），发射极串联电阻越低但死层越厚；第二次剂量占比越高，剖面越"尖峰"（浓度从表面向体内单调快速下降），死层较薄但发射极电阻较高。量产中通过 DOE（Design of Experiments）确定最优剂量比，平衡 Voc 和 FF 的取舍。

5 温度变化的等效时间换算

Arrhenius 公式的定积分形式可用于不同温度间的等效扩散时间换算：

t₂ / t₁ = exp[(Ea/k) · (1/T₂ − 1/T₁)]

例如：950°C（1223 K）扩散 30 min 等效于 1050°C（1323 K）扩散多长时间？

t₂ = 30 × exp[(3.12/8.617×10⁻⁵) · (1/1323 − 1/1223)] ≈ 30 × exp(−2.19) ≈ 3.4 min

即 1050°C 扩散 3.4 min 即可达到 950°C 扩散 30 min 的相同扩散深度。这一换算可用于工艺转移和热预算优化。

5.1 多段温度曲线的等效积分

实际扩散工艺包含升温和降温的过渡阶段，并非全程恒温。完整的热预算需对时间-温度曲线进行积分：

(Dt)_effective = ∫ D(T(t)) dt = D₀ · ∫ exp(−Ea/kT(t)) dt

对于控温精度 ±5°C 的现代扩散炉，升降温阶段的积分贡献通常 < 5% 总热预算，简化为仅计算恒温段时间的 D×t 积的误差可接受。但对于 RTA（升降温速率 > 50°C/s），恒温段时间仅 1 ~ 30 s，升降温阶段的积分贡献可达总热预算的 20% ~ 40%，使用等效恒温时间模型会产生显著偏差，须使用完整积分公式计算。

6 方阻均匀性的产线控制

6.1 片内均匀性

片内方阻均匀性（within-wafer uniformity）定义为：（Rs_max − Rs_min）/（Rs_max + Rs_min）× 100%，或以标准差/均值（1σ）表示。典型 25 点或 49 点方阻 mapping 可揭示系统性非均匀特征：

• 对称型（中心高/低、边缘高/低）：来源于炉管径向温度梯度和气流分布，通过调整炉管加热区功率配比和气体流量均流器（baffle）优化
• 非对称型（一侧高/低）：来源于硅片在石英舟中的倾斜、贴壁或气流遮挡，通过调整石英舟间距和装载方式纠正

6.2 片间和批间均匀性

片间均匀性反映同一批次内不同硅片位置（炉口 vs 炉中 vs 炉尾）的方阻差异，通常源于轴向温度梯度和硼源浓度沿气流方向的消耗递减。批间均匀性反映不同扩散批次之间的方阻漂移，主要来源为热电偶老化（导致温度测量偏离真实值）、石英管壁 BSG 累积（改变管内热辐射特性）和硼源纯度波动。

产线统计过程控制（SPC）方案：每批次抽取 3 ~ 5 片（炉前、炉中、炉后各 1 ~ 2 片），每片测量 5 ~ 9 点方阻，记录均值和控制图。控制上限（UCL）/下限（LCL）通常设为目标值 ±（15 ~ 20）Ω/□。连续 3 点超出 ±1σ 线或单点超出 ±3σ 线触发停线排查。

7 总结

TOPCon 硼扩散的方阻与温度之间遵循Rs ∝ exp(Ea/2kT) 的 Arrhenius 指数关系。硼在硅中的扩散激活能（3.12 ~ 3.50 eV）决定了扩散系数对温度的极度敏感性——10°C 温度漂移导致 35% ~ 50% 的扩散系数变化。工艺实践中，方阻的控制通过温度（粗调）+ 时间（微调）+ 多步沉积预淀积（剂量调控）的组合实现。氧化增强扩散（OED）效应可导致方阻偏离预测值 5% ~ 15%，工艺模型中需显式修正。SE 激光重掺技术解耦了方阻与 Voc 的固有矛盾，使高方阻（> 130 Ω/□）浅结发射极成为 TOPCon 高效工艺的标准选择。产线方阻均匀性的持续监控依赖四探针 mapping 和 SPC 控制图，温度校准和气流优化是降低批间漂移的两个主攻方向。