[論文レビュー] A genetic algorithm approach to fitting interferometric data of post-AGB objects: I. the case of the Ant nebula

Daniel Macdonald, Orsola De Marco+3 more|arXiv (Cornell University)|Apr 29, 2017

Stellar, planetary, and galactic studies参考文献 1被引用数 5

ひとこと要約

本論文は、光干渉計測定データに放射線輸送モデルと幾何学的モデルをフィットさせることで、後-AGB惑星状星雲環境をモデル化する遺伝的アルゴリズムベースのPythonモジュールGADRADを紹介する。これは不偏な多次元パラメータ確率密度関数を生成し、相関関係やデゲネラシーを明らかにする。アンツ星雲（Mz3）に適用したところ、$3.5^{+7.5}_{-2.2} \times 10^{-5}$ M$_\odot$ のほこりを含む惑星環が得られ、先行研究と整合的であるが、信頼区間が改善されており、Mz3の磁気的共通エンvelope噴出起源を示唆する。

ABSTRACT

We present GADRAD, a Python module that adopts heuristic search techniques in the form of genetic algorithms, to efficiently model post-asymptotic giant branch (post- AGB) disc environments. GADRAD systematically constructs the multi-dimensional pa- parameter probability density functions that arise from the fitting of radiative transfer and geometric models to optical interferometric data products. The result provides unbiased descriptions of the object's potential morphology, component luminosities and temperatures, dust composition, disc density profiles and mass. Correlation in the estimated parameters as well as potential degeneracies are revealed. Estimated probability distributions of the post-AGB environment parameters provide insight into the shaping processes that may occur in the transition from the post-AGB to the planetary nebula phase. We test parameter recovery on simulated artificial data products of a typical post-AGB environment. We then use GADRAD to model the mid-infrared spectrum and visibilities of the Ant nebula (Mz3), taken with the Very Large Telescope Interferometer's instrument MIDI. Our result is consistent with a large dusty disc with similar parameter values to those previously found by Chesneau et al., except for a larger dust mass of $3.5^{+7.5}_{-2.2} imes10^{-5}$ M$_{\odot}$. The parameter confidence intervals determined by GADRAD, can however be relied upon to impose additional constraints on all disc and system parameters. Based on our analysis and other considerations, we tentatively suggest that Mz3 is a pre-PN ejected during a magnetic (polar) common envelope interaction, where the binary may or may not have survived at the core of the nebula.

研究の動機と目的

• 複雑な惑星環形状を示す後-AGB天体の干渉計測定データに適した、頑健でヒューリスティックに基づくフィッティング手法の開発。
• パラメータの不確実性と相関関係を考慮した、多次元パラメータ確率密度関数を体系的に構築すること。
• 干渉計測定のVisibilityおよびスペクトルデータを用いて、後-AGB系におけるほこり質量、温度、組成、密度プロファイルの制約を強化すること。
• 実観測（例：アンツ星雲）への適用前に、シミュレートされたデータを用いてパラメータ回復の信頼性をテストすること。
• パラメータ制約とデゲネラシー解析を通じて、特に惑星状星雲への移行を引き起こす物理的メカニズムを探索すること。

提案手法

• GADRADは、高次元パラメータ空間におけるヒューリスティックなグローバル最適化を実行する遺伝的アルゴリズムを用い、複雑な尤度関数の表面を効率的に探索する。
• 放射線輸送モデルと幾何学的惑星環モデルを組み合わせ、非常に大型望遠鏡干渉計（VLT-I）のMIDI機器からの中赤外可視度およびスペクトルをシミュレートする。
• 遺伝的アルゴリズムの枠組み内において、マルコフ連鎖モンテカルロ法を用いてモデルパラメータの事後確率分布を推定する。
• パラメータ空間の探索には、ほこりの組成、半径方向密度プロファイル、惑星環のスケールヘイト、全放射光度、温度、全ほこり質量を含む。
• 多次元確率密度推定を用いて、パラメータのデゲネラシーと相関関係を体系的に同定する。
• 実データ（例：Mz3）への適用の前に、パラメータ回復性能の妥当性を確認するために、シミュレートされた干渉計測定データを用いる。

実験結果

リサーチクエスチョン

• RQ1 遺伝的アルゴリズムに基づくフィッティング手法は、干渉計測定データから後-AGB系の主要な惑星環パラメータをどの程度正確に回復できるか？
• RQ2 中赤外可視度およびスペクトルのフィッティングにおいて、惑星環の形状、ほこり質量、温度に関する主要なパラメータのデゲネラシーと相関関係は何か？
• RQ3 フィットされたモデルは、アンツ星雲（Mz3）の進化歴にどのような物理的制約を課えるか？
• RQ4 推定された惑星環パラメータは、Mz3の形成に磁気的共通エンベロープ噴出シナリオを支持するか？
• RQ5 フィットされたパラメータの信頼区間は、先行研究と比較して、システムの物理的性質の制約をどのように改善するか？

主な発見

• GADRADは、シミュレートされた干渉計測定データにおいて入力パラメータを正確に回復でき、実世界への応用における頑健さと信頼性を検証した。
• アンツ星雲に対して、$3.5^{+7.5}_{-2.2} \times 10^{-5}$ M$_\odot$ のほこり質量が得られ、先行研究と整合的であるが、不確実性の評価が改善された。
• この手法は、ほこり質量、温度、密度プロファイルパラメータの間で顕著なパラメータの相関関係とデゲネラシーを明らかにした。
• 推定された惑星環パラメータは、数百年AUにわたる広大でほこりの多い惑星環で、フラレート（膨らんだ）幾何形状を示している。
• GADRADによる信頼区間は、先行研究よりも、惑星環およびシステムのパラメータに対する制約をきめ細かく与えることで、物理的解釈可能性を向上させた。
• モデル結果および追加的考察に基づき、本研究は、Mz3が磁気的（極性）共通エンベロープ相互作用によって形成された可能性を一時的に提案する。この場合、二重星は星雲の中心に残存している可能性がある。