Multiwavelength variability of the high-energy neutrino candidate PKS 0735+178 over three decades

技術要約：30年間にわたる高エネルギーニュートリノ候補天体 PKS 0735+178 のマルチ波長変動

問題提起
PKS 0735+178 は、高エネルギーニュートリノ・イベント IC211208A (171 TeV) と関連付けられている典型的な BL Lac 天体である。この天体は極端なマルチ波長変動と、小さな角度で観測される相対論的ジェットで知られているが、異なる周波数における放射領域間の物理的な結合、長期的な準周期的挙動の性質、およびニュートリノ・イベントを取り巻く活動を駆動する具体的なメカニズムについては、依然として十分に理解されていない。先行研究は主に短期間のスペクトル特性や孤立したフレアに焦点を当てており、ジェットの幾何学構造や擾乱の伝播を制約するための、マルチデケード（数十年間）のベースラインにわたる光源の時系列的進化に関する統一的な解析を欠いている。

手法
著者らは、以下のデータを統合した約30年間（1995年–2026年）の包括的なマルチ波長解析を提示している：

• 電波データ： RATAN-600 電波望遠鏡による 1–22 GHz の準同時観測、および Submillimeter Array (SMA) による 230 GHz データ。
• 光学データ： Zeiss-1000、AS-500/2、KAIT、ZTF、および CRTS の R バンド測光。
• ガンマ線データ： フレアリング状態における微細構造を解像するために適応的ビンニングを用いた Fermi-LAT ライトカーブ (1.5–300 GeV)。
• VLBI 制約： ジェットの運動学および物理パラメータを導出するための、MOJAVE および VLBA-BU-BLAZAR データを用いた 15 および 43 GHz でのコアシフト測定。

解析には、以下の統計的および物理的モデリング手法を用いた：

1. 相互相関： 離散相関関数 (DCF) を用いて、バンド間のタイムラグを測定。
2. スペクトル解析： 電波スペクトルを3つの主要なフレア・エピソードにグループ化し、スペクトル指数 ($\alpha$) の進化を追跡。
3. 準周期性探索： 重み付きウェーブレット Z 変換 (WWZ) および 一般化ロンスガード (GLS) ピリオドグラムを適用して特徴的なタイムスケールを検出し、モンテカルロ・シミュレーションを用いて赤ノイズ背景に対する有意性を評価。
4. 物理モデリング：
   • 衝撃波伝播： タイムディレイをべき乗則依存性 ($\Delta t \propto \nu^{-1/k_r}$) に適合させ、不透明度で層状化されたジェットモデルを検証。
   • ジェット歳差運動： ベイジアン最適化モデル (UltraNest) を用い、幾何学的効果が長期的な変調を説明できるかどうかを検証するために、ライトカーブを歳差運動シナリオに適合。
   • VLBI コアシフト： ジェットの視線角、見かけの速度、磁場強度、および磁気力線に対する磁束（磁気阻止円盤 [MAD] 極限に対する比率）を推定するために使用。

主な結果

• 周波数依存の時間遅延： 電波、光学、および $\gamma$ 線バンドのライトカーブは相関しているが、電波放射は高周波に対して遅れて現れる。遅延は低周波数に向かって系統的に増加し、$\gamma$ 線に対して 1–2 GHz では $\sim$1200 日に達する。これは $k_r = 1.96 \pm 0.24$ となるべき乗則に従っており、擾乱が下流へと伝播する不透明度で層状化されたジェットと一致している。
• ニュートリノ・イベントのフレア： IC211208A イベントの数日以内に、光学および $\gamma$ 線バンドで短く明るいフレアが発生した。光学が先行する 12 日間のラグ が光学と $\gamma$ 線ピークの間で検出されており、これは非共空間的な放射領域、または外部コンプトン散乱のようなプロセスの寄与を示唆している。第3の主要な電波フレア (2019–2026) は、低周波数において逆転または平坦なスペクトルを示しており、粒子加速が強化されたコンパクトで光学的厚みを持つ放射領域を示している。
• 準周期的な変動： 電波およびミリ波 (5–230 GHz) において、特徴的な変動タイムスケール $\sim$10–11 年 がロバストに検出された ($\ge 3\sigma$)。この変調はこれらのバンド間でコヒーレントであるが、厳密な周期性ではない。対照的に、光学および $\gamma$ 線データは、時間的に局在化しており、コヒーレンスが低い $\sim$5–9 年のより短い周期の信号を示している。
• ジェットの幾何学と物理：
  • VLBI コアシフト測定により、ジェットの視線角 $\theta \approx 3.6^\circ$ および見かけの速度 $\beta_{app} \approx 0.8c$ が得られ、これはタイムディレイから導出された伝播速度と一致している。
  • 1 pc 地点での磁場強度は $\sim$0.36 G と推定される。磁束と MAD 極限の比は、現在の光度レベルでは降着円盤が MAD 状態にはないことを示唆している。
• 歳差運動モデルの限界： ジェットの歳差運動モデルはライトカーブの一般的な形状を再現できるが、導出されたパラメータ（視線角 $\sim$27.4°）は VLBI の制約 ($\theta \sim 1^\circ-3^\circ$) と矛盾する。これは、歳差運動単独では観測された変動を説明できないことを示している。

意義および主張
本論文は、PKS 0735+178 の変動が、異なるタイムスケールで作用する2つの異なる物理プロセスによって支配されていると主張している：

1. 長期的な変調： ジェット基部のエネルギー放出における緩やかな変化が、$\sim$11 年の準周期的な変調を駆動しており、これは降着流またはジェット基部のダイナミクスに関連し、擾乱の注入率を制御していると考えられる。
2. 短期的フレア： ニュートリノ・イベントに関連する個々のフレアは、不透明度で層状化されたジェット内を伝播する衝撃波から生じる。周波数依存の遅延とスペクトル進化は、衝撃波が異なる距離でシンクロトロン光学的厚さの表面（フォトスフィア）を通過することによって自然に説明される。

本研究は、この天体における光学と $\gamma$ 線放射の間の 12 日間のラグを初めて信頼性高く測定したものであり、単純な一ゾーン・シンクロトロン・セルフ・コンプトン（SSC）シナリオに疑問を投げかけるものである。さらに、VLBI コアシフトと電波タイムディレイからの独立した制約の一致は、層状流における衝撃波伝播という統一的なシナリオを強く支持している。ニュートリノ・イベントのマルチ波長フレアへの時間的な集積は、PKS 0735+178 と高エネルギーニュートリノ生成との関連性を補強しており、その生成はジェット内のより内側の、光学的厚みを持つ領域で行われている可能性が高い。