Ibn型超新星「SN 2023uqf」の可視光観測から、宇宙線をPeV級まで加速できる短い期間が浮かび上がった。理化学研究所などの研究チームは、光度曲線を再現した放射流体モデルから衝撃波の半径、速度、周囲のガス密度を取り出し、時間依存のニュートリノ生成計算へつないだ。2023年にIceCubeが捉えた442 TeVのニュートリノ事象「IC-231004A」は、モデルが予測する時期とエネルギーに収まる。ただし、同じモデルが見積もる検出数は1天体あたり最大でも約1万分の1にとどまる。今回確かめたのは発生源であることではなく、この超新星が候補から物理的に脱落しないことだ。成果は2026年7月8日付の『The Astrophysical Journal Letters』に掲載され、理研が同月21日に発表した。
可視光から衝撃波の履歴を逆算
SN 2023uqfは、IC-231004Aが2023年10月4日に到来した約21時間後、Zwicky Transient Facility（ZTF）の追跡観測で見つかった。地球からの距離は723 Mpc、約24億光年である。水素に乏しくヘリウムに富む星周物質へ超新星の放出物質が衝突するIbn型に分類され、明るさは絶対g等級で少なくともマイナス19.7に達した。半値以上の明るさを保った期間は8日と短く、ZTFが捉えたIbn型の中で最も速く変化した天体だった。
先行研究は、半解析的なMOSFiTモデルで光度曲線を合わせ、爆発時刻をZTF初検出の2.4日前と推定していた。これに対して今回の研究は1次元放射流体コード「STELLA」を使い、光の放射と衝撃波の運動を同時に解いた。狙いは光度曲線の再現そのものではない。可視光を生んだ衝撃波が、いつ、どれほど速く、どの密度のガスを通過したかを得ることにある。
基準モデルは、放出物質を2太陽質量、爆発の運動エネルギーを10^51 ergに固定した。星周物質はヘリウム90%で、密度が半径の3乗に反比例して急減すると仮定する。密度規格化をD'=50とした計算は、ZTFのg、r、iバンドの変化を初検出後約30日までおおむね再現した。そこから推定した爆発時刻は初検出の約10日前で、MOSFiTの値より早い。どちらも直接測定ではなく、採用した物理モデルに依存する。
この差はニュートリノとの時間関係を変える。新しい時間軸では、IC-231004Aが届いた時点で衝撃波はすでに濃い星周物質の中を進み、粒子加速が可能になる時期へ入っていたことになる。可視光を時刻合わせの手掛かりに使うのではなく、粒子加速を決める動力学へ変換した点が今回の前進である。
5.74日後に開くペバトロンの窓
爆発直後の衝撃波は光子に支配されるため、荷電粒子を往復させて加速する仕組みが働きにくい。基準モデルでは、衝撃波の前方にある物質の光学的厚さが下がり、無衝突衝撃波として粒子を加速できるのはコア崩壊から5.74日後だった。前方衝撃波はこの時期に毎秒約10^9 cmで進む。
研究チームは、加速される宇宙線をヘリウム原子核とし、衝撃波のエネルギーの10%が宇宙線へ、1%が磁場へ回る基準値を置いた。加速粒子が領域から逃げるまでの時間で上限を決めると、最大エネルギーは複数PeVの領域へ届く。IC-231004Aの442 TeVニュートリノを作るには、親となるハドロンが少なくとも数PeVのエネルギーを持たなければならない。基準モデルはこの条件を満たす。
濃いガスは別の役割も担う。加速された原子核が星周物質と衝突するとパイ中間子が生まれ、その崩壊でニュートリノが放出される。密度が高ければ衝突は増えるが、爆発直後は衝撃波が放射媒介となり加速を妨げる。時間がたちすぎればガスが薄くなり、衝突効率が落ちる。SN 2023uqfでは、無衝突衝撃波への移行と十分なガス密度が重なる短い期間が、過渡的なペバトロンの窓になる。
期待検出数は2.0×10^-5〜9.7×10^-5
物理的に作れることと、実際にIceCubeで1件捉えることの間には大きな隔たりがある。研究チームがモデルのニュートリノ放射にIceCubeの有効面積を畳み込むと、723 Mpc離れたSN 2023uqfから期待されるトラック状事象は、選別条件に応じて2.0×10^-5〜9.7×10^-5件だった。最も楽観的な条件でも、同条件の天体約1万個から1件に相当する。
この低計数域では、モデルの時間分布を「関連する確率」として読むことはできない。1件を検出したと条件づけたうえで、どの時刻に出やすいかを比べる重みである。IC-231004Aの検出時刻は重みの高い期間に入り、442 TeVというエネルギーも逃走制限で決めたスペクトルと矛盾しなかった。それでも、単一事象からSN 2023uqfを発生源と特定する統計的な根拠にはならない。
先行する観測研究は、43回のZTFニュートリノ追跡で同様のIbn型が偶然見つかる確率を0.28±0.06%と推定していた。これは「追跡領域に珍しいIbn型が現れる頻度」を測る値だ。一方、2.0×10^-5〜9.7×10^-5件という値は「基準モデルの1天体がIceCubeへ何件もたらすか」を表す。母集団や選択効果を組み込まないまま、二つの数字を同じ確率として比べることはできない。今回の研究も偶然一致の評価は更新していない。
近傍天体と損失込みの計算が判定を進める
基準モデルには、ペバトロン候補としての強さを左右する未解決点がある。最大エネルギーは粒子の逃走時間で制限しており、光中間子生成や光分解、非弾性原子核衝突などの損失を明示的に計算していない。損失をすべて入れれば、高エネルギー側の上限が下がる可能性がある。ヘリウム原子核の衝突も核子の重ね合わせで近似しており、著者らは数倍規模の系統誤差を見込む。
光度曲線の再現にも一意性はない。放出物質の質量と爆発エネルギーを固定し、星周物質を一つのべき乗則で表したため、約10日以降には観測とのずれが残る。SN 2023uqfは明るく急速に変化するIbn型の中でも、濃い星周物質を持つ部分集合の例とみるべきで、Ibn型全体を同じ効率のニュートリノ源とは扱えない。
観測面では距離が効く。同様の光度と衝撃波を持つ超新星なら、期待検出数はおおむね距離の2乗に反比例する。SN 2023uqfより近いIbn型を爆発直後から可視光、X線、電波で追い、IceCubeをはじめとするニュートリノ観測と重ねれば、星周物質の密度と粒子加速の開始時刻を厳しく絞れる。損失込みの計算が442 TeVまで届き、近傍例で複数事象が積み上がるか。それが超新星ペバトロンを候補から発生源へ進める条件になる。
