A HINSA view of cosmic-ray ionization in IC 348 and NGC 1333: evidence for a strong low-energy cosmic-ray disparity

技術要約：IC 348およびNGC 1333におけるHINSAによる宇宙線電離の観測

問題提起
宇宙線電離率（CRIR, $\zeta_{ion}$）は、特に磁場とガスの結合を制御することによって、分子雲の化学的および力学的進化を支配する基本的なパラメータである。近年の観測により、大質量星形成領域や原始 stars（プロトスター）付近において上昇したCRIRの値が明らかになっているが、低エネルギー宇宙線（LECRs, $E < 1$ GeV）の物理的な起源および加速メカニズムについては依然として不明である。さらに、近傍のペルセウス分子雲複合体において不一致が存在する。NGC 1333は拡散的なサブGeVガンマ線残差の欠損を示している一方で、その隣にあるIC 348は正の残差とより硬い宇宙線エネルギースペクトルを示すが、NGC 1333の方が星形成活動が高い。局所的な宇宙線源がどの程度まで領域内の銀河背景放射を圧倒しているのか、またその格差の性質がどのようなものであるのかを解明する必要がある。

手法
本研究では、H I narrow self-absorption（HINSA）法を用いて、IC 348およびNGC 1333に対する新たなCRIRの推定値を提示する。手法は以下の通りである：

1. 観測： FAST望遠鏡（NGC 1333）からのH IデータおよびGALFA-H Iサーベイ（IC 348）を利用し、これらをCOMPLETEサーベイおよびJCMTからのアーカイブ$^{13}$CO ($J=1-0$ および $3-2$) データ、ならびにHerschel Gould Belt Surveyからのダスト連続波データと組み合わせる。
2. HINSA解析： $^{13}$COの速度場に対するH I吸収特徴の二次微分フィッティングを通じて、原子状水素カラム密度 $N(\text{HINSA})$ を導出する。HINSA分率（$f_{\text{HINSA}} = N(\text{HINSA})/N(\text{H}_2)$）は、光解離の影響を最小限に抑えるため、$N(\text{H}_2) > 3 \times 10^{21}$ cm$^{-2}$ の領域に対して計算される。
3. 密度制約： $^{13}$COの輝線強度によって制約された非LTE放射伝達モデリング（RADEX）およびマルコフ連鎖モンテカルロ（MCMC）サンプリングを用いて、ガス体積密度 $n(\text{H}_2)$ を推定する。
4. CRIRの算出： ダスト粒子上でのH$_2$生成と、宇宙線による破壊（電離、解離、および宇宙線誘起UVプロセスを含む）の定常状態のバランスから $\zeta_{ion}$ を導出する。
5. スペクトルモデリング： 有限スラブ減衰フレームワーク内でLECRスペクトルをモデリングする。このモデルは、観測された $\zeta_{ion}$ プロファイルを適合させるために、Fermiガンマ線観測から推測される高エネルギー宇宙線スペクトルの制約を取り入れつつ、低エネルギースペクトルの傾き（$\alpha$）と運動量ブレイク（$p_0$）を変化させる。

主な結果

• CRIRの格差： 本研究は、2つの雲の間で顕著なCRIRの差を見出した。IC 348における加重平均 $\zeta_{ion}$ は、$(5.35 \pm 1.49) \times 10^{-17}$ s$^{-1}$ から $(2.33 \pm 0.55) \times 10^{-16}$ s$^{-1}$ の範囲にある。対照的に、NGC 1333はそれよりも約1桁低く、$(3.94 \pm 2.46) \times 10^{-18}$ s$^{-1}$ から $(1.68 \pm 0.25) \times 10^{-17}$ s$^{-1}$ の範囲にある。
• カラム密度の傾向： 両方の雲において、$\zeta_{ion}$ はH$_2$カラム密度の増加とともに減少しており、これは内部加速源を欠いた減衰モデルと一致している。
• スペクトルモデリング： 最良適合の低エネルギー宇宙線スペクトルは、観測されたプロファイルを再現している。モデルは、$\sim$10倍の差が生じている要因が、低エネルギー（$\sim$10 MeV）陽子フラックスの対応する差によって駆動されていることを示している。IC 348のスペクトルはNGC 1333よりも著しく硬く、かつ強度が大きく、NGC 1333はGeVエネルギー領域において背景放射に対して欠損を示している。
• 密度および構造解析： 著者らは、ガス密度推定および雲の構造（3D-PDRシミュレーションによる）からのバイアスの可能性を調査した。$n(\text{H}_2)$ の不確実性やフラクタルな雲構造が絶対的な $\zeta_{ion}$ の値にバイアスを与える可能性があるが（低 $N(\text{H}_2)$ では過小評価、高 $N(\text{H}_2)$ では過大評価の可能性）、解析の結果、これらの要因が両領域間の1桁の格差を説明することはできないと結論付けた。

意義および主張
本論文は、観測されたCRIRの格差は、進化のタイムスケールや測定バイアスの産物ではなく、雲の物理的特性に固有のものであると主張している。

• 外部加速源： 高いCRIRを持つIC 348と、その硬いガンマ線スペクトルは、高密度雲コアの外側に存在する局所的なLECR加速源の存在を示唆しており、それはNGC 1333で支配的な原始星活動とは異なるものである可能性が高い。著者らは、より進化した、より埋没度の低いIC 348クラスター（より多くのB型星を含む）における恒星風終端衝撃が、妥当なソースであると提案している。
• 原始星の寄与： 結果は、原始星ジェットが小規模なスケールでCRIRを増強する可能性があるものの（NGC 1333のマッピングで見られるように）、これらがこれらの環境における広域的なCRIRを駆動するものではないことを示唆している。特に、NGC 1333の方が高い星形成活動を持っているにもかかわらず、広域的なCRIRが低いという事実がこれを裏付けている。
• 星形成への影響： IC 348における全域的に高いCRIRは、より高い電離度を意味し、ガスと磁場の結合を強める。著者らは、これがコア崩壊中の磁気ブレーキを強化し、IC 348の原始星円盤がNGC 1333のものよりも平均して小さくなるという観測結果を説明できる可能性があると述べている。
• 手法の妥当性： 本研究は、低質量星形成領域におけるCRIRのプローブとしてのHINSA手法の妥当性を検証しており、絶対値を精緻化するための将来的な3Dモデリングの必要性を強調すると同時に、相対的な格差の堅牢性を確認している。