STEPSIC: Initial condition generator for stereographic cosmological simulations

ステレオグラフィック宇宙論シミュレーションの初期条件生成器 stepsic の技術的概要

問題提起
従来の宇宙論的初期条件（IC）生成器は、完全に周期的な立方体領域（$T^3$ トポロジー）向けにのみ設計されている。標準的な構造形成研究では成功を収めているものの、$T^3$ トポロジーはボックスサイズに匹敵するスケールで人工的な相関を導入し、完全回転群 $SO(3)$に対する不変性の欠如により角運動量を保存しない。これらの限界により、標準的なツールは非周期的で観測者中心のシミュレーションや、回転不変性を必要とする現象には不適切となる。具体的には、逆ステレオグラフィック射影を用いて$R^3$多様体を、半径方向に変化する粒子分解能を持つ有限の球状（$R^3$）および円筒状（$S^1 \times R^2$）領域へ写像する StePS（STEreographically Projected cosmological Simulations）コードには、互換性のある公開 IC 生成器が存在しない。既存のコード（例：2lptic、music、monofonic）は、ステレオグラフィック射影に固有の非標準的な幾何学形状を生成したり、内在する質量分解能の勾配を処理したりすることができない。

手法
本論文では、このギャップを埋めるために設計されたオープンソースの初期条件生成器 stepsic を提示する。このコードは、ラグランジュ摂動論（LPT）を拡張し、球状、円筒状、および任意のアスペクト比を持つ直方体領域をサポートする。技術的なワークフローは以下の通りである：

1. 粒子負荷生成：stepsic は 4 つの負荷戦略をサポートする：

   • 規則格子：周期的 $T^3$ 領域用の立方体ボクセル。
   • ポアソン分布：テストやグラス生成のための均一なランダムサンプリング。
   • 同心シェル/円環：StePS 幾何学に固有のもので、体積を固定粒子数を含むシェル（球状）または円環（円筒状）で充填する。これにより、ステレオグラフィック射影に固有の半径方向に変化する質量分解能が作成される。2 つのビンニングモードがサポートされている：一定のステレオグラフィック角度（$\Delta\omega$）と一定のコンパクト空間体積。巨大な質量差を持つ粒子間の相互作用に起因する力精度の劣化を緩和するため、一定分解能の内部領域が定義されている。
   • グラス入力：外部で生成されたグラス構成を読み込み、スケーリングする機能。

2. 場構築：

   • フーリエ格子：異方性のない立方体ボクセルを持つフーリエ格子が構築される。非立方体領域の場合、格子はボックス幾何学によって条件付けられ、単に大きな等方性実現を切り取るのではなく、アクセス可能なモードのアンサンブルを変更する（短軸に沿った分散を抑制する）。
   • 密度場：白色ノイズを用いてガウスランダム場を生成し、初期赤方偏移における目標線形物質パワースペクトル $P(k, z)$ にスケーリングする。このコードは、標準的な実現、ペアリングされた実現（位相シフト付き）、および宇宙論的分散を低減するための固定振幅実現をサポートする。
   • LPT 実装：コードは 1 次（ツェルドビッチ近似）および 2 次（2LPT）のラグランジュ摂動論を実装する。変位場と速度場はフーリエ空間で計算され、物理空間へ変換される。実装は 2lptic および monofonic によって確立された規約に従い、標準的なワークフローとの互換性を保ちつつ、非標準的なトポロジーへ拡張している。

3. 補間と多分解能：

   • 変位場と速度場は、平滑化を補正するためのフーリエ空間デコンボリューションを伴う B スプライン質量割り当てカーネル（NGP、CIC、TSC）を用いて粒子位置へ補間される。
   • 多分解能スキームは、半径方向に変化する質量分解能全体にわたって変位場をマッピングする。粒子は、その質量に対応する様々な分解能の格子から変位場を割り当てられ、より大きな体積を表す高質量粒子が、ノイズとして作用する可能性のある未解決の小規模変動を受け取らないようにしている。

主要な結果と検証
著者は、再現性、数値的正確性、理論および他のコードとの統計的整合性という 3 つのレベルのテストを通じて stepsic を検証する。

• パワースペクトルの回復：周期的な立方体ボックスにおいて、回復された物質パワースペクトルは、ナイキスト波数の半分まで、入力線形理論予測と 0.5% 以内で一致する。この精度は、ボックスのアスペクト比（10:1 までテスト）に依存しない。
• 相互検証：同一の白色ノイズ場と宇宙論的パラメータを用いて確立された monofonic コードと比較した際、stepsic はパーセント未満のパワースペクトル整合性を達成した。わずかな残留オフセット（約 0.6%）は、線形成長因子計算の実装の違い（stepsic では colossus ライブラリを使用、monofonic では内部ソルバーを使用）に起因すると考えられる。
• 過渡現象の抑制：一致する円筒状 StePS 実行の完全な N 体進化により、準線形スケールにおいて 1 次初期条件に存在する 2–3% の過渡的なパワースペクトル過剰 が、2 次 LPT によって正しく抑制されることが確認された。
• 幾何学的整合性：スラブ幾何学に対するテストは、アスペクト比が 10:1 までにおいて、立方体ボクセルフーリエ格子の構築が、より小さな体積におけるサンプル分散の期待される減少を除き、幾何学依存の系統的バイアスを導入しないことを示している。

意義と主張
本論文は、stepsic を、StePS コードが必要とする非標準的な幾何学（$R^3$ および $S^1 \times R^2$）と半径方向に変化する質量分解能をサポートできる最初の公開初期条件生成器として位置づけている。これらの領域へ LPT を成功裡に拡張し、出力を理論的期待値および既存のコードと比較して検証することにより、著者らは、周期的 $T^3$ トポロジーに内在する人工的な相関なしに、正確な観測者中心の宇宙論シミュレーションを実行可能であることを実証している。このコードは、StePS フレームワーク内でのフィラメント構造および異方性宇宙論モデルの研究に必要なツールとして提示されている。著者は控えめに、3 次 LPT および原始非ガウス性は将来のリリースで予定されているものの、現在の実装はステレオグラフィックシミュレーションのための堅牢な 2 次基盤の提供に焦点を当てていると指摘している。