X-ray Fourier lag-frequency spectra modulated by stochastic turbulent acceleration in the jets of high-frequency-peaked BL Lac

技術要約：HBLジェットにおける確率論的乱流加速によって変調されるX線フーリエ・ラグ周波数スペクトル

問題提起
X線帯域間のタイムラグは、ジェット物理学の重要な診断指標である。観測では、高エネルギーの変化が低エネルギーの変化に対して遅れるハード（正）ラグや、その逆が起こるソフト（負）ラグを含む、多様なラグのシグネチャが頻繁に示される。フーリエ変換された電子輸送方程式に基づく従来の理論モデルは、ソフトラグを再現したり、ハードラグを説明するために一般化された輸送方程式を必要としたりしたが、非線形なシンクロトロン自己コンプトン（SSC）冷却を無視していることが多い。また、特に、これらと連動してラグ周波数スペクトルの形状を決定する確率論的乱流加速（STA）の役割については、十分に探索されてこなかった。本研究は、HBLのフレアで見られる複雑なタイムラグのシグネチャを解釈するために、STA、放射冷却、粒子脱出、および非線形SSC冷却を考慮した統一的な定量的フレームワークの必要性に対処するものである。

手法
著者らは、相対論的ローレンツ因子を持つ、空間的に一様で光学的厚さが薄い球状のプラズマ塊を仮定した、一般的な1ゾーン・レプトニック・モデルを用いている。電子集団の固有の時間進化は、時間依存のフォッカー・プランク（FP）方程式によって支配される。この方程式には以下が含まれる：

1. 注入： 衝撃的に注入されるべき電子の冪乗分布。
2. 加速： エネルギー空間における拡散係数（$D_p \propto \gamma'^q$）を介してモデル化された確率論的乱流加速（STA）。ここでは「ハードスフィア」型乱流（$q=2$）を主要なケースとしている。
3. 冷却： シンクロトロン放射およびSSC散乱による放射損失。後者には非線形効果とクライン＝ニシナ補正が含まれる。
4. 脱出： エネルギーに依存しない脱出タイムスケール。

FP方程式は、インプリシット（隠解法）なクランク・ニコルソン・スキームを用いて数値的に解かれる。得られた時間依存の電子分布を用いて、シンクロトロンおよびSSC放射スペクトルを計算する。因果律を考慮するため、放出領域のスライスからのフラックス寄与を積分することによって、光学的通過効果（light-crossing effects）を組み込む。

ソフトX線（0.5–2 keV）とハードX線（2–10 keV）のライトカーブ（LC）間のタイムラグは、以下の3つの手法を用いて分析される：

• 離散フーリエ変換（DFT）： フーリエ・ラグ周波数スペクトルを導出するため。
• 正規化フラックス・マッチング遅延（NFMD）： 瞬時的な時間遅延を追跡するため。
• 相互相関関数（CCF）： ライトカーブ全体にわたるセントロイド・ラグを測定するため。

本研究では、最大注入ローレンツ因子（$\gamma'_{i,max}$）、スペクトル指数（$p$）、拡散係数（$D_0$）、および注入率を系統的に変化させ、ハードラグ、ソフトラグ、および遷移状態の3つの異なるレジームをシミュレートしている。

主な貢献と結果

1. ラグ・レジームの統一： 本モデルは、以下の3つの異なるタイムラグ・レジームを再現することに成功した：

   • ハード／正のラグ： 最大注入エネルギーが平衡エネルギーを大幅に下回る場合（$\gamma'_{i,max} \ll \gamma'_{eq}$）。これは加速支配のダイナミクスを示している。
   • ソフト／負のラグ： $\gamma'_{i,max} \gg \gamma'_{eq}$ の場合であり、冷却が支配的となる。
   • 遷移レジーム： $\gamma'_{i,max} \approx \gamma'_{eq}$ のとき、周波数に対して正から負への滑らかな遷移が観察される。この遷移は、加速支配の低エネルギー電子と冷却支配の高エネルギー電子の共存を反映している。

2. STAの役割： STA、放射冷却、および脱出プロセスの競合が、ラグ周波数スペクトルを形成する。具体的には、STAは、低エネルギーの電子を継続的に再加速することにより、高エネルギー電子の放射冷却を抑制する。この抑制により、純粋な冷却モデルと比較して、遷移レジームおよびソフトラグ・レジームにおけるタイムラグが増幅される。

3. 非線形SSC冷却の影響：

   • SSC冷却は、特に遷移レジームおよびソフトラグ・レジームにおいて、タイムラグをさらに増幅させる。
   • この増幅は、STAが（再加速された低エネルギー電子を介して）ソフトX線放射を維持する一方で、SSC冷却がハードX線放射の減衰を加速させることに起因する。
   • モデルは、フレアの継続時間とラグの振幅の間に正の相関を予測している：フレアが長くなると電子の蓄積が進み、SSC冷却が強まるため、より大きなラグが生じる。これは、TeV輝度の高いフレアで見られる大きなラグを説明するものである。

4. 手法の比較： 本研究は、ライトカーブ全体を平均化するCCFラグが、フーリエ・ラグ周波数スペクトルによって明確に解明される遷移的挙動や周波数依存の変動を覆い隠してしまう可能性があることを強調している。非対称なライトカーブ（ソフトラグのケース）では、CCFラグは低周波のフーリエ・ラグよりも著しく小さくなることがある。

5. 堅牢性： ラグ周波数スペクトルの定性的な特徴は、乱流指数を「ハードスフィア」型（$q=2$）からコルモゴロフ型（$q=5/3$）へと変化させても一貫しているが、ラグの振幅や特性周波数には定量的な差異が生じる。フェルミI型（系統的）加速の導入は、ハードラグのケースでは無視できる影響であったが、他のレジームではスペクトルを修正する可能性がある。しかし、STAおよびSSCに関する主要な結論は堅牢である。

意義
本論文は、HBLのX線フレアで見られる多様なタイムラグのシグネチャを解釈するための、統一的な定量的フレームワークを提供すると主張している。STAと非線形SSC冷却を統合することで、本モデルは以下を説明する：

• ハードラグとソフトラグの両方の存在、およびそれらの間の遷移。
• TeV輝度の高いフレアにおけるラグの増幅。
• より大きなフレアの継続時間が、より大きなタイムラグに対応するという観察された傾向。

著者らは、正から負のラグへの滑らかなシフトを特徴とする遷移レジームが、相対論的電子の再加速の証拠となり得ることを示唆している。結論として、フーリエ・ラグ周波数スペクトルは、スペクトルエネルギー分布（SED）モデリングと組み合わせることで、相対論的ジェットのコンパクトな放射領域における極限的な物理条件（加速効率、磁場強度、および冷却率）を制約するための強力な診断ツールとなる。本研究は、優れた時間分解能とエネルギー分解能を持つAthenaのような将来のミッションが、これらの理論的予測を検証するための精密なスペクトル測定を可能にすることを期待している。