Time-dependent cosmic-ray escape from wind bubbles: hard spectra formation

技術要約：風によるバブルからの時間依存的な宇宙線脱出：硬いスペクトルの形成

問題提起
若い大質量星団（YMSC）、ウォルフ・レート星、あるいは超エディントン降着コンパクト天体などの源から生成される風駆動型のバブルは、宇宙線をPeVエネルギーまで加速する候補（PeVatron）として有望である。内側風終端衝撃波（termination shock）における拡散衝撃波加速（DSA）が主要な加速メカニズムとして広く受け入れられているが、これらの粒子が終端衝撃波から前方衝撃波（ISMへの脱出境界）へと輸送される時間依存的なプロセスについては、依然として十分に解明されていない。従来の研究の多くは、定常的な記述に依存するか、あるいは輸送スケープタイムを推定するにとどまっており、フル・タイム依存の輸送方程式を解くまでには至っていなかった。本研究は、バブルの動的な膨張と、終端衝撃波と前方衝撃波の間の距離の時間的変化が、脱出する宇宙線のエネルギー・スペクトルにどのような影響を与えるかという理解における空白を埋めるものである。

手法
著者らは、風駆動型バブルを、終端衝撃波（$R_{ts} \propto t^{2/5}$）と前方衝撃波（$R_{fs} \propto t^{3/5}$、減速期において）の両方が時間依存の速度で膨張する球対称系としてモデル化している。粒子の輸送は、衝撃波後方の風の流れによる移流と拡散を含む、1次元球形時間依存輸送方程式（フォッカー・プランク方程式）によって記述される。

この方程式を解くために、著者らは伊藤の補題（Itô's lemma）を用いて、方程式を確率微分方程式（SDE）の集合へと変換した。これらのSDEは、オープンソースのCRPropa3.2フレームワークの修正版を用いて積分されている。主な手法の特徴は以下の通りである：

• 注入： 粒子は、事前加速された冪乗則スペクトル（$dN/dE \propto E^{-2}$）と、ヒラス限界および加速時間によって決定される時間依存の高エネルギー・カットオフを持つ状態で、終端衝撃波（TS）に注入される。
• 移流： 下流領域では、放射状速度プロファイル $u(r) \propto r^{-2}$ が仮定され、衝撃波後方の速度に正規化される。
• 拡散： 本研究では、ボーム（$\delta=1$）、クライコネン（$\delta=3/2$）、コルモゴロフ（$\delta=5/3$）の3つの乱流モデルを調査している。拡散係数は、局所的な磁場強度と乱流パワースペクトルに依存する。
• 磁場： 磁場はラム圧の一定割合に基づいて推定される。研究では、下流磁場の半径方向の依存性（$B \propto r^{-\alpha}$、ここで $\alpha = 0, -1, -2$）および磁化レベル（$\epsilon_B$）をテストしている。
• シミュレーション・パラメータ： シミュレーションは、超新星爆発がエネルギー収支を支配する前の風駆動フェーズをカバーする、5 Myr（500万年）まで実行される。

主な結果
時間依存的な取り扱いにより、標準的な期待値からの顕著な逸脱が明らかになった：

1. スペクトルの硬化： 特定の物理的構成において、脱出する宇宙線スペクトルは、注入された $E^{-2}$ スペクトルよりも硬くなることがある。この硬化は、より強いエネルギー依存性を持つ乱流モデル（例：ボーム拡散）や、より弱い磁場依存性において最も顕著になる。
2. 低エネルギー抑制： 極めて重要な知見は、脱出フラックスにおける低エネルギー粒子の著しい抑制である。前方衝撃波（脱出境界）が時間の経過とともに源から遠ざかるため、拡散係数が小さい低エネルギー粒子は、脱出境界に追いつくことができず、下流領域に閉じ込められる。
   • ボーム拡散の場合、低エネルギー粒子（$E < 1$ TeV）は脱出フラックスからほぼ完全に消失する。
   • クライコネン拡散の場合、中程度の硬化（$\Delta \gamma \approx 0.1-0.2$）と部分的な低エネルギー抑制が観察される。
   • コルモゴロフ拡散の場合、減少は最小限であり、脱出スペクトルは注入スペクトルに酷似しているが、最大エネルギーは依然として低い（$\sim 20$ TeV）。
3. パラメータ依存性： 低エネルギーの減少度合いとスペクトルの硬化は、乱流モデル、磁場の半径方向の依存性（依存性が弱いほどスペクトルが硬くなる）、および磁化レベルに依存する。
4. 観測的含意： 地球への伝播をシミュレートした場合（3 kpcに位置する源を想定）、低エネルギー抑制は持続するが、長期間の銀河系内拡散によって硬化は部分的に平滑化される。

意義と主張
本論文は、風によるバブルの膨張という時間依存的な性質が、エネルギー損失プロセスとは独立して、観測される宇宙線スペクトルを形成する決定的な要因であることを主張している。著者らは以下を強調している：

• 硬いスペクトルの形成： 時間依存の輸送が行われる限り、風によるバブルは、非標準的な加速物理を必要とすることなく、自然に硬い宇宙線スペクトル（$E^{-2}$ より硬いもの）を生成することができる。
• 低エネルギー・カットオフ： バブル内での低エネルギー粒子の「トラッピング（捕捉）」は、この時間依存モデルの堅牢な予測である。これは、バブル内部からのマルチメッセンジャー放射（ガンマ線およびニュートリノ）に観測可能な特徴をもたらし、宇宙線の累積グラマージュ（物質量）にも影響を与える可能性がある。
• 超新星残骸（SNR）との区別： 得られる硬いスペクトルと低エネルギー・カットオフは、超新星残骸（SNR）から通常期待されるものとは異なり、これが「膝（knee）」付近の宇宙線スペクトルや、特定の星団（例：ウェステルンド1）のガンマ線観測を説明する助けとなる可能性がある。

著者らは、本モデルは風駆動フェーズに焦点を当てているものの、捕捉された低エネルギー粒子は後に超新星衝撃波によって再加速されるか、あるいは前方衝撃波が停止した際に脱出し、修正されたスペクトル形状を持つ低エネルギー宇宙線集団に寄与する可能性があると結論づけている。今後の研究では、時間依存のガンマ線およびニュートリノ放射を予測するために、ハドロン相互作用を組み込む予定である。