等离子体状态与特征尺度

Learning objectives / 学习目标

完成本课后，你应当能够：

1. 区分电离与等离子体集体行为，说明“完全电离”不是等离子体的必要条件。
2. 写出多组分等离子体的准中性条件，并解释“准中性”为什么不等于“处处电荷密度严格为零”。
3. 从 Boltzmann 电子响应与 Poisson 方程推导电子 Debye 长度及 Yukawa 屏蔽势。
4. 用电子薄片位移模型推导电子等离子体频率，并区分角频率、频率和周期。
5. 计算 Debye 球内粒子数 \(N_D\) 与 Coulomb 耦合参数 \(\Gamma\)，判断经典弱耦合近似是否合理。
6. 用 \(n\)、\(T\) 和系统尺度 \(L\) 对一个给定系统作数量级体检，并说清结论的适用边界。

Scope and conventions / 范围与约定

本课研究经典、非相对论、近 Maxwell 分布的多组分等离子体，并把系统尺度 \(L\)、可用时间 \(\tau_{\rm avail}\) 与 Debye/response scales 明确比较。默认 SI；体区准中性不外推到 sheath、double layer 或 \(L\sim\lambda_D\) 的结构。强耦合、量子简并、强磁化和高碰撞性只作为模型边界，不用弱耦合公式强行覆盖。

Prerequisites / 先修知识

• Coulomb force、Gauss/Poisson equation 与 electric potential；
• simple harmonic motion、Boltzmann factor 与 Taylor expansion；
• 数密度、温度、量纲与数量级估计；
• 明确 system size、boundary 与 response time 的建模习惯。

Notation and modeling contract / 符号与建模契约

0.1 Temperature has two common meanings

等离子体教材经常把温度直接写成能量。为了避免少一个 \(k_B\)，本笔记统一采用：

\[1\ \mathrm{eV}=1.602\,176\,634\times10^{-19}\ \mathrm J, \qquad 1\ \mathrm{eV}/k_B\approx1.1605\times10^4\ \mathrm K.\]

0.2 A working definition

判断一个给定过程是否处在常规等离子体区间，可以先问：

\[\boxed{L\gg\lambda_D} \qquad \boxed{\omega_{pe}\tau_{\rm avail}\gg1} \qquad \boxed{N_D=\frac{4\pi}{3}n\lambda_D^3\gg1}\]

Knowledge 1 — Ionization and collective behavior / 电离与集体行为

Knowledge ID: P13-U01-L01-K01

1.1 Ionization is a composition statement

对最简单的单电离原子气体，定义电离度（ionization fraction）

\[\boxed{\alpha\equiv\frac{n_i}{n_n+n_i}}\]

其中 \(n_n\) 是中性原子密度，\(n_i\) 是正离子密度：

• \(\alpha\ll1\)：弱电离（weakly ionized）；
• \(\alpha\sim1\)：高度或近完全电离；
• 若存在多价离子，应分别记录各电荷态，而不是把所有离子简单合并。

电离可由高温碰撞、强电场、光电离或高能粒子触发；复合过程则把自由电子重新束缚。一个系统的稳态电离度取决于产生与损失的竞争，并不只由“温度够不够高”决定。低温放电、D 层电离层和尘埃等离子体都说明：**等离子体不必完全电离，也不必处于热平衡。**¹³

1.2 Collective behavior is a dynamical statement

中性气体中，分子主要在短程碰撞时显著改变彼此运动；等离子体中，一小部分电荷位移就能产生跨越许多粒子间距的电场。每个粒子既是场的源，又受总场的作用：

\[\{f_s\} \longrightarrow (\rho,\mathbf J) \longrightarrow (\mathbf E,\mathbf B) \longrightarrow \{f_s\}.\]

这里 \(f_s(\mathbf r,\mathbf v,t)\) 表示物种 \(s\) 的分布函数。这个闭环称为自洽（self-consistent）响应，也是等离子体波、屏蔽、不稳定性和输运的共同起点。

1.3 Ionized gas versus plasma: a scale-dependent distinction

设自由电子刚刚出现：若电子在建立集体电场前便撞到中性粒子、逃出系统或被复合，那么“有电离”仍不足以保证典型等离子体行为。反之，即使电离度很低，只要集体响应比这些破坏过程更快、更长程，系统仍可按等离子体处理。

一个常用的附加比较是

\[\omega_{pe}\gtrsim\nu_{\mathrm{loss}},\]

其中 \(\nu_{\mathrm{loss}}\) 可代表有效碰撞、附着或逃逸速率。它不是普适定义，而是提醒我们：必须比较竞争时间尺度。

Knowledge 2 — Quasineutrality / 准中性

Knowledge ID: P13-U01-L01-K02

2.1 Core statement

对任意多组分等离子体，局域电荷密度为

\[\rho(\mathbf r,t)=\sum_s q_sn_s(\mathbf r,t).\]

在远大于 Debye 长度的等离子体体区（bulk），通常有

\[\boxed{\sum_s q_sn_s\approx0}.\]

对电子与单一正离子物种（电荷 \(+Ze\)）：

\[\boxed{n_e\approx Zn_i}.\]

2.2 Why a tiny imbalance is enough

Poisson 方程给出

\[\nabla^2\phi=-\frac{\rho}{\varepsilon_0}.\]

若电势在尺度 \(L\) 上改变约 \(\phi\)，则数量级上

\[|\rho|\sim\varepsilon_0\frac{|\phi|}{L^2}.\]

再取常见的热电势尺度 \(e|\phi|\sim k_BT_e\)：

\[\frac{|\rho|}{en_e} \sim \frac{\varepsilon_0k_BT_e}{n_ee^2L^2} = \boxed{\left(\frac{\lambda_{De}}{L}\right)^2}.\]

因此，只要 \(L\gg\lambda_{De}\)，维持有限电势所需的相对电荷分离就极小。准中性不是额外假定出来的魔法，而是长尺度极限的结果。

Worked Example 1 — How neutral is “quasi”?

2.3 Bulk and boundary layer

在等离子体体区，\(n_e\) 与 \(Zn_i\) 几乎重合；靠近吸收电子和离子的材料壁面，二者响应速度不同，形成厚度通常为若干 \(\lambda_D\) 的鞘层（sheath）。鞘层、电双层、局域带电束以及高频电荷密度波都是准中性近似需要谨慎或直接失效的地方。⁶⁴

2.4 Scale hierarchy

常规弱耦合、宏观体区常出现

\[\boxed{a\ll\lambda_D\ll L},\]

其中 \(a\) 是平均粒子间距的量级。左边保证屏蔽云包含许多粒子；右边保证屏蔽层相对系统很薄。

Knowledge 3 — Debye shielding and Debye length / 德拜屏蔽与德拜长度

Knowledge ID: P13-U01-L01-K03

3.1 Physical picture: the plasma polarizes

在均匀等离子体中放入一个正试探电荷 \(Q\)：

1. 电子被吸引，正离子被排斥；
2. 试探电荷周围形成净负的屏蔽云（shielding cloud）；
3. 距离很近时仍看见近似 Coulomb 势；
4. 超过某个特征距离后，试探电荷与屏蔽云的总电荷几乎抵消，势快速衰减。

这个特征距离就是 Debye 长度。它不是硬边界，也不是“粒子之间的距离”；它是指数屏蔽的 e-folding length。¹⁷

3.2 Derivation — fixed ions and Boltzmann electrons

电子电荷为 \(-e\)，其静电势能是 \(U_e=-e\phi\)。Boltzmann 因子给出

\[n_e(\mathbf r) =n_0\exp\!\left(-\frac{U_e}{k_BT_e}\right) =n_0\exp\!\left(\frac{e\phi}{k_BT_e}\right).\]

令无量纲电势

\[u\equiv\frac{e\phi}{k_BT_e}.\]

当 \(|u|\ll1\) 时，

\[e^u=1+u+\frac{u^2}{2!}+\cdots\approx1+u,\]

所以

\[n_e\approx n_0\left(1+\frac{e\phi}{k_BT_e}\right).\]

Poisson 方程必须包含外加扰动电荷密度 \(\rho_{\rm ext}\)：

\[\nabla^2\phi =-\frac{\rho_{\rm ext}+en_i-en_e}{\varepsilon_0}.\]

代入 \(n_i=n_0\) 和线性化电子密度：

\[\nabla^2\phi =-\frac{\rho_{\rm ext}}{\varepsilon_0} +\frac{n_0e^2}{\varepsilon_0k_BT_e}\phi.\]

定义电子 Debye 长度

\[\boxed{ \lambda_{De} \equiv \sqrt{\frac{\varepsilon_0k_BT_e}{n_ee^2}} }\]

（在线性化精度内可令 \(n_e=n_0\)），得到线性化 Poisson–Boltzmann 方程：

\[\boxed{ \left(\nabla^2-\frac{1}{\lambda_{De}^2}\right)\phi =-\frac{\rho_{\rm ext}}{\varepsilon_0} }.\]

3.3 Point charge: the Yukawa or screened Coulomb potential

对位于原点的点电荷 \(Q\)，

\[\rho_{\rm ext}=Q\,\delta^{(3)}(\mathbf r),\]

球对称解为

\[\boxed{ \phi(r)=\frac{Q}{4\pi\varepsilon_0r} \exp\!\left(-\frac{r}{\lambda_{De}}\right) }.\]

它由两部分组成：

\[\underbrace{\frac{Q}{4\pi\varepsilon_0r}}_{\text{vacuum Coulomb potential}} \times \underbrace{e^{-r/\lambda_{De}}}_{\text{screening factor}}.\]

3.4 Electron Debye length versus total Debye length

若多个物种都能在所研究的时间尺度上作 Boltzmann 响应，线性化后各物种的屏蔽贡献相加：

\[\boxed{ \frac{1}{\lambda_D^2} =\sum_s\frac{n_{s0}q_s^2}{\varepsilon_0k_BT_s} }.\]

对单电离、\(n_e=n_i=n_0\) 的电子—离子等离子体：

\[\frac{1}{\lambda_D^2} =\frac{n_0e^2}{\varepsilon_0k_B} \left(\frac1{T_e}+\frac1{T_i}\right).\]

若 \(T_e=T_i\)，则 \(\lambda_D=\lambda_{De}/\sqrt2\)。这解释了为什么不同教材的 Yukawa 指数有时相差 \(\sqrt2\)：有的把 \(\lambda_D\) 定义为单物种电子长度，却同时让电子和离子都响应；有的直接把所有物种并入有效总长度。本笔记的原则是：

3.5 Scaling laws and practical formulas

电子 Debye 长度的核心标度是

\[\boxed{\lambda_{De}\propto\sqrt{\frac{T_e}{n_e}}}.\]

• 温度升高：电子更难被小电势束缚，屏蔽云更大；
• 密度升高：可供重排的电荷更多，屏蔽距离更短。

将 \(T_e\) 用 eV、\(n_e\) 用 \(\mathrm{m^{-3}}\) 表示：

\[\boxed{ \lambda_{De}[\mathrm m] =7.43\times10^3 \sqrt{\frac{T_e[\mathrm{eV}]}{n_e[\mathrm{m^{-3}}]}} }.\]

等价的常见 cgs 实用形式是

\[\lambda_{De}[\mathrm{cm}] =7.43\times10^2 \sqrt{\frac{T_e[\mathrm{eV}]}{n_e[\mathrm{cm^{-3}}]}}.\]

Worked Example 2 — Debye length of a laboratory plasma

给定

\[n_e=1.0\times10^{16}\ \mathrm{m^{-3}}, \qquad T_e=2.0\ \mathrm{eV}.\]

用 SI 基础公式时，\(k_BT_e=(2.0\ \mathrm{eV})e=3.204\times10^{-19}\ \mathrm J\)：

\[\lambda_{De} =\sqrt{ \frac{(8.854\times10^{-12})(3.204\times10^{-19})} {(1.0\times10^{16})(1.602\times10^{-19})^2} }.\]

\[\boxed{ \lambda_{De}=1.05\times10^{-4}\ \mathrm m =0.105\ \mathrm{mm} }.\]

若容器尺度 \(L=1.0\ \mathrm{cm}\)，则

\[\frac{L}{\lambda_{De}}\approx95.\]

因此容器体区可以很好地准中性；但距材料边界几倍 \(\lambda_{De}\) 的区域仍必须保留 Poisson 方程。

3.6 What the derivation does not cover

Knowledge 4 — Plasma frequency / 等离子体频率

Knowledge ID: P13-U01-L01-K04

4.1 Physical picture: charge separation creates a restoring force

考虑均匀电子流体与静止正离子背景。令电子整体向 \(+x\) 方向微移 \(\xi(t)>0\)：

• 左端露出正离子，形成 \(+\sigma\)；
• 右端电子过剩，形成 \(-\sigma\)；
• 两面之间的电场从左向右；
• 电子受力 \(-e\mathbf E\) 向左，指向平衡位置。

4.2 Derivation — electron slab model

电子位移形成的面电荷密度大小为

\[\sigma=en_e\xi.\]

两张相反无限面电荷之间的合电场为

\[E=\frac{\sigma}{\varepsilon_0} =\frac{en_e}{\varepsilon_0}\xi.\]

电子运动方程为

\[m_e\ddot\xi=-eE =-\frac{n_ee^2}{\varepsilon_0}\xi.\]

于是

\[\boxed{\ddot\xi+\omega_{pe}^2\xi=0},\]

其中

\[\boxed{ \omega_{pe} =\sqrt{\frac{n_ee^2}{\varepsilon_0m_e}} }\]

称为电子等离子体角频率（electron plasma angular frequency）。这是简谐振子：

\[\xi(t)=\xi_0\cos(\omega_{pe}t+\delta).\]

该推导显示了集体行为的核心：恢复力不是某个指定离子对某个电子的“弹簧”，而是微小电荷分离产生的整体自洽电场。⁴⁶

4.3 Do not mix \(\omega\), \(f\), and the period

\[\boxed{f_{pe}=\frac{\omega_{pe}}{2\pi}}, \qquad \boxed{\tau_{pe}=\frac1{f_{pe}}=\frac{2\pi}{\omega_{pe}}}.\]

4.4 Phase and energy exchange

若取 \(\delta=0\)：

\[\frac{\xi}{\xi_0}=\cos\omega_{pe}t, \qquad \frac{v}{\omega_{pe}\xi_0}=-\sin\omega_{pe}t.\]

最大位移时，电荷分离和场能最大、速度为零；穿过平衡位置时，场能最小、电子集体动能最大。

4.5 Every mobile species has a plasma frequency

一般物种 \(s\) 的特征频率为

\[\boxed{ \omega_{ps}=\sqrt{\frac{n_sq_s^2}{\varepsilon_0m_s}} }.\]

对电荷 \(+Ze\) 的离子：

\[\omega_{pi}=\sqrt{\frac{n_iZ^2e^2}{\varepsilon_0m_i}}.\]

若 \(n_e=Zn_i\)，则

\[\frac{\omega_{pi}}{\omega_{pe}} =\sqrt{\frac{Zm_e}{m_i}}.\]

氢等离子体中 \(Z=1\)、\(m_i=m_p\)：

\[\frac{\omega_{pi}}{\omega_{pe}}\approx\frac1{42.85}.\]

因此没有特别说明时，“the plasma frequency”通常指电子等离子体频率。电子能在离子几乎未动时先完成快速响应。

4.6 The bridge between space and time scales

定义本笔记采用的电子热速度

\[\boxed{v_{te}\equiv\sqrt{\frac{k_BT_e}{m_e}}}.\]

把 \(\lambda_{De}\) 与 \(\omega_{pe}\) 相乘：

\[\omega_{pe}\lambda_{De} =\sqrt{\frac{n_ee^2}{\varepsilon_0m_e}} \sqrt{\frac{\varepsilon_0k_BT_e}{n_ee^2}} =\sqrt{\frac{k_BT_e}{m_e}} =v_{te}.\]

所以

\[\boxed{\lambda_{De}=\frac{v_{te}}{\omega_{pe}}}.\]

它把本课的空间尺度与时间尺度连接起来：电子以热速度在响应时间 \(\omega_{pe}^{-1}\) 内移动的距离，正是 \(\lambda_{De}\) 的量级。

Worked Example 3 — Plasma frequency and response time

把密度用 SI 的 \(\mathrm{m^{-3}}\) 表示：

\[\boxed{ \omega_{pe}[\mathrm{rad\,s^{-1}}] =56.4\sqrt{n_e[\mathrm{m^{-3}}]} }\]

\[\boxed{ f_{pe}[\mathrm{Hz}] =8.98\sqrt{n_e[\mathrm{m^{-3}}]} }.\]

对前例 \(n_e=1.0\times10^{16}\ \mathrm{m^{-3}}\)：

\[\omega_{pe}=5.64\times10^9\ \mathrm{rad\,s^{-1}},\]

\[f_{pe}=8.98\times10^8\ \mathrm{Hz}=0.898\ \mathrm{GHz},\]

\[\boxed{\tau_{pe}=1.11\times10^{-9}\ \mathrm s=1.11\ \mathrm{ns}}.\]

同时，\(T_e=2.0\ \mathrm{eV}\) 给出

\[v_{te}=5.93\times10^5\ \mathrm{m\,s^{-1}},\]

且

\[\frac{v_{te}}{\omega_{pe}} =1.05\times10^{-4}\ \mathrm m =\lambda_{De},\]

与 Debye 长度的独立计算完全一致。

4.7 Model boundary

Knowledge 5 — Plasma parameter and weak/strong coupling / 等离子体参数与弱强耦合

Knowledge ID: P13-U01-L01-K05

5.1 How many particles build one shielding cloud?

以 \(\lambda_D\) 为半径的 Debye 球体积为

\[V_D=\frac{4\pi}{3}\lambda_D^3.\]

其中平均粒子数称为 Debye number 或一种常见定义下的 plasma parameter：

\[\boxed{ N_D\equiv\frac{4\pi}{3}n\lambda_D^3 }.\]

当 \(N_D\gg1\) 时：

• 一个屏蔽云由许多粒子共同构成，连续密度近似合理；
• 粒子数的 Poisson 相对涨落量级 \(\delta N/N\sim N_D^{-1/2}\) 很小；
• “平均场 + 小涨落”的弱耦合动理学图像有希望成立。

当 \(N_D\) 不大时，单个邻近粒子、离散涨落和粒子相关性不能再被平滑平均掉。

5.2 Interparticle spacing and Coulomb coupling

定义 Wigner–Seitz 半径 \(a\)，使每个粒子平均占据一个半径为 \(a\) 的球形体积：

\[\frac{4\pi}{3}a^3=\frac1n, \qquad \boxed{a=\left(\frac{3}{4\pi n}\right)^{1/3}}.\]

相距 \(a\) 的同号电荷，其典型 Coulomb 势能与热能之比为

\[\boxed{ \Gamma \equiv \frac{q^2/(4\pi\varepsilon_0a)}{k_BT} }.\]

对电荷数为 \(Z\) 的粒子，应使用 \(q=Ze\)；对多组分等离子体，\(\Gamma\) 是物种与相互作用相关的，电子、离子和尘埃可能处在完全不同的耦合区。

5.3 Exact relation between \(N_D\) and \(\Gamma\)—under matched conventions

对单组分经典系统，若 \(\lambda_D\)、\(a\)、\(\Gamma\) 都使用同一个 \(n,T,q\)：

\[N_D =\frac{4\pi}{3}n\lambda_D^3 =\left(\frac{\lambda_D}{a}\right)^3.\]

又因为

\[\lambda_D^2=\frac{\varepsilon_0k_BT}{nq^2}, \qquad n=\frac{3}{4\pi a^3},\]

所以

\[\frac{\lambda_D^2}{a^2} =\frac{1}{3\Gamma}.\]

最终得到

\[\boxed{ N_D=(3\Gamma)^{-3/2} }, \qquad \boxed{ \Gamma=\frac13N_D^{-2/3} }.\]

因此在这一组严格一致的定义下：

\[N_D\gg1 \quad\Longleftrightarrow\quad \Gamma\ll1.\]

5.4 Scaling laws

在固定粒子种类、经典电子定义下：

\[\lambda_D\propto T^{1/2}n^{-1/2}, \qquad \omega_p\propto n^{1/2},\]

\[N_D\propto T^{3/2}n^{-1/2}, \qquad a\propto n^{-1/3}, \qquad \Gamma\propto n^{1/3}T^{-1}.\]

5.5 Density–temperature parameter map

电子耦合参数的实用形式为

\[\boxed{ \Gamma_e =2.32\times10^{-9} \frac{n_e[\mathrm{m^{-3}}]^{1/3}}{T_e[\mathrm{eV}]} }.\]

因此 \(\Gamma_e=1\) 在 \(\log n\)–\(\log T\) 图上是一条斜率 \(1/3\) 的线。线下方更冷、更密，经典强耦合倾向更强；线上方更热、更稀，弱耦合倾向更强。

Worked Example 4 — Is the laboratory plasma weakly coupled?

沿用

\[n_e=1.0\times10^{16}\ \mathrm{m^{-3}}, \quad T_e=2.0\ \mathrm{eV}, \quad \lambda_{De}=1.051\times10^{-4}\ \mathrm m.\]

Step 1 — Debye number

\[N_D=\frac{4\pi}{3}(1.0\times10^{16})(1.051\times10^{-4})^3 =\boxed{4.87\times10^4}.\]

实用公式给出同样结果：

\[\boxed{ N_D =1.72\times10^{12} \frac{T_e[\mathrm{eV}]^{3/2}} {\sqrt{n_e[\mathrm{m^{-3}}]}} }.\]

Step 2 — Wigner–Seitz radius

\[a=\left(\frac{3}{4\pi(1.0\times10^{16})}\right)^{1/3} =\boxed{2.88\times10^{-6}\ \mathrm m}.\]

所以 \(\lambda_D/a\approx36.5\)，一个 Debye 球包含约 \(36.5^3\) 个粒子。

Step 3 — Coupling parameter

\[\Gamma_e =\frac{e^2}{4\pi\varepsilon_0ak_BT_e} =\boxed{2.50\times10^{-4}}.\]

Step 4 — Independent consistency check

\[(3\Gamma_e)^{-3/2}=4.87\times10^4=N_D.\]

5.6 Comparison across plasma environments

下表全部按电子独立 Debye 长度、SI 制和给定的数量级 \(n_e,T_e\) 重新计算；环境参数本身只代表典型尺度，不是固定常数。⁵⁸

6. Lesson synthesis / 全课整合

6.1 One input set, four characteristic outputs

给定一个电子组分的 \((n_e,T_e)\)，本课的四个核心输出相互关联：

6.2 Core equation sheet

6.3 A seven-step diagnostic workflow

拿到一个新系统时，不要先问“它叫什么”，而应按以下顺序量化：

1. Identify species. 列出电子、各价态离子、中性粒子及其 \(n_s,T_s,q_s,m_s\)。
2. Choose the process. 明确系统尺度 \(L\)、过程时间 \(\tau\) 或频率 \(\omega\)；“是不是等离子体”依赖研究尺度。
3. Choose responding species. 高频时离子可能固定；静态慢响应时多个物种都可能参与屏蔽。
4. Compute \(\lambda_D\). 比较 \(L/\lambda_D\)；若不大，保留 Poisson 方程并期待 sheath/space charge。
5. Compute \(\omega_{ps}\). 比较外部驱动、碰撞与损失速率；判断粒子是否来得及建立集体响应。
6. Compute \(N_D\) and \(\Gamma\). 检查连续、弱耦合、线性 Debye 模型是否自洽。
7. State the boundary. 检查强场、非 Maxwell、磁化、碰撞、相对论或量子简并，并用一句话限定结论。

Formula-and-decision table / 公式与决策表

Model boundaries / 模型边界

• Debye–Hückel/Yukawa 推导要求经典、近 Maxwell、静态、小振幅且弱耦合；\(\lvert q\phi\rvert\ll k_BT\) 失效时需非线性 Poisson–Boltzmann 或 kinetic treatment。
• 准中性只适用于 \(L\gg\lambda_D\) 的 bulk；sheath、double layer、probe 附近和高频 charge-separation mode 必须保留 Poisson equation。
• 冷电子 slab 模型给出 \(\omega_{pe}\) 的基本尺度，但真实 dispersion、thermal pressure、magnetic field、collisions 与 boundaries 会改变 mode。
• \(N_D\gg1\) 与 \(\Gamma\ll1\) 检验的是经典连续弱耦合描述；强耦合等离子体并不因此“不再是等离子体”。
• 量子简并、相对论温度、多电荷态和非平衡 distribution 需要相应统计或 kinetic extensions。

7. Misconception clinic / 易错概念诊所

8. Practice and retrieval / 练习与检索

1. [K01 概念] 为什么 ionization fraction \(\alpha\) 不能单独判定一个系统是否适合 plasma description？

2. [K02 符号] 正试探电荷附近，electron Boltzmann factor 为什么含 \(+e\phi/(k_BT_e)\)？

3. [K02 尺度] 若 \(L=20\lambda_D\) 且 \(e\phi/(k_BT_e)\sim0.20\)，估算 \(\lvert\rho\rvert/(en)\)。

4. [K03 函数] 在 \(r=2\lambda_D\) 处，Yukawa potential 与 bare Coulomb potential 的比是多少？

5. [K03 多组分] 单电离 plasma 中若 \(T_i=T_e\) 且 electron、ion 都能静态响应，总 Debye length 与 electron-only \(\lambda_{De}\) 的关系？

6. [K04 频率] 为什么 \(1/\omega_{pe}\) 不是完整 oscillation period？

7. [K04 species] hydrogen plasma 中 \(\omega_{pi}/\omega_{pe}\) 约多少？

8. [K05 scaling] 密度变为 9 倍、温度变为 4 倍时，\(\lambda_D,\omega_p,N_D,\Gamma\) 各乘多少？

9. [K05 边界] “\(N_D=0.2\)，所以这个系统绝不可能是 plasma。”批判此说法。

10. [综合推导] 对 \(\phi=Ae^{-r/\lambda_D}/r\) 求 \(E_r\)，并给出 \(r=\lambda_D\) 时相对 bare Coulomb field 的比。

11. [K02–K04 概念桥接] 宏观准中性的 plasma 为什么仍能支持 electrostatic waves？

12. [K03 尺度判断] 某系统的 \(\lambda_D\) 绝对值很大，能否因此直接判定它不是 plasma？

13. [K05 弱耦合] 为什么 \(N_D\gg1\) 与 \(\Gamma\ll1\) 在匹配定义下描述同一个经典弱耦合极限？

14. [K03–K05 综合定量] 对 \(n_e=4.0\times10^{14}\ \mathrm{m^{-3}}\)、\(T_e=1.0\ \mathrm{eV}\) 的经典 electron plasma，计算 \(\lambda_{De}\)、\(f_{pe}\) 和 \(N_D=(4\pi/3)n_e\lambda_{De}^3\)。

Exit ticket

不用查表，用三句话分别说明 \(\lambda_D\)、\(\omega_{pe}\)、\(N_D\) 与 \(\Gamma\) 检验的物理问题。

9. References and source notes / 参考资料与使用说明

9.1 Source-to-content map

9.2 Bibliographic notes

1. Hutchinson, I. H. “Introduction to Plasma Physics I: Chapter 1 — Introduction.” MIT OpenCourseWare 22.611J, Fall 2003. 本课以此作为定义、准中性与 Debye 推导的主干来源。
2. Fitzpatrick, R. Plasma Physics: An Introduction. University of Texas at Austin online text. 本课采用其 Parameter → Frequency → Shielding → Plasma Parameter 的本科教学逻辑，并显式化解不同 \(\lambda_D\) 约定的差异。
3. Duarte, V. N. “Introduction to Plasma Physics.” Princeton Plasma Physics Laboratory SULI course, 2024. 用于核对当前入门课程对经典三判据与电子薄片模型的呈现。
4. MIT 16.55 course staff. “Fundamentals of Plasma Physics, Lecture 1.” MIT OpenCourseWare, 2014. 用于核对弱电离、碰撞性与准中性例外。
5. Beresnyak, A. 2023 NRL Plasma Formulary. U.S. Naval Research Laboratory. 所有实用系数先按原 cgs 约定读取，再独立转换并代回 SI 基础公式检查。
6. Militzer, B. “The High-Temperature Phase Diagram.” UC Berkeley dissertation web notes. 用于给出 \(\Gamma\) 与量子边界的模型限制，而非把 \(n\)–\(T\) 图误当成普适相图。

9.3 Reading route after this lesson

• 想继续单粒子动力学：Single-Particle Motion。
• 想继续碰撞与平均自由程：Collisions and Transport。
• 想理解 Boltzmann 响应何时被 Vlasov 动力学替代：Kinetic Theory and Vlasov Equation。
• 想理解 \(\omega_{pe}\) 如何成为真实波的色散关系：后续 Plasma Waves and Instabilities 单元。