[论文解读] Numerical simulations of stellar SiO maser variability. Investigation of the effect of shocks

E. M. L. Humphreys, M. D. Gray+4 more|Research Explorer (The University of Manchester)|Feb 22, 2002

Astrophysics and Star Formation Studies参考文献 45被引用 75

一句话总结

本研究将恒星流体动力学脉动模型与SiO脉泽辐射转移代码相结合，模拟了AGB星体中脉泽的变异性，重点关注激波引起的物理条件变化。尽管模型再现了环状脉泽结构和自行运动等关键特征，但由于冷却过程过于简化、缺少时变红外辐射场以及缺乏化学演化，导致模型表现欠佳，表明仅靠激波特性无法完全解释观测到的脉泽变异性。

ABSTRACT

A stellar hydrodynamic pulsation model has been combined with a SiO maser model in an attempt to calculate the temporal variability of SiO maser emission in the circumstellar envelope (CE) of a model AGB star. This study investigates whether the variations in local physical conditions brought about by shocks are the predominant contributing factor to SiO maser variability because, in this work, the radiative part of the pump is constant. We find that some aspects of the variability are not consistent with a pump provided by shock-enhanced collisions alone. In these simulations, gas parcels of relatively enhanced SiO abundance are distributed in a model CE by a Monte Carlo method, at a single epoch of the stellar cycle. From this epoch on, Lagrangian motions of individual parcels are calculated according to the velocity fields encountered in the model CE during the stellar pulsation cycle. The potentially masing gas parcels therefore experience different densities and temperatures, and have varying line-of-sight velocity gradients throughout the stellar cycle, which may or may not be suitable to produce maser emission. At each epoch (separated by 16.6 days), emission lines from the parcels are combined to produce synthetic spectra and VLBI-type images. We report here the results for v=1, J=1-0 (43-GHz) and J=2-1 (86-GHz) masers.

研究动机与目标

• 探究激波引起的密度、温度和速度梯度变化是否为AGB星体中SiO脉泽变异性主要驱动力。
• 检验在结合动态激波结构时，恒定辐射泵浦机制是否能再现观测到的脉泽光曲线与运动学特征。
• 评估非绝热激波、冷却 timescales 和化学不均一性在塑造脉泽发射轮廓与空间形态中的作用。
• 识别当前模型中的不足之处——特别是辐射转移、激波冷却和化学演化方面的缺陷——这些缺陷导致无法准确再现观测到的脉泽行为。
• 通过识别关键改进方向，为未来建模提供指导：时变红外场、准确的激波冷却模型，以及SiO和CO、H2O等冷却剂物种的动态化学过程。

提出的方法

• 采用拉格朗日方法追踪具有增强SiO丰度的气团，沿脉动周期运动，利用流体动力学恒星模型提供的速度场。
• 每隔16.6天使用辐射转移代码计算脉泽发射，该代码基于局部密度、温度和速度梯度计算脉泽增益。
• 生成$ v=1 $, $ J=1-0 $（43-GHz）和$ J=2-1 $（86-GHz）跃迁的合成光谱与VLBI类图像，以与观测结果对比。
• 模型假设恒定辐射泵浦，以隔离由激波增强的碰撞泵浦的作用。
• 激波冷却通过流体动力学timescale（60–100天）建模，但由于缺少CO等分子冷却剂，导致冷却时间过长。
• 采用蒙特卡洛方法在单一初始相位分布SiO富集气团，以模拟致密、不均一的包层环境。

实验结果

研究问题

• RQ1 仅靠激波引起的密度、温度和速度梯度变化，在多大程度上能解释观测到的SiO脉泽时间变异性？
• RQ2 尽管激波结构合理，为何模拟仍无法再现观测脉泽光曲线的完整动态范围？
• RQ3 缺少时变红外辐射场和准确的激波冷却如何影响脉泽放大与发射持续时间？
• RQ4 化学不均一性——特别是SiO及CO等冷却剂分子的生成与破坏——在脉泽变异性中扮演何种角色？
• RQ5 球对称的脉动模型能否再现观测中脉泽环在不同恒星相位下的持久性？还是非球对称结构才是必需的？

主要发现

• 强SiO脉泽发射在恒星光球面几倍恒星半径范围内形成环状空心结构，与TX Cam等恒星观测到的切向放大型态一致。
• 模型中合成脉泽环的直径随恒星脉动相位变化，与TX Cam及其他Mira变星的观测结果相符。
• 模型中明亮脉泽气团的自行速度与R Aqr和TX Cam中实测值相当，但并非所有观测到的运动都能被再现。
• 合成光谱显示线轮廓在恒星径向速度处达到峰值，典型Mira变星的线宽与强度变化与恒星相位相关。
• 模型低估了脉泽亮度，且无法再现完整的动态范围，主要原因是激波冷却timescale过长（60–100天）以及缺少CO等分子冷却剂。
• 模型预测脉泽发射偶尔会出现在恒星圆面之上，该特征已被某些恒星的VLBI观测所证实。