1.2 统计描述与德拜球
1.2 统计描述与德拜球
本节摘要：一立方厘米的托卡马克等离子体里有约 10^14 个粒子，逐个追踪既不可能也不必要。本节建立分布函数与速度空间的统计语言，并用统计观点重新审视德拜球——为什么"德拜球内粒子数足够多"是集体行为的判据，以及麦克斯韦分布、温度、压强这些热力学量如何从统计描述中自然涌现。
驯服上亿个囚徒，逐个点名是点不过来的。可行的办法是人口普查式的：不问"3 号电子现在在哪"，而问"速度在某个区间内的电子占多大比例"。这个比例函数就是分布函数，它是第 3 章动理学理论的主角，也是本章把热力学语言接入等离子体物理的桥梁。
分布函数：从点名到普查
定义 f(r, v, t) 使得 f·d³r·d³v 是时刻 t 位于位置 r 附近体积元、速度在 v 附近速度体积元中的粒子数。六个坐标（三个位置 + 三个速度）构成相空间。几个基本量由 f 的积分给出：
- 密度：n(r,t) = ∫ f d³v；
- 平均速度：u = (1/n)∫ v f d³v；
- 温度与压强：由速度对平均值的二阶矩决定，对麦克斯韦分布有 p = n·T。
这种描述的威力在于降维：10^14 个粒子的 6N 维问题，压成一个六维函数。代价是丢掉了粒子间的关联信息——这正是第 1 节德拜球判据要守的门：只有当德拜球内粒子数 N_D 足够大、每个粒子同时与大量同伴发生弱作用时，两两关联才不重要，单粒子分布函数才是充分的描述。N_D 太小的系统（如某些量子点、尘埃等离子体的极端参数）必须用更精细的多体或分子动力学方法。
麦克斯韦分布与"温度"的含义
碰撞足够的体系会弛豫到麦克斯韦速度分布。需要特别强调等离子体语境的两个要点。
其一，电子与离子各有各的温度。库仑碰撞中，同种粒子交换动量快、异种交换慢（电子-离子质量比悬殊，一次碰撞只能传递很小份额的能量），于是典型放电等离子体中电子温度可以是离子温度的一百倍——这直接决定了第 11 章低温等离子体"电子温度高、气体近室温"的工艺特性。
其二，磁化等离子体是各向异性的：沿磁力线与垂直磁力线的温度可以不同，压强变成并矢 p∥ 与 p⊥。这个各向异性是不稳定性（第 6 章的火蛇模等）的自由能来源之一。
数值演练：从分布函数提取宏观量
import numpy as np # 用 5 万个抽样粒子代替 1e14 个真实粒子做统计 rng = np.random.default_rng(42) T_eV, m_ratio = 10.0, 1.0 # 一维麦克斯韦(高斯)抽样：速度标准差 sigma = sqrt(T/m) def sample(T_eV, mass_in_eV=1.0, N=50_000): sigma = np.sqrt(T_eV / mass_in_eV) # 单位: sqrt(eV/eV)~光速无关的自然单位 return rng.normal(0.0, sigma, N) v = sample(T_eV) v2_mean = np.mean(v**2) T_from_moment = m_ratio * v2_mean # 一维: T = m<v^2> print(f"输入 T = {T_eV} eV, 二阶矩反推 T = {T_from_moment:.2f} eV") # 输出: 输入 T = 10.0 eV, 二阶矩反推 T = 10.02 eV (统计涨落 ~1/sqrt(N)) # 演示电子/离子温度悬殊的后果：能量弛豫时间与质量的联系 for r in (100, 1836, 3672): # 模拟质量比 frac = 4*m_ratio*r/(1+r)**2 # 单次碰撞能量转移份额的上限系数 print(f"质量比 {r:5d}: 单次碰撞最大能量转移系数 ~ {frac:.3f}")
输出示例：
输入 T = 10.0 eV, 二阶矩反推 T = 10.02 eV 质量比 100: 单次碰撞最大能量转移系数 ~ 3.846 质量比 1836: 单次碰撞最大能量转移系数 ~ 2.174 质量比 3672: 单次碰撞最大能量转移系数 ~ 1.086
注意最后一行的系数归一化方式不同质量比不可直接比大小，真正要看的结论是：电子对离子的能量弛豫速率正比于电子质量与离子质量之比，对氘离子这个比约 1/3672，所以电子温度与离子温度可以长时间脱钩——这就是"双温度"不是近似玄学而是物理现实的定量根源。
用统计观点重看德拜球
静态推导德拜长度用的是玻尔兹曼关系：电子密度在势场中按 n ∝ exp(eφ/T) 分布，与泊松方程联立、线性化后得到屏蔽势 φ ∝ exp(-r/λ_D)/r。这里值得停下来体会统计观点带来的三个新认识：
- 屏蔽是热运动的功劳：若 T→0，λ_D→0，屏蔽层无限薄——冷电子会被势阱吸死，形成点状的经典"电子跪"。有限温度让电子有动能"讨价还价"，屏蔽层才有有限厚度；
- 屏蔽是线性的前提是 eφ ≪ T：势能远小于热能时玻尔兹曼因子才能线性化。在强场区（鞘层、激光近临界等离子体）线性屏蔽失效，需要泊松方程全量求解——第 11 章的鞘层物理会正面遇到这一点；
- N_D 判据是"统计描述合法性"判据：N_D 大意味着一个粒子感受到的是大量同伴的平均场而非少数几个的近距离遭遇，平均场近似（第 3 章弗拉索夫方程）因此成立。可以算出库仑对数的典型值 ln Λ ≈ ln(λ_D/b_min) ≈ 10–20，这个对数在第 7 章输运系数里反复出现。
案例展开：给一份电离层档案补全参数
背景：电离层测站报告 F 区电子密度 1e12 m⁻³、电子温度 0.2 eV（约 2300 K）。操作：按第 1.1 节公式算得 λ_D ≈ 3.3 mm，德拜球粒子数 N_D ≈ 1.5e5，远超判据；德拜球体积相对探测体积（公里级）完全可忽略。结果：准中性极好，可用流体或动理学描述。解读：在这个参数区，"离子质量大、跟不上电子热运动"使得屏蔽几乎全由电子完成——德拜长度公式里也确实只出现电子温度与密度。变式：若在极区等离子体密度坑（n 掉到 1e9 m⁻³）处，λ_D 增大到约 10 cm，N_D 掉到 1e6 以下两个量级，但判据仍通过；真正让判据失效的是尘埃等离子体——毫米级尘埃带上千个元电荷，N_D 可以降到个位数，统计描述崩塌，必须回到"颗粒物理"。
易错点
⚠️ 把电子温度 10 eV 读成"10 度"或把 eV 直接当 K：差 11600 倍，所有后续估算全错。
⚠️ 以为分布函数必须是麦克斯韦的：驱动中的等离子体（束流、射频加热、逃逸电子尾）分布可以严重非麦克斯韦，此时"温度"只是二阶矩的代名词，第 3 章的朗道阻尼恰恰依赖非麦克斯韦分布的细节。
💡 衔接提示：本节的 f(r,v,t) 是静态快照，第 3 章要问的是它如何随时间演化——答案是动理学方程组。而在那之前，第 2 章先看最简单的情形：单个粒子在给定场里的运动，那是所有统计图像的"原子事实"。