功率器件的基本公式概念
目录
- 载流子输运与电阻
- 静电方程
- 载流子浓度与费米能级
- 电流密度方程
- 爱因斯坦关系式
- 连续性方程
- 电导率与电阻率
- 临界电场与击穿特性:碰撞电离系数
- 泊松方程
- 性能折衷与优值因子
- Baliga静态优值因子
- 藤平静态优值因子(Fujihira's Static Figure of Merit)
- 开关器件的总功率损耗
- Baliga动态优值因子
- Kim的动态品质因数
- Huang的动态品质因数
载流子输运与电阻
静电方程
著名的第一麦克斯韦方程组可表述为:
∇ ⋅ D = ρ ( x , y , z ) (1.1) \nabla \cdot D = \rho\left( x, y, z \right) \tag{1.1} ∇⋅D=ρ(x,y,z)(1.1)
其中D表示电通量密度,ρ为电荷密度。其简化形式即为高斯定律(又称泊松方程):
d 2 ψ d x 2 = − d E d x = − ρ ε S = q ( n − p + N A − N D ) ε S (1.2) \frac{d^2\psi}{dx^2} = -\frac{dE}{dx} = -\frac{\rho}{\varepsilon_S} = \frac{q(n - p + N_A - N_D)}{\varepsilon_S} \tag{1.2} dx2d2ψ=−dxdE=−εSρ=εSq(n−p+NA−ND)(1.2)
式中ψ、E、 ε S ε_S εS、q分别代表电势、电场强度、材料介电常数和单位电荷量。n、p、 N A N_A NA、 N D N_D ND分别对应自由电子浓度、自由空穴浓度、电离受主浓度和电离施主浓度。
载流子浓度与费米能级
在本征半导体中,自由载流子浓度n可通过电子态密度N(E)与费米-狄拉克分布函数F(E)的乘积积分求得:
n = ∫ E C m ∞ N ( E ) F ( E ) d E , (1.3) n = \int_{E_{C\,m}}^{\infty} N(E)F(E)\:dE,\tag{1.3} n=∫ECm∞N(E)F(E)dE,(1.3)
其中 E C m E_{Cm} ECm表示导带底能量最小值。根据费米-狄拉克统计,载流子占据概率F(E)遵循如下关系式:
F ( E ) = 1 1 + exp ( E − E F k T ) (1.4) F(E) = \frac{1}{1 + \exp\left( \frac{E - E_F}{kT} \right)} \tag{1.4} F(E)=1+exp(kTE−EF)1(1.4)
其中 E F E_F EF 为费米能级, k 为玻尔兹曼常数,T 为热力学温度。
将玻尔兹曼统计应用于式(1.4)并求解式(1.3),可得自由电子与空穴浓度方程:
n = N C exp ( − E C m − E F k T ) 或 E C m − E F = k T ln ( N C n ) (1.5) n = N_C \exp\left( -\frac{E_{C_m} - E_F}{kT} \right) \quad \text{或} \quad E_{C_m} - E_F = kT \ln\left( \frac{N_C}{n} \right) \tag{1.5} n=NCexp(−kTECm−EF)或ECm−EF=kTln(nNC)(1.5)
p = N V exp ( − E F − E V m k T ) 或 E F − E V m = k T ln ( N V p ) (1.6) p = N_V \exp\left( -\frac{E_F - E_{V_m}}{kT} \right) \quad \text{或} \quad E_F - E_{V_m} = kT \ln\left( \frac{N_V}{p} \right) \tag{1.6} p=NVexp(−kTEF−EVm)或EF−EVm=kTln(pNV)(1.6)
式中:
- N C N_C NC 和 N V N_V NV 分别为导带和价带的有效态密度(单位:cm − 3 ^{-3} −3)
- E C m E_{C_m} ECm 和 E V m E_{V_m} EVm 分别为导带底和价带顶的能量(单位:eV)
电流密度方程
电子电流密度 J n J_n Jn 和空穴电流密度 J p J_p Jp 由扩散电流与漂移电流叠加构成:
J n = q μ n n E + q D n d n d x (1.7) J_n = q \mu_n n E + q D_n \frac{dn}{dx} \tag{1.7} Jn=qμnnE+qDndxdn(1.7)
J p = q μ p p E − q D p d p d x (1.8) J_p = q \mu_p p E - q D_p \frac{dp}{dx} \tag{1.8} Jp=qμppE−qDpdxdp(1.8)
其中:
- μ n , μ p \mu_n, \mu_p μn,μp:电子与空穴的迁移率(单位:cm 2 ^2 2/(V·s))
- D n , D p D_n, D_p Dn,Dp:电子与空穴的扩散系数(单位:cm 2 ^2 2/s)
- E:电场强度(单位:V/cm)
- 负号表示空穴扩散电流方向与空穴漂移电流方向相反
总电流密度 J T J_T JT 为电子与空穴电流密度之和:
J T = J n + J p (1.9) J_T = J_n + J_p \tag{1.9} JT=Jn+Jp(1.9)
爱因斯坦关系式
对于掺杂非均匀的n型或p型半导体,在无外加电压时总电流为零。通过改写方程(1.7)和(1.8)
J n = q μ n n E + q D n d n d x = 0 (1.10) J_n = q\mu_n nE + qD_n\frac{dn}{dx} = 0 \tag{1.10} Jn=qμnnE+qDndxdn=0(1.10)
J p = q μ p p E − q D p d p d x = 0 (1.11) J_p = q\mu_p pE - qD_p\frac{dp}{dx} = 0 \tag{1.11} Jp=qμppE−qDpdxdp=0(1.11)
将式(1.5)和(1.6)代入电流密度方程(1.7)与(1.8),通过求解可得扩散系数与迁移率的定量关系:
D n = μ n k T q (1.12) D_n = \mu_n \frac{kT}{q} \tag{1.12} Dn=μnqkT(1.12)
D p = μ p k T q (1.13) D_p = \mu_p \frac{kT}{q} \tag{1.13} Dp=μpqkT(1.13)
此即著名的爱因斯坦关系式,仅适用于非简并半导体。在室温(T=300 K)下,热电压kT/q的值为0.0259 V,且电子与空穴的迁移率(mobility)与其扩散系数(diffusivities)高度相关。例如,若载流子迁移率较高,其扩散系数也必然较高。因此,对于需要电子和空穴具有高扩散系数的双极型器件(如双极晶体管、IGBT等),采用高迁移率的半导体材料更有利于实现良好的导电性。正因如此,爱因斯坦关系式(Einstein relationships)被广泛应用于分析双极导电机制。需注意的是,载流子的扩散系数还与载流子寿命(carrier’s life times) τ n \tau_{n} τn和 τ p \tau_{p} τp、扩散长度(diffusion lengths)Ln和Lp直接相关,
L n = D n τ n (1.14) L_{n}=\sqrt{D_{n}\tau_{n}} \tag{1.14} Ln=Dnτn(1.14)
L p = D p τ p (1.15) L_{p}=\sqrt{D_{p}\tau_{p}} \tag{1.15} Lp=Dpτp(1.15)
连续性方程
方程(1.10)和(1.11)给出的电流密度方程仅适用于稳态条件。连续性方程则用于描述半导体中载流子在瞬态行为下的特性,适用于存在低注入水平、载流子产生与复合的场景。
∂ n ∂ t = G n − U n + 1 q ∇ ⋅ J n , ∂ p ∂ t = G p − U p + 1 q ∇ ⋅ J p \begin{align*} \frac{\partial n}{\partial t} &= G_{n} - U_{n} + \frac{1}{q} \nabla \cdot J_{n} ,\tag{1.16}\\ \frac{\partial p}{\partial t} &= G_{p} - U_{p} + \frac{1}{q} \nabla \cdot J_{p} \tag{1.17} \end{align*} ∂t∂n∂t∂p=Gn−Un+q1∇⋅Jn,=Gp−Up+q1∇⋅Jp(1.16)(1.17)
其中 G n G_n Gn 和 G p G_p Gp 分别是电子和空穴的产生率, U n U_n Un 和 U p U_p Up 分别是电子和空穴的复合率。将方程 (1.10) 与 (1.12) 代入方程 (1.16) 与 (1.17) 后,方程变为:
∂ n p ∂ t = G n − n p − n p 0 τ n + n p μ n ∂ E ∂ x + μ n E ∂ n p ∂ x + D n ∂ 2 n p ∂ x 2 (1.18) \frac{\partial n_{p}}{\partial t} = G_{n} - \frac{n_{p} - n_{p0}}{\tau_{n}} + n_{p}\mu_{n}\frac{\partial E}{\partial x} + \mu_{n} E\frac{\partial n_{p}}{\partial x} + D_{n}\frac{\partial^{2}n_{p}}{\partial x^{2}} \tag{1.18} ∂t∂np=Gn−τnnp−np0+npμn∂x∂E+μnE∂x∂np+Dn∂x2∂2np(1.18)
∂ p n ∂ t = G n − p n − p n 0 τ p + p n μ p ∂ E ∂ x − μ p E ∂ p n ∂ x + D p ∂ 2 p n ∂ x 2 (1.19) \frac{\partial p_{n}}{\partial t} = G_{n} - \frac{p_{n} - p_{n0}}{\tau_{p}} + p_{n}\mu_{p}\frac{\partial E}{\partial x} - \mu_{p} E\frac{\partial p_{n}}{\partial x} + D_{p}\frac{\partial^{2}p_{n}}{\partial x^{2}} \tag{1.19} ∂t∂pn=Gn−τppn−pn0+pnμp∂x∂E−μpE∂x∂pn+Dp∂x2∂2pn(1.19)
其中 n p n_p np、 p n p_n pn、 n p 0 n_{p0} np0 和 p n 0 p_{n0} pn0 分别表示 p 型半导体中的电子浓度、n 型半导体中的空穴浓度、p 型半导体中的初始电子浓度以及 n 型半导体中的初始空穴浓度。
电导率与电阻率
在均匀掺杂的 n 型或 p 型半导体中,导电仅由漂移电流贡献。漂移电流与材料参数 q n μ n q n \mu_n qnμn 和 q p μ p q p \mu_p qpμp 以及外加电场 E E E 成正比,如方程 (1.10) 和 (1.11) 所示。因此, q n μ n q n \mu_n qnμn 和 q p μ p q p \mu_p qpμp 即成为该材料的电子电导率 σ n \sigma_n σn 和空穴电导率 σ p \sigma_p σp,
σ n = q n μ n σ p = q p μ p \begin{align} \sigma_{n} &= q n \mu_{n} \tag{1.20}\\ \sigma_{p} &= q p \mu_{p} \tag{1.21} \end{align} σnσp=qnμn=qpμp(1.20)(1.21)
电子和空穴的电阻率, ρ n \rho_n ρn 和 ρ p \rho_p ρp,分别是电导率的倒数形式。
ρ n = 1 q n μ n , ρ p = 1 q p μ p \begin{align} \rho_n &= \frac{1}{q\,n\,\mu_n}, \tag{1.22}\\ \rho_p &= \frac{1}{q\,p\,\mu_p} \tag{1.23} \end{align} ρnρp=qnμn1,=qpμp1(1.22)(1.23)
n 型和 p 型半导体的电阻 R n R_n Rn 和 R p R_p Rp 分别与导电面积 A A A 成反比,与长度 L L L 成正比,
R n = 1 q n μ n L A = ρ n L A , R p = 1 q p μ p L A = ρ p L A \begin{align} R_n &= \frac{1}{q\,n\,\mu_n}\,\frac{L}{A} = \rho_n\,\frac{L}{A}, \tag{1.24}\\ R_p &= \frac{1}{q\,p\,\mu_p}\,\frac{L}{A} = \rho_p\,\frac{L}{A} \tag{1.25} \end{align} RnRp=qnμn1AL=ρnAL,=qpμp1AL=ρpAL(1.24)(1.25)
由于电阻随导电面积的增大而减小,在功率电子领域,为了评估器件性能,需要考察电阻与面积的乘积。分别将面积与电阻相乘,即可得到 n 型和 p 型器件的比导通电阻,记作 R sp.n R_{\text{sp.n}} Rsp.n 和 R sp.p R_{\text{sp.p}} Rsp.p,
R s p , n = L q n μ n = ρ n L , R s p , p = L q p μ p = ρ p L \begin{align} R_{\mathrm{sp},n} &= \frac{L}{q n \mu_{n}} = \rho_{n} L,\tag{1.26} \\ R_{\mathrm{sp},p} &= \frac{L}{q p \mu_{p}} = \rho_{p} L \tag{1.27} \end{align} Rsp,nRsp,p=qnμnL=ρnL,=qpμpL=ρpL(1.26)(1.27)
临界电场与击穿特性:碰撞电离系数
当自由载流子在电场作用下加速运动时,会与晶格原子发生持续碰撞,引发晶格散射。若载流子在足够强的电场(通常高于该材料饱和速度对应的电场强度)中获得超额能量,其与原子碰撞时将通过碰撞电离过程产生电子-空穴对。这些新生载流子被强电场进一步加速,最终导致雪崩击穿。雪崩击穿条件取决于碰撞电离系数。硅(Si)与4H-SiC材料的碰撞电离系数可分别通过Fulop模型(αF)与Baliga模型(αB)近似表征。
α F ( S i ) = 1.8 × 10 − 35 E 7 (1.28) \alpha_{F}(\mathrm{Si}) = 1.8 \times 10^{-35} E^{7} \tag{1.28} αF(Si)=1.8×10−35E7(1.28)
α B ( 4 H - S i C ) = 3.9 × 10 − 42 E 7 (1.29) \alpha_{B}(4H\text{-}\mathrm{SiC}) = 3.9 \times 10^{-42} E^{7} \tag{1.29} αB(4H-SiC)=3.9×10−42E7(1.29)
泊松方程
半导体中触发雪崩击穿的电场强度值称为临界电场(EC)。该参数可通过碰撞电离系数结合一维P+/N结泊松方程求解获得。对于均匀掺杂的N型半导体,描述其电场分布的一维泊松方程见式(1.30),
d 2 V d x 2 = − d E d x = − q N D ε S (1.30) \frac{\mathrm{d}^{2} V}{\mathrm{d} x^{2}} = -\frac{\mathrm{d} E}{\mathrm{d} x} = -\frac{q N_{D}}{\varepsilon_{S}} \tag{1.30} dx2d2V=−dxdE=−εSqND(1.30)
其中,V 表示 P+/N 结两端维持的静电势,ND 为 N 型半导体的掺杂浓度。通过求解方程 (1.13),可导出电场强度与电势随距离变化的函数关系:
E ( x ) = − q N D ε S ( W D − x ) (1.31) E(x) = -\frac{qN_{D}}{\varepsilon_{S}}(W_{D} - x) \tag{1.31} E(x)=−εSqND(WD−x)(1.31)
V ( x ) = − q N D ε S ( W D x − x 2 2 ) (1.32) V(x) = -\frac{qN_{D}}{\varepsilon_{S}}\left(W_{D}x - \frac{x^{2}}{2}\right) \tag{1.32} V(x)=−εSqND(WDx−2x2)(1.32)
P⁺/N结中电场强度与静电势随距离变化的分布示意图
通过施加边界条件 V(WD)=VA,可解得:结界面处最大电场强度 E(0)=EM,耗尽区宽度 WD
W D = 2 ε S V a q N D (1.33) W_{D} = \sqrt{\frac{2 \varepsilon_{S} V_{a}}{q N_{D}}} \tag{1.33} WD=qND2εSVa(1.33)
E M = 2 q N D V a ε S (1.34) E_{M} = \sqrt{\frac{2 q N_{D} V_{a}}{\varepsilon_{S}}} \tag{1.34} EM=εS2qNDVa(1.34)
临界电场强度由碰撞电离系数的积分表达式确定,具体形式见方程(1.28)与(1.29):
∫ 0 W D α d x = 1 (1.35) \int_{0}^{W_{D}} \alpha dx = 1 \tag{1.35} ∫0WDαdx=1(1.35)
根据Sze与Baliga的研究,硅(Si)和碳化硅(SiC)的临界电场强度与掺杂浓度ND呈正比关系,该现象源于载流子与掺杂离子间增强的散射效应,
E C ( S i ) = 4010 N D 1 / 8 (1.36) E_C(\mathrm{Si}) = 4010 N_D^{1/8} \tag{1.36} EC(Si)=4010ND1/8(1.36)
E C ( 4 H - S i C ) = 33000 N D 1 / 8 (1.37) E_C(4H\text{-}\mathrm{SiC}) = 33000 N_D^{1/8} \tag{1.37} EC(4H-SiC)=33000ND1/8(1.37)
由方程(1.36)与(1.37)可导出硅(Si)和碳化硅(SiC)的击穿电压(VB) 随掺杂浓度(ND) 变化的函数关系:
V B ( S i ) = 5.34 × 10 13 N D − 3 / 4 (1.38) V_{B}(\mathrm{Si}) = 5.34 \times 10^{13} \, N_{D}^{-3/4} \tag{1.38} VB(Si)=5.34×1013ND−3/4(1.38)
V B ( 4 H - S i C ) = 3.0 × 10 15 N D − 3 / 4 (1.39) V_{B}(4H\text{-}\mathrm{SiC}) = 3.0 \times 10^{15} \, N_{D}^{-3/4} \tag{1.39} VB(4H-SiC)=3.0×1015ND−3/4(1.39)
性能折衷与优值因子
Baliga静态优值因子
联立比导通电阻方程(1.26)、耗尽区宽度方程(1.33)及临界电场方程(1.34),可将比导通电阻 Rsp 表示为材料本征参数的函数
R s p = 4 V B 2 ε S μ n E C 3 (1.40) R_{\mathrm{sp}} = \frac{4 V_{B}^{2}}{\varepsilon_{S} \mu_{n} E_{C}^{3}} \tag{1.40} Rsp=εSμnEC34VB2(1.40)
方程(1.40)的分母定义为Baliga优值因子(BFOM),该参数表征击穿电压平方(VB²)与比导通电阻(Rsp)的比值关系。
B F O M = ε S μ n E C 3 (1.41) \mathrm{BFOM} = \varepsilon_{S} \, \mu_{n} \, E_{C}^{3} \tag{1.41} BFOM=εSμnEC3(1.41)
需特别指出:Baliga优值因子(BFOM) 仅适用于标准平面型功率MOSFET,而无法用于评估超级结(Superjunction)等先进功率器件。
藤平静态优值因子(Fujihira’s Static Figure of Merit)
对于具有柱区长度L的超级结功率MOSFET,其理想比导通电阻可表述为:
R s p = 2 L q μ n N D (1.42) R_{\mathrm{sp}} = \frac{2L}{q \mu_{n} N_{D}} \tag{1.42} Rsp=qμnND2L(1.42)
方程(1.42)的推导基于对称超级结结构假设(即n型柱区与p型柱区宽度相等)。根据Fujihira的研究,超级结MOSFET的理想比导通电阻 Rsp 可表示为材料本征参数的函数:
R s p = 4 V B ε S μ n E C 2 d (1.43) R_{\mathrm{sp}} = \frac{4 V_{B}}{\varepsilon_{S} \mu_{n} E_{C}^{2}} d \tag{1.43} Rsp=εSμnEC24VBd(1.43)
其中 d 是超结金属氧化物半导体场效应晶体管(MOSFET)的元胞间距。该方程的问题是:超结 MOSFET 的比导通电阻可通过缩小元胞间距无限降低。方程 (1.43) 的分母被定义为藤平品质因数(Fujihira’s Figure of Merit, FFOM),
F F O M = ε S μ n E C 2 (1.44) \mathrm{FFOM} = \varepsilon_{S} \, \mu_{n} \, E_{C}^{2} \tag{1.44} FFOM=εSμnEC2(1.44)
开关器件的总功率损耗
如图所示,在硬开关应用中,器件总功率损耗包含四个分量:
- 导通损耗: D ⋅ P ON D \cdot P_{\text{ON}} D⋅PON
- 关断损耗: ( 1 − D ) ⋅ P OFF (1-D) \cdot P_{\text{OFF}} (1−D)⋅POFF
- 栅极驱动损耗: P GS P_{\text{GS}} PGS
- 开关损耗: P SW P_{\text{SW}} PSW
其中 D D D表示工作周期的占空比:
P T = D P O N + ( 1 − D ) P O F F + P G S + P S W (1.45) P_{T} = D P_{\mathrm{ON}} + (1 - D) P_{\mathrm{OFF}} + P_{\mathrm{GS}} + P_{\mathrm{SW}} \tag{1.45} PT=DPON+(1−D)POFF+PGS+PSW(1.45)
硬开关模式下功率器件的功率损耗
Baliga动态优值因子
功率MOSFET在开关应用中,导通电阻与开关损耗共同构成总功率损耗。1989年Baliga提出高频优值因子(BHFFOM)
B H F F O M = 1 R s p × C i s s . s p = f B (1.46) \mathrm{BHF}\;\mathrm{FOM} = \frac{1}{R_{\mathrm{sp}} \times C_{\mathrm{iss.sp}}} = f_{B} \tag{1.46} BHFFOM=Rsp×Ciss.sp1=fB(1.46)
其中:
- C i s s . s p C_{iss.sp} Ciss.sp:器件单位面积输入电容
- f B f_B fB:Baliga特征频率上限
需注意的是,式(1.46)定义的BHFFOM未计入输出电容相关损耗。
Kim的动态品质因数
为解决BHFFOM的局限性,Kim提出了一种新的品质因数NHFFOM,该指标在开关过程中考虑了输出电容 C o s s . s p C_{oss.sp} Coss.sp的影响。
N H F F O M = 1 R s p × C o s s . s p (1.47) \mathrm{NHF\ FOM} = \frac{1}{R_{sp} \times C_{oss.sp}} \tag{1.47} NHF FOM=Rsp×Coss.sp1(1.47)
方程(1.47)的局限性在于未纳入输入开关损耗 P G S P_{GS} PGS的影响。
Huang的动态品质因数
栅极-漏极电容( C G D C_{GD} CGD或 C r s s C_{rss} Crss)在充放电(关断与导通)过程中消耗的功率可间接视为开关功率损耗。尽管Huang提出的品质因数HDFOM考虑了其他电容分量(如 C G S C_{GS} CGS和 C D S C_{DS} CDS),但它为降低功率器件损耗提供了指导意义。
H D F O M = 1 R s p × C r s s , s p (1.48) \mathrm{HDFOM} = \frac{1}{R_{\mathrm{sp}} \times C_{\mathrm{rss},\mathrm{sp}}} \tag{1.48} HDFOM=Rsp×Crss,sp1(1.48)
然而,与Kim的品质因数类似,HDFOM也未纳入输入开关损耗 P G S P_{GS} PGS。目前仍缺乏一个全面的品质因数,能够通过方程(1.45)完整描述功率器件的总功率损耗。