[论文解读] 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，可保存和重启完整的 3D 场配置、尺度因子演化和模拟参数。

实验结果

研究问题

• RQ1 如何使在膨胀宇宙中对标量场和规范场的晶格模拟既具有高度准确性，又能扩展到长时间尺度？
• RQ2 在保持规范不变性和约束满足的前提下，实现高阶时间积分器的最有效方法是什么？
• RQ3 如何利用符号场语言在不牺牲性能或数值精度的前提下简化模型定义？
• RQ4 MPI 与基于 PFFT 的并行化在多分辨率宇宙学模拟中能多大程度上提升性能？
• RQ5 如何在不损害数据完整性的前提下，将稳健的检查点与重启机制集成到大规模晶格模拟中？

主要发现

• 在基准测试中，该代码实现了约 99% 的并行效率（α ≈ 0.99），表明在多核处理器上具有近乎理想的可扩展性。
• CosmoLattice 中所有规范理论算法均以机器精度强制执行高斯约束，确保模拟全程保持规范不变性。
• 符号场语言允许用户以接近连续数学的语法定义复杂的场方程和算子。
• 基于 HDF5 的 I/O 支持可实现模拟的可靠保存与重启，包括完整的 3D 场分布和时间演化数据。
• 该框架支持自洽膨胀的模型模拟，其中尺度因子因场的反作用而动态演化。
• 该代码在 http://www.cosmolattice.net 公开提供，适用于笔记本电脑、台式机和高性能计算集群。