费米能级理论计算

一、项目背景与科学问题

费米能级是半导体物理中决定载流子统计分布与器件电学行为的核心概念。在金属中费米能级位于导带内，是电子占据与未占据的分界线；在半导体中它落在带隙内，其精确位置由掺杂浓度与温度共同决定。费米能级的位置直接控制着pn结势垒高度、MOSFET阈值电压、欧姆接触的形成以及热平衡下载流子浓度的数量级。本项目的科学问题是：硅在300 K下的本征费米能级位于带隙何处？当掺入施主或受主杂质后，费米能级如何随掺杂浓度移动？通过费米-狄拉克统计与质量作用定律，可以对这些工程关键量进行精确计算。

二、计算方法与理论背景

对非简并半导体，导带电子浓度与价带空穴浓度分别表示为n=Nc·exp[-(Ec-EF)/kT]与p=Nv·exp[-(EF-Ev)/kT]，其中Nc与Nv为导带与价带有效态密度。由电中性条件n=p=ni可得本征费米能级Ei=(Ec+Ev)/2+(kT/2)·ln(Nv/Nc)，即本征费米能级并不严格位于带隙中央，而是略偏向有效态密度较小的一侧。掺杂时电中性条件变为n+N_a=p+N_d，在非简并近似下n型半导体费米能级移动为EF-Ei=kT·ln(Nd/ni)，p型则为EF-Ei=-kT·ln(Na/ni)。当掺杂浓度接近或超过有效态密度时进入简并状态，费米能级进入导带或价带，需改用完整的费米积分描述。

三、计算设置与模型构建

选取硅为研究对象，300 K下取带隙Eg=1.12 eV，导带有效态密度Nc=2.8E19 cm⁻³，价带有效态密度Nv=1.04E19 cm⁻³。首先根据上式求解本征载流子浓度ni与Ei的位置。随后在n型一侧使掺杂浓度Nd从1E13 cm⁻³按对数间隔增加到1E20 cm⁻³，在p型一侧同步扫描Na，计算EF-Ei随掺杂的变化并绘制半对数曲线。为完整呈现费米能级的统计含义，还计算在固定温度下电子与空穴浓度随费米能级位置的变化关系，标出本征点与载流子浓度相等的交叉位置。

四、结果分析与案例图解读

图5的左侧面板（a）给出300 K下电子浓度n与空穴浓度p随费米能级位置的半对数曲线。蓝色实线为电子浓度，当费米能级靠近导带底时n急剧上升；红色虚线为空穴浓度，当费米能级靠近价带顶时p急剧上升。两条曲线在带隙内某处相交，交点即本征费米能级Ei，对应载流子浓度约为1E10 cm⁻³量级。由于Nv略小于Nc，Ei略偏向价带一侧，图中绿色点线标出Ei的位置。此图直观说明费米能级在带隙内移动约一个kT，载流子浓度便改变约一个数量级，体现了费米能级对载流子占据的指数级调控能力。

图5的右侧面板（b）给出EF-Ei随掺杂浓度的半对数曲线。蓝色实线为n型，费米能级随掺杂浓度对数增大而上移；红色虚线为p型，随掺杂浓度增大而下移。在低掺杂区，曲线保持线性，斜率为kT，对应非简并状态下的经典关系EF-Ei=kT·ln(Nd/ni)；当掺杂浓度升高至约1E19 cm⁻³、接近导带有效态密度Nc时，费米能级趋近导带底（图中灰色点线标出Ec-Ei=0.56 eV），曲线出现饱和迹象，进入简并状态。此时经典对数近似失效，掺杂浓度的继续增加对费米能级位置的提升作用减弱，半导体呈现类金属导电特性。该曲线对设计掺杂工艺、预估接触势垒与器件阈值具有直接参考价值。

两面板结合，左面板从统计角度揭示了费米能级对载流子浓度的控制，右面板给出了费米能级随掺杂的定量演化，共同完成费米能级理论计算的完整分析。

五、结论与展望

本项目基于费米-狄拉克统计与质量作用定律计算了硅中费米能级的位置及其随掺杂的移动规律，300 K下本征费米能级略偏价带一侧，n型与p型掺杂分别使EF按kT·ln(Nd/ni)上移与下移，高掺杂时进入简并状态。案例图以载流子浓度与EF-Ei双面板呈现。后续可拓展：(1)考虑温度依赖的Nc、Nv与带隙，计算宽温区费米能级；(2)引入简并情形的费米积分求解；(3)结合泊松方程求解实际器件中的费米能级空间分布，支撑器件仿真。