[論文レビュー] The Multi-Mission Maximum Likelihood framework (3ML)

G. Vianello, R. Lauer+8 more|arXiv (Cornell University)|Jul 29, 2015

Astrophysics and Cosmic Phenomena参考文献 1被引用数 94

ひとこと要約

Multi-Mission Maximum Likelihoodフレームワーク（3ML）は、オープンソースでプラグインベースのソフトウェアフレームワークであり、複数の機器および波長にわたる天体源の広帯域スペクトルおよび空間的モデリングを一貫性を持って統合的に実行できる。これは、各機器の尤度に基づく解析パイプラインを統合することで実現される。3MLは、機器固有の応答関数とデータ形式を用いて、点源と拡張源を同時にフィッティング可能であり、余計なパラメータや相互キャリブレーションを透明に処理するため、ユーザーの干渉を最小限に抑えつつ、マルチメッセンジャーおよびマルチウェーブレングス解析を可能にする。

ABSTRACT

Astrophysical sources are now observed by many different instruments at different wavelengths, from radio to high-energy gamma-rays, with an unprecedented quality. Putting all these data together to form a coherent view, however, is a very difficult task. Each instrument has its own data format, software and analysis procedure, which are difficult to combine. It is for example very challenging to perform a broadband fit of the energy spectrum of the source. The Multi-Mission Maximum Likelihood framework (3ML) aims to solve this issue, providing a common framework which allows for a coherent modeling of sources using all the available data, independent of their origin. At the same time, thanks to its architecture based on plug-ins, 3ML uses the existing official software of each instrument for the corresponding data in a way which is transparent to the user. 3ML is based on the likelihood formalism, in which a model summarizing our knowledge about a particular region of the sky is convolved with the instrument response and compared to the corresponding data. The user can choose between a frequentist analysis, and a Bayesian analysis. In the former, parameters of the model are optimized in order to obtain the best match to the data (i.e., the maximum of the likelihood). In the latter, the priors specified by the user are used to build the posterior distribution, which is then sampled with Markov Chain Monte Carlo or Multinest. Our implementation of this idea is very flexible, allowing the study of point sources as well as extended sources with arbitrary spectra. We will review the problem we aim to solve, the 3ML concepts and its innovative potential.

研究の動機と目的

• 異なる望遠鏡から得られる多機器・多波長天体データを、互換性のないデータ形式、ソフトウェア、分析手順を有する状況でも、一貫した形で統合的に分析する課題に対処すること。
• 機器の出どころやデータ形式にかかわらず、一貫した広帯域スペクトルフィッティングを可能にする共通のフレームワークを提供すること。
• 各機器の公式ソフトウェア（OS）の整合性と正確性を保つために、それを変更せずにプラグインを通じて透明に使用することで、再実装を回避すること。
• 1つの統合されたモデリング環境内で、頻度主義的およびベイズ的統計的分析を両立させること。
• 背景の正規化や相互キャリブレーション係数などの余計なパラメータを、プロファイル尤度法を用いて体系的に取り扱うこと。

提案手法

• 3MLは、スカイモデル S(E, P) が、機器固有の応答関数 R(E, e, P, p) と畳み込まれ、観測されたカウントを予測する尤度形式を用いる。
• フレームワークは、各機器の公式ソフトウェア（OS）に直接インターフェースするプラグインアーキテクチャを採用しており、変更を加えずに機器固有の尤度を計算する。
• 総尤度は、各機器からの個別の尤度の合計として構築され、複数のデータセットにわたる同時フィッティングを可能にする。
• ベイズ的分析では、事前分布を尤度と組み合わせて事後分布を形成し、MCMC や MultiNest を用いてサンプリングする。頻度主義的分析では、パラメータは MIGRAD を用いて最適化される。
• 余計なパラメータ（例：背景の正規化）は、プロファイル尤度法を用いてプロファイルアウトされ、システム的な不確実性を考慮しつつ、源のパラメータに焦点を当てる。
• フレームワークは、任意のスペクトル形状を持つ点源および拡張源のモデリングをサポートしており、ニュートリノまたは重力波の放射率予測を含める形で拡張することで、マルチメッセンジャーアナリシスを可能にする。

実験結果

リサーチクエスチョン

• RQ1 異なるデータ形式やソフトウェアパイプラインを持つ多機器・多波長天体データを、互換性のない状況下でも、一貫したスペクトルおよび空間的分析に統合的に統合する方法は何か？
• RQ2 既存の機器固有の公式ソフトウェアを、変更や再実装を加えずに、統一された分析フレームワークに統合する最適な方法は何か？
• RQ3 背景の正規化や相互キャリブレーション係数などの余計なパラメータを、統合尤度フレームワーク内で体系的に取り扱う方法は何か？
• RQ4 1つのモデリングフレームワークが、電波からガンマ線に至る広帯域にわたるエネルギー範囲で、点源と拡張源の両方の解析をサポートできるか？
• RQ5 3MLは、ニュートリノや重力波などの非電磁気的メッセンジャーからのデータ統合を通じて、マルチメッセンジャーアストロノミーをどの程度促進できるか？

主な発見

• 3MLは、Fermi/LAT、HAWC、VERITAS、Swift、HESS といった複数の機器からのデータを用いて、広帯域スペクトルフィッティングを一貫して実行でき、多様なデータタイプにおいても一貫した結果を得た。
• プラグインアーキテクチャにより、各機器の公式ソフトウェアへの完全な透明性と忠実性が確保され、精度を損なわず、再実装の必要がない正確なデータ処理が可能になった。
• フレームワークは頻度主義的およびベイズ的推論を両方サポートしており、MCMC や MultiNest を用いてパラメータの不確実性を完全にサンプリング可能である。
• Fermi/LAT や HAWC における背景正規化などの余計なパラメータは、効果的にプロファイルアウトされ、バイアスのない源パラメータ推定が可能になった。
• 3MLは、広視野と高解像度望遠鏡を組み合わせた研究を可能にする、複数の機器にわたる拡張源の最初の統合的モデリングを実現した。
• フレームワークはマルチメッセンジャーアストロノミーに拡張可能であり、ニュートリノおよび重力波探査機からのデータ統合への明確な道筋が、互換性のある尤度プラグインを通じて確保されている。