Accretion history dependence of the halo depletion radius

以下は、Jiale Zhou と Jiaxin Han による論文「Accretion history dependence of the halo depletion radius（ハロー枯渇半径の降着履歴依存性）」の技術的な要約です。

1. 問題提起 (Problem)

ダークマターハローの物理的な外縁（境界）を定義することは、宇宙構造形成や銀河進化の理解において重要な課題ですが、未解決の問題です。

• 従来の定義の限界: 一般的に用いられる「球対称崩壊モデル」に基づく過密度半径（例：$R_{200}$ や $R_{vir}$）は、境界の定義が恣意的であり、擬似進化（pseudo-evolution）の影響を受けやすく、動的に関連する物質が定義された半径を超えて存在するなどの問題があります。
• 新しい境界の必要性: 動的な状態や成長をより直接的に反映する境界として、「スプラッシュバック半径（$R_{sp}$）」が提案されましたが、これとは別に、密度プロファイルの時間進化から定義される「枯渇半径（Depletion Radius, $R_{id}$）」が注目されています。
• 未解明の点: 枯渇半径がハローの最近の質量降着率（MAR）、ハロー質量、赤方偏移、およびハローの二次的な性質（形状、スピンなど）にどのように依存するか、またスプラッシュバック半径との関係性が体系的に解明されていません。特に、$R_{id}/R_{vir} \simeq 2$ というほぼ一定の比率が観測・シミュレーションで報告されていますが、その物理的メカニズムは不明確でした。

2. 手法 (Methodology)

• シミュレーション: 「CosmicGrowth」スイートに属する大規模な宇宙論的 N 体シミュレーション（3072^3 個のダークマター粒子、ボックスサイズ 600 $h^{-1}$ Mpc）を使用しました。
• ハロー選定: 質量 $10^{12} \sim 10^{14.5} h^{-1} M_{\odot}$ のハローを、赤方偏移 $z=0$ から $z=3$ までの 9 つの出力で選定しました。環境汚染を避けるため、より大きなハローの近傍にあるものは除外し、中心サブハローの安定した進化履歴を確保する選定基準を適用しました。
• 境界の定義と測定:
  • 枯渇半径 ($R_{id}$): 連続の方程式に基づき、局所的な密度変化率 $\partial \rho / \partial t = 0$ となる点、すなわち質量流率（MFR: Mass Flow Rate）プロファイルの極小値として定義しました。
  • スプラッシュバック半径 ($R_{sp}$): 密度プロファイルの対数勾配の極小値として定義しました。
  • 質量降着率 ($\Gamma_{dyn}$): 現在のダイナミカルタイムスケール（$t_{dyn} = 2R_{vir}/V_{vir}$）における質量の対数成長率として定義しました。
• 解析手法: 異なる $\Gamma_{dyn}$、ハロー質量、赤方偏移、およびハローの二次的性質（形状パラメータ、スピン、濃度、形成時間）でハローを層別化し、スタックされた MFR プロファイルや密度プロファイルから $R_{id}$ と $R_{sp}$ を測定・比較しました。また、最近の主要合併（質量比 1:9 以上）を経験したハローを除外したサンプルを用いて、合併の影響を評価しました。

3. 主要な貢献と結果 (Key Contributions & Results)

A. 降着率と赤方偏移への依存性

• 降着率依存性: $R_{id}/R_{vir}$ はハロー質量にはほとんど依存しませんが、最近の質量降着率 $\Gamma_{dyn}$ に強く依存します。$\Gamma_{dyn}$ が高い（急速に成長している）ハローほど $R_{id}/R_{vir}$ は小さくなります（よりコンパクトな成長層）。
• スプラッシュバック半径との相関: $R_{id}$ と $R_{sp}$ は強く相関しており、ほぼ線形関係 ($R_{id}/R_{vir} \approx -0.50 + 1.72 R_{sp}/R_{vir}$) にあります。両者は同じ物理的プロセス（最初のカオスティック）を反映していますが、$R_{id}$ は密度と速度の両方を含む MFR で定義されるため、より大きな半径を示します。
• 赤方偏移依存性の違い: $R_{sp}/R_{vir}$ は赤方偏移に対して比較的弱い依存性を示しますが、$R_{id}/R_{vir}$ は明確な赤方偏移依存性を示します。低赤方偏移では $\Gamma_{dyn}$ に対する $R_{id}/R_{vir}$ の傾きが急になります。これは自己相似的な球対称崩壊モデルの予測と一致しますが、$R_{sp}$ の振る舞いとは異なります。

B. 主要合併の影響

• 合併の役割: 最近の主要合併を経験したハローを除外すると、$R_{id}/R_{vir}$ の赤方偏移依存性が大幅に弱まり、単一の曲線（$\Omega_m \approx 1$ の高赤方偏移の振る舞いに近い）で記述できるようになります。
• メカニズム: 合併はハローの動的構造を乱し、一貫した流入（inflow）を強化して外向きの成分（outflow）の形成を抑制します。その結果、MFR プロファイルの極小値が内側にシフトし、$R_{id}$ が小さくなります。特に高速降着（$\Gamma_{dyn} \gtrsim 1$）の領域でこの効果が顕著です。

C. 二次的なハロー性質への依存性

• 緩やかな降着 ($\Gamma_{dyn} \lesssim 1$): $R_{id}$ は主に $\Gamma_{dyn}$ によって決定され、ハローの形状、スピン、濃度、形成時間への依存性はほとんど見られません。
• 急速な降着 ($\Gamma_{dyn} \gtrsim 1$): 二次的な性質への依存性が顕在化します。特に、中心サブハローの半質量形成時間 ($a_{1/2}$) との依存性が強く、形成が早い（$a_{1/2}$ が小さい）ハローほど $R_{id}$ が小さくなります。これは、形成が早いハローが現在の急速な質量増加を「平滑な降着」ではなく「合併（特に衛星の降着）」によって達成していることを示唆しており、合併駆動型の成長が $R_{id}$ に影響を与えることを意味します。

D. $R_{id}/R_{vir} \simeq 2$ という普遍性の説明

• 以前に報告された $R_{id}/R_{vir} \simeq 2$ というほぼ一定の比率は、物理定数によるものではなく、統計的な平均効果です。
• 赤方偏移が進むにつれて、典型的なハローの降着率 $\Gamma_{dyn}$ は増加し（$R_{id}/R_{vir}$ を減少させる要因）、一方で同じ $\Gamma_{dyn}$ における $R_{id}/R_{vir}$ の関係自体が赤方偏移とともに上昇します（$R_{id}/R_{vir}$ を増加させる要因）。
• この 2 つの相反するトレンドが互いに相殺し合うことで、赤方偏移全体を通じて $R_{id}/R_{vir}$ がほぼ一定に保たれるという結果が導かれました。

4. 意義 (Significance)

• 新しい物理的プローブ: 枯渇半径 ($R_{id}$) は、ハローの成長率だけでなく、詳細な降着モード（平滑な降着か、合併駆動か、異方的か） に対する感度が高いことを示しました。これは、スプラッシュバック半径を補完する、ハローの動的状態と成長履歴を調べる強力な手段となります。
• 理論モデルへの示唆: 球対称崩壊モデルの予測とシミュレーション結果の一致・不一致を明らかにし、特に合併や非球対称性がハロー境界に与える影響を定量化しました。
• 観測への応用: 弱重力レンズや衛星の運動学を通じて観測可能な枯渇半径の特性を理論的に裏付けることで、将来の観測データからハローの降着履歴や合併履歴を推測する可能性を開きました。

この研究は、ダークマターハローの境界が単なる幾何学的な定義ではなく、ハローの複雑な成長履歴（特に最近の降着モード）を刻み込んだ物理的な量であることを確立しました。