論文要約：原始星包囲雲におけるイオン化率の増大の証拠

論文タイトル: Evidence of Enhanced Ionization in Protostellar Envelopes
著者: K. R. Schwarz ら
掲載日: 2026 年 4 月 6 日（ドラフト版）

1. 研究の背景と問題意識

星形成過程におけるイオン化レベルは、ガスと磁場の結合（MHD 効果）や化学進化の両方に決定的な影響を及ぼします。特に、原始惑星系円盤（Class II 段階）では、ガス相の CO などが塵に吸着したり、化学反応によって減少したりする現象が観測されています。化学モデルでは、これらの変化を説明するには、宇宙線によるイオン化率（$\zeta$）が $10^{-17} \sim 10^{-18} \text{ s}^{-1}$ 程度であることが必要とされています。

しかし、より若い Class 0 段階（原始星がまだ包囲雲に埋没している段階）におけるイオン化環境は十分に解明されていません。過去の研究では、原始星の衝撃波や表面での宇宙線加速により、局所的なイオン化率が大幅に上昇する可能性が示唆されていましたが、Class 0 段階の包囲雲全体でのイオン化率を直接制約した研究は限られていました。本研究は、**「Class 0 段階の原始星包囲雲において、イオン化率が拡散星間物質（ISM）やより進化した円盤と比較して著しく増大しているか」**を検証することを目的としています。

2. 研究方法と手法

2.1 観測対象とデータ

研究対象は、3 つの Class 0 原始星（NGC 1333 IRAS4A, L1448-C, L1157）です。以下の観測データを統合して使用しました。

• NOEMA (Northern Extended Millimeter Array): 2016-2018 年に実施された新しい観測。$H^{13}CO^+ \ J=1-0$ 線（86.75 GHz）をターゲットとした。
• IRAM 30m 望遠鏡: 2022 年 7 月に実施された単一電波望遠鏡観測。同じく $H^{13}CO^+ \ 1-0$ 線を観測し、干渉計データでは失われるゼロ間隔（large-scale）のフラックスを補完した。
• CALYPSO サイバーデータ: 既存のアーカイブデータから、$C^{18}O \ J=2-1$ 線を使用。

2.2 データ処理と解析手法

1. データの統合: 単一電波望遠鏡（30m）と干渉計（NOEMA）のデータを結合し、空間分解能とフラックス保存の両方を確保したデータキューブを作成。
2. 吸収補正: IRAS4A の $C^{18}O$ スペクトルには、落下（infall）を示す逆 P シグニ型プロファイルが見られたため、ガウスプロファイルを用いて吸収成分を補正し、正確な柱密度を導出した。
3. 柱密度の算出:
   • $H^{13}CO^+$ と $C^{18}O$ の輝線強度から、局所熱平衡（LTE）を仮定して各分子の柱密度（$N$）を算出。
   • 励起温度には、各源の既存モデル（S. Anderl et al. 2016）から得られた値（26-28 K）を使用。
   • 同位体比（$^{12}C/^{13}C = 69$, $^{16}O/^{18}O = 557$）を仮定して、全 CO と全 $HCO^+$ の柱密度に変換。
4. イオン化率（$\zeta$）の導出:
   • 定常状態の化学反応モデルに基づき、$HCO^+$ と CO の存在量比からイオン化率を算出する式を再導出した（E. Redaelli et al. 2024 の手法を適用）。
   • 電子密度 $n(e)$ については、$n(e) \geq n(HCO^+)$ という保守的な仮定を置いた。
   • 水素分子密度 $n(H_2)$ は、既存の塵連続放射モデル（L. E. Kristensen et al. 2012）に基づく密度プロファイルを使用。

3. 主要な結果

3.1 イオン化率の増大

3 つのすべての源において、導出されたイオン化率は、拡散星間物質の平均値（$\approx 6 \times 10^{-17} \text{ s}^{-1}$）や、宇宙船観測に基づく太陽系周辺の値（$\approx 1.3 \times 10^{-17} \text{ s}^{-1}$）を数桁上回ることが判明しました。

• 導出されたイオン化率の範囲: $\zeta = 10^{-16} \sim 10^{-13} \text{ s}^{-1}$
• 源ごとの傾向:
  • IRAS4A, L1448-C: $\zeta \approx 10^{-16} \sim 10^{-15} \text{ s}^{-1}$
  • L1157: 他の源よりも顕著に高く、$\zeta \approx 10^{-14} \sim 10^{-13} \text{ s}^{-1}$
  • L1157 の高い値は、$H^{13}CO^+$ と $C^{18}O$ のフラックス比が他の源に比べて高いため導出されました。

3.2 空間分布

イオン化率は、$C^{18}O$ のフラックスが最大となる中心付近でピークを示す傾向がありますが、これは $C^{18}O$ が光学的に厚くなることによる過大評価の可能性も指摘されています。しかし、中心部から離れた領域（包囲雲の外側）においても、依然として増大したイオン化率（$10^{-16} \sim 10^{-13} \text{ s}^{-1}$）が観測されており、これは局所的な現象ではなく、包囲雲全体にわたる現象であることを示唆しています。

3.3 既存研究との比較

L1157 のアウトフロー内の弓状衝撃波（L1157-B1）における過去の研究（L. Podio et al. 2014）では、$\zeta \approx 3 \times 10^{-16} \text{ s}^{-1}$ と報告されていました。本研究で導出した包囲雲全体のイオン化率（特に L1157）はそれよりも高い値を示しましたが、これは解析手法の違い（化学モデル vs 解析的アプローチ）や、衝撃波前のガスの化学状態の違いによるものと考えられます。

4. 考察と意義

4.1 イオン化の源

Class 0 段階で観測される高いイオン化率の源として、以下の可能性が挙げられます。

• 宇宙線: 原始星のジェットや衝撃波、あるいは星風による加速により、宇宙線フラックスが増大している可能性。
• X 線: 原始星からの X 線放射（ただし、包囲雲による吸収が強く、深部への到達は限定的）。
• 紫外線 (UV): 衝撃波や原始星からの UV 放射（ただし、包囲雲は UV に対して不透明であり、深部への影響は小さい）。
  特に、宇宙線は包囲雲の深部や円盤内部まで浸透できるため、化学反応を駆動する主要な要因である可能性が高いです。

4.2 星形成化学への示唆

本研究の結果は、**「原始星進化の初期段階（Class 0 段階、年齢 $< 10^5$ 年）において、イオン化駆動型の化学反応が極めて効率的に進行している」**という仮説を強く支持します。

• Class II 段階（円盤が形成された後）で観測される揮発性物質の減少や化学的再処理は、Class 0 段階の初期に高いイオン化率によって引き起こされた化学変化の蓄積の結果である可能性があります。
• 従来の化学モデルが想定していた低いイオン化率では、Class II 段階で観測されるような化学的変化を説明できないため、初期段階でのイオン化率の見直しが必要です。

5. 結論

NOEMA と IRAM 30m 望遠鏡を用いた高感度観測により、3 つの Class 0 原始星の包囲雲において、イオン化率が拡散星間物質の値を数桁上回る（$10^{-16} \sim 10^{-13} \text{ s}^{-1}$）ことを初めて体系的に示しました。この「増大したイオン化」は、原始星進化の初期段階における化学進化の主要な駆動力であり、後の円盤段階での化学組成を決定づける重要なプロセスであると考えられます。今後のより詳細な化学モデルや、より多くの源を対象とした調査により、イオン化の具体的なメカニズム（宇宙線加速の有無など）の解明が期待されます。