轫致辐射的理论计算

轫致辐射是等离子体物理中极为重要的一种辐射过程。当带电粒子(通常是电子)在离子或其他带电粒子的库仑场中发生偏转时,由于其速度发生变化,会向外辐射电磁波。在高温等离子体中,轫致辐射是不可避免的能量损失机制,也是诊断等离子体参数的重要手段。本文将系统介绍轫致辐射的理论计算方法,从经典理论出发,过渡到量子力学修正,并给出适用于等离子体辐射功率的实用计算公式。
在非相对论近似下,电子在离子库仑场中的运动可以视为偶极辐射问题。根据经典电动力学,加速运动的带电粒子辐射总功率由拉莫尔公式给出:

$$ P = \frac{e^2 a^2}{6\pi \varepsilon_0 c^3} $$

其中,$e$ 为电子电荷,$a$ 为电子的加速度,$\varepsilon_0$ 为真空介电常数,$c$ 为光速。

在距离为 $r$ 处,电子受离子(电荷为 $Ze$)库仑力产生的加速度为 $a = Ze^2/(4\pi\varepsilon_0 m_e r^2)$。将其代入拉莫尔公式,可以看出辐射功率与距离的四次方成反比,这意味着电子越靠近离子,偏转越剧烈,辐射的能量也越高。

量子力学修正与 Gaunt 因子

经典理论在处理电子与离子的近距离碰撞时存在局限性。当电子的德布罗意波长与碰撞参数相当时,必须引入量子力学修正。为此,理论计算中引入了Gaunt 因子 $g_{ff}$ 来修正经典结果与量子结果之间的偏差。

量子力学给出的单电子-离子碰撞微分截面由贝特-海特勒公式描述,其辐射能量在 $\hbar\omega$ 到 $\hbar\omega + d(\hbar\omega)$ 范围内的截面为:

$$ d\sigma = \frac{16\pi Z^2 e^6}{3\sqrt{3} (4\pi\varepsilon_0)^3 c^3 m_e^2 v^2} \frac{g_{ff}}{\hbar\omega} d(\hbar\omega) $$

其中 $v$ 为电子速度,$\hbar\omega$ 为辐射光子能量。Gaunt 因子 $g_{ff}$ 是一个量级为 1 的无量纲数,通常在 1 到 2 之间,具体数值取决于电子温度和辐射频率。在高温等离子体的粗略估算中,常取 $g_{ff} \approx \sqrt{3}$。

等离子体总辐射功率计算

在实际等离子体物理研究中,我们更关心的是单位体积内所有粒子碰撞产生的总轫致辐射功率,即辐射功率密度。这需要对电子的速度分布(通常为麦克斯韦分布)和所有可能的碰撞参数进行积分。

对于由电子和单电荷离子组成的氢等离子体,假设电子数密度为 $n_e$,离子数密度为 $n_i$,电子温度为 $T_e$,其轫致辐射功率密度 $P_{br}$ 的计算公式为:

$$ P_{br} = 1.69 \times 10^{-32} n_e^2 T_e^{1/2} \quad (\text{W/m}^3) $$

如果等离子体中包含多种离子,或者离子带有多个正电荷 $Z$,则公式需推广为:

$$ P_{br} = 1.69 \times 10^{-32} n_e n_i Z_{eff}^2 T_e^{1/2} \quad (\text{W/m}^3) $$

其中 $Z_{eff}$ 为有效电荷数。从公式中可以得出轫致辐射的两个重要特征:

  • 与密度的平方成正比:辐射功率随密度增加而急剧上升,因为碰撞频率与粒子数密度的乘积相关。
  • 与温度的平方根成正比:温度升高虽然增加了单个电子的动能,但也增加了电子速度,从而减少了电子在离子附近停留的时间(即相互作用时间缩短),综合表现为 $T_e^{1/2}$ 的弱依赖关系。

计算示例:托卡马克等离子体辐射损失

为了更直观地理解轫致辐射的能量损失规模,我们以典型的托卡马克聚变等离子体参数为例进行估算。

假设等离子体参数如下:

  • 电子数密度 $n_e = 1 \times 10^{20} , \text{m}^{-3}$
  • 离子数密度 $n_i = 1 \times 10^{20} , \text{m}^{-3}$
  • 有效电荷数 $Z_{eff} = 1.5$(考虑了少量杂质)
  • 电子温度 $T_e = 10 , \text{keV}$ (换算为开尔文温度约为 $1.16 \times 10^8 , \text{K}$)

将这些参数代入轫致辐射功率密度公式:

$$ P_{br} = 1.69 \times 10^{-32} \times (1 \times 10^{20}) \times (1 \times 10^{20}) \times (1.5)^2 \times (1.16 \times 10^8)^{1/2} $$

计算得:

$$ P_{br} \approx 1.69 \times 10^{-32} \times 10^{40} \times 2.25 \times 1.08 \times 10^4 \approx 4.1 \times 10^5 , \text{W/m}^3 $$

即每立方米等离子体由于轫致辐射将损失约 410 千瓦的功率。在磁约束聚变中,这部分能量以 X 射线形式逃逸出等离子体,是决定聚变点火条件(劳逊判据)的关键因素之一。如果等离子体中混入高 Z 杂质(如钨或钼),由于 $Z_{eff}^2$ 的急剧增大,轫致辐射损失将成倍增加,甚至导致等离子体冷却熄灭。

总结

轫致辐射的理论计算从经典的拉莫尔公式出发,经过量子力学的 Gaunt 因子修正,最终推导出适用于高温等离子体的宏观辐射功率密度公式。掌握轫致辐射的计算方法,不仅有助于评估聚变反应堆中的辐射能量损失,也为利用 X 射线谱线进行等离子体诊断提供了坚实的理论基础。在实际工程与物理分析中,必须严格控制高 Z 杂质的含量,以抑制轫致辐射带来的不利影响。