[論文レビュー] Characterising ultra-high-redshift dark matter halo demographics and assembly histories with the GUREFT simulations

L. Y. Aaron Yung, Rachel S. Somerville+4 more|arXiv (Cornell University)|Sep 25, 2023

Galaxies: Formation, Evolution, Phenomena被引用数 4

ひとこと要約

本稿では、赤方偏移 z ≈ 20 から 6 の範囲で、1024³ 個の粒子を有する 4 つのボックスサイズ（5–90 Mpc/h）の高解像度、ダークマターのみの宇宙論的シミュレーション・スイート、GUREFT を提示する。このシミュレーションは、超高赤方偏移銀河団の特徴と形成歴を特徴づけるために用いられる。z ≳ 10 において、従来の解析的モデルとは顕著に異なる、更新された正確な銀河団質量関数、濃縮度、スピン、降着率のフィッティング関数が得られている。

ABSTRACT

Dark matter halo demographics and assembly histories are a manifestation of cosmological structure formation and have profound implications for the formation and evolution of galaxies. In particular, merger trees provide fundamental input for several modelling techniques, such as semi-analytic models (SAMs), sub-halo abundance matching (SHAM), and decorated halo occupation distribution models (HODs). Motivated by the new ultra-high-redshift (z > 10) frontier enabled by JWST, we present a new suite of Gadget at Ultrahigh Redshift with Extra-Fine Timesteps (GUREFT) dark matter-only cosmological simulations that are carefully designed to capture halo merger histories and structural properties in the ultra-z universe. The simulation suite consists of four 1024^3-particle simulations with box sizes of 5, 15, 35, and 90 Mpc h-1, each with 170 snapshots stored between 40 > z > 6. With the unprecedented number of available snapshots and strategically chosen dynamic range covered by these boxes, gureft uncovers the emerging dark matter halo populations and their assembly histories in the earliest epochs of cosmic history. In this work, we present the halo mass functions between z ~ 20 to 6 down to log(Mvir/Msun) ~ 5, and show that at high redshift, these robust halo mass functions can differ substantially from commonly used analytic approximations or older fitting functions in the literature. We also present key physical properties of the ultra-z halo population, such as concentration and spin, as well as their mass growth and merger rates, and again provide updated fitting functions.

研究の動機と目的

• ジェイ・ウィルビーエス・テレスコープ（JWST）による観測が、宇宙の初期を探索し始めている超高赤方偏移（z ≳ 10）におけるダークマター団の特徴と形成歴をモデル化すること。
• 宇宙初期の時代（z ≈ 20 から 6）をカバーする高解像度・高スナップショット密度のシミュレーションが不足している現状を是正すること。このようなシミュレーションは、正確な銀河形成モデルの構築に不可欠である。
• 広い質量および赤方偏移範囲にわたって有効である、質量関数、濃縮度、スピン、降着率といった銀河団の性質の、信頼性の高いシミュレーションベースのフィッティング関数を提供すること。
• 最初の銀河の高精度な合体木および構造的データを供給することで、半解析的モデル、準銀河団数マッチング、銀河団占有分布技術の精度を向上させること。
• 極端な銀河団性質（例：高濃縮度や高スピン）に及ぼす質量解像度の影響を評価する収束性の検証を実施すること。

提案手法

• 1024³ 個の粒子を有する 4 つの高解像度、ダークマターのみのシミュレーションを実行し、ボックスサイズは 5、15、35、90 Mpc/h であり、z = 40 から z = 6 の間で 170 スナップショットを取得する。
• Gadget-2 を用い、超高赤方偏移領域における銀河団形成および合体歴の正確な追跡を保証するため、超微細なタイムステップを適用する。
• 標準的な銀河団検出アルゴリズム（例：位相空間フォローアップ・フォリオ、FoF）を用いて銀河団を同定し、半径質量、濃縮度、スピン、円運動速度を測定する。
• 銀河団質量関数（HMF）を導出し、全赤方偏移および質量範囲に適合する新しい関数形（累乗則および指数補正を含む）を用いてフィッティングする。
• 合体率、降着率（dM/dt および ḡM/M）、およびスケーリング関係（例：Mvir–Vmax、Mvir–c、Mvir–λ）を赤方偏移および質量にわたって計算する。
• 隣接するシミュレーションボックス間の重複する質量および赤方偏移領域を比較することで、極端な銀河団性質に及ぼす解像度効果を評価する収束性の検証を実施する。

実験結果

リサーチクエスチョン

• RQ1 z ≳ 10 における銀河団質量関数は、文献に既存の解析的近似や古いフィッティング関数と比較してどのように異なるか？
• RQ2 超高赤方偏移宇宙における銀河団の濃縮度、スピン、最大円運動速度の赤方偏移傾向は何か？
• RQ3 銀河団の降着率および合体歴は、z ≈ 6 から z ≈ 20 にかけてどのように変化するか？また、低赤方偏移シミュレーションからの外挿と比較するとどうなるか？
• RQ4 Mvir–Vmax、Mvir–c、および Mvir–dM/dt の既存のフィッティング関数は、高赤方偏移および低質量領域でどの程度破綻するか？新たなフィッティングは必要か？
• RQ5 質量解像度は、初期宇宙における極端な銀河団（例：高濃縮度または高スピン）の集団にどのような影響を及ぼすか？

主な発見

• GUREFT の銀河団質量関数（HMF）は、log(Mvir/M☉) ≈ 5 から 12 の範囲、z ≈ 6 から 20 の範囲をカバーし、広い動的範囲を正確に記述する新しいフィッティング関数が得られた。
• 一般的に用いられる解析的 HMF モデル（例：拡張プレス＝シェューティング形式に基づくもの）は、z ≳ 10 において GUREFT の結果と最大で 1 デキシタル・オクテット（1 dex）の差を示し、超高赤方偏移領域ではその不正確さが顕著に現れる。
• Mvir–Vmax 関係は、z = 0 から z = 20 にかけて約 0.25 デキシタル・オクテット増加しており、これは同じ質量に対して高赤方偏移の銀河団がより高い円運動速度を持つことを示している。
• 銀河団スピン（λ）は赤方偏移の増加に伴い減少する。λB はやや弱い進化を示し、λP はより強い進化を示す。z ≳ 10 におけるこれらの分布のフィッティング関数は、本稿で初めて提供される。
• Mvir–cVir 関係は、z ≈ 6 から 20 の間で平坦なまま維持され、初期の時点でわずかに正規化値が増加する。これは、同じ質量に対して z ≈ 20 の銀河団が低赤方偏移時よりもより濃縮されていることを意味する。
• 特定降着率（Ṁ/M）は、z ≈ 6 から 20 にかけて単調に増加する。Ṁ/M–z 関係の勾配は高赤方偏移域でわずかに平坦化するが、トレンドは継続する。既存の低赤方偏移フィットは、この領域への外挿に不適切である。