推导
在均匀等离子体附近线性化
写 f=f0(v)+f1,f0 是均匀分布,f1 很小。只保留 f1 和电场的一阶项。E∂f1/∂v 是二阶小量,可以略去,剩下一个线性方程,驱动它的是电场作用在背景分布的斜率上。
设模式 f1,E∝ei(kx−ωt)。Vlasov 方程用 E 表出 f1,泊松方程又把 f1 的密度变回 E。要有非零电场,介电函数必须为零。
波速处的极点
被积函数在电子速度等于波速 v=ω/k 处有极点。朗道指出,初值问题决定了怎样绕过它:从复 ω 平面的上半平面延拓积分,积分路径从极点下方经过。这样会多出一项,正比于波速处的斜率 f0′。
对麦克斯韦分布,这个积分用等离子体色散函数 Z 表示。程序通过 Faddeeva 函数计算 Z,再用牛顿法找根。式中 ζ=ω/(2k);根落在下半平面,波就衰减。
完整例题
例题:三个波数
解 ε=0:k=0.3 时 ω=1.1598−0.0126i,k=0.5 时 1.4157−0.1534i,k=0.8 时 1.7999−0.5346i。短波衰减得快得多,因为它的相速度更接近热速度,能跟着波走的电子多得多。
Lab 只用电场的历史就测出同样的速率:在 t=4 之后、非线性俘获或网格复现可能出现之前,对电场峰值的对数做直线拟合。拟合结果是 −0.0127、−0.1536、−0.5358,实部也吻合到 0.2% 以内。
长波时,一个简短的公式就能抓住大部分结果。k=0.5 时它给出 γ≈−0.151,接近真根。它适合估算,不适合做对照:换个波数,误差会很快变大。
这个速率说明了什么
阻尼率描述的是电场和密度,也就是对速度平均之后的量。分布本身并没有弛豫:能量离开电场,以越来越细的速度条纹留在电子里。这个速率还假定波一直很小。波强到能抓住它想交出能量的那些电子时,情况会变,见第 5 章。
带着问题动手
先预测,再改变,最后解释
- 用上面的三个根,预测 k=0.8 的电场衰减到 e−5 需要多少时间单位。
- 在相空间 Lab 里把波数滑块从 0.3 拉到 0.8,对照“测得的衰减率”和“介电函数根给出的衰减率”。
- 解释为什么把波纹深度调大一点时仍吻合,调大很多时就对不上了。
三个相空间 Lab 用等离子体单位求解一维电子模型,每个 Lab 都能导出参数、模型版本和电场历史。
试着自己解释
1. k=0.8 时,电场衰减 e5 倍要多久?
展开答案与解释
衰减率 0.5346,时间为 5/0.5346≈9.4 个单位,比两个振荡周期(每个 2π/1.80≈3.5)略长。
2. 线性计算为什么能丢掉 E∂f1/∂v?什么时候这样做不安全?
展开答案与解释
它是两个小量的乘积。电场强到能在衰减时间之内弯折电子轨道时,这一项就要紧了;第 5 章用弹跳时间来衡量这一点。
3. 如果分布在波速处的斜率为正,根的虚部应是什么符号?
展开答案与解释
为正:斜率项变号,波转而增长。第 5 章用尾部鼓包检验这一点。