简并与非简并半导体：物理判据与工程实践

半导体物理这门课最容易被跳过的一页，往往就是"简并与非简并"这六个字。教材上通常一句话带过：掺杂浓度高了就是简并，用玻尔兹曼分布算不准了就换成费米-狄拉克分布。可真到了器件仿真、霍尔测试数据拟合、或者接触电阻调不下来的深夜，你会发现所有的麻烦都堆在这个"换公式"的动作上——载流子浓度被高估、本征浓度被低估、迁移率曲线对不上、带隙变窄没算进去，最后仿真和实测差出百分之几十，还找不到原因。这篇就围绕简并半导体和非简并半导体这条分界线，把判据的物理来源、数学工具的落地方案、材料性质在越线之后的连锁变化，以及器件里哪些地方必须主动把自己推进简并区，一条一条拆开讲。适合已经会算载流子浓度、但一碰到重掺杂就开始心里没底的读者，也适合做器件建模、需要判断模型开关该不该打开的人。

1. 从一条在高掺杂下悄悄失准的公式说起

1.1 那个让仿真结果偏离实测三成的假设

多数人第一次算载流子浓度，用的都是n = N_C · exp(−(E_C − E_F)/kT)这一条。低掺杂下它非常好用，误差可以忽略，算出来的数和霍尔测试对得上，于是很多人就把它当成普适公式一直用下去。问题出在掺杂浓度往上走的时候。

假设你在硅里掺了 1×10¹⁹ cm⁻³ 的施主，室温全电离，那 n 大约就是这个数，N_C 在 300 K 下是 2.8×10¹⁹ cm⁻³，于是 n/N_C ≈ 0.357。如果继续用玻尔兹曼公式反推，会得到 E_C − E_F = kT·ln(1/0.357) ≈ 1.03kT，也就是费米能级在导带底下约 26 meV，看起来还挺"正常"，离简并还远。但真实情况是，费米-狄拉克积分给出的解对应 η = (E_F − E_C)/kT ≈ −0.85，费米能级只低于导带边 22 meV 左右，而且这还是不考虑带隙变窄的结果。

这个差别在低掺杂时无关痛痒，但到了 3×10¹⁹ cm⁻³ 以上就会翻车。按玻尔兹曼算，n/N_C = 1.07 时公式直接给不出解——指数项大于 1 意味着 E_F 已经跑到导带里面去了，而玻尔兹曼分布在这个区域根本没有物理意义。更要命的是，很多人不会意识到公式失效，而是拿这个"算不出来"的结果去调别的参数，把错误一路传下去。

我自己踩过的最典型一次坑，是给一个重掺杂源漏的短沟道器件做标定。实测的方块电阻比仿真低了将近 25%，我先后怀疑了迁移率模型、接触电阻、甚至网格划分，折腾了两天才反应过来：仿真里的载流子统计用的是玻尔兹曼近似，而源漏区的掺杂已经到 5×10¹⁹ 量级，属于明确简并。换成费米-狄拉克统计之后，同样的掺杂浓度算出的电子浓度更低、费米能级位置更高，配合带隙变窄模型，方块电阻一下就对上了。这件事之后我养成一个习惯：只要看到掺杂浓度进入 10¹⁹ cm⁻³ 量级，第一件事就是把统计模型从玻尔兹曼切到费米-狄拉克，别等对不上数再回头找。

1.2 简并与非简并，本质是统计分布的选择题

把话说透一点：简并和非简并不是两种材料，而是同一种材料在不同载流子填充程度下呈现的两种统计状态。半导体里的电子是费米子，服从费米-狄拉克分布，占据概率是

f(E) = 1 / (1 + exp((E − E_F)/kT))

当 E − E_F 远大于 kT 时，分母里的 1 可以忽略，分布退化成经典的玻尔兹曼尾巴exp(−(E − E_F)/kT)，这就是非简并情形——载流子稀稀拉拉地分布在能带边缘，彼此之间几乎不"抢座位"，泡利不相容原理不起作用。当 E_F 逼近甚至进入能带内部时，1 这一项不能再扔，能带底附近的量子态被大量填满，载流子开始"挤"在一起，这就是简并。

打个生活化的比方：非简并就像深夜的高速公路，车少，每辆车想开多快开多快，互相不干扰；简并就像早高峰，车道被占满，你想加速也没空间，整体行为从"个体自由"变成了"集体拥堵"，表现出来反而更像金属——电导率高、温度系数小、费米能级钉在能带里不动。

这个类比能帮你记住一个关键点：简并的核心不是"掺杂多"，而是"量子态被填满"。掺杂只是实现填充的手段之一，低温、量子限域、强光注入、甚至电场导致的能带弯曲，都能把系统推进简并区。我见过有人在 77 K 下测一个中等掺杂的样品，发现数据怎么都对不上室温模型，最后发现是低温下有效状态密度按 T^1.5 缩水了七八倍，材料在低温下已经进入简并——掺杂浓度一点没变，状态变了。

2. 判据的物理含义与工程上真正好用的那条线

2.1 费米能级和能带边的相对位置才是唯一标准

判断简并与否，只看一件事：费米能级 E_F 相对能带边的位置。n 型看 E_C，p 型看 E_V。E_F 在禁带里、离带边足够远，就是非简并；E_F 贴近带边甚至进入带内，就是简并。掺杂浓度、温度、有效质量，全都只是影响 E_F 位置的输入参数，不是判据本身。

严格来说，简并程度用一个无量纲量描述最清楚：

η_n = (E_F − E_C) / kT （n 型） η_p = (E_V − E_F) / kT （p 型）

• η < −3 到 −4：深度非简并，玻尔兹曼公式误差可以接受
• −3 < η < 0：弱简并过渡区，玻尔兹曼开始明显高估载流子浓度
• 0 < η < 2：中度简并，E_F 已进入能带
• η > 2：强简并，载流子气接近金属行为

工程上常用的说法"E_C − E_F > 3kT 为非简并"，是把误差压到可接受范围的经验线，不是严格的数学边界。严格讲，玻尔兹曼近似和真实费米-狄拉克积分之间的相对偏差是随 η 单调变化的，η 越负偏差越小，但收敛得并不快。所以在精度要求高的场合（比如做少子寿命或者本征浓度反演），我会把这条线推到 4kT 甚至 5kT，宁可多算一步费米积分，也不要在边界上赌。

2.2 3kT 这条线到底怎么来的

这条线的来源其实很朴素。载流子浓度严格写是

n = N_C · (2/√π) · F_1/2(η)

其中 F_1/2 是 1/2 阶费米-狄拉克积分。当 η 很负时，F_1/2(η) → exp(η)，两者重合；当 η 往正走时，F_1/2 的增长速度明显慢于指数，因为它被态密度"卡住"了。两者的比值在 η = 0 处大约只有 0.76 左右，也就是说在带边位置，用玻尔兹曼公式会把载流子浓度高估三成左右。要让这个偏差压到百分之几的量级，η 需要退到 −3 以下。这就是 3kT 这条经验线的由来——它对应的是"公式误差可接受"，而不是"物理上开始简并"。

这一点特别值得强调，因为很多人把它理解反了。物理意义上的简并起点是 η ≈ 0，即费米能级到达带边；而工程上用 3kT 是因为再靠近就不好用近似了。两个概念差了 3 个 kT，在室温下差不多 78 meV，不是小数目。

2.3 不同材料的简并门槛差得很远

这是我特别想强调的一点：不存在一个通用的"掺杂浓度阈值"。门槛由有效状态密度 N_C、N_V 决定，而后者跟有效质量、能谷数、温度直接挂钩。同样是 300 K，不同材料的差距能到一个数量级以上。

看这张表能读出两件事。第一，p 型普遍比 n 型更容易简并，因为价带的有效状态密度通常更小（Si 里 N_V 只有 N_C 的三分之一多一点），同样的掺杂浓度在 p 型侧更容易把费米能级推到带边。做 p+ 接触的时候如果按 n+ 的经验去定掺杂，很容易发现实际早就简并了。

第二，宽禁带和化合物半导体的门槛低得多。GaAs 的 N_C 只有 4.7×10¹⁷，掺杂到 5×10¹⁷ 就已经进入简并区了。这解释了一个常见困惑：为什么 GaAs 的欧姆接触相对好做，而宽禁带材料即使掺到很高浓度也难做低阻接触——不是掺杂不到位，而是这些材料的掺杂固溶度、电离能、自补偿效应限制了实际能达到的载流子浓度，很多时候根本推不到简并区。

还有一层容易忽略：温度会改变门槛。N_C ∝ T^{3/2}，从 300 K 降到 77 K，N_C 缩水约 7.7 倍。一个在室温下稳稳处于非简并的样品，降到液氮温度就可能变成简并。低温霍尔测试如果还套用室温的玻尔兹曼公式去反推载流子浓度，结果会系统性偏高。

3. 把费米-狄拉克积分算明白：从状态密度到反解费米能级

3.1 有效状态密度的量级估算

N_C 的解析式是

N_C = 2 · (2π m_de* kT / h²)^(3/2)

m_de* 是态密度有效质量。这个式子的实用价值不在于精确计算（因为 m* 本身就有温度依赖和能谷修正），而在于让你能快速估量级。记住 300 K 下 Si 的 N_C ≈ 2.8×10¹⁹ cm⁻³ 这个锚点，其他材料按有效质量和温度缩放就行：有效质量翻倍，N_C 涨 2^1.5 ≈ 2.83 倍；温度翻倍，N_C 涨 2.83 倍。

实操中我一般不去手算 N_C，而是直接从材料参数表里取，然后把注意力放在后面的积分上。因为 N_C 的取值在不同手册之间能差 20% 以上（主要差在有效质量的选取和能谷简并度的处理），纠结这个不如把统计模型选对。

3.2 费米-狄拉克积分的数值处理

F_1/2(η) 的标准写法是

F_1/2(η) = (2/√π) · ∫₀^∞ x^(1/2) / (1 + exp(x − η)) dx

这个积分没有初等闭式解，实际处理有三条路：查表插值、解析近似、数值积分。我的建议是数值积分优先，因为现在算力根本不是瓶颈，而近似公式的适用范围各有各的坑。

几个常用的近似值得知道：

• 玻尔兹曼：F_1/2 ≈ exp(η)，只在 η < −3 时靠谱
• Nilsson 近似：F_1/2 ≈ 1/(exp(−η) + 0.27)，在 η = 0 附近给 0.787（真值约 0.765），跨整个过渡区都能用，误差在几个百分点
• 强简并渐近：η > 1.5 时 F_1/2 ≈ (4/(3√π))·η^(3/2)，这就是金属自由电子气的行为，n ∝ (E_F − E_C)^{3/2}

如果你只需要快速判断是否简并，Nilsson 近似足够了。但如果要做参数提取或者仿真对标，老老实实数值积分。

下面这段是我平时用的实现，直接可以跑：

import numpy as np from scipy.integrate import quad from scipy.optimize import brentq def fermi_half(eta): """1/2 阶费米-狄拉克积分 F_{1/2}(eta)""" f = lambda x: np.sqrt(x) / (1.0 + np.exp(x - eta)) val, err = quad(f, 0, np.inf, limit=200) return 2.0 / np.sqrt(np.pi) * val def solve_eta(n_over_Nc): """给定 n/N_C，反解 eta = (E_F - E_C)/kT""" target = n_over_Nc func = lambda eta: fermi_half(eta) - target # 搜索区间覆盖从深度非简并到强简并 return brentq(func, -30.0, 60.0, xtol=1e-10) # 硅 300K 试算 Nc_Si = 2.8e19 for Nd in [1e16, 1e18, 1e19, 3e19, 1e20]: ratio = Nd / Nc_Si eta = solve_eta(ratio) print(f"N_D={Nd:.1e}, n/N_C={ratio:.4f}, eta={eta:+.3f}, " f"E_C-E_F={-eta*0.02585*1000:+.1f} meV")

跑出来的结果是：1×10¹⁶ 时 η ≈ −8.4，深度非简并；1×10¹⁹ 时 η ≈ −0.85；3×10¹⁹ 时 η ≈ +0.34，费米能级进入导带约 9 meV；1×10²⁰ 时 η ≈ +2.1，高出导带约 54 meV，属于强简并。

这组数字非常值得记住，它给你一个直觉刻度：硅在室温下，掺杂到 1×10¹⁹ 还只是"逼近"简并，到 3×10¹⁹ 才算真正跨过去，而 1×10²⁰ 已经是强简并了。这和很多人口口相传的"1e19 就简并"有出入，差别就在于用的是什么判据。按物理定义（η > 0）是 3e19 左右；按"玻尔兹曼误差 5%"是 1e19 左右。看你问的是哪个问题。

3.3 别忘了电离率这一步

上面所有计算都默认"掺杂进去的杂质全部电离"，掺杂浓度等于载流子浓度。这个假设在重掺杂下反而更成立，因为杂质能级在重掺杂时会展宽成杂质带，与能带重叠，电离能趋近于零。所以简并区用 n ≈ N_D 是合理的。

但在中等掺杂加低温的情况下，这个假设会崩。比如硅里掺 1×10¹⁷ 的磷，室温全电离没问题，但降到 77 K 时电离率可能只有百分之几十，载流子浓度远低于掺杂浓度，此时再用 n = N_D 去套简并判据会得出完全错误的结论。这时候需要联立电中性方程：

n + N_A⁻ = p + N_D⁺ N_D⁺ = N_D / (1 + g_D · exp((E_F − E_D)/kT))

g_D 是基态简并因子，施主通常取 2，受主取 4。这个方程得数值求解，但逻辑很简单：给一个 E_F，算出 n、p、电离杂质浓度，代入电中性方程看是否平衡，二分或者牛顿迭代都行。

我吃过一次亏：一个样品在室温下测出载流子浓度 8×10¹⁶，按 77 K 的 N_C 折算觉得应该刚好进简并区，结果低温霍尔测出来只有 3×10¹⁶，差了不止一倍。查了半天以为是测试问题，最后是电离率没考虑——那个掺杂浓度下，77 K 的电离率只有 40% 左右。判断简并之前，先确认你的载流子浓度是真的载流子浓度，不是掺杂浓度。

4. 越过那条线之后，材料发生的连锁变化

4.1 带隙变窄与被"放大"的本征载流子浓度

这是简并区最容易被漏掉、影响又最大的一项。重掺杂之后，杂质之间、杂质与载流子之间的多体相互作用会让导带底和价带边分别向禁带内移动，宏观表现就是带隙变窄（Bandgap Narrowing，BGN）。

量级上，硅在 N = 1×10¹⁹ cm⁻³ 时带隙收缩大约 50 到 60 meV，到 1×10²⁰ 时能到 100 meV 以上。这个数字看起来不大，但对本征载流子浓度是指数级的影响：

n_i² = N_C · N_V · exp(−E_g / kT)

E_g 减小 60 meV，在室温下 exp(60/25.85) ≈ 10，也就是 n_i² 放大十倍，n_i 放大三倍多。少子浓度直接跟着涨。

这个效应在双极器件里是实打实的性能损失。BJT 的发射区重掺杂，带隙比基区窄，形成所谓的"带隙差"，导致发射效率下降、电流增益被压缩——这是重掺杂发射区的固有代价，不是工艺缺陷。太阳能电池的背面场、LED 的有源区重掺杂，也都要为这个付出代价。做器件仿真时如果模型里 BGN 的开关没打开，少子浓度会被严重低估，漏电流、增益、开路电压全都不对。

常用的 BGN 模型有 Slotboom 和 del Alamo 两套，形式都是

ΔE_g = A · [ln(N/N_ref) + sqrt((ln(N/N_ref))² + 0.5)]

硅的常用参数是 A ≈ 9.0 meV，N_ref ≈ 1×10¹⁷ cm⁻³。不同模型在 1×10²⁰ 以上的外推差异可以到 30 meV，这个量级在仿真里已经能造成明显的曲线偏移，选模型时要和你对标的数据源保持一致。

4.2 迁移率退化：散射机制的切换

简并区的迁移率会明显下降，原因有三层。

第一层是电离杂质散射增强。掺杂浓度上去，库仑散射中心密度跟着上去，迁移率按 N 的幂次下降。硅里电子迁移率从低掺杂的 1400 cm²/(V·s) 掉到 1×10¹⁹ 时的 100 多，掉了十倍不止。

第二层是载流子-载流子散射。载流子浓度高到一定程度，电子之间互相碰撞也成了有效的动量弛豫机制。这一项在低掺杂时可以忽略，在 1×10²⁰ 以上不能。

第三层是屏蔽效应。高载流子浓度会屏蔽电离杂质的库仑势，反而减弱了电离杂质散射——这是个反向效应，和第一层部分抵消。所以在中等掺杂区间迁移率下降得快，到了极高掺杂反而趋缓，曲线不是简单的幂律。

工程上描述这段曲线，Masetti 模型和 Arora 模型都常用，形式是多个指数项叠加。我自己的经验是：如果只是要个量级估算，记住硅在 1×10¹⁹ 时 μ_n ≈ 100–150，1×10²⁰ 时 ≈ 50–80，p 型还要再低一些。精确值必须查模型或者实测，因为工艺路径（是扩散还是外延、有没有补偿掺杂）对结果影响很大。

顺带说一句简并区里迁移率的物理含义。强简并时费米面附近的态才是导电的主力，迁移率的概念本身开始向金属的 Drude 模型靠拢，霍尔因子也不再是 1。这时候用低掺杂的迁移率-浓度关系去外推，误差会累积到接触电阻的估算上，而接触电阻又是对迁移率最敏感的参数之一。

4.3 光学性质：吸收边移动与复合路径改写

重掺杂对光学性质的影响，做光电器件的人绕不开。

Burstein-Moss 位移：n 型重掺杂时费米能级进入导带，导带底附近的态被填满，最低能量的光学跃迁不再是 E_g，而是要跳到 E_F 以上，表现为吸收边向高能方向蓝移。位移量约等于 E_F − E_C。这个效应在透明导电氧化物里是常态，也是判断材料是否简并的一个直接证据——测个透射谱，看吸收边是不是比本征带隙高出一截。

俄歇复合增强：俄歇复合率正比于 n²p 或 np²，在重掺杂区随浓度三次方增长，很快成为主导复合机制。LED 和激光器在重掺杂有源区效率掉得快，主因之一就是这个。做高效率发光器件时，简并是有代价的——你想要高电导，就得接受俄歇复合。

自由载流子吸收：简并半导体的载流子浓度接近金属，红外波段的自由载流子吸收变得显著，这在透明导电膜的设计里是个必须权衡的点：方阻要低就得掺杂高，掺杂高了红外透过率就下降。

5. 器件里主动把自己推进简并的几类场合

5.1 欧姆接触：简并不是可选项，是必需品

金属-半导体接触要想做成欧姆接触，核心是把势垒做得足够薄，让载流子能隧穿过去。隧穿概率对势垒宽度是指数敏感的，而势垒宽度反比于 √N。所以路径很明确：把半导体一侧重掺杂到简并区，势垒宽度压到几纳米，隧穿成为主导输运机制。

比接触电阻的经验量级：掺杂 1×10¹⁸ 时大概是 10⁻⁴ Ω·cm² 量级，1×10¹⁹ 时能到 10⁻⁶，1×10²⁰ 时进 10⁻⁷ 甚至更低。三个数量级的差距，就是"能用"和"性能瓶颈"的差距。

这里有个实操细节值得说：接触区的简并是局部的，不需要整个有源区都重掺杂。常见的做法是用离子注入在接触窗口做局部重掺，退火激活后再做金属化。这样做的好处是既拿到了低接触电阻，又避免有源区整体的少子寿命被重掺杂毁掉。我见过有人为了降接触电阻把整片衬底都重掺，结果器件的漏电流和噪声全崩了——接触到的地方简并就够了。

还有一个坑：重掺杂区的实际激活浓度往往达不到注入剂量。固溶度、团簇、自补偿都会吃掉一部分。设计时按"名义注入浓度"去估接触电阻，实测经常差一个数量级。所以接触工艺的窗口要按激活后的实际载流子浓度来定，这个必须用电化学电容电压或者霍尔测试去标定，不能拍脑袋。

5.2 隧穿结：两侧都得简并才有意义

多结叠层电池、隧穿场效应晶体管、隧道二极管，这些器件的共同点是依赖带间隧穿，而带间隧穿要求两侧都重掺杂到简并，耗尽区宽度压到几纳米，能带在极短距离内急剧弯曲。

以多结电池的隧穿结为例，典型结构是 n++/p++ 各 1×10²⁰ 量级，厚度十几纳米。设计上有两个相互拉扯的要求：掺杂要足够高以形成薄势垒，但厚度也要足够薄以减少吸收损耗。隧穿结的电阻如果做不好，整个叠层电池的填充因子就被它拖死，而且这个损失是串联进去的，前面几层做得再好也补不回来。

这里的简并还有一个副作用：带隙变窄会降低隧穿结的有效势垒高度，反过来又影响隧穿概率。所以隧穿结的建模不能只算能带图，要把 BGN 一起带进去，否则电流密度会算偏。

5.3 透明导电膜与简并半导体激光器

透明导电氧化物是简并半导体最商业化的一类应用。这类材料（比如氧化铟锡体系）通过高浓度掺杂实现简并，费米能级进入导带，因此既有接近金属的电导率，又在可见光波段保持高透过率。设计的核心矛盾是：载流子浓度越高，方阻越低，但等离子体边往可见光方向移动，短波透过率下降。工程上的平衡点通常在载流子浓度 10²⁰ 到 10²¹ 量级，迁移率尽可能高以在同等电导下降低所需的载流子浓度。这也是为什么这类材料的研究重点一直在"高迁移率"而不是"高浓度"上——同样方阻下，迁移率高的方案光学损失更小。

简并半导体激光器的逻辑不太一样。激光要满足粒子数反转条件，Bernard-Duraffourg 条件要求准费米能级差大于光子能量，这天然要求有源区处于简并或准简并状态。所以对激光器来说，简并不是负担，是工作前提。真正要权衡的是俄歇复合和自由载流子吸收带来的内损耗，这两个都和掺杂浓度正相关，于是设计上要找一个平衡点：掺杂太低达不到反转条件，太高内损耗吃掉增益。

6. 参数提取与仿真建模里那些容易判断错的地方

6.1 霍尔测试在简并区的系统性偏差

霍尔测试是提取载流子浓度和迁移率最常用的手段，但它在简并区有几个系统性的坑。

霍尔因子偏离 1。低掺杂时霍尔迁移率和漂移迁移率的比值 r_H 接近 1，简并时理论上趋近于 1 附近但散射机制变了，实际取值在 1 到 2 之间浮动，取决于主导散射。用 r_H = 1 去反推，浓度和迁移率会同时带上系统误差。

表面耗尽层的影响。重掺杂样品的载流子浓度高，但表面态引起的耗尽层如果相对薄膜厚度不可忽略，测出来的有效厚度偏小，浓度被高估。做薄膜样品时这一项尤其要注意，薄膜越薄影响越大。

接触质量。霍尔电极如果做不成欧姆接触，测出来的电压里混进了接触电阻的贡献，迁移率会被严重低估。这一点在宽禁带材料上尤其突出——前面说过它们的简并门槛低但掺杂难，接触往往是瓶颈。

我的做法是：简并样品至少用两种方法交叉验证，霍尔配一个电容电压法或者范德堡法，看一下浓度是否自洽。差得超过 20% 就说明有系统性偏差，得先把测试本身搞清楚，再谈材料参数。

6.2 仿真模型里几个必须确认的开关

做器件仿真，进入重掺杂区之前，我一般会按这张清单过一遍：

这张表里最容易被跳过的是"带隙变窄"和"俄歇复合"，因为它们不像统计模型那样有一个显眼的开关，很多人根本不知道要开。我建议的做法是：只要仿真结构里出现 1×10¹⁹ 以上的掺杂区，就把这两个默认打开，然后通过对比开关前后的结果差异来判断影响量级，而不是凭感觉决定。

6.3 几个反直觉但很实用的结论

第一，"掺杂高"不等于"简并"。前面算过，硅在 1×10¹⁹ 时按物理定义还没跨过去。宽禁带材料更难，因为掺杂本身就有固溶度和自补偿的限制。见到高掺杂就先假定简并，是另一种形式的偷懒。

第二，简并不一定是坏事。接触区、源漏、透明电极、隧穿结，全都靠简并吃饭。判断简并"好不好"要看它在器件里的位置，而不是一概而论。

第三，温度是会翻盘的变量。低温下 N_C 缩水，中等掺杂也可能简并；高温下 N_C 膨胀，重掺杂也可能退回非简并。做宽温区应用的器件，简并判据必须按工作温度逐点算，不能拿室温的结论一路套到底。

第四，简并区的"载流子浓度"和"掺杂浓度"必须分开谈。前者是能带电学量，后者是化学量，两者在重掺杂下的差距来自电离率和固溶度，在低温下的差距来自不完全电离。混淆这两个概念是重掺杂建模出错的高频原因。

第五，玻尔兹曼公式不会报错，它只会悄悄给你一个错数。这是我觉得最值得警惕的一点——非简并公式在简并区不会崩，它会返回一个看起来合理的值。所以判据这东西必须前置，不能等结果不合理了再回头查。

7. 我个人的几条实操习惯

写到这里，把前面散落的经验收一下。第一条是先算 η，再谈其他。拿到任何一个掺杂浓度，先花三十秒算一下 n/N_C，反解 η，看它在哪个区间。这个动作成本极低，但能挡掉后面百分之八十的方向性错误。第二条是任何进入 10¹⁹ 量级的设计，都把统计模型、带隙变窄、俄歇复合三个开关作为默认配置，用开关前后对比确认影响量级，而不是凭经验判断"应该不大"。

第三条是关于测试的：简并区的数据一定要交叉验证，霍尔、电容电压、范德堡法至少两种对上，再拿去和模型比。单点数据在重掺杂区的不确定性比低掺杂大得多，拿一条曲线去拟合往往是在拟合测试误差。第四条是关于建模心态的：简并不是一个"高级话题"，它是器件做到一定程度必然撞上的墙。早一点把这套东西理顺，比后面在仿真里瞎调参数省时间得多。

最后一个我自己常用来做快速判断的口诀：看费米能级，别看掺杂数；看工作温度，别看室温表；看实际载流子，别看名义剂量。这三句话挡掉过我至少一半的返工。