A dense dark matter core of the subhalo in the strong lensing system JVAS B1938+666

技術要約：強重力レンズ系 JVAS B1938+666 におけるサブハロの密なダークマターコア

問題提起
ダークマターの根本的な性質は依然として未解決であり、特に $\Lambda$CDM パラダイムにおける小スケールの不一致に関する問題が残されている。冷たいダークマター（CDM）は大規模構造をうまく説明できる一方で、ナバロ・フレンク・ホワイト（NFW）密度プロファイルのような中心部の「カスプ（尖点）」を予測するが、これは矮小銀河で観測される平坦な「コア（核）」と矛盾する（コア・カスプ問題）。代替的な枠組み、具体的には自己相互作用ダークマター（SIDM）や、フェジー／波動ダークマター（FDM または $\psi$DM）は、コアを形成するメカニズムを提案している。しかし、これらのモデルを区別するには、従来の運動学的手法では高赤方偏移での観測が困難な低質量サブハロ（$10^7-10^9 M_\odot$）の密度プロファイルを精密に制約する必要がある。

手法
著者らは、強重力レンズ系 JVAS B1938+666（$z_{lens}=0.881, z_{source}=2.059$）におけるダークサブハロの質量分布の非パラメトリック再構成を提示している。

• データ・リダクション： 本研究では、Keck II 望遠鏡（NIRC2 装置、レーザーガイド星適応制御）によるアーカイブされた近赤外 $K_p$ バンド観測データを使用している。利用可能な 93 フレームから、66 枚の高品質フレームを選択し、サブピクセル単位のシフトを用いて整列させ、スタックした。中心のレンズ銀河の光は、2 つの Sersic プロファイルを用いてモデル化および減算され、低 SNR ピクセルはレンズ弧を孤立させるためにマスクされた。
• レンズモデリング： LensCharm ソフトウェアを用いた 2 ステップのモデリング手法を実装した。
  1. ステップ 1： 主要なレンズ銀河（擬似等温楕円体）と光源（パラメトリックな Sersic モデル）のパラメトリックモデルをデータに適合させた。
  2. ステップ 2： サブハロの質量分布をモデル化するために、非パラメトリックな収束摂動場（$\Delta\kappa_{sub}$）を導入した。この場は、特定の密度プロファイルを事前には仮定せずに、サブハロの質量分布を適応的かつモデルに依存しない形で再構成できるよう、積分ウィーナー過程によって制御されたガウス過程を利用している。
• プロファイル・フィッティング： 再構成されたサブハロ密度プロファイルに対し、競合する 3 つの理論モデルへのフィッティングを行った：標準的な NFW プロファイル（CDM）、シミュレーションから導出されたパラメトリックな SIDM プロファイル、およびソリトニック $\psi$DM プロファイルである。ベイズ推論によるマルコフ連鎖モンテカルロ法（MCMC）を用いて、対数尤度、$\Delta\chi^2$、およびベイズ因子を計算し、モデル間の比較を行った。

主な貢献

• 非パラメトリック再構成： 本研究は、検出段階において特定の関数形式（例：冪乗則）を仮定することによるバイアスを回避し、JVAS B138+666 におけるサブハロの密度プロファイルの堅牢でモデルに依存しない再構成を提供している。
• 定量的モデル比較： 高解像度の Keck データを用い、単なる傾き（$\gamma$）の制約を超えて、NFW、SIDM、$\psi$DM の間の直接的な統計的比較を行う。
• コアの特性評価： 再構成の結果、サブハロに特徴的なキロパーセック規模のコアが存在することが明らかになった。これは、$10^8 M_\odot$ の質量領域におけるダークマター物理学をテストするための稀な経験的データポイントを提供するものである。

結果

• 密度プロファイル： 再構成されたサブハロは、半径 $r_c \approx 0.5$ kpc、中心密度 $\rho_c \approx 2.5 \times 10^7 M_\odot \text{kpc}^{-3}$ の特徴的なコアを示す。
• モデルの選好性：
  • 標準的な NFW プロファイルはデータに対して適合が悪く（$\psi$DM に対して $\Delta\text{BIC} = -186.3$）、観測された平坦なコアとの間に重大な不一致を示している。
  • SIDM と $\psi$DM はともに NFW よりも有意に優れた適合を示す。
  • SIDM は、ベイズ因子 $BF = 14.44$（$\Delta\text{BIC} = 2.67$）をもって $\psi$DM よりも好まれる。
• パラメータ制約：
  • 最も好ましい $\psi$DM モデルについて、粒子質量は $m_\psi = 1.32^{+0.22}_{-0.31} \times 10^{-22}$ eV と制約される。
  • 再構成されたプロファイルは、局所銀河群の矮小銀河におけるコア・ハロー質量スケーリング関係と一致している。
• バリオン・フィードバック： 本研究では、有意な恒星質量への寄与は見られない（$M_* < 10^3 M_\odot$）。バリオン・フィードバック（超新星駆動のガス流出）は理論的にコアを作り出すことができるが、得られた光度限界からの低すぎる恒星対ハロー質量比（$M_*/M_{halo}$）は、フィードバックだけでは観測されたコアを説明するには不十分である可能性を示唆しており、SIDM や $\psi$DM のような素粒子物理学的な解を支持している。

意義と主張
本論文は、JVAS B1938+666 における高密度のダークマターコアが、この質量スケールのサブハロに対する標準的な NFW 予測に対する強力な反証となることを主張している。結果は、この観測された密度プロファイルを説明するために、SIDM や $\psi$DM といった代替的なダークマターの枠組みが必要であることを示唆している。著者らは、データは $\psi$DM よりも SIDM を緩やかに支持しているものの、この特定の系においては両者とも CDM より優れていると述べている。

本研究は、摂動体が（以前の HST による解析で $z \sim 1.4$ にある可能性が示唆された）視線方向のインターローパーではなく、レンズ面（$z \approx 0.881$）にあるサブハロであるという仮定に基づいているという制限事項を認めている。もし対象が視線方向のハローであれば、質量と集中度の制約は変化するが、平坦なコア自体は、強力なバリオン・フィードバックがない限り、標準的な CDM 予測に依然として挑戦するものとなる。著者らは、ダークマター候補を決定的に区別し、バリオン・フィードバック効果をさらに制約するためには、将来の高感度観測装置（ELT、JWST など）とより大規模なデータセットが必要であると結論づけている。