[論文レビュー] New mass estimates for massive binary systems: a probabilistic approach using polarimetric radiative transfer

Andrew G. Fullard, Jack O’Brien+5 more|arXiv (Cornell University)|Feb 21, 2022

Gamma-ray bursts and supernovae参考文献 67被引用数 5

ひとこと要約

本論文は、偏光的放射移動を用いた、質量の大きなWolf-Rayet + O連星星の質量を推定する新しい確率的手法を提示する。この手法では、モンテカルロ放射移動モデルをニューラルネットワークでエミュレートすることで、解析的モデルに起因するバイアスを克服する。このアプローチにより、より正確な軌道傾き角の推定が可能となり、特にWR 153について顕著な既存推定値との不一致が生じ、質量制約が再評価された。これは、質量の大きな連星の進化モデルを精緻化するため、偏光およびスペクトル観測データの改善が不可欠であることを示している。

ABSTRACT

Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths.

研究の動機と目的

• 質量の大きな連星系、特にWolf-Rayet + O星の質量推定を改善すること。従来の手法はバイアスと不確実性に苦しんでいる。
• 解析的偏光モデルの限界に対処すること。これらは高軌道傾き角にバイアスがかかる上、光学的に薄い散乱領域に限定される。
• 位相偏光と放射移動を用いて、軌道傾き角と質量を推定する計算効率が高く、正確な手法を開発すること。
• 軌道傾き角、光学的厚さ、フラックス比といったモデルパラメータの信頼性の高い事後分布を提供し、不確実性を明確にすること。
• 公開済みの偏光データを用いて4つのWR + O連星系を再分析することで、本手法が従来の解析的手法を上回ることを実証すること。

提案手法

• 光度が厚い風を有する連星系からの偏光光をシミュレートするために、モンテカルロ放射移動（MCRT）コードであるslipが使用される。
• MCRTシミュレーションは、フル放射移動モデルを近似するニューラルネットワークエミュレータを訓練することで高速化され、迅速なサンプリングが可能になる。
• エミュレータを用いてベイズ推論が実施され、軌道傾き角、光学的厚さ、フラックス比といったモデルパラメータの事後分布が探索される。
• システムパラメータに関する事前知識が組み込まれ、マルコフ連鎖モンテカルロ（MCMC）サンプリングを用いて信頼区間が推定される。
• トレーニング範囲内での正確性を保証するため、エミュレータはフルMCRTシミュレーションと照合されて検証される。
• 光学的厚さ（τ）、風の非対称性（zell）、フラックス比（Ifrac）といったモデルパラメータが同時に制約され、強いデゲネラシーが明らかになった。

実験結果

リサーチクエスチョン

• RQ1 ニューラルネットワークでエミュレートされた放射移動モデルは、質量の大きな連星系において、従来の解析的モデルよりも正確でバイアスのない軌道傾き角推定を可能にするか？
• RQ2 偏光データからの再評価された軌道傾き角推定値は、Wolf-Rayet + O連星系の導出される質量にどのように影響するか？
• RQ3 光学的厚さ、風の非対称性、フラックス比といった主要パラメータの事後分布は何か？また、それらの相関関係は？
• RQ4 新しい質量推定値は、既存の推定値とどの程度異なるか、特にWR 153においては？そしてこれは進化モデルにどのような含意をもたらすか？
• RQ5 この手法は、離心率を持つ軌道やBe星を伴う連星系を含む、他の質量の大きな連星系へ一般化可能か？

主な発見

• 本手法は、位相偏光とエミュレートされた放射移動を用いることで、M sin iのデゲネラシーを効果的に解消し、軌道傾き角を高精度に推定できた。
• WR 153において、新しい質量推定値は、以前の偏光推定値と顕著に異なる結果となり、従来のモデルがバイアスを有していたことを示唆している。
• 全4系統のシステムにおいて、昇交点線の経度（Ω）の事後分布は、不確実性の適切な評価がなかった過去の推定値と完全に食い違う結果となった。
• 光学的厚さ（τ）は通常低く（約0.2）、これは、バンド幅偏光がWolf-Rayet風の外側、密度が低い領域で発生していることを示している。
• τとzellの間、およびIfracとτの間に強い相関関係が確認され、解釈の際にこれらデゲネラシーを考慮しなければならないことが明らかになった。
• 本手法は、すべてのシステムにおいて質量の信頼性の高い上限を提供し、不確実性を明確に定量化しており、解析的モデルに対する統計的に妥当な代替手法を提供している。