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.831.8112.279.85爆炸能量(10^51erg)1.01.01.01.0预期超新星类型IIP-likeIb/c-likeIIb-likeIb/c-like2.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} $$ 其中n11反映致密前身星特征核心半径$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取$χ_e0.1$时各模型在典型时标的DM值如下表模型t10年t30年t100年t300年11M⊙单星450180652511M⊙双星1204515630M⊙单星380200904030M⊙双星35017075353.3 与观测暴的对比分析将模拟结果与FRB 20190520BdDM/dt-12.4 pc cm⁻³ yr⁻¹和FRB 20121102dDM/dt-3.93 pc cm⁻³ yr⁻¹匹配得到关键时标FRB 20190520B11M⊙单星t≈27年DM≈184 pc cm⁻³30M⊙双星t≈23年DM≈168 pc cm⁻³FRB 2012110211M⊙单星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} $$ 基准模型取$ε_B0.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}$强烈依赖未受扰动区电离度高电离$χ_e0.1$$t_{esc}≲70$年低电离$χ_e0.01$早期即可透明这一结果支持年轻SNR作为FRB源区环境的合理性同时解释了为何仅部分FRB表现出明显DM时变特征。6. 讨论与展望6.1 模型局限性当前1D模拟存在以下简化未考虑磁星诞生 kick≲500年内位移可忽略忽略脉冲星风星云(PWN)对电离度的额外贡献磁场处理为参数化形式未自洽演化6.2 未来方向多维模拟研究流体不稳定性导致的混合效应辐射转移自洽计算电离平衡多信使约束结合X射线/光学观测校验模型这些改进将有助于更精确地利用FRB作为探针研究极端天体物理环境的等离子体特性。