[論文レビュー] A dynamics-based density profile for dark haloes. I. Algorithm and basic results

Benedikt Diemer|arXiv (Cornell University)|Dec 7, 2021

Galaxies: Formation, Evolution, Phenomena参考文献 166被引用数 36

ひとこと要約

本論文は、N体シミュレーションの粒子の近点通過を追跡することで、ダークマター銀河団の密度プロファイルを軌道回り成分と降着成分に分離する、動力学的アプローチを提案する。その結果、軌道回り項はスプラッシュバック半径で急峻な切り捨てを示し、主に質量降着率によって支配され、散らばりが非常に小さい（0.1–0.4 dex）一方、降着項は環境要因や降着歴に強く依存し、単純な理論的モデルに挑戦するものである。

ABSTRACT

The density profiles of dark matter haloes can potentially probe dynamics, fundamental physics, and cosmology, but some of the most promising signals reside near or beyond the virial radius. While these scales have recently become observable, the profiles at large radii are still poorly understood theoretically, chiefly because the distribution of orbiting matter (the one-halo term) is partially concealed by particles falling into halos for the first time. We present an algorithm to dynamically disentangle the orbiting and infalling contributions by counting the pericentric passages of billions of simulation particles. We analyse dynamically split profiles out to 10 R200m across a wide range of halo mass, redshift, and cosmology. We show that the orbiting term experiences a sharp truncation at the edge of the orbit distribution. Its sharpness and position are mostly determined by the mass accretion rate, confirming that the entire profile shape primarily depends on halo dynamics and secondarily on mass, redshift, and cosmology. The infalling term also depends on the accretion rate for fast-accreting haloes but is mostly set by the environment for slowly accreting haloes, leading to a diverse array of shapes that does not conform to simple theoretical models. While the resulting scatter in the infalling term reaches 1 dex, the scatter in the orbiting term is only between 0.1 and 0.4 dex and almost independent of radius. We demonstrate a tight correspondence between the redshift evolution in LCDM and the slope of the matter power spectrum. Our code and data are publicly available.

研究の動機と目的

• 銀河団の外縁部における軌道回り成分と降着成分を分離する理論的課題に取り組む。
• 内側銀河団と大規模構造の間の遷移領域におけるプロファイル形状の起源を理解する。
• ハロー質量、降着率、赤方偏移、および宇宙論的パラメータに依存する外縁部密度プロファイルの定量的評価を行う。
• N体シミュレーションデータの動的分解に耐久的かつ公開可能なフレームワークを提供する。

提案手法

• N体シミュレーションにおける粒子の近点通過回数を数えることで、動的歴史に基づいて軌道回りまたは降着と分類するアルゴリズム。
• 任意の半径カットを避ける位相空間ベースの粒子分類を用い、1ハロー（軌道回り）項と降着（降着中）項の正確な分離を可能にする。
• Erebusシミュレーション・スイートにおいて、多様な質量、赤方偏移、宇宙論的条件下の10,000個のハローにこの手法を適用。
• 10 R200mまでの球殻でプロファイルを計算し、質量、力、時間分解能における数値収束性に注意を払う。
• ハロー境界定義や合体除外の感度にかかわらず、収束テストと感度解析により動的分離の妥当性を検証。
• コード名はSpartaで、オープンソースであり、分解済みプロファイルの読み取りと解析に使えるツールを含む。

実験結果

リサーチクエスチョン

• RQ1 外縁部における軌道回り密度プロファイルの形状は、ハローの降着率と宇宙論にどのように依存するか？
• RQ2 軌道回りおよび降着成分の密度プロファイルにおける散らばりの原因は何か？
• RQ3 ΛCDMにおけるハローの赤方偏移進化は、異なる物理的サイズのハローが探査する物質パワースペクトルの有効勾配によって説明可能か？
• RQ4 降着成分と軌道回り成分は、ハローの環境および成長歴にどのように異なる依存性を示すか？
• RQ5 NFW や Einasto といった標準的なフィッティング関数は、スプラッシュバック半径付近の遷移領域をどの程度うまく記述できないか？

主な発見

• 軌道回り項は、アポセントル分布の端に相当するスプラッシュバック半径で急峻な切り捨てを示し、これは必ずしも全プロファイルの急峻な点とは一致しない。
• 軌道回りプロファイルの形状は主に質量降着率によって決定され、ゆっくり降着するハローでは徐々に傾きがきつくなり、急速な降着者では累乗則的挙動に近づく。
• 軌道回りプロファイルの対数的散らばりは小さく（0.1–0.4 dex）、半径にほぼ依存せず、ハロー間で高い動的整合性を示す。
• 降着項の散らばりは大きく（最大1 dex）、主にゆっくり降着するハローでは環境要因、急速な降着者では降着率が主因であり、普遍的なフィッティングモデルの構築を困難にしている。
• ΛCDMにおけるハローの赤方偏移進化は、ハローの物理的サイズに応じて探査する物質パワースペクトルの有効勾配によって完全に説明可能である。
• 動的分離アルゴリズムは解像度スケールにわたって収束し、ハロー境界定義や合体除外の取り扱いに対してもロバストである。コードとデータは公開済みである。