模型移植实战：BSIM模型向PSP模型的SPICE网表转换 - 第3节
3.1 PSP模型概述：从表面电势到紧凑模型
PSP模型（Philips Surface-Potential-based compact model）是业界公认的、基于表面电势的MOSFET紧凑模型标准。与BSIM系列基于阈值电压的分段建模思路不同，PSP模型从器件物理的根本——沟道表面电势出发，实现了从亚阈值到强反型、从线性区到饱和区的统一描述。
核心思想：PSP模型以泊松方程和电流连续性方程为基础，通过求解表面电势（ψs）作为沟道位置的函数，直接计算漏极电流、电荷和电容。这种“从头算起”的物理方法，使得PSP模型在模拟/射频电路设计中具有极高的精度和连续性。
3.2 核心方程：表面电势的求解
PSP模型的核心是求解沟道中任意位置x处的表面电势ψs(x)。该电势由以下隐式方程决定：
VGB - VFB - ψs(x) = γ · √[ψs(x) + φt · exp((ψs(x) - 2·φF - Vch(x)) / φt)]其中：
- VGB：栅-体偏压
- VFB：平带电压
- γ：体效应系数（γ = √(2·q·εsi·Nsub) / Cox）
- φt：热电压（kT/q）
- φF：费米势
- Vch(x)：沟道准费米势（从源端VS到漏端VD变化）
PSP vs BSIM 求解差异：BSIM使用显式的阈值电压公式分段描述电流，而PSP需要迭代求解上述隐式方程。虽然计算量略大，但PSP在亚阈值区、过渡区和强反型区之间实现了完全平滑的过渡，这对模拟电路（如运算放大器、比较器）的跨导和输出电阻建模至关重要。
3.3 参数体系总览：四大参数组
PSP模型的参数体系庞大而精细，按功能可分为四大参数组。下表列出了各组的关键参数及其物理含义：
| 参数组 | 典型参数 | 物理含义 | 典型值范围 |
|---|---|---|---|
| 几何/工艺参数 | TOXE, NSUB, GAMMA, VFB | 栅氧厚度、衬底掺杂浓度、体效应系数、平带电压 | TOXE: 1~5 nm; NSUB: 1e17~1e18 cm-3 |
| 物理/输运参数 | U0, VSAT, THESAT, MUE | 低场迁移率、饱和速度、速度饱和温度系数、迁移率退化 | U0: 200~600 cm²/Vs; VSAT: 1e7 cm/s |
| 二阶效应参数 | RDSW, CIT, DIBL, CLM | 源漏电阻、界面态电容、漏致势垒降低、沟道长度调制 | RDSW: 50~500 Ω·μm; CIT: 1e-5~1e-4 F/m² |
| 温度/噪声参数 | TNOM, TC1, TC2, NOIA, NOIB | 标称温度、迁移率温度系数、饱和速度温度系数、1/f噪声参数 | TNOM: 25°C; NOIA: 1e-30~1e-28 |
3.4 关键物理效应建模
PSP模型对以下关键物理效应提供了精确的物理描述，这些效应在BSIM模型中通常通过经验公式或分段函数处理：
3.4.1 迁移率退化效应
PSP使用统一的迁移率模型，同时考虑垂直电场和横向电场的影响：
μeff = μ0 / [1 + (MUE · Eeff)MUEEXP + (CS · Elat)CS]其中Eeff为有效垂直电场，Elat为横向电场。该公式在低场和高场下均保持物理连续性。
3.4.2 沟道长度调制（CLM）与DIBL
PSP通过表面电势的二维求解自然包含了DIBL效应，无需像BSIM那样引入额外的Vth修正项。CLM效应通过沟道夹断点的位置计算实现：
Leff = L - ΔLCLM - ΔLSCE其中ΔLCLM和ΔLSCE分别由横向电场和短沟道效应决定。
3.4.3 量子效应与多晶硅耗尽
PSP模型内置了量子力学修正（QM）和多晶硅栅耗尽效应（PDE），这些在先进工艺节点（≤45nm）中变得不可忽略：
- QM效应：通过增加等效氧化层厚度（ΔTox）来模拟载流子量子化导致的能带分裂
- PDE效应：通过栅电压降（ΔVpoly）修正实际栅压
移植注意事项：BSIM模型中的许多经验参数（如DIBL系数ETA0、CLM系数PCLM等）在PSP中没有直接对应项。移植时需要将这些效应“反解”为PSP的物理参数，而非简单的一一映射。例如，BSIM的DIBL效应需要转化为PSP中的沟道掺杂分布和栅氧厚度参数调整。
3.5 PSP模型与BSIM模型的核心差异对比
下表从建模方法论、参数体系、精度和计算效率四个维度对比PSP与BSIM模型：
| 对比维度 | PSP模型 | BSIM模型（BSIM4/6） |
|---|---|---|
| 建模基础 | 表面电势（ψs）统一求解 | 阈值电压（Vth）分段函数 |
| 区域连续性 | 全域连续，无分段点 | 需平滑函数拼接，存在残余不连续 |
| 参数数量 | 约300+（核心50~80） | 约400+（核心60~100） |
| 物理可解释性 | 高，参数直接对应物理量 | 中，部分参数为经验拟合系数 |
| 计算效率 | 中等（需迭代求解） | 高（显式公式） |
| 模拟/射频精度 | 极高（跨导、输出电阻连续） | 高（但高阶导数可能不连续） |
| 噪声模型 | 物理噪声模型（1/f, 热噪声） | 经验噪声模型 |
| 温度缩放 | 物理温度模型（Arrhenius型） | 经验温度系数 |
3.6 PSP模型在SPICE网表中的实例化
在SPICE仿真中，PSP模型通过.MODEL语句实例化，并指定模型类型为PSP或PSP103（取决于版本）。以下是一个典型的PSP模型定义示例：
****************************************************************
* PSP 103 Level 3 模型定义示例（28nm 工艺节点）
****************************************************************
.MODEL NMOS_28NM PSP103 LEVEL=3
+ TOXE = 1.2E-9 * 栅氧物理厚度 (m)
+ NSUB = 5.0E17 * 衬底掺杂浓度 (cm^-3)
+ GAMMA = 0.45 * 体效应系数 (V^0.5)
+ VFB = -0.95 * 平带电压 (V)
+ U0 = 350.0 * 低场迁移率 (cm²/Vs)
+ VSAT = 1.2E7 * 饱和速度 (cm/s)
+ THESAT = 1.5 * 速度饱和温度系数
+ MUE = 2.0 * 迁移率退化指数
+ RDSW = 200.0 * 源漏电阻 (Ω·μm)
+ CIT = 2.5E-5 * 界面态电容 (F/m²)
+ TNOM = 25.0 * 标称温度 (°C)
+ TC1 = -1.5E-3 * 迁移率温度系数 (1/°C)
+ NOIA = 1.0E-29 * 1/f 噪声参数 A
+ NOIB = 5.0E-25 * 1/f 噪声参数 B
****************************************************************
LEVEL参数说明：PSP103模型支持LEVEL=1（基础）、LEVEL=2（标准）和LEVEL=3（完整，包含所有二阶效应和噪声模型）。在模型移植过程中，建议从LEVEL=2开始验证，逐步开启LEVEL=3的高级特性。
3.7 本章小结
PSP模型作为基于表面电势的紧凑模型标准，其核心优势在于：
- 物理一致性：从泊松方程出发，统一描述所有工作区域
- 高精度连续性：跨导、输出电阻等关键指标在区域边界处完全平滑
- 参数可解释性：每个参数都有明确的物理含义，便于工艺关联和缩放
- 射频友好：物理噪声模型和NQS效应使其在RF设计中表现优异