[論文レビュー] Origin of the Galactic 511 keV emission from positrons produced in irregular supernovae

Hagai B. Perets|arXiv (Cornell University)|Jul 8, 2014

Gamma-ray bursts and supernovae参考文献 30被引用数 12

ひとこと要約

本論文は、ヘリウムの爆発が白色矮星で発生する弱いカルシウム豊富なType Ib超新星（2005Eに類似したSNe）が、$^{44}$Tiの崩壊によって十分な量の陽電子を生成し、銀河系の511 keVガンマ線放射を説明できると提唱する。これらの超新星は主に古く、星族IIの領域に存在するため、観測された銀河中心部に偏った分布と全放射能を自然に再現でき、ダークマターの存在や精密に調整された天体的源の必要性を排除する。

ABSTRACT

Gamma ray emission of 511 keV lines arising from electron-positron annihilation has been detected from the Galaxy since the 70s. Spatially resolved observations using the INTEGRAL satellite have shown its full sky distribution to be strongly concentrated in the Galactic bulge, with a smaller contribution from the disk, unlike the situation in any other wavelength. The puzzling distribution of the positrons gave rise to various suggestions, including stellar nucleosynthesis in core-collapse (CC) and type Ia thermonuclear supernovae (SNe), accreting compact objects, and more "exotic" explanations of annihilation of dark-matter particles. However, such models encounter difficulties in reproducing the total Galactic 511 keV flux as well as its peculiar bulge-centered distribution. Theoretical models of SNe from thermonuclear Helium detonations on white dwarfs (WDs) were suggested as potential additional sources of positrons, contributing to the Galactic Gamma-ray emission. Here, we show that the recently discovered class of faint, calcium-rich type Ib SNe (with the prototype SN 2005E), thought to arise from such explosions, also called ".Ia" SNe, can explain both the 511 keV flux and its distribution, and can eliminate the need for non-astrophysical (e.g., dark matter annihilation) sources. Such SNe comprise a fraction of only ~2% of all SNe, but they currently inject hundreds of times more positrons (from $^{44}$Ti decay) to the ISM than injected by CC SNe, enough to reproduce the total Galactic 511 keV line emission. They exclusively explode in old environments, and their contribution to the 511 keV emission therefore follows the old (>10 Gyrs) Galactic stellar population, dominated by the bulge, thereby reproducing the observed large bulge-to-disk ratio.

研究の動機と目的

• 銀河系の511 keVガンマ線放射の起源を説明すること。これは数十年にわたる観測にもかかわらず、未だに解明されていない。
• 既知の源から予想される円盤支配の陽電子消失と対照的に、陽電子消失が銀河中心部に強く偏っているという謎を解明すること。
• かつて見過ごされてきた弱くカルシウムが豊富なType Ib超新星というクラスが、放射の全放射能と空間的分布の両方を説明できることを示すこと。
• ダークマターの消失など非天体的説明を排除し、観測的制約と整合的で、精密調整の必要のない天体的源を提示すること。
• 過剰な調整や特異な物理学を必要とせず、観測された陽電子生成率と銀河中心部対円盤比を再現できる物理的に妥当なメカニズムを提供すること。

提案手法

• $^{44}$Ti崩壊による陽電子生成を2005Eに類似した超新星に対してモデル化し、観測された噴出物質量とカルシウムが豊富な組成を用いて$^{44}$Ti生成量を推定する。
• 星族合成法を採用し、楕円銀河や銀河のハローでの観測環境に基づき、2005Eに類似した超新星を古く（100億年以上）の星族に割り当てる。
• Picardら（2003年）に従い、銀河の構造に「銀河中心部＋厚い円盤」モデルを採用し、空間的に分解された陽電子注入率を計算する。
• INTEGRAL観測データから得られた値（$\dot{N}_{e^+} = 19.6_{-2}^{+2.3} \times 10^{42} \, \text{s}^{-1}$、$B/D \approx 1.42$）と、モデル予測の陽電子生成率を比較する。
• 超新星発生率の推定誤差（検出バイアスによる低減）と星族質量関数の不確実性を組み込む。
• コアコリプス型および標準的Type Ia超新星からの$^{44}$Ti寄与を既知の値と比較し、その合計寄与が支配的でないことを確認する。

実験結果

リサーチクエスチョン

• RQ1 2005Eに類似した超新星は、銀河系における観測された511 keV放射全流量を説明するのに十分な量の陽電子を生成できるか？
• RQ2 2005Eに類似した超新星の空間的分布（古く、星族IIに存在）は、511 keV放射の観測された銀河中心部対円盤比を再現できるか？
• RQ3 このメカニズムは、過剰な調整や特異な物理学を必要とせず、観測された陽電子生成率と（数MeVの）注入エネルギーを再現できるか？
• RQ4 2005Eに類似した超新星が511 keV放射の厚い円盤成分にどの程度寄与するか？また、薄い円盤成分と比べてどうなるか？
• RQ5 現在の望遠鏡（例：NuSTAR）で検出可能な$^{44}$Ti線（68/78 keV）が、銀河中心部やM31の超新星残骸で観測可能であり、これにより本モデルが確認できる予測は可能か？

主な発見

• 2005Eに類似した超新星モデルは、全陽電子生成率が$12.7 \pm 3.7 \times 10^{42} \, \text{s}^{-1}$と予測され、観測値$19.6_{-2}^{+2.3} \times 10^{42} \, \text{s}^{-1}$と整合的である。
• モデルは、観測された銀河中心部対円盤陽電子生成比$B/D = 1.5 \pm 0.6$を再現しており、観測値$1.42_{-0.30}^{+0.34}$と一致する。
• 厚い円盤への寄与は$1.95 \pm 0.73 \times 10^{42} \, \text{s}^{-1}$であり、厚い円盤対薄い円盤比が0.3に相当する。これは、中心ブラックホールやダークマター源とは整合しない。
• モデルは、全Type Ia超新星の約10％が2005Eに類似していると予測しており、これは銀河中心部での検出バイアスにより、実際の発生率が低く見積もられている可能性がある。
• 今後の観測で、銀河中心部やM31の超新星残骸（例：SN 1885A）において$^{44}$Ti崩壊の68/78 keV線が検出可能であり、これにより本モデルの妥当性を確認できる重要なテストが可能となる。
• コアコリプス型超新星からの$^{44}$Ti寄与（最大$3 \times 10^{42} \, \text{s}^{-1}$）を含めても、モデルの整合性は損なわれず、全銀河系の$\dot{N}_{e^+}^{\text{MW}} \approx 19 \times 10^{42} \, \text{s}^{-1}$および$B/D \approx 1.1$が得られ、観測誤差範囲内に収まる。