Revisiting neutrino event epochs for the blazar PKS 0735+178 with TESS

技術要約：TESSを用いたブレイザーPKS 0735+178のニュートリノ・イベント・エポックの再検討

問題と動機
超高エネルギー宇宙線の起源、およびブレイザーにおける高エネルギーニュートリノ放射を駆動する具体的なメカニズムは、依然として高エネルギー天体物理学における未解決の課題である。IceCubeニュートリノ観測所は、特定のブレイザー（例：TXS 0506+056）に関連するニュートリノ事象を特定してきたが、光学的変動性とニュートリノ生成の間の物理的な関連性を理解するには、高密度かつ高頻度の時間データが必要である。ブレイザーPKS 0735+178（$z = 0.45$）は、2021年12月の主要なマルチ波長フレアの期間中、ニュートリノ事象（IceCube-2121208A、$\sim$171 TeV）に対して空間的に近接していることが観測された。これまでの研究では様々な観測装置のデータが利用されてきたが、ニュートリノ検出の前・中・後の各エポックをカバーする、連続的で高頻度な光学ライトカーブ（LC）が欠けていた。本論文は、TESS（Transiting Exoplanet Survey Satellite）のデータを用いてPKS 0735+178の光学的変動を解析することで、この空白を埋め、ニュートリノ放射ウィンドウにおける光源の挙動を特徴付ける。

手法
著者らは、約75日間のフレアを捉えた3つの連続したセクター（44、45、46）におけるTESSの観測データを解析した。

• データ簡約: 本研究では、デトレンド処理されたPDCSAP_FLUXではなく、Simple Aperture Photometry Flux (SAP_FLUX) を使用した。著者らは、デトレンド処理によって、ブレイザーに特有の真の確率論的なフラックス変動が人工的に変化してしまう可能性があると主張している。近傍の微光天体による汚染を軽減するため、Zwicky Transient Facility (ZTF) のデータを用いた差分測光を行い、近傍の天体が観測された変動に寄与していないことを確認した。
• フレアのセグメンテーション: 全局的な最小値と最大値に基づく視覚的検査に基づき、ライトカーブを上昇相（BTJD $\sim$2504.3 から 2558.1）と下降相（BTJD $\sim$2558.1 から 2578.7）に分割した。
• 統計解析: 以下の4つの主要な解析手法を用いた：
  1. フラックス分布: 完全なフレア、上昇相、および下降セグメントのフラックス値のヒストグラムを作成し、正規分布または対数正規分布をテストした。
  2. 過剰分散 ($F_{var}$): 測定誤差を考慮しつつ、固有の変動振幅を定量化するために計算した。
  3. 変動タイムスケール: 上昇相と下降相の特性的なタイムスケールを推定するために、倍増/半減時間のアプローチを適用した。
  4. パワースペクトル密度 (PSD): 両方のフェーズに対して一般化ローソン・スカーグ・ピリオドグラムを計算し、3つのモデル（Simple Power Law (SPL)、Bending Power Law (BPL1)、BPL2）でフィッティングを行った。モデル選択はベイズ情報量基準 (BIC) を用いて決定された。

主な結果

• ライトカーブの形態: この光源は、$\sim$120 $e^-s^{-1}$ からピークの $\sim$375 $e^-s^{-1}$ へのフラックス増加を伴う強いフレアを示し、その後に急速な減衰が続いた。全体の変動振幅は $25.16 \pm 0.01\%$ と算出された。
• フラックス分布: 完全なフレアは多峰性のフラックス分布を示した。上昇相は、明確でほぼ等間隔なガウス型プロファイルを示し、コヒーレントなフラックス状態を示唆した。対照的に、下降相は、特に高フラックス値（$\sim$350 $e^-s^{-1}$）において、歪んだ不規則な分布を示した。
• タイムスケール: 変動タイムスケールの解析により、顕著な非対称性が明らかになった。上昇相は $43 \pm 6$ 日のタイムスケールによって支配されていたが、下降相は有意に速く、$17 \pm 3$ 日のタイムスケールであった。
• スペクトル傾斜: 両方のフェーズ（上昇および下降）のパワースペクトル密度 (PSD) は、Simple Power Law (モデルM1) によって最も良く記述された。スペクトル指数は、上昇相で $2.00 \pm 0.01$、下降相で $1.95 \pm 0.01$ であり、2つのフェーズ間でスペクトル傾斜の公称的な変化があることを示した。
• ニュートリノの相関: Baksan、IceCube、Baikal、および KM3NeT からのニュートリノアラートは、BTJD 2553 から 2564 の間に発生した。IceCube および Baikal のアラートは光学最大値（$\sim$BTJD 2558.3）と密接に一致しており、KM3NeT のアラートは下降相の間に発生した。

物理的解釈と意義
著者らは、以前の超長基線干渉計（VLBI）によるラジオ観測結果（Y.-S. Kim & J.-Y. Kim 2025）——新しいジェット成分（C2）がより遅い成分（C1）を横切って放出されたことを特定したもの——に基づき、光学的な挙動とニュートリノ事象を結びつける物理的シナリオを提案している。

• 上昇相: 多峰性のフラックス分布と上昇相における明確な倍増タイムスケールは、コアからの規則的かつ離散的なプラズマの注入を示唆している。これらの注入が集合して、ジェット内を伝播するC2成分を形成している。
• 下降相: 急速な減衰と不規則なフラックス分布は、C2とC1の相互作用に起因すると考えられる。著者らは、この相互作用がコンパクトな成分を、ランダムに移動する様々なサイズのプラズマ断片のアンサンブルへと崩壊させると示唆している。この断片化は、（視線方向の変化による）ドップラー因子の減少と急激なフラックス減少をもたらす。
• ニュートリノとの接続: ニュートリノアラートの光学最大値およびC2-C1相互作用との時間的一致は、この特定の相互作用の窓において、粒子加速とターゲットフォトン密度の増加により、ニュートリノ放射が促進される可能性を示唆している。

主張と限界
本論文は、TESSのデータがこのフレアに対する最も高密度な光学サンプリングを提供し、ニュートリノ放射エポックに対する重要な時間的制約を与えると主張している。著者らは、観測された光学フラックスの挙動が、ニュートリノ生成に資する条件（例：強化された粒子加速）と矛盾しないことを挙げ、光学変動とニュートリノ事象の間の有力な関連性を提唱している。

しかし、著者らは控えめな姿勢を維持している。彼らは以下の点を明示している：

1. 統計的な有意性が限定的であり、同時期のマルチメッセンジャー制約が欠如しているため、この関連性は依然として暫定的なものである。
2. 既知のニュートリノ放射ブレイザーのサンプルは、現在非常に少ない。
3. 偶然の一致を含む代替の説明を排除することはできない。
4. ハドロンシナリオはスペクトルエネルギー分布を再現できるが、予測されるニュートリノフラックスはモデル依存であり、レプトンモデルも電磁放射を説明できる可能性がある。

したがって、本論文は、より大きなサンプルとなるブレイザーのさらなる観測的証拠を待つ間、ニュートリノ放射のリンクは「確定的」ではなく「示唆的」なものとして扱うべきであると結論付けている。