Flux-ratio anomalies in cusp quasars reveal dark matter beyond CDM

技術要約：カスプ・クエーサーにおけるフラックス比アノマリーは、CDMを超えたダークマターを明らかにする

問題提起
ダークマターの根本的な性質は、宇宙論における主要な問いである。冷たいダークマター（CDM）の枠組みは大規模構造を記述することには成功しているが、銀河スケールおよびサブ銀河スケールにおいては、内側ハローのプロファイル、矮小銀河の回転曲線、および小規模なクラスタリングに関する持続的な課題に直面している。強重力レンズ効果は、$\sim 10^6 M_\odot$ までのサブハローを探索するためのユニークなプローブを提供している。これは、CDM、自己相互作用ダークマター（SIDM）、およびファジーダークマター（FDM）が明確に異なる予測を行う領域である。具体的には、「カスプ関係（cusp relation）」、すなわち、カウスティック・カスプ付近の3つの合体像の符号付き倍率がほぼ相殺される（$R_{cusp} \approx 0$）という理論的期待は、小規模な摂動に対して非常に敏感である。しかし、これまでの解析は、サンプルサイズの小ささ、ステラー・マイクロレンズ効果による汚染、およびマクロなレンズ幾何学（例：ディスクや合体による高次マルチポール）による摂動とダークマター・サブストラクチャーの信号を系統的に区別できない点に制限されてきた。

手法
著者らは、観測データと広範なモンテカルロ・シミュレーションを組み合わせた包括的な統計フレームワークを提示し、ダークマターの微物理学を検証している。

1. 観測サンプル： 本研究では、マイクロレンズ効果の影響を受けないフラックス比を持つ、利用可能な最大規模の17個の四重像カスプ・クエーサーのセットを利用している。データは、GO-2046 JWSTプログラム（マイクロレンズ効果の影響を受けにくい温かい塵のフラックス比を提供する）から、電波観測、およびHSTのナローライン測定から抽出されている。サンプルは、カスプ・トリプレットのレンズ銀河の主軸に対する位置に基づき、「メジャー軸（major-axis）」および「マイナー軸（minor-axis）」のカスプ構成として幾何学的に分類されている。
2. シミュレーション・フレームワーク： 著者らは、PICSレンズ・フレームワークのJAX加速GPU実装を用いて、$\sim 10^8$ 個の模擬レンズ実現体を生成している。
   • マクロモデル： ベースライン・モデルには、外部シア（external shear）を伴う特異等温楕円体（SIE）が含まれ、ディスク、非対称性、または合体駆動の構造を考慮するために、高次マルチポール摂動（$m=3, 4$）へと拡張されている。
   • サブストラクチャー： 以下の3つのシナリオに対して、滑らかなポテンシャル上にサブストラクチャー摂動が重畳されている：
     • CDM： NFWプロファイルに従うサブハローおよび視線方向（LOS）のハロー。
     • SIDM： 低質量ハロー（$<10^8 M_\odot$）が熱的重力崩壊を起こし、高密度なコアを形成する速度依存型散乱モデル。
     • FDM： 波による密度揺らぎ（グラニュール）および中心部のソリトニック・コア。これは結合したサブハローを抑制する一方で、連続的な干渉パターンを導入する。
3. 統計解析： 著者らは、ソース面（光源面）の構成が観測されたカスプ幾何学（開き角 $\phi$）と一致するように、2段階のマルコフ連鎖モンテカルロ（MCMC）スキームを用いてサンプリングを行っている。$R_{cusp}$ の理論分布は、$\phi$ のビン（階級）にわたるカーネル密度推定（KDE）を通じて構築される。ベイズモデル比較を用い、各ダークマター・シナリオに対するデータの証拠を定量化するためにベイズ因子を算出している。

主な貢献

• 幾何学的識別： 本論文は、マルチポール自由度を持つグローバルにパラメータ化されたマクロモデルを用いても、フラックス比アノマリーが持続する特定のパラメータ空間領域（マイナー軸および狭いメジャー軸のカスプ・レンズで構成される領域）を特定している。この領域は、ダークマターの微物理学の信号を、マクロモデルの縮退から効果的に分離する。
• マイクロレンズフリーのデータセット： JWSTの温かい塵および電波データを利用することで、ステラー・マイクロレンズ効果を混乱変数として排除しており、これは従来のフラックス比アノマリー研究に対する決定的な改善である。
• マルチポールの頑健性： 著者らは、高次マルチポールがメジャー軸構成における一部のアノマリーを説明できる一方で、マイナー軸構成で見られる特定のアノマリーや、極端な外れ値であるJ1042+1641を再現することはできないことを示している。

結果

• CDMおよびSIDMとの不適合性： システムJ1042+1641は、マルチポール摂動を許容した場合でも、CDMおよびSIDMの両方のシナリオに対して$3\sigma$以上の不適合であることが判明した。
• FDMの優位性： FDMシナリオは、特にマイナー軸構成において明確なシグネチャを生み出す。ここでは、波による揺らぎが $R_{cusp}$ の分布をより高い値へとシフトさせ、観測された外れ値と一致する。
• ベイズ的証拠： 解析の結果、最も楽観的なCDMシナリオに対するFDMのベイズ因子は100を超え（CDMに対して $BF \sim 10^6$、SIDMに対して $BF \sim 10^4$）、FDMがCDMおよびSIDMよりも優れているという「非常に強い」証拠が得られた。これはジェフリーズ尺度における強力な証拠である。
• マクロモデルからの独立性： 本研究は、サブストラクチャーのない滑らかなレンズモデルが決定的に拒絶されることを確認した（$BF < 10^{-35}$）。さらに、FDMモデルにマルチポールの自由度を加えることはむしろ支持されない。これは、FDMが追加のマクロモデルの複雑さを必要とせず、データに対してより直接的かつ自然な説明を提供するものであることを示唆している。

意義
本論文は、現在のカスプ・レンズ・データに基づき、CDMやSIDMよりもFDMに対して非常に強い統計的証拠を提供することを主張している。本研究は、詳細なシステムごとのマクロモデリングに依存しない、フラックス比アノマリーを解釈するための堅牢で統計的に駆動されたフレームワークを確立した。著者らは、診断領域内にある現在の11個のレンズのサンプルサイズは小さいものの、本手法は将来のサーベイ（Euclid, Roman, CSST）がより大きなデータセットを用いてこれらの知見を決定づけるための明確な道筋を示すものであると述べている。本研究は、伝統的な結合サブハローが存在しない場合でも、波動的なダークマターが顕著で区別可能なレンズ・シグネチャを生み出し得ることを強調している。