[论文解读] 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¹⁰) 的辛积分器，能够以机器精度保持高斯约束，从而实现对 3+1D 中非线性动力学（包括自洽膨胀）的精确模拟。其主要贡献是一个稳健且公开可用的 C++ 软件包 CosmoLattice，用于研究如 SU(2)×U(1) 重加热等暴胀后动力学。

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）宇宙中标量场和规范场动力学的数值稳定且高精度的格点积分技术。
• 解决在数值模拟中保持规范不变性和约束保持的挑战，特别是在场能量和压强驱动的自洽膨胀情况下。
• 系统性地分类并实现从二阶到十阶精度的数值算法，用于经典格点场论。
• 为 CosmoLattice 软件包建立理论和计算基础，以实现对非线性早期宇宙现象（如重加热）的高保真模拟。

提出的方法

• 采用自然单位制和平坦 FLRW 度规，作者推导了在膨胀时空中包含 U(1) 和 SU(2) 相互作用的标量场和规范场的连续场方程。
• 在空间格点上离散化场方程，通过离散傅里叶变换实现格点导数算符和动量表示，利用叠铺和威尔逊型形式确保规范不变性。
• 实现了一套分层的数值积分器：从二阶叠铺跳跃和 Verlet 方法，到四阶显式 Runge-Kutta 方法，以及最高达 O(δt¹⁰) 的隐式 Gauss-Legendre 和 Yoshida 方法。
• 对于规范理论，通过构建在机器精度下保持约束的辛积分器，数值强制执行高斯约束，即使在自洽模拟中亦然。
• 该框架包括能量密度分量、场谱和哈勃约束等可观测量，通过格点正则化期望值计算。
• 通过一个工作示例（带电暴胀子的 SU(2)×U(1) 重加热）验证了该方法，展示了在多个积分器阶次下均保持高精度和约束保持性。

实验结果

研究问题

• RQ1 如何设计并实现高阶辛积分器，以在最小数值误差下模拟膨胀宇宙中非线性标量场和规范场的动力学？
• RQ2 在非阿贝尔规范理论的格点模拟中，尤其当膨胀由场本身驱动时，高斯约束违反能在多大程度上被抑制？
• RQ3 在 3+1D 格点场论模拟中，显式与隐式积分器（如 Runge-Kutta 与 Gauss-Legendre）之间的性能与精度权衡如何？
• RQ4 能否系统性地构建相互作用标量场和规范场的格点形式，以保持能量守恒和规范不变性等关键物理约束？
• RQ5 从 O(δt²) 到 O(δt¹⁰) 的不同积分器阶次如何影响暴胀后动力学中场谱和能量密度的演化？

主要发现

• 作者成功实现了在 3+1D 格点模拟中针对标量场和规范场的最高达 O(δt¹⁰) 精度的辛积分器，并在表 1–3 中提供了时间步长权重的显式公式。
• 所有规范场理论积分器均以机器精度保持高斯约束，即使在自洽膨胀模拟中亦然，这是维持物理一致性的关键成就。
• Gauss-Legendre 和 Yoshida 方法实现了 O(δt¹⁰) 精度，采用稳定且对称的时间步长方案，通过 s = 2 到 5 阶段的 Butcher 表格表示得到验证。
• 对带电暴胀子的 SU(2)×U(1) 重加热模拟显示，不同积分器阶次下能量和谱演化保持一致，证实了数值稳定性和精度。
• 基于该框架实现的 CosmoLattice 软件包，使公众能够进行高精度的早期宇宙动力学模拟，包括非阿贝尔规范场和自洽膨胀。
• 该方法在计算场谱和能量密度分量方面实现了高精度，误差控制通过在多个积分器阶次下的收敛性测试得到验证。