半导体器件建模仿真-边界条件

边界条件

• 欧姆接触:
  源极: \(\phi = 0 + V_{bi}\) 或 \(\phi = 0\) (若以本征费米能级为参考)
  漏极: \(\phi = V_d + V_{bi}\) 或 \(\phi = V_d\)

• 肖特基接触:
  \(\phi = V_{applied} + \phi_B\)
  \(\phi_B\)是金属-半导体势垒。载流子浓度不再钉扎，电流由边界处载流子偏离平衡的浓度差驱动。

绝缘/对称边界: \(\partial\phi/\partial n = 0\) (诺伊曼条件), 电场法向分量为零，无载流子净流动穿过该边界。用于简化模型。
半导体-绝缘体界面: \(\epsilon_s(\partial\phi/\partial n)_s = \epsilon_{ox}(\partial\phi/\partial n)_{ox}\), 通常仅对多数载流子设表面复合边界.

NEGF

量子输运边界: 通过“自能 \(\Sigma\)”实现
自能 \(\Sigma\) 描述载流子如何从器件流入/流出半无限的源漏电极. 在哈密顿量 \(H\) 所描述的器件区域边界上， \(\Sigma_L\) 和 \(\Sigma_R\) 以算符形式“镶嵌”进去，数学上代表了开放系统的边界。它们决定了电极的能带结构、态密度以及最重要的载流子注入的费米分布 \(f(E-\mu_L)\), \(f(E-\mu_R)\)。

\(\mu_L\) 和 \(\mu_R\) 是左右电极的化学势（费米能级）。对于PMOS器件，在零偏压下，\(\mu_L = \mu_R = E_F\) 系统费米能级. 加偏压时，例如 \(V_ds > 0\), 则 \(\mu_L = E_F - q*V_{ds}\), \(\mu_R = E_F\) (假设漏端接地).

经典静电边界: 通过“电势（\(\phi\)）”实现
这是Poisson方程需要的. 采用第一类（狄利克雷）边界条件，即在源漏接触的边界网格点上，直接指定静电势 \(\phi_{contact}\) 的值. 对于理想欧姆接触，一个核心物理假设是：金属电极的费米能级\(\mu_{L/R}\)与半导体接触区的费米能级对齐。由此可推导出：

\[\phi_{contact} = \frac{\mu_{L/R}}{q} + \chi + \frac{E_g}{2q} + \phi_F \]

这里, 对于P+型源漏 (高掺杂), \(\chi\) 是半导体电子亲和能 (从真空能级到导带底), \(E_g / 2q\) 表示本征半导体的费米能级大约在带隙中央. \(\phi_F\) 是费米势, 由于高掺杂导致的费米能级相对本征位置的偏移。对于P型掺杂（\(N_A\)），\(\phi_F = -(k_B T / q) * ln(N_A / n_i)\)，

注意, 对于极高掺杂（\(N_A > 10^{19}\) cm^{-3}），需使用费米-狄拉克统计精确计算 \(\phi_F\)，并考虑带隙变窄效应，这会修正 \(E_g\) 和 \(\chi\)。若更为严格的话，在接触区，通过额外增加一个局域的自洽循环，使得该区域的静电势和电荷密度不仅满足泊松方程，也满足平衡状态的NEGF关系 (此时 \(G^<\)仅由单个费米函数决定). 通过这个循环自洽地确定接触区的有效费米能级，再将其作为电极的 \(\mu\)。

待整理

这就像测量山高，是选海平面为零点，还是选山脚下的村庄为零点，山顶的相对高度差（电势差）才是决定水流（电流）的关键。

两种参考系的详细解释

1. 参考系A：本征费米能级参考

• 核心：将未加偏压、未掺杂的本征半导体的费米能级 E_i 处的电势定义为 0。
• 边界条件：
  • 在n型重掺杂的源极接触，金属费米能级 E_fs 比 E_i 高 qV_bi。因此，静电势 φ_s = V_bi。
  • 漏极同理：φ_d = V_d + V_bi。
• 优点：非常符合器件物理的直觉，能清晰地看到金属-半导体接触处的能带弯曲和内建势垒。
• 缺点：在数值计算中，边界值随掺杂变化（V_bi 是掺杂浓度的函数），且整个计算域有一个恒定的电势背景值。

2. 参考系B：源极费米能级参考

• 核心：将源极金属的费米能级 E_fs 处的电势定义为 0。这是绝大多数现代商业和学术仿真工具的默认做法。
• 边界条件：
  • 源极接触：φ_s = 0。
  • 漏极接触：φ_d = V_d。
• 优点：
  1. 计算简洁：边界条件直接由外加电压 V_d 和 V_g 决定，与材料的具体掺杂参数解耦。
  2. 物理清晰：计算得到的电势 φ(x) 直接反映了相对于源极的电压降。
  3. 数值稳定：避免了在边界上处理一个可能很大的 V_bi 值，有利于迭代收敛。
• 内建电势去哪了？
  内建电势并没有消失。在参考系B中，V_bi 是通过在计算域内设置正确的掺杂分布来体现的。在源/漏重掺杂n型区域，泊松方程的电荷源项 ρ 中包含正的电离施主电荷。求解泊松方程时，系统为了满足边界条件（φ_s=0, φ_d=V_d）和内部的电荷分布，会自动地在从源端到沟道轻掺杂区的过渡区域产生一个电势变化，这个变化正是内建电势差。换句话说，V_bi 是计算的结果，而不是输入的边界条件。

关键结论与操作建议

1. 二者等价：只要仿真设置是自洽的，两种参考系最终计算出的电场分布（-∇φ）、能带图、载流子浓度和电流都完全一致。它们只相差一个全局的常数（V_bi）。
2. 现代仿真通行做法：在绝大多数科研论文和仿真软件（如Sentaurus, COMSOL, NanoTCAD ViDES等）的示例中，都采用参考系B（φ_s=0）。因为它将复杂的物理细节（功函数差、掺杂）封装在材料参数设置中，让使用者更专注于偏压的设置。
3. 你的操作指南：
   • 如果你在阅读文献或代码：看到 φ_s=V_bi，要意识到作者采用了绝对参考系。你需要检查他对费米能级的定义是否一致。
   • 如果你自己在设置仿真：强烈建议采用参考系B（φ_s=0, φ_d=V_d）。这更通用，更不易出错。
   • 确保自洽：无论采用哪种，都必须确保整个仿真模型的参数系统一致。例如，在参考系B中，如果你将 φ_s 设为0，那么：
     • 左引线（源）电子库的费米能级 E_Fs 也应设为 0。
     • 右引线（漏）电子库的费米能级 E_Fd 应设为 -qV_d（因为 E_F - qφ = constant）。
     • 材料的电子亲和能、功函数等参数也应按此参考系进行设置。