Learning objectives
学完本课后，你应能：
- 用碰撞电离、Townsend 雪崩、击穿条件和鞘层解释直流气体放电的建立；
- 由 Debye 长度与 Bohm 速度估算壁面附近的空间尺度和粒子通量尺度；
- 区分电子温度、离子温度与中性气体温度，并用约化电场 \(E/N\) 与电子能量分布函数组织非平衡反应；
- 建立离子通量、离子能量、自由基通量与表面反应之间的加工模型；
- 对照明、Hall 推进器和环境等离子体完成能量与动量账本；
- 判断局域热平衡、碰撞性、准中性与连续介质近似何时失效。
Prerequisites
- Plasma State and Characteristic Scales：准中性、Debye 屏蔽、等离子体参数；
- 电势能 \(q\phi\)、碰撞截面、平均自由程；
- 指数增长、常微分方程和数量级估算；
- 功率、动量流率与效率。
Scope and conventions
- 采用 SI；电子温度常以电子伏特表示，并约定 \(k_{\mathrm B}T_e\) 的数值写成 eV；
- 一维平板鞘层坐标 \(x\) 从准中性体区指向壁面，壁面通常相对体区为负；
- 放电电压 \(V\) 指电极间电势差的大小；电流密度 \(J\) 的正方向按常规电流；
- 主要讨论弱电离、低温、碰撞性可调的实验室等离子体，不把所有低温等离子体等同于同一种放电；
- 表面模型默认稳态、均匀束流；真实反应器中的射频波形、三维输运与复杂化学需要数值模型。
Notation and modeling contract
| Symbol | Meaning | SI unit / convention |
|---|---|---|
| \(n_e,n_i,n_n\) | 电子、离子与中性粒子数密度 | \(\mathrm{m^{-3}}\) |
| \(\alpha\) | Townsend 第一电离系数 | \(\mathrm{m^{-1}}\) |
| \(\gamma\) | 二次电子发射系数 | 无量纲 |
| \(p,d\) | 气压、电极间距 | Pa, m |
| \(\lambda_D\) | 电子 Debye 长度 | m |
| \(u_B\) | Bohm 速度 | \(\mathrm{m\,s^{-1}}\) |
| \(E/N\) | 约化电场 | \(\mathrm{V\,m^2}\)；\(1\ \mathrm{Td}=10^{-21}\ \mathrm{V\,m^2}\) |
| \(f_e(\varepsilon)\) | 电子能量分布函数 EEDF | 按所用归一化 |
| \(k_r=\langle\sigma_r v\rangle\) | 反应速率系数 | \(\mathrm{m^3\,s^{-1}}\) |
| \(\Gamma_i\) | 离子粒子通量 | \(\mathrm{m^{-2}\,s^{-1}}\) |
1. Gas discharges and sheaths / 气体放电与鞘层
Knowledge ID: P13-U04-L03-K01
1.1 从种子电子到自持放电
外加电场使少量种子电子增能；当电子与中性粒子的碰撞能超过电离阈值，就产生新电子和正离子。在均匀场、忽略附着与空间电荷时，
\[\frac{\mathrm dn_e}{\mathrm dx}=\alpha n_e \quad\Longrightarrow\quad n_e(x)=n_e(0)e^{\alpha x}.\]
若正离子或光子轰击阴极引起二次电子发射，自持 Townsend 放电的理想条件是
\[\gamma\left(e^{\alpha d}-1\right)=1.\]
在经验关系 \(\alpha/p=A\exp(-Bp/E)\) 下，可得到 Paschen 形式
\[V_b= \frac{Bpd} {\ln(Apd)-\ln\!\left[\ln\!\left(1+\gamma^{-1}\right)\right]}.\]
它说明击穿电压主要受 \(pd\) 控制，但只在气体成分、温度、电极状态与近似均匀场固定时成立。
图 1 把指数倍增与阴极反馈画在同一账本里：没有反馈时雪崩会随种子电子消失，满足闭环增益条件后才可能自持。
图 2 强调放电的伏安特性可出现负微分电阻；电源与限流阻抗也是工作点模型的一部分。
1.2 为什么壁面前会形成鞘层
电子热速度远大于离子热速度，未带偏压的壁面起初更快收集电子，因而变负并排斥后续电子。准中性体区与壁面之间形成非准中性 鞘层（sheath）。尺度先由
\[\lambda_D= \sqrt{\frac{\varepsilon_0 k_{\mathrm B}T_e}{n_e e^2}}\]
控制；鞘层边缘的冷离子必须至少以 Bohm 速度进入鞘层：
\[u_i\ge u_B=\sqrt{\frac{k_{\mathrm B}T_e}{m_i}}.\]
该 Bohm 条件是理想平面、无碰撞、Maxwellian 电子模型的结果；碰撞、磁场、负离子或非 Maxwellian EEDF 会修正它。
图 3 中电势朝壁面降低；正离子因此加速，电子被势垒筛除。
图 4 是粒子流而非力图：红箭头表示离子收集，细蓝箭头分别表示能越过势垒的电子流向壁面和被势垒反射的电子。忽略 secondary-electron emission 与位移电流时，浮置壁的稳态净电流为零，即 \(\lvert J_e\rvert\simeq\lvert J_i\rvert\)；负壁势降低电子收集通量，却不会把它降为零。
2. Non-equilibrium electron and ion temperatures / 非平衡电子与离子温度
Knowledge ID: P13-U04-L03-K02
2.1 “低温”不等于所有组分同温
电子以速度撞击初始近静止、质量为 \(M\) 的重粒子时，一维正碰可转移的最大动能分数为
\[\left(\frac{\Delta\varepsilon_e}{\varepsilon_e}\right)_{\max} =\frac{4m_eM}{(m_e+M)^2} \simeq\frac{4m_e}{M}.\]
若再对各向散射角作适当平均，常用的代表性尺度约为 \(2m_e/M\)；该系数依微分截面和角权重，不能称为普适最大值。无论取最大或角平均，\(m_e/M\ll1\) 才是决定性能量交换缓慢的小参数。
因此电场可以优先加热电子，而离子和中性气体仍接近环境温度，常见排序为
\[T_e\gg T_i\approx T_g.\]
这就是 非平衡等离子体（nonequilibrium plasma） 的核心，不意味着系统违反热力学，而是不同自由度之间的能量交换时间比外部驱动与损失时间更长。
图 5 用三个能量库区分温度；它们不是同一条热平衡温标。
2.2 约化电场与 EEDF
电子在两次碰撞之间从场获得能量。平均自由程随中性密度 \(N\) 改变，所以控制电子动力学的常用相似参数是
\[\frac{E}{N}.\]
电子驱动反应 \(r\) 的速率系数不是“由平均温度直接代入”：
\[k_r=\langle\sigma_r v\rangle =\int_0^\infty \sigma_r(\varepsilon)v(\varepsilon) f_e(\varepsilon)\,\mathrm d\varepsilon.\]
当 EEDF 明显非 Maxwellian 时，单一 \(T_e\) 只能是约定的有效参数；阈值反应尤其由高能尾部决定。
图 6 说明提高 \(E/N\) 会重分配电子功率，但曲线不是所有气体通用的经验定律。
图 7 中两种 EEDF 即使平均能量相近，阈值右侧面积也可差很多。
3. Plasma processing and surface modification / 等离子体加工与表面改性
Knowledge ID: P13-U04-L03-K03
3.1 把“化学选择性”和“方向性”分开建模
体区电子碰撞产生自由基、激发态与离子。自由基通常以近各向同性热运动到达表面，负责选择性化学；正离子跨越鞘层获得近法向动量，帮助断键、去除钝化层或溅射。简化能量估算为
\[\mathcal E_i\simeq e\lvert V_s\rvert\]
但射频鞘层中离子能量分布取决于离子渡越时间、波形、碰撞与偏置，不能只用一个电压代表。
若 \(R\) 表示每单位面积、每单位时间去除的粒子数，量纲必须与入射粒子通量一致。一个不预设协同函数形式的概念模型是
\[R=\Gamma_iY_i+\Gamma_rY_r +R_{ir}(\Gamma_i,\Gamma_r,\theta,T_s),\]
其中 \(Y_i,Y_r\) 是无量纲的每入射粒子有效产额，而 \(R_{ir}\) 本身必须具有 \(\mathrm{m^{-2}\,s^{-1}}\)，表示离子—自由基协同去除通量。若希望把协同写成“每个离子额外去除多少”，可令 \(R_{ir}=\Gamma_iY_{ir}(\Gamma_r/\Gamma_{r0},\theta,T_s)\)，此时 \(Y_{ir}\) 无量纲，\(\Gamma_{r0}\) 是明示的参考通量。直接写 \(\Gamma_i\Gamma_rY_{ir}\) 且宣称 \(Y_{ir}\) 无量纲会得到错误的 \(\mathrm{m^{-4}\,s^{-2}}\)。
图 8 把体区生成、鞘层加速与表面反应分成三个子模型。
图 9 显示为什么法向离子通量可形成各向异性刻蚀；若侧壁也被强烈轰击，几何就会改变。
3.2 表面不是被动边界
壁面会吸附、复合、充电、发射二次电子并改变局部电场。绝缘结构可累积电荷；高深宽比沟槽会出现阴影效应与充电偏转。因而可靠工艺必须同时控制：
- 物种通量与能量分布；
- 基片温度与表面覆盖；
- 偏压、频率和脉冲占空比；
- 选择比、均匀性、损伤和颗粒污染。
4. Lighting, propulsion, and environmental applications / 照明、推进与环境应用
Knowledge ID: P13-U04-L03-K04
4.1 照明：把电子能量转成辐射
气体放电灯中，电子碰撞激发原子或分子，随后辐射退激。荧光灯还用紫外光激发荧光粉。能量链为
\[P_{\rm electric}\rightarrow P_{\rm electron}\rightarrow P_{\rm excitation}\rightarrow P_{\rm radiation}\rightarrow P_{\rm visible}+P_{\rm loss}.\]
光效受碰撞分支、辐射俘获、镇流器、荧光转换和壁面损失共同限制。
图 10 是能量账本，不把“电子温度高”误写成灯管整体很热。
4.2 电推进：动量来自高速喷流
Hall 推进器用磁场限制电子、用轴向电场加速离子。若离子束速度 \(v_i\) 且推进剂质量流率 \(\dot m\)，忽略压强项，
\[F\simeq \dot m v_i, \qquad P_{\rm beam}=\frac12\dot m v_i^2, \qquad I_{\rm sp}=\frac{v_i}{g_0}.\]
高比冲不等于高推力；在功率受限下，提高喷速会降低可承受的质量流率。
4.3 环境应用：反应性物种而非“万能消毒”
大气压非热等离子体可生成 O、OH、臭氧、激发态与带电粒子，用于表面活化、废气和水处理。效能必须以剂量、传质、选择性、副产物和能耗评估；“有等离子体”本身不是安全或处理完成的证据。
Worked Examples
Worked Example 1 — Paschen 击穿估算
已知： 某气体经验常数 \(A=15\ \mathrm{cm^{-1}\,Torr^{-1}}\)、\(B=365\ \mathrm{V\,(cm\,Torr)^{-1}}\)，\(\gamma=0.010\)，且 \(pd=1.00\ \mathrm{Torr\,cm}\)。
系统与模型： 均匀平板电极、Townsend—Paschen 经验模型；忽略电极粗糙与附着。
求： 击穿电压。
\[V_b=\frac{365} {\ln 15-\ln[\ln(1+100)]}\ \mathrm V.\]
数值为
\[\ln15=2.708,\qquad \ln(\ln101)=1.529,\]
\[V_b=3.10\times10^2\ \mathrm V.\]
检查： 分母无量纲，\(Bpd\) 给出伏特；结果为典型几百伏数量级。若 \(pd\) 远离适用范围或场强不均匀，本结果不能外推。
Worked Example 2 — 氩等离子体的 Debye 长度与 Bohm 速度
已知： \(n_e=1.0\times10^{16}\ \mathrm{m^{-3}}\)，\(k_{\mathrm B}T_e=3.0\ \mathrm{eV}\)，氩离子质量 \(m_i=6.63\times10^{-26}\ \mathrm{kg}\)。
参考系： 实验室系；冷离子、平面无碰撞鞘层。
求： \(\lambda_D\)、\(u_B\)。
\[\lambda_D= \sqrt{\frac{\varepsilon_0(3.0e)}{n_e e^2}} =1.29\times10^{-4}\ \mathrm m =0.129\ \mathrm{mm},\]
\[u_B=\sqrt{\frac{3.0e}{6.63\times10^{-26}}} =2.69\times10^3\ \mathrm{m\,s^{-1}}.\]
检查： \(\lambda_D\) 远小于厘米级反应器，支持体区准中性；\(u_B\) 是鞘边界速度尺度，不是整个体区所有离子的固定速度。
Worked Example 3 — 最大与角平均能量交换尺度
已知： 电子与氩原子弹性碰撞，\(m_e=9.11\times10^{-31}\ \mathrm{kg}\)，\(M_{\rm Ar}=6.63\times10^{-26}\ \mathrm{kg}\)。
模型： 重粒子近静止、非相对论二体弹性碰撞；正碰给最大值，\(2m_e/M\) 仅作常见角平均数量级。
求： 单次正碰最大能量交换分数，并与角平均尺度比较。
\[\frac{4m_eM_{\rm Ar}}{(m_e+M_{\rm Ar})^2} =5.50\times10^{-5}, \qquad \frac{2m_e}{M_{\rm Ar}}=2.75\times10^{-5}.\]
检查与解释： 两个量都远小于 1，且最大值约为角平均尺度的两倍；因此需很多次弹性碰撞才能把电子能量转给中性气体。非弹性激发与电离仍可高效消耗电子能量。数量级支持 \(T_e\gg T_g\)，但不代表气体永远不会因功率沉积而升温。
Worked Example 4 — Hall 推进器动量与功率账本
已知： \(\dot m=5.0\ \mathrm{mg\,s^{-1}}=5.0\times10^{-6}\ \mathrm{kg\,s^{-1}}\)，\(v_i=2.0\times10^4\ \mathrm{m\,s^{-1}}\)。
系统边界： 推进器与喷流；忽略中和电子的动量、羽流发散和压强项。正方向沿喷流反向的推力。
求： 推力、束流功率和理想比冲。
\[F=\dot m v_i=0.100\ \mathrm N,\]
\[P_{\rm beam}=\frac12\dot m v_i^2 =1.00\times10^3\ \mathrm W,\]
\[I_{\rm sp}=\frac{v_i}{g_0} =2.04\times10^3\ \mathrm s.\]
检查： \(P_{\rm beam}=Fv_i/2\) 一致。若输入电功率为 \(2.0\ \mathrm{kW}\)，该简化束流效率为 \(50\%\)；真实总效率还受电离、发散、壁损与中和限制。
Lesson synthesis: a transferable modeling workflow
- 辨认驱动： 直流、射频、微波还是脉冲；确定 \(E/N\) 与能量输入位置；
- 建立体区粒子账本： 电离、附着、复合、扩散与壁损；
- 建立电子能量账本： EEDF 决定 \(\langle\sigma v\rangle\)，不要只凭平均温度；
- 连接边界： 鞘层把体区物种转换成壁面的能量—角度分布；
- 连接应用指标： 工艺速率、光效、推力、污染物去除率；
- 做边界检查： 准中性、碰撞性、稳态、局域场与表面均匀性是否成立。
Formula-and-decision table
| Decision | Relation | Use when | Main limitation |
|---|---|---|---|
| 雪崩增长 | \(n_e=n_0e^{\alpha x}\) | 弱空间电荷、均匀场 | 不含附着、扩散与场畸变 |
| 自持条件 | \(\gamma(e^{\alpha d}-1)=1\) | Townsend 击穿 | \(\gamma\) 与表面状态有关 |
| 鞘层尺度 | \(\lambda_D=\sqrt{\varepsilon_0k_BT_e/(n_ee^2)}\) | 近 Maxwellian 电子 | 多组分需广义屏蔽长度 |
| 鞘边离子速度 | \(u_B=\sqrt{k_BT_e/m_i}\) | 冷离子、无碰撞平面鞘 | 磁化、碰撞、负离子会修正 |
| 反应速率 | \(k_r=\langle\sigma_rv\rangle\) | 已知 EEDF 与截面 | 截面与 EEDF 均有不确定度 |
| 表面离子能量 | \(\mathcal E_i\simeq e\lvert V_s\rvert\) | 稳态无碰撞鞘层估算 | 射频与碰撞使其成为分布 |
| 推力 | \(F\simeq\dot m v_i\) | 准直稳态喷流 | 忽略发散、压强与多价离子 |
Model boundaries
- 低温不是低电子能量： 常指重粒子温度低；电子可达数 eV；
- 准中性不覆盖鞘层： 体区 \(n_e\simeq n_i\)，壁面附近必须允许空间电荷；
- Paschen 定律非普适： 微间隙、强场发射、高压多体碰撞和复杂混合气体会偏离；
- 温度参数非总是热力学温度： 非 Maxwellian EEDF 下要注明定义；
- 表面反应非仅由入射能量决定： 覆盖率、温度、晶向与历史均重要；
- 环境应用需副产物核查： 臭氧、NOx 与活性物种可能有安全边界；
- 连续模型与动理学模型选择： 若平均自由程或 Debye 长度不远小于几何尺度，应使用动理学或混合模型。
Misconception clinic
Practice and retrieval
- [K01 概念] 为什么单次电子雪崩不足以保证自持放电？
- [K01 推导] 从 \(\mathrm dn/\mathrm dx=\alpha n\) 推导 \(n(d)/n(0)\)。
- [K01 定量] \(T_e=2.0\ \mathrm{eV}\)、\(n_e=4.0\times10^{16}\ \mathrm{m^{-3}}\)，求 \(\lambda_D\)。
- [K02 概念] 为什么 \(E/N\) 比单独的 \(E\) 更适合比较不同气压放电？
- [K02 图像] 若阈值右侧 EEDF 面积增加而 \(n_e\) 不变，阈值反应速率如何变化？
- [K02 定量] 氙原子质量 \(2.18\times10^{-25}\ \mathrm{kg}\)，分别估算正碰最大分数 \(4m_e/M\) 与角平均尺度 \(2m_e/M\)。
- [K03 建模] 分别指出自由基与离子对各向异性刻蚀的主要作用。
- [K03 定量] 刻蚀速率 \(180\ \mathrm{nm\,min^{-1}}\)，运行 \(12.0\ \mathrm{min}\)，去除厚度是多少？
- [K04 定量] 推进器 \(\dot m=2.0\ \mathrm{mg\,s^{-1}}\)、\(v_i=30\ \mathrm{km\,s^{-1}}\)，求 \(F\) 与 \(P_{\rm beam}\)。
- [综合边界] 反应器尺寸 \(L=0.50\ \mathrm{mm}\)、\(\lambda_D=0.20\ \mathrm{mm}\)，能否把整个反应器视为“薄鞘层加准中性体区”？
1. Answer
自持还需要阴极二次发射、光电离或其他反馈不断补充新一代种子电子；否则雪崩粒子被电极收集后终止。
2. Answer
分离变量得 \(\mathrm dn/n=\alpha\,\mathrm dx\)，积分为 \(\ln[n(d)/n(0)]=\alpha d\)，故 \(n(d)/n(0)=e^{\alpha d}\)。
3. Answer
\[\lambda_D=\sqrt{\frac{\varepsilon_0(2.0e)}{(4.0\times10^{16})e^2}} =5.26\times10^{-5}\ \mathrm m =52.6\ \mu\mathrm m.\]
4. Answer
平均自由程近似随 \(1/N\) 变化，电子在碰撞间获得的能量尺度与 \(E\lambda\) 有关，因此由 \(E/N\) 控制更能保持碰撞动力学相似。
5. Answer
\(k_r=\int\sigma v f_e\,\mathrm d\varepsilon\) 的阈值以上贡献增大，所以在其他条件不变时反应速率 \(n_en_tk_r\) 增大。
6. Answer
\[\frac{4m_e}{M}=1.67\times10^{-5}, \qquad \frac{2m_e}{M}=8.36\times10^{-6}.\]
7. Answer
自由基提供材料选择性化学与挥发性产物；鞘层定向离子提供法向动量、断键和去钝化，因此控制侧壁与底部差异。
8. Answer
\(d=Rt=(180)(12.0)=2160\ \mathrm{nm}=2.16\ \mu\mathrm m\)。有效数字为三位。
9. Answer
\[F=(2.0\times10^{-6})(3.0\times10^4) =6.0\times10^{-2}\ \mathrm N,\] \[P_{\rm beam}=\frac12(2.0\times10^{-6})(3.0\times10^4)^2 =9.0\times10^2\ \mathrm W.\]
10. Answer
不宜。\(\lambda_D/L=0.40\) 并非小参数，鞘层可能占据系统显著部分，准中性体区与薄边界层分离失效，应考虑全域空间电荷或动理学模型。
Exit ticket
- 设计一个区分“高电子温度”和“高气体温度”的双诊断方案，并说明各自测量对象。
- 给定一种新应用，按“体区反应—鞘层—表面或喷流—性能指标—副作用”写出最小模型。
References and source notes
以下资料均于 2026-08-31 实际打开核验；正文为重新组织的教学表达，图均依方程原创：
- MIT OpenCourseWare, Sheaths — 鞘层、Bohm 条件与空间推进课程背景。
- Princeton Plasma Physics Laboratory SULI, Physics of Low Temperature Plasmas — 非平衡、放电与 EEDF 教学框架。
- NIST, Electricity Division: Programs, Activities, and Accomplishments (NISTIR 6587) — “Plasma Processing and Gaseous Dielectrics”（pp. 68–69）的等离子体加工数据、反应与模型需求。
- NASA Technology Transfer, Annular Hall Thruster — Hall 推进器结构与应用。
- U.S. Department of Energy, General Plasma Science — 低温等离子体研究范围与应用。
- NASA Small Spacecraft Systems Virtual Institute, In-Space Propulsion — 电推进性能指标与工程边界。
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| MIT OCW Sheaths | 鞘层与 Bohm 进入条件 | K01 §1.2 | 冷离子、平面鞘层 |
| PPPL SULI LTP | 放电、非平衡与 EEDF | K01–K02 | 低温弱电离等离子体 |
| NIST IR 6587 | 工艺反应、表面与数据需求 | K03 | 工业加工需耦合化学与输运 |
| NASA Hall Thruster | Hall 推进器物理图景 | K04 §4.2 | 工程装置并非理想一维束流 |
| DOE General Plasma Science | 应用范围 | K03–K04 | 分类性来源，不替代具体装置数据 |
| NASA In-Space Propulsion | 推力、功率与系统限制 | K04 / Example 4 | 需计入总效率和子系统 |
复习与拓展
- 后续可加入 Boltzmann 方程求 EEDF、global model 与 particle-in-cell/Monte Carlo collision 的层级比较；
- 与材料科学接口：溅射产额、表面吸附动力学、原子层刻蚀；
- 与 Unit 14.1 接口：射频匹配网络、光谱诊断、数据采集与闭环工艺控制。