A Novel Lensed Point Source Modeling Pipeline using GIGA-Lens with Application to SN Zwicky and SN iPTF16geu

技術要約：GIGA-Lensを用いた新しいレンズ点光源モデリング・パイプライン

問題提起
本論文は、初期宇宙と後期宇宙のハッブル定数（$H_0$）の測定値の間にある顕著な不一致（約5$\sigma$）である「ハッブル・テンション」に取り組んでいる。強重力レンズによる時間変動イベント（クエーサーや超新星）は、タイムディレイ（時間遅延）測定を通じてこのテンションを解決する経路を提供するが、レンズ超新星（SNe）は稀少である。点光源に対する既存のモデリング手法には、以下の特有の課題が存在する：

1. ピクセルベース手法の適用不能性： 従来のレンズモデリング（拡張光源用のGIGA-Lensなど）は、観測画像とモデル画像の間のピクセルごとの差を最小化することに依存している。しかし、点光源の場合、最適化の開始時にフォワードモデル化された画像位置が観測位置と重ならないことがあり、その結果、勾配がゼロになり最適化が機能しなくなる。
2. 損失関数におけるバイアス： 画像位置のみを使用する従来の手法は、距離指標を最小化するものの物理的制約に反する非物理的なパラメータ領域（例：質量勾配 $\gamma \sim 1$）へとサンプラーを誘導してしまう。
3. データの希少性： 多くのレンズ超新星は、拡張したホスト銀河の光が検出できないため、拡張光源モデリング技術を利用することができない。
4. 系統的な不確実性： マスシート縮退（MSD）や、精密な$H_0$制約を得るためのフルフォワードモデリング（全パラメータの同時推論）の必要性は、依然として重要な障壁である。

手法
著者らは、4枚のA100 GPU上で動作するJAXプラットフォームを活用した、GIGA-Lensフレームワークに基づく新しいモデリング・パイプラインを導入している。核心となる革新は、拡張したホスト銀河のデータを必要とせず、重心位置、観測されたフラックス（または倍率）、およびタイムディレイのみを利用する、点光源専用に設計された微分可能な損失関数である。

尤度関数 $\log \mathcal{L}$ は、経験的に決定された相対的な重み（$\omega_D, \omega_F, \omega_{TD}$）を用いて3つの項を組み合わせている：

1. ソースプレーンのコンパクト性（$\mathcal{L}_D$）： ピクセル比較の代わりに、パイプラインは各観測された画像位置を、レンズモデルパラメータを用いてソースプレーンへと「デレンズ（逆レンズ効果による復元）」する。損失関数は、これらのデレンズされた位置とその平均との間の二乗距離を最小化し、すべての画像が単一の点光源に由来することを強制する。
2. フラックス／倍率（$\mathcal{L}_F$）： 質量勾配を制約し縮退を打破するために、パイプラインは観測された逆二乗倍率をモデル化する。ここでの重要な技術的貢献は、非微分的な倍率関数（クリティカルカーブ付近で特異点を持つもの）を、滑らかで微分可能な成分へと分解することである。これにより、倍率の計算に**自動微分（AD）**を用いることが可能となり、有限差分近似を回避して、計算速度を20〜40%向上させている。
3. タイムディレイ（$\mathcal{L}_{TD}$）：パイプラインは、フェルマポテンシャル（または$H_0$を自由パラメータとして扱う直接的なタイムディレイ）の観測値と予測値の差を最小化することで、タイムディレイを組み込んでいる。

技術的実装の詳細：

• 微分可能性とビヨクター（Bijectors）： パラメータ境界（例：質量勾配 $\gamma > 1.5$ の制限）を、数値的不安定性を導入せずに扱うために、著者らは制約のないパラメータ空間を制約された物理的空間に写像するビヨクターを使用している。
• 重み付けスキーム： 損失項間の相対的な重みは、観測の不確かさのみによって設定することはできないことを著者らは示している。彼らは、ハミルトニアン・モンテカルロ（HMC）サンプリング中に事後分布の各成分が同等に寄与するように（非物理的な解へのバイアスを防ぐために）、重みを増大させる（例：フラックスに対して）という経験的な重み付けスキームを採用している。
• サンプリング： パイプラインは、3段階の最適化を使用する：最大事後確率（MAP）のためのマルチスタート勾配降下法、ガウス近似のための確率的変分推論（SVI）、および完全な事後分布サンプリングのためのHMCである。

主な結果

• シミュレーション系： パイプラインは、5つの典型的なレンズ構成（クロス、ロングカスプ、ショートカスプ、フォールド、ダブル）でテストされた。

  • 画像位置、フラックス、およびタイムディレイの回収に成功した。
  • 点光源データのみを用いて質量勾配 $\gamma$ を制約した。
  • $H_0$を自由パラメータとしたシミュレーション系において、フルフォワードモデリングを用いた単一のシステムから、$H_0$に対して$\sim 3.6\%$（$\sim 2.5$ km s$^{-1}$ Mpc$^{-1}$）の統計的不確定度を達成した。
  • ほとんどのシステムで収束指標（$\hat{R} < 1.01$ および $\text{ESS} \gtrsim 10^3$）を達成し、4枚のA100 GPUを用いたモデリング時間は数分から約23分であった。

• SN iPTF16geuへの適用：

  • パイプラインは、拡張したホスト銀河を無視し、4つのレンズ化された超新星画像（位置と倍率）のみを使用してこのシステムをモデル化した。
  • 推定された質量勾配（$\gamma \approx 1.92$）およびその他のパラメータは、速度分散やホスト銀河のデータを利用した先行研究（例：Mörtsell et al. 2020）と$1\sigma$以内で一致した。
  • 結果は、滑らかなレンズモデル単独ではフラックスのアノマリーを説明できないことを裏付け、マイクロレンズ効果の仮説を支持した。

• SN Zwickyへの適用：

  • パイプラインは、先行研究（Goobar et al. 2023; Pierel et al. 2023）と比較して代替的な最良適合モデルを生成した。
  • 従来のモデルは、低い楕円率を示唆し、クリティカルカーブに対して対称な画像を示していたが、これは観測された明るさの非対称性（画像AとCがBとDよりも著しく明るいこと）を説明できなかった。
  • 本モデルは、より高い楕円率（$q \approx 0.38$）と質量勾配 $\gamma \approx 2.11$ を推定した。この幾何学的配置により、画像AとCはクリティカルカーブ上に直接配置され（これが明るさを説明する）、一方でBとDはそこから離れた位置になる。
  • モデルは、観測された位置、フラックス、および改訂されたタイムディレイの順序（画像Aが最初に到着する）を正確に予測しており、これは最近の精緻化された測定値（Larison et al. 2024）と一致している。

意義と主張
論文は以下のことを主張している：

1. 点光源データの十分性： 拡張したホスト銀河の情報がなくても、レンズ化された点光源データ（位置、フラックス、タイムディレイ）のみを用いて、複雑なレンズ質量分布（質量密度勾配 $\gamma$ を含む）を推論することが可能である。
2. フルフォワードモデリング： すべてのレンズパラメータと $H_0$ を同時に推論することは可能であり、システムが十分に大きなアインシュタイン半径（$\theta_E \gtrsim 1''$）を持つ場合、単一のシステムから精密な制約（3.6%の不確定度）が得られる。
3. 特異性への堅牢性： 倍率の分解に基づくアプローチにより、クリティカルカーブ付近の画像に対するロバストなモデリングが可能となり、これは従来の手法が失敗したり手動介入を必要とした領域である。
4. パラメータ空間の探索： SN Zwickyのケーススタディは、パラメータ空間を完全に探索することの重要性を強調している。著者らのモデルは、画像対称性と明るさ非対称性の間の矛盾を解決する物理的に妥当な構成を明らかにしており、これは従来のモデルが見落としていたものである。
5. 効率性： JAXとGPU加速の使用により、迅速かつ統計的に厳密なモデリング（$\hat{R} < 1.01$ という厳格な収束基準の達成）が可能となり、発見が増加しているレンズ超新星に対する実用的なツールとなっている。

著者らは、現在の焦点はSN Iaにあるが、このフレームワークはフラックス比を通じてクエーサーや非Ia型超新星にも容易に拡張可能であり、ホスト銀河のモデリングと組み合わせることで、さらに強力な宇宙論的制約を得ることができると結論付けている。