[論文レビュー] A Comprehensive Statistical Description of Radio-Through-$γ$-Ray Spectral Energy Distributions of All Known Blazars

Peiyuan Mao, C. M. Urry+4 more|arXiv (Cornell University)|Apr 13, 2016

Astrophysics and Cosmic Phenomena被引用数 3

ひとこと要約

本研究では、2214個の既知のブラザールに対して、モデルに依存しないアプローチを用いて、電波からγ線に至るスペクトルエネルギー分布（SED）の包括的な統計的分析を実施した。主成分分析（PCA）を用い、欠損したフラックスをk近傍法で補完することで、電波の放射度とコヒーレントピーク周波数の間に顕著な逆相関があることが確認された。SEDの形状は、選択効果とは無関係に、主に1.4 GHz電波放射度と赤方偏移に依存していることが示された。

ABSTRACT

We combined multi-wavelength data for blazars from the Roma-BZCAT catalog and analyzed hundreds of X-ray spectra. We present the fluxes and Spectral Energy Distributions (SEDs), in 12 frequency bands from radio to $γ$-rays, for a final sample of 2214 blazars. Using a model-independent statistical approach, we looked for systematic trends in the SEDs; the most significant trends involved the radio luminosities and X-ray spectral indices of the blazars. We used a Principal Component Analysis (PCA), to determine the basis vectors of the blazar SEDs and, in order to maximize the size of the sample, imputed missing fluxes using the K-nearest neighbors method. Using more than an order of magnitude more data than was available when Fossati et al. (1997, 1998) first reported trends of SED shape with blazar luminosity, we confirmed the anti-correlation between radio luminosity and synchrotron peak frequency, although with greater scatter than was seen in the smaller sample. The same trend can be seen between bolometric luminosity and synchrotron peak frequency. Finally, we used all available blazar data to determine an empirical SED description that depends only on the radio luminosity at 1.4~GHz and the redshift. We verified that this statistically significant relation was not a result of the luminosity-luminosity correlations that are natural in flux-limited samples (i.e., where the correlation is actually caused by the redshift rather than the luminosity).

研究の動機と目的

• 電波からγ線に至る周波数範囲において、既知のすべてのブラザールのSEDを、統計的に信頼性が高く、モデルに依存しない記述を提供すること。
• パラメトリックな仮定を排除して、SED形状の系統的傾向が放射度と赤方偏移に関連しているかどうかを特定すること。
• SED形状と放射度の間の観察された相関が物理的要因によるものか、それとも赤方偏移に依存するフラックス限界などの選択効果に起因するのかを検証すること。
• サンプルの完全性と選択バイアスを検証した上で、1.4 GHz電波放射度と赤方偏移のみに依存する経験的SED推定器を開発すること。
• 赤方偏移に起因する相関を除去した上で、全放射度とコヒーレントピーク周波数の間の相関が本質的であるかどうかを評価することにより、ブラザールシーケンスの本質的性質を検証すること。

提案手法

• ローマ-BZCATカタログを用いて、2214個のブラザールの電波からγ線に至る12の周波数帯域における多波長フラックスデータを収集した。
• 関数形を仮定せずに、SED変動を記述する主要な基底ベクトルを特定するために主成分分析（PCA）を用いた。
• サンプルサイズと統計的パワーを最大化するために、欠損フラックス値をk近傍法（k-NN）を用いて補完した。
• 観測された相関が赤方偏移に起因する選択効果によるものかどうかを評価するために、広範なブートストラップ法と赤方偏移の入れ替えテストを実施した。
• 宇宙論的パrameter H₀ = 70 km s⁻¹ Mpc⁻¹ および ΩΛ = 0.72 を用いて、各周波数での放射度を計算し、log₁₀ベースの対数を用いた。
• フラックス限られたサンプルに内在する放射度-放射度相関に依存しないことを検証することで、経験的SEDモデルの妥当性を検証した。

実験結果

リサーチクエスチョン

• RQ1 電波放射度とコヒーレントピーク周波数の間の逆相関は、物理的傾向であるのか、それとも赤方偏移に起因する選択効果によるものなのか？
• RQ2 SED形状が1.4 GHz電波放射度と赤方偏移に依存する程度は、他のパラメータとは無関係にどの程度高いのか？
• RQ3 赤方偏移に起因する放射度相関を除去した後でも、全放射度とコヒーレントピーク周波数の間の相関が維持されるのか？
• RQ4 高コヒーレントピークを持つブラザール（HSP）は、主な調査で系aticallyに不足しているのか。その結果、観測されたSED傾向に影響を与えているのか？
• RQ5 他のSEDパラメータに依存せず、1.4 GHz放射度と赤方偏移のみに依存する、統計的に有意な経験的SEDモデルを構築できるのか？

主な発見

• 電波放射度とコヒーレントピーク周波数の間の逆相関は、以前の研究よりも統計的パワーが向上したが、より大きな不均一なサンプルのため、散乱が大きくなった。
• 全放射度とコヒーレントピーク周波数の間の相関は、統計的に有意であり、赤方偏移の入れ替えとブートストラップテストを厳密に実施した後でも持続した。これは、物理的根拠があることを示している。
• 赤方偏移の入れ替え後に、測定された相関係数と勾配が偶然に得られる確率は1%未満であり、相関が単に赤方偏移バイアスによるものであるという仮説は強く否定された。
• 1.4 GHz電波放射度と赤方偏移にのみ依存する経験的SEDモデルが導出され、フラックス限られたサンプルに内在する放射度-放射度相関とは独立した統計的有意性が検証された。
• 本研究は、SED形状と放射度を結ぶブラザールシーケンスが、調査の選択バイアスの産物ではなく、本質的な物理的起源を持つことを示した。
• これらの結果は、将来のモンテカルロシミュレーションの基盤を提供し、HSPブラザールがLSPブラザールよりも本質的に多数存在するかどうかを検証するための鍵となるテストを可能にする。