Constraining r-process nucleosynthesis with multi-objective Galactic chemical evolution models

技術的概要：多目的銀河化学進化モデルによる r 過程核合成の制約

問題提起
急速中性子捕獲過程（r 過程）の天体物理学的起源は依然として不確実であり、中性子星合体（NSM）と磁気回転超新星（MR-SNe）を含む競合するシナリオが存在する。マルチメッセンジャー観測は NSM を実在可能なサイトとして確認しているが、その遅延した発生と低い発生率は、最も初期の低金属量星における r 過程の増殖を説明する能力に疑問を投げかけている。一方、MR-SNe のような即時的な源は理論的に妥当であるが、その発生頻度と収量に関する直接的な観測的制約は欠けている。銀河化学進化（GCE）モデリングにおける重要な課題は、「r 過程普遍性」の仮定、すなわち太陽系 r 過程分布に基づき、r 過程元素の相対的豊度パターンがユーロピウム（Eu）から他の元素へ均一にスケーリングするという仮説である。低金属量星の観測は、この普遍性が重元素（$Z \ge 56$）に対しては成り立つ可能性があるが、r 過程に乏しい星において Eu に対してしばしば増強された比率を示す、より軽い中性子捕獲元素（$38 < Z < 50$）に対しては成り立たないことを示唆している。現在のモデルはしばしば単一目的最適化（Eu へのみ適合させること）に依存しており、複数の元素の化学進化を同時に再現するために必要なトレードオフを潜在的に隠蔽している可能性がある。

手法
著者らは、$\sim 1.5 \times 10^5$ の 1 回降着・均質銀河化学進化モデルのグリッドを用いた、体系的かつデータ駆動型のアプローチを採用している。r 過程を特定の天体物理サイト（例えば NSM や MR-SNe）に固定するのではなく、4 つの主要パラメータを変化させる柔軟なパラメトリック枠組みを採用している：

1. イベントあたりの Eu 収量（$Y_{Eu,r}$）： $[10^{-7}, 10^{-5}] M_\odot$ の範囲で対数的にサンプリング。
2. イベント率割合（$\alpha$）： r 過程物質を生成する親星の割合。$[0, 1]$ の範囲で線形にサンプリング。
3. 遅延時間（$\tau$）： 恒星の死から増殖までの時間。$[0, 14]$ Gyr の範囲で準対数的にサンプリング。
4. 親星の質量範囲（$[M_l, M_u]$）： $[8, 10)$ から $[95, 100] M_\odot$ まで $5 M_\odot$ 刻みで変化。

他の中性子捕獲元素の収量は、太陽系 r 過程パターン（Simmerer et al. 2004）に従って Eu 収量をスケーリングすることで導出される。モデルは、Gaia-ESO サーチ、MINCE、および SAGA データベースのデータを組み合わせて構築された、広範な金属量範囲（$-2.7 \le [\text{Fe/H}] \le \text{太陽}$）にわたる観測データセットと比較される。

モデルの性能を評価するため、著者らはパレートフロンティア分析に基づく多目的最適化戦略を採用している。単一の「最良適合」モデルを求めるのではなく、ある元素への適合を改善することが必然的に他の元素への適合を悪化させるような、パレート最適解の集合を特定する。これにより、異なる元素制約間のトレードオフ（例：Eu 対 Zr、または軽元素群 Sr-Y-Zr 対重元素群 Ba-La-Ce）のマッピングが可能となる。

主要な貢献と結果

• パラメータ制約： 最良のパフォーマンスを示すモデルは、一貫して短い遅延時間（$\tau \le 30$ Myr）、低質量の親星（中心は $20-25 M_\odot$ 付近）、およびイベントあたり $\sim 2 \times 10^{-7} M_\odot$ の実効的な Eu 注入率を支持する。$\sim 80 M_\odot$ よりも質量の大きい親星は、高い収量であっても、観測された r 過程増殖を単独で説明するにはあまりに希少であることが判明した。
• 軽元素に対する太陽系スケーリングの失敗： 最良の Eu モデルを太陽系 r 過程パターンにスケーリングすると、中程度から低金属量において、軽元素（Sr, Zr）の豊度を体系的に過小評価する一方、重元素（Ba, La, Ce）はより成功裡に再現する。
• パレートトレードオフ：
  • Eu 対 Zr： 2 次元最適化において、Eu と Zr の両方の傾向を同時に再現する単一のモデルは存在しない。Zr に良く適合するモデルは体系的に Eu を過大評価し、その逆もまた然りである。
  • 軽元素群対重元素群： 軽元素群（Sr, Y, Zr）に対する 3 次元最適化は、大きな矛盾を明らかにしており、パレートフロンティアは理想解から大きく離れている。対照的に、重元素群（Ba, La, Ce）は、Eu と同時に、はるかに低いトレードオフと理想点へのより小さな距離で再現可能である。
  • グローバル比較： Eu と軽元素群（Sr-Y-Zr）を組み合わせると、Eu と重元素群を組み合わせる場合（$D \approx 50$）と比較して、目的空間における理想点までの距離が著しく大きくなる（$D \approx 200$）。これは、太陽系スケーリング関係がより軽い中性子捕獲元素に向かうにつれて、徐々に妥当性を失うことを確認するものである。

意義と結論
本論文は、太陽系スケーリング収量の仮定の下で動作する単一の r 過程イベントのクラスでは、天の川銀河における軽および重中性子捕獲元素の両方の化学進化を同時に説明することはできないと結論付けている。観測された緊張関係は、r 過程普遍性が質量範囲全体にわたって拡張しないことを示唆している。

むしろ、結果は少なくとも 2 つの異なるスケーリング成分の必要性を示唆している：

1. 主要 r 過程成分（H 成分）： 太陽系パターンに類似し、r 過程に富む星で観測されるものであり、重元素（$Z \ge 56$）の生成を支配する。
2. 弱い r 過程成分（L 成分）： ランタノイドに対して軽元素（Sr, Y, Zr）の生成が強化されることを特徴とし、r 過程に乏しい星で観測されるパターンに類似する。

著者らは、これらの成分が異なる天体物理的条件またはイベントに由来する可能性があると仮定しているが、特定の物理的サイトを特定しているわけではない。本研究は、化学進化モデルにおいて普遍的な r 過程パターンを仮定することの限界を浮き彫りにし、複雑な核合成寄与を解きほぐすための多目的最適化の有用性を示している。今後の研究としては、GCE モデル内、あるいはシミュレートされた豊度パターンを通じて、弱い成分を直接制約することが提案されている。