Semiconductor Fundamentals这个系列写到这里前面几篇花了不少篇幅把能带、态密度、掺杂这些“静态”概念讲清楚了。说实话能带图画得再漂亮掺杂浓度算得再准如果不知道载流子到底怎么运动、怎么产生、怎么消失那在面对二极管、MOSFET、IGBT这些真实器件时还是会一头雾水。Part 4我要做的就是把这几块拼图补上载流子的漂移、扩散、产生复合以及PN结如何把这些机制统合起来。这篇内容适合微电子专业的在校生、刚入行做器件设计的工程师以及做材料工艺想回头补电路物理基础的朋友——读完以后你再看器件手册里的耐压、恢复时间、导通电阻这些参数会明显感觉它们不再是孤立的数据。我需要先把话放前面这一篇不是泛泛的概念罗列而是会讲到可以直接用在估算里的公式、数值案例以及我自己在实验和仿真中踩过的真实坑。尤其数值估算那块建议你拿纸笔跟着算一遍算完对半导体物理的感觉会完全不一样。1. 载流子从哪来重新理解本征激发和掺杂1.1 本征激发为什么ni随温度飙升半导体里的载流子不是凭空出现的。理想的本征半导体里电子-空穴对只能靠热激发跨越禁带产生这一过程叫本征激发。本征载流子浓度ni是所有后续计算的基石它的表达式是ni² Nc · Nv · exp(-Eg / kT)注意这里有个平方因为电子和空穴是成对产生的n p ni。为什么是指数关系因为要跨越禁带必须从晶格振动中获得至少Eg的能量而晶格振动能量分布服从玻尔兹曼分布所以概率的指数项就是exp(-Eg/kT)。我在实际估算里经常会用300K下硅的ni ≈ 1.5×10¹⁰ cm⁻³这个值。但你要记住这个数只在室温附近成立。温度从300K升到400K指数项会发生剧烈变化——具体来说ni会从10¹⁰量级跳到10¹²量级整整两个数量级。这就是为什么功率器件在高温下漏电会显著增加因为它本质上就是你挡不住的热激发。1.2 掺杂后的载流子从“完全电离”开始本征半导体的载流子浓度太低没法直接做器件所以必须掺杂。掺入五价元素磷、砷、锑形成N型半导体掺入三价元素硼、镓、铟形成P型半导体。核心假设是“完全电离”在室温下施主杂质几乎全部给出电子受主杂质几乎全部接受电子。这个假设的物理基础在于施主能级通常位于导带底下方几十毫电子伏。比如磷在硅中的电离能大约0.045 eV而室温kBT ≈ 0.026 eV热运动能量完全足够把束缚电子“踢”进导带。所以工程计算里我们常用n ≈ ND N型p ≈ NA P型低温下这个假设会失效出现“载流子冻结”现象这也是为什么很多器件在低温下电阻反而变大。但在常温设计里完全电离是足够好的近似。1.3 补偿掺杂不是只有净浓度那么简单很多初学者不知道实际芯片里往往同时存在施主和受主杂质。浓度高的那一种决定导电类型这种工艺叫补偿掺杂。补偿掺杂最经典的应用就是半绝缘衬底或者在MEMS工艺里做电阻率精准调控。补偿掺杂的关键是漂移区的电阻率由净浓度决定但载流子寿命和迁移率会受到总杂质浓度的影响。比如N型区里Na和Nd都是10¹⁶ cm⁻³净浓度是0看上去像本征材料但实际上散射中心数量是2×10¹⁶迁移率比真正的本征材料低不少。这个坑我在早期做器件仿真时踩过后面第7章还会细说。2. 载流子怎样动漂移和扩散的两套规则2.1 漂移电场下的定向运动迁移率μ怎么定义载流子在半导体里首先会做无规则热运动就像教室里的学生下课铃一响到处乱跑。但一旦加上电场整个群体就会沿着电场方向产生一个微小的净位移。电子逆电场方向走空穴顺电场方向走这个定向运动就是漂移。漂移速度v和电场E的关系是v_d μE迁移率μ的单位是cm²/(V·s)它衡量的是“单位电场下载流子能跑多快”。对于低掺杂硅室温下电子迁移率大约1400 cm²/(V·s)空穴大约470 cm²/(V·s)。电子之所以比空穴“跑得快”是因为导带底的有效质量小于价带顶的有效质量同样的电场下加速度更大。但要特别注意v μE只在低电场下成立。电场超过大约10⁴ V/cm时漂移速度开始饱和硅的饱和速度大约10⁷ cm/s。这意味着你加再大的电场载流子也快不到哪去这是所有高频器件绕不过去的物理上限。2.2 扩散浓度梯度驱动的“扎堆散开”如果半导体里载流子浓度不均匀它们就会自动从高浓度区向低浓度区扩散这和墨水滴进清水的原理完全一样。扩散不依赖电场依赖的是浓度梯度。电子和空穴的扩散电流密度可以写成Jn_diff qDn · dn/dx Jp_diff -qDp · dp/dx符号不一样因为电子带负电而电流方向定义为正电荷运动方向。扩散系数D的单位是cm²/s它和迁移率之间有一个非常重要的关系叫爱因斯坦关系D / μ kBT / q ≈ 0.026 V室温这个关系说明了什么说明扩散和漂移本质上是同一个微观机制的两面——都是载流子与晶格散射相互作用的宏观结果。你测量出迁移率就能直接推出扩散系数不用单独测两遍。这个关系在后面的PN结计算里会反复用到。2.3 合到一起看总电流方程和费米能级视角当电场和浓度梯度同时存在时总的电子电流就是漂移和扩散两项的代数和Jn qμn nE qDn · dn/dx空穴的表达式也类似。这里有个很实用的技巧如果引入“准费米能级”EFn和EFp电流公式可以改写为更简洁的形式Jn μn · n · dEFn/dx Jp μp · p · dEFp/dx言下之意载流子永远倾向于从高费米能级流向低费米能级无论是电场引起的能带倾斜还是浓度差引起的费米能级位置变化造成的驱动力在数学上是统一的。我后面算PN结时你会发现用费米能级的视角看问题很多复杂的边界条件会一下子清楚非常多。3. 迁移率不是常数背后是散射在作祟3.1 两种主要散射机制迁移率之所以不是常数本质上是因为载流子在运动过程中会不断被散射散射越多净漂移速度越小迁移率越低。最主要的两种散射机制分别是声学声子散射和电离杂质散射。声子散射来自晶格原子的热振动温度越高振动越剧烈散射越强所以这种散射导致的迁移率随温度升高而下降大致按T⁻³/²规律变化。电离杂质散射则来自掺杂原子失去电子后的带电位点掺杂浓度越高散射中心越多但这种散射在温度升高时反而减弱因为载流子跑得快经过带电中心时被偏转的时间更短。实际硅材料中两种机制叠加起来迁移率随温度的变化会出现一个峰值。通常在中低温区100K到200K电离杂质散射占主导迁移率随温度升高而上升到了室温以上声子散射占主导迁移率随温度升高而下降。记住这个趋势对理解器件的负温度系数非常重要。3.2 掺杂浓度对迁移率的影响与电阻率估算掺杂浓度直接影响迁移率这点在做器件设计时怎么强调都不为过。硅的电子迁移率在极低掺杂时可以到1400 cm²/(V·s)以上但当掺杂浓度提高到10¹⁹ cm⁻³时迁移率会掉到100到200 cm²/(V·s)左右降了一个数量级还要多。工程上常用Caughey-Thomas经验公式来近似描述这个依赖关系μ(N) μ_min (μ_max - μ_min) / [1 (N / N_ref)^α]这个公式的物理思想很直接低掺杂时迁移率趋近μ_max高掺杂时逐渐降到μ_min中间有个平滑过渡。为什么这个细节那么重要因为电阻率的公式是ρ 1 / (q · n · μ)如果你只记得n ND却拿低掺杂时的迁移率去算高掺杂区的电阻率误差能到好几倍。比如漂移区浓度1e16时迁移率约1200电阻率约0.5 Ω·cm而源漏区浓度1e19时迁移率只有200电阻率约0.003 Ω·cm。如果拿1400去算结果会偏小不少反映不了真实的导通损耗。3.3 实测迁移率霍尔测量和四点探针理论归理论做芯片的人还是要测实际数值。最经典的方法是霍尔效应测量给样品通电流同时加一个垂直磁场载流子会在洛伦兹力作用下偏向一侧形成横向霍尔电压。霍尔电压的极性决定载流子类型大小可以算出浓度结合电导率就能反推迁移率。四点探针则是测薄层电阻的标配工具。四根探针排成一条线外面两根通电流里面两根测电压这样能消除接触电阻的影响测出来的薄层电阻Rs更可靠结合已知结深就能得到方块电阻率。我在实验室里见过不少测量翻车的情况最常见的是没有避光——室温下光生载流子会让电导率显著升高测出来的电阻率偏低数据完全失真。另外样品表面如果有氧化层或者沾污四探针读数也会漂移所以测量前用稀HF去氧化层是常规操作。4. 产生与复合非平衡载流子的“生死循环”4.1 直接复合、SRH复合、俄歇复合半导体器件在工作时注入或者光生出来的非平衡载流子不可能永远存在它们最终会通过复合回到平衡态。复合路径主要有三条。第一条是直接复合导带电子直接跃迁回价带与空穴复合多余能量以光子形式释放。这件事在直接带隙半导体比如砷化镓里发生得很顺利所以LED和激光器都用直接带隙材料。但在硅这种间接带隙半导体里直接复合需要声子参与同时满足动量守恒概率非常低基本可以忽略。第二条是SRH复合也就是通过禁带中的缺陷能级做“中转站”。缺陷能级先俘获一种载流子再俘获另一种完成复合。这条路径是硅里最主要的复合机制因为硅晶圆里或多或少总有金属杂质、位错、界面态。工程上甚至会主动利用这点——比如在快恢复二极管里掺入金或铂人为制造深能级复合中心把少子寿命从微秒级压到几十纳秒从而缩短反向恢复时间。第三条是俄歇复合高浓度下载流子直接三体碰撞把多余能量交给另一个载流子而不是光子。在重掺杂区和高注入条件下这一机制变得重要。功率器件在大电流密度下效率骤降往往就是俄歇复合在作祟。4.2 少子寿命和扩散长度器件的“内存时间”少子的寿命τ是一个极端重要的参数它决定了非平衡载流子能存在多久。少子寿命不仅影响双极晶体管的电流增益和开关速度还影响太阳能电池的收集效率。扩散长度L则描述载流子在被复合前能扩散多远L √(D · τ)这个量简单好用。如果器件的特征尺寸远大于扩散长度那注入的载流子还没到达收集端就复合掉了效率自然上不去。所以做功率器件的都追求大少子寿命做高频开关器件的反而希望寿命短一点这就是为什么同样一种二极管的制造工艺会因应用方向不同而截然不同。4.3 测量少子寿命的常见方法少子寿命的实验测量最常用的是微波光电导衰减法用激光脉冲在样品表面激发非平衡载流子然后用微波反射探针监测电导率随时间衰减的曲线衰减时间常数就是有效寿命。准稳态光电导法则通过不同光照强度下的稳态电导来反推注入依赖的寿命曲线。这里有个容易忽略的问题测到的寿命是“有效寿命”它包含了表面复合成份。如果表面复合速度很高测出来的寿命会远低于体寿命。所以在测量前必须做表面钝化比如用热氧化硅片或者用碘溶液钝化否则数据没什么参考意义。我在带学生做测量时每次都会先确认表面状态不然数据离散到根本没法用。5. 第一次把N型和P型拼起来PN结如何工作5.1 平衡状态扩散、漂移与内建电场PN结是所有半导体器件的母体想理解它就得先看两种半导体刚接触时会发生什么。N区电子浓度高P区空穴浓度高一接触浓度差就出现了电子往P区扩散空穴往N区扩散。但是这些载流子并不是无限扩散下去——它们离开后会在结附近留下带正电的施主离子和带负电的受主离子形成一个空间电荷区也叫耗尽区。这个电荷区会产生一个从N指向P的电场电场推动的漂移电流恰好和浓度差推动的扩散电流方向相反。平衡时就达到一种动态平衡扩散过去多少漂移拉回来多少净电流为零。这时候整个结构里费米能级处处相等能带在结区发生弯曲。内建电势的大小由两边的掺杂浓度和本征浓度共同决定Vbi (kT / q) · ln(NA · ND / ni²)举个例子两边浓度都是1×10¹⁶ cm⁻³室温下Vbi大约0.75 V到0.8 V和硅PN结的经验值0.7 V对得上。5.2 正向和反向偏压非平衡下的电流给PN结加正向偏压相当于外电场削弱内建电场耗尽区变窄势垒降低于是扩散电流大幅增加。这个过程的物理图景是N区的电子被注入到P区P区的空穴被注入到N区结附近出现大量非平衡少子它们边扩散边复合形成正向电流。加反向偏压则相反势垒升高扩散电流被压到几乎为零剩下的就是少子被耗尽区边缘“扫”过去形成的反向饱和电流I0。理想二极管方程I I0 · [exp(qV / kT) - 1]能精确描述这个行为。反向饱和电流I0非常小和少子寿命、扩散系数以及本征载流子浓度都有关系而由于ni²随温度指数上升I0对温度极其敏感——这就是为什么二极管的反向漏电在高温下会失控。5.3 温度对PN结的影响PN结最让人头疼的工程特性就是温度稳定性。正向导通电压V_F随温度升高大约以-2 mV/°C的速率下降温度越高同样的电流需要的正向电压越低。反向漏电则相反随温度升高指数上升。这两个趋势叠加起来会导致双极型功率器件在高温下很容易产生热失控。我做功率模块测试时经常要处理“正温度系数还是负温度系数”这个问题。合理的器件设计会让导通压降和电流呈正温度系数让电流自动均匀分布避免局部热点。单纯从PN结物理来看负温度系数是天然倾向所以必须靠工艺和结构去纠正。6. 三个手算实例把公式变成具体数字6.1 实例1掺杂后载流子浓度和费米能级假设有一块N型硅掺杂浓度ND 1×10¹⁶ cm⁻³室温300Kni取1.5×10¹⁰ cm⁻³。根据完全电离假设n ≈ ND 1×10¹⁶ cm⁻³然后由质量作用定律p ni² / n (1.5×10¹⁰)² / 1×10¹⁶ 2.25×10⁴ cm⁻³这个空穴浓度非常低约是电子浓度的十亿分之一。费米能级相对本征费米能级的位置EF - Ei kT · ln(ND / ni)代入数值ln(1×10¹⁶ / 1.5×10¹⁰) ln(6.67×10⁵) ≈ 13.41乘kT约0.0259 eV得到0.347 eV。也就是说费米能级上移到了带隙中点上方约0.35 eV处距离导带底约0.21 eV。6.2 实例2电阻率估算迁移率用对了吗继续用ND 1×10¹⁶ cm⁻³的N型硅来算电阻率。查实验迁移率数据这个掺杂浓度下电子迁移率约1200 cm²/(V·s)而不是本征值1400。代入ρ 1 / (q · n · μ) 1 / (1.6×10⁻¹⁹ × 1×10¹⁶ × 1200) ≈ 0.52 Ω·cm如果拿1400去算结果是0.45 Ω·cm偏差看着不大。但把浓度提到1×10¹⁹ cm⁻³后迁移率掉到约200电阻率约0.0031 Ω·cm如果还用1400就是0.00045 Ω·cm差了近7倍。这种量级误差在重掺杂区和外延层设计里是致命的。6.3 实例3不对称PN结的内建电势与耗尽层考虑一个不对称PN结P侧NA 1×10¹⁷ cm⁻³N侧ND 1×10¹⁵ cm⁻³。先算内建电势Vbi 0.0259 × ln(1×10¹⁷ × 1×10¹⁵ / (1.5×10¹⁰)²)括号里是1×10³² / 2.25×10²⁰ ≈ 4.44×10¹¹取自然对数为26.82Vbi ≈ 0.695 V。然后算零偏耗尽层宽度用公式W sqrt(2ε_s · Vbi / q · (NA ND) / (NA · ND))硅的介电常数约1.04×10⁻¹² F/cm算出来W约0.95 μm。注意这个宽度几乎全部展在轻掺的N侧因为两侧耗尽层宽度之比等于掺杂浓度反比这里约1:100。轻掺侧承担了几乎所有的耐压这就是功率器件中漂移区要用低掺杂的原因——只有让耗尽层尽量展宽才能承受更高的反向电压。7. 我在学习和器件仿真中踩过的几个坑7.1 把ni当成常温常数结果一算高温就变形最开始做器件仿真我想当然地认为ni永远是1.5×10¹⁰ cm⁻³于是算一个125°C环境下的漏电时怎么都对不上实测数据。后来才意识到125°C下ni已经涨到接近10¹² cm⁻³量级反向饱和电流因为含有ni²项会放大到室温的几千倍。做温度特性设计一定要用完整的ni(T)表达式不能只看常温值。7.2 算电阻率时拿纯本征迁移率直接套这个是仿真报告里最常出现的错误。有一次审核同学的报告发现他算源漏区电阻率用了电子迁移率1400我一看注入浓度是1×10¹⁹ cm⁻³告诉他得用Caughey-Thomas公式或实验曲线查迁移率至少应该是200左右。他用错误迁移率算出来的导通电阻比实际小了近7倍整个方案结论都变了。7.3 漂移和扩散的符号问题在PN结电流推导里扩散电流和漂移电流的符号是最容易看走眼的。我早期算一个PN结构的正向电流方程里Jp那一项少了个负号结果电流方向反了仿出来完全是“负电阻”的假象。后来养成了一个习惯每写一项就停下来思考物理方向——空穴从P区到N区是正的它的浓度梯度方向是负的所以fp的扩散项必须带负号才能让电流为正。7.4 小注入假设被打破时SRH公式不能用SRH复合的标准公式建立在小注入条件下也就是非平衡少子浓度远小于平衡多子浓度。但我在做高压大电流仿真时发现晶体管基区在大电流下会进入大注入少子浓度甚至可以超过多子浓度此时复合率对注入水平的依赖完全变了寿命会随着注入浓度增加而上升。生搬硬套小注入公式算出来的增益和开关损耗都偏得离谱。这个坑我踩过两回后来再做器件仿真第一件事就是检查工作区是不是落在了小注入区间不是的话赶紧把SRH模型换成注入依赖的形式。半导体基础这部分看着不难但只要想把它用到真实器件里处处都是细节。