基于Maxwell–Bloch框架对G358.93?0.03中6.7 GHz甲醇脉泽耀发现象的建模研究

《New Astronomy》：Modelling the 6.7 GHz methanol flaring in G358.93-0.03 using the Maxwell–Bloch framework

　　

研究人员采用Maxwell–Bloch框架对G358.93?0.03中2019年非周期性6.7 GHz甲醇脉泽耀发期间若干代表性速度通道的单天线光变曲线进行建模。将每个速度通道的光变曲线表示为多个有效一维Maxwell–Bloch发射系统强度响应的加权求和，每个系统由表征其局部激发历史的唯象泵浦函数驱动。环境响应参数（T₁=144.4 d、T₂=6.9 d、T_p=1.7 d、L=2.51×10¹⁵ cm）在所有拟合中保持固定，而有效泵浦幅度、到达时间和标度系数因组分而异。结果表明，所选单天线轮廓可由多个有效一维Maxwell–Bloch响应之和再现，耀发组分与超辐射（superradiance） regime下的快瞬变Maxwell–Bloch响应一致；由拟合背景泵浦率推得的初始反转柱密度n₀L?10³–10⁴ cm^?2对应特征超辐射时标T_R,0?0.3–1.4 d，短于所采用的退相与弛豫时标，满足相干瞬变发射的基本时标条件。具有缓慢衰减尾部的组分通过有效长尾泵浦再现，应解释为缓慢变化局部激发的唯象表示而非吸积暴光度曲线的直接测量。单速度通道内使用多个有效样品系由VLBI观测所推断的空间结构化脉泽环境所驱动，多组分拟合可能不唯一，但为在统一Maxwell–Bloch响应框架内描述复杂单天线耀发轮廓提供了有用方式。

研究背景：Class II甲醇脉泽与高质量年轻星体密切关联，由暖尘埃的红外辐射抽运。6.7 GHz跃迁是巨质量恒星形成区监测最广泛的脉泽谱线之一，其光变可灵敏探测局部激发环境的变化。G358.93?0.03（简称G358）是伴随便携吸积暴的非周期性6.7 GHz甲醇耀发的著名实例，2019年1月iMet探测到耀发 onset，流量密度在两周内由约5 Jy增至99 Jy，约两个月内峰值达103 Jy量级，红外观测确认该事件与巨质量年轻星体的吸积暴有关。多历元VLBI观测追踪了热辐射热波在脉泽区的传播，并推断存在开普勒盘与四臂螺旋结构。以往研究虽提出吸积暴提供了增强红外激发的可能来源，但其本身无法决定耀发的精细时间形态；观测到的单天线光变曲线并非简单孤立脉冲，相邻速度通道显示陡升、多峰、次级结构和延展衰减。因此本研究要解决的核心问题是耀发介质如何响应2019年事件相关的激发，而非吸积暴本身的触发机制。Dicke超辐射此前已被提出作为脉泽宿主区中快速瞬变辐射机制，Maxwell–Bloch框架已被用于阐明准稳态脉泽与快瞬变超辐射 regime之间的联系并应用于多个源。研究人员将这一框架应用于G358 2019年耀发的选定6.7 GHz速度分量，检验复杂耀发形态能否被理解为多个有效发射系统对变化激发环境的组合强度响应。

主要关键技术方法：采用一维Maxwell–Bloch方程（MBEs）数值求解，使用四阶Runge–Kutta方法，将每个发射区近似为含二能级甲醇分子的圆柱系统，辐射沿对称轴传播并于推迟时间τ=t?z/c中演化，Fresnel数设为1；观测数据取自Ibaraki 6.7 GHz甲醇脉泽监测数据库中使用Hitachi 32-m单天线望远镜在2019年耀发及后续衰减期间对G358的观测（约MJD 58500至59000），速度分辨率经三通道平滑约0.13 km s^?1，典型1σ rms噪声约0.3 Jy；采用唯象泵浦函数（单脉冲及非对称长尾泵浦）驱动各有效系统，并以多组分强度响应加权求和拟合单天线光变曲线；VLBI图仅用于提供空间背景而非拟合约束。

研究结果：
4.1 接近?18 km s^?1的速度组：该组在相邻通道间包含相似的耀发组分序列，最早窄组分再现陡峭初始上升，较晚组分解释次级结构与延展衰减；在v_LSR=?18.05和?17.95 km s^?1的profile 3中再现出类振铃次级结构，与快瞬变Maxwell–Bloch响应一致，该特征在?17.85 km s^?1变得不明显并与相邻组分混合；较晚组分相对强度亦随速度变化，profile 6在所有情形下先对较宽泵浦组分显示初始超辐射瞬变响应随后进入饱和脉泽 regime。这表明相邻速度通道并非同一耀发的简单缩放，而包含同一组分的不同混合。
4.2 接近?15.6 km s^?1的速度组：该组所需有效组分少于?18 km s^?1组；profile 1和2出现在耀发早期且在该组内相对不变；中心组分profile 3从?15.75至?15.53 km s^?1逐渐占主导，profile 4较弱且在?15.53 km s^?1重要性下降；该组低强度profile 5似对应于?18 km s^?1组的profile 6。这些变化被解释为有效组分相对贡献的改变而非与单个物理发射区的唯一映射。
4.3 时标与模型解释：由拟合背景泵浦率、T₁=144.4 d和L=2.51×10¹⁵ cm推得初始反转柱密度n₀L?10³–10⁴ cm^?2，对应初始特征超辐射时标约T_R,0?0.3–1.4 d，短于所采用的退相与弛豫时标T₂=6.9 d和T₁=144.4 d，故所有拟合组分满足快瞬变超辐射响应的基本时标条件；系数C_v,j将各有效一维系统模型强度标度至观测单天线流量密度，吸收有效发射立体角、聚束及具有相似时间响应的未分辨样品数目的不确定性。
4.4 VLBI背景与局限：使用多个有效系统系由G358的空间复杂性所驱动，VLBI图追踪了耀发 onset后多个历元6.7 GHz甲醇发射的演化，发射区由致密分布演化为更延展结构，被解释为热辐射波传播而非气体物理运动；这些图非拟合数据，但确立单天线速度通道混合多个空间上不同发射区的事实，故拟合的Maxwell–Bloch组分不应被解释为对单个VLBI脉泽斑点的逐一对应。?18和?15.6 km s^?1附近速度组形成空间上不同的簇，可解释光变曲线中不同特征与结构。有效Maxwell–Bloch系统的特征尺度远小于VLBI分辨的空间尺度，一个VLBI脉泽斑点或区域可能包含许多未分辨有效样品而非对应单个Maxwell–Bloch样品。将所采用退相时标与简单碰撞时标估计比较，对τ_coll=T₂=6.9 d、σ_g～10^?16 cm²、v?～10⁵ cm s^?1，得n_H2?1.68×10⁵ cm^?3，处于Cragg et al.（2005）所报6.7 GHz甲醇跃迁10⁴–10⁹ cm^?3密度范围之内但偏向低端。

讨论与结论：研究人员采用一维Maxwell–Bloch框架对G358中选定6.7 GHz甲醇耀发轮廓进行建模，所分析光变曲线对应伴随便携吸积暴的2019年非周期性耀发事件；研究未对吸积暴本身建模，而将泵浦函数唯象处理为有效发射系统所经历的局部激发历史。所选单天线轮廓可表示为有效一维Maxwell–Bloch响应之和，在代表性拟合中环境响应参数保持固定而有效泵浦幅度、到达时间和标度系数因组分而异，该方法使包括陡升、次级结构和延展衰减在内的复杂耀发形态得以按发射气体的多个有效响应加以组织。耀发组分与超辐射 regime下的快瞬变Maxwell–Bloch响应一致，对当前拟合参数集推得的特征超辐射时标短于所采用的弛豫与退相时标，满足相干瞬变发射所需基本条件；具有缓慢衰减尾部的组分通过有效长尾泵浦再现，应解释为缓慢变化局部激发的唯象表示而非吸积暴光度曲线的直接测量。单速度通道内使用多个有效样品系由G358.93?0.03-MM1的VLBI观测所推断的空间结构化脉泽环境所驱动，多组分拟合可能不唯一，但为在统一Maxwell–Bloch响应框架内描述复杂单天线耀发轮廓提供了有用方式。未来建模将受益于同步高 cadence射电与红外监测以及耀发演化期间的空间分辨观测，此类数据有助于更直接约束局部激发历史并检验单个耀发组分能否与发射环境中不同空间区域相关联。