[論文レビュー] Radiative transfer with POLARIS: I. Analysis of magnetic fields through synthetic dust continuum polarization measurements

Stefan Reißl, S. Wolf+1 more|arXiv (Cornell University)|Apr 18, 2016

Astrophysics and Star Formation Studies参考文献 80被引用数 3

ひとこと要約

この論文は、星形成領域における磁場構造の分析を可能にする、複数の波長でダスト連続スペクトルの偏光をシミュレートすることを目的とした3次元モンテカルロ放射移動コードPOLARISを紹介する。放射移動と先進的な粒子配列理論を統合することで、POLARISは、磁気力学的（MHD）モデルに基づく円周星周囲ディスクや雲構造のような複雑な系においても、理論的予測やMC3Dなどのベンチマークコードと一致するダスト温度分布と合成偏光マップを正確に再現する。

ABSTRACT

Aims: We present POLARIS (POLArized RadIation Simulator), a newly developed three-dimensional Monte-Carlo radiative transfer code. POLARIS was designed to calculate dust temperature, polarization maps, and spectral energy distributions. It is optimized to handle data that results from sophisticated magneto-hydrodynamic simulations. The main purpose of the code is to prepare and analyze multi-wavelength continuum polarization measurements in the context of magnetic field studies in the interstellar medium. An exemplary application is the investigation of the role of magnetic fields in star formation processes. Methods: We combine currently discussed state-of-the-art grain alignment theories with existing dust heating and polarization algorithms. We test the POLARIS code on multiple scales in complex astrophysical systems that are associated with different stages of star formation. POLARIS uses the full spectrum of dust polarization mechanisms to trace the underlying magnetic field morphology. Results: Resulting temperature distributions are consistent with the density and position of radiation sources resulting from magneto-hydrodynamic (MHD) - collapse simulations. The calculated layers of aligned dust grains in the considered cirumstellar disk models are in excellent agreement with theoretical predictions. Finally, we compute unique patterns in synthetic multi-wavelength polarization maps that are dependent on applied dust-model and grain-alignment theory in analytical cloud models.

研究の動機と目的

• 星間空間における磁場解析を目的とした、複数波長観測におけるダスト偏光をシミュレートできる3次元モンテカルロ放射移動コードの開発。
• 内部配列、パラ磁性緩和、放射場配列といった最新の粒子配列理論を統合した1つの計算フレームワークへの統合。
• 特に円周星周囲ディスクや磁気力学的（MHD）崩壊モデルのような複雑な天体物理学的環境において、既存のシミュレーションおよびベンチマーク結果との照合。
• ダストモデルおよび粒子配列メカニズムに依存する独自のパターンを示す合成偏光マッピングを可能にし、観測解釈を支援。
• 星形成に関連する非球形で不均一かつ磁化された環境における放射移動シミュレーションの精度向上。

提案手法

• 円周星周囲ディスクやMHDでシミュレートされた雲構造を含む、複雑な天体物理学的構造をモデル化するため、適応的解像度を有する3次元直交座標系オクトリー格子を採用。
• 光子パッケージの追跡にモンテカルロ法を適用し、確率的放射輸送を用いて加熱、散乱、偏光を計算。
• IDG、RAT、Goldの複数の粒子配列メカニズムを、解析的および数値近似による内部配列関数（複素誤差関数を用いた効率的計算を含む）を用いて統合。
• 偏光マップの信号対ノイズ比を向上させ、計算効率を向上させるために、強制的1回目の散乱アルゴリズムとピールオフ技術を採用。
• 複数の粒子配列メカニズムの重ね合わせを近似するため、修正されたレイリー還元係数の公式（R = (R1+R2+R3+R1R2+R1R3+R2R3+3R1R2R3)/(1+2R1R2+2R1R3+2R2R3+2R1R2R3)）を用い、合成配列効果をモデル化。
• 同一のMHDシミュレーション入力を使ってMC3Dコードと照合し、異なるディスクおよび雲モデルにおけるダスト温度および偏光フラックスマップを比較。

実験結果

リサーチクエスチョン

• RQ1 POLARISは、既存のコードと比較して、MHDでシミュレートされた星形成環境におけるダスト温度分布をどの程度正確に再現できるか？
• RQ2 POLARISが生成する合成偏光マップは、円周星周囲ディスクや高密度雲核における磁場の姿かたちをどの程度正確に反映しているか？
• RQ3 RAT、IDG、Goldといった粒子配列理論の変更が、複数波長の合成観測における偏光パターンにどのような影響を及えるか？
• RQ4 グリッド構造（直交座標系オクトリー対球座標系）の違いが、高密度領域における温度および偏光予測の正確さに与える影響は何か？
• RQ5 結合レイリー還元係数の公式は、1回のシミュレーション内で複数の粒子配列メカニズムの重ね合わせを実用的かつ近似的にモデル化できるか？

主な発見

• POLARISは、MHD崩壊シミュレーションからの放射場および密度構造と整合するダスト温度分布を正確に再現しており、内側ディスク縁付近でのわずかな偏差はグリッド幾何の違いに起因するにとどまる。
• コードは、円周星周囲ディスクモデルにおける配列したダスト粒子層を正確にモデル化しており、さまざまなメカニズム下での粒子配列理論的予測と優れた一致を示している。
• POLARISが計算する複数波長の合成偏光マップは、選択されたダストモデルおよび粒子配列理論に応じて独自で明確なパターンを示し、競合する物理的モデルの観測的区別を可能にする。
• グリッド構造の違い（直交座標系オクトリー対球座標系）があるにもかかわらず、POLARISの偏光フラックス予測は、ベンチマークコードMC3Dと期待される範囲内で一致しており、最適化されたアルゴリズムのおかげで信号対ノイズ比がより高い。
• 複素誤差関数（erfi）を用いた内部配列の数値近似により、粒子配列度の効率的かつ正確な計算が可能となり、超熱的領域では1に収束することが確認された。
• 結合レイリー還元係数の公式は、複数の粒子配列メカニズムの重ね合わせを実用的かつ近似的にモデル化する有効な手法を提供するが、完全な正確性を保証するにはさらなる照合が必要である。