Cylindrical cosmological simulations with StePS

技術要約：StePSを用いた円柱形宇宙論シミュレーション

問題提起
標準的な宇宙論的$N$体シミュレーションは、通常、周期的な立方体ボックス内で物質を進化させる。これは三次元トーラス（$T^3 = S^1 \times S^1 \times S^1$）のトポロジーを課すものである。このトポロジーは、数値計算上の利便性が高く統計的一様性を保持する一方で、幾何学が連続回転群$SO(3)$ではなく八面体群（$O_h$）に対してのみ不変であるため、統計的な等方性を破ってしまう。その結果、周期的なボックス内では角運動量が一般に保存されず、立方体の対称性が重力場に微妙な異方性を導入し、特にシミュレーション・ボリュームが関心のスケールに対して十分に大きくない場合、大規模構造の統計量にバイアスを与える可能性がある。さらに、標準的な周期境界条件は、宇宙のフィラメントや特定の異方的な宇宙モデルのような、固有の円柱対称性を持つ系のシミュレーションには適していない。

手法
著者らは、$S^1 \times \mathbb{R}^2$（スラブ）トポロジカル多様体に基づいたコンパクト化されたシミュレーション・フレームワークを導入している。この幾何学において、シミュレーション領域は単一の軸（$z$軸、周期$L_z$）に沿って周期的であるが、半径方向（$\mathbb{R}^2$）には開いている。

• 幾何学的コンパクト化: 著者らは、非コンパクトな$\mathbb{R}^2$平面をコンパクト化するために、半径方向への逆ステレオ投影を用いている。これにより、無限の半径座標を有限の領域へと写像し、「ズームイン」戦略が可能となる。これにより、中心軸付近では空間および質量分解能が高くなり、外側の境界に向かって滑らかに減少する。これは、自由な外側境界におけるエッジ・アーティファクトを抑制する。
• 運動方程式: 宇宙論的共動座標系においてニュートン運動方程式を解く。重力計算は、$z$軸に沿った周期的な像（イメージ）を足し合わせることで、$S^1 \times \mathbb{R}^2$のトポロジーを考慮する。直接和の収束の遅さを処理するため、著者らは円柱幾何学に適応させたエバルト（Ewald）総和法を採用し、力を実空間成分と逆空間成分に分割している。
• 力計算アルゴリズム: StePS（STEreographically Projected cosmological Simulations）フレームワーク内の実装は、2つの力計算法をサポートしている。
  1. 直接 $O(N^2)$ 法: GPUへのマッピングが効率的であり、より少ない粒子数やより小さな円柱状のボリュームに適している。
  2. オクトリー $O(N \log N)$ 法: CPUクラスター向けに実装されている。円柱領域におけるオクトリー・グリッドの立方体的な性質から生じる系統誤差を軽減するため、コードは毎ステップごとに領域に対してランダムなシフトと回転を適用する。重要な較正ステップとして、背景密度の安定性を確保し、偽の半径方向の崩壊または膨張を防ぐために、半径方向の力誤差を補正する。
• 初期条件: 初期条件は、$S^1 \times \mathbb{R}^2$トポロジーにおける斥力重力によって生成された「グラス（ガラス）」粒子分布を用いて生成される。摂動は、円柱幾何学およびズームイン設定の非一様なサンプリングに適応させたゼルドビッチ近似または二次ラグランジュ摂動論（2LPT）を用いて刻み込まれる。パワースペクトルは、非自明なウィンドウ関数を考慮するために、フェルドマン・カイザー・ピーケイ（FKP）推定法を用いて推定される。

主要な貢献

1. 新規なトポロジーの実装: 本論文は、$S^1 \times \mathbb{R}^2$多様体における宇宙論的$N$体シミュレーションの初の実装を提示しており、StePSコードを従来の$T^3$および球形$\mathbb{R}^3$の能力を超えて拡張している。
2. 保存則: 著者らは、このトポロジーにおいて、ラグランジアンが円柱軸に沿った並進および横平面内での回転（$SO(2)$）に対して不変であることを示している。その結果、$z$軸に沿った線運動量および同軸に沿った角運動量が保存され、$T^3$トポロジーにおける対称性の破れの問題に対処している。
3. アルゴリズムの適応: 力計算（円柱形エバルト総和）、初期条件生成（円柱形グラス）、およびハロー探索（$S^1$の周期性と半径方向の分解能変化に適応させたカスタム球面過密度ファインダー）に必要な修正の詳細を述べている。
4. 検証: 手法の妥当性を検証するために、フル$\Lambda$CDMシミュレーションが実行され、標準的な宇宙論的期待と一致する宇宙網（コスミック・ウェブ）の形成が示された。

結果

• シミュレーション実行: $L_z = 1.0$ Gpc、半径 $R_{sim} = 500$ Mpc の円柱内で、Planck 2018 $\Lambda$CDMパラメータを用いて、$2.4 \times 10^7$ 個の粒子によるシミュレーションを実行した。
• 構造形成: シミュレーションは、非線形な宇宙網（ボイド、壁、フィラメント、およびハロー）を正常に再現した。ズームイン戦略により、中心軸付近の低質量ハローを高精度に解像できる一方で、より大きな体積では、より大きな半径にあるより質量の大きいハローを解像できる。
• パワースペクトル: 推定された物質パワースペクトルは、大規模スケールでは線形CAMB予測と、小規模スケールでは非線形 Mead et al. (2021) の「halofit」エミュレータとよく一致しており、構造成長を捉える手法の正確性を裏付けている。
• ハロー特性: 著者らは $z=0$ において約 $1.1 \times 10^5$ 個のダークマターハローを特定し、各ハローについて質量、集中度、角運動量を含む68個のパラメータを算出した。
• 異方性テスト: 銀河のスピン異方性に関する最近の主張に基づき、著者らはハローのスピン軸と周期的な$z$軸との分布をテストした。コルモゴロフ・スミルノフ検定を用いた結果、$L_z = 1.0$ Gpcにおいて、スピン分布は統計的等方性と一致していることが分かった。彼らは、このスケールにおいては、このトポロジーが報告されている銀河スピンの非対称性を自然に説明するものではないと結論付けている。

意義と主張
本論文は、$S^1 \times \mathbb{R}^2$トポロジーが、標準的な立方体シミュレーションに対する特化した代替案を提供すると主張している。これは特に、円柱対称性を持つ系（例：フィラメント）や、特定の軸に沿った角運動量の保存が物理的に関連する場合の研究に適している。著者らは、この手法が独特の系統誤差（半径方向に依存する分解能や選択効果など）を導入するものの、$T^3$ボックスの立方体的な異方性から解放された自己整合的なフレームワークを提供することを強調している。

著者らは控えめな主張をしており、このアルゴリズムが一般的な宇宙論のための確立された立方体$T^3$シミュレーションに取って代わることは期待していない。むしろ、立方体の対称性が制限要因となる特定の科学的問いに対する補完的なツールとして機能する。また、トポロジーの異方的な力場は、周期長が小さすぎる場合（$L_z \lesssim 200$ Mpc）にアーティファクトを刻み込む可能性があることを指摘しており、これは今後の詳細な特性評価に委ねられている。コードは、これらの幾何学のさらなる探求を促進するために、オープンソース（StePS v2）として公開されている。