[論文レビュー] The art of simulating the early Universe -- Part I

Daniel G. Figueroa, Adrien Florio+2 more|arXiv (Cornell University)|Jun 26, 2020

Cosmology and Gravitation Theories参考文献 275被引用数 13

ひとこと要約

本稿では、高精度な格子技術を用いて、膨張する初期宇宙におけるスカラー場およびゲージ場のダイナミクスをシミュレートする包括的なフレームワークを提示する。O(δt²) から O(δt¹⁰) のシンプレクティック積分法を導入し、マシン精度でゲージ制約（Gauss制約）を保存する。これにより、自己自己一致的な膨張を伴う3+1次元における非線形ダイナミクスの正確なシミュレーションが可能になる。主な貢献は、SU(2)×U(1)プリヒーティングなどのインフレーション後ダイナミクスを研究するための、信頼性が高く公開可能なC++パッケージ、CosmoLatticeの確立である。

ABSTRACT

We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of $O(dt^2)$ accuracy like $staggered~leapfrog$ and $Verlet$ integration, to $Runge-Kutta$ methods up to $O(dt^4)$ accuracy, and the $Yoshida$ and $Gauss-Legendre$ higher-order integrators, accurate up to $O(dt^{10})$. We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in $3+1$ dimensions, including the case of `self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: $i)$ Interacting scalar fields, $ii)$ Abelian $U(1)$ gauge theories, and $iii)$ Non-Abelian $SU(2)$ gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from $O(dt^2)$ up to $O(dt^{10})$. For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. Remarkably, all our algorithms for gauge theories respect the Gauss constraint to machine precision, including when `self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under $SU(2) imes U(1)$. The present manuscript is meant as part of the theoretical basis for $CosmoLattice$, a modern C++ MPI-based package for simulating the non-linear dynamics of scalar-gauge field theories in an expanding universe, publicly available at www.cosmolattice.net

研究の動機と目的

• 膨張するフレリッド・レマトゥル・ロバートソン・ウォーカー（FLRW）宇宙におけるスカラー場およびゲージ場のダイナミクスをシミュレートするための、数値的に安定かつ高精度な格子積分技術の開発。
• 特に、場のエネルギーと圧力によって駆動される自己自己一致的な膨張下でも、ゲージ不変性と制約保存を維持する挑戦への対処。
• 古典的格子場理論への応用を想定し、2次から10次までの精度を有する数値アルゴリズムの体系的分類と実装。
• CosmoLatticeパッケージの理論的・計算的基盤を確立し、プリヒーティングなどの非線形的初期宇宙現象を高精度にシミュレート可能にする。

提案手法

• 自然単位と平坦なFLRW計量を採用し、膨張する時空におけるスカラー場およびゲージ場（U(1)およびSU(2)相互作用を含む）の連続場方程式を導出する。
• 空間格子上での場の方程式の離散化を実施し、離散フーリエ変換を用いた格子微分作用素および運動量表現を実装。スタガードおよびウィルソン型形式により、ゲージ不変性を保証する。
• 積分法の階層を構築：2次ステップのスタガードリープフロッグおよびバーレット法から、4次精度の明示的ルンゲ＝クッタ法、および10次精度の陰的ガウス＝レジェンドルおよびヨシダ法まで、O(δt¹⁰) までの多様な方法を含む。
• ゲージ理論の文脈では、マシン精度で制約を保存するシンプレクティック積分法を構築し、自己自己一致的シミュレーションにおいてもゲージ制約（Gauss制約）を厳密に満たす。
• エネルギー密度成分、場スペクトル、ハッブル制約などの観測量を、格子正則化された期待値を用いて計算する。
• 実用例として、電荷を有するインフレートンを伴うSU(2)×U(1)プリヒーティングを提示。複数の積分法の次数において、精度と制約保存の両方が確認された。

実験結果

リサーチクエスチョン

• RQ1 非線形的スカラー場およびゲージ場のダイナミクスを、最小限の数値誤差で、膨張宇宙において10次精度までのシンプレクティック積分法をどのように設計・実装できるか？
• RQ2 特に、場自身が膨張を駆動する非アーベルゲージ理論の格子シミュレーションにおいて、Gauss制約の破れをどの程度抑制できるか？
• RQ3 3+1次元格子場理論シミュレーションにおいて、明示的（例：ルンゲ＝クッタ）と陰的（例：ガウス＝レジェンドル）積分法の間で、性能と精度のトレードオフはどのようなものか？
• RQ4 相互作用するスカラー場およびゲージ場の格子形式を、エネルギー保存則やゲージ不変性といった重要な物理的制約を保持するように体系的に構築できるか？
• RQ5 2次から10次までの異なる積分法次数（O(δt²) から O(δt¹⁰)）が、インフレーション後ダイナミクスにおける場スペクトルおよびエネルギー密度の進化にどのように影響を与えるか？

主な発見

• 著者らは、3+1次元格子シミュレーションにおいて、スカラー場およびゲージ場に対してO(δt¹⁰)までのシンプレクティック積分法を成功裏に実装した。時間ステップ重みの明示的式は、表1～3に提示された。
• すべてのゲージ理論用積分法は、自己自己一致的膨張を伴うシミュレーションにおいても、マシン精度でGauss制約を保存する。これは物理的整合性を保つために極めて重要な成果である。
• ガウス＝レジェンドル法およびヨシダ法は、s = 2～5段階のバーサー図表現を用いて検証された安定的かつ対称的な時間ステッピングスキームを実現し、O(δt¹⁰)の精度を達成した。
• 電荷を有するインフレートンを伴うSU(2)×U(1)プリヒーティングのシミュレーションにおいて、複数の積分法次数において一貫したエネルギーおよびスペクトルの進化が確認され、数値的安定性と精度が裏付けられた。
• 本フレームワークに基づくCosmoLatticeパッケージの実装により、非アーベルゲージ場および自己自己一致的膨張を伴う初期宇宙ダイナミクスの公開・高精度なシミュレーションが可能になった。
• 場スペクトルおよびエネルギー密度成分の計算において高い精度が達成され、複数の積分法次数にわたる収束テストにより誤差制御が検証された。