[論文レビュー] 3D chemical structure of diffuse turbulent ISM. I. Statistics of the HI-to-H$_2$ transition

Elena Bellomi, B. Godard+5 more|SPIRE - Sciences Po Institutional REpository|Sep 11, 2020

Astrophysics and Star Formation Studies参考文献 176被引用数 4

ひとこと要約

本研究では、詳細な化学反応を組み込んだ3次元磁気流体力学的シミュレーションを用いて、拡散乱した乱流状態の銀河間空間媒体（ISM）をモデル化し、標準的なISM設定（200 pcのボックス、n̄_H = 1–2 cm⁻³、G₀ = 1）が観測されたHIからH₂への遷移統計を驚くほどよく再現することを示している。キングフィッシャー図（HおよびH₂の2次元密度分布ヒストグラム）は、解像度、磁場（≤4 μG）、重力、乱流駆動力の変化に対して安定しており、H₂は約3–10 pcの冷たい中性媒体（CNM）構造で形成され、HIは約10–20 pcの密度相関構造で蓄積されることを示している。

ABSTRACT

We studied the statistical properties of the HI-to-H$_2$ transition observed in absorption in the local diffuse and multiphase ISM to identify the physical processes controlling the probability of occurrence of any line of sight. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chemistry of H$_2$. The impacts of the UV radiation field, the mean density, the turbulent forcing, the integral scale, the magnetic field, and the gravity on the molecular content of the gas are explored through a parametric study covering a wide range of physical conditions. The statistics of the HI-to-H$_2$ transition are interpreted through analytical prescriptions and compared with the observations using a modified and robust version of the Kolmogorov-Smirnov test. The results of one simulation, convolved with the distribution of distances of the observational sample, are able to explain most of the statistical properties of the HI-to-H$_2$ transition observed in the local ISM. The tightest agreement is obtained for a neutral diffuse gas modeled over ~200 pc, with a mean density of $1-2$ cm$^{-3}$, illuminated by the standard interstellar UV radiation field, and stirred up by a large-scale compressive turbulent forcing. Within this configuration, the 2D probability histogram of the column densities of H and H$_2$ is remarkably stable and is almost unaltered by gravity, the strength of the turbulent forcing, the resolution of the simulation, or the strength of the magnetic field $B_x$. The weak effect of the resolution and our analytical prescription suggest that the column densities of HI are likely built up in large-scale WNM and CNM structures correlated in density over ~20 pc and ~10 pc, respectively, while those of H$_2$ are built up in CNM structures between ~3 pc and ~10 pc.

研究の動機と目的

• 観測データと数値的データを統合して、局所的な拡散・多相状態のISMにおけるHIからH₂への遷移の統計的性質を理解すること。
• 異なる視線方向における出現確率や、視線方向にわたる分子分率の変動を制御する物理的メカニズムを特定すること。
• 乱流、UV放射、磁場、重力、密度の揺らぎが、原子状および分子状水素の分布に与える影響を評価すること。
• 修正されたコルモゴロフ＝スミルノフ検定を用いて、3次元数値モデルを観測された密度分布と照合すること。
• キングフィッシャー図を、拡散状態のISMの物理的状態を探査する安定な統計的ツールとして確立すること。

提案手法

• RAMSESコードを用いて、MHD、熱的進化、H₂化学を組み込んだ、乱流状態で磁場が作用する化学反応を伴うISMをシミュレートする。
• UV放射（G₀）、平均密度（n̄_H）、乱流駆動力（圧縮型対比非圧縮型）、ボックスサイズ（L）、磁場（Bₓ）、重力のパラメータを変化させたパラメトリックスタディを実施する。
• 観測されたHIおよびH₂の密度分布と比較するため、修正されたコルモゴロフ＝スミルノフ検定を適用する。
• 観測された背景源の距離分布と合成することで、シミュレーション結果を観測統計と一致させる。
• N(H)対N(H₂)の2次元確率ヒストグラム（キングフィッシャー図）を分析し、物理的および数値的変化に対する安定性を評価する。
• 解析的パラダイムを用いて、密度分布構造の起源とそのスケール依存性を解釈する。

実験結果

リサーチクエスチョン

• RQ1 どのような物理的条件が、局所的な拡散状態のISMにおける観測されたHIからH₂への遷移の統計的分布を再現するか？
• RQ2 乱流駆動力、磁場、重力、UV放射が、拡散状態のISMにおけるH₂の形成と分布にどのように影響を与えるか？
• RQ3 キングフィッシャー図（N(H) 対 N(H₂) の2次元ヒストグラム）は、シミュレーションの解像度、磁場強度、乱流駆動力の変化に対してどの程度安定しているか？
• RQ4 シミュレートされた拡散状態のISMにおける、HIおよびH₂の密度分布構造の特徴的な空間スケールは何か？
• RQ5 1つの3次元シミュレーションを、観測された距離分布と合成することで、観測されたHIからH₂への遷移の位置、幅、分散、形状を再現できるか？

主な発見

• 観測とシミュレーションの間で最も良好な一致が得られるのは、200 pcのボックス、平均密度 n̄_H = 1–2 cm⁻³、標準的なUV放射場（G₀ = 1）の設定である。
• キングフィッシャー図は、解像度、磁場（≤4 μG）、重力、乱流駆動力タイプの変化に対して顕著に安定しており、駆動力が圧縮型である限り成立する。
• HIの密度は、約20 pc（WNM）および約10 pc（CNM）のスケールで相関する大規模構造に蓄積されるが、H₂は約3–10 pcスケールの小さなCNM構造で形成される。
• キングフィッシャー図はドーピラー幅パラメータおよび磁場強度に対して弱く依存しており、これらのパラメータの診断的有用性は限定的である。
• シミュレートされた視線長さの分布は、約100 pcから約3 kpcの範囲で観測された平坦な対数分布と一致している。
• 統計的検定の結果、シミュレーションで3倍以上に過小または過大に予測される視線方向は全体の約0.5にとどまり、全体として強い一致が示されている。