1. 快速射电暴与超新星遗迹环境研究概述
快速射电暴（FRB）作为宇宙中最神秘的射电瞬变现象之一，其色散测量（DM）和旋转测量（RM）已成为揭示其源区环境物理特性的关键探针。当FRB信号穿过电离介质时，高频分量比低频分量传播得更快，这种频散效应形成的DM值直接反映了信号路径上的电子柱密度。而RM则揭示了磁场沿视线方向的分量与电子密度的积分关系，为研究FRB源区的磁化环境提供了独特窗口。
近年来，重复暴FRB 20121102和FRB 20190520B的观测显示，它们的DM值呈现出明显的时变特征——前者表现为先升后降的两阶段演化，后者则呈现近乎线性的持续下降趋势。这些现象强烈暗示着FRB源区存在动态演化的致密等离子体环境。理论研究表明，年轻磁星嵌入超新星遗迹（SNR）的模型能够自然解释这类时变特征：随着遗迹膨胀，电子柱密度降低导致DM减小，而激波放大的磁场则贡献显著的RM。
2. 研究方法与数值模拟框架
2.1 前身星模型构建
本研究采用两种典型的前身星演化通道：
- 单星（SS）通道：11M⊙和30M⊙的孤立恒星，通过星风损失部分包层后发生II型超新星爆发
- 双星剥离（BS）通道：通过洛希瓣溢流失绝大部分氢包层的剥离氦星，最终产生Ib/c型超新星
关键参数对比如下表所示：
| 参数 | 11M⊙单星 | 11M⊙双星 | 30M⊙单星 | 30M⊙双星 |
|---|---|---|---|---|
| 抛射物质量(M⊙) | 7.83 | 1.81 | 12.27 | 9.85 |
| 爆炸能量(10^51erg) | 1.0 | 1.0 | 1.0 | 1.0 |
| 预期超新星类型 | IIP-like | Ib/c-like | IIb-like | Ib/c-like |
2.2 超新星遗迹模拟设置
采用一维辐射流体力学代码结合非平衡电离(NEI)计算，自洽追踪激波区热力学状态和电离演化。模拟从爆发后3年开始，持续追踪500年演化，时间步长采用自适应策略：
- 3-300年：1年分辨率
- 300-500年：4年分辨率
抛射物密度剖面采用分段幂律分布： $$ ρ_{ej}(r) = \begin{cases} ρ_0 & r \leq R_{core} \ ρ_0(r/R_{core})^{-n} & r > R_{core} \end{cases} $$ 其中n=11反映致密前身星特征，核心半径$R_{core}$由抛射物质量$M_{ej}$和爆炸能量$E_{SN}$决定。
3. 色散测量(DM)的演化特征
3.1 激波区贡献分析
在仅考虑激波区（反向激波到前向激波之间）的情况下，DM演化呈现以下特点：
- 单星模型：峰值DM可达8.3 pc cm⁻³（11M⊙）和5.1 pc cm⁻³（30M⊙），出现在爆发后约10年
- 双星模型：DM始终低于0.1 pc cm⁻³，反映剥离前身星抛射物质量较小的特性
激波区DM的时间演化近似遵循幂律衰减$DM \propto t^{-α}$，其中α≈1.0-1.5。这种衰减主要源于：
- 抛射物膨胀导致的电子密度降低（$n_e \propto t^{-3}$）
- 积分路径长度增长（$\Delta R \propto t$）的部分抵消效应
3.2 全电离区贡献分析
当计入未受扰动抛射物和星际介质的贡献时，DM演化呈现更丰富的特征：
3.2.1 主导机制转变
早期抛射物主导阶段（t≲100年）：
- DM以未受扰动抛射物贡献为主
- 遵循$DM \propto t^{-2}$的陡峭衰减
- 11M⊙单星模型早期DM高于30M⊙模型（尽管后者抛射物质量更大），这是由于：
- 30M⊙抛射物的平均分子量更高（氢包层较薄）
- 11M⊙模型具有更高的特征膨胀速度
晚期CSM主导阶段：
- 星际介质贡献逐渐显现
- 衰减斜率趋缓至$DM \propto t^{-1}$以下
- 单双星模型的DM排序可能反转
3.2.2 电离度参数化影响
未受扰动区域的电离分数$χ_e$对DM有显著影响。我们采用以下参数范围：
- 未受扰动抛射物：$0.01 ≤ χ_{e,unej} ≤ 1$
- 星际介质：$10^{-4} ≤ χ_{e,ISM} ≤ 1$
基准模型（HH）取$χ_e=0.1$时，各模型在典型时标的DM值如下表：
| 模型 | t=10年 | t=30年 | t=100年 | t=300年 |
|---|---|---|---|---|
| 11M⊙单星 | 450 | 180 | 65 | 25 |
| 11M⊙双星 | 120 | 45 | 15 | 6 |
| 30M⊙单星 | 380 | 200 | 90 | 40 |
| 30M⊙双星 | 350 | 170 | 75 | 35 |
3.3 与观测暴的对比分析
将模拟结果与FRB 20190520B（dDM/dt=-12.4 pc cm⁻³ yr⁻¹）和FRB 20121102（dDM/dt=-3.93 pc cm⁻³ yr⁻¹）匹配，得到关键时标：
FRB 20190520B：
- 11M⊙单星：t≈27年，DM≈184 pc cm⁻³
- 30M⊙双星：t≈23年，DM≈168 pc cm⁻³
FRB 20121102：
- 11M⊙单星：t≈40年，DM≈89 pc cm⁻³
- 30M⊙单星：t≈42年，DM≈98 pc cm⁻³
二阶导数分析显示，DM衰减曲率在早期阶段约为1-4 pc cm⁻³ yr⁻²，与观测估计值相符。
4. 旋转测量(RM)与磁场结构
4.1 计算方法
RM由以下积分决定： $$ RM \simeq 0.81 \int n_e B_\parallel dl \quad (\text{rad m}^{-2}) $$ 其中磁场强度通过激波区压强分数$ε_B$参数化： $$ B = \sqrt{4πε_Bρv^2} $$ 基准模型取$ε_B=0.1$，覆盖典型超新星遗迹中湍流放大磁场的能量占比（1%-30%）。
4.2 与FRB 20121102的匹配
通过定义对数斜率β≡dlogRM/dlogt，将模拟RM演化与观测匹配。对于FRB 20121102的后峰衰减阶段（β≈-0.36），发现：
- 11M⊙单星模型在t≈9.54年时最佳匹配观测斜率
- 30M⊙单星模型需t≈9.99年
这表明虽然激波区对DM贡献有限，但其放大的磁场足以产生观测到的RM演化特征。
5. GHz辐射逃逸与光学深度
5.1 不透明度来源
总光学深度包含两项贡献： $$ τ_{tot} = τ_{es} + τ_{ff} $$
- 汤姆孙散射：早期主导（t≲10年） $$ τ_{es} = \int n_eσ_T dl $$
- 自由-自由吸收：晚期主导 $$ τ_{ff}(ν) = 1.9\times10^{-2} T_e^{-3/2}ν^{-2}n_e \langle g_{ff} \rangle \langle Z^2 \rangle dl $$
5.2 逃逸时标分析
当$τ_{tot}≤1$时，GHz频段FRB辐射可逃逸。全区域计算显示：
- 激波区：始终透明（$τ_{tot}≪1$）
- 全电离区：逃逸时标$t_{esc}$强烈依赖未受扰动区电离度
- 高电离（$χ_e=0.1$）：$t_{esc}≲70$年
- 低电离（$χ_e=0.01$）：早期即可透明
这一结果支持年轻SNR作为FRB源区环境的合理性，同时解释了为何仅部分FRB表现出明显DM时变特征。
6. 讨论与展望
6.1 模型局限性
当前1D模拟存在以下简化：
- 未考虑磁星诞生 kick（≲500年内位移可忽略）
- 忽略脉冲星风星云(PWN)对电离度的额外贡献
- 磁场处理为参数化形式，未自洽演化
6.2 未来方向
- 多维模拟：研究流体不稳定性导致的混合效应
- 辐射转移：自洽计算电离平衡
- 多信使约束：结合X射线/光学观测校验模型
这些改进将有助于更精确地利用FRB作为探针，研究极端天体物理环境的等离子体特性。
