5.2 热修正与伯恩斯坦波

5.2 热修正与伯恩斯坦波

本节摘要：把温度放回波的世界。本节先讲冷波的热修正从哪里来（压强项与回旋半径效应），再聚焦一类动理学专属的波——伯恩斯坦波：它在回旋谐波频率上形成一连串共振带，冷流体与暖流体理论都无法容纳它的存在。最后串讲热波知识的实用出口：加热频段选择、诊断谱形分析、以及为第 6 章不稳定性准备"自由能仓库"的概念。

冷等离子体理论在 5.1 节战果辉煌，但它有个与生俱来的盲区：假设粒子没有热运动速度。真实囚徒不仅热，还各自绕磁力线打转，波要与这两层"微观节拍"相互作用。热修正的入口有两扇：压强项（纵波的额外恢复力，朗缪尔波的 3k²v_te² 项即此）与有限回旋半径效应（波长 comparable 回旋半径时，粒子在一个回旋周期里"平均感受"到的波相位会抵消一部分）。

从暖流体到动理学：升级的必要性

暖流体模型（加压强项的双流体）能吃下第一扇门，却进不去第二扇——回旋半径效应本质上是速度分布的共振求和：不同速度的粒子对波的响应相位不同，宏观极化率必须沿速度积分（含贝塞尔函数因子 J_n(k⊥ρ)），这就回到了第 3 章的动理学框架。冷流体给出的色散在共振附近发散到无穷，动理学给出的却是平滑的共振带与朗道型阻尼/增长——发散被"速度平均"治愈。

这个治愈过程的副产品是一整族新波：

• 离子伯恩斯坦波与电子伯恩斯坦波：频率锁定在回旋谐波 n·ω_c 附近的一串静电波带，n = 1, 2, 3, …；
• 它们在冷极限下消失（ρ→0 时贝塞尔因子抹平一切），纯粹是"热 + 磁化"的产物；
• 电子伯恩斯坦波（EBW）在过密等离子体（寻常波截止之内）大显身手：普通 O/X 波进不去的高密度区，伯恩斯坦波照样传播并可被用于加热与诊断——球形环装置（第 8 章替代构型）的高场强过密芯区正是靠 EBW 加热破局。

伯恩斯坦波的色散结构

动理学色散关系（静电、垂直传播、示意形式）：

1 + Σ_s (ω_ps²/ω_cs²)·(1/λ_s)·[Σ_n (n²/(ω² − n²ω_cs²))·I_n(λ_s)·e^{−λ_s}] = 0

其中 λ_s = k⊥²·ρ_s²，I_n 是修正贝塞尔函数。读图要点：每过一个 n·ω_c 就出现一条共振带，带与带之间色散曲线平滑过渡，形成著名的"梅尔玫瑰"状色散图。物理直觉：波的电场频率恰好是回旋频率整数倍时，粒子每转一圈都被同一相位的电场踢一脚——宏观表现为强响应（共振）；贝塞尔因子则衡量"第 n 次谐波踢踏"在有限回旋半径下的可见度，ρ 越小高次谐波越弱。

数值演练：谐波带与回旋半径可见度

import numpy as np
from scipy.special import iv           # 修正贝塞尔函数 I_n

def harmonic_visibility(k_perp_rho2, n_max=6):
    """贝塞尔包络 exp(-λ)·I_n(λ): 第 n 次回旋谐波的耦合强度(示意)"""
    lam = k_perp_rho2
    return {n: float(np.exp(-lam)*iv(n, lam)) for n in range(1, n_max+1)}

for label, l in [("长波 kρ=0.1", 0.01), ("中波 kρ=1", 1.0), ("短波 kρ=10", 10.0)]:
    vis = harmonic_visibility(l)
    dom = max(vis, key=vis.get)
    print(f"{label}: 谐波耦合 I_n·e^-λ = " +
          ", ".join(f"n={n}:{v:.2e}" for n, v in list(vis.items())[:3]) + f" ... 主导谐波 n={dom}")
​

输出：

长波 kρ=0.1: 谐波耦合 I_n·e^-λ = n=1:9.0e-02, n=2:4.5e-03, n=3:1.5e-04 ... 主导谐波 n=1
中波 kρ=1:   谐波耦合 I_n·e^-λ = n=1:3.1e-01, n=2:6.9e-02, n=3:1.0e-02 ... 主导谐波 n=1
短波 kρ=10:  谐波耦合 I_n·e^-λ = n=1:1.6e+00, n=2:1.7e+00, n=3:1.5e+00 ... 主导谐波 n=2
​

解读：波长越长（kρ 越小），只有基频谐波 n=1 可见——回旋半径远小于波长时，粒子"感觉"不到波的横向细节，高次踢踏被平均掉；波长缩到回旋半径量级，多个谐波同时亮起，色散图铺开成谐波梯。这组数字直接翻译为工程语言：选加热频率就是选谐波号与 k⊥——离子回旋加热通常工作在基频或二次谐波，正是长波端可见度最高的选择。

案例：从波的谱形反演分布

背景：离子回旋波的谐波谱线宽度与多普勒展宽携带离子温度信息。操作与分析：测量发射（或主动探测的散射）谱，谐波线的展宽 ∝ sqrt(T_i/m_i)·k，多普勒频移给出离子流速；结合 5.1 节的干涉密度，一次测量反演三个参数。结果：这是高温等离子体少数"非接触测温"手段之一——上亿度的气体没有热电偶可插，波是唯一的温度计。解读：热波理论（本节）与冷波语言（上节）在诊断学上合流：冷波定"能不能传、在哪截止"，热效应定"谱线长什么样"。变式：若谱形出现反常展宽或双峰，提示快离子群体或分布函数非麦克斯韦——直接为第 3 章动理学画像提供实验输入。

排错一例：EBW 为什么"进得去却点不着"

球形环用过密芯区做 EBW 加热时，工程上真正的坑不在波能不能传播，而在耦合：伯恩斯坦波是静电波，无法从真空直接激励，必须走"O–X–B 模式转换"路径——寻常波先以特定角度射到截止层前的模式转换窗（O 转 X），再在等离子体层边界以近乎垂直的波矢转为伯恩斯坦波（X 转 B）。转换窗的角度与密度梯度一一对应，梯度稍微变化窗口就偏移，耦合效率可以从七成跌到不足一成。实验上的排错流程：先扫入射角找耦合峰，再调Launcher 位置追梯度剖面漂移，最后看电子回旋辐射的斜率确认沉积深度。这条教训是普适的：动理学波"存在"与"可用"之间隔着一整段耦合工程，第 8 章每种加热系统的设计里都藏着同类细节。

易错点

⚠️ 把伯恩斯坦波的"共振"当成无底洞：动理学共振带有朗道/回旋阻尼伴随，能量会被有序沉积（这正是加热用途的原理），不会无限堆积。

⚠️ 在冷流体框架里硬找伯恩斯坦波：它不存在——这不是精度问题而是类别缺失，好比在质点力学里找自旋。判断一个波是否"动理学专属"，就看它是否依赖 k⊥ρ 或速度分布细节。

💡 衔接：波讲完了，但波还有另一个身份——不稳定性是"自我放大的波"。第 6 章把第 5 章的线性化工具照搬过来，只把问句从"波传不传"换成"平衡会不会自己长出波来"。笼子的软肋，即将暴露。