[论文解读] The 2.35 year itch of Cygnus OB2#9 III. X-ray and radio emission analysis based on three dimensional hydrodynamical modelling

E. R. Parkin, J. M. Pittard+2 more|arXiv (Cornell University)|Jun 22, 2014

High-pressure geophysics and materials参考文献 2被引用 11

一句话总结

本研究利用三维自适应网格细化（AMR）流体动力学模拟，对天鹅座OB2 #9双星系统的风-风碰撞激波进行建模，以研究其X射线与非热射电辐射。结果表明，为匹配X射线观测，需采用较低的电子-离子温度比（~0.1）和降低的质量损失率（~6.5–7.5×10⁻⁷ M☉ yr⁻¹）；而射电辐射主要由同步辐射主导，最佳拟合模型显示磁场能量密度与相对论性能量密度之间存在非平衡状态，暗示其粒子加速物理机制超出标准模型的复杂性。

ABSTRACT

X-ray and radio data recently acquired as part of a project to study Cyg OB2#9 are used to constrain physical models of the binary system, providing in-depth knowledge about the wind-wind collision and the thermal, and non-thermal, emission arising from the shocks. We use a three-dimensional, adaptive mesh refinement simulation (including wind acceleration, radiative cooling, and the orbital motion of the stars) to model the gas dynamics of the wind-wind collision. The simulation output is used as the basis for radiative transfer calculations considering the thermal X-ray emission and the thermal/non-thermal radio emission. To obtain good agreement with the X-ray observations, our initial mass-loss rate estimates require a down-shift by a factor of roughly 7.7 to $6.5 imes10^{-7}$ and $7.5 imes10^{-7}$ solar mass per year for the primary and secondary star, respectively. Furthermore, the low gas densities and high shock velocities in Cyg OB2#9 are suggestive of unequal electron and ion temperatures, and the X-ray analysis indicates that an (immediately post-shock) electron-ion temperature ratio of $\simeq 0.1$ is also required. The radio emission is dominated by (non-thermal) synchrotron emission. A parameter space exploration provides evidence against models assuming equipartition between magnetic and relativistic energy densities. However, fits of comparable quality can be attained with models having stark contrasts in the ratio of magnetic-to-relativistic energy densities. The radio models also reveal a subtle effect whereby inverse Compton cooling leads to an increase in emissivity as a result of the synchrotron characteristic frequency being significantly reduced. Finally, using the results of the radio analysis, we estimate the surface magnetic field strengths to be $\approx 0.3-52\;$G. (Abridged)

研究动机与目标

• 利用多波段观测约束大质量双星天鹅座OB2 #9风-风碰撞区的物理条件。
• 在激波加热与粒子加速的背景下，探究X射线与非热射电辐射的起源。
• 检验标准假设（如磁场与相对论性能量密度之间的能量均分）在该系统中是否成立。
• 确定辐射冷却、风加速及视线方向对观测光曲线形状的影响。
• 基于射电辐射建模估算表面磁场强度，并与磁化恒星的观测值进行比较。

提出的方法

• 采用包含风加速、辐射冷却及双星轨道运动的三维自适应网格细化（AMR）流体动力学模拟。
• 将模拟输出作为输入，用于辐射转移计算，以推导出射的X射线与射电辐射谱。
• 利用流体动力学模拟中获得的激波后气体温度与密度结构，建模热X射线辐射。
• 通过同步辐射模拟非热射电辐射，引入相对论电子能量分布与磁场强度，参数化注入效率（ζ_rel, ζ_B）。
• 在ζ_rel、ζ_B与谱指数p的参数空间中进行搜索，以找到能匹配观测X射线与射电光曲线及谱的模型。
• 在射电辐射建模中考虑自由-自由吸收、Razin效应及逆康普顿冷却等消光机制。

实验结果

研究问题

• RQ1 为重现天鹅座OB2 #9的观测X射线谱，激波后区域所需的电子-离子温度比是多少？
• RQ2 推断的主星与伴星质量损失率与以往估计相比如何？为匹配X射线观测需作何调整？
• RQ3 观测到的非热射电辐射是否与激波区磁场与相对论性能量密度之间的能量均分一致？
• RQ4 预激波风速度与双星间距的变化如何调制整个轨道相位中X射线与射电光曲线？
• RQ5 何种物理机制可解释观测到的谱指数p ≪ 2？其与标准扩散激波加速理论预测相比如何？

主要发现

• X射线数据要求激波后电子-离子温度比约为0.1，表明电子加热效率极低。
• 所需的质量损失率相比以往估计显著降低，约低7.7倍，主星为6.5×10⁻⁷ M☉ yr⁻¹，伴星为7.5×10⁻⁷ M☉ yr⁻¹。
• 射电辐射主要由非热同步辐射主导，热辐射贡献可忽略不计，与早期射电观测结果一致。
• 即使允许谱指数p在1至2之间变化，假设磁场与相对论性能量密度均分的模型仍无法重现观测结果。
• 最佳拟合模型要求能量密度比存在显著差异：要么ζ_rel ≫ ζ_B（如ζ_rel = 0.15，ζ_B = 5×10⁻⁵，p = 1.6），要么ζ_rel ≪ ζ_B（如ζ_rel = 3×10⁻⁴，ζ_B = 0.5，p = 1.2），表明存在非标准的加速物理机制。
• 估算的表面磁场强度约为0.3–52 G，低于强磁化单星的典型值，且观测到的p ≪ 2值与上游粒子对激波特性的修正不一致，提示可能存在替代加速机制，如在多个弱激波处的再加速。