跳到主要内容

半导体物理知识总结

固体物理

参考:

固体类型

  • 无定形:在原子/分子的尺度内有序
  • 单晶:晶粒,由晶界彼此分离,几何周期性好
  • 多晶:在多个原子/分子的尺度上有序,称为单晶区域,有不同的大小方向。

空间晶格

  • 晶格:晶体中原子的周期性排列
  • 晶胞/单胞:可以复制出整个晶体的一小部分晶体
  • 初基原胞:最小晶胞
  • 三维晶体中的每一个等效格点都可用矢量 r=pa+qb+sc\overline{r}=p\overline{a}+q\overline{b}+s\overline{c} 表示, pp , qq , ss 是正整数, a\overline{a} , b\overline{b} , c\overline{c} 是晶格常数。

晶体结构

rr 是原子半径

参数简单立方晶胞体心立方晶胞面心立方晶胞金刚石晶胞
晶胞边长2r2r43r\frac{4}{\sqrt{ 3 }}r22r2\sqrt{ 2 }r83r\frac{8}{\sqrt{ 3 }}r
配位数66881212 注意44
单个晶胞中的原子个数11224488
密排指数π6=52%\frac{\pi}{6}=52\%3π8=68%\frac{\sqrt{ 3 }\pi}{8}=68\%26π=74%\frac{\sqrt{ 2 }}{6}\pi=74\%3π16=34%\frac{\sqrt{ 3 }\pi}{16}=34\%(注意并不是最高)
六方最密堆积的配位数也是 12 ,密排指数同样为 74%74\%

密勒指数

晶向和晶面都有密勒指数。

  • 晶向:描述晶体内部的一个方向:分解矢量。
    • 晶向:[hkl][hkl],如果为负数要写为 h\overline{h}
    • 等价晶向:<hkl><hkl>
  • 晶面:描述晶体内部的一个平面:截距→倒数→取最小公倍数
    • 晶面:(hkl)(hkl),负数同上
    • 等价晶面:{hkl}\{hkl\}
  • 特性
    • 互质
    • 都是整数
    • 立方晶体:晶面和晶向的垂直面密勒指数相同
    • 其他晶体:若晶面和晶向密勒指数相同,晶向/晶面满足矢量分解/合成的关系。
    • 对称性
      • 是否能把晶向/晶面转化为等价晶面/晶向与空间晶格的对称性有关
      • 简单立方晶体:由于对称性,八种晶向可以等价为 <111><111>

其他

  • 1A˚=108cm=1010m1\mathring{A}=10^{-8}cm=10^{-10}m
  • 原子面密度
面密度=每晶面的原子数每晶面的面积面密度=\frac{每晶面的原子数}{每晶面的面积}

量子物理

高数知识Recall

  • 对于形如 y+yp+qy=0y''+y'p+qy=0 的微分方程,直接解特征根方程 r2+rp+q=0r^2+rp+q=0 ,解得 y=C1er1x+C2er2xy=C_{1}e^{ r_{1}x }+C_{2}e^{ r_{2}x }

相关参数

  • hh 普朗克常数,h=6.625×1034Jsh=6.625\times10^{-34}J\cdot s
  • hwˉ\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}} 约化普朗克常数, hwˉ=h2π=1.054×1034Js\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}=\frac{h}{2\pi}=1.054\times 10^{-34}J\cdot s
  • m0=9.1×1031kgm_{0}=9.1\times 10^{-31}kg ,电子质量
  • q=1.6×1019Cq=1.6\times 10^{-19}C ,电子电荷
  • ϵ0=8.85×1012F/m\epsilon_{0}=8.85\times 10^{-12} F/m ,真空介电常数
  • 1eV=1.6×1019J1eV=1.6\times 10^{-19} J ,一个电子伏特

基本原理

能量量子化

  • 量子能量(其中 ν\nu 是频率)
E=hνE=h\nu
  • 光电效应中光电子最大动能
T=12mv2=hνhν0T=\frac{1}{2}mv^2=h\nu-h\nu_{0}

其中 hν0h\nu_{0} 为逸出功。

  • 光电效应
    • 高频率的光→更高能量的电子
    • 光强大的光→发射更多电子

波粒二象性与物质波

  • 粒子动量/德布罗意公式/普朗克-爱因斯坦方程式
p=hλ=hwˉkp=\frac{h}{\lambda}=\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}k

(其中:k=2πλk=\frac{2\pi}{\lambda}

  • 粒子能量
E=hν=hwˉwE=h\nu=\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}w

(其中:w=2πνw=2\pi \nu

不确定原理

  • 动量不确定度
ΔpΔxhwˉ2\Delta p\Delta x\geq \frac{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}}{2}
  • 能量不确定度
ΔEΔthwˉ2\Delta E\Delta t\geq \frac{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}}{2}

薛定谔方程

  • 一维非相对论薛定谔波动方程
22m2Ψ(x,t)x2+V(x)Ψ(x,t)=jΨ(x,t)t{\frac{-\hbar^{2}} {2 m}} \cdot{\frac{\partial^{2} \Psi( x, \, t )} {\partial x^{2}}}+V ( x ) \Psi( x, \, t )=\mathrm{j} \hbar{\frac{\partial\Psi( x, \, t )} {\partial t}}

其中 Ψ(x,t)\Psi(x,t) 是一个复值函数,是孤立系统量子态的数学表示。

  • 令波函数 Ψ(x,t)=ψ(x)ϕ(t)\Psi(x,t)=\psi(x)\phi(t) ,可以将薛定谔波动方程化为一边为 xx 的函数,一边为 tt 的函数的形式,并通过设出 EE 这个参量求得 ϕ(t)=ejwt\phi(t)=e^{ -jwt } 。则薛定谔波动方程可以被简化为:
2ψ(x)x2+2m2(EV(x))ψ(x)=0\frac{\partial^{2} \psi( x )} {\partial x^{2}}+\frac{2 m} {\hbar^{2}} \, ( E-V ( x ) ) \psi( x ) \,=\, 0
  • 由于平方可以消去函数中的复数值,故而将其平方,可以得到与时间无关的 Ψ(x,t)2\mid \Psi(x,t)\mid^2 ,也即概率密度函数, 其中 Ψ(x,t)2dx\mid \Psi(x,t)\mid^2dx 表示在 tt 时刻, xxdxdx 位置找到粒子的概率。
  • 波函数归一化条件
Ψ(x)2dx=1\int ^\infty_{-\infty}\mid \Psi(x)\mid^2 \, dx=1
  • 概率密度的条件
    • 根据概率密度的定义, ψ(x)\psi(x) ==有限==、==单值==
    • 如果粒子总能量 EE 和势函数 V(x)V(x) 在任何位置==均为有限值==,则二阶函数必须有限,做那一点的泰勒展开,可以得出==一阶导数连续==。由于一阶导数跟粒子动量有关,也必须==有限单值==,从而得出函数==连续==
    • 如果势函数 V(x)V(x) 在某处是无限值,由于上面的推导,可以知道一阶导数不一定连续

薛定谔方程的应用

自由空间的电子

  • 条件 V(x)=0V(x)=0
  • 解微分方程2ψ(x)x2+2mE2ψ(x)=0\frac{\partial^{2} {\psi} ( x )} {\partial x^{2}}+\frac{2 m E} {\hbar^{2}} \, {\psi} ( x )=0 ,得 ψ(x)=Aejkx+Bejkx\psi(x)=Ae^{ jkx }+Be^{ -jkx } 。我们使用这种形式的解(行波解)是因为它在自由空间中。 ^977c93
  • 其中:k=2mEhwˉ2k=\sqrt{ \frac{2mE}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2} } 称为波数, 由于 2mE=p22mE=p^2 ,于是有 k=phwˉk=\frac{p}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}} 。此处可以联动前面 [[#波粒二象性与物质波) 的内容。
  • 代入时间方程:Ψ(x,t)=Aej(kxwt)+Bej(kx+wt)\Psi(x,t)=Ae^{ j(kx-wt) }+Be^{ -j(kx+wt) } ,表现为一个方向为 +x+x 的波和一个方向为 x-x 的波的和。
  • 其中, AA / BB 的值可以通过边界条件确定。
  • 如果有一个沿 +x+x 方向运动的电子,则 B=0B=0 ,此时概率密度函数是一个常量,物理意义是:有明确动量定义的自由粒子在空间任意位置出现的概率相等。

无限深势阱

  • 由薛定谔方程可知,如果 EE 有限,在 VV 为无限的区域, ψ(x)\psi(x) 必须为 00
  • VV00 的区域,同自由空间的电子,可以解微分方程。然而此时粒子处于束缚态,所以解的形式需要写成驻波解: ψ(x)=A1coskx+A2sinkx\psi(x)=A_{1}\cos kx+A_{2}\sin kx ,其中 A1A_{1} / A2A_{2} 为复数。 ^916c31
  • 考虑薛定谔方程概率密度条件:
    • 在边界( x=0x=0x=ax=a )位置波函数必须连续:ψ(x=0)=ψ(x=a)=0\psi(x=0)=\psi(x=a)=0
      • A1=0A_{1}=0
      • ka=nπka=n\pi 时成立,其中参数 nn 为正整数(负数不考虑,在求概率密度的时候,实际与正整数相同)
  • 考虑归一化条件 0aA22sin2kxdx=1\int ^a_{0}A_{2}^2\sin^2kx \, dx=1 ,解得 A2=2aA_{2}=\sqrt{ \frac{2}{a} }
  • 得到:ψ(x)=2asin(nπxa)\psi(x)=\sqrt{ \frac{2}{a} }\sin\left( \frac{n\pi x}{a} \right) ,其中 nn 为正整数。
  • 由于 kk 的值需满足 k=2mEhwˉ2k=\sqrt{ \frac{2mE}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2} } 以及 k=nπak=\frac{n\pi}{a} 两个条件,联立可得粒子的总能量
  • 无限深势阱中粒子的波函数为
ψ(x)=2asinknx\psi(x)=\sqrt{ \frac{2}{a} }\sin k_{n}x

由于 knk_{n} 必须是分立值,相应的总能量也只能是分立值。这里就体现了量子力学中能量的量子化。

阶跃势函数

  • 粒子流入射至阶跃势垒,流动方向为 +x+x ,起始点 -\infty 。粒子的总能量小于势垒的高度 E<V0E<V_{0}
  • V=0V=0 的区域,方程的通解为 ψ1(x)=A1ejk1x+B1ejk1x\psi_{1}(x)=A_{1}e^{ jk_{1}x }+B_{1}e^{ -jk_{1}x }
  • 在势函数 V=V0V=V_{0} 的区域,假设 E<V0E<V_{0} ,则有 ψ2(x)=A2ek2x+B2e+k2x\psi_{2}(x)=A_{2}e^{ -k_{2}x }+B_{2}e^{ +k_{2}x }注意这里e指数上是没有虚数标记的,原因是 E<V0E<V_{0}) 其中 k2=2m(V0E)hwˉ2k_{2}=\sqrt{ \frac{2m(V_{0}-E)}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2} }
  • 边界条件:
    • ψ2(x)\psi_{2}(x) 必须保持有限, B2=0B_{2}=0
    • x=0x=0 处保持连续, ψ1(0)=ψ2(0)\psi_{1}(0)=\psi_{2}(0)
    • 势函数有限——波函数的一阶导 ψ1xx=0=ψ2xx=0\left. \frac{\partial\psi_{1}} {\partial x} \right|_{x=0}=\left. \frac{\partial\psi_{2}} {\partial x} \right|_{x=0} 连续
  • 解得系数:
B1=(k22+2jk1k2k12)k22+k12A1B_{1}=\frac{-( k_{2}^{2}+2 \mathrm{j} k_{1} k_{2}-k_{1}^{2} )} {k_{2}^{2}+k_{1}^{2}} \cdot A_{1} A2=2k1(k1jk2 )k22+k12A1A_{2}=\frac{2 k_{1} ( k_{1}-\mathrm{ j k_{2} ~} )} {k_{2}^{2}+k_{1}^{2}} \cdot A_{1}
  • 驻波的概率密度函数同样是取共轭相乘。例如,入射粒子流的流量 viA1A1v_{i} \cdot A_{1}\cdot A_{1}^* ,反射粒子流的流量 vrB1B1v_{r}\cdot B_{1}\cdot B_{1}^*
  • k1=2mEhwˉ2k_{1}=\sqrt{ \frac{2mE}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2} } 以及 E=12mv2E=\frac{1}{2}mv^2 可得(此时由于 V(x)=0V(x)=0 粒子的总能量与动能相等), vi=hwˉ2mk1v_{i}=\frac{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2}{m}\cdot k_{1} ,也可以知道 vr=hwˉ2mk1v_{r}=\frac{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2}{m}\cdot k_{1} (我觉得是从第二项得到的,方法同上,因为二者的 kk 都是 k1k_{1}
  • 由此可以得到反射系数 R=νrB1B1ν1A1A1R=\frac{\nu_{r} \cdot B_{1} \cdot B_{1}^{*}} {\nu_{1} \cdot A_{1} \cdot A_{1}^{*}}1.01.0
  • 存在有限概率的粒子穿越势垒,但最后粒子会完全返回到势垒为 00 的区域中。

矩形势垒

  • 同样是 E<V0E<V_{0} 的情况,则三个薛定谔方程的解为
    • ψ1(x)=A1ejk1x+B1ejk1x\psi_{1}(x)=A_{1}e^{ jk_{1}x }+B_{1}e^{ -jk_{1}x }
    • ψ2(x)=A2ek2x+B2ek2x\psi_{2}(x)=A_{2}e^{ k_{2}x }+B_{2}e^{ -k_{2}x }
    • ψ3(x)=A3ejk1x+B3ejk1x\psi_{3}(x)=A_{3}e^{ jk_{1}x }+B_{3}e^{ -jk_{1}x }
    • k1=2mEhwˉ2k_{1}=\sqrt{ \frac{2mE}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2} }
    • k2=2mhwˉ2(V0E)k_{2}=\sqrt{ \frac{2m}{\mathrm{h}\kern{-7pt}\bar{\small\phantom{w}}^2}(V_{0}-E) }
  • 边界条件:
    • 连续:x=0x=0x=ax=a 处连续
    • 一阶导数连续
    • B3=0B_{3}=0 因为在那个区域没有反射的粒子
  • 透射系数:
T=νtA3A3νiA1Ai=A3A3A1A1T=\frac{\nu_{t} \cdot A_{3} \cdot A_{3}^{*}} {\nu_{i} \cdot A_{1} \cdot A_{i}^{*}}=\frac{A_{3} \cdot A_{3}^{*}} {A_{1} \cdot A_{1}^{*}}
  • EV0E\ll V_{0} 时,T16(EV0)(1EV0)e2k2aT\approx 16\left( \frac{E}{V_{0}} \right)\left( 1-\frac{E}{V_{0}} \right)e^{ -2k_{2}a }
  • 存在有限概率穿过势垒到达对面的区域。

单电子原子

  • 三维的与时间无关的薛定谔波动方程
2ψ(r,θ,ϕ)+2m02(EV(r))ψ(r,θ,ϕ)=0\nabla^{2} \psi( r, \theta, \phi)+{\frac{2 m_{0}} {\hbar^{2}}} ( E-V ( r ) ) \psi( r, \theta, \phi)=0
  • 球坐标系的
1r2r(r2ψr)+1r2sin2 ⁣θ2ψϕ2+1r2sin ⁣θθ(sinθ ⁣ ⁣ψθ)+2m02(EV(r))ψ=0{\frac{1} {r^{2}}} \cdot{\frac\partial {\partial r}} \left( r^{2} \, {\frac{\partial\psi} {\partial r}} \right)+{\frac{1} {r^{2} \sin^{2} \! \theta}} \cdot{\frac{\partial^{2} \psi} {\partial\phi^{2}}}+{\frac{1} {r^{2} \sin\! \theta}} \cdot{\frac\partial {\partial\theta}} \left( \sin\theta\! \cdot\! {\frac{\partial\psi} {\partial\theta}} \right) + \, \frac{2 m_{0}} {\hbar^{2}} \, ( E-V ( r ) ) \psi=0
  • 就像前面解一维形式的薛定谔方程一样:ψ(r,θ,ϕ)=R(r)Θ(θ)Φ(ϕ)\psi(r,\theta ,\phi)=R(r)\cdot\Theta(\theta)\cdot \Phi(\phi) ,同样可以看到第二项仅是 ϕ\phi 的函数,同样令第二项为 m2-m^2 ,可以解得 Φ=ejmϕ\Phi=e^{ jm\phi } 。我们知道,在球坐标系中, ϕ\phi 是表示极角的量。ϕ\phi 可以取任意大于 2π2\pi 的值,但对于一个球坐标系,相同的角度只能有一个函数值。因而为了保证波函数的单值性, mm 需要是整数。我们把 mm 称为==量子数==。
  • 势函数
V(r)=e24πϵ0rV(r)=-\frac{e^2}{4\pi\epsilon_{0}r}
  • 继续解薛定谔方程,可以得出:
    • 主量子数:n=1,2,3n=1,2,3…
    • 角量子数:l=n1,n2,,0l=n-1,n-2,…,0
    • 磁量子数:m=l,l1,,0\mid m\mid=l,l-1,…,0
  • 电子的能量(其中 nn 是能量量子数/轨道数) En=m0e4(4πϵ0)222n2=13.6n2eVE_{n}=\frac{-m_{0} \mathrm{e}^{4}} {( 4 \pi\epsilon_{0} )^{2} 2 \hbar^{2} n^{2}}=-\frac{13.6}{n^2}eV 来源于对方程
m0v2r=14πϵ0q2r2\frac{m_{0}v^2}{r}=\frac{1}{4\pi\epsilon_{0}} \frac{q^2}{r^2}

以及方程

L=n=m0vrL=n\hbar=m_{0}vr

的联立,解得

r=4πϵ0n22q2mr=\frac{4\pi\epsilon_{0}n^2\hbar^2}{q^2m}

代入 E=Ek+Ep=1214πϵ0q2rE=E_{k}+E_{p}=-\frac{1}{2} \frac{1}{4\pi\epsilon_{0}} \frac{q^2}{r} 之后解出电子能量。 其中, nn 为主量子数。这体现了能量的量子化。

  • 最小的能量状态,波函数表示为
ψ100=1π(1a0)3/2 er/a0\psi_{1 0 0}=\frac{1} {\sqrt{\pi}} \cdot\left( \frac{1} {a_{0}} \right)^{3 / 2} \mathrm{~ e}^{-r / a_{0}}

其中 a0a_{0}

a0=4πϵ02m0e2=0.529A˚a_{0}=\frac{4 \pi\epsilon_{0} \, \hbar^{2}} {m_{0} e^{2}}=0. 5 2 9\mathring{A}

将该函数取共轭后与自己相乘,可以得到径向概率密度函数,最大概率出现在半径 r=a0r=a_{0} 处。

元素周期表

  • 增加了电子的本征角动量(自旋)值为 s=±12s=\pm \frac{1}{2}

载流子模型与特性

符号

  • EcE_{c} 导带底
  • EvE_{v} 价带顶
  • EGE_{G} 带隙/禁带宽度
  • nn 本征半导体内部的载流子数
  • pp 本征半导体内部的空穴数
  • EDE_{D} 施主态能级(虚线)
  • EAE_{A} 受主态能级(虚线)
  • EFE_{F} 费米能级/费米能量
  • k=8.617×105eV/Kk=8.617\times 10^{-5} eV /K 玻尔兹曼常数
  • n0n_{0} 导带中的电子浓度值
  • p0p_{0} 价带中的电子浓度值
  • nin_{i} 本征半导体中的电子浓度
  • pip_{i} 本征半导体中的空穴浓度
  • EFiE_{Fi}EiE_{i} 本征半导体的费米能级(本征费米能级)
  • kT=0.0259eVkT=0.0259eV

载流子模型

价键模型-与空间相关

  • 类似于高中化学的模型,硅原子在金刚石晶体中,每个原子与四个最近邻原子相互吸引,形成共价结合。每个原子贡献出四个共享的电子,也接受四个从其他原子共享的电子。
  • 一种理想模型
    • 点缺陷/缺少原子
    • 价键断裂/释放出自由电子
  • 在温度接近 0K0K ,半导体中没有缺陷和杂质原子时,价键模型对半导体的描述是等效的。

能带模型-与能量有关

思想脉络:独立原子→晶体

  • 当两个原子距离变近时,原子中电子的势函数相互交叠,于是两个原子互相影响。这种影响使离散的量子化能级分裂为两个分离的能级。下图是一个示意图,参数 r0r_{0} 代表晶体中平衡原子间的距离。

  • 分离的能级中,能量仍然是离散的。根据泡利不相容原理,量子数不能改变,故而每个电子须占据独立不同的电子态。由于晶体中的原子所带的电子很多,且允带宽度不大,故而每个能级的间距是很小的,可以认为允带处于准连续能量分布。
  • 就硅原子而言,我们在讨论硅的能带时只需考虑 n=3n=3 能级上的价电子,因为剩下两个电子壳层受到核的束缚。 n=3n=3 能级上有俩能量状态, 3s3s 包括俩量子态, 3p3p 包含六个量子态。
  • 随着原子间距的减小,每个原子的其中 44 个量子态处于较低能带,另外 44 个处于较高能带。 0K0K 时,电子都处于最低能量状态,较低能带(价带)的所有状态都满,较高能带(导带)的所有状态都空,如下图所示。

  • 我们令 EcE_{c} 为最低可能的导带能量, EvE_{v} 为最高可能的价带能量, EG=EcEvE_{G}=E_{c}-E_{v} ,称为带隙的能量。
  • 通过克勒尼希-彭尼模型(Neamen那本书3.1.2)解薛定谔方程可以得到粒子能量 EE 和对应波数 kk 的关系如图所示。

2π2\pi 为周期进行平移,得到

  • 基本能带模型图,一条线表示导带中最低的能量,另一条表示价带顶的能量, yy 是电子的能量轴, xx 方向是位置轴。(如图所示)

能带模型的三维扩展

  • 把一维的 EkE-k 图扩展到三维,kk 坐标的变化不仅代表 kk 空间的变化,还代表晶向的变化(由于在特定晶向上, EkE-k 关系是对称的,故而我们可以在负半轴上画更多的东西),因此,从三维空间的 EkE-k 关系图上,我们可以看到不同晶向的 EkE-k 关系。
  • 直接带隙半导体:导带最小能量与价带最大能量具有相同的 kk 坐标,两个允带之间电子的跃迁不会对动量产生影响。例子: GaAsGaAs ,如图所示为砷化镓的能带结构:
  • 间接带隙半导体:价带能量最大值和导带能量最小值的 kk 坐标不同的半导体,电子在价带和导带中跃迁时,需要使用动量守恒定律(包含与晶体的相互作用),通过吸收/释放声子的方式对其进行补偿。例子: SiGeSi、Ge ,如图所示是硅的能带结构:

载流子-电流传输的性质

  • 电子与空穴倾向于向能级低的地方运动。
  • 不存在载流子/电流的情况
    • 价键模型没有断裂的键
    • 能带模型中价带被电子填满而导带是空的(电子的动量也是量子化的,每个可能的动量状态有一个大小相等、方向相反的动量状态),因为被填满了,所以电子的净动量在能带中恒等于 00
  • 存在载流子的情况
    • 价键模型:出现键的断裂,价电子成为自由电子/载流子
    • 能带模型:价带电子受到激发进入导带,产生载流子
  • 空穴
    • 共价键断裂时,共价键中产生缺键/空位,电子整体的移动造成空位的移动,称为空穴。

带隙与材料分类

材料性质能带情况电阻率
绝缘体热能只能将很少的电子从价带激发到导带,很难用加热的方法产生电子和空穴,不良导体能带全满或全空,带隙能量 EgE_{g} 通常为 3.56eV3.5\sim 6eV ,带隙很宽ρ=10141018Ωcm\rho=10^{14}\sim 10^{18}\Omega\cdot cm
金属良导体,其中有大量的载流子带隙很小,或者是相互重合没有带隙,导带和价带相互重合ρ=105106Ωcm\rho=10^{-5}\sim 10^{-6}\Omega\cdot cm ,电阻率随温度升高而升高
半导体室温下,价带中的一些电子受到激发可以得到足够的热能,到达导带,产生一定数量的载流子带隙介于上述两者之间。近似空带或近似满带,外加电场可使得电荷运动并产生电流, EgE_{g} 通常为 1eV1eVρ=103103Ωcm\rho=10^{-3}\sim 10^3\Omega \cdot cm ,电阻率随温度降低而升高

载流子特性

  • 电荷粒子:导电电子与空穴;导电电子带负电荷,空穴带正电荷

有效质量

  • 电子在晶格中的运动与在自由空间中不同,既受外力影响也受内力影响。因为 很难一一考虑粒子所受的内力,即
Fext=maF_{ext}=m^*a

mm^* 称为有效质量,概括了粒子的质量以及内力的作用。

  • 如何求出粒子的有效质量:
    • E=p22m=2k22mE=\frac{p^2}{2m}=\frac{\hbar^2k^2}{2m} 开始,之后对 kk 求导,得到 dEdk=2km=pm=v{\frac{\mathrm{d} E} {\mathrm{d} k}}={\frac{\hbar^{2} k} {m}}={\frac{\hbar p} {m}}=\hbar v ,可知 EEkk 的一阶导数与粒子速度有关。
    • kk 求二阶导数:得到 12d2Edk2=1m\frac{1} {\hbar^{2}} \, \frac{\mathrm{d}^{2} E} {\mathrm{d} k^{2}}=\frac{1} {m} ,我们可以通过这个方程计算有效质量或者判断有效质量的正负。
  • 利用有效质量解题
    • 对电子而言,运动方程为 F=qE=mndvdtF=-qE=m_{n}^* \frac{dv}{dt}
    • 对空穴而言,运动方程为 F=qE=mpdvdtF=qE=m_{p}^* \frac{dv}{dt}
    • 注意电子和空穴之间的区别

载流子数

  • 本征半导体:一种半导体材料,杂质含量小于热激发的电子数和空穴数,或者本征半导体是没有添加其他杂质的天然材料
  • 热平衡条件下,半导体内部的本征半导体
n=p=nin=p=n_{i}

其中硅的 nin_{i}1×1010/cm31\times 10^{10}/cm^3

载流子数的控制(掺杂)

  • 施主元素
    • 元素周期表中的 VV 族元素,多出一个施主电子,不能进入共价键,只能在施主电子周围形成较弱束缚,容易挣脱晶格的束缚,在晶格中自由运动称为载流子。
    • 我们可以把施主离子看作氢原子的原子核,而把施主电子看作氢原子的电子,于是可以得到施主电子的结合能 EBmnm01KS2EHn=10.1eVE_{B}\simeq \frac{m_{n}^*}{m_{0}} \frac{1}{K_{S}^2}E_{H\mid n=1} \simeq -0.1eV ,其中 mnm^*_{n} 是电子的有效质量, KS=11.8K_{S}=11.8 是硅的介电常数。
    • 束缚电子占据的电子能级 ED=EcEBE_{D}=E_{c}-\mid E_{B}\midEB\mid E_{B}\mid 是束缚电子的结合能),能带图如下图 注意这里的能级表现为一组非连续的线,因为束缚电子和施主之间的距离是局域化的,也就是图中 Δx\Delta x 的范围
    • 带正电的施主原子不移动
    • 不增加空穴的浓度
  • 受主元素
    • IIIII I 族元素,缺少一个电子,需要从临近的硅原子中夺取一个价电子,于是产生一个空穴,称为载流子
    • 可以类比施主元素得到受主电子的结合能(使用有效质量与介电常数)与受主能级 EAE_{A}
    • 带负电的受主离子不移动
    • 不增加电子的浓度
  • 施主能级与受主能级在不同温度下的填充效果可以见下图

  • 非元素半导体的掺杂:类似前述掺杂,施主/受主元素的确定来源于不同元素的施主/受主的状态。
  • 如果杂质既有施主的作用又有受主的作用,则称为两性杂质。

相关术语

  • 掺杂:为了增加半导体内电子或空穴的浓度,将一定数量的特殊杂质原子掺入到半导体的体内。
  • 本征半导体:没有掺杂的半导体;非常纯净的半导体(样品内的杂质原子数量可以忽略不计);具有材料固有性质的半导体。
  • 非本征半导体:掺杂半导体;加入的杂质原子控制半导体性质的半导体。
  • 施主:能增加电子浓度的杂质原子;n型掺杂。
  • 受主:能增加空穴浓度的杂质原子;p型掺杂
  • n型材料:掺有施主的材料;半导体内的电子浓度大于空穴浓度。
  • p型材料:掺有受主的材料;半导体内的空穴浓度大于电子浓度。

热平衡态下的载流子浓度

  • 掺杂半导体内载流子浓度的精确数值
  • 掺杂半导体内载流子的分布(是能量的函数)

态密度

  • 研究导带/价带中任意给定的能级有多少种状态,单位为 (Numbers of States/cm3)/eV(Numbers\ of\ States/cm^3) /eV
  • 当电子能量为 EE 且距离带边不远时:
gc(E) ⁣= ⁣8π2h3(mn) ⁣3/2(EEc)    for  E ⁣ ⁣Ecg_{c} ( E ) \!=\! \frac{8 \pi\sqrt{2}} {h^{3}} \left( m_{n}^{*} \right)^{\! 3 / 2} \sqrt{\left( E-E_{c} \right)} \; \; \mathrm{f o r} \; E \! \geq\! E_{c} gv(E) ⁣= ⁣8π2h3(mp) ⁣3/2(Ev ⁣ ⁣E)    for  E ⁣ ⁣Evg_{v} ( E ) \!=\! \frac{8 \pi\sqrt{2}} {h^{3}} \left( m_{p}^{*} \right)^{\! 3 / 2} \sqrt{\left( E_{v} \!-\! E \right)} \; \; \mathrm{f o r} \; E \! \le\! E_{v}

or

gc(E)=mn2mn(EEc)π23  ,EEcgv(E)=mp2mp(EvE)π23  ,EEv\begin{array} {c c} {{g_{\mathrm{c}} ( E ) \,=\, \displaystyle{\frac{m_{\mathrm{n}}^{\star} \sqrt{2 m_{\mathrm{n}}^{\star} ( E \,-\, E_{\mathrm{c}} )}} {\pi^{2} \hbar^{3}}} \;, \qquad E \, \geqslant\, E_{\mathrm{c}}}} \\ {{g_{\mathrm{v}} ( E ) \,=\, \displaystyle{\frac{m_{\mathrm{p}}^{\star} \sqrt{2 m_{\mathrm{p}}^{\star} ( E_{\mathrm{v}} \,-\, E )}} {\pi^{2} \hbar^{3}}} \;, \qquad E \, \leqslant\, E_{\mathrm{v}}}} \end{array}
  • 物理意义举例: gc(E)dEg_{c}(E)dE 描述能量 EEE+dEE+dE 之间每立方厘米内导带的状态数(条件 EEcE\geq E_{c}
  • 相关图如下图所示

费米分布函数

  • 在==热平衡==条件下,求得能量为 EE 的状态被电子占据的几率
  • 费米分布函数 ( TT 的单位是 KK
f(E)=11 + e(EEF)/kTf ( E ) \,=\, {\frac{1} {1 ~+~ e^{( E-E_{F} ) / k T}}}
  • 性质 ^af2f6b
    • T0KT\to0K
      • E<EFE<E_{F}f(E)1f(E)\to 1 ,费米能级之下的量子态全被电子占据
      • E>EFE>E_{F}f(E)0f(E)\to{0} ,费米能级之上的量子态几乎全空
      • 费米能级是绝对零度时,粒子所能达到的最高能级
    • E=EFE=E_{F}f(EF)=12f(E_{F})=\frac{1}{2}
    • EEF+3kTE\geq E_{F}+3kTf(E)e(EEF)/kTf(E)\simeq e^{ -(E-E_{F})/kT } ,在费米能级之上 3kT3kT 或者更多的量子态是空的
    • EEF3kTE\leq E_{F}-3kTe(EEF)/kT1e^{ (E-E_{F})/kT }\ll 1 ,因此在费米能级低 3kT3kT 或者更多的量子态几乎是满的。
    • 示意图:温度等于/大于 0K0K 时费米分布函数随能量变化的特性
    • 1fF(E)=11+e(EFE)/kT1-f_{F}(E)=\frac{1}{1+e^{ (E_{F}-E )/kT}} ,与 fF(E)f_{F}(E) 关于费米能级对称,也就是 EFE_{F} 以上 dEdE 距离处与 EFE_{F} 以下 dEdE 距离处空状态的概率相等。
    • 联系前面有关施主能级与受主能级的知识,我们可以知道: T=0KT=0K 的时候,对于 nn 型半导体,费米能级必定高于施主能级;对于 pp 型半导体,费米能级必定低于受主能级。

平衡载流子分布

  • 不同能带中载流子的分布: 态密度×占有因子态密度\times 占有因子
    • n=EcEtopgc(E)f(E)dEn=\int_{E_{c}}^{E_{t o p}} g_{c} ( E ) f ( E ) d E
    • 导/价带:导/价带的态密度;
    • 电子/空穴:选取相应的占有因子( fF(E)f_{F}(E) 或者 1fF(E)1-f_{F}(E)
  • 当费米能级处于不同的位置时,载流子的分布情况也不同(注意,下图中的横轴分别为态密度和费米分布函数。
  • ==导带电子==的热平衡浓度由上一节的相乘后积分(从导带往上积)以及一堆复杂的数学推导得出(注意: 推导的前提条件是:(EcEF)3kT(E_{c}-E_{F})\gg 3kT ,此时费米概率分布函数简化为玻尔兹曼近似,也即 f(E)e(EEF)/kTf ( E ) \, \cong\, e^{-( E-E_{F} ) \, / \, k T}
n0=Ncexp[(EcEF)kT]n_{0}=N_{c} \exp\left[ {\frac{-( E_{c}-E_{F} )} {k T}} \right]

其中: Nc=2(2πmnkTπh2)3/2N_{c}=2 \Bigl( {\frac{2 \pi m_{n}^{*} \, k T} {\pi\, h^{2}}} \Bigr)^{3 / 2} ,是导带的状态密度

  • 价带空穴的热平衡浓度,从价带往下积,经过数学推导得出(注意: 推导的前提条件是:(EFEv)3kT(E_{F}-E_{v})\gg 3kT ,此时费米概率分布函数简化为玻尔兹曼近似
p0=Nvexp[(EFEv)kT].p_{0}=N_{v} \exp\left[ {\frac{-( E_{F}-E_{v} )} {k T}} \right].

其中: Nv=2(2πmpkTh2)3/2N_{v}=2 \left( \frac{2 \pi m_{p}^{*} k T} {h^{2}} \right)^{3 / 2} ,是价带的状态密度

  • 对比上述两式,有所变化的是不同载流子的有效质量以及导带/价带的顺序,其他没什么变化。
  • 从上面的推导条件,我们可以得出简并半导体和非简并半导体的概念:
    • Ev+3kTEFEc3kTE_{v}+3kT\leq E_{F}\leq E_{c}-3kT ,相关载流子的热平衡浓度可以简化成上面的式子,称之为“非简并半导体”(也就是 EFE_{F} 大概在能带的中间)
    • 反之:称为简并半导体。
    • 那么,什么时候 EFE_{F} 离施主/受主能级很近呢?我们考虑杂质浓度增加的情形。
      • 如果杂质浓度增加到载流子开始互相作用的程度,分立的施主/受主能级就会分裂成能带,而能带可能与导带底相交叠。
      • 如果导带中的电子浓度超过 NcN_{c} /价带中的空穴浓度超过 NvN_{v} 时,费米能级就位于导带/价带内部(看上面的两个方程)
      • 此时材料部分表现金属特性
    • 因此,简并半导体/非简并半导体这组概念可以从两个角度来看待:
      • EFE_{F} 也即费米能级的位置
      • 半导体中杂质原子的浓度
      • 而且上述两个角度是相互关联的
本征半导体
本征载流子浓度
  • 我们将上述式子应用于本征半导体,于是有:
n0=ni=Ncexp((EcEEi)kT)p0=pi=ni=Nvexp((EEiEv)kT)\begin{aligned} n_0 &= n_i = N_c \exp\left( \frac{-(E_c - E_{E_i})}{kT} \right) \\ p_0 &= p_i = n_i = N_v \exp\left( \frac{-(E_{E_i} - E_v)}{kT} \right) \end{aligned}
  • 两式相乘:
ni2=NcNvexp[(EcEv)kT]=NcNvexp[EgkT]n_{i}^{2}=N_{c} N_{v} \exp{\left[ {\frac{-( E_{c}-E_{v} )} {k T}} \right]}=N_{c} N_{v} \exp{\left[ {\frac{-E_{g}} {k T}} \right]}

or

ni=NcNveEG/2kTn_{i}=\sqrt{ N_{c}N_{v} }e^{ -E_{G}/2kT }
  • T=300KT=300K 时,对 SiSi 而言 ni=1010cm3n_{i}=10^{10}cm^{-3}
  • 影响 nin_{i} 的因素:
    • TT 越大, nin_{i} 越大
    • EgE_{g} 越大, nin_{i} 越小
本征费米能级位置
  • 利用电子浓度和空穴浓度相等进行计算:
Ncexp[(EcEFi)kT]=Nvexp[(EFiEv)kT]N_{c} \exp\left[ {\frac{-( E_{c}-E_{F i} )} {k T}} \right]=N_{v} \exp\left[ {\frac{-( E_{F i}-E_{v} )} {k T}} \right]
  • 解得
EFi=12(Ec+Ev)+12kTln(NvNc)E_{F i}={\frac{1} {2}} \left( E_{c}+E_{v} \right)+{\frac{1} {2}} \, k T \ln\left( {\frac{N_{v}} {N_{c}}} \right)
  • NvN_{v}NcN_{c} 代入:
EFi=12(Ec+Ev)+34kTln(mpmn)E_{F i}=\frac{1} {2} \left( E_{c}+E_{v} \right)+\frac{3} {4} \, k T \ln\left( \frac{m_{p}^{*}} {m_{n}^{*}} \right)
  • 由于电子和空穴的有效质量与状态密度函数相关,因此也跟本征费米能级(电子浓度与空穴浓度相等)有关。但在实际的计算中,上式的后一项因为太小,常常可以忽略。故而有:
Ei=Ec+Ev2 ⁣+ ⁣3kT4ln ⁣(mpmn) ⁣ ⁣Ec+Ev2E_{i}=\frac{E_{c}+E_{v}} {2} \!+\! \frac{3 k T} {4} \ln\! \left( \frac{m_{p}^{\ast}} {m_{n}^{\ast}} \right) \! \cong\! \frac{E_{c}+E_{v}} {2}
非本征半导体
载流子浓度
  • 我们使用本征半导体的公式推理出 NcN_{c}NvN_{v} 的表达式,之后把它代入非本征半导体的公式中。因为我们想通过本征求得非本征。(非本征的这些量看上去就不怎么好求的样子)于是有:
NC=nie(EcEi)/kTNV=nie(EiEv)/kT\begin{array} {c} {{N_{\mathrm{C}} \,=\, n_{\mathrm{i}} e^{( E_{\mathrm{c}}-E_{\mathrm{i}} ) / k T}}} \\ {{N_{\mathrm{V}} \,=\, n_{\mathrm{i}} e^{( E_{\mathrm{i}}-E_{\mathrm{v}} ) / k T}}} \end{array}

代入到非本征/一般半导体中:

n=nie(EFEi)/kTn \,=\, n_{\mathrm{i}} e^{( E_{\mathrm{F}}-E_{\mathrm{i}} ) / k T} p=nie(EiEF)/kTp \,=\, n_{\mathrm{i}} e^{( E_{\mathrm{i}}-E_{\mathrm{F}} ) / k T}

(这两个式子中的 EiE_{i} 是本征半导体的费米能级)

  • 之后我们把两式相乘,于是得到
pn=ni2pn=n_{i}^2

注意:由于上述式子在推理过程中均使用了玻尔兹曼近似,所以上述公式既适用于本征半导体材料,也适用于非简并的杂质半导体材料

掺杂浓度与载流子浓度的关系
  • 通过热平衡条件我们知道,半导体中不可能存在电场。因为如果存在电场,就会出现载流子的运动,也就不符合热平衡条件的假定了。既然不存在电场,那么半导体中也就没有净电荷。故而我们可以得到==均匀掺杂半导体材料==的内部电中性条件:
p    n  +  ND+    NA  =  0p \;-\; n \;+\; N_{\mathrm{D}}^{+} \;-\; N_{\mathrm{A}}^{-} \;=\; 0

其中 ND+N^+_{D} 是单位体积内电离施主数, NAN_{A}^- 是单位体积内电离受主数。假设在室温条件下,半导体中受主和施主全部电离,也即:

ND+=NDNA=NA\begin{matrix} N_{\mathrm{D}}^{+}=N_{\mathrm{D}} \\ N_{\mathrm{A}}^{-}=N_{\mathrm{A}} \end{matrix}

于是我们得到

p    n  +  ND    NA  =  0p \;-\; n \;+\; N_{\mathrm{D}} \;-\; N_{\mathrm{A}} \;=\; 0
  • 将这个方程与 np=ni2np=n_{i}^2 联立,并且解方程,可以得到:
n=NDNA2+[(NDNA2)2+ni2]1/2n \,=\, \frac{N_{\mathrm{D}} \,-\, N_{\mathrm{A}}} {2} \,+\, \left[ \left( \frac{N_{\mathrm{D}} \,-\, N_{\mathrm{A}}} {2} \right)^{2} \,+\, n_{\mathrm{i}}^{2} \, \right]^{1 / 2} p=ni2n=NA    ND2+[(NA    ND2)2+  ni2]1/2p \,=\, {\frac{n_{\mathrm{{i}}}^{2}} {n}}={\frac{N_{\mathrm{{A}}} \;-\; N_{\mathrm{{D}}}} {2}} \,+\, \left[ \left( {\frac{N_{\mathrm{{A}}} \;-\; N_{\mathrm{{D}}}} {2}} \right)^{2} \,+\; n_{\mathrm{{i}}}^{2} \, \right]^{1 / 2}
  • 性质
    • 本征半导体 NA=0N_{A}=0ND=0N_{D}=0 ,简化为 n=p=nin=p=n_{i}
    • NDNANDniN_{D}-N_{A}\simeq N_{D}\gg n_{i} 或者相反的掺杂半导体。此时 nNDn\simeq N_{D}p=ni2np=\frac{n_{i}^2}{n}
    • niNDNAn_{i}\gg\mid N_{D}-N_{A}\mid ,当温度升到足够高时,本征激发占主导地位,此时 npnin\simeq p \simeq n_{i} 。在温度足够高,niNDNAn_{i}\gg\mid N_{D}-N_{A}\mid 时,所有半导体都将成为本征半导体。
    • NDNA=0N_{D}-N_{A}=0 ,类似本征半导体材料,如果 NAN_{A}NDN_{D} 二者是可比的但不相等时,这种材料称为补偿半导体,这时不能舍去二者之中的任何一个。
载流子浓度与温度的关系
  • 两种影响平衡载流子数的机制(以施主掺杂半导体为例,导带中的电子来自施主电子价带电子):
    • T=0KT=0K ,不论是施主电子还是价带电子都无法激发
    • 低温 TT ,施主电子不断激发
    • 中等温度 TT ,此时施主状态弱束缚内的全部电子都被释放, nNDn\simeq N_{D} 并且进入非本征区,此时价带电子不断被激发,但没有达到 NDN_{D} 的数量级
    • 最终电子被激发穿过带隙数达到 NDN_{D} 的数量级,此时来自施主的电子数基本等于杂质的浓度
    • 示意图
费米能级的位置(对于掺杂/非简并/全部电离的半导体)
  • 直接通过上面方程解出(这几个量都是知道了其中的一个才能求另一个)
EF    Ei=  kTln(n/ni)  =  kTln(p/ni)E_{\mathrm{F}} \;-\; E_{\mathrm{i}} \,=\; k T \, \ln( n / n_{\mathrm{i}} ) \;=\;-\, k T \, \ln( p / n_{\mathrm{i}} )
  • 有特殊情况,比如在合适的温度下
EFEi=kTln(NDni)EiEF=kTln(NAni)If NDNA,NDniIf NAND,NAni\begin{aligned} E_{\mathrm{F}} - E_{\mathrm{i}} &= k T \ln\left( \frac{N_{\mathrm{D}}}{n_{\mathrm{i}}} \right) \\ E_{\mathrm{i}} - E_{\mathrm{F}} &= k T \ln\left( \frac{N_{\mathrm{A}}}{n_{\mathrm{i}}} \right) \end{aligned} \qquad \begin{aligned} & \text{If } N_{\mathrm{D}} \gg N_{\mathrm{A}}, N_{\mathrm{D}} \gg n_{\mathrm{i}} \\ & \text{If } N_{\mathrm{A}} \gg N_{\mathrm{D}}, N_{\mathrm{A}} \gg n_{\mathrm{i}} \end{aligned}
  • 非简并情况下,随着施/受主杂质的增加,费米能级从 EiE_{i} 能量有规律地向上/向下移动。如图所示:
  • TT 增大, EFE_{F}EiE_{i} 离得越近,因为 TT 增大 nin_{i} 也增大。

这个费米能级到底TM的是什么东西?

  • 费米能级表征了某个能级,比它小的能级大多数被电子填满,而它上面的能级大多数为空。
  • 当我们把两块半导体材料放在一起的时候,我们如何判断载流子的浓度呢?通过观察费米能级的大小,我们可以判断载流子的流向。这种载流子的流动以热平衡状态为止。
  • 当达到热平衡状态时,两块材料的费米能级相同

载流子的输运

相关数值:

  • 300K300K 时, kTq=0.0259V\frac{kT}{q}=0.0259V
  • NTN_{T} 是 R-G 中心的数量

能带图与费米能级

  • 当我们考虑载流子沿着 xx 轴输运的时候,能带能量就变为了位置的函数(能带弯曲)
  • 能带图的本质是能带标出了半导体中允许的电子能量与位置的关系。

电子/空穴的动能和势能以及有外电场时能带的弯曲情况

^397caf

  • 考虑某个电子吸收能量超过带隙宽度 EGE_{G} ,则其动能为:EEcE-E_{c} 。如果把空穴视作电子的反面(也即电子获得动能时,同时产生了一个空穴,二者能量是守恒(?)的),那么空穴的动能就是 EvEE_{v}-E
  • 电子的势能: EcErefE_{c}-E_{ref}
  • 那么静电势 VV 与势能的关系就是: qV=EcEref-qV=E_{c}-E_{ref}注意这里的负号!
  • 场强则是电势下降的梯度: ε=dVdx\varepsilon=-\frac{d V} {d x}注意这里的负号!
  • 联立两式,则有:(由于上面那两个式子中的负号可以相互抵消,所以我们可以直接从能带图的斜率近似得到电场和位置的关系
E=1qdEcdx=1qdEvdx=1qdEidx{\mathcal{E}}={\frac{1} {q}} {\frac{d E_{\mathrm{c}}} {d x}}={\frac{1} {q}} {\frac{d E_{\mathrm{v}}} {d x}}={\frac{1} {q}} {\frac{d E_{\mathrm{i}}} {d x}}

可以这样表示的原因:EcE_{c}EvE_{v}EiE_{i} 之间值相差一个常数。

  • 注意此时费米能级的弯曲程度与导带和价带相同。

在扩散情况下费米能级的稳定性

  • 如果非均匀掺杂的半导体的平衡条件成立,那么它的费米能级是与位置无关的。
  • 更深入地来看,费米能级表征了某个状态被电子占据的几率。在平衡条件下,能量的填充态应该是相同的,否则载流子会在不同几率的状态下迁移,由此产生静电流——这就违反了平衡的条件。
  • 下面的示意图是非均匀掺杂半导体的相关能带图像。

载流子的漂移

  • 带电粒子在外电场作用下的运动,并且与电离的杂质原子和热振动的晶格原子不断地与载流子发生碰撞(载流子的加速度不断改变)(当然,还有它本身的热运动)
  • 宏观观测到的是固定的漂移速度项( vdv_{d}
  • 散射机制
    • 晶格散射
    • 电离杂质(施主/受主)散射
    • 晶体中电子与原子的内部场——有效质量

漂移电流

  • 漂移电流(单位时间内流过垂直于电流流动方向任意平面的电荷数):
IPdrift=qpvdAI_{P\mid drift}=qpv_{d}A

(以空穴漂移电流为例)

  • 电流密度 JJ :(以空穴漂移电流为例)
JPdrift=qpvdJ_{P\mid drift}=qpv_{d}
  • 可以看出,得到漂移速度跟电场的关系很重要。我们可以得到 vdv_{d} 跟电场的函数关系。(其中 E\mathcal{E} 是电场,在硅中,空穴 β1\beta \cong 1,电子 β2\beta \cong 2 )(以空穴漂移电荷为例)
υd=μ0E[1+(μ0Eυsat)β]1/β={μ0Eif E0υsatif E\upsilon_{\mathrm{d}} = \frac{\mu_{0} \mathcal{E}}{\left[ 1 + \left( \frac{\mu_{0} \mathcal{E}}{\upsilon_{\mathrm{s a t}}} \right)^\beta \right]^{1 / \beta}} = \begin{cases} \mu_{0} \mathcal{E} & \quad \text{if } \mathcal{E} \to 0 \\ \upsilon_{\mathrm{s a t}} & \quad \text{if } \mathcal{E} \to \infty \end{cases}
  • 在低电场下,μ0\mu_{0} 作为比例常数,同时 vdv_{d}E\mathcal{E} 成比例;在高电场下,vdv_{d} 接近饱和速度 vsatv_{sat} ,漂移电流不显示对电场的依赖性。
  • 在低电场中:(以空穴漂移电荷为例)
Jpdrift=qμppE{\bf J}_{\mathrm{p \mid d r i f t}} \,=\, q \mu_{p} p \mathcal{E}

电子的漂移电流公式也类似,除了将 pp 变成 nn 。虽然最终的公式相同,但我们也不能忽略对于电子而言: vd=μ0Ev_{d}=-\mu_{0}\mathcal{E} ,并且 JNdrift=qnvdJ_{N\mid drift} =-qnv_{d} 。(扩散电流与电子运动方向相反)

迁移率

  • μn\mu _n 或者 μp\mu_{p} ,指移动的自由度,迁移率越高,载流子在半导体内受到的碰撞越少
  • 单位 cm2/Vscm^2 / V \cdot s 或者 m2/Vsm^2 / V \cdot s
  • 与散射的关系:从动量角度列式 vd=qτEmv_{d}=\frac{q\tau E}{m^*}τ\tau 是两次碰撞的平均自由时间。我们假设低电场,于是有 vd=μ0Ev_{d}=\mu_{0}{E} ,于是有 μ0=qτm\mu_{0}=\frac{q\tau}{m^*}
  • 与掺杂的关系:对硅而言,在低掺杂浓度下,迁移率与掺杂浓度无关,超过 1015/cm310^{15}/cm^3 ,随掺杂浓度增大而单调减少。
  • 与温度的关系:
    • 晶格散射/声子散射:与原子的热运动有关。设 μL\mu_{L} 是只有晶格散射存在时的迁移率,则有 μLT3/2\mu_{L} \propto T^{-3 / 2}
    • 电离杂质散射:与掺杂浓度有关。室温下电子/空穴与电离杂质之间存在库伦作用。引起的碰撞与散射会改变载流子的速度特性。设 μI\mu_{I} 是只有电离杂质散射存在时的迁移率。则有 μIT+3/2NI\mu_{I} \propto\frac{T^{+3 / 2}} {N_{I}}
    • 温度升高时可能存在的过程
      • 低温区:温度升高,电子被电离,浓度上升,迁移率上升
      • 中温区(非本征区):温度升高,但杂质已经被完全电离,电子浓度不变。但此时随温度升高晶格散射上升,因此迁移率下降(当然此时电离杂质散射下降,但仍然是晶格散射占主导)
      • 高温区:本征载流子浓度增加(见上一章),开始主导电子浓度,随温度升高,迁移率增加。
      • 对于高掺杂样品,电离杂质散射比较显著。在低温区,电离杂质散射在全部散射中占的比例较大,因此在这个区域内迁移率与温度的依赖关系的斜率减小。
  • 两种散射过程以及迁移率的关系:
    • 微分时间 dtdt 中受到散射的概率:dtτ=dtτI+dtτL\frac{\mathrm{d} t} {\tau}=\frac{\mathrm{d} t} {\tau_{I}}+\frac{\mathrm{d} t} {\tau_{L}}
    • 迁移率(马西森定则): 1μ=1μI+1μL\frac{1} {\mu}=\frac{1} {\mu_{I}}+\frac{1} {\mu_{L}}

电阻率

  • E=ρJE=\rho J ρ\rho 称为电阻率; J=σEJ=\sigma Eσ=1ρ\sigma=\frac{1}{\rho} ,称为电导率。
  • 均匀材料中电流密度与电场的关系
Jdrift=JNdrift+Jpdrift=q(μnn+μpp)E{\bf J}_{\sf d r i f t} \,=\, {\bf J}_{\sf N \mid d r i f t} \,+\, {\bf J}_{\sf p \mid d r i f t} \,=\, q ( \mu_{n} \, n \,+\, \mu_{p} p ) {\bf\cal E}
  • 于是得到电阻率为
ρ=1qμnNDn型半导体\rho = \frac{1}{q \mu_{\mathrm{n}} N_{\mathrm{D}}} \quad \dots \quad \text{n型半导体} ρ=1qμpNAp型半导体\rho = \frac{1}{q \mu_{\mathrm{p}} N_{\mathrm{A}}} \quad \dots \quad \text{p型半导体}

测量电阻率

  • 四探针法,省去了接触电阻对测量的影响。
  • 示意图:
  • 相关计算公式
ρ=VI2πsf(d,t)\rho=\frac{V} {I} \cdot2 \pi s \cdot f ( d, t )

其中 f(d,t)f(d,t) 称为校正因子。

霍尔效应

  • 原理图 相关公式:
qEy=qvxBzq \mathrm{E}_{\mathrm{y}}=q v_{x} B_{z} VH=+EHWV_{H}=+\mathrm{E}_{H} W

Ey=EHE_{y}=E_{H} ,于是:

VH=vxWBzV_{H}=v_{x}WB_{z}

其中, pp 型半导体中的空穴漂移速度(注意, nn 型半导体中这个式子是有负号的):

vx=Jxqp=Ixqp(Wd)v_{x}=\frac{J_{x}}{qp}=\frac{I_{x}}{qp(Wd)}

求解霍尔电压 VHV_{H} 与电流 IxI_{x} ,磁场 BzB_{z} 的关系:

VH=IxBzqpdV_{H}={\frac{I_{x} \, B_{z}} {q p d}}

空穴浓度:

p=IxBzqdVHp=\frac{I_{x} B_{z}} {q d V_{H}}

电子浓度:

n=IxBzqdVHn=-\frac{I_{x} B_{z}} {q d V_{H}}

求出载流子浓度之后,利用空穴漂移速度 vx=μpEx=μpVxLv_{x}=\mu_{p}E_{x}=\mu_{p} \frac{V_{x}}{L} ,将其代入上面的式子。(条件:弱电场) 得到空穴迁移率:

μp=IxLqpVxWd\mu_{p}=\frac{I_{x} \, L} {q p V_{x} \, \, W d}

电子迁移率:

μn=IxLqnVxWd\mu_{n}=\frac{I_{x} \, L} {q n V_{x} \, W d}

载流子的扩散

  • 由于无规则的热运动,粒子从浓度高的区域向浓度低的区域在宏观尺度上的移动,结果使得粒子重新分布。
  • 扩散与粒子电荷的变化、热运动以及粒子间的斥力无关。

扩散电流

  • 菲克第一定律
Jn,diff = qDn dndxJ_{\mathrm{n, d i f f}} ~=~ q D_{n} ~ {\frac{d n} {d x}} J p,diff = qDp dpdxJ_{\mathrm{\tiny~ p, d i f f}} \ =\ -q D_{p} \ {\frac{d p} {d x}}

其中 DnD_{n}DpD_{p} 是电子与空穴的扩散系数,单位是 cm2/scm^2 /s

  • 总电流:
JP=JPdrift+JPdiffJ_{P}=J_{P\mid drift}+J_{P|diff} JN=JNdrift+JNdiffJ_{N}=J_{N\mid drift}+J_{N\mid diff} J=JN+JPJ=J_{N}+J_{P}

扩散系数和迁移率的关系

  • 前提:在平衡条件下的非均匀掺杂半导体
  • 由于是平衡条件,所以总电流为零。又因为半导体是非均匀的,所以此时半导体中存在内建电场,并且内建电场不为零。
  • 我们令:
JNdnft+JNdiff=qμnnE+qDNdndx=0J_{\mathrm{N | d n f t}} \,+\, J_{\mathrm{N | d i f f}} \,=\, q \mu_{\mathrm{n}} \, n {\mathcal{E}} \,+\, q D_{\mathrm{N}} {\frac{d n} {d x}}=\, 0

我们首先看电场那一项:

E=1qdEcdx=1qdEvdx=1qdEidx{\mathcal{E}}={\frac{1} {q}} {\frac{d E_{\mathrm{c}}} {d x}}={\frac{1} {q}} {\frac{d E_{\mathrm{v}}} {d x}}={\frac{1} {q}} {\frac{d E_{\mathrm{i}}} {d x}}

再代入:n=nie(EFEi)/kTn \,=\, n_{\mathrm{i}} e^{( E_{\mathrm{F}}-E_{\mathrm{i}} ) / k T}dndx=nikTe(EFEi)/kTdEidx=nkTdEidx\frac{dn}{dx}=-\frac{n_{i}}{kT}e^{ (E_{F}-E_{i})/kT } \frac{dE_{i}}{dx}=-\frac{n}{kT} \frac{dE_{i}}{dx} 。 将上述两式代入,得到对于电子的爱因斯坦关系式:

DNμn=kTq\frac{D_{\bf N}} {\mu_{\bf n}}=\frac{k T} {q}

对于空穴

Dpμp=kTq{\frac{D_{\mathrm{p}}} {\mu_{\mathrm{p}}}}={\frac{k T} {q}}

上述关系式在非平衡/非简并状态下也有效。

载流子的产生-复合

  • 当半导体的平衡状态受到微扰时,半导体内部的载流子浓度会涨落。复合-产生是一种恢复机制。
  • 分析工作时非平衡的器件。
  • 二者是一组相关的过程:
    • 复合:电子和空穴被湮灭/消失的过程
    • 产生:电子/空穴被创建的过程

几种R-G过程

  • 直接复合/产生:导带电子和价带空穴的直接湮灭/电子被激发后,直接从价带进入导带;这个过程会产生过剩的能量释放/吸收热能或者光能
  • 间接复合/产生(R-G中心复合/产生):通过中介物进行的复合过程,只在半导体内复合中心(晶格缺陷/特殊的杂质原子,相比受主和施主浓度很低)的特殊位置产生。在靠近带隙中心引入了允许电子能级 ETE_{T} ,不同杂质产生的能级不尽相同。复合:首先一种载流子漂移进入R-G中心附近,被相关势阱捕获,失去能量;之后出现另一种载流子,二者互相吸引,失去能量,在中心内与电子一同湮灭。也有可能,电子失去能量后,在价带与空穴一起湮灭。产生过程也类似,热能作为中间媒介,在R-G中心产生载流子。
  • 俄歇过程:两个同种类型的载流子发生碰撞,发生直接复合,释放的能量传递给经过碰撞保存下来的那个,后者与晶格碰撞失去能量,并产生热量。
  • 碰撞电离:高能载流子与晶格碰撞,释放能量并产生一对电子和空穴。(常常在外电场 E\mathcal{E} 比较大的地方产生)

不同R-G过程在不同半导体中是否占据主导

  • 直接/间接R-G过程:室温下非简并掺杂半导体的低 E\mathcal{E} 区域
    • 发射光子:能量+很少量的动量,光助跃迁在 EkE-k 图中基本垂直
    • 发射声子:很少量的能量+动量,声助跃迁在 EkE-k 图中基本水平
    • 直接带隙半导体:产生/复合时,能量变化+很小的动量变化
    • 间接带隙半导体:产生/复合时,能量变化+较大的动量变化
    • 在直接带隙半导体中,直接复合以很快的速度产生。
    • 所以,如果在间接带隙半导体中进行直接复合,必须伴有声子+光子的吸收与发射。事实上,间接带隙半导体中,R-G中心复合(间接复合)机制是主要的复合过程。
  • 俄歇复合:器件结构的高掺杂区域
  • 碰撞电离:器件的高 E\mathcal{E} 区域

对R-G过程的数学描述

  • R-G过程不导致电荷的传输,但会改变载流子浓度。
光产生
ntlight=ptlight= GL(x,λ)\left. \frac{\partial n} {\partial t} \right|_{\mathrm{l i g h t}}=\left. \frac{\partial p} {\partial t} \right|_{\mathrm{l i g h t}}=\ G_{\mathrm{L}} ( x, \lambda)

GLG_{L} 光产生率:(1/cm3s)(1个/cm^3\cdot s)

GL(x,λ)=GL0eαxG_{\mathrm{L}} ( x, \lambda) \,=\, G_{\mathrm{L 0}} e^{-\alpha x}

GL0G_{L{0}}x=0x=0 处的光产生率,α\alpha 是吸收系数,与材料和波长有关。

特殊情况下的R-G过程(间接热复合-产生)
  • 产生的微扰必须是少数载流子“小注入”
    • Δpn0\Delta p\ll n_{0} ,且 nn0n\simeq n_{0} (在 nn 型材料内)
    • Δnp0\Delta n\ll p_{0} ,且 pp0p \simeq p_{0} (在 pp 型材料内)
    • 注意,其中少数载流子可以比平衡浓度下的少数载流子的浓度高。
  • 考虑 R-G 过程下的小注入导致的微扰的过剩导致的少数载流子的湮灭
    • 由于是小注入,少数载流子远小于多数载流子,所以电子总能迅速到 R-G 能级上与空穴结合。所以单位时间内少数载流子的变化与 R-G 中心数量 NTN_{T} 成比例。
    • 同时,少数载流子的变化也与少数载流子本身的浓度成比例。
    • 引入比例常数 cpc_{p}
    • 考虑到少数载流子不断减小,所以在公式前面加负号。
    • 故有公式:ptR=cpNTp\left. \frac{\partial p} {\partial t} \right|_{\mathrm{R}}=-c_{p} N_{\mathrm{T}} p
  • 考虑少数载流子的产生
    • 在实验上这个数量只与R-G中心的数量有关。
    • 故有公式:ptG=cpNTp0\left. \frac{\partial p} {\partial t} \right|_{\mathrm{G}} \,=\, c_{\mathrm{p}} \, N_{\mathrm{T}} p_{0}
  • 两者相加: nn 型材料中的空穴:
ptRGithermal= cpNTΔp\frac{\partial p} {\partial t} \Bigg|_{\stackrel{\mathrm{i-t h e r m a l}} {\mathrm{R-G}}}=\ -c_{\mathrm{p}} N_{\mathrm{T}} \Delta p

pp 型材料中的电子:

ntRGi ithermal =  cnN TΔn\frac{\partial n} {\partial t} \Bigg|_{\stackrel{\scriptstyle i \mathrm{\tiny~ i-t h e r m a l}} {\scriptstyle R-G}} \ =\ -\ c_{n} N_{\mathrm{\tiny~ T}} \Delta n
  • 时间常数
τp=1cpNT\tau_{\mathrm{p}}={\frac{1} {c_{\mathrm{p}} N_{\mathrm{T}}}} τn=1cnNT\tau_{\mathrm{n}}={\frac{1} {c_{\mathrm{n}} N_{\mathrm{T}}}}
  • 代入
ptRGithermal  =  Δpτp\frac{\partial p} {\partial t} \Bigg|_{\stackrel{\mathrm{i-t h e r m a l}} {\mathrm{R-G}}} \;=\;-\frac{\Delta p} {\tau_{p}} ntRGithermal  =  Δnτn\frac{\partial n} {\partial t} \Bigg|_{\stackrel{\mathrm{i-t h e r m a l}} {\mathrm{R-G}}} \;=\;-\frac{\Delta n} {\tau_{n}}
  • 在更一般的情况下
ptRG=ntRG=ni2npτp(n+n1)+τn(p+p1)\left. \frac{\partial p} {\partial t} \right|_{R-G}=\left.\frac{\partial n} {\partial t} \right|_{R-G}=-\frac{n_{i}^{\, 2}-n p} {\tau_{p} ( n+n_{1} )+\tau_{n} ( p+p_{1} )}

其中

n1=niexp[(ETEi)/kT]p1=niexp[(EiET)/kT]\begin{array} {c} {n_{1}=n_{i} e x p [ ( E_{T}-E_{i} ) / k T ]} \\ {p_{1}=n_{i} e x p [ ( E_{i}-E_{T} ) / k T ]} \\ \end{array}
少子寿命
  • τn\tau_{n}τp\tau_{p} 是少数载流子的平均寿命,也就是:少数载流子在多子海洋中存活的平均时间。
  • 寿命测量:时域光电导,连一块半导体到电路里,照光,测电导的变化
υLVA2(1+Δσσdark)=VL+υ0et/τp\upsilon_{\mathrm{L}} \, \equiv\, {\frac{V_{\mathrm{A}}} {2}} \left( 1 \,+\, {\frac{\Delta\sigma} {\sigma_{\mathrm{d a r k}}}} \right) \,=\, V_{\mathrm{L}} \,+\, \upsilon_{0} \, e^{-\, t / \tau_{\mathrm{p}}}

后者是瞬态衰减项。

状态方程

  • 连续性方程
nt=ntdrift+ntdiff+ntthermalRG+nt(light.etc)other\frac{\partial n} {\partial\, t} \,=\, \frac{\partial n} {\partial\, t} \bigg|_{d r i f t} \,+\, \frac{\partial n} {\partial\, t} \bigg|_{d i f f} \,+\, \frac{\partial n} {\partial\, t} \bigg|_{t h e r m a l \atop R-G}+\left. \frac{\partial n} {\partial\, t} \right|_{\stackrel{o t h e r} {(li g h t. e t c )}} pt=ptdrift+ptdiff+ptRGthermal+pt(light.etc)other\frac{\partial p} {\partial t} \,=\, \frac{\partial p} {\partial t} \bigg|_{d r i f t} \,+\, \frac{\partial p} {\partial t} \bigg|_{d i f f} \,+\, \frac{\partial p} {\partial t} \bigg|_{\stackrel{t h e r m a l} {R-G}} \,+\, \frac{\partial p} {\partial t} \bigg|_{\stackrel{o t h e r} {{( l i g h t. e t c )}}}

代入后有:

nt=1qJn(x)xΔnτn+GL\frac{\partial n} {\partial t}=\frac{1} {q} \, \frac{\partial J_{n} ( x )} {\partial x}-\frac{\Delta n} {\tau_{n}}+G_{L} pt=1qJp(x)xΔpτp+GL{\frac{\partial p} {\partial t}} \,=\,-\, {\frac{1} {q}} \, {\frac{\partial J_{p} ( x )} {\partial x}}-{\frac{\Delta p} {\tau_{p}}} \,+\, G_{L}

少子扩散方程

  • 条件
    • 一维
    • 少子
    • 电场约为 00 (舍去漂移项)
    • 平衡少子浓度不是位置的函数
    • 小注入条件
    • 间接热复合-产生(间接带隙半导体)
    • 排除光产生
  • 少子扩散方程
Δnpt=DN2Δnpx2Δnpτn+GLΔpnt=Dp2Δpnx2Δpnτp+GL\begin{aligned} \frac{\partial \Delta n_{\mathrm{p}}}{\partial t} &= D_{\mathrm{N}} \frac{\partial^2 \Delta n_{\mathrm{p}}}{\partial x^2} - \frac{\Delta n_{\mathrm{p}}}{\tau_{\mathrm{n}}} + G_{\mathrm{L}} \\ \frac{\partial \Delta p_{\mathrm{n}}}{\partial t} &= D_{\mathrm{p}} \frac{\partial^2 \Delta p_{\mathrm{n}}}{\partial x^2} - \frac{\Delta p_{\mathrm{n}}}{\tau_{\mathrm{p}}} + G_{\mathrm{L}} \end{aligned}
  • 条件与简化
    • 稳态: Δnpt0 (Δpnt0) {\frac{\partial\Delta n_{\mathrm{p}}} {\partial t}} \to0 \quad{\mathrm{~ \left( {\frac{\partial\Delta p_{\mathrm{n}}} {\partial t}} \to0 \right) ~}}
    • 无浓度梯度/扩散电流: DN2Δnpx20  (Dp2Δpnx20)D_{\mathrm{N}} {\frac{\partial^{2} \Delta n_{\mathrm{p}}} {\partial x^{2}}} \rightarrow0~ ~\left( D_{\mathrm{p}} {\frac{\partial^{2} \Delta p_{\mathrm{n}}} {\partial x^{2}}} \rightarrow0 \right)
    • 无热 R-G 时: Δnpτn0(Δpnτp0)\frac{\Delta n_{\mathrm{p}}} {\tau_{\mathrm{n}}} \rightarrow0 \quad\quad\left( \frac{\Delta p_{\mathrm{n}}} {\tau_{\mathrm{p}}} \rightarrow0 \right)
    • 无光照时:GL0G_{L}\to0 (或瞬时光照)

扩散长度

  • 少子在被湮灭前能够在大量多子内扩散的平均距离
  • nn 型材料的少子空穴
LP=DPτpL_{P}=\sqrt{ D_{P}\tau_{p} }
  • pp 型材料的少子电子
LN=DNτnL_{N}=\sqrt{ D_{N}\tau_{n} }
  • 计算的前提:稳态,无光照

准费米能级-非平衡态描述载流子浓度

  • 将能带图推广到非平衡态,也就是“==非热平衡状态==”
  • 非平衡态的非简并半导体:
nnie(FNEi)/kT或   FNEi+kTln ⁣(nni)pnie(EiFP)/kT或   FP=EikTln ⁣(pni)\begin{array} {l l} {n \, \equiv\, n_{\mathrm{i}} e^{( F_{\mathrm{N}}-E_{\mathrm{i}} ) / k T}} & {或 \ \ \ F_{\mathrm{N}} \, \equiv\, E_{\mathrm{i}} \,+\, \, k T \, \operatorname{l n} \! \left( \frac{n} {n_{\mathrm{i}}} \right)}\\ {p \, \equiv\, n_{\mathrm{i}} e^{( E_{\mathrm{i}}-F_{\mathrm{P}} ) / k T}} & {或 \ \ \ F_{\mathrm{P}} \,=\, \, E_{\mathrm{i}} \,-\, \, k T \, \operatorname{l n} \! \left( \frac{p} {n_{\mathrm{i}}} \right)} \\ \end{array}
  • 非平衡到平衡态,也就是 FNEFF_{N}\to E_{F}FPEFF_{P}\to E_{F} 的过程。
  • 准费米能级随位置变化,表示在半导体中有电流流动。

其他输运

  • 量子隧穿
    • 磁隧道结
    • 扫描隧道显微镜
  • 弹道输运
  • 热离子发射

PN结的特性

  • 基本知识
    • P型半导体和N型半导体放在一起的时候,电子/空穴形成扩散,扩散之后形成电场,电场与电子/空穴的扩散作用发生抵消。
    • P区与N区的分界线:冶金结
    • 两边叫做准中性区,中间叫做耗尽层
  • 杂质分布-净杂质浓度 NDNAN_{D}-N_{A} ,实际工艺具有复杂度,进行理想化处理
    • 突变结分布(离子注入,轻掺杂)
    • 线性缓变结分布(中等掺杂,重掺杂)
    • 示意图

零偏

  • 假设
    • 玻尔兹曼分布(非简并半导体)推导出不同能级的载流子浓度
    • 完全电离(室温下完全电离)

内建电压

  • 由于PN结处于热平衡状态,故而半导体的费米能级处处相等,并且是恒定的值。然而P型半导体与N型半导体中,费米能级的位置并不相同。于是,空间电荷区所在位置的导带和价带就要发生弯曲。示意图
  • 内建电势差 VbiV_{bi} ,维持多子/少子的平衡,故而它在半导体内不产生电流。外加探针与半导体之间也会产生相应的电势差(电子/空穴的运动),故而伏安表无法测出。
  • 载流子浓度 → ϕFp\phi_{F_{p}}ϕFn\phi_{F_{n}} 的大小 → VbiV_{bi} 的大小
Vbi=kTqln(NaNdni2)V_{b i}=\frac{k T} {q} \ln\left( \frac{N_{a} N_{d}} {n_{i}^{2}} \right)

其中, NdN_{d}nn 区的施主浓度;NaN_{a}pp 区的受主浓度。

  • 另一种计算 VbiV_{bi} 的方法(利用平衡条件下 J=0J=0 得出)
E=dVdxE=-\frac{d V} {d x}

积分: V(xn)V(xp)=VbiV(x_{n})-V(-x_{p})=V_{bi}

E = DNμndn/dxn= kTqdn/dxn{{E}} \ =\ -\, \frac{D_{\mathrm{N}}} {\mu_{n}} \, \frac{d n / d x} {n}=\ -\, \frac{k T} {q} \, \frac{d n / d x} {n} Vbi=xpxnEdx=kTqn(xp)n(xn)dnn=kTqln(n(xn)n(xp))V_{\mathrm{b i}} = - \int_{-x_{\mathrm{p}}}^{x_{\mathrm{n}}} E \, dx = \frac{k T}{q} \int_{n(-x_{\mathrm{p}})}^{n(x_{\mathrm{n}})} \frac{dn}{n} = \frac{k T}{q} \ln\left( \frac{n(x_{\mathrm{n}})}{n(-x_{\mathrm{p}})} \right)
  • 这个式子相比上式更加普适。上式也是把 xnx_{n}xpx_{p} 处的数值代入后得出的。

耗尽近似

  • 一维泊松方程:
dEdx=ρKSε0=qKSε0(p  n + ND  NA)\frac{d {\cal E}} {d x}=\frac{\rho} {K_{\mathrm{S}} \varepsilon_{0}}=\frac{q} {K_{\mathrm{S}} \varepsilon_{0}} ( p \ -\ n \ +\ N_{\mathrm{D}} \ -\ N_{\mathrm{A}} )
  • 为了方便地得到 nn / pp 的值,需要对实际电荷的分布进行理想化近似
    • 冶金结附近( xpxxn-x_{p}\leq x\leq x_{n} ),与杂质浓度相比,载流子浓度近似忽略不计
    • 耗尽区之外的电荷密度处处为 00
  • 对于一维突变结,引入耗尽近似,电荷密度: (注意,下面的积分也要用耗尽近似:也即 E(xp)E(-x_{p})E(xn)E(x_{n}) 均为 00
ρ={qNAif xpx0qNDif 0xxn0if xxp or xxn\rho = \begin{cases} -q N_{\mathrm{A}} & \quad \text{if } -x_{\mathrm{p}} \leq x \leq 0 \\ q N_{\mathrm{D}} & \quad \text{if } 0 \leq x \leq x_{\mathrm{n}} \\ 0 & \quad \text{if } x \leq -x_{\mathrm{p}} \text{ or } x \geq x_{\mathrm{n}} \end{cases}

两边分别积分:

E(x)=qNAKSε0(xp+x)...xpx0{\mathcal E} ( x ) \,=\,-\, \frac{q N_{\mathrm A}} {K_{\mathrm S} \, \varepsilon_{0}} ( x_{\mathrm p} \,+\, x ) \qquad. \,. \,. \, \,-x_{\mathrm p} \leqslant\, x \leqslant\, 0 E(x)=qNDKSε0(xnx):...0xxn{\mathcal E} ( x ) \,=\,-\, \frac{q N_{\mathrm{D}}} {K_{\mathrm{S}} \varepsilon_{0}} ( x_{\mathrm{n}} \,-\, x ) : \quad. \,. \,. \, 0 \, \le\, x \, \leqslant\, x_{\mathrm{n}}
  • x=0x=0 处电场连续: qNAxp=qNDxnqN_{A}x_{p}=qN_{D}x_{n}
  • 通过 E=dV/dx{\mathcal{E}}=-d V / d x ,可以有方程(注意上下界的选取
dVdx={qNAKSε0(xp  +  x).  .  .  xp  x  0qNDKSε0(xn    x).  .  0xxn\frac{d V} {d x}=\left\{\begin{array} {l l} {{\displaystyle\frac{q N_{\mathrm{A}}} {K_{\mathrm{S}} \varepsilon_{0}} ( x_{\mathrm{p}} \;+\; x ) \quad}} & {{. \;. \;. \;-x_{\mathrm{p}} \leqslant\; x \leqslant\; 0}} \\ {{\displaystyle\frac{q N_{\mathrm{D}}} {K_{\mathrm{S}} \varepsilon_{0}} ( x_{\mathrm{n}} \;-\; x ) \quad}} & {{. \;. \; 0 \, \leqslant\, x \, \leqslant\, x_{\mathrm{n}}}} \end{array} \right.
  • 边界条件:
V=0atx=xpV=Vpiatx=xn\begin{array}{c} V = 0 \quad \text{at} \quad x = -x_{\text{p}} \\ V = V_{\text{pi}} \quad \text{at} \quad x = x_{\text{n}} \end{array}
  • 积分之后:
V(x)=qNA2KSε0(xp+x)2 . ...xpx0V ( x ) \,=\, \frac{q N_{\mathrm{A}}} {2 K_{\mathrm{S}} \varepsilon_{0}} ( x_{p} \,+\, x )^{2} \quad\mathrm{~. ~}. \,. \,. \,-x_{p} \, \leqslant\, x \, \leqslant\, 0 V(x)=VbiqND2KSε0(xnx)2..0xxnV ( x ) \,=\, V_{\mathrm{b i}} \,-\, \frac{q N_{\mathrm{D}}} {2 K_{S} \varepsilon_{0}} ( x_{\mathrm{n}} \,-\, x )^{2} \qquad. \, \,. \, \, 0 \, \leqslant\, x \, \leqslant\, x_{\mathrm{n}}
  • V(0)V(0) 相等,并且使用上面的电荷守恒的条件,可以推出 xpx_{p}xnx_{n} ,以及 W=xn+xpW=x_{n}+x_{p} 的值
Wxn+xp=[2Ksε0q(NA+NDNAND)Vbi]1/2W \equiv\, x_{\mathrm{n}} \,+\, x_{\mathrm{p}} \,=\, \left[ \frac{2 K_{\mathrm{s}} \varepsilon_{0}} {q} \left( \frac{N_{\mathrm{A}} \,+\, N_{\mathrm{D}}} {N_{\mathrm{A}} N_{\mathrm{D}}} \right) V_{\mathrm{b i}} \right]^{1 / 2}

其中:

xp=NDNA+NDWx_{p}=\frac{N_{D}}{N_{A}+N_{D}} W xn=NANA+NDWx_{n}=\frac{N_{A}}{N_{A}+N_{D}} W

偏置电压下的PN结

  • 对于小电流注入,接触电极到耗尽层边界的准中性区上电阻压降忽略不计。电压压降在耗尽层上。
  • 设正偏电压为 VAV_{A} ,耗尽层的电压为: VbiVAV_{bi}-V_{A} (可以把整个电路看作两个电源)
  • 我们可以根据上面的推论对零偏时的电场、电压与耗尽层的长度进行修正。 对 0xxn0\leq x\leq x_{n}
V(x) = Vbi  VA  qND2KSε0(xn  x)2V ( x ) ~=~ V_{\mathrm{b i}} ~-~ V_{\mathrm{A}} ~-~ \frac{q N_{\mathrm{D}}} {2 K_{\mathrm{S}} \varepsilon_{0}} ( x_{\mathrm{n}} ~-~ x )^{2} W=[2KSε0q(NA + NDNAND) ⁣(Vbi  VA)]1/2W=\left[ \frac{2 K_{\mathrm{S}} \varepsilon_{0}} {q} \left( \frac{N_{\mathrm{A}} ~+~ N_{\mathrm{D}}} {N_{\mathrm{A}} N_{\mathrm{D}}} \right) \! ( V_{\mathrm{b i}} ~-~ V_{\mathrm{A}} ) \right]^{1 / 2}

前提:VAVbiV_{A}\leq V_{bi} ,当 VAV_{A} 趋近 VbiV_{bi} 的时候,不能忽略准中性区的压降。

  • 结果分析:
    • VA>0V_{A}>0 耗尽层宽度减少,结内电荷减少,电场降低,电势减少。
    • VA<0V_{A}<0 耗尽层宽度增大,结内电荷增多,电场增大,电势增大。
  • 费米能带图: 由上面的推测,可以得到
EFpEFn=qVAE_{F_{p}}-E_{F_{n}}=-qV_{A}

PN结的I-V特性

  • V=0V=0 时,载流子通过扩散,只有较高能量的电子可以克服势垒,进入结的另一侧。而另一侧,少数载流子进入耗尽层后,却会被电场很快地扫到结的另一侧。平衡时:对电子而言: pp 型到 nn 型的漂移电流抵消 nn 型到 pp 型的扩散电流。
  • VA>0V_{A}>0 ,势垒出现下降,更多载流子可以通过扩散的形式到耗尽层的另外一侧,导致电流 I=IN+IPI=I_{N}+I_{P} (电子电流+空穴电流)。
    • 示意图:
    • 外加电压和电流的关系:外加电压与能级变化成线性,能级变化与载流子浓度成指数关系,而载流子浓度又与电流有线性关系,所以外加电压与电流呈指数关系。
    • 注入的载流子如何湮灭:多子扩散到耗尽层的另外一侧后,借助复合-产生中心发生复合。
  • VA<0V_{A}<0 ,势垒上升,载流子很难通过扩散越过势垒,而少数载流子更容易经由电场通过耗尽区到对面的准中性区。因此此时,反向偏置电流与少数载流子有关,且非常小。
    • 示意图:
    • 外加电压与电流的关系:少数载流子的漂移电流不受势垒高度的影响。与电流大小有关的是每秒进入耗尽层的少数载流子数目。
    • 注入和抽取的载流子如何得到补充:少子进入耗尽层后,借助复合-产生中心产生出的载流子会替代离开的载流子。而过剩的多数载流子会导致一个局部电场,将载流子推向接触电极,使得多余的载流子进入外部电路。
  • 由上述推理,得到 IVI-V 关系的通用形式:
I=I0(eVA/Vref1)I \,=\, I_{0} ( e^{V_{A} / V_{\mathrm{r e f}}} \,-\, 1 )

或者

I=I0(eqVA/kT1)I=I_{0}(e^{ qV_{A}/kT }-1)

定量求解

假设
  • 二极管在稳态
  • 非简并掺杂,突变结
  • 一维
  • 准中性区域,小电流注入
  • GL=0G_{L}=0 ,二极管内除去漂移/扩散/热复合-产生之外没有其他过程。
  • 二极管的准中性区满足少数载流子扩散方程的假设条件, E0\mathcal{E}\cong0 和小电流注入。
耗尽层
  • 假设:耗尽层内载流子没有热复合-产生过程。于是代入连续性方程,有:
dJNdx=dJPdx=0\frac{dJ_{N}}{dx}=\frac{dJ_{P}}{dx}=0

也即:耗尽层中的 JNJ_{N}JPJ_{P} 为常数。可以得到

J=JN(xp)+Jp(xn)J \,=\, J_{\mathrm{N}} (-\, x_{\mathrm{p}} ) \,+\, J_{\mathrm{p}} ( x_{\mathrm{n}} )
准中性区
  • 由于有假设:稳态/ GL=0G_{L}=0 ,于是有少数载流子扩散方程:
0=DNd2Δnpdx2Δnpτn   xxp0=Dpd2Δpndx2Δpnτp   xxn\begin{array} {c} {{0 \,=\, D_{\mathrm{N}} \displaystyle\frac{d^{2} \Delta n_{\mathrm{p}}} {d x^{2}} \,-\, \displaystyle\frac{\Delta n_{\mathrm{p}}} {\tau_{\mathrm{n}}} \ \ \ \, \, \, \, \, x \, \leqslant\,-\, x_{\mathrm{p}}}} \\ {{0 \,=\, D_{\mathrm{p}} \displaystyle\frac{d^{2} \Delta p_{\mathrm{n}}} {d x^{2}} \,-\, \displaystyle\frac{\Delta p_{\mathrm{n}}} {\tau_{\mathrm{p}}} \ \ \ \, \, \, \, \, \, \, x \, \geqslant\, x_{\mathrm{n}}}} \end{array}
  • 由于有假设: E0\mathcal{E}\cong0 以及 dn0dx=dp0dx=0\frac{dn_{0}}{dx}=\frac{dp_{0}}{dx}=0 。则在中性区的载流子电流密度可以化简为
JN=qDNdΔnpdx..xxpJp=qDpdΔpndx...xxn\begin{array} {c l} {{J_{\mathrm{N}} \,=\, q D_{\mathrm{N}} \displaystyle\frac{d \Delta n_{\mathrm{p}}} {d x}}} & {{\quad. \,. \, \, x \, \leqslant\,-\, x_{\mathrm{p}}}} \\ {{{}}} & {{{}}} \\ {{J_{\mathrm{p}} \,=\,-\, q D_{\mathrm{p}} \displaystyle\frac{d \Delta p_{\mathrm{n}}} {d x}}} & {{\quad. \,. \, \,. \, \, x \, \geqslant\, x_{\mathrm{n}}}} \end{array}
  • 边界条件:“宽基层”二极管,二极管的接触电极与耗尽层边界之间的距离是少数载流子扩散长度的几倍或更多
Δnp(x)=0 \Delta n_{\mathrm{p}} ( x \rightarrow-\infty) \,=\, 0 \ \, Δpp(x+)=0\Delta p_{\mathrm{p}} ( x \rightarrow+\infty) \,=\, 0
  • 边界条件:考虑耗尽层边界,并认为边界处准费米能级与 pp / nn 区域的费米能级等同(FNFP=qVAF_{N}-F_{P}=qV_{A}),由准费米能级的两个式子,可以得出PN结定律
np=ni2eqVA/kT.  .  .  xp  x  xnn p \,=\, n_{\mathrm{i}}^{2} e^{q V_{\mathrm{A}} / k T} \qquad. \;. \;. \;-x_{p} \leqslant\; x \; \leqslant x_{n}

由此可以得出:

Δnp(xp)=ni2NA(eqVA/kT1)\Delta n_{\mathrm{p}} (-x_{\mathrm{p}} ) \,=\, \frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{A}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, \, 1 ) Δpn(xn)=ni2ND(eqVA/kT1)\Delta p_{\mathrm{{n}}} ( x_{\mathrm{{n}}} ) \,=\, \frac{n_{\mathrm{{i}}}^{2}} {N_{\mathrm{{D}}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, 1 )

(注意是 Δpn\Delta p_{n} 哦,不是 pnp_{n}

  • 将边界条件代入到准中性区的少数载流子扩散方程,同时设定新的坐标系如图:
  • 可以得到:
Δpn(x)=ni2ND(eqVA/kT1)ex/LP\Delta p_{n} ( x^{\prime} ) \,=\, \frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, 1 ) e^{\,-\, x^{\prime} / L_{\mathrm{P}}} Δnp(x)=ni2NA(eqVA/kT1)ex/LN\Delta n_{\mathrm{p}} ( x^{\prime\prime} ) \,=\, {\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{A}}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, 1 ) e^{\,-\, x' / L_{\mathrm{N}}}
  • 这里也可以严格推导出PN结定律,并求出 I0I_{0} 的值
I0qA(DNLNni2NA+DPLPni2ND)I_{0} \, \equiv\, q A \left( {\frac{D_{\mathrm{N}}} {L_{\mathrm{N}}}} \, {\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{A}}}} \,+\, {\frac{D_{\mathrm{P}}} {L_{\mathrm{P}}}} \, {\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}}} \right)

(注意:在计算这个值的时候,仅需要考虑轻掺杂一侧)

结论

  • 少数载流子扩散电流:在准中性区,电流随着距离耗尽层边界处的距离加大而衰减。(正向偏置)
  • 载流子浓度:
    • 正偏时,少数载流子浓度高于其对应的热平衡值,随着与耗尽层边界的距离的增加,非平衡载流子浓度呈指数形式衰减,在无限远处回归其平衡值。示意图
    • 反偏时,少数载流子浓度低于其平衡值,其余同上。

与理想情况的误差

假设的打破:没有其他过程——反向击穿
  • 击穿是可逆的
  • 击穿电压:(对于单边结), x=0x=0 处的电场( ECRE_{CR} )为 2q(VbiVA)NKsϵ0\sqrt{ \frac{2q(V_{bi}-V_{A})N}{K_{s}\epsilon_{0}} } 。如果反偏电压超过了某个特定的值(反偏电压 VBRV_{BR} ),那么PN结就会被击穿。相关公式
ECR = 2qN ⁣(Vbi + VBR)Ksϵ0E_{C R} ~=~ {\sqrt{\frac{2 q N \! \left( V_{b i} ~+~ V_{B R} \right)} {K_{s} \epsilon_{0}}}} VBR = Ksε0ECR22qNVbiV_{\it B R} ~=~ {\frac{K_{s} \varepsilon_{0} {E_{\it C R}}^{2}} {2 q N}}-V_{b i}
  • 对单边结而言, NN 增大 VBRV_{BR} 降低;带隙宽度越大, VBRV_{BR} 降低。
  • 雪崩击穿:穿过耗尽层时,载流子一边加速一边与半导体晶格发生碰撞,并损失能量。当反向偏压不断增加并接近击穿电压时,每次碰撞传递的能量可以使得半导体原子“电离化”。(使原子释放出一个价带电子,产生一个电子-空穴对,称为碰撞电离)额外产生的载流子与原来的载流子一起在耗尽层中加速,产生出更多载流子,反向电流急剧增加至无穷大。在击穿前的几个伏特处,就存在明显的递增现象。
    • VBR1NBV_{\mathrm{B R}} \propto{\frac{1} {N_{\mathrm{B}}}}
    • VBRV_{BR} 会随温度的上升而增加。因为温度上升导致晶格散射的作用增强,平均自由程减小,需要更大能量使雪崩发生。
  • 齐纳击穿:反向偏置二极管中出现量子隧穿现象。要求势垒另一侧的能级位置存在着允态,势垒宽度要很薄( <10nm<10nm )。
    • 结两侧都是重掺杂的二极管中,齐纳过程会比较显著
    • 击穿电压 VBRV_{BR} 随温度上升而降低。温度升高,价电子更容易摆脱共价键的束缚,可以隧穿的价带电子增多,因此更容易击穿。
假设的打破:假设价带中有电流的复合-产生
  • 在推上面的公式的时候,忽略了耗尽层中的热复合-产生过程。
    • 二极管处于反向偏置的时候,耗尽层中载流子浓度会下降,高于其热平衡值(因为其中只有少数载流子),导致耗尽层中电子和空穴的产生,所以反向电流增大
    • 二极管处于正向偏置的时候,耗尽层中载流子浓度上升,低于其热平衡值(其中流过多数载流子),导致耗尽层中出现电流的复合,因此正向电流减少。
  • 使用R-G过程相关的公式进行计算:
IRG=qAxpxnnp    n12τp(n  +  n1)  +  τn(p  +  p1)dxI_{\mathrm{R-G}} \,=\, q A \, \int_{-x_{p}}^{x_{n}} {\frac{n p \;-\; n_{1}^{2}} {\tau_{p} ( n \;+\; n_{1} ) \;+\; \tau_{n} ( p \;+\; p_{1} )}} d x
  • 反向偏压(耗尽区载流子忽略不计,且反向偏压 >> 几个 kTq\frac{kT}{q}):
IRG =  qAni2τ0WI_{\mathrm{R-G}} ~=~-~ {\frac{q A n_{i}} {2 \tau_{0}}} W

其中:

τ012 ⁣(τpn1ni+τnp1ni)=12 ⁣(τpe(ETEi)/kT+τne(EiET)/kT)\tau_{0} \, \equiv\, {\frac{1} {2}} \! \left( \tau_{\mathrm{p}} {\frac{n_{1}} {n_{\mathrm{i}}}} \,+\, \tau_{\mathrm{n}} {\frac{p_{1}} {n_{\mathrm{i}}}} \right) \,=\, {\frac{1} {2}} \! \left( \tau_{\mathrm{p}} e^{( E_{\mathrm{T}}-E_{\mathrm{i}} ) / k T} \,+\, \tau_{\mathrm{n}} e^{( E_{\mathrm{i}}-E_{\mathrm{T}} ) / k T} \right)
  • 正向偏压(耗尽区载流子不能忽略不计,正向偏压大于 >> 几个 kTq\frac{kT}{q} ):
IRGqAniWeqVA/2kTI_{R-G} \propto q A n_{i} W e^{q V_{A} / 2 k T}
  • 在分析流过二极管的总电流时有公式
I=IDIFF+IRGI=I_{DIFF}+I_{R-G}
  • 反向偏压时, VAV_{A} 越大, WW 越大,因此 IRGI_{R-G} 也就越大。
假设的打破:大电流时,准中性区的压降相比耗尽层的压降不能忽略
  • 需要考虑准中性区上的压降。因此耗尽层上的压降(结电压 VJV_{J} )与 VAV_{A} 有如下关系
VJ = VA  IRSV_{\mathrm{J}} ~=~ V_{\mathrm{A}} ~-~ I R_{\mathrm{S}}
  • 修正串联电阻的影响:用 VJV_{J} 代替 VAV_{A}
I=I0eqVJ/kT=I0eq(VAIRS)/kT{ I} \,=\, {I}_{0} e^{q V_{J} / k T} \,=\, {I}_{0} e^{q ( V_{\mathrm{A}}-I R_{S} ) / k T}
  • 计算:假设已知 VJV_{J} ,可以计算出 II ,之后可以计算出 VAV_{A}
  • 相当于消耗掉了一部分外电压,需要更大的外电压来达到相同的电流。
  • 确定串联电阻:
假设的打破:准中性区域边界处有可能不是小电流注入
  • 正偏时,当电压很大的时候,耗尽层边界处少数载流子的浓度接近杂质浓度,少数载流子和多数载流子的浓度都出现波动,于是产生了反向电流。
  • 室温下,外加电压低于 VbiV_{bi} 零点几个伏特时出现。
  • 示意图:
  • 常常与串联电阻效应混同。
总结

电荷控制模型

  • 描述准中性区内过剩/不足的少数载流子电荷的行为表达式。
  • p+np^+-n 结二极管 nn 型一侧的过剩少数载流子:
QP=qAxn ⁣ ⁣ ⁣Δpn(x,t)dxQ_{\mathrm{P}} \,=\, q A \, \int_{x_{\mathrm{n}}}^{\infty} \! \! \! \Delta p_{\mathrm{n}} ( x, \, t ) \, \, d x
  • 使用少数载流子扩散方程,得到:
dQPdt= iDIFF QPτp{\frac{d Q_{\mathrm{P}}} {d t}}=~ i_{\mathrm{DI F F}}-~ {\frac{Q_{\mathrm{P}}} {\tau_{\mathrm{p}}}}
  • 两种方式可以改变区域的过剩空穴电荷:空穴流入流出/R-G机制。
  • 通过上述方式,我们也可以求出准中性区的电流。考虑稳态下: dQPdt=0\frac{dQ_{P}}{dt}=0 ,所以有:
IDIFF=QPτpI_{\mathrm{D I F F}}={\frac{Q_{\mathrm{P}}} {\tau_{\mathrm{p}}}}
  • 也有
IN(xp)=QNτnI_{N} (-x_{p} )=\frac{-Q_{N}} {\tau_{n}}
  • 我们可以通过估计 QPQ_{P} 的值简化求出 IDIFFI_{DIFF} 的步骤。认为 QPQ_{P} 服从长方形分布,于是有:
QPq(ALP)Δpn(x ⁣= ⁣0)=qALPni2ND(eqVA/kT1)Q_{\mathrm{P}} \, \cong\, q ( A L_{\mathrm{P}} ) \Delta p_{\mathrm{n}} ( x^{\prime} \!=\! 0 ) \,=\, q A L_{\mathrm{P}} {\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, \, 1 ) IDIFF=QPτp=qALPτpni2ND(eqVA/kT1)I_{\mathrm{D I F F}}={\frac{Q_{\mathrm{P}}} {\tau_{\mathrm{p}}}}=\, q A {\frac{L_{\mathrm{P}}} {\tau_{\mathrm{p}}}} \, {\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, \, 1 )

窄基区二极管

  • 仍然是对前述假设的修正。这里我们考虑准中性区域宽度等于或者小于少数载流子扩散长度的二极管。对于 p+np^+ -n 结二极管,有如下的边界条件:
Δpn(x=0)=ni2ND(eqVAkT1)Δpn(x=xc)=0\begin{array}{c} \Delta p_{n} ( x^{\prime}=0 ) = \displaystyle\frac{n_{\mathrm{i}}^{2}}{N_{\mathrm{D}}} \left( e^{\frac{q V_{\mathrm{A}}}{k T}} - 1 \right) \\ \Delta p_{n} ( x^{\prime}=x_{c}^{\prime} ) = 0 \end{array}

其他推导同理想二极管,可以解出:

Δpn(x)=Δpn(0)sinh(xcxLp)sinh(xcLp)\Delta p_{\mathrm{n}} ( x^{\prime} ) = \Delta p_{\mathrm{n}} ( 0 ) \frac{\sinh \left( \frac{x_{\mathrm{c}}^{\prime}-x^{\prime}}{L_{\mathrm{p}}} \right)}{\sinh \left( \frac{x_{\mathrm{c}}^{\prime}}{L_{\mathrm{p}}} \right)} IDIFF=I0(eqVAkT1)I0qADPLPni2NDcosh(xcLP)sinh(xcLP)\begin{array}{l} I_{\mathrm{DIFF}} = I_{0}^{\prime} \left( e^{\frac{q V_{A}}{k T}} - 1 \right) \\ I_{0}^{\prime} \equiv q A \frac{D_{\mathrm{P}}}{L_{\mathrm{P}}} \frac{n_{\mathrm{i}}^{2}}{N_{\mathrm{D}}} \frac{\cosh \left( \frac{x_{\mathrm{c}}^{\prime}}{L_{\mathrm{P}}} \right)}{\sinh \left( \frac{x_{\mathrm{c}}^{\prime}}{L_{\mathrm{P}}} \right)} \end{array}

其中:

sinh(x)=exex2\sinh(x) = \frac{e^{x} - e^{-x}}{2} cosh(x)=ex+ex2\cosh(x) = \frac{e^{x} + e^{-x}}{2} sinh(ξ){ξξ0eξ2ξ\sinh(\xi) \to \begin{cases} \xi & \cdots \xi \to 0 \\ \frac{e^{\xi}}{2} & \cdots \xi \to \infty \end{cases} cosh(ξ){1+ξ22ξ0eξ2ξ\cosh(\xi) \to \begin{cases} 1 + \frac{\xi^{2}}{2} & \cdots \xi \to 0 \\ \frac{e^{\xi}}{2} & \cdots \xi \to \infty \end{cases}
  • 极限情况:
    • xcx_{c}'\to \infty ,该式化简为理想二极管公式
    • xc0x_{c}'\to0 or xcLP1\frac{x_{c}'}{L_{P}}\ll 1 ,可以将该式化为: Δpn(x)=Δpn(0)(1xxc)\Delta p_{\mathrm{n}} ( x^{\prime} ) \,=\, \Delta p_{\mathrm{n}} ( 0 ) \left( 1 \, \,-\, {\frac{x^{\prime}} {x_{\mathrm{c}}^{\prime}}} \right) 这时, Δpn\Delta p_{n} 是一个线性的方程。
      • JpJ_{p} 是一个常数。
      • 忽略了热复合-产生项(因为准中性区太短了,可以忽略不计)
      • I0qA(ni2ND)(DPxc)I_{0}'\to qA\left( \frac{n_{i}^2}{N_{D}} \right)\left( \frac{D_{P}}{x_{c}'} \right) 。反向偏压升高,耗尽层宽度增加,准中性区域的宽度 xc=xcxnx_{c}'=x_c-x_{n} 减小。由上面近似的式子可以知道, IDIFFI_{DIFF} 会随着反向偏压增加而增大。
      • xc0x_{c}'\to{0} 的时候,称其为穿通IDIFF(VA<0)I_{DIFF}(V_{A}<0)\to-\infty 。此时理论公式无效。
  • 推广到双边的窄基区二极管
    • WPW_{P}'WNW_{N}' 分别是两侧准中性区的长度,且 WPLPW_{P}'\ll L_{P}WNLNW'_{N}\ll L_{N} ,可以得到
I=qAni2[DPWNND+DNWPNA](eqVAkT1)=I0(eqVAkT1)I = q A n_{i}^{2} \left[ \frac{D_{P}}{W_{N}^{\prime} N_{D}} + \frac{D_{N}}{W_{P}^{\prime} N_{A}} \right] \left( e^{\frac{q V_{A}}{k T}} - 1 \right) = I_{0} \left( e^{\frac{q V_{A}}{k T}} - 1 \right)
  • 窄基区二极管的电荷控制模型:用少数载流子运输时间 τtr\tau_{tr} 代替 τp\tau_{p}τn\tau_{n}
  • 于是有
IDIFF=QPτPI_{DIFF}=\frac{Q_{P}}{\tau_{P}}

分别计算 QPQ_{P}IPI_{P} ,由于 Δpn\Delta p_{n} 是线性的,所以在这里:

QP=qA ⁣ ⁣ ⁣xnWNΔ ⁣pn(x)dx ⁣= ⁣qA ⁣ ⁣ 12Δ pn(xn)WNQ_{P}=q A \! \! \! \int_{x_{n}}^{W_{N}} \Delta\! p_{n} ( x ) d x \!=\! q A \! \!\ \frac{1} {2} \Delta\ p_{n} ( x_{n} ) W_{N}^{\prime} IP=AJP=qADpdΔ ⁣pndx=qADpΔ ⁣pn(xn)WNI_{P}=A J_{P}=-q A D_{p} \, \frac{d \Delta\! p_{n}} {d x}=q A D_{p} \, \frac{\Delta\! p_{n} ( x_{n} )} {W_{N}^{\prime}}

二者相除:

τtr,p=QPIp ⁣= ⁣(WN)22Dp\tau_{t r, p}=\frac{Q_{P}} {I_{p}} \!=\! \frac{\left( W_{N}^{\prime} \right)^{2}} {2 D_{p}} τtr,n=(WP)22DN\tau_{t r, n}={\frac{( W_{P}^{\prime} )^{2}} {2 D_{N}}}

小信号导纳

  • 小信号:一个正弦电压 vBv_{B} 叠加在外加直流偏压上。
    • 小信号导纳 Y=ivgY=\frac{i}{v_{g}} 表征无源器件的交流响应。
    • 小信号等效电路表示响应情况
    • pnpn 结二极管结区的响应包括电容和电导分量,导纳为 Y=G+jwCY=G+jwCww 为交流信号的角频率)
  • 前提:多数载流子需要快速流入/流出相关区域,也即器件随交流信号准静态变化(载流子可以对交流信号做出及时的响应,就如该信号是直流偏压的情况)。

反向偏置

反向偏置电容

  • 电导很小, Y=jwCY=jwC ;只与器件内的多数载流子振荡有关
  • 小交流信号会导致耗尽层宽度在其稳态值附近进行小幅的摆动,因此有电容的效应。类似于平行板电容器,二极管的电容:
CJ=KSεoAWC_{\mathrm{J}}=\frac{K_{\mathrm{S}} \varepsilon_{\mathrm{o}} A} {W}
  • 反向偏置中 C=CJC=C_{J} ,由于 WW 随反向偏压升高而增加,所以当反向偏压增大的时候 CC 会降低。
  • 我们会发现 WWVAV_{A} 的关系随杂质分布而变化,这非常不方便,我们希望获得一个通用的表达式。由于这种变化主要来自杂质浓度分布的不同,因此我们用一个公式来表征上述分布 NB(x)=bxmx>0N_{B}(x)=bx^m \dots x>0 (其中 b>0b>0 ,且 mm 为常数)(这个式子表征了单边结的指数分布,也即一边重掺杂,一边轻掺杂的杂质浓度)
    • m=0m=0 对应非对称突变结
    • m=1m=1 对应单边线性缓变结
    • m<0m<0 称为超突变分布(超突变结)
  • 代入后, m>2m>-2 时耗尽层宽度关系式:
W=[(m  +  2)KSε0qb(Vbi    VA)]1/(m+2)W=\left[ \frac{( m \;+\; 2 ) K_{\mathrm{S}} \varepsilon_{0}} {q b} ( V_{\mathrm{b i}} \;-\; V_{\mathrm{A}} ) \right]^{1 / ( m+2 )} CJ=KSε0A[(m + 2)KSε0qb(Vbi  VA)]1/(m+2)C_{J}=\frac{K_{\mathrm{S}} \varepsilon_{0} A} {\left[ \frac{( m \ +\ 2 ) K_{\mathrm{S}} \varepsilon_{0}} {q b} ( V_{\mathrm{b i}} \ -\ V_{\mathrm{A}} ) \right]^{1 / ( m+2 )}}

也可以用 CJ0C_{J_{0}}CJC_{J} 来描述上述电容关系

CJ0CJVA=0=KSε0A[(m+2)KSε0qbVbi]1/(m+2)C_{J 0} \, \equiv\, C_{J} \big|_{V_{\mathrm{A}}=0} \,=\, {\frac{K_{S} \varepsilon_{0} A} {\left[ {\frac{( m \,+\, 2 ) K_{S} \varepsilon_{0}} {q b}} \, V_{\mathrm{b i}} \right]^{1 / ( m+2 )}}} CJ=CJ0(1  VAVbi)1/(m+2)C_{\mathrm{J}}={\frac{C_{\mathrm{J 0}}} {\left( 1 ~-~ {\frac{V_{\mathrm{A}}} {V_{\mathrm{b i}}}} \right)^{1 / ( m+2 )}}}
  • 有上述关系的 pnpn 二极管称为变容二极管。该二极管的一个参数称为调制率 TRTR ,表征给定电压范围内二极管表现出的最大的电容比。
TR=C1(VA1)C1(VA2)(VA2VA1)1/(m+2)T R=\frac{C_{1} ( V_{\mathrm{A 1}} )} {C_{1} ( V_{\mathrm{A 2}} )} \cong\left( \frac{V_{\mathrm{A 2}}} {V_{\mathrm{A 1}}} \right)^{1 / ( m+2 )}

(反偏时, VAVbiV_{A}\gg V_{bi}mm 越小, TRTR 越大。

  • 计算相关参数:假设测试器件是非对称掺杂的突变结( m0m\to{0}bNBb\to N_{B} ),采用上述电容公式并近似,有:
1CJ2=2qNBKSε0A2(Vbi  VA)\frac{1} {C_{\mathrm{J}}^{2}}=\frac{2} {q N_{\mathrm{B}} K_{\mathrm{S}} \varepsilon_{0} A^{2}} ( V_{\mathrm{b i}} ~-~ V_{\mathrm{A}} )

通过斜率和截距的关系,可以得到 NBN_{B}VbiV_{bi} 等的数值。

反向偏置电导

  • 假设二极管准静态地响应一个交流信号,那么交流电导就等于微分直流电导。
G0=dIdVAG_{0}=\frac{d I} {d V_{\mathrm{A}}}
  • 理想二极管的式子 I=I0(eqVA/kT1)I=I_{0}(e^{ qV_{A}/kT }-1)VAV_{A} 求导。
G0=qkTI0eqVA/kT=qkT(I+ I0)G_{0}={\frac{q} {k T}} I_{0} e^{q V_{A} / k T}={\frac{q} {k T}} ( I+\ I_{0} )
  • 如果考虑复合-产生电流占主导地位(参见前面的修正)
G0=ddVA(qAni2τ0W)=qAniW/2τ0(m+2)(VbiVA)G_{0} \,=\, {\frac{d} {d V_{\mathrm{A}}}} \bigg(-\, {\frac{q A n_{\mathrm{i}}} {2 \tau_{0}}} W \bigg) \,=\, {\frac{q A n_{\mathrm{i}} W / 2 \tau_{0}} {( m \,+\, 2 ) ( V_{\mathrm{b i}} \,-\, V_{\mathrm{A}} )}}

正向偏置

  • 反向偏置导纳分析中引入结电容是多数载流子进入和离开稳态耗尽层引起的,此时少数载流子浓度的涨落对导纳的影响忽略不计。但在正向偏置中,交流信号引起的少数载流子电荷也会对结电容有贡献。
  • 与多数载流子不同,当角频率接近少数载流子寿命的倒数时,少数载流子电荷的涨落很难与交流信号保持同步,而是异步的变化——增加测量电导,减小测量电容,并且依赖频率。
  • 少数载流子在耗尽层边界的积累由扩散电流引起,导纳称为扩散导纳。
YD = GD + jωCD Y_{\mathrm{D}} ~=~ G_{\mathrm{D}} ~+~ j \omega C_{\mathrm{D}} ~

计算过程

  • p+np^+ -n 结二极管:
  • 要考虑少数载流子,就要利用少数载流子扩散方程
Δpn(x,t)t=Dp2Δpn(x,t)x2Δpn(x,t)τp\frac{\partial\Delta p_{\mathrm{n}} ( x, t )} {\partial t} \,=\, D_{\mathrm{p}} \, \frac{\partial^{2} \Delta p_{\mathrm{n}} ( x, t )} {\partial x^{2}} \,-\, \frac{\Delta p_{\mathrm{n}} ( x, t )} {\tau_{\mathrm{p}}}
  • 交流信号带来的 Δpn(x,t)\Delta p_{n}(x,t) 的改变:
Δpn(x,t)=Δpn(x)+p~n(x,ω)ejωt\Delta\, p_{\mathrm{n}} ( x, t ) \,=\, \overline{{{\Delta\, p_{\mathrm{n}}}}} ( x ) \,+\, \tilde{p}_{\mathrm{n}} ( x, \omega) e^{j \omega t}
  • 分别计算直流项和交流项:
0=Dpd2Δpndx2Δpnτp0 \,=\, D_{\mathrm{p}} \, \frac{d^{2} \overline{{{\Delta p_{\mathrm{n}}}}}} {d x^{2}} \,-\, \frac{\overline{{{\Delta p_{\mathrm{n}}}}}} {\tau_{\mathrm{p}}} 0=Dpd2p~ndx2p~nτp/(1+jωτp)0 \,=\, D_{\mathrm{p}} \, \frac{d^{2} \tilde{p}_{\mathrm{n}}} {d x^{2}} \,-\, \frac{\tilde{p}_{\mathrm{n}}} {\tau_{\mathrm{p}} / ( 1 \,+\, j \omega\tau_{\mathrm{p}} )}
  • 边界条件:
Δpn()=p~n()=0\overline{{{\Delta p_{\mathrm{n}}}}} ( \infty)=\widetilde{p}_{\mathrm{n}} \left( \infty\right)=0

直接把交流的电压代入前面推出的边界处的过剩载流子公式:

Δpn(x ⁣= ⁣xn)=Δpn(xn)+p~n(xn)=n12ND(eq(VA+vn)/kT1)\Delta p_{n} ( x \!=\! x_{n} ) \,=\, \overline{{ {\Delta p_{n}} }} ( x_{n} ) \,+\, \tilde{p}_{n} ( x_{n} ) \,=\, \frac{n_{1}^{2}} {N_{\mathrm{D}}} ( e^{q ( V_{\mathrm{A}}+v_{n} ) / k T} \,-\, \, 1 )

已知直流形式与之前推出的相同:

Δpn(xn)=ni2ND(eqVA/kT1)\overline{{ {\Delta p_{n}} }} ( x_{n} ) \,=\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}} ( e^{q V_{\mathrm{A}} / k T} \,-\, 1 )

和上式相减,在采用无穷小时的近似:

p~n(xn)=ni2NDeqVA/kT(eqva/kT1)ni2ND ⁣ ⁣(qvakTeqVA/kT)vakT/q\tilde{p}_{\mathrm{n}} ( x_{\mathrm{n}} ) \,=\, \frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}} e^{q V_{\mathrm{A}} / k T} ( e^{q v_{\mathrm{a}} / k T} \,-\, \, 1 )\cong\frac{n_{\mathrm{i}}^{2}} {N_{\mathrm{D}}} \! \! \left( \frac{q v_{a}} {k T} e^{q V_{\mathrm{A}} / k T} \right) \qquad\dots v_{a} \ll k T / q
  • 利用边界条件,以及 τpτp1+jwτp\tau_{p}\to \frac{\tau_{p}}{1+jw\tau_{p}} 的变化进行修正。得到:
idiff=G01+jωτpvai_{\mathrm{d i f f}} \,=\, G_{0} \sqrt{1 \,+\, j \omega\tau_{\mathrm{p}}} \, v_{\mathrm{a}}
  • 正向偏置时的电导:
YD=idiffυa=G01+jωτpY_{\mathrm{D}} \,=\, {\frac{i_{\mathrm{d i f f}}} {\upsilon_{\mathrm{a}}}} \,=\, G_{0} \sqrt{1 \,+\, j \omega\tau_{\mathrm{p}}}
  • 如果是 n+pn^+ -p 结的二极管: τp\tau_{p} 换成 τn\tau_{n}
  • 分解实虚部
    • 扩散电导(实部) GD=G02(1+ω2τp2  +  1).1/2G_{\mathrm{D}}={\frac{G_{0}} {\sqrt{2}}} \left( \sqrt{1+\omega^{2} \tau_{\mathrm{p}}^{2}} \;+\; 1 \right)_{.}^{1 / 2}
    • 扩散电容(虚部) CD=G0ω2(1+ω2τp21)1/2C_{\mathrm{D}}=\frac{G_{0}} {\omega\sqrt{2}} \left( \sqrt{1 \,+\, \omega^{2} \tau_{\mathrm{p}}^{2}} \,-\, 1 \right)^{1 / 2}
  • 低频( wτp1w \tau_{p}\ll 11+w2τp21+w2τp22\sqrt{ 1+w^2\tau^2_{p} } \cong 1+w^2 \frac{\tau_{p}^2}{2}
GDG0G_{D}\to G_0 CDG0τp2C_{D}\to \frac{G_{0 }\tau _{p}}{2}

瞬态响应

  • 关心 pnpn 结二极管在开关过程中达到的速度,在瞬态关断过程中,速度受限最为显著。
  • 前提:与二极管的最大正向压降 VONV_{ON} 相比,电源电压(无论正偏还是反偏)足够大。

关断特性

电路特性

  • 什么是瞬态关断: t=0t=0 之前在二极管两端加正向偏置电压并达到稳态, t=0t=0 时刻开关打向反偏。
  • 电流特性:开关转换后的瞬间,反向电流的大小与正向电流相当,在衰减到稳态值之前的时间内,基本保持常数。
    • 存贮延迟时间 tst_{s} :反向电流保持为常数的时间
    • 反向恢复时间 trrt_{rr} :反向电流衰减到它最大数值 10%10\% 所需的时间
    • 恢复时间 trt_{r} :反向电流衰减的时间,前两者之差。
  • 瞬时电压变化:
    • 在存贮延迟时间仍保持正偏,在 t=tst=t_{s} 时等于 00

物理机制

  • 二极管正向偏置时,使紧邻耗尽区的准中性区内产生过剩少数载流子的积累;反偏时,紧邻耗尽区的少数载流子显著减少。过剩少数载流子需要一定的时间从准中性区移走,这时间就是存贮延迟时间。
  • 由电荷控制模型可知,由两种方法可以移走过剩少数载流子的电荷:
    • 复合:如果是移走载流子的唯一方式,那么从开态到关态的时间与少数载流子的寿命相同。
    • 净载流子漂移:流回结的另一边,非常快。此时二极管几乎等效为短路,最大反向电流由外电压 VRV_{R} 与外电阻(保护电阻) RRR_{R} 之比 IRI_{R} 决定。
    • 以上两种机理中没有一种可以以足够快的速度把电荷完全消除,因此会观察到延迟时间。
  • 0ts0\sim t_{s} 时,二极管如何保持正偏:对 p+np^+ -n 结二极管的过剩少数载流子而言,它与二极管两端的电压有如下的关系(这里我们用表示交流电压的 vAv_{A} 代替 VAV_{A}
Δpn(xn)=ni2ND(eqvA/kT1)\Delta p_{\mathrm{{n}}} ( x_{\mathrm{{n}}} ) \,=\, \frac{n_{\mathrm{{i}}}^{2}} {N_{\mathrm{{D}}}} ( e^{q v_{\mathrm{A}} / k T} \,-\, 1 )

从这个式子中我们可以观察到:当 vAv_{A} 大于 00 时,紧邻耗尽层边界的区域中有过剩的少数载流子。实际上,是耗尽区边界与内部过剩载流子的存在才使结保持正偏。 考虑公式: i=AJP(xn)=qADPdΔpndxx=x0i=A J_{P}(x_{n})=-qA D_{P} \frac{d\Delta p_{n}}{dx}\big|_{x=x_{0}} ,可知:

dΔpndxx=xn=iqADP\frac{d\Delta p_{n}}{dx}\big|_{x=x_{n}}=-\frac{i}{qAD_{P}}

这也是 pn(x,t)xp_{n}(x,t)-x 图上 x=xnx=x_{n} 处曲线的斜率。由于 0ts0\sim t_{s} 时, i=IRi=-I_{R} 为常数,所以对于所有 t>0t>0 的曲线,在 x=xnx=x_{n} 处的斜率必定相同。

存贮延迟时间

  • 推算存贮延迟时间 tst_{s} 的定量关系式:考虑 p+np^+ -n 二极管,使用电荷控制模型。(由于是单边结,相比 nn 型一侧的空穴存贮电荷 QPQ_{P} 而言,存贮于 pp 型一侧的电子电荷可以忽略;理想二极管,忽略耗尽区的R-G过程,故而 i=iDIFFi=i_{DIFF}
  • 0ts0\sim t_{s}i=IRi=-I_{R} 为常数,代入电荷控制模型:
dQPdt=  (IR+QPτp)0+ tts\frac{d Q_{\mathrm{P}}} {d t}=\ -\ \left( I_{\mathrm{R}} \,+\, \frac{Q_{\mathrm{P}}} {\tau_{\mathrm{p}}} \right) \quad\ldots\, 0^{+} \leqslant\ t \, \leqslant\, t_{\mathrm{s}}
  • 分离 QPQ_{P}tt ,从 t=0+t=0^+t=tst=t_{s} 积分,得到
ts=τpln(IR+QPτP)QP(0+)QP(ts)=τpln[IR+QP(0+)τpIR+QP(ts)τp]t_{s} = -\tau_{p} \ln \left( I_{R} + \frac{Q_{P}}{\tau_{P}} \right) \bigg|^{Q_{P}(t_{s})}_{Q_{P}(0^+)} = \tau_{p} \ln \left[ \frac{I_{R} + \frac{Q_{P}(0^+)}{\tau_{p}}}{I_{R} + \frac{Q_{P}(t_{s})}{\tau_{p}}} \right]
  • 边界条件:
QP(0+)=QP(0)Q_{P}(0^+)=Q_{P}(0^-)

关断时: dQPdt=0\frac{dQ_{P}}{dt}=0i=IFi=I_{F} 。代回前面列出的电荷控制模型方程,得到:

IF=Qp(0)τp=Qp(0+)τpI_{\mathrm{F}}={\frac{Q_{\mathrm{p}} ( 0^{-} )} {\tau_{\mathrm{p}}}}={\frac{Q_{\mathrm{p}} ( 0^{+} )} {\tau_{\mathrm{p}}}}

同时假设 QP(ts)Q_{P}(t_{s}) 近似为 00 。将这些条件代入前面得到的有关 tst_{s} 的式子中:

ts=τpln(1+IFIR)t_{\mathrm{s}} \,=\, \tau_{\mathrm{p}} \, \operatorname{l n} \Biggl( 1 \,+\, \frac{I_{\mathrm{F}}} {I_{\mathrm{R}}} \Biggr)
  • 从公式以及理解的角度,可以得到(当其他变量不变时):
    • IFI_{F} 变大,使得二极管的存贮电荷增加,需要更长的时间来移走存贮电荷。
    • IRI_{R} 变大,意味着反向电流移走存贮电荷的速率增加,移走存贮电荷需要更短的时间。
    • τP\tau_{P} 减少,则少数载流子的复合速率增加,因此 tst_{s} 会减少。
  • 阶跃恢复二极管 trt_{r} 不分很短,反向电流在达到 t=tst=t_{s} 之后快速地返回稳态值。

瞬态开启

  • 什么是瞬态开启:二极管从反偏的关态切换到正偏的开态。可以通过电流脉冲或电压脉冲或两种脉冲的组合完成,这里只考虑电流脉冲。

电路特性

  • 维持反偏电压的少子迅速注入耗尽区,多子迅速重分布,结电压迅速从反偏移动到 0V0V 。之后电压缓慢从 0V0V 上升到 VONV_{ON}

分析

  • 认为它是理想的 p+np^+ -n 突变结二极管,并采用电荷控制方法。显示预测的 nn 型一侧空穴存贮电荷随时间的增加。
  • 边界条件
    • t=0t=0QP=0Q_{P}=0
    • 瞬态开启过程中 i=IFi=I_{F}
  • 电荷控制方法
dQPdt = IF  QPτp{\frac{d Q_{\mathrm{P}}} {d t}} ~=~ I_{\mathrm{F}} ~-~ {\frac{Q_{\mathrm{P}}} {\tau_{\mathrm{p}}}}
  • 分离变量并积分
t=τpln(IFQPτp)0QP(t)=τpln[1Qp(t)IFτp]t \,=\,-\tau_{\mathrm{p}} \, \ln\Biggl( I_{\mathrm{F}} \,-\, \frac{Q_{\mathrm{P}}} {\tau_{\mathrm{p}}} \Biggr) \Biggr|_{0}^{Q_{\mathrm{P}} ( t )} \,=\,-\tau_{\mathrm{p}} \, \ln\Biggl[ 1 \, \,-\, \frac{Q_{\mathrm{p}} ( t )} {I_{\mathrm{F}} \tau_{\mathrm{p}}} \Biggr]
  • 求解 QPQ_{P}
Qp(t)=IFτp(1et/τp)Q_{\mathrm{p}} ( t ) \,=\, I_{\mathrm{F}} \tau_{\mathrm{p}} ( 1 \, \,-\, \, e^{\,-t / \tau_{\mathrm{p}}} )
  • 稳态条件下,空穴存贮电荷(注意跟上面的不一样)
QP=IDIFFτp=I0τp(eqVA/kT1)Q_{\mathrm{P}} \,=\, I_{\mathrm{D I F F}} \tau_{\mathrm{p}} \,=\, I_{0} \tau_{\mathrm{p}} ( e^{q V_{A} / k T} \,-\, 1 )

由于是瞬态过程,用 vAv_{A} 表示 QP(t)Q_{P}(t) ,与上式联立,得到:

υA(t)=kTqln[1+IFI0(1et/τp)]\upsilon_{\mathrm{A}} ( t ) \,=\, \frac{k T} {q} \, \ln\left[ 1 \,+\, \frac{I_{\mathrm{F}}} {I_{0}} ( 1 \,-\, \, e^{\,-t / \tau_{\mathrm{p}}} ) \right]
  • 结论
    • QP()=IFτPQ_{P}(\infty)=I_{F}\tau_{P}
    • 瞬态开启的开始阶段有一个快速的上升,后期变得缓慢。

金属半导体接触

  • 金属半导体从物理上和功能上,都与非对称的突变二极管 p+np^+-nn+pn^+ -p 很相似。

MS分析

  • 理想的MS接触
    • 金属与半导体在原子尺度上紧密地接触,两者之间不存在任何夹层
    • 金属和半导体之间不存在相互扩散或混合
    • MS界面没有吸附的杂质或表面电荷。
  • 能带图
  • E0E_{0} 真空能级:电子完全脱离材料本身必须具有的最小能量
  • Φ\Phi 功函数:真空能级与费米能级之间的能量差。
    • ΦM\Phi_{M} 金属功函数,不变的基本参数,接触形成过程中保持不变
    • ΦS\Phi_{S} 半导体功函数,公式:ΦS=χ+(EcEF)FB\Phi_{\mathrm{S}} \,=\, \chi\,+\, ( E_{\mathrm{c}} \,-\, E_{\mathrm{F}} )_{\mathrm{F B}}
      • χ\chi 电子亲和势, χ=(E0Ec)surface\chi=(E_{0}-E_{c})|_{surface} ,对于一定的半导体材料是恒定的基本参数,在接触形成过程中保持不变
      • (EcEF)FB(E_{c}-E_{F})_{FB} 是在平带和零电场条件下 EcE_{c}EFE_{F} 之间的能量差。
        • (EcEF)FB=EG2+EFEi(E_{c}-E_{F})_{FB}=\frac{E_{G}}{2}+E_{F}-E_{i}
        • EF    Ei=  kTln(n/ni)  =  kTln(p/ni)E_{\mathrm{F}} \;-\; E_{\mathrm{i}} \,=\; k T \, \ln( n / n_{\mathrm{i}} ) \;=\;-\, k T \, \ln( p / n_{\mathrm{i}} )
  • MS接触的过程
    • 接触形成之后的瞬间:金属和半导体的能带图以公共的 E0E_{0} 参考能级为基准垂直地对齐,并在界面处对接。此时,由于费米能级不相等,故而MS接触明显处于非平衡态。
    • 之后,载流子会根据半导体的类型与势垒的大小流动,达到平衡态。

n型半导体

ΦM>ΦS\Phi_{M}>\Phi_{S} 的情况
  • 刚接触时:
  • 之后,电子从半导体向金属内转移(因为 nn 型半导体的费米能级高于金属的 EFME_{FM} ),而电子损失导致电子势垒的增加,阻止电子从半导体向金属内转移,直到电子的转移和势垒的阻止互相抵消,此时 EFE_{F} 在整个结构中都一样。
  • 注意在边界处的能带不发生变化,但在无限远处 EFE_{F} 与金属的 EFME_{FM} 齐平。其中 ΦB=ΦMχ\Phi_{B}=\Phi_{M}-\chi 是在金属中具有 E=EFE=E_{F} 能量的电子遇到的表面势垒。
  • 肖特基势垒实际上既阻止电子从金属到半导体,也阻止电子从半导体到金属,只不过后者更容易一些,而且会受到外加电压的影响。
添加偏压 VA>0V_{A}>0
  • EFM<EFSE_{FM}<E_{FS} ,对半导体中的电子而言势垒降低,有电流从半导体流向金属。 VAV_{A} 增加,正向偏置电流快速上升,可以越过表面势垒的电子数量呈指数增加。
  • 示意图:
添加偏压 VA<0V_{A}<0
  • EFM>EFSE_{FM}>E_{FS} ,阻止电子从半导体向金属的流动;同时金属中只有很少的电子可以越过 ΦB\Phi_{B} 势垒,反向偏置电流相对较小。
  • 示意图:
总结
  • 类似于 pnpn 结二极管,具有整流特性or整流接触。
  • 金属与 pp 型半导体接触,其中 ΦM<ΦS\Phi_{M}<\Phi_{S} 的情况与这种情况类似,示意图:
ΦM<ΦS\Phi_{M}<\Phi_{S} 的情况
  • 刚接触时:
  • 之后,由于 EFME_{FM}EFSE_{FS} 之间的差距比较小,所以电子在半导体和金属之间自由流动,达到平衡状态。
  • 这里的 ΦB\Phi_{B} 不是特别大。
添加偏压 VA>0V_{A}>0
  • 没有任何势垒,很小的电压就会造成很大的正向偏置电流。
添加偏压 VA<0V_{A}<0
  • 从金属到半导体有一个很小的势垒
  • 反偏超过零点几个伏特时,势垒就变为 00
  • 较小的反偏电压下就有比较大的反偏电流,而且反偏电流不会饱和。
总结
  • 非整流型的/欧姆型的接触。
  • 金属与 pp 型半导体接触时,其中 ΦM>ΦS\Phi_{M}>\Phi_{S} 的情况与这种情况类似,示意图:

肖特基二极管

静电特性

  • 回顾一下内建电压的概念:在分析理想 pnpn 结的时候,内建电压来源于势垒,也就是当两个半导体拼在一起的时候,它们费米能级的差值再除以 qq 。ok,我们把这个思路迁移到整流MS二极管上:
  • ΦBn\Phi_{Bn} 是肖特基势垒。
  • 从这个图中我们把金属和半导体拼在一起,明显看到 EFME_{FM}EFSE_{FS} 之间的差值是 ΦMΦS\Phi_{M}-\Phi_{S} 。但 ΦS\Phi_{S} 实际上是算出来的,它等于 χ+(ECEFS)FB\chi+(E_{C}-E_{FS})_{FB} ,我们把这个式子代入,并且用 ΦB\Phi_{B} 代替 ΦMχ\Phi_{M}-\chi ,所以就可以得到:
Vbi=1q[ΦB(EcEF)FB]V_{\mathrm{b i}} \,=\, {\frac{1} {q}} \left[ \Phi_{\mathrm{B}} \,-\, \left( E_{\mathrm{c}} \,-\, E_{\mathrm{F}} \right)_{\mathrm{F B}} \right]
  • 考虑载流子堆积的情况。对 nn 型MS接触(整流),由于半导体中的电子流向金属,半导体中会出现一个耗尽区,宽度为 xnx_{n} ,耗尽区中具有施主离子产生的净正电荷。不过,在金属一侧,不会出现负的受主电荷了。所以,MS接触中,正电荷以如下方式进行堆积
ρ{qND0xW0x>W\rho \cong \begin{cases} q N_{\mathrm{D}} & \quad \dots 0 \leqslant x \leqslant W \\ 0 & \quad \dots x > W \end{cases}

在交界处,负电荷以冲激函数的形式存在,半导体一侧的空间正电荷把金属中的空间负电荷吸引到结区的边缘,且金属中的自由负电荷非常多,因此形成了冲激函数。

  • 然后是一套丝滑小连招:泊松方程:
dEdx=ρKsε0qNDKsε0 0xW\frac{d {\mathcal E}} {d x}=\frac{\rho} {K_{s} \varepsilon_{0}} \cong\frac{q N_{\mathrm{D}}} {K_{s} \varepsilon_{0}} \qquad\dots \ 0 \, \leqslant\, x \, \leqslant\, W

得到:

E(x)=qNDKSε0  (Wx)      .  .  .  0xW{\mathcal E} ( x ) \,=\,-\, \frac{q N_{\mathrm{D}}} {K_{S} \varepsilon_{0}} \; ( W \,-\, x ) \; \; \;. \;. \;. \; 0 \, \leqslant\, x \, \leqslant\, W
  • 继续,静电势积一下
V(x)=qND2KSε0(Wx)20xWV ( x ) \,=\,-\, \frac{q N_{\mathrm{D}}} {2 K_{\mathrm{S}} \varepsilon_{0}} \, ( W \,-\, x )^{2} \quad\dots \, \, 0 \, \leqslant\, x \, \leqslant\, W
  • 等效于 pnpn 结分析中的结果,只是电势不同(在这里,我们为了分析方便,将电势参考点选在半导体表面,方便对半导体进行分析),因此在 x=0x=0 处,电势为 Vbi-V_{bi} (此处刚好是金属的电势,但我们设半导体的电势为 00 点,同时要考虑电子电荷为 q-q 对能带的影响),加入并考虑外加电场 VAV_{A} ,可以假设 x=0x=0V=(VbiVA)V=-(V_{bi}-V_{A})
  • 代入上述 VV 可以解出耗尽区宽度
W=[2KSε0qND(Vbi  VA)]1/2W=\left[ {\frac{2 K_{\mathrm{S}} \varepsilon_{0}} {q N_{\mathrm{D}}}} \left( V_{\mathrm{b i}} \ -\ V_{\mathrm{A}} \right) \right]^{1 / 2}

I-V特性

  • 对MS二极管( nn 型)来说,正偏时电流主要来源于半导体进入金属电子注入(势垒非常低),耗尽层中的复合电流与空穴扩散可以忽略(我很好奇一个金属哪来的空穴),示意图
  • 反偏时电流主要来自金属流向半导体的电子流,反偏空穴扩散电流与耗尽区的载流子相关的R-G电流一般可以忽略。
  • MS二极管称为“多数载流子器件”
正偏
  • 正偏时,电子越过势垒进行注入,可看做热电子发射电流
KEx=12 mnυx2q(Vbi  VA)\mathrm{K E}_{\mathrm{x}}=\frac{1} {2} ~ m_{\mathrm{n}}^{\star} \upsilon_{\mathrm{x}}^{2} \geq q ( V_{\mathrm{b i}} ~-~ V_{\mathrm{A}} )
  • 可以通过解出电子速度 vxv_{x} 从而得到与电子相关的电流
ISM,νx  =  qAυxn(υx)I_{\mathrm{S \bullet} \to{\bf M}, \nu_{x}} \;=\;-\, q A \upsilon_{x} n ( \upsilon_{x} )

(此处省略一个 n(vx)n(v_{x}) 的推导式子,感觉贴了也没啥用,想看看书P352页罢,半导体器件基础那本(雾),总而言之,可以简化得到:

ISM=AAT2eΦB/kTeqVA/kTI_{\mathbf{S \bullet} \to \mathbf{M}} = \mathcal{A} \mathcal{A}^{*} T^{2} e^{-\Phi_{\mathbf{B}} / k T} e^{q V_{\mathbf{A}} / k T} A(mnm0)A\mathcal{A}^{*} \equiv \left( \frac{m_{\mathrm{n}}^{*}}{m_{\mathrm{0}}^{*}} \right) \mathcal{A}

理查森常数:

A4πqm0k2h3=120amps/(cm2K2)\mathcal{A} \equiv \frac{4 \pi q m_{0} k^{2}}{h^{3}} = 120 \, \mathrm{amps} / \mathrm{(cm^{2} \cdot K^{2})}

电子从金属到半导体,穿过界面时遇到的势垒高度 ΦB\Phi_{B} 始终不变,并且基本不受到外加电压的影响,于是有:

IMS(VA)=IMS(VA=0)I_{\mathrm{M \bullet} \to \mathrm{S}} (V_{\mathrm{A}}) = I_{\mathrm{M \bullet} \to \mathrm{S}} (V_{\mathrm{A}} = 0)

且,平衡时无电流,所以可以得出:

IMS(VA=0)=ISM(VA=0)=AAT2eΦB/kTI_{\mathrm{M \bullet} \to \mathrm{S}} (V_{\mathrm{A}} = 0) = -I_{\mathrm{S \bullet} \to \mathrm{M}} (V_{\mathrm{A}} = 0) = -\mathcal{A} \mathcal{A}^{*} T^{2} e^{-\Phi_{\mathrm{B}} / k T}
  • 总电流
I=ISM+IMS=ISM+IMS(VA=0)I = I_{\mathbf{S \bullet} \to \mathbf{M}} + I_{\mathbf{M \bullet} \to \mathbf{S}} = I_{\mathbf{S \bullet} \to \mathbf{M}} + I_{\mathbf{M \bullet} \to \mathbf{S}} (V_{\mathbf{A}} = 0)
  • 代入上面的公式:(IsI_{s} 比较大)
I=Is(eqVA/kT1)I = I_{s} \left( e^{q V_{A} / k T} - 1 \right) Is=AAT2eΦB/kTI_{s} = \mathcal{A} \mathcal{A}^{*} T^{2} e^{-\Phi_{B} / k T}

交流响应

  • 交流小信号加在直流反偏电压上,MS界面处 δ\delta 函数电荷发生波动,多数载流子从半导体中快速地移进移出,耗尽区宽度也发生波动。示意图
  • 电容类似于突变结且一边重掺杂的半导体的电容,有:
1C2=2qNDKSε0A2(VbiVA)\frac{1}{C^{2}} = \frac{2}{q N_{\mathrm{D}} K_{\mathrm{S}} \varepsilon_{0} A^{2}} (V_{\mathrm{bi}} - V_{\mathrm{A}})

如果 VbiV_{\mathrm{bi}}NDN_{\mathrm{D}} 已知,即可算出 ΦBn\Phi_{\mathrm{Bn}}

qVbi=ΦBn(EcEF)FB=ΦBnkTln(NcND)q V_{\mathrm{bi}} = \Phi_{\mathrm{Bn}} - (E_{\mathrm{c}} - E_{\mathrm{F}})_{\mathrm{FB}} = \Phi_{\mathrm{Bn}} - k T \ln \left( \frac{N_{\mathrm{c}}}{N_{\mathrm{D}}} \right)

欧姆接触

  • 由于费米能级钉扎效应的存在,很难满足欧姆接触的条件。
  • 实际上,要达到欧姆接触,需要高掺杂率的半导体,这样载流子可以隧穿过势垒。
  • 比接触电阻率(每单位面积的接触电阻)
R=ρcAR = \frac{\rho_{c}}{A}

就实际的欧姆接触而言:

ρceHΦB/ND\rho_{c} \propto e^{H \Phi_{B} / \sqrt{N_{\mathrm{D}}}}
  • 理想情况下欧姆接触电阻很小,压降可以忽略不计。