[論文レビュー] DiskFit: a code to fit simple non-axisymmetric galaxy models either to photometric images or to kinematic maps

J. A. Sellwood, Kristine Spekkens|arXiv (Cornell University)|Sep 23, 2015

Astronomy and Astrophysical Research被引用数 24

ひとこと要約

DiskFit は、非軸対称な銀河モデル（たとえば棒やねじれ）を、全領域のカイ二乗最小化手法を用いて、光度的画像または運動学的マップにフィットする公開ソフトウェアパッケージである。ブートストラップリサンプリングを用いて統計的に妥当な不確実性を提供し、回転曲線やgalfitが仮定するような傾いたリングやパラメトリックなプロファイルを避けることで、データ駆動型の、頑健なディスク銀河構造および運動学のモデリングを可能にする。

ABSTRACT

This posting announces public availability of version 1.2 of the DiskFit software package developed by the authors, which may be used to fit simple non-axisymmetric models either to images or to velocity fields of disk galaxies. Here we give an outline of the capability of the code and provide the link to downloading executables, the source code, and a comprehensive on-line manual. We argue that in important respects the code is superior to rotcur for fitting kinematic maps and to galfit for fitting multi-component models to photometric images.

研究の動機と目的

• 非軸対称な特徴（たとえば棒）を、光度的および運動学的データの両方に対して柔軟かつオープンソースでフィットできるツールの開発。
• rotcur（傾いたリング） や galfit（パラメトリックなプロファイル）といった既存ツールの制限を克服するため、非パラメトリックで全領域のフィッティングアプローチを採用すること。
• ブートストラップ手法によって空間相関を持つ残差を保持することで、フィットパラメータの統計的に妥当な不確実性推定を提供すること。
• 棒、ねじれ、バルジといった複雑な特徴をフィット可能にしつつ、計算の柔軟性と正確さを維持すること。
• ユーザー定義の半径サンプリングと見かけのぼやけ補正を許容することで、銀河構造および運動学の高精度モデリングを支援すること。

提案手法

• 2次元観測データ（画像または速度場）と、表形式で与えられる非パラメトリックモデル（たとえば回転曲線または光度プロファイル）を、全領域のカイ二乗最小化で比較する。
• 空間相関を持つ残差を保持するブートストラップリサンプリング技術を用いて、フィットパラメータの現実的な不確実性推定を実現する。
• 径方向および方位角方向の速度成分における m=2（二重対称）歪みの重ね合わせにより、非軸対称な特徴（たとえば棒）をサポートする。
• 軸対称モデルのフィッティングを可能とし、必要に応じてねじれを含める。ねじれの振幅は半径の2乗に比例し、固定された線のノードの位置角によって制約される。
• 非軸対称特徴が投影された主要軸とほぼ一致する場合に、フィットを安定化させるためにスムージングペナルティをオプションで導入する。
• ユーザー指定の FWHM を用いたガウス畳み込みを適用して、観測のぼやけ効果をモデル化する。最大 FWHM は 3 ピクセルまで許容される。

実験結果

リサーチクエスチョン

• RQ1 パラメトリックなプロファイルを仮定せずに、非軸対称な特徴（たとえば棒）を銀河の運動学的および光度的データに頑健にフィットする方法は何か？
• RQ2 全領域モデルフィッティングは、rotcur のような傾いたリング法に比べて、回転曲線および非円運動の推定においてどのように優れているか？
• RQ3 残差が空間的に相関している場合、フィットパラメータの不確実性推定を統計的に妥当にするにはどうすればよいか？
• RQ4 投影された主要軸とほぼ一致する棒をフィットする際の制限とは何か。また、それらはどのように緩和できるか？
• RQ5 DiskFit は、ねじれやバルジといった複雑な特徴を、正確さと柔軟性を保ちつつどの程度モデル化できるか？

主な発見

• DiskFit は、全領域最小化フレームワークを用いて、光度的および運動学的データの両方に対して、非軸対称な特徴（たとえば棒やねじれ）を成功裏にフィットした。
• ブートストラップリサンプリングによって残差の空間相関を保持することで、標準的手法よりも現実的であるとされる誤差バーを提供する、統計的に妥当な不確実性推定を実現した。
• 運動学的データに対しては、関数形を仮定せずに、ユーザー定義の半径で平均円運動速度および非円運動速度（径方向および方位角方向）を推定した。
• 光度的ブランチでは、パラメトリックなプロファイルを仮定せず、ディスク、棒、バルジ成分の光度割合を推定でき、galfit よりも柔軟性が向上した。
• バージョン 1.2 では、動的配列割り当てによりデータサイズのソフトウェア的制限を撤廃し、重大なバグ（特に見かけのぼやけ補正モジュールのバグ）を修正した。
• スムージングペナルティは、主軸とほぼ一致する棒のフィットにおいて、正規化されていないフィットが物理的に非現実的な振動を生じるのを防ぐのに有効であった。