[論文レビュー] 3D Grid-Based Monte Carlo Code for Radiative Transfer through Raman and Rayleigh Scattering with Atomic Hydrogen -- STaRS

Seok-Jun Chang, Hee‐Won Lee|arXiv (Cornell University)|Dec 7, 2020

Spectroscopy and Laser Applications参考文献 21被引用数 7

ひとこと要約

本稿では、シンボリック星や活動銀河核などの天体物理学的環境において、原子水素内のラマン散乱およびレイリー散乱をモデル化するための3次元グリッドベースのモンテカルロ放射移動コードSTaRSを提示する。このコードは、線幅や偏光パターンを含む多次元分光偏光的特徴を正確にシミュレートでき、幾何学的切断がない完全な空間的カバレッジのおかげで、先行研究と定量的に一致し、解像度が向上している。

ABSTRACT

Emission features formed through Raman scattering with atomic hydrogen provide unique and crucial information to probe the distribution and kinematics of a thick neutral region illuminated by a strong far UV emission source. We introduce a new 3 dimensional Monte-Carlo code in order to describe the radiative transfer of line photons that are subject to Raman and Rayleigh scattering with atomic hydrogen. In this code entitled "${\bf S}$ejong Radiative ${\bf T}$r${\bf a}$nsfer through ${\bf R}$aman and Rayleigh ${\bf S}$cattering (${\it STaRS}$), each photon is traced until escape with a tag attached carrying information including the position, direction, wavelength, and polarization. The thick neutral scattering region is divided into numerous cells with each cell being characterized by its velocity and density, which ensures huge flexibility of the code in analyzing Raman-scattered features formed in a neutral region with complicated kinematics and density distribution. As a test of the code, we revisit the formation of Balmer wings through Raman scattering of far UV continuum near Ly$β$ and Ly$γ$ in a static neutral region. An additional check is made to investigate Raman scattering of O VI in an expanding neutral medium. We find fairly good agreement of our results with previous works, demonstrating the capability of dealing with radiative transfer modeling that can be applied to spectropolarimetric imaging observations of various objects including symbiotic stars, young planetary nebulae, and active galactic nuclei.

研究の動機と目的

• 厚い中性水素領域におけるラマン散乱およびレイリー散乱を通る放射移動をモデリングするための柔軟で3次元グリッドベースのモンテカルロコードの開発。
• シンボリック星や活動銀河核のような天体系において、線幅や偏光方向を含む分光偏光的特徴の正確なシミュレーションを可能にする。
• 特にLee & Lee (1997)のモデルと比較して、6825 Åにおけるラマン散乱O VI特徴を再現することで、コードのテストおよび検証を行う。
• 画像化応用のための表面輝度、偏光度、偏光方向マップの生成能力を示す。
• 天体物理学コミュニティが複雑で非等方的な散乱幾何学におけるラマン散乱を研究できるように、GitHub経由でオープンソースで提供するツールを提供する。

提案手法

• 空間領域を3次元デカルトグリッドで離散化することで、密度および速度構造が変化する散乱領域の高解像度モデリングを可能にする。
• 原子水素におけるレイリーおよびラマン過程の波長依存散乱断面積を用いたモンテカルロ光子トレースを採用する。
• 完全な散乱相関関数（特にトンプソンに類似した角度分布）を組み込むことで、非等方的幾何学における偏光発展を正確にモデリングする。
• 速度場を球対称な拡散流とし、最終速度 $v_{\infty}$ および密度プロファイル $n(\mathbf{r}) \propto (R_*/r)^2(1 - R_*/r)^{-1}$ を仮定する。
• ストークスパラメータ $I$, $Q$, および $U$ を計算して強度および偏光を追跡し、$P = Q/I$ で線形偏光の度合いと方向を示す。
• 吸収および再放射はエネルギー保存遷移により処理され、ラマン散乱によって遠紫外光子（例：O VI λ1032）が可視光Hα光子に変換される。

実験結果

リサーチクエスチョン

• RQ1 3次元グリッドベースのモンテカルロコードは、拡張するH I領域におけるラマン散乱O VI放射の線幅および偏光特徴をどの程度正確に再現できるか？
• RQ2 従来のモデルでは幾何学的切断が生じるが、その影響は線幅および偏光にどのような影響を及ぼし、完全な空間的カバレッジが結果をどのように改善するか？
• RQ3 コードは表面輝度、偏光度、偏光方向を含む多次元分光偏光画像化をどの程度正確にシミュレートできるか？
• RQ4 コードは、厚い中性媒体内における遠紫外源からのHα光子のラマン変換効率およびフラックス分布をどの程度正確にモデリングできるか？
• RQ5 コードは、シンボリック星における複数ピーク線幅や偏光反転といった既知の観測的特徴を信頼性高く再現できるか？

主な発見

• STaRSコードは、Lee & Lee (1997)の結果と良好に一致し、特に偏光挙動および線幅形状において、6825 Åにおけるラマン散乱O VI特徴を正確に再現した。
• 幾何学的切断がないため、STaRSでは線幅の赤ピークがLee & Lee (1997)と比較して10–15%強化された。
• 観測者方向へ脱出するラマン散乱光子の割合は、光学濃度 $\tau_0 = 0.5$, 1, 10 のそれぞれで0.11%、0.21%、1.82%であり、光学濃度に強く依存することが示された。
• コードは赤翼における偏光反転を正しく捉えており、青い光子は対称軸に対して垂直に偏光し、赤い光子は拡張領域での散乱により弱く平行に偏光する。
• 偏光度は観測者方向の視線（x–y平面）に沿って最大に達し、軸対称系における理論的期待と一致する。
• STaRSは、表面輝度、偏光度、偏光方向マップを含む高精度な分光偏光画像データを生成でき、グリッド解像度が向上するにつれて解析解に収束することが示された。