利用斯皮策太空望远镜研究M型米拉变星大气中CO2的相位依赖性中红外光谱分析

本刊推荐：为揭示脉动对米拉变星周星环境的影响，研究人员利用斯皮策太空望远镜的红外光谱仪（IRS）对13颗M型米拉变星进行中红外光谱观测，结合帕洛马测试床干涉仪（PTI）的角直径测量，通过RADEX辐射转移模型首次定位CO2气体发射区。研究发现CO2分布高度延伸且温度/柱密度跨度大，其14.98μm特征在吸收/发射间的转换揭示了"制冷区"的存在，为星尘形成机制提供新视角。该成果发表于《皇家天文学会月报》。

在浩瀚宇宙中，老年恒星如同星际化学工厂，通过星风将富含新元素的物质抛洒至星际空间。特别是渐近巨星支（AGB）恒星，作为低中质量恒星演化的终末阶段，被认为是星系中尘埃颗粒的主要贡献者。其中脉动周期为100-500天的M型米拉变星，因其动态大气中丰富的分子化学过程而备受关注。然而，由于地球大气层对特定红外波段的吸收，地面望远镜难以观测这些恒星的关键分子特征，阻碍了对其化学演化机制的深入理解。

早期红外空间天文台（ISO）的观测虽首次在富氧AGB星中探测到二氧化碳（CO₂），但其光谱特征的复杂变化行为（如14.98μm波段在吸收与发射间转换）始终未得到系统解释。传统理论认为碳原子主要被锁定在CO分子中，且基于热平衡假设的模型无法解释脉动对化学环境的影响。Duari等人提出CO₂可能通过激波化学在光球层附近形成，但需要更多时间分辨观测验证这一假说。

为解开这些谜团，Baylis-Aguirre等人利用斯皮策太空望远镜红外光谱仪（IRS）对13颗M型米拉变星进行相位跟踪观测，结合帕洛马测试床干涉仪（PTI）的K波段干涉测量数据，首次实现分子层物理条件与空间定位的同步分析。所有样本均显示13.87μm、14.98μm和16.18μm处CO₂的未分辨转动振动Q支谱线，其中14.98μm特征表现出显著的相位依赖性变化。

研究团队创新性地构建包含8个振动能级（最高J=50转动能级）的CO₂分子文件，集成818条辐射跃迁和20,000个碰撞速率系数，利用RADEX辐射转移代码建立双层模型。通过将发射/吸收层计算叠加，同时拟合三个主要波段，并结合PTI测量的角直径反演发射区位置。关键技术方法包括：斯皮策IRS中红外光谱观测（10-37μm）、帕洛马测试床干涉仪角直径测量、RADEX非局部热动平衡辐射转移建模、以及针对13μm特征的高斯拟合分析。

3.1 CO₂特征

样本光谱显示出13.87μm、14.98μm和16.18μm处明显的Q支谱线，其中14.98μm基频v₂带在吸收与发射间动态转换，而13.87μm/16.18μm作为费米简并对始终呈发射态。以RT Cyg为例，其七次观测覆盖近完整脉动周期，前五次观测显示14.98μm为发射态，后两次转为吸收态，对应发射区从12R_★收缩至3-5R_★，表明CO₂分布跨越恒星大气层不同区域。

3.2 13μm特征

七颗米拉变星呈现13μm宽发射特征，其形态随相位变化。如S Peg在φ=0.631时特征尖锐（线心13.11μm），至φ=0.098时展宽至0.61μm且线心蓝移。高斯拟合显示该特征与CO₂谱线强度存在协同变化，暗示尘埃与气体间的相互作用。

4.1 RADEX建模结果

模型显示吸收组恒星（如R Tri、S Peg）的14.98μm层对应600-1000K高温区，位于2-5R_★；发射组（如SS Cas、W Her）则该层温度降至400-500K，延伸至5-15R_★。当14.98μm特征呈吸收态时，其激发温度较动力学温度高约50K，表明存在非局部热动平衡过程。

4.2 脉动周期与谱线行为

短周期（143-273天）米拉变星倾向于显示14.98μm发射特征，而长周期（≥308天）样本多表现为吸收态。具有硅酸盐发射SE7/SE8分类的恒星（如R Cet、Z Cyg）均伴随脉泽活动，且CO₂谱线光学厚度显著高于其他样本。

4.3 大气结构

CO₂在米拉大气中呈现高度延伸分布，吸收组中探测到距恒星2-4R_★的"制冷区"，温度低于辐射平衡预测值。该低温环境使硅酸盐尘粒可在更接近恒星处凝结，突破传统模型对尘埃形成位置的限制。发射组则揭示距星体5-15R_★的延展气体层，柱密度随径向距离增加而增大，与Cherchneff激波化学模型预测一致。

本研究通过多波段协同观测揭示米拉变星大气中CO₂的动态分布特征，首次为"制冷区"理论提供观测证据。双层RADEX模型证明恒星脉动驱动大气层周期性膨胀/收缩，导致分子激发条件改变谱线形态。13μm特征与CO₂谱线的关联性暗示尘埃-气体相互作用在星风驱动中的关键作用。该研究方法为解析詹姆斯·韦伯太空望远镜（JWST）中红外光谱数据建立新范式，对理解恒星演化末期的物质循环机制具有里程碑意义。