[论文解读] Mid-infrared interferometric monitoring of evolved stars - The dust shell around the Mira variable RR Aql at 13 epochs

I. Karovicova, M. Wittkowski+4 more|arXiv (Cornell University)|Jul 5, 2011

Stellar, planetary, and galactic studies被引用 13

一句话总结

本研究利用VLTI/MIDI仪器对Mira变星RR Aql进行了13个周期的中红外干涉观测，采用基于脉动的动态大气模型与辐射转移方法对尘埃壳层进行建模。数据最佳拟合为硅酸盐尘埃壳层，光学深度τ_V = 2.8 ± 0.8，内半径R_in = 4.1 ± 0.7 R_Photosphere，幂律指数p = 2.6 ± 0.3，表明脉动对中红外通量和可见度的影响分别低于25%和20%，因此在测量不确定度范围内无法检测到。

ABSTRACT

We obtained 13 epochs of mid-infrared interferometry with the MIDI instrument at the VLTI between April 2004 and July 2007, covering pulsation phases 0.45-0.85 within four cycles. The data are modeled with a radiative transfer model of the dust shell where the central stellar intensity profile is described by a series of dust-free dynamic model atmospheres based on self-excited pulsation models. We examined two dust species, silicate and Al2O3 grains. We performed model simulations using variations in model phase and dust shell parameters to investigate the expected variability of our photometric and interferometric data. The observed visibility spectra do not show any indication of variations as a function of pulsation phase and cycle. The observed photometry spectra may indicate intracycle and cycle-to-cycle variations at the level of 1-2 standard deviations. The best-fitting model for our average pulsation phase of 0.64+/-0.15 includes the dynamic model atmosphere M21n (T_model=2550 K) with a photospheric angular diameter of 7.6+/-0.6 mas, and a silicate dust shell with an optical depth of 2.8+/-0.8, an inner radius of 4.1+/-0.7 R_Phot, and a power-law index of the density distribution of 2.6+/-0.3. The addition of an Al2O3 dust shell did not improve the model fit. The photospheric angular diameter corresponds to a radius of 520^+230_-140 R_sun and an effective temperature of ~ 2420+/-200 K. Our modeling simulations confirm that significant visibility variations are not expected for RR Aql at mid-infrared wavelengths within our uncertainties. We conclude that our RR Aql data can be described by a pulsating atmosphere surrounded by a silicate dust shell. The effects of the pulsation on the mid-infrared flux and visibility values are expected to be less than about 25% and 20%, respectively, and are too low to be detected within our measurement uncertainties.

研究动机与目标

• 利用高分辨率中红外干涉测量，研究富氧Mira变星中脉动、尘埃凝结与质量损失之间的关联。
• 确定中红外可见度与通量变化是否能追踪RR Aql的脉动相位与周期。
• 利用动态自激脉动模型与辐射转移方法，对尘埃壳层结构与中心恒星大气进行建模。
• 评估Al2O3尘埃相对于硅酸盐在尘埃壳层中的贡献及其可检测性。
• 利用多 epoch 干涉数据，约束尘埃壳层与光球层的几何与物理参数。

提出的方法

• 在2004至2007年间，利用甚大望远镜干涉仪（VLTI）的MIDI仪器，对4个脉动周期内的RR Aql进行了多 epoch 中红外干涉观测。
• 采用基于自激脉动模拟的动态无尘模型大气（例如，T_model = 2550 K 的M21n模型）对中心恒星的强度分布进行建模。
• 应用尘埃壳层的辐射转移模型，假设密度分布为幂律形式（ρ ∝ r^−p），并变化光学深度（τ_V）与内半径（R_in）。
• 在模型中测试了两种尘埃成分——硅酸盐与Al2O3，以评估其对观测SED与可见度谱的相对贡献。
• 通过蒙特卡洛模拟传播不确定性，比较模拟与观测的通量与可见度谱在8–13 μm范围内的拟合优度，评估模型拟合效果。
• 通过模型模拟脉动相位与尘埃壳层参数的变化，预测预期的可见度与光度变异性水平。

实验结果

研究问题

• RQ1 RR Aql的中红外可见度与通量变化是否与脉动相位或周期相关？
• RQ2 能否通过围绕脉动动态大气的尘埃壳层辐射转移模型再现观测到的中红外SED与可见度曲线？
• RQ3 数据中是否能检测到Al2O3尘埃壳层的存在，还是与纯硅酸盐壳层一致？
• RQ4 哪些物理参数（光学深度、内半径、幂律指数）最能拟合多 epoch 干涉数据？
• RQ5 由脉动驱动的恒星大气变化在短时标与长时标上对中红外通量与可见度的影响程度如何？

主要发现

• 未检测到显著的周期内或周期间可见度变化，模型预测的变异性低于5–20%的测量不确定度。
• 光谱通量显示可能存在1–2σ水平的周期内与周期间变化，提示存在通量调制的微弱证据，但尚未得到确认。
• 最佳拟合模型包含一个硅酸盐尘埃壳层，光学深度τ_V = 2.8 ± 0.8，内半径R_in = 4.1 ± 0.7 R_Photosphere，幂律指数p = 2.6 ± 0.3。
• 中心动态模型大气M21n（T_model = 2550 K，Φ_model = 0.1）与光球角直径θ_Photosphere = 7.6 ± 0.6 mas 最佳拟合数据。
• 光球半径为R_phot = 520⁺²³⁰₋₁₄₀ R☉，有效温度T_eff ≈ 2420 ± 200 K，与模型大气一致。
• 尽管Al2O3尘埃未改善拟合效果，但模型模拟表明，低光学深度的内层Al2O3壳层无法被排除，可能对尘埃凝结序列有贡献。