A Boundary-Consistent Two-Zone Electron Kernel for Distant Pulsar Contributions to Positron Flux and Anisotropy

技術要約：遠方のパルサーによる陽電子フラックスおよび異方性への寄与に関する、境界整合的な二層電子カーネル

問題提起
局所的な宇宙線電子および陽電子のスペクトルは、近傍の高エネルギー加速器、特にパルサーを特定するための精密なプローブとして機能する。GemingaやMonogemのような成熟したパルサーは、TeVスケールのレプトン・フラックスの自然な候補であるが、ここで重要な問いが残る。すなわち、局所的なレプトン・フラックスは最も近いパルサーによって支配される必要があるのか、それとも、より遠方のパルサーが集合的に寄与し得るのか？この問いは、2つの要因によって複雑化している。(1) 急激な放射冷却（シンクロトロンおよび逆コンプトン散乱）が、高エネルギー粒子の伝播ホライゾンを制限すること、(2) Fermi-LATおよびAMS-02による異方性の制限が、単一の近傍ソースが広範囲にわたって支配するモデルを制約していることである。先行研究では、遠方のソースが低エネルギーの電子を再現できることが示唆されているが、近傍ソースにおける低速拡散（slow diffusion）が存在する場合の幾何学的なソースカウント・ファクターを定量化するための、安定した半解析的なフレームワークは欠如していた。さらに、銀河系の内側数百パーセク（ローカルバブル）は、分子ガス密度の低さやパルサーの誕生時のキックを考慮すると、滑らかな連続的なパルサー密度を想定するには物理的な不連続点が存在する。

手法
著者らは、球状の二層拡散モデルにおける電子および陽電子伝播のための、半解析的な級数解を開発した。

• 二層カーネル： このモデルは、ソースを取り囲む半径 $R_1$ の低速拡散領域（拡散係数 $D_1(E)$）と、標準的な星間領域（$R_1 < r < R_2$、拡散係数 $D_2(E)$）を定義する。各係数は同じエネルギー依存性（$D \propto f_D(E)$）を共有しており、これにより解析的な固有関数解が可能となる。
• 境界条件： 解は、$R_1$ における密度と拡散フラックスの連続性を、および $R_2$（20 kpcに設定）における吸収境界を強制する。この定式化は、二層モデルの直接的な積分評価に内在する振動的な打ち消しを回避し、鋭い放射冷却境界を保持する。
• 冷却： 放射損失（シンクロトロンおよびCMB、赤外、可視光場に対するクライン＝ニシナ逆コンプトン散乱）は、注入エネルギーとソースの年齢を関連付けるエネルギー特性を用いて処理される。
• 集団モデリング： カーネルは、連続的な円盤状のパルサー集団に適用される。著者らは、ローカルバブルとパルサーの固有運動を考慮するため、$R < 0.1$ kpcに対して空洞（empty cavity）を設けている。
• 統計解析： AMS-02の陽電子フラックス・データ、CALETの全電子ダイポール異方性測定値、およびFermi-LATの異方性上限値にモデルを適合させるため、「ソルバー・イン・ザ・ループ」型のマルコフ連鎖モンテカルロ（MCMC）解析が行われた。このフィットには、連続的なスピンダウン注入履歴、破断冪乗則を持つ対生成注入スペクトル、および二次成分的な陽電子基底と滑らかな冪乗則背景の非負線形振幅上のプロファイルが含まれる。

主要な貢献

1. 安定した級数解： 本論文は、放射冷却を伴う二層拡散問題に対する、数値的に安定した半解析的な級数解を提供する。この手法は計算効率が高く、特に中距離領域において、ブルートフォース的な有限体積法や振動積分が抱える数値的不安定性を回避する。
2. 幾何学的定量化： 本研究は、「シェル・ファクター」（距離に対するソースカウント）が伝播カーネルとどのように相互作用するかを定量化している。3Dの一様分布では10–100 GeVにおいてキロパーセック規模のシェルが支配的であるが、円盤形状（$h=0.2$ kpc）ではその寄与が内側へとシフトするものの、1 kpcを超えるソースが依然として10–100 GeVのフラックスの37–47%を供給していることを示している。
3. 低速ハローの役割： 著者らは、低速拡散ハローを含めることは、一層モデルと比較して遠方のソースからのフラックスを「増強」するのではなく、むしろエネルギー依存的な冷却ペナルティを導入し、すべての距離からの絶対的なフラックスを減少させ、高エネルギーにおいてより強い抑制をもたらすことを明らかにしている。
4. 異方性制約： 本研究は、CALETのダイポール振幅とFermi-LATの上限値を用いて、共通のハロー・パラメータの妥当性をテストすることにより、異方性制約を集団フィットに直接統合している。

結果

• フラックス分解： フィットされた連続的なパルサー集団（$R < 0.1$ kpcを除く）において、0.3 kpcを超えるソースからの寄与は、10 GeVからTeVスケールに至るすべてのエネルギーにおいて、局所的な0.1–0.3 kpcのアニュラス（環状領域）を支配している。10 GeVにおいて、$R > 0.3$ kpcの寄与は、局所アニュラスよりも約1桁大きい。
• ハロー・パラメータ： MCMC解析により、低速ハロー半径 $R_1$ の事後分布は $0.0709^{+0.0510}_{-0.0407}$ kpc（$\approx 71$ pc）を中心とし、拡散係数は $D_1(100 \text{ TeV}) \approx 2.08 \times 10^{28}$ cm$^2$ s$^{-1}$ となった。この結果は、GemingaおよびMonogemのTeVハロー形態の研究と一致している。
• 異方性の適合性： 導出されたハロー・パラメータを用いたフィットモデルは、現在の異方性制限（CALETおよびFermi-LAT）と両立している。ダイポール振幅の予測値は95%の上限値内に収まっている。
• エネルギー依存性： 遠方のソースはサブTeVのフラックスを支配するが、放射冷却によってその寄与はTeVスケール付近に局在化する。1 TeVにおいて、円盤モデルのフラックスの約90%は1 kpc以内から来ている。

意義と主張
本論文は、遠方のパルサー（直近のローカルバブル近傍を超えたもの）からの集合的な寄与が、数百GeV以下の陽電子および電子フラックスの、妥当かつ支配的な寄与因子となり得ることを結論づけている。フィットされたパルサー成分は、最も近いソースではなく、0.3 kpc以上に広がる集団によって支配されている。

しかし、著者らは自らの知見の独自性については慎重な立場を維持している：

• フラックスおよび異方性のデータのみでは、ハローのサイズ（$R_1$）や拡散係数（$D_1$）を一意に決定することはできない。注入形状、背景モデル、放出履歴、およびカタログの不完全性との間に縮退が存在する。
• 導出された「Gemingaスケール」のハローサイズは、使用された特定の集団モデルとデータセットに基づく条件付きの制約であり、物理的な低速拡散ハローの確定的な測定値ではない。
• 共通のハロー・パラメータの検証から、低速拡散領域の決定へと進むためには、個々のハローのTeV逆コンプトン表面輝度プロファイル、詳細なパルサー集団データ（固有運動、誕生場所）、および完全なベクトル異方性尤度が、追加の観測量として必要である。

本研究は、伝播制御された実証実験として、現在のデータが共通ハローのパラメータ空間をテストできることを示しており、将来の集団推論のための安定したカーネルを提供すると同時に、完全な解決にはマルチメッセンジャーの入力が必要であることを強調している。