Validating Digital Twins of the Local Universe with the Thermal Sunyaev-Zel'dovich Signal

技術要約：熱的スニヤエフ・ゼルドビッチ信号を用いた局所宇宙のデジタルツインの検証

問題提起
熱的スニヤエフ・ゼルドビッチ（tSZ）効果は、銀河団や宇宙網における電離ガスの熱圧力を探るための主要なプローブである。標準的な宇宙論的解析では、tSZデータは統計的な制約（例：パワースペクトル、クラスター数）のために利用されることが多いが、これらの手法はオブジェクトごとの比較に必要な空間位相情報に欠けることが多い。一方、局所宇宙を再現することを目的とした制約付きシミュレーションは、速度場との比較を超えた、頑健で独立した検証を必要とする。本論文は、これらの「デジタルツイン」（データに基づき制約された事後分布シミュレーション）が、局所的な宇宙構造の位置、質量、および熱圧力分布の再構成における忠実度を評価するために、tSZ信号を用いて検証する必要性に取り組んでいる。

手法
著者らは、最新のベイズ的銀河起源再構成（BORG）再構成に基づくデジタルツインの新しいスイートであるCSiBORG-Manticore (CBM)を導入している。これらのシミュレーションは、2M++銀河カタログから推論された局所密度場（$\sim155 h^{-1}$ Mpc以内）の事後分布をサンプリングしており、高解像度のダークマター・ズームインによって拡張されている。著者らは、CBMを前世代のCSiBORG (CB2)、およびCLUES初期条件に基づく流体シミュレーションSLOWと比較している。

検証フレームワークは、Planckコンプトン$y$マップ（McCarthy and Hill 2024）およびeROSITA X線クラスターカタログ（弱重力レンズによる較正済み）を用いた3つの補完的な指標を用いている：

1. 角度的位置精度 ($p_{tSZ}$): 特定の観測されたクラスターに対し、シミュレーションされたハロー位置におけるtSZ信号をアパーチャフォトメトリで測定する。観測されたtSZホットスポットと信号を比較することで、$p$値（$p_{tSZ}$）を算出する。低い$p_{tSZ}$値は、シミュレートされたハローが観測されたtSZホットスポットと一致していることを示す。
2. スタックされた半径方向プロファイル: 質量依存性を評価するため、著者らは質量によってランク付けされたハロー（上位10、50、100）に対して、全リアリゼーションにわたる1次元の半径方向tSZプロファイルをスタックする。これにより、シミュレーションが質量増加に伴うtSZ信号の増加という期待される傾向を再現しているかどうかをテストする。
3. 質量スケーリング関係: 質量の不確実性を周辺化するベイズ線形回帰を用いて、BORGハロー質量と観測量との間のスケーリング関係を推論する：
   • $Y_{500c}^{tSZ}$ vs. ハロー質量（自己相似的な傾き $m \approx 5/3$ のテスト）。
   • BORG質量 vs. Planck tSZ質量。
   • BORG質量 vs. 弱重力レンズで較正されたeROSITA X線質量。
     マッチングは、厳格な基準（$\Delta\theta < 60'$, $\Delta cz < 300$ km s$^{-1}$）および、観測されたクラスターとシミュレートされたハローの関連付けを行うための局所宇宙モデル（LUM）のマッチング手順を用いて行われる。

主な貢献

• CSiBORG-Manticore (CBM) の導入: 最新のBORG再構成（Manticore-Local初期条件）を表す50個の事後サンプルからなる新しいデジタルツイン・スイート。これは、前世代と比較して改善された初期条件と銀河バイアスモデリングを提供する。
• 新しい検証フレームワークの開発: 速度場との比較を超え、宇宙網の熱圧力を直接テストすることによる、制約付きシミュレーションのためのスコアリングシステムの開発。
• 質量・観測量関係の較正: ダークマター分布の既知の分布を利用することで、外部の質量プロキシに頼ることなく、高精度なデジタルツインがいかに質量–観測量関係（特に$Y$–Mass関係）を較正できるかを実証する。

結果

• 位置精度: 選択された局所クラスターの大部分（例：Coma, Shapley A3558, Perseus）は、CBMにおいて高い忠実度で再現されており（$p_{tSZ} < 0.05$）、CBMはSLOWシミュレーションを大幅に上回っている。SLOWは2つのクラスター（PerseusとShapley A3558）のみが $p_{tSZ} < 0.05$ を示し、大きな角度オフセット（Comaで$\sim15^\circ$）を持つ。対照的に、CBMはほとんどのターゲットに対して高い有意性を達成し、オフセットは通常$<1^\circ$である。
• クラスターごとの再構成: CB2には存在しなかったクールコア・クラスターであるPerseus (A426) は、CBMにおいてロバストに再現されており（$p_{tSZ} \approx 2.9 \times 10^{-4}$）、これは基礎となるBORG推論の改善を強調している。逆に、Hercules (A2151) はCBMにおいて依然として不十分な再現となっている。
• 質量スケーリング:
  • $Y$–Mass 関係: CBMは、$Y_{500c}^{tSZ}$–Mass関係に対して $m = 1.79 \pm 0.28$ という傾きを与えており、これは自己相似的な期待値（$5/3$）と一致し、CB2（$m = 1.43 \pm 0.31$）よりも理論値に近い。また、CBMはCB2（$\sigma_{int} = 0.09 \pm 0.07$ dex）と比較して、より低い固有の散らばりを示す。
  • 質量較正: BORG質量をPlanck tSZ質量と比較すると、CBMは系統的なオフセット（$c \approx -0.20$ dex）を示しており、これはBORG質量が$\sim1.6$倍高いことを意味する。これは、Planck tSZ質量における既知の$\sim30\%$の静水圧質量バイアスと一致している。
  • 弱重力レンズによる検証: 決定的なことに、弱重力レンズで較正されたeROSITA X線質量と比較した場合、CBMは一対一の関係（$m = 0.84 \pm 0.19, c = 0.02 \pm 0.06, 0.5\sigma$ のテンション）と極めて良好に一致している。これに対し、CB2は大きく逸脱している（$2.1\sigma$ のテンション）。このことは、CBMのハロー質量が、前世代の再構成よりも真のクラスター質量をより正確に反映していることを示唆している。

意義
本論文は、tSZ効果が、速度場テストを補完する、制約付きシミュレーションのための強力かつ独立した検証観測量を提供すると主張している。結果は、最新のBORGベースのデジタルツイン（CBM）が、前世代や代替の速度制約シミュレーション（SLOW）よりも大幅に高い位置および質量の忠実度を達成していることを示している。

著者らは、これらの高精度なデジタルツインが、tSZデータから追加の物理情報を抽出するための実用的な経路を提供すると考えている。tSZ信号を事後分布の広域構造のリアリゼーションに固定することで、研究者は、統計的解析において通常は摂動パラメータとして吸収される投影効果や視線方向の寄与を明示的にモデル化できる。この研究は、制約付きシミュレーションに内在する大規模構造の情報を、CMB二次異方性の研究へと統合する道を開くものであり、将来的にCMB一次および二次信号と大規模構造を自己整合的なフィールドレベルの尤度フレームワーク内で結合解析することを可能にする。