Detection of a dark matter subhalo in the strongly lensed system PJ011646

技術的サマリー：強重力レンズ系 PJ011646 におけるダークマターサブハロの検出

問題提起
ダークマターの性質は、宇宙論における未解決の根本的な課題のままである。冷たいダークマター（CDM）モデルは、非常に低い質量まで及ぶ自己束縛ハロの階層的な形成を予測する一方、温かいダークマター（WDM）モデルは、特徴的なカットオフスケール以下の低質量ハロの抑制を予測する。これらのモデルを区別するには、特に予測が分かれる質量領域（$M_{200} \sim 10^7 - 10^9 M_\odot$）におけるサブハロ集団の統計的測定が必要である。強重力レンズを介してこれらのサブ構造を検出する際の主要な課題は、真のサブハロ信号とレンズ銀河の滑らかな質量分布におけるモデル化されていない角度的複雑さとの間の「縮退」である。単純な楕円形のパワーロー（PL）モデルは、早期型銀河（ETG）の偏差を捉えきれず、誤ったサブハロ検出を引き起こしたり、実際のものを隠蔽したりすることが多い。さらに、過去の検出は光学望遠鏡（例：HST）の角分解能と、レンズ銀河の光による汚染に制限されてきた。

手法
著者らは、アタカマ大型ミリ波サブミリ波干渉計（ALMA）のバンド 7 における高解像度（約 0.1 秒角）のダスト連続放射観測を用いて、強重力レンズ系 PJ011646（レンズ $z=0.555$、光源 $z=2.125$）の解析を提示する。解析には、オープンソースソフトウェアスイートである PyAutoLens、PyAutoFit、および PyAutoGalaxy が用いられた。

手法は主に 3 つの段階で進行する：

1. マクロモデリング：レンズ質量分布は、外部せん断を伴う楕円形パワーロー（EPL）プロファイルを用いてモデル化される。単純な楕円形モデルの限界に対処するため、著者らは質量密度の角度的偏差を捉えるために、高次のフーリエ多重極（$m=3$ および $m=4$、さらに $m=1$）を組み込んだ。これらの多重極の収束は、周期的位相に関連する縮退を避けるためにパラメータ化された。
2. ソース再構成：ソースの表面輝度は、ボロノイメッシュ画素化と自然近傍補間を用いたピクセル化ベイズ逆変換法（Suyu et al. 2006）によって再構成される。これは完全な干渉計データを活用するため、直接 $uv$ 平面で行われる。
3. サブハロ探索：画像平面でグリッドベースの探索を行い、擾乱を検出する。各グリッドセルに球形のナバーロ・フレンク・ホワイト（NFW）プロファイルを加擾乱体として追加し、最良適合マクロモデルに対するベイズ対数尤度差（$\Delta \ln Z$）を計算する。堅牢な候補を識別するために、$\Delta \ln Z > 10$ という閾値を採用する。また、標準的な CDM プロファイルからの偏差を検証するため、自由な内部傾き（$\gamma$）を持つ一般化 NFW（gNFW）プロファイルもテストされる。

主要な結果

• マクロモデルの改善：3 次および 4 次多重極（$m=3, 4$）の組み込みは、データへの適合を著しく改善し、純粋な PL モデルに対して $\Delta \ln Z = 52.1$ のベイズ証拠の増加をもたらした。これらの多重極の振幅は収束の約 1.5% であり、近傍の ETG の観測と一致する。1 次多重極（$m=1$）を含むモデルは、わずかな追加の改善（$\Delta \ln Z = 6.5$）を提供し、物理的な伴天体の可能性を示唆するが、主要な信号を駆動するものではない。
• サブハロ検出：最良適合の PL+MP マクロモデルを用いて、画像平面で $\Delta \ln Z > 10$ を示す 2 つの領域をグリッド探索が特定した。しかし、フォローアップ解析により、これら 2 つの領域がソース平面では同じ位置にマッピングされ、縮退を示していることが明らかになった。2 つのサブハロモデルは、単一サブハロモデルよりも証拠を改善せず、2 番目の候補の事後分布は非物理的な濃度（$c_{200} > 1000$）をもたらす。したがって、著者らは、北側の弧の右側に位置する 1 つのサブハロのみが物理的に実在すると結論づけた。
• サブハロの性質：検出されたサブハロの質量は $M_{200} = 2.78^{+0.43}_{-0.66} \times 10^{10} M_\odot$、濃度は $c_{200} = 30^{+5}_{-7}$ である。この濃度は、親ハロ内の潮汐剥離された NFW サブハロに対する期待値と一致する。検出の有意性は、PL+MP マクロモデルに対して約 5.8$\sigma$ である。
• 堅牢な質量制限：NFW スケール半径と正規化の間の縮退により、総質量 $M_{200}$ は厳密に制限されていない。代わりに、著者らは投影された閉込め質量の分数的不確かさが最小化される「堅牢半径」（$R_{sub} \approx 2$ kpc）を定義する。この半径内での閉込め質量は $M_{sub} = 3.57^{+0.16}_{-0.14} \times 10^9 M_\odot$ である。
• 感度限界：サブハロが検出されなかった（$\Delta \ln Z < 10$）グリッドセルにおいて、著者らは検出可能なサブハロ質量の上限を導出した。最も感度の高い領域（明るいレンズ化された弧に沿った領域）において、Ludlow et al. (2016) の質量 - 濃度関係を仮定した場合、3$\sigma$ 有意性における最小検出可能質量は $M_{200} \approx 8 \times 10^8 M_\odot$ である。

意義と主張
本論文は、サブ秒角分解能を持つ ALMA 連続放射イメージングが、CDM モデルと WDM モデルが異なる存在量を予測する質量領域（$M_{200} \sim 10^9 M_\odot$）におけるダークマターサブ構造を探査する能力を実証していると主張している。多重極を介してレンズ銀河の角度的複雑さを明示的にモデル化することで、著者らは歴史的にサブハロ探索を悩ませてきた系統的不確かさを軽減した。CDM 予測と一致する性質を持つサブハロの検出と、非検出領域における厳格な上限の導出は、小規模構造形成の堅牢なテストを提供する。著者らは、角分解能が HST と同等である一方で、この系の高い増幅率と ALMA データにおける前景光の汚染の欠如により、以前光学データで達成可能だったよりも低い質量を探査できることに言及している。本研究は、ALMA で観測される強重力レンズ化された塵埃星形成銀河（DSFGs）の増加するサンプルを用いた、体系的で統計的に意味のあるダークマターサブ構造の探索への一歩として位置づけられている。