[論文レビュー] A fast direct method of mass reconstruction for gravitational lenses

M. Lombardi, G. Bertin|arXiv (Cornell University)|Jun 7, 1999

Astrophysical Phenomena and Observations参考文献 1被引用数 3

ひとこと要約

本論文は、変分原理を効率的に解くことで、重力レンズの質量分布を再構成する高速な直接法を提示している。従来の手法と比較して計算時間を100〜1000倍短縮し、N = 400までのグリッドに対して標準的なワークステーションでも1秒未塔の再構成が可能になった。これにより、大規模な統計的・宇宙論的シミュレーションが、実用的な速度で初めて可能となった。

ABSTRACT

Statistical analyses of observed galaxy distortions are often used to reconstruct the mass distribution of an intervening cluster responsible for gravitational lensing. In current projects, distortions of thousands of source galaxies have to be handled efficiently; much larger data bases and more massive investigations are envisaged for new major observational initiatives. In this article we present an efficient mass reconstruction procedure, a direct method that solves a variational principle noted in an earlier paper, which, for rectangular fields, turns out to reduce the relevant execution time by a factor from 100 to 1000 with respect to the fastest methods currently used, so that for grid numbers N = 400 the required CPU time on a good workstation can be kept within the order of 1 second. The acquired speed also opens the way to some long-term projects based on simulated observations (addressing statistical or cosmological questions) that would be, at present, practically not viable for intrinsically slow reconstruction methods.

研究の動機と目的

• 重力レンズの銀河のせん断データから質量分布を再構成する際の計算ボトルネックを解消すること。
• 近い将来の観測プロジェクトで求められる、数千の源銀河を含む大規模データセットの効率的処理を可能にすること。
• 質量再構成の精度を維持しつつ、計算時間を著しく短縮する手法を開発すること。
• かつて再構成が遅すぎて現実的でなかった長期間のレンズ観測シミュレーションを、実用的に可能にする。

提案手法

• 文献に先行して指摘された変分原理を解くことで、質量再構成を制約付き最小化問題に再定式化する。
• 長方形のフィールドでは、基礎となる線形方程式系の解法を高速化するために、高速フーリエ変換（FFT）に基づくアルゴリズムを活用する。
• 反復解法を回避するため、変分原理から導かれた閉形式の式を直接計算することで解を求める。
• 定期的なグリッドを想定して最適化されており、FFTを用いた効率的な行列・ベクトル演算が可能である。
• 標準的なワークステーションに実装され、N = 400グリッドに対して1秒未塔の性能を達成している。

実験結果

リサーチクエスチョン

• RQ1 現在の反復的手法よりも高速に質量分布を再構成できる直接法を開発することは可能か？
• RQ2 変分原理の直接解法を用いることで、質量再構成の最大速度向上はどの程度達成可能か？
• RQ3 数千の源銀河を含む大規模データセットに対しても、この手法はリアルタイム性能を維持できるか？
• RQ4 この再構成速度のおかげで、長期間のレンズ観測シミュレーションが今や現実的に行えるようになるか？

主な発見

• 本手法は、質量再構成のための最も高速な既存手法と比較して、CPU時間を100〜1000倍短縮した。
• N = 400までのグリッドサイズでは、優れたワークステーション上でも必要な計算時間が約1秒に抑えられた。
• 変分原理を介した直接解法により、統計的および宇宙論的調査に向けたシミュレートされたレンズ観測の実用的利用が可能になった。
• FFTに基づく高速化のおかげで、本手法は大規模データ処理に適している。