[论文解读] Kinetic theory of cosmic ray and gamma-ray production in supernova remnants expanding into wind bubbles

Evgeny Berezhko, H. J. Voêlk|arXiv (Cornell University)|Feb 22, 2000

Astrophysics and Cosmic Phenomena被引用 11

一句话总结

该论文构建了一个在恒星风吹拂形成的泡状结构中膨胀的超新星遗迹（SNR）的宇宙射线（CR）与伽马射线产生动力学模型，考虑了非线性CR反作用和激波特性的效应。研究发现，CR可有效吸收20%至40%的爆发能量，产生硬幂律谱（γ ≈ 2.0–2.1），最大能量可达约10¹⁴ eV；而在低密度风泡中，伽马射线通量显著被抑制，除非前身星为红超巨星（SN II），此时在距离达30 kpc范围内，TeV辐射可能在数百年内可被探测到。

ABSTRACT

Abbreviated Abstract: A kinetic model of particle acceleration in supernova remnants (SNRs) is extended to study the cosmic ray (CR) and associated high energy gamma-ray production during SN shock propagation through the inhomogeneous circumstellar medium of a progenitor star that emits a strong wind. The wind forms a low-density bubble surrounded by a swept-up shell of interstellar matter. Gamma rays are produced as a result of decay of pions which in turn are the result of collisions of CRs with nuclei of the thermal plasma. The time evolution of the SNRs is followed numerically, taking into account the nonlinear backreaction of the accelerated CRs. Examples typical for SN type Ib and SN type II explosions are considered. Apart from the confirmation of the known fact that acceleration is extremely rapid and that the upper momentum cutoff is reached almost immediately after the explosion due to the high wind magnetic field, it is also shown that the CRs are accelerated with high efficiency. The CR energy spectrum, ultimately produced in the SNRs, has a power law form with index = 2.0 to 2.1 in a wide energy range up to a maximum energy which is about 10^14 eV if the CR diffusion coefficient is as small as the Bohm limiting value. However, the expected pi^0-decay gamma-ray flux is considerably lower than in the case of a uniform interstellar medium; a relatively high gamma-ray luminosity in the TeV band, detectable at distances of several kpc, is only expected in the case of a relatively dense ISM with a number density above 10/cc.

研究动机与目标

• 模拟超新星遗迹（SNR）在前身星恒星风形成的非均匀周介质中膨胀时的宇宙射线（CR）加速及其相关伽马射线发射。
• 研究低密度风泡与致密壳层的存在如何通过pion衰变影响CR加速效率与伽马射线产生。
• 考察非线性CR反作用、激波特性的改变以及磁流体力学波耗散在塑造CR谱与伽马射线光度中的作用。
• 比较SN II型（红超巨星前身星）与SN Ib型（沃尔夫-拉叶星前身星）在可探测高能伽马射线发射方面的预测。

提出的方法

• 采用动力学方法求解包含能量依赖扩散系数的宇宙射线分布函数输运方程，结合扩散性激波特加速机制。
• 模型自洽地处理通过CR反作用引起的激波特性改变，包含时间依赖的激波速度与演化中的预区结构。
• 预区中磁流体力学波的耗散加热热等离子体，使超热粒子能够被注入加速过程。
• 根据CR与壳层中热等离子体核素发生的非弹性核碰撞后产生的π⁰衰变，计算伽马射线通量。
• 对系统方程应用相似性解法，假设激波速度恒定，以在相似变量中导出时间无关的分布形态。
• 数值模拟追踪SNR随时间的演化，包括自由膨胀相到绝热相的过渡，重点关注CR能量预算与谱形变化。

实验结果

研究问题

• RQ1 风吹拂泡的存在如何影响SNR中宇宙射线加速的效率与谱形？
• RQ2 非线性CR反作用与激波特性的改变在决定CR最大能量与谱形中起何作用？
• RQ3 与均匀星际介质相比，风吹拂泡中的π⁰衰变伽马射线光度有何差异？
• RQ4 在何种条件下，风泡中的SNR可产生可探测的TeV伽马射线辐射，持续时间多长？
• RQ5 CR谱与成分如何随能量与时间演化，其特征如何反映前身星风的化学组成？

主要发现

• 当注入分数η = 10⁻³时，宇宙射线在演化过程中可吸收20%至40%的超新星爆发能量，具体取决于周介质参数。
• CR动量谱在宽广能量范围内呈幂律分布，谱指数γ = 2.0至2.1，若扩散系数处于玻姆极限，则最大能量可达约10¹⁴ eV。
• 与均匀星际介质相比，风吹拂泡中的π⁰衰变伽马射线通量显著降低，尤其在低密度星际介质中更为明显。
• 对于SN II型（红超巨星前身星），在1 kpc距离处，TeV伽马射线通量在约300年内可超过10⁻⁹ cm⁻² s⁻¹，使得此类SNR在当前仪器灵敏度下，最远可达30 kpc范围内被探测到。
• 对于SN Ib型，若星际介质密度低于0.3 cm⁻³，则TeV伽马射线通量始终低于10⁻¹¹ cm⁻² s⁻¹，这解释了高阈值阵列未探测到的原因。
• 在玻姆扩散条件下，伽马射线谱的截止能量约为10¹³ eV，低于许多当前伽马射线望远镜的探测阈值，部分解释了未被探测到的现象。