[論文レビュー] CosmoLattice

Daniel G. Figueroa, Adrien Florio+2 more|arXiv (Cornell University)|Feb 1, 2021

Educator Training and Historical Pedagogy被引用数 7

ひとこと要約

CosmoLattice は、膨張宇宙におけるスカラー場およびゲージ場のダイナミクスをモデル化するための現代的な C++ ベースで MPI 並列化された格子シミュレーションフレームワークであり、記号的場言語、高次時間積分法（最大 O(δt¹⁰)）、ゲージ理論におけるガウスの法則の正確な強制を特徴とする。これは、プリヘーティング や相転移などの初期宇宙現象の高分解能・長時間シミュレーションを可能にし、再起動機能および HDF5 を用いた場の保存機能を完全にサポートする。

ABSTRACT

This is the user manual for CosmoLattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoLattice incorporates a series of features that makes it very versatile and powerful: $i)$ it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, $ii)$ it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well-separated scales, running very high resolution simulations, or simply very long ones, $iii)$ it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, $iv)$ it includes a library of numerical algorithms, ranging from $O(δt^2)$ to $O(δt^{10})$ methods, suitable for simulating global and gauge theories in an expanding grid, including the case of `self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g.~energy densities, field spectra, lattice snapshots) and we note that remarkably all our algorithms for gauge theories always respect the Gauss constraint to machine precision. In this manual we explain how to obtain and run CosmoLattice in a computer (let it be your laptop, desktop or a cluster). We introduce the general structure of the code and describe in detail the basic files that any user needs to handle. We explain how to implement any model characterized by a scalar potential and a set of scalar fields, either singlets or interacting with $U(1)$ and/or $SU(2)$ gauge fields. CosmoLattice is publicly available at www.cosmolattice.net.

研究の動機と目的

• 初期宇宙におけるスカラー場およびゲージ場の非摂動的ダイナミクスを研究するための柔軟で高性能なシミュレーションツールを提供すること。
• 広範なスケールの隔たりを持つ宇宙論的シナリオの高分解能・長時間格子シミュレーションを可能にすること。
• アルゴリズム的設計により、ゲージ理論における正確な制約満足（例：ガウスの法則）を保証すること。
• 記号的場言語とモジュラーなコード構造を通じたユーザーフレンドリーなモデル定義を可能にすること。
• MPI 並列化と HDF5 を用いたチェックポイントおよび再起動機能を通じて、再現可能性とスケーラビリティを促進すること。

提案手法

• コードはオブジェクト指向設計を用いて C++ で実装されており、物理的記述と数値的実装の明確な分離がなされている。
• 効率的な多次元空間並列化のため、MPI と並列化された離散フーリエ変換（PFFT）が使用されている。
• 記号的場言語により、連続的表記に近く、ユーザーが場や微分作用素を定義できる。
• O(δt²) から O(δt¹⁰) までの数値的時間発展アルゴリズムのライブラリが実装されており、グローバル理論およびゲージ理論の両方に適している。
• U(1) および SU(2) ゲージ理論におけるガウス制約は、すべてのアルゴリズムで機械精度まで強制されている。
• HDF5 I/O サポートにより、3次元場の全配置、スケール因子の時間発展、シミュレーションパラメータの保存と再起動が可能である。

実験結果

リサーチクエスチョン

• RQ1 どのようにして、膨張宇宙におけるスカラー場およびゲージ場の格子シミュレーションを、高精度かつ長時間スケールにまでスケーラブルに実現できるか？
• RQ2 ゲージ不変性および制約の満足を保ちながら、高次時間積分法を効果的に実装する最良の方法は何か？
• RQ3 記号的場言語をどのように用いることで、パフォーマンスや数値的正確性を損なわず、モデル定義を簡素化できるか？
• RQ4 MPI および PFFT を用いた並列化は、高分解能の宇宙論的シミュレーションのパフォーマンスをどの程度向上できるか？
• RQ5 大規模な格子シミュレーションに、データ整合性を損なわずに、堅牢なチェックポイントおよび再起動メカニズムを統合するにはどうすればよいか？

主な発見

• ベンチマークでは、並列効率が約 99%（α ≈ 0.99）に達しており、コア数に応じたほぼ理想のスケーリングが実現されている。
• CosmoLattice 内のすべてのゲージ理論アルゴリズムが、機械精度までガウス制約を強制しており、シミュレーション全体を通じてゲージ不変性が保たれている。
• 記号的場言語により、ユーザーは連続数学に近い表記で複雑な場の式や作用素を定義できる。
• HDF5 を用いた I/O 機能により、3次元場の分布や時間発展データを含む、信頼性の高いシミュレーションの保存と再起動が可能である。
• フレームワークは、場の反作用によるスケール因子の動的変化を伴う自己自己一致的な拡張を伴うモデルのシミュレーションをサポートしている。
• コードは http://www.cosmolattice.net で公開されており、ラップトップ、デスクトップ、ハイパフォーマンスクラスタのすべてで使用可能である。