Nuclear Isomers and Their Impact on Gamma-Ray Emission in Binary Neutron Star Mergers

技術要約：核異性体と連星中性子星合体におけるガンマ線放出への影響

問題提起
連星中性子星（BNS）合体からのガンマ線（$\gamma$線）放出に関する現在のモデルは、特定の原子核に対して不正確である可能性のある2つの暗黙的な仮定に依存している。(1) 核反応ネットワークが基底状態の種のみを追跡していること、および (2) 異性体遷移が基底状態への瞬時な脱励起として扱われ、放出されるすべての$\gamma$線が親核の$\beta$崩壊のタイムスケールに割り当てられていることである。これらの仮定は、長寿命の核異性体（メタステーブル状態）が関与する場合に成立しない。このような異性体は、崩壊のタイムスケールを変化させ、異なる崩壊チャネルを開き（基底状態をバイパスする）、温度依存的な遷移率を示す可能性がある。その結果、標準的なパイプラインは、$\gamma$線放出のタイミング、強度、およびスペクトル特性を誤って特徴付け、rプロセス核合成の直接的なプローブを覆い隠してしまう可能性がある。

手法
著者らは、核異性体をrプロセス反応ネットワーク内の独立した種として動的にモデル化するフレームワークを開発した。この手法は、主に以下の4つの段階で進行する。

1. 候補の選定: 原子番号のチャートから、著者らは、rプロセスでアクセス可能であり、半減期が分単位から1年であり、$\gamma$線エネルギーが$>0.3$ MeV（将来のMeV観測装置の感度範囲内）であり、かつ異性体遷移（IT）による分岐比が$>99\%$であるという、事前フィルタリングを適用した。これにより、初期の50個の候補から12個へと絞り込まれた。さらに、30のポスト合体噴出物軌道（一般相対論的磁気流体力学シミュレーションによる）における核合成収量に基づく事後的な選定により、有意な存在量を持つ6つの異性体にまで絞り込まれた：$^{91\text{m}}\text{Y}$, $^{97\text{m}}\text{Nb}$, $^{90\text{m}}\text{Zr}$, $^{175\text{m}}\text{Yb}$, $^{179\text{m}}\text{Lu}$, および $^{207\text{m}}\text{Tl}$ である。

2. 遷移率の理論的枠組み: 著者らは、基底状態と異性体状態間の温度依存的な有効遷移率（$\Lambda_{AB}$）を計算した。熱的な光子バス（プランク分布）と準位の縮退を考慮し、中間核準位を経由した遷移確率を解くために、2つの確立された定式化（M1およびM2）を利用した。

   • M2 (厳密解): 中間準位を介したランダムウォークによる最終状態への到達確率を決定するために、線形方程式系を解く。
   • M1 (近似解): 部分和行列アプローチを用い、M2が数値的に不安定になる可能性のある低温域において堅牢である。
   • $\beta$フィード: 高温では娘核の中間励起状態が占有されやすくなり、それらが基底状態よりも異性体へ優先的にカスケードすること（$\beta$フィード）を認め、温度依存的な$\beta$フィード確率を組み込んでいる。
   • データの補完: 実験的な核データ（スピン、パリティ、強度）が不完全な場合、著者らは遷移率の計算に理論的なワイスコップ推定値を用いた。

3. ネットワークの実装: PRISM（Portable Routines for Integrated nucleosynthesis Modeling）反応ネットワークを拡張し、異性体を独立した種（異なる$(Z, N)$ペアとしてエンコード）として扱うようにした。ネットワークは、$T \approx 10$ GKから$\sim 60$日間にわたって進化する、BNS合体シミュレーションからの30の代表的な噴出物軌道に対して実行された。

4. 観測可能性の分析: 選択された異性体について$\gamma$線フラックスを計算し、それを2つの背景放射、すなわち噴出物の熱的ブラックボディ連続成分および他の不安定なrプロセス核からの$\beta$崩壊連続成分と比較した。光子の脱出を、一様なスラブ（悲観的）と均質な球体（オスターブロック・モデル）という2つの限定的な幾何学的形状を用いてモデル化し、放射輸送の仮定に対する感度の幅を決定した。フラックスは、15 kpcの典型的な距離において、現在および提案されているMeV観測装置（COSI, AMEGO, LOX）の感度に対して評価された。

主な貢献

• 動的な異性体モデリング: 基底状態の静的な近似を超え、温度依存的な有効遷移率および$\beta$フィード確率の堅牢な実装を提供している。
• ネットワークの拡張: 異性体を独立した種として明示的に追跡するようにPRISMネットワークを拡張したことを示し、$^{26\text{m}}\text{Al}$のようなよく研究されている事例によって検証されている。
• スペクトル特性: $\beta$崩壊の広範な連続成分や熱的背景とは異なる、異性体遷移によって生成される特定の鋭い$\gamma$線ラインを特定した。

結果

• 存在量: 6つの異性体が、シミュレートされた噴出物中で有意な存在量（$Y \ge 10^{-12}$）に達した。
• 支配的なライン: $^{97\text{m}}\text{Nb}$ (743.3 keV) と $^{91\text{m}}\text{Y}$ (555.6 keV) の遷移が、それぞれのエネルギーにおける$\gamma$線フラックスの主要な寄与者として特定された。
• 検出可能性:
  • $^{97\text{m}}\text{Nb}$ の743.3 keV ラインと $^{91\text{m}}\text{Y}$ の555.6 keV ラインは、銀河系内イベント（15 kpc）において、スラブおよび球体両方の脱出幾何学において、COSI、AMEGO、およびLOXによって検出可能であると予測されている。
  • $^{175\text{m}}\text{Yb}$ と $^{179\text{m}}\text{Lu}$ は、AMEGOおよびLOXの球体幾何学において限界的な検出可能性を示したが、スラブ幾何学では閾値を下回った。
  • その他の候補（$^{90\text{m}}\text{Zr}$, $^{207\text{m}}\text{Tl}$）は検出閾値を下回った。
• タイミング: 「熱的出現時間」（ラインフラックスが熱的背景を超える時間）は急速に（$\sim 0.3 - 0.7$ 秒）起こるが、主要なラインのピークフラックスは数時間から数日のタイムスケール（例：$^{91\text{m}}\text{Y}$については約20時間、$^{97\text{m}}\text{Nb}$については約1日）で発生し、光学的に厚い状態から薄い状態への転移と一致する。
• 幾何学依存性: ピークフラックスとタイミングは、仮定された脱出幾何学に強く依存しており、球体モデルは、光学的に厚い領域においてスラブモデルよりも大幅に高いフラックスを予測している。

意義と主張
著者らは、標準的なモデリング・パイプラインは、rプロセス核合成からの$\gamma$線放出を正確に特徴付けるには不十分であると主張している。彼らの主張は以下の通りである：

1. 異性体は観測可能である: 核異性体は、特定の核合成された核の直接的な指紋として機能する、明確で鋭い$\gamma$線ラインを生成することができ、キロノバの光度曲線に対する補完的なプローブを提供する。
2. 動的な取り扱いが不可欠である: これらの信号を正確に予測するには、反応ネットワークへの核異性体の堅牢な組み込みと、結果として生じる分光特性（特に温度依存的な遷移率および$\beta$フィード）の注意深いモデリングが必要である。
3. 限界: 著者らは、自らの結果が、厳密な最終計算を提供するものではなく、異性体の寄与の重要性の可能性を示すものであることを明示している。彼らは、決定的な予測を行うには、代表的なサブセットだけでなく噴出物全体をモデル化する必要があること、および詳細な放射輸送計算が必要であることを述べている。また、$\beta$フィード分岐比や準位構造に関して、不完全な核データが不確実性を導入することも認めている。

要約すると、本論文は、BNS合体における核異性体の検出のための理論的および計算的な基礎を確立しており、$^{97\text{m}}\text{Nb}$ と $^{91\text{m}}\text{Y}$ を、将来のMeVガンマ線天文学ミッションにおける主要な候補として強調している。