简介围绕半导体掺杂原理整理的入门级知识文档面向微电子、集成电路及材料相关专业学生与初入行的技术人员。内容从本征半导体导电能力低讲起说明杂质对载流子数目的影响再依次展开掺磷形成N型硅、掺硼形成P型硅的价电子变化讲解donor、acceptor、多数与少数载流子等术语并延伸到counterdoping、III-V、II-VI化合物半导体及砷化镓器件应用帮助读者建立掺杂工艺与器件原理的整体认识。资源包内含1个PDF文件约68KB为图文结合的章节式讲义附掺磷、掺硼硅简化晶体结构示意图与半导体术语汇总表篇幅紧凑适合课堂补充或复习查阅。已有2750人学习下载可快速梳理概念脉络为后续集成电路工艺学习打基础。1. 掺杂到底改了什么从本征半导体到可控导电室温下高纯硅的本征载流子浓度只有约 1×10¹⁰ cm⁻³对应本征电阻率在 10⁵ Ω·cm 量级比常见金属差了十几个数量级这种材料直接做器件几乎没法用。掺杂做的事情很集中把特定杂质原子按可控浓度放进晶格替位让载流子的类型和浓度都能被设计。掺磷给出电子得到 n 型掺硼给出空穴得到 p 型浓度从 10¹⁴ 到 10²⁰ cm⁻³ 跨六个数量级电阻率随之能压到毫欧量级。做器件、跑工艺、做 TCAD 仿真和材料建模的人都会反复回到这条链路掺杂浓度决定费米能级位置费米能级位置决定载流子浓度迁移率再决定电阻率。抓住这条链路掺杂就不再是一堆数字而是一组可反推、可验证的工程参数。2. 施主与受主掺杂的能带模型与载流子浓度计算2.1 替位掺杂如何改变晶格与化学键硅是金刚石立方结构晶格常数约 0.543 nm每个硅原子与四个最近邻形成 sp³ 共价键四个价电子全部参与成键晶格在绝对零度附近表现为绝缘体。把一个 V 族原子磷、砷、锑放到硅原子位置上它用四个价电子参与共价键第五个价电子只受到原子核和周围晶格的弱库仑束缚很小的能量就能把它激发到导带成为自由电子。这类杂质叫施主束缚能叫施主电离能。III 族原子硼、铝、镓、只有三个价电子替位后共价键缺一个电子相当于产生一个可以移动的空穴这类杂质叫受主。受主电离后空穴进入价带参与导电。两类杂质改变的是同一个东西本征材料里电子和空穴的平衡被打破多子浓度由杂质浓度主导少子浓度被质量作用定律 n·p n_i² 约束。这里有个容易忽略的前提掺杂原子必须占据替位位置。间隙位置的杂质、团簇或者沉淀相通常不贡献自由载流子所以后面要区分化学浓度和电活性浓度两者在重掺杂下能差出近一个数量级。2.2 电离能与杂质能级浅能级和深能级的分界杂质在禁带里引入的能级位置决定了它在室温下能不能有效贡献载流子。判断尺子是室温热能 kT ≈ 0.026 eV如果电离能远小于这个值杂质在室温下基本全部电离如果接近禁带中央它更可能充当复合中心而不是掺杂剂。杂质类型电离能 (eV)室温电离情况在硅中的角色磷 P施主0.045基本全电离常用 n 型掺杂砷 As施主0.054基本全电离离子注入常用锑 Sb施主0.039基本全电离特殊 n 型硼 B受主0.045基本全电离常用 p 型掺杂铝 Al受主0.057基本全电离少见受主金 Au深能级0.54基本不电离少子寿命控制铜 Cu深能级0.24基本不电离污染源、复合中心浅能级杂质和深能级杂质的选择性很明显想控制载流子浓度就选浅能级想控制载流子寿命就故意引入深能级。固溶度是另一条硬边界硼在硅里的固溶度上限大约在 5×10²⁰ cm⁻³ 量级超过这个值多余杂质会以团簇或沉淀形式存在不再贡献载流子这也是重掺杂源漏区浓度通常卡在这个量级附近的原因。2.3 用 Python 从掺杂浓度算费米能级与载流子浓度掺杂浓度和费米能级的关系可以用有效态密度直接算出来。n 型完全电离近似下 n ≈ N_D再用 n N_C·exp(-(E_C-E_F)/kT) 反解 E_F和本征费米能级 E_i 的差值就是掺杂把费米能级推离禁带中央的距离。import math # 物理常数与硅材料参数 k_B 8.617333e-5 # 玻尔兹曼常数, eV/K T 300.0 # 温度, K N_C 2.8e19 # 导带有效态密度, cm^-3 N_V 1.1e19 # 价带有效态密度, cm^-3 E_g 1.12 # 硅室温禁带宽度, eV n_i 1.0e10 # 本征载流子浓度, cm^-3 kT k_B * T def fermi_shift_n(N_D): n 型半导体费米能级相对本征能级的偏移 (eV)完全电离近似 return kT * math.log(N_D / n_i) def fermi_shift_p(N_A): p 型半导体费米能级相对本征能级的偏移 (eV) return -kT * math.log(N_A / n_i) # 扫描几个典型掺杂浓度 for N in [1e14, 1e16, 1e18, 1e20]: dE fermi_shift_n(N) n N p n_i**2 / n # 质量作用定律给出少子浓度 print(fN_D{N:.0e} cm^-3 n{n:.2e} p{p:.2e} E_F-E_i{dE:.4f} eV)这段代码做三件事用对数关系把掺杂浓度映射成费米能级偏移用质量作用定律从多子浓度推少子浓度再打印出四个数量级上的变化趋势。参数上要注意两点N_C 和 N_V 取的是硅在 300 K 的有效态密度换材料Ge、GaAs、SiC要重新取值n_i 对温度非常敏感300 K 大约 1×10¹⁰ cm⁻³到 400 K 会涨到 10¹² 量级高温下本征载流子开始压倒掺杂浓度器件会失稳。跑出来的结果是 N_D1×10¹⁴ 时 E_F-E_i 约 0.24 eV1×10¹⁶ 时约 0.35 eV到 1×10²⁰ 时已经接近导带边 0.55 eV。费米能级一旦逼近导带边玻尔兹曼近似开始失效必须换成费米-狄拉克积分同时出现带隙变窄效应这是重掺杂模型和轻掺杂模型的分界点。3. 半导体掺杂的主流工艺路线扩散、离子注入与外延原位掺杂3.1 高温扩散固溶度限制下的推进与再分布扩散是较早出现的掺杂方法把硅片放在含杂质源的高温气氛里杂质从表面进入晶格并向体内推进。常见做法分两步先做恒定表面源扩散再关掉源做有限源推进。结深和表面浓度由温度和时间的组合决定常见的经验关系是结深 x_j ≈ 2√(D·t)其中扩散系数 D 服从 Arrhenius 形式。import math k_B 8.617e-5 # eV/K D0 0.037 # 硼在硅中的扩散前置因子, cm^2/s Ea 3.46 # 硼扩散激活能, eV T 1273.0 # 扩散温度, K (约 1000 摄氏度) t 1800.0 # 扩散时间, s (30 分钟) D D0 * math.exp(-Ea / (k_B * T)) # 扩散系数, cm^2/s x_j 2.0 * math.sqrt(D * t) # 特征结深估计, cm print(fD {D:.3e} cm^2/s) print(fx_j {x_j*1e4:.3f} um)参数含义D0 和 Ea 是材料对硼在硅中的这组值决定扩散快慢T 和 t 是工艺窗口的可调量温度每升高约 50 K扩散系数大致翻倍所以扩散炉的温度均匀性直接决定片内结深一致性。逻辑上这是平方根关系想加倍结深就要四倍扩散时间实际产线更倾向用温度微调而不是拉长时间避免热预算过大影响已形成的浅结。扩散的边界同样明显受固溶度限制表面浓度很难突破 5×10²⁰ cm⁻³横向扩散和纵向推进同量级做深亚微米浅结几乎无能为力需要高温长时间热预算大。这也是为什么在微米级及更小节点之后扩散逐步被离子注入取代只在一些特殊结构里保留。3.2 离子注入剂量、能量、射程与退火激活离子注入把杂质离子加速到几十到几百 keV直接轰进硅片注入剂量决定浓度注入能量决定射程和峰值深度两者解耦控制这是它取代扩散的核心优势。注入后的杂质落在间隙位置必须经过退火才能进入替位、恢复晶格损伤并实现电激活。参数典型范围控制的目标注入能量10 keV ~ 500 keV射程与峰值深度注入剂量10¹² ~ 10¹⁶ cm⁻²积分浓度倾斜角0° ~ 7°抑制沟道效应退火温度900 ~ 1100 ℃激活率与损伤恢复退火时间毫秒级尖峰到数十秒控制扩散量两个必须注意的点。倾斜角如果不加离子会沿晶格沟道穿得很深形成拖尾分布器件一致性会变差常见做法是偏 7°。退火窗口要同时兼顾激活和扩散尖峰退火或毫秒退火能在高激活的同时把杂质扩散压到最小这是浅结掺杂的关键工艺。3.3 原位掺杂与外延界面陡峭度的取舍外延生长时同步通入掺杂气体杂质在薄膜生长过程中直接进入晶格这叫原位掺杂。它和注入退火的区别在于掺杂原子在生长阶段就处在替位位置不需要高剂量注入和高强度退火因此能把掺杂剖面控制得很陡界面过渡区可以做到纳米量级。代价是只能在能外延的层上做浓度上限受生长气氛和固溶度约束厚层外延的热预算也不低。选择时可以按三条线判断要浅结和陡界面优先原位掺杂外延要局部选区掺杂、对片内均匀性要求高优先离子注入要低成本深结、对结深精度不敏感才考虑扩散。4. 掺杂结构建模与电学参数标定4.1 Materials Studio 里建掺杂超胞与 Make P1 的取舍用 Materials Studio 做掺杂结构优化绕不开两个动作建足够大的超胞避免掺杂原子和它的周期性镜像相互作用以及在放好掺杂原子后处理对称性。原始硅的晶胞是 Fd-3m 空间群默认对称操作会把一个掺杂位点复制成多个等价位置或者在你替换原子后强行把结构拉回高对称导致掺杂原子实际被放到了不等价位置。常见做法是替换原子之前或之后执行 Make P1把对称性降到 P1只保留平移对称这样掺杂位点和局部弛豫就能被正确描述。是否需要 Make P1 取决于模型目的只想看单个替位原子的局域结构Make P1 后放开全部原子坐标做几何优化更稳想做高对称路径或声子计算则要权衡对称性带来的计算量收益和掺杂引入的实际破缺。判断标准很简单替换后如果发现多个位置被同时改变或者能量异常低得不像掺杂体系就检查一下对称性有没有被保留。超胞尺寸和掺杂浓度的换算要提前算清楚不然建模浓度和实验浓度对不上。a 0.543e-7 # 硅晶格常数, cm atoms_per_cell 8 # 金刚石结构每晶胞 8 个原子 nx, ny, nz 3, 3, 3 # 超胞沿三个方向的晶胞数 n_cells nx * ny * nz n_atoms n_cells * atoms_per_cell V n_cells * a**3 # 超胞体积, cm^3 density n_atoms / V # 原子数密度, cm^-3 doping density / n_atoms # 替换一个原子对应的体浓度, cm^-3 print(f总原子数 {n_atoms}) print(f原子数密度 {density:.3e} cm^-3) print(f替换 1 个原子对应掺杂 {doping:.3e} cm^-3)逻辑说明原子数密度由晶格常数和每晶胞原子数直接得到替换一个原子相当于掺杂浓度等于原子数密度除以总原子数。3×3×3 超胞的 216 个原子里换一个对应约 2.3×10²⁰ cm⁻³正好落在重掺杂量级想模拟 10¹⁸ cm⁻³ 的轻掺杂就要把超胞扩到几十个晶胞边长第一性原理计算的成本会迅速上升这也是很多掺杂建模只做趋势对照而不直接对齐产线浓度的原因。4.2 掺杂浓度、电阻率与方块电阻的换算电学参数之间有一条可直接用的关系链。电阻率 ρ 1/(q·n·μ_n)方块电阻 R_sh ρ/x_j其中 x_j 是结深。迁移率不是常数它随掺杂浓度上升而下降硅中电子迁移率常用 Caughey-Thomas 形式描述。import math q 1.602e-19 # 元电荷, C mu_max 1400.0 # 低掺杂电子迁移率, cm^2/(V·s) mu_min 65.0 # 高掺杂渐近电子迁移率, cm^2/(V·s) N_ref 8.5e16 # 参考浓度, cm^-3 alpha 0.72 # 拟合指数 def mu_n(N): 硅中电子迁移率随掺杂浓度的经验模型 return mu_min (mu_max - mu_min) / (1 (N / N_ref)**alpha) def sheet_resistance(N_D, x_j_um): 由掺杂浓度和结深估算方块电阻, 单位 Ω/sq n N_D mu mu_n(n) x_j x_j_um * 1e-4 # um - cm rho 1.0 / (q * n * mu) # Ω·cm return rho / x_j for N in [1e16, 1e18, 1e19, 1e20]: print(fN_D{N:.0e} cm^-3 mu{mu_n(N):.1f} cm^2/Vs R_sh{sheet_resistance(N, 0.2):.1f} Ω/sq)参数说明N_ref 和 alpha 决定迁移率随浓度衰减的快慢不同材料、不同载流子类型要换不同的拟合参数x_j 单位换算成厘米后才能和电阻率一起给出方块电阻。这段计算的意义是把掺杂浓度、迁移率退化和工艺可测的方块电阻连起来实测方块电阻偏高时可能来自激活率不足、结深偏浅或者补偿度过大而不是掺杂浓度本身写错了。4.3 Hall、CV、SIMS 各测什么别把化学浓度当电活性浓度掺杂表征方法各管一段信息混用会得出错误结论。方法直接测量量给出的是化学浓度还是电活性浓度主要限制霍尔 Hall载流子浓度、类型、迁移率电活性浓度面密度需要良好欧姆接触受表面态影响CV载流子浓度-深度分布电活性浓度可分辨深度需要肖特基结或 PN 结靠耗尽近似SIMS元素浓度-深度分布化学浓度含未激活破坏性标样依赖强扩展电阻 SRP载流子浓度-深度分布电活性浓度探针接触、结深附近分辨率下降四探针方块电阻电学综合量需配合结深才能反推浓度一条实用原则SIMS 剖面给出的是总杂质浓度霍尔和 CV 给出的是电活性浓度两者相除就是激活率。如果 SIMS 显示杂质很多但霍尔测到的载流子明显偏少问题通常出在退火不充分、杂质落在间隙位置或者形成了团簇而不是扩散或注入的剂量算错了。5. 用实测方块电阻反推掺杂浓度与补偿比工程上更常见的场景是只有电学测试数据要从方块电阻和结深反推电活性浓度。区块电阻和结深、迁移率的关系可以直接反解得到的是有效电活性面浓度。q 1.602e-19 R_sh 120.0 # 实测方块电阻, Ω/sq x_j 0.35 # 结深, um mu 200.0 # 由霍尔测试得到的迁移率, cm^2/V·s x_j_cm x_j * 1e-4 N_eff 1.0 / (q * mu * R_sh * x_j_cm) # 有效电活性体浓度, cm^-3 print(f有效电活性浓度 ≈ {N_eff:.3e} cm^-3)这里的 N_eff 是平均到结深上的有效浓度实际情况并非均匀分布高浓度集中在表面附近时这个平均值会低于峰值浓度。更精细的做法是把 SRP 或 CV 得到的浓度剖面分段积分每段用对应浓度下的迁移率再求和得到面浓度和方块电阻能较好地还原剖面形状。补偿比是排错时最容易被忽略的一项。n 型区的电活性浓度不是 N_D而是有效值 N_D − N_A。当注入 n 型杂质的区域本身有较高背景 p 型掺杂或者工艺引入了反型杂质实测载流子浓度会明显低于预期但 SIMS 上的施主总量却是对的这种偏差会随着补偿比接近 1 而急剧放大。遇到激活率看起来很低的情况先算一遍背景掺杂浓度再判断是不是补偿问题。温度相关参数也要一起看。迁移率随温度下降而上升本征载流子浓度随温度上升而快速上升两条曲线方向相反导致高温下用固定迁移率反推的浓度会偏大。可靠做法是在器件工作温度下重测霍尔用实测迁移率而不是查表值去反推浓度这一步能把中高温下的浓度偏差压下来。最后一点落在剖面验证上反推得到的浓度用 SIMS 剖面交叉核对确认峰值位置和积分剂量两块数据对得上掺杂工艺窗口才算真正锁定。本文还有配套的精品资源点击获取 半导体掺杂浓度与费米能级：载流子计算、工艺与方块电阻
简介围绕半导体掺杂原理整理的入门级知识文档面向微电子、集成电路及材料相关专业学生与初入行的技术人员。内容从本征半导体导电能力低讲起说明杂质对载流子数目的影响再依次展开掺磷形成N型硅、掺硼形成P型硅的价电子变化讲解donor、acceptor、多数与少数载流子等术语并延伸到counterdoping、III-V、II-VI化合物半导体及砷化镓器件应用帮助读者建立掺杂工艺与器件原理的整体认识。资源包内含1个PDF文件约68KB为图文结合的章节式讲义附掺磷、掺硼硅简化晶体结构示意图与半导体术语汇总表篇幅紧凑适合课堂补充或复习查阅。已有2750人学习下载可快速梳理概念脉络为后续集成电路工艺学习打基础。1. 掺杂到底改了什么从本征半导体到可控导电室温下高纯硅的本征载流子浓度只有约 1×10¹⁰ cm⁻³对应本征电阻率在 10⁵ Ω·cm 量级比常见金属差了十几个数量级这种材料直接做器件几乎没法用。掺杂做的事情很集中把特定杂质原子按可控浓度放进晶格替位让载流子的类型和浓度都能被设计。掺磷给出电子得到 n 型掺硼给出空穴得到 p 型浓度从 10¹⁴ 到 10²⁰ cm⁻³ 跨六个数量级电阻率随之能压到毫欧量级。做器件、跑工艺、做 TCAD 仿真和材料建模的人都会反复回到这条链路掺杂浓度决定费米能级位置费米能级位置决定载流子浓度迁移率再决定电阻率。抓住这条链路掺杂就不再是一堆数字而是一组可反推、可验证的工程参数。2. 施主与受主掺杂的能带模型与载流子浓度计算2.1 替位掺杂如何改变晶格与化学键硅是金刚石立方结构晶格常数约 0.543 nm每个硅原子与四个最近邻形成 sp³ 共价键四个价电子全部参与成键晶格在绝对零度附近表现为绝缘体。把一个 V 族原子磷、砷、锑放到硅原子位置上它用四个价电子参与共价键第五个价电子只受到原子核和周围晶格的弱库仑束缚很小的能量就能把它激发到导带成为自由电子。这类杂质叫施主束缚能叫施主电离能。III 族原子硼、铝、镓、只有三个价电子替位后共价键缺一个电子相当于产生一个可以移动的空穴这类杂质叫受主。受主电离后空穴进入价带参与导电。两类杂质改变的是同一个东西本征材料里电子和空穴的平衡被打破多子浓度由杂质浓度主导少子浓度被质量作用定律 n·p n_i² 约束。这里有个容易忽略的前提掺杂原子必须占据替位位置。间隙位置的杂质、团簇或者沉淀相通常不贡献自由载流子所以后面要区分化学浓度和电活性浓度两者在重掺杂下能差出近一个数量级。2.2 电离能与杂质能级浅能级和深能级的分界杂质在禁带里引入的能级位置决定了它在室温下能不能有效贡献载流子。判断尺子是室温热能 kT ≈ 0.026 eV如果电离能远小于这个值杂质在室温下基本全部电离如果接近禁带中央它更可能充当复合中心而不是掺杂剂。杂质类型电离能 (eV)室温电离情况在硅中的角色磷 P施主0.045基本全电离常用 n 型掺杂砷 As施主0.054基本全电离离子注入常用锑 Sb施主0.039基本全电离特殊 n 型硼 B受主0.045基本全电离常用 p 型掺杂铝 Al受主0.057基本全电离少见受主金 Au深能级0.54基本不电离少子寿命控制铜 Cu深能级0.24基本不电离污染源、复合中心浅能级杂质和深能级杂质的选择性很明显想控制载流子浓度就选浅能级想控制载流子寿命就故意引入深能级。固溶度是另一条硬边界硼在硅里的固溶度上限大约在 5×10²⁰ cm⁻³ 量级超过这个值多余杂质会以团簇或沉淀形式存在不再贡献载流子这也是重掺杂源漏区浓度通常卡在这个量级附近的原因。2.3 用 Python 从掺杂浓度算费米能级与载流子浓度掺杂浓度和费米能级的关系可以用有效态密度直接算出来。n 型完全电离近似下 n ≈ N_D再用 n N_C·exp(-(E_C-E_F)/kT) 反解 E_F和本征费米能级 E_i 的差值就是掺杂把费米能级推离禁带中央的距离。import math # 物理常数与硅材料参数 k_B 8.617333e-5 # 玻尔兹曼常数, eV/K T 300.0 # 温度, K N_C 2.8e19 # 导带有效态密度, cm^-3 N_V 1.1e19 # 价带有效态密度, cm^-3 E_g 1.12 # 硅室温禁带宽度, eV n_i 1.0e10 # 本征载流子浓度, cm^-3 kT k_B * T def fermi_shift_n(N_D): n 型半导体费米能级相对本征能级的偏移 (eV)完全电离近似 return kT * math.log(N_D / n_i) def fermi_shift_p(N_A): p 型半导体费米能级相对本征能级的偏移 (eV) return -kT * math.log(N_A / n_i) # 扫描几个典型掺杂浓度 for N in [1e14, 1e16, 1e18, 1e20]: dE fermi_shift_n(N) n N p n_i**2 / n # 质量作用定律给出少子浓度 print(fN_D{N:.0e} cm^-3 n{n:.2e} p{p:.2e} E_F-E_i{dE:.4f} eV)这段代码做三件事用对数关系把掺杂浓度映射成费米能级偏移用质量作用定律从多子浓度推少子浓度再打印出四个数量级上的变化趋势。参数上要注意两点N_C 和 N_V 取的是硅在 300 K 的有效态密度换材料Ge、GaAs、SiC要重新取值n_i 对温度非常敏感300 K 大约 1×10¹⁰ cm⁻³到 400 K 会涨到 10¹² 量级高温下本征载流子开始压倒掺杂浓度器件会失稳。跑出来的结果是 N_D1×10¹⁴ 时 E_F-E_i 约 0.24 eV1×10¹⁶ 时约 0.35 eV到 1×10²⁰ 时已经接近导带边 0.55 eV。费米能级一旦逼近导带边玻尔兹曼近似开始失效必须换成费米-狄拉克积分同时出现带隙变窄效应这是重掺杂模型和轻掺杂模型的分界点。3. 半导体掺杂的主流工艺路线扩散、离子注入与外延原位掺杂3.1 高温扩散固溶度限制下的推进与再分布扩散是较早出现的掺杂方法把硅片放在含杂质源的高温气氛里杂质从表面进入晶格并向体内推进。常见做法分两步先做恒定表面源扩散再关掉源做有限源推进。结深和表面浓度由温度和时间的组合决定常见的经验关系是结深 x_j ≈ 2√(D·t)其中扩散系数 D 服从 Arrhenius 形式。import math k_B 8.617e-5 # eV/K D0 0.037 # 硼在硅中的扩散前置因子, cm^2/s Ea 3.46 # 硼扩散激活能, eV T 1273.0 # 扩散温度, K (约 1000 摄氏度) t 1800.0 # 扩散时间, s (30 分钟) D D0 * math.exp(-Ea / (k_B * T)) # 扩散系数, cm^2/s x_j 2.0 * math.sqrt(D * t) # 特征结深估计, cm print(fD {D:.3e} cm^2/s) print(fx_j {x_j*1e4:.3f} um)参数含义D0 和 Ea 是材料对硼在硅中的这组值决定扩散快慢T 和 t 是工艺窗口的可调量温度每升高约 50 K扩散系数大致翻倍所以扩散炉的温度均匀性直接决定片内结深一致性。逻辑上这是平方根关系想加倍结深就要四倍扩散时间实际产线更倾向用温度微调而不是拉长时间避免热预算过大影响已形成的浅结。扩散的边界同样明显受固溶度限制表面浓度很难突破 5×10²⁰ cm⁻³横向扩散和纵向推进同量级做深亚微米浅结几乎无能为力需要高温长时间热预算大。这也是为什么在微米级及更小节点之后扩散逐步被离子注入取代只在一些特殊结构里保留。3.2 离子注入剂量、能量、射程与退火激活离子注入把杂质离子加速到几十到几百 keV直接轰进硅片注入剂量决定浓度注入能量决定射程和峰值深度两者解耦控制这是它取代扩散的核心优势。注入后的杂质落在间隙位置必须经过退火才能进入替位、恢复晶格损伤并实现电激活。参数典型范围控制的目标注入能量10 keV ~ 500 keV射程与峰值深度注入剂量10¹² ~ 10¹⁶ cm⁻²积分浓度倾斜角0° ~ 7°抑制沟道效应退火温度900 ~ 1100 ℃激活率与损伤恢复退火时间毫秒级尖峰到数十秒控制扩散量两个必须注意的点。倾斜角如果不加离子会沿晶格沟道穿得很深形成拖尾分布器件一致性会变差常见做法是偏 7°。退火窗口要同时兼顾激活和扩散尖峰退火或毫秒退火能在高激活的同时把杂质扩散压到最小这是浅结掺杂的关键工艺。3.3 原位掺杂与外延界面陡峭度的取舍外延生长时同步通入掺杂气体杂质在薄膜生长过程中直接进入晶格这叫原位掺杂。它和注入退火的区别在于掺杂原子在生长阶段就处在替位位置不需要高剂量注入和高强度退火因此能把掺杂剖面控制得很陡界面过渡区可以做到纳米量级。代价是只能在能外延的层上做浓度上限受生长气氛和固溶度约束厚层外延的热预算也不低。选择时可以按三条线判断要浅结和陡界面优先原位掺杂外延要局部选区掺杂、对片内均匀性要求高优先离子注入要低成本深结、对结深精度不敏感才考虑扩散。4. 掺杂结构建模与电学参数标定4.1 Materials Studio 里建掺杂超胞与 Make P1 的取舍用 Materials Studio 做掺杂结构优化绕不开两个动作建足够大的超胞避免掺杂原子和它的周期性镜像相互作用以及在放好掺杂原子后处理对称性。原始硅的晶胞是 Fd-3m 空间群默认对称操作会把一个掺杂位点复制成多个等价位置或者在你替换原子后强行把结构拉回高对称导致掺杂原子实际被放到了不等价位置。常见做法是替换原子之前或之后执行 Make P1把对称性降到 P1只保留平移对称这样掺杂位点和局部弛豫就能被正确描述。是否需要 Make P1 取决于模型目的只想看单个替位原子的局域结构Make P1 后放开全部原子坐标做几何优化更稳想做高对称路径或声子计算则要权衡对称性带来的计算量收益和掺杂引入的实际破缺。判断标准很简单替换后如果发现多个位置被同时改变或者能量异常低得不像掺杂体系就检查一下对称性有没有被保留。超胞尺寸和掺杂浓度的换算要提前算清楚不然建模浓度和实验浓度对不上。a 0.543e-7 # 硅晶格常数, cm atoms_per_cell 8 # 金刚石结构每晶胞 8 个原子 nx, ny, nz 3, 3, 3 # 超胞沿三个方向的晶胞数 n_cells nx * ny * nz n_atoms n_cells * atoms_per_cell V n_cells * a**3 # 超胞体积, cm^3 density n_atoms / V # 原子数密度, cm^-3 doping density / n_atoms # 替换一个原子对应的体浓度, cm^-3 print(f总原子数 {n_atoms}) print(f原子数密度 {density:.3e} cm^-3) print(f替换 1 个原子对应掺杂 {doping:.3e} cm^-3)逻辑说明原子数密度由晶格常数和每晶胞原子数直接得到替换一个原子相当于掺杂浓度等于原子数密度除以总原子数。3×3×3 超胞的 216 个原子里换一个对应约 2.3×10²⁰ cm⁻³正好落在重掺杂量级想模拟 10¹⁸ cm⁻³ 的轻掺杂就要把超胞扩到几十个晶胞边长第一性原理计算的成本会迅速上升这也是很多掺杂建模只做趋势对照而不直接对齐产线浓度的原因。4.2 掺杂浓度、电阻率与方块电阻的换算电学参数之间有一条可直接用的关系链。电阻率 ρ 1/(q·n·μ_n)方块电阻 R_sh ρ/x_j其中 x_j 是结深。迁移率不是常数它随掺杂浓度上升而下降硅中电子迁移率常用 Caughey-Thomas 形式描述。import math q 1.602e-19 # 元电荷, C mu_max 1400.0 # 低掺杂电子迁移率, cm^2/(V·s) mu_min 65.0 # 高掺杂渐近电子迁移率, cm^2/(V·s) N_ref 8.5e16 # 参考浓度, cm^-3 alpha 0.72 # 拟合指数 def mu_n(N): 硅中电子迁移率随掺杂浓度的经验模型 return mu_min (mu_max - mu_min) / (1 (N / N_ref)**alpha) def sheet_resistance(N_D, x_j_um): 由掺杂浓度和结深估算方块电阻, 单位 Ω/sq n N_D mu mu_n(n) x_j x_j_um * 1e-4 # um - cm rho 1.0 / (q * n * mu) # Ω·cm return rho / x_j for N in [1e16, 1e18, 1e19, 1e20]: print(fN_D{N:.0e} cm^-3 mu{mu_n(N):.1f} cm^2/Vs R_sh{sheet_resistance(N, 0.2):.1f} Ω/sq)参数说明N_ref 和 alpha 决定迁移率随浓度衰减的快慢不同材料、不同载流子类型要换不同的拟合参数x_j 单位换算成厘米后才能和电阻率一起给出方块电阻。这段计算的意义是把掺杂浓度、迁移率退化和工艺可测的方块电阻连起来实测方块电阻偏高时可能来自激活率不足、结深偏浅或者补偿度过大而不是掺杂浓度本身写错了。4.3 Hall、CV、SIMS 各测什么别把化学浓度当电活性浓度掺杂表征方法各管一段信息混用会得出错误结论。方法直接测量量给出的是化学浓度还是电活性浓度主要限制霍尔 Hall载流子浓度、类型、迁移率电活性浓度面密度需要良好欧姆接触受表面态影响CV载流子浓度-深度分布电活性浓度可分辨深度需要肖特基结或 PN 结靠耗尽近似SIMS元素浓度-深度分布化学浓度含未激活破坏性标样依赖强扩展电阻 SRP载流子浓度-深度分布电活性浓度探针接触、结深附近分辨率下降四探针方块电阻电学综合量需配合结深才能反推浓度一条实用原则SIMS 剖面给出的是总杂质浓度霍尔和 CV 给出的是电活性浓度两者相除就是激活率。如果 SIMS 显示杂质很多但霍尔测到的载流子明显偏少问题通常出在退火不充分、杂质落在间隙位置或者形成了团簇而不是扩散或注入的剂量算错了。5. 用实测方块电阻反推掺杂浓度与补偿比工程上更常见的场景是只有电学测试数据要从方块电阻和结深反推电活性浓度。区块电阻和结深、迁移率的关系可以直接反解得到的是有效电活性面浓度。q 1.602e-19 R_sh 120.0 # 实测方块电阻, Ω/sq x_j 0.35 # 结深, um mu 200.0 # 由霍尔测试得到的迁移率, cm^2/V·s x_j_cm x_j * 1e-4 N_eff 1.0 / (q * mu * R_sh * x_j_cm) # 有效电活性体浓度, cm^-3 print(f有效电活性浓度 ≈ {N_eff:.3e} cm^-3)这里的 N_eff 是平均到结深上的有效浓度实际情况并非均匀分布高浓度集中在表面附近时这个平均值会低于峰值浓度。更精细的做法是把 SRP 或 CV 得到的浓度剖面分段积分每段用对应浓度下的迁移率再求和得到面浓度和方块电阻能较好地还原剖面形状。补偿比是排错时最容易被忽略的一项。n 型区的电活性浓度不是 N_D而是有效值 N_D − N_A。当注入 n 型杂质的区域本身有较高背景 p 型掺杂或者工艺引入了反型杂质实测载流子浓度会明显低于预期但 SIMS 上的施主总量却是对的这种偏差会随着补偿比接近 1 而急剧放大。遇到激活率看起来很低的情况先算一遍背景掺杂浓度再判断是不是补偿问题。温度相关参数也要一起看。迁移率随温度下降而上升本征载流子浓度随温度上升而快速上升两条曲线方向相反导致高温下用固定迁移率反推的浓度会偏大。可靠做法是在器件工作温度下重测霍尔用实测迁移率而不是查表值去反推浓度这一步能把中高温下的浓度偏差压下来。最后一点落在剖面验证上反推得到的浓度用 SIMS 剖面交叉核对确认峰值位置和积分剂量两块数据对得上掺杂工艺窗口才算真正锁定。本文还有配套的精品资源点击获取