Multi-messengers from the radioactive decay of $r$-process nuclei

本論文は、$r$過程原子核の$\beta$崩壊に由来する電子、ガンマ線、ニュートリノ、および中性子の時間依存する非熱的放射スペクトルを第一原理計算により提示し、キロノバ観測と並んで重元素形成の直接的かつ相補的なプローブとしてこれらのマルチメッセンジャー信号を確立する。

要約

宇宙を、金、白金、ウランなどの最重元素が調理されている巨大で混沌としたキッチンだと想像してみてください。この調理プロセスはr 過程と呼ばれ、2 つの中性子星の衝突のような極限の宇宙現象で起こります。

長年、科学者たちは、これらの現象が放つ光（「キロノバ」）を観測することで、この「調理」がどのように行われるかを正確に理解しようと試んできました。しかし、光を見ることは、完成したケーキだけを眺めてレシピを理解しようとするようなものです。個々の材料やオーブンの熱は見えません。

本論文は、オーブンの扉を開け、放射性の材料そのものから放出される熱と蒸気を直接観察することについて述べています。

以下に、著者たちが行ったことと発見したことを簡潔にまとめます。

1. レシピ：放射性崩壊を「粒子シェーカー」として

重元素が生成されると、それらは不安定になります。安定化するためには、余分なエネルギーを「振り落とす」必要があります。不安定な原子核を、強く振りすぎた炭酸飲料のボトルだと考えてください。蓋を開けると、中身が飛び散ります。

• 飛び散るもの： 炭酸飲料の代わりに、これらの原子は 4 種類の粒子を噴射します。電子（微小な帯電粒子）、ニュートリノ（ほとんど何とも相互作用しない幽霊のような粒子）、ガンマ線（高エネルギーの光）、そして中性子です。
• 目的： 著者たちは、時間経過に伴い、何が、どれくらい、そしてどの速さで放出されるかを正確に計算することを意図しました。

2. 手法：デジタルシミュレーション

実際の宇宙爆発（稀で遠く離れた現象）を待つ代わりに、科学者たちは超高精度のコンピュータシミュレーションを構築しました。

• 彼らは「核反応ネットワーク」を使用しました。これは、数百万種類もの異なる原子の材料を追跡する巨大なスプレッドシートのようなものです。
• これを詳細な物理モデルと組み合わせることで、各原子がどのように崩壊するかを正確に予測しました。
• 結果： 彼らは放出物の「メニュー」を作成し、最初の 1 秒から 1 年後までの、電子、ニュートリノ、ガンマ線、中性子のエネルギーと粒子数を示しました。

3. 大きな驚き：穏やかな暖房ではない

著者たちは、これらの爆発から放出されるエネルギーは、科学者が以前に想定していたものとは非常に異なることを発見しました。

• 「熱的」ではない： 通常、熱を考えると、均一で滑らかな分布（暖かいオーブンのようなもの）を想像します。しかし、著者たちはここにはそれが当てはまらないことを発見しました。粒子は「非熱的」であり、巨大で混沌としたエネルギーの奔流を伴って噴射されます。
  • 比喩： 焚き火を想像してください。「熱的」な火は、一定で穏やかな暖かい光を放ちます。一方、これらの核爆発は、高速で飛び散る巨大な火花に続き、長い間小さな火花が尾を引くような花火のショーに似ています。
• 「幽霊」粒子の勝利： 時間の大部分において、ニュートリノ（幽霊粒子）が全エネルギーの約 40% から 50% を運び去ります。残りを電子とガンマ線が分け合います。
• ガンマ線の「指紋」：
  • 初期： 原子の寿命が短く、変化が速すぎて特定のパターンが見えないため、ガンマ線はごちゃ混ぜのぼやけとして現れます。
  • 後期（数日〜数週間後）： 塵が落ち着くと、特定の「線」が現れます。これらはバーコードのようです。著者たちは、特定の原子（例えばタリウム -208）が明確な痕跡（2.6 MeV の線）を残すことを発見しました。これらの線を見ることができれば、どの重元素が生成されたかを正確に知ることができます。

4. 見えるか？（「聴く」部分）

この論文は問いかけます。「実際にこれらの粒子を検出できるのか？」と。

• 電子と中性子： いいえ。周囲の破片にすぐに捕捉されてしまいます。厚い霧の中で懐中電灯を見ようとするようなものです。
• ニュートリノ： はい、しかし困難です。幽霊のような性質のため、容易に逃げ出します。著者たちは計算により、もし我々の銀河（約 15,000 光年先）で巨大な爆発が起きた場合、ハイパー・カミオカンデ（巨大な水タンク）のような巨大な検出器が約2 個のニュートリノ事象を検出する可能性があるとしました。これは微小な信号ですが、確かに存在します。
• ガンマ線： はい、そしてここがエキサイティングな部分です。当初、破片が厚すぎてガンマ線は逃げ出せません。しかし、数日〜数週間後には霧が晴れます。著者たちは、将来のガンマ線望遠鏡で銀河を観測すれば、これらの特定の「バーコード」線を数週間、あるいは数ヶ月にわたって観測できる可能性があると提案しています。

結論

この論文は、重元素の生成から放出されるエネルギーに関する、新しく極めて詳細な「地図」を提供します。

• 重要性： 現在の宇宙爆発のモデルは、エネルギーの分布を推測することが多々あります。この論文は、その推測を精密な計算に置き換えます。
• 成果： これらの粒子がどのように放出されるかを正確に理解することで、天文学者はこれらの現象からの光をより良く解釈できるようになります。それ以上に、爆発の輝きに基づいて推測するのではなく、宇宙が最も重い元素をどのように生成するかを証明するために、核の「煙」（ニュートリノとガンマ線）を直接観測する扉を開くことになります。

技術概要

技術的サマリー：r 過程原子核の放射性崩壊からのマルチメッセンジャー

問題提起
急速中性子捕獲過程（r 過程）は、中性子星合体や潜在的にはガンマ線バーストのココンのような極限的な天体物理環境において発生し、重元素の合成を担っている。その結果生じる「キロノバ」遷移現象は、これらの r 過程原子核の放射性崩壊によって駆動されるが、現在のキロノバ光曲線のモデル化は、崩壊生成物からのエネルギー沈着に関する簡略化された仮定に依存している。具体的には、既存のモデルはしばしば崩壊 Q 値を通じてエネルギー生成率を計算し、放出される粒子（電子、ニュートリノ、$\gamma$線、中性子）に対して単純な熱化分布を仮定している。このアプローチは、熱化およびその結果生じる光曲線の進化を理解する上で決定的に重要な、崩壊生成物の詳細な非熱的スペクトル特徴を無視している。さらに、中性子過剰同位体に関する実験データの限界により、r 過程収量に対する異なる天体物理サイトからの相対的寄与は依然として不確かである。

手法
著者らは、詳細な原子核構造データと原子核反応ネットワークを結合させた、r 過程原子核の$\beta$崩壊に起因する放出スペクトルの第一原理計算を提示する。

• スペクトル計算: 本研究は、電子、反ニュートリノ、$\gamma$線、および$\beta$遅延中性子に対する放出スペクトルを計算するために、統計的・多段階フレームワーク（M. Mumpower ら 2025b に基づく）を利用する。$\beta$強度関数については、利用可能な場合は実験データを使用し、そうでない場合は準粒子ランダム位相近似（QRPA）を採用する。電子とニュートリノのスペクトルは、総崩壊エネルギーがこれら間で分配されると仮定して導出され、一方、$\gamma$線および中性子の放出は、励起状態の緩和および分岐比をモデル化するハウザー・フェッシュバック統計的アプローチを用いて計算される。
• 反応ネットワーク: これらの個別のスペクトルは、PRISM（Portable Routines for Integrated nucleoSynthesis Modeling）反応ネットワーク（v1.6.0）と結合される。このネットワークは、原子核入力に有限範囲液滴モデルを、放射捕獲および核分裂率に CoH3 コードを用いて、核合成をシミュレートする。
• 天体物理的軌跡: シミュレーションは、ガンマ線バーストのココンを表す特定の軌跡（M. R. Mumpower ら 2025a の軌跡 b）の進化を追跡する。この軌跡は、特にアクチノイド領域において太陽の r 過程残留物を再現する能力が選定理由である。物理的条件は、初期密度$3.2 \times 10^4$ g/cm$^3$、初期温度 2 GK、電子分率$Y_e = 0.034$の断熱膨張に従う。
• フラックス計算: 総放出スペクトル$S^{(i)}(E, t)$は、反応流$F(t)$（単位質量あたりの秒あたりの崩壊数）と崩壊あたりのスペクトル放出$S(E)$に因数分解される。これにより、時間関数としての総粒子数、エネルギーフラックス、および平均エネルギーの計算が可能となる。

主要な結果
本研究は、放出分布が時間依存性が強く、非熱的であり、実質的な平均エネルギーを有することを明らかにしている。

• 非熱的性質: 放出粒子の平均エネルギーは、熱平衡から期待されるもの（keV 範囲）よりも著しく高い。最初の数秒間、電子とニュートリノの平均エネルギーは約 5 MeV であり、10 秒後には 1.5〜2.5 MeV に低下し、中性子捕獲が停止するまで（約 100 秒）1 MeV 以上を維持する。対照的に、$\gamma$線と中性子はより低い平均エネルギー（$\lesssim$1 MeV）を維持するが、それでも非熱的である。
• スペクトル進化:
  • 電子とニュートリノ: これらのスペクトルは形状がほぼ同一であり、初期には高エネルギー（約 10〜20 MeV）まで伸びる。滑らかな高エネルギーの尾部は、核分裂生成物に起因し、長期的な時間スケールまで持続する。
  • $\gamma$線: 初期には低エネルギーピーク（約 1 MeV）と高エネルギー尾部を持つ二峰性であるが、約 1 秒後には単峰性分布へと進化する。初期段階では、短寿命原子核の理論的モデリングによりスペクトルは滑らかである。約 1 日後には、より長寿命の原子核（例えば$^{208}$Tl からの 2.6 MeV 線や$^{140}$La からの線など）に由来する特定のスペクトル線が現れる。
  • 中性子: 中性子スペクトルは初期には顕著であるが、$\beta$遅延中性子放出の減少と中性子捕獲の終了に伴い、約 100 秒後に急速に減少する。
• エネルギー分配: 崩壊エネルギーの配分は一定ではない。ニュートリノは進化を通じて最大の割合（40〜50%）を担う。初期段階（$\lesssim$1 秒）では、$\gamma$線の寄与は小さく、電子の寄与が支配的である。中性子のエネルギー寄与は、ごく初期の段階でのみ有意である。

意義と主張
著者らは、これらの計算された放出スペクトルが、将来のキロノバの放射輸送モデルにとって必要不可欠な自己整合的な入力データを提供すると主張する。簡略化された仮定を詳細な非熱的スペクトルに置き換えることで、これらの結果はキロノバ光曲線のより精密な解釈を可能にし、モデル化における縮退を軽減することを目指している。

さらに、本論文は銀河系内でのこれらの信号の直接観測の可能性を推定している。

• ニュートリノ: 銀河系内の事象（放出質量$0.01 M_\odot$、距離 15 kpc）の場合、最初の 10 秒間の放出は、逆$\beta$崩壊を介してハイパー・カミオカンデサイズの検出器で約 2.2 個の検出可能な事象をもたらす可能性がある。
• $\gamma$線: 初期には不透明であるが、放出物質の不透明度が低下するにつれ、$\gamma$線は事象発生から数週間から数ヶ月後に観測可能になると予想される。特定のスペクトル線（例えば$^{208}$Tl からの線など）は直接観測可能であり、キロノバ光曲線のモデル化に依存せずに核合成を探る手段を提供する。

本研究は、正確な観測可能性の時間スケールを決定するには完全な放射輸送が必要であるものの、銀河系内の r 過程事象からのニュートリノと$\gamma$線の直接検出は、モデル化の縮退を回避し、重元素形成を直接探る有望な手段であると結論付けている。