半导体器件物理习题答案怎么用：从推导链到参数速查脚本

简介《半导体器件物理》第三版课后习题答案PDF面向微电子学、半导体物理方向的本科生、研究生及课程复习备考者针对施敏教材课后习题缺少系统解答的问题提供可直接对照的解题参考。压缩包内仅含1个PDF文件约3.13MB便于在电脑或移动端查阅。已有10417人学习说明其在相关课程复习中具有一定参考价值。内容围绕晶体结构、能带理论、杂质掺杂、载流子迁移率与扩散以及pn结等核心章节展开涵盖原胞划分、金刚石结构最近邻距离、体心立方基矢与倒格子、无限深势阱量子态、施主杂质电离能、金与硼掺杂效应、电阻率与迁移率关系以及pn结耗尽层电容、耗尽层宽度、内建电势、电场分布和击穿电压等典型题目并给出公式代入与计算过程适合用于核对思路、巩固概念和考前查漏补缺。1. 一份“答案”真正的用法从半导体器件物理习题到可复算的推导链很多人下载课后习题答案第一反应是对数字第 3 题 0.25a³、第 4 题 d√3a/4、第 12 题 E_C−E_F0.046 eV对上了就合上文档。但这份《半导体器件物理》第三版施敏著的习题解答价值不在答案那一行而在每道题前面那几步建模动作——原胞怎么切、势阱边界怎么设、耗尽区往哪边扩。它面向的是微电子、集成电路方向的本科生和研究生也适合已经工作、但pn结电容和击穿电压只在仿真器里见过的人。第 1 章晶体结构、第 2 章 pn 结、第 4 章 MOS 电容、第 5 章双极晶体管四块内容串起来正好是一条从材料参数到器件特性的推导链能当复习提纲也能当仿真前的验算底稿。2. 晶体结构与能带把几何题还原成参数题第 1 章的题目看起来像立体几何实际上是让你把“晶格常数 a”这个唯一输入换算成原子间距、原胞体积、倒格子基矢、态密度有效质量这些后续章节要反复调用的中间量。答案里那些 0.25a³、√3a/4、b4π/a如果只是抄下来到第 2 章算耗尽层宽度时你会发现 N_C、N_V 的来源是断的。这一章的正确读法是把每道题的输出当成一个可以写进脚本的常量。2.1 原胞体积与最近邻距离别急着代 a5.43Å第 3 题的解法是先切分原胞把它拆成一个边长 √2a/2 的正十二面体加两个正四面体再分别算体积求和得到 0.25a³。这个思路的关键不是拆分本身而是告诉你金刚石结构的原胞里实际只含 8 个原子所以每原子占体积是 0.25a³/8 a³/32。第 4 题则直接用立方体对角线推最近邻距离 d√3a/4室温下硅的 a5.43 Å 代进去得 d≈2.35 Å这个数字和后面算原子体密度、做紧束缚近似时的跃迁积分都直接相关。我一般会在文档边上留一张这样的对照表把散在各题里的常量统一到同一套单位物理量表达式硅的取值a5.43 Å后续用在哪原胞体积0.25a³约 4.00×10⁻²³ cm³原子体密度换算最近邻距离√3a/4约 0.235 nm键长、跃迁积分倒格子基矢模4π/a约 4.63×10⁸ cm⁻¹布里渊区边界原子体密度8/(0.25a³)约 5.0×10²² cm⁻³掺杂浓度量级参照注意表里第四行掺杂浓度 10¹⁵~10¹⁹ cm⁻³ 相对原子密度 10²² 只是千分之一到百万分之一这就是“轻掺杂近似”为什么在多数题里成立。把这句话记住后面第 2 章做单边突变结近似时就不会心虚。2.2 倒格子与态密度第 6 题和第 9 题是一条线第 6 题给出体心立方基矢 a₁0.5a(ij−k)、a₂0.5a(−ijk)、a₃0.5a(i−jk)按定义求倒格子基矢得 b₁0.5b(ij)、b₂0.5b(jk)、b₃0.5b(ik)结论是体心立方的倒格子是晶格常数为 b4π/a 的面心立方。这个结论本身要背但更重要的是它解释了第 9 题电子被限制在边长 a 的立方体内k 空间一个量子态占据的体积是 (π/a)³再对 4πk²dk 积分并乘上能谷数 M_C才得到 N(E) ∝ M_C·√(2m₀)··· 那个形式的态密度。把这两题连起来看会发现一个常被忽略的点态密度公式里的 m₀ 在严格推导中应换成有效质量 m*而 m* 的数值恰恰来自能带在极值点附近的曲率。习题答案里保留 m₀ 是因为原书在该章做了简化处理真正做器件仿真时我一般用一维泊松—漂移扩散自写求解器或商用 TCAD必须换成各能谷的有效质量并乘上能谷简并度。如果你的仿真结果和手算差一个量级八成就是这一步没换。2.3 施主电离能与载流子浓度0.044 eV 怎么变成 4.1×10¹³ cm⁻³第 12 题给磷的施主电离能 E_C−E_D0.044 eV在 77 K 下本征激发基本为零电离的施主电子数近似等于导带电子数联立电中性条件和费米分布求得 E_C−E_F0.046 eV导带电子浓度约 4.1×10¹³ cm⁻³。这道题的实操价值在于它演示了低温下“杂质未完全电离”的处理流程而这个流程在分析低温器件特性、红外探测器工作点时天天用。复算时最容易翻车的地方是费米分布的近似。77 K 时 kT≈0.0066 eVE_C−E_F0.046 eV 约为 7 个 kT用玻尔兹曼近似误差还能接受但如果温度降到 20 K 甚至更低E_C−E_F 与 kT 同量级必须回到费米—狄拉克积分此时简并因子 1/2 也要计入。文档里没有写这一步我建议自己补一个判断当 (E_C−E_F)/kT 3 时就别用指数近似了。用 Python 把这道题的迭代过程复现一遍顺便留个可改参数的模板比对着答案看十遍有用import numpy as np from scipy.optimize import brentq # 物理常量SI 与 cm 混用注意单位标注 k 8.617e-5 # eV/K q 1.602e-19 # C T 77.0 # K, 低温工况 Nd 1e16 # cm^-3, 施主浓度题中可改 Nc300 2.8e19 # cm^-3, 硅 300K 导带有效态密度 Ed 0.044 # eV, 施主电离能 # 态密度按 T^(3/2) 缩放 Nc Nc300 * (T/300)**1.5 # 电中性n Nd Nd / (1 2*exp((Ef-Ed)/kT))n 用玻尔兹曼 def neutrality(Ef): n Nc * np.exp(-(Ef - 0.0)/k/T) # 以 Ec0 为参考 Nd_ion Nd / (1 2*np.exp((Ef Ed)/k/T)) return n - Nd_ion Ef brentq(neutrality, -0.2, 0.2) n Nc * np.exp(-Ef/k/T) print(fEc-Ef {-Ef*1000:.1f} meV, n {n:.2e} cm^-3)这段代码做了三件事用 T^1.5 缩放 N_C用含简并因子 2 的电离率公式算电离施主数用 brentq 求电中性方程的根。参数说明Nd改成题给值就能对上 4.1×10¹³ 这个量级Ed换成其他杂质如砷 0.049 eV、硼受主 0.045 eV时要把符号和简并因子一起改T降到 30 K 以下时玻尔兹曼近似会失效需要把np.exp换成费米积分。第 16 题讲金和硼掺杂对电子、空穴数目的影响本质是同一套电中性方程多写两个杂质能级把neutrality函数里加一项受主电离项即可结构不用动。3. pn 结从耗尽层宽度到击穿电压的完整算例链第 2 章的题量最大覆盖耗尽电容、内建电势、电场分布、产生—复合电流、雪崩击穿答案里大量出现 W、E_m、V_bi、J_gen 这些符号。这一章如果只挑几道题对答案很容易漏掉一条隐藏的主线所有题目都建立在“耗尽区近似”和“突变结/线性缓变结”这两个假设上假设一换公式全变。读懂这一章等于拿到了一张判断“什么条件下能用哪个公式”的地图。3.1 内建电势与耗尽层宽度0.8351 V 是怎么来的第 4 题给 ψ_bi (kT/q)·ln(N_A·N_D/n_i²)代入数值得 0.8351 V。这个式子看着简单但 N_A、N_D 的取值直接决定后面 W 的量级。第 5 题进一步用泊松方程分两段求解电场边界条件 E₁(0)−E_m、E₂(W_D)0、两段在 x₁ 处连续再联立电压积分式得 W_D1.61 μm、E_m5.217×10⁴ V/cm。这套推导流程是 pn 结章节的标准动作建议至少手推一遍因为仿真软件里看电场分布时你就是靠这套直觉判断曲线形状对不对。第 5 题还有一个值得单独拎出来的细节反向电流由势垒区产生电流和扩散电流两部分组成答案算得 J_gen3.96×10⁻⁷ A/cm²而扩散电流 J_diff≈5.372×10⁻¹¹ A/cm²前者大四个数量级。这解释了为什么硅二极管在室温下的反向饱和电流往往由产生—复合主导而不是少子扩散。如果你的器件模型里只写扩散电流反向特性会偏小几个量级。复算这里最容易出错的是单位。ε_s 用 11.9×8.854×10⁻¹⁴ F/cmq 用 1.6×10⁻¹⁹ C浓度用 cm⁻³那么 W 算出来自然就是 cm。我见过不少人把 ε_s 代成 8.854×10⁻¹² F/m结果 W 差了两个数量级还以为是公式错。下面是把第 5 题前半段写成函数的样子import numpy as np eps0 8.854e-14 # F/cm eps_s 11.9 * eps0 # 硅的介电常数, F/cm q 1.6e-19 # C kT 0.0259 # eV, 300K ni 1.5e10 # cm^-3 def pn_junction(NA, ND, Vr0.0): 单边近似不成立时用完整式Vr 为反向偏压正值 Vbi kT * np.log(NA * ND / ni**2) V Vbi Vr W np.sqrt(2 * eps_s * V * (NA ND) / (q * NA * ND)) # 峰值电场在冶金结处 Em np.sqrt(2 * q * V * NA * ND / (eps_s * (NA ND))) return Vbi, W, Em Vbi, W, Em pn_junction(1e16, 1e16, 0) print(fVbi{Vbi:.3f} V, W{W*1e4:.3f} um, Em{Em:.2e} V/cm)参数说明NA、ND相等时是对称结第 5 题的数值大致落在这个区间要复现单边突变结如 p⁺n 结就把NA调到 1e19 以上此时 W 主要由轻掺一侧决定公式可退化成 W√(2ε_sV/(qN_D))。Vr改成 30 或 100 就是第 11、13 题的反偏情形注意反偏时 W 随 √V 增长、E_m 同步上升而击穿就发生在 E_m 达到临界场硅约 3×10⁵ V/cm的那一刻。3.2 产生电流与扩散电流反向电流到底由谁主导第 5 题后半段、第 6 题、第 7 题都在算反向电流答案反复出现 J_genq·n_i·W/τ_g 这个形式。它的物理含义是反偏时耗尽区里的产生中心不断发射载流子产生的电流正比于耗尽区体积也就是正比于 W。这与扩散电流正比于 n_i²/N 的规律完全不同所以两者的温度依赖也不一样——产生电流随 n_i 变化∝exp(−E_g/2kT)扩散电流随 n_i² 变化∝exp(−E_g/kT)高温下扩散电流才会追上并超过产生电流。第 6 题用 E_t−E_i0.02 eV 的深能级算产生率得 J_gen 与 τ 的关系再结合第 7 题反偏下 W 的扩展得到 J_gen≈2.758×10⁻⁹ A/cm² 这个量级。这几道题串起来其实在教你一件事反向漏电不是“一个数”而是温度、偏压、少子寿命三个参数的函数。做低漏电设计时降低 W提高掺杂和增大 τ改善工艺洁净度是两条独立路径这道题就是给这两条路径做定量估算的。3.3 雪崩击穿V_BD 与掺杂浓度的对应关系第 8、9、11、13 题都在算击穿答案里出现 V_BD≈60×(E_g/1.1)^{3/2}×(N/10¹⁶)^{−3/4} 这类经验式以及电离积分 ∫αdx1 的判据。这里有一个必须建立的概念击穿电压不是材料常数而是掺杂浓度的函数轻掺杂一侧的浓度决定 V_BD因为耗尽区主要向轻掺侧扩展、电场峰值也由它承受。第 13 题给 V_BD32 V反推 E_m4.409×10⁵ V/cm、W_D1.452 μm并画出电场随 W 线性下降的三角形分布再对三角形面积积分得击穿电压 45.76 V。这个“三角形面积等于电压”的做法非常实用拿到任何突变结的电场分布图肉眼估面积就能判断 V_BD 量级和仿真结果对不上时也能快速定位是电场峰值错了还是耗尽区宽度错了。第 12 题则给出另一条判据当 N 型外延层减薄到 5 μm 时外延层在击穿前就完全耗尽这就是穿通punch-through此时击穿电压由外延层厚度而非掺杂决定。实际功率器件设计里外延层厚度和掺杂是一对要同时定的参数只按雪崩公式算 V_BD 会高估耐压这一题的提醒价值就在这。4. MOS 电容与双极晶体管把 C-V 曲线和基区输运算明白第 4 章 MOS 电容和第 5 章双极晶体管是这份答案里计算量最重的两部分。MOS 部分从平带电压、阈值电压一路算到高频/低频 C-V 曲线的特征点BJT 部分则围绕基区渡越、发射效率、穿通与雪崩的竞争展开。这两章的题如果只对最终数值会错过最有用的东西——它们其实是一套“给定工艺参数预测器件电学特性”的手算流程和 TCAD 仿真前的 sanity check 完全对应。4.1 平带电压与阈值电压氧化层电荷怎么进公式第 4 章第 1 题算理想 MOS 电容的金属—半导体功函数差、氧化层电荷量进而求平带电压 V_FB 和阈值电压 V_T。第 10 题给 d15 nm、V_T0.7 V反推衬底掺杂第 16 题给 C2.5×10⁻⁸ F/cm²反推氧化层厚度和掺杂。这一串题的核心公式只有两个V_FBΦ_ms−Q_ox/C_oxV_TV_FB2ψ_B√(2ε_s qN_A·2ψ_B)/C_ox。容易忽略的是 Q_ox 的分布位置。第 15 题专门讨论了氧化层内薄层电荷对平带电压的贡献 dV_FB−xρ(x)dx/ε_s结论是同样总量的电荷离半导体界面越近影响越大。这解释了工艺上为什么对界面态和可动离子那么敏感——不是总量问题是位置问题。做 C-V 仿真时如果 V_FB 偏移和实测对不上先查电荷分布假设再查掺杂。第 6、10 题都用 C_min/C_i≈0.2~0.5 这个比值反推掺杂这个比值来自耗尽层电容与氧化层电容的串联。实际测量中 C_min 受频率影响明显低频时反型层能响应C_min 会升高高频时反型层来不及响应C_min 由耗尽层最大宽度决定。答案是按理想高频情形算的实测曲线如果 C_min 偏高先确认测试频率而不是怀疑掺杂算错。4.2 基区渡越与穿通电压13.6 V 这个数说明了什么第 5 章第 2 题文档中“小 2”那部分算共基极连接下的集电结耐压结论是先发生穿通穿通电压 13.6 V而不是雪崩击穿。推导分两步先按突变结算发射结、集电结在基区一侧的耗尽层宽度0.2 μm 和 0.068 μm再令两者之和等于中性基区宽度反解出穿通电压然后用 N_C5×10¹⁵ cm⁻³ 查雪崩击穿曲线得约 90 V比较后确认穿通先到。这道题的工程意义比数值本身大基区越薄、掺杂越轻穿通电压越低而且穿通是“软”失效——电压再升高电流迅速增大但不像雪崩那样陡峭。设计高频 BJT 时基区必须做薄所以穿通电压往往就是耐压瓶颈这一题给的就是这个矛盾的定量版本。第 1 题“靳 1”还算了两结耗尽层在基区的总占用 W0.714 μm以及少子电荷 Q2.89×10⁻¹³ C这俩量直接对应基区存储电荷和开关速度。基区少子分布那道题文档里“唐 3”给出了一个很实用的近似把中性基区掺杂看成线性分布少子面密度就是三角形面积 ½·p_n(0)·W再乘 qA 得电荷。这个近似在掺杂梯度不大时误差可接受但如果基区掺杂从发射结到集电结变化超过一个数量级最好还是解扩散方程。我一般先用三角形近似估量级再用数值解校正两步都做一遍心里才有底。5. 避坑与排查这份答案最容易用错的五个地方文档是扫描整理来的公式排版有大量上下标丢失、符号错位加上原书本身对不同章节做了不同近似直接照抄很容易在关键处翻车。下面五条是我在复算过程中真实踩过的按“现象→原因→解决”写遇到类似情况可以对着排。现象一第 12 题按答案代公式算出的 n 和 4.1×10¹³ 差两个量级。原因答案里 E_C−E_F0.046 eV 是在 77 K、且假定杂质完全电离或部分电离的特定条件下得到的如果直接套 300 K 的 N_C 值指数项 exp(−0.046/0.0259) 和 exp(−0.046/0.0066) 差了 e⁶ 倍。解决先确认题目温度N_C 必须按 T^{3/2} 缩放kT 用对应温度的值再检查是否考虑了简并因子。现象二第 5 题耗尽层宽度算出来比答案大一个数量级。原因ε_s 单位混用或 N_A、N_D 用了 m⁻³ 而 q 用 C、ε 用 F/cm量纲不闭合。解决统一到 cm-g-s 制ε_s11.9×8.854×10⁻¹⁴ F/cm浓度 cm⁻³算完 W 的单位就是 cm用 Python 写函数时把单位写进变量名如 eps_s_cm比事后检查省事。现象三MOS C-V 曲线的 C_min 和答案对不上。原因答案按理想高频、无界面态计算实际曲线受测试频率、界面态密度、串联电阻影响另外文档第 6 题的 C_min/C_i≈0.5 和第 10 题的 0.2 是不同掺杂下的取值不能混用。解决先确认掺杂对应的理论 C_min再看测试频率是否足够高一般 1 MHz最后才怀疑工艺参数。现象四击穿电压手算和仿真差 30% 以上。原因手算用的是突变结加经验电离率公式仿真里掺杂分布是缓变的电场峰值位置和形状都不同另外击穿判据用 ∫αdx1 时α 的经验系数在不同文献里能差一倍。解决手算只用来定量级和趋势精确值交给仿真如果非要手算对齐把掺杂分布按实际外延参数改为线性缓变结公式要换一套。现象五BJT 穿通电压算出来是 13.6 V但手册上同结构器件耐压 30 V。原因答案里基区宽度 W_B1.0 μm、N_B10¹⁶ cm⁻³ 是题设值实际器件基区更厚或掺杂更高而且穿通电压对 W_B 是平方依赖W_B 差 20% 穿通电压就差 40%。解决先核对题设参数和实际结构是否一致再检查是否漏掉了基区两侧耗尽层的温度依赖V_bi 随 T 下降高温下穿通电压更低。提示文档中第 2 章和第 4 章有多处公式因扫描识别出现符号缺失如根号、下标复算时以原书公式为准答案只用来核对最终数值量级。6. 进阶用法把答案改造成自己的参数速查脚本把这份答案用出价值的方式不是把它当标准答案对而是把每类题提炼成一个带默认参数的函数形成自己的“手算库”。我现在的习惯是晶体结构类题目留下几何常量和密度换算pn 结类留下 W、E_m、V_BD 三个函数MOS 类留下 V_FB、V_T、C_min 三个函数BJT 类留下穿通电压和基区电荷两个函数。每学一章就往里加半年下来这个脚本比任何速查表都好用因为每个参数都是自己推过一遍的。下面这段是把前面几章的公式收拢成一个可调用模块的骨架重点看参数默认值和单位约定具体数值按你的题设改import numpy as np EPS0 8.854e-14 # F/cm EPS_S 11.9 * EPS0 # 硅 Q 1.6e-19 # C KT300 0.0259 # eV NI 1.5e10 # cm^-3, 300K def builtin_potential(NA, ND, T300.0): 内建电势T 可改以支持低温/高温工况 kT KT300 * T / 300.0 return kT * np.log(NA * ND / NI**2) def depletion_width(NA, ND, Vr0.0, T300.0): 突变结耗尽层总宽度Vr 为反向偏压正值 V builtin_potential(NA, ND, T) Vr return np.sqrt(2 * EPS_S * V * (NA ND) / (Q * NA * ND)) def breakdown_rough(ND): 单边突变结雪崩击穿电压粗估ND 为轻掺侧浓度 cm^-3 return 60.0 * (1.12 / 1.1)**1.5 * (ND / 1e16)**(-0.75) def mos_threshold(NA, Cox, Qox0.0, phi_ms0.0, T300.0): N 型衬底简化阈值电压Cox 单位 F/cm^2 kT KT300 * T / 300.0 psi_B kT * np.log(NA / NI) Qdep np.sqrt(2 * EPS_S * Q * NA * 2 * psi_B) return phi_ms 2 * psi_B Qdep / Cox - Qox / Cox参数说明几个关键点depletion_width的Vr在正向偏置时要传负值函数内部不做符号保护调用前自己确认breakdown_rough是经验式只对单边突变结、掺杂 10¹⁴~10¹⁷ cm⁻³ 范围有效超出范围误差会很大mos_threshold里的Qox符号取决于电荷极性固定电荷为正时平带电压向负方向偏移这一点和第 15 题的推导一致别把符号弄反。这套脚本的用法不是替代仿真而是给仿真结果做第一道过滤仿真出来的 W 和depletion_width差一个数量级先查单位V_BD 和breakdown_rough差两倍以上先查掺杂分布是不是缓变V_T 和mos_threshold对不上先查氧化层厚度和 Qox 假设。最后说一个我自己的教训。早些年做一道 MOS 电容的题C_min 怎么算都和答案差 20%折腾一晚上最后发现是把氧化层厚度 15 nm 当成了 15 Å 代进去差了十倍只是串在其他近似误差里不明显。从那以后我每次写这类计算脚本第一行都强制加单位注释变量名带后缀函数入口先 assert 一遍量级。血泪经验就是半导体器件物理的题错基本不在物理在单位和近似条件的适用边界。希望帮到你。本文还有配套的精品资源点击获取