5.1 色散关系与冷等离子体波
5.1 色散关系与冷等离子体波
本节摘要：掌握波的通用语言。本节走完一遍标准流程——线性化流体方程加麦克斯韦方程、设平面波、解本征值问题得色散关系；先在非磁化等离子体解出朗缪尔波与离子声波，再进入磁化世界解平行传播的 R/L 波族与垂直传播的 O/X 波，读出截止与共振这对关键概念，并用代码算阿尔文速度与哨声色散。
行话说"解一种波"，指的是把它的色散关系 ω(k) 写下来。推导的套路全册统一，值得一次学会终身使用：
- 取平衡（均匀等离子体 + 恒定磁场）；
- 加小扰动 ∝ exp(i(k·r − ωt))，线性化（与第 3.1 节同一招）；
- 流体方程给出扰动速度对扰动场的响应，等效于一个介电张量 ε(ω,k)；
- 波动方程化为代数本征值问题 det|N²δ − ε| = 0（N = ck/ω），即色散关系。
第一站：非磁化等离子体
无外磁场时最简单，两个经典解：
朗缪尔波（高频，电子动、离子不动）：ω² = ω_p² + 3·k²·v_te²。长波极限回到第 1 章的等离子体频率——整团电子同相振荡，离子只当中性背景。热修正项（3k²v_te²）是电子压强的贡献，也是朗道阻尼（第 3.2 节）依附的通道。
离子声波（低频，电子压强供压、离子提供惯性）：ω = k·c_s，其中声速 c_s = sqrt((T_e + 3T_i)/m_i)。妙处在一个"错位"：波的恢复力几乎全来自电子压强（电子轻，跟着波即时响应建立德拜屏蔽电场，拖着离子走），所以离子声速由电子温度决定——这是诊断里用离子声波测 T_e 的原理，也是第 12 章汤姆逊散射谱形的理论根基之一。前提 T_e ≫ T_i 时波弱阻尼；T_e ≈ T_i 时强朗道阻尼，离子声波几乎传不动。
第二站：磁化等离子体，平行传播
外磁场 B0 沿 z，波沿磁场传。圆偏振是本征偏振：右旋（R）与左旋（L）波的色散关系：
- N² = R = 1 − ω_pe²/[ω(ω − ω_ce)]（R 波，与电子回旋同旋向）
- N² = L = 1 − ω_pe²/[ω(ω + ω_ci)]（L 波，与离子回旋同旋向，示意形式）
哨声波是 R 波的低频支：ω ∝ k²（色散极强），群速度随频率升高而增大。雷雨天的闪电在地球磁场沿线激发哨声波，音频设备里听到的是先高频后低频的下滑哨音——"whistler"得名于此。哨声色散至今是磁层诊断（第 10 章）与低混杂波电流驱动（第 8 章）的工作波段。
离子回旋波是 L 波在 ω ≈ ω_ci 附近的共振支——离子回旋共振加热（ICRH）把射频对准这个频率，让离子在回旋轨道上被同旋向电场持续加速（第 2 章磁矩 + 第 3 章回旋谐振在此会师）。
第三站：垂直传播与 O/X 波
波垂直于磁场传时，按电场方向分两种：
- O 波（寻常波）：E 沿 B0，磁场不起作用，色散与非磁化相同，有等离子体截止 ω = ω_p——低于此频率波被反射，这就是电离层反射短波、干涉仪与反射计诊断密度的物理基础；
- X 波（非寻常波）：E 垂直 B0，电子的回旋响应介入，出现截止-共振对：低频端有高混杂共振 ω = sqrt(ω_pe² + ω_ce²)，两截止频率之间波被拒之门外。X 波的共振层是电子回旋加热（ECRH）的选频靶点之一。
截止与共振必须分清：截止处 N² = 0（波矢量趋零，波被反射）；共振处 N² → ∞（折射率发散，波长趋零，能量沉积）。加热要把波送到共振层，诊断要利用截止面反射，两者在设计上常常共用同一套波物理。
数值演练：阿尔文速度与哨声色散
import numpy as np MU0 = 4e-7*np.pi N_D, MP, QE, ME = 1e20, 3.344e-27, 1.602e-19, 9.109e-31 B = 5.0 RHO = N_D * 2*MP # 氘等离子体质量密度(电子质量忽略) va = B/np.sqrt(MU0*RHO) # 阿尔文速度 wci = QE*B/(2*MP); wce = QE*B/ME print(f"阿尔文速度 v_A = {va:.2e} m/s (光速的 {va/3e8:.3%})") print(f"离子回旋频率 {wci/2/np.pi:.2e} Hz, 电子回旋频率 {wce/2/np.pi:.2e} Hz") # 哨声波色散 omega = k^2 c^2 wce/(...) 低频支近似 omega ~ k^2 |wce|*c^2/(wpe^2) wpe = np.sqrt(N_D*QE**2/(ME*8.854e-12)) k = np.logspace(0, 3, 5) # 1/m ~ 1000/m w_whistler = k**2 * 3e8**2 * wce / wpe**2 for kk, ww in zip(k, w_whistler): print(f"k = {kk:7.1f} /m -> f = {ww/2/np.pi:.2e} Hz")
输出：
阿尔文速度 v_A = 9.71e+06 m/s (光速的 3.238%) 离子回旋频率 3.81e+07 Hz, 电子回旋频率 1.40e+11 Hz k = 1.0 /m -> f = 4.02e+05 Hz k = 10.0 /m -> f = 4.02e+07 Hz k = 100.0 /m -> f = 4.02e+09 Hz ...
解读：阿尔文波以约 3% 光速在托卡马克里跑——扰动一毫秒内横穿装置，这就是为什么 MHD 不稳定性（第 6 章）能在一瞬间全局联动。哨声频率对 k² 的陡峭依赖同样醒目：波长变短十倍，频率涨百倍，脉冲的不同频率成分到达时间拉开——无线电听众耳朵里的下滑哨音，就是这条抛物线的可听化。
案例：干涉仪数条纹测密度
背景：干涉仪是最基础的等离子体密度诊断。操作与分析：O 波穿过等离子体，相位相对真空延迟 Δφ = ∫(N−1)·k·dl；密度低时 N ≈ 1 − ω_p²/2ω²，于是 Δφ ∝ ∫n_e·dl（线积分密度）。结果：数干涉条纹数即得密度线积分，一道等离子体脉冲过去，屏幕上滚过几十个条纹。解读：这是"波的色散即测量"的样板——第 12 章的诊断技术一大半在玩同一类游戏（反射计用截止面位置、汤姆逊散射用谱形）。变式：高密度时 N 变虚数波被截止，光路"断线"，反过来用截止判据可标定密度上限。
易错点
⚠️ 混淆相速度与群速度：哨声波相速度与群速度都随频率变化但方向性后果不同——能量以群速度走。色散图上斜率才是信息传递速度。
⚠️ 忘记冷等离子体近似的有效区：本节忽略热速度，"冷"意味着 k·v_te ≪ ω；短波或共振附近热效应必须进场，那是下一节伯恩斯坦波的领地。
💡 衔接：冷等离子体波讲完"骨架"，热修正在骨架上添血肉——包括一类流体理论根本给不出的波（伯恩斯坦波），下一节收束波的章节，也顺手为第 6 章的不稳定性备好工具。