A Consistent Implementation of Cluster Strong Lensing in Cosmological Simulation Light Cones

技術要約：宇宙論的シミュレーションのライトコーンにおける、一貫した銀河団強重力レンズの実装

問題提起
銀河団の強重力レンズ現象は精密宇宙論の礎石であるが、堅牢な理論的予測は歴史的に、高品質な観測結果に比べて遅れをとってきた。このギャップは、主に2つの問題に起因している。一つは、標準的な立方体状のシミュレーションボックスと、ライトコーン（光錐）に求められる細長い幾何学的形状との間の不一致であり、もう一つは、巨大なクラスターを収容するために必要な広大な体積と、背景にあるソース銀河やサブ構造を解像するために必要な高解像度との間の根本的な緊張関係である。

現在の予測手法は、シミュレーションから個々のレンズを抽出し、それらを解析的または観測されたソース集団と組み合わせるハイブリッドなアプローチに依存することが多い。これらの手法は、赤方偏移にわたる相関関係を保持できないことが頻繁にあり、あるいは、相関のある、および相関のない視線方向（LoS）の構造を省略してしまう。その結果、主要なレンズ、ソース、および介在するすべての物質が、標準的な冷たい暗黒物質（CDM）パラダイムに忠実な形で同一の宇宙論的体積から描かれている、シミュレーション出力から一貫して生成された強重力レンズ像の生成手法が欠如している。

手法
著者らは、IllustrisTNG宇宙論的流体シミュレーションの粒子データから直接、強重力レンズ像を生成するための、完全にシミュレーションに基づいた手順を提示している。この手法は、主に以下の3つのコンポーネ Zeit 構成要素からなる。

1. 体積の再マッピング（Volume Remapping）: 幾何学的な制約に対処するため、著者らは BoxRemap コードを用いて、立方体のシミュレーション体積を直方体に変換するボックス再マッピング技術を採用している。この幾何学的形状は、2つの短い軸（空の面を表す）と1つの長い軸（視線方向を表す）を備えている。この変換は、総体積と局所的な構造を保持しつつ、単一のライトコーン内に重複した構造を導入することなく（フィールド・オブ・ビューがボックスサイズに対して適切に選択されている場合）、再マッピングを適用する。この再マッピングは、シミュレーションのすべてのスナップショットに対して一貫して適用される。
2. 平面の抽出と平滑化（Plane Extraction and Smoothing）: 再マッピングされたボックスから、質量および光の情報を抽出し、シミュレーションのスナップショットに対応する一連のレンズ面およびソース面を作成する。
   • 質量面（Mass Planes）: 粒子をスナップショットの赤方偏移に基づく動径方向のシェルに分割する。質量マップは、近傍粒子（32個の近傍、Wendland-C^2 カーネル）を用いた平滑化によって、ダークマター、恒星、およびガスの粒子を用いて生成される。希薄な領域における拡散質量を考慮し、ショットノイズを避けるために、シェルの総質量を一様なシートとして加え、宇宙の平均物質密度を差し引くことで質量コントラストマップを作成する。
   • ソース面（Source Planes）: 恒星粒子を、特定の観測フィルター（例：JWST f200w）に対するフラックスマップへと平滑化する。フラックスは、恒星の年齢、初期質量、および金属量に基づいて推定され、ソース面の赤方偏移に合わせて計算され、見かけの等級に変換される。粒子数が不十分なサブハロは、完全性を確保するために解析的なセルシック（Sérsic）プロファイルで置き換えられる。
3. マルチプレーン・レイトレーシング（Multi-Plane Ray Tracing）: 抽出された平面は、マルチプレーン・レイトレーシング・アルゴリズム（GLAMER コードを使用）によって用いられる。赤方偏移 $z_s$ にあるソース面からの光は、すべてのレンズ面 $z_l < z_s$ を通って逆方向に追跡される。得られたレンズ像はスタックされ、最終的な理論的予測が生成される。このパイプラインは、ソースの強制的な注入や解析的なプロファイルを必要とせず、宇宙論および銀河形成モデルが事前分布として機能する。

主な貢献と結果
著者らは TNG300-1 シミュレーションを用い、レンズ赤方偏移 $z_G \approx 0.38$ において質量的な銀河群（$M_{200m} \ge 5 \times 10^{14} M_\odot$）を中心とし、$z_{max} \approx 4.18$ まで延びるライトコーンを生成した。

• 画像生成: パイプラインは、シミュレーションデータから直接、解析的なソース注入なしに、強重力レンズの特徴（多重像の構成、弧状像、および臨界曲線を含む）を正常に生成した。生成された画像は、クラスターの合体やサブ構造によって駆動される、接線方向の弧（タンジェンシャル・アーク）や複雑な臨界曲線の形態といった典型的な特徴を示している。
• 視線方向（LoS）構造の影響: 著者らは、フル・ライトコーン・レイトレーシングと、主要なレンズ面のみを通るレイトレーシングを比較することで、相関のない LoS 構造の影響を定量化した。
  • 画像の配置: 相関のない LoS 構造は、レンズ像の相対的な位置を数秒角（極端なケースでは最大 $\sim 10''$ まで）シフトさせる可能性があり、これは典型的な強重力レンズの質量モデリングの残差（$\lesssim 1''$）を大幅に上回る。
  • 臨界曲線の形態: LoS 物質の包含は、クラスターの主要な臨界曲線の面積に $\sim 6\%$ の散らつきをもたらす。いくつかのケースでは、主要な臨界曲線の面積が30%増加、あるいは50%減少した。
  • 総臨界面積: クラスターのポテンシャルの最小値から $100''$ 以内において、$z_s = 4$ のときの総臨界面積は、LoS 構造が含まれることで $16^{+20}_{-14}\%$ 変化する。この増大は、LoS オブジェクトによる新しい臨界曲線の生成と、既存の曲線の修正の両方によって生じる。

意義と限界
本論文は、宇宙論的流体シミュレーションのライトコーンから自己整合的に生成された、クラスター強重力レンズ像の最初の実証であると主張している。このアプローチは、モック観測の膨大なカタログを生成するための、より自動化された「ハンズオフ」な手法を提供し、これらはレンズモデリングの系統誤差、質量推論のバイアス、および銀河―銀河強重力レンズ（GGSL）の断面積の緊張を研究するために使用できる。

著者らは以下の限界を指摘している：

• 大規模構造のパワー: この手法は、初期のシミュレーションボックスに欠けている大規模なパワーモードを導入しないため、大規模構造のパワースペクトルのサンプリングにバイアスを与える可能性がある。
• 解像度: この手法はシミュレーションの質量解像度に制限されており、巨大なクラスターと低質量のサブ構造／ソースを同時に解像することに苦慮している。
• 拡散質量: 拡散質量は、その特定の構造を保持するのではなく一様なシートとして近似されているが、これはショットノイズを避けるために必要である。
• 計算コスト: 単一のライトコーンの生成には、各スナップショットにおけるシミュレーションボックス内のすべての粒子を処理する必要があり、計算コストが非常に高く、膨大なデータストレージを必要とする。
• 外部シア（External Shear）: 現在の実装には、ライトコーンの外側からの大規模構造による外部シアは含まれていないが、期待される効果（$\gamma_{ext} \sim 0.02$）は推定されている。

著者らは、フル・ライトコーンの情報を取り入れることは理論的予測にとって極めて重要であり、観測を解釈する際には、相関のない構造が画像構成や臨界曲線の特性を著しく変化させ得るため、LoS 構造を慎重に考慮することが必要であると結論付けている。