[论文解读] Simulating the inflationary Universe: from single-field to the axion-U(1) model

Angelo Caravano|arXiv (Cornell University)|Sep 27, 2022

Computational Physics and Python Applications被引用 4

一句话总结

本论文提出了一套晶格模拟框架，用于研究单场及轴子-U(1)模型中的暴胀过程，通过在离散时空网格上求解非线性运动方程的数值解。主要贡献在于表明，由于反作用效应，曲率扰动中的非高斯性可能被强烈抑制，原初密度场的统计特性显示出超出三谱之外的累积量级联增长，揭示了标准n点函数描述的局限性。

ABSTRACT

We present a nonlinear study of the inflationary epoch based on numerical lattice simulations. Lattice simulations are a well-known tool in primordial cosmology, and they have been extensively used to study the reheating epoch after inflation. We generalize this known machinery to the inflationary epoch. Being this the first simulation of the inflationary epoch much before the end of inflation, the first part of the thesis focuses on the minimal single-field model of inflation. We discuss the conceptual and technical ingredients needed to simulate inflation on a lattice. The simulation is used to reproduce the nearly scale-invariant spectrum of scalar perturbations, as well as the oscillations in the power spectrum caused by a step in the potential. In the second part, we focus on the more complicated axion-U(1) model of inflation and present the first lattice simulation of this model during the deep inflationary epoch. We use the simulation to discover new properties of primordial scalar perturbations from this model. In the linear regime of the theory, we find high-order non-Gaussianity (beyond trispectrum) to be key to describing the statistical properties of scalar perturbations. Conversely, we find perturbations to be nearly Gaussian in the nonlinear regime of the theory. This relaxes existing constraints from the overproduction of primordial black holes, allowing for a gravitational waves signal in the observable range of upcoming experiments such as LISA. Our results show that lattice simulations can be a powerful tool to study the inflationary epoch and its observational signatures.

研究动机与目标

• 开发一种数值晶格框架，用于在强 gauge 场粒子产生模型中，超越微扰理论模拟暴胀动力学。
• 研究在轴子-U(1)暴胀模型中，反作用对标量和引力波扰动的影响。
• 分析曲率扰动 ζ 的完整非高斯统计特性，包括高阶累积量，以评估其对原初黑洞（PBH）形成的潜在影响。
• 将该方法推广至具有类似非微扰动力学但复杂度更高的非交换 SU(2) 规范场模型。
• 提供一种独立于标准n点函数形式体系的实空间、非微扰的暴胀前宇宙描述。

提出的方法

• 采用基于晶格的数值方法，在离散时空网格上求解标量场与规范场的非线性运动方程，使用半经典近似。
• 对空间导数采用有限差分离散化方案，并对有效动量进行细致处理，以考虑色散关系中的晶格效应。
• 时间演化采用四阶龙格-库塔（RK4）积分器，并通过监测能量守恒来评估数值精度。
• 初始条件通过标量场与规范场的量子真空涨落设定，并在U(1)理论中应用规范固定以保持一致性。
• 从晶格数据中计算统计输出，如功率谱、三谱以及 ζ 的完整概率密度函数（PDF）。
• 在慢滚和阶跃势能的单场暴胀模型上验证该方法后，将其应用于具有强规范场产生的轴子-U(1)模型。

实验结果

研究问题

• RQ1 在轴子-U(1)暴胀模型中，规范场产生引起的反作用如何影响曲率扰动的统计特性？
• RQ2 在强反作用存在下，标准n点函数（如功率谱、三谱）在多大程度上无法完整描述 ζ 的非高斯统计特性？
• RQ3 晶格模拟能否解析轴子-U(1)暴胀中规范场产生与能量守恒的非微扰动力学？
• RQ4 ζ 的完整PDF对本模型中原初黑洞（PBH）形成具有何种影响？
• RQ5 该晶格框架能否推广至具有类似动力学但复杂度更高的非交换 SU(2) 规范场？

主要发现

• 在强反作用的轴子-U(1)模型中，曲率扰动 ζ 的功率谱在大尺度被抑制，表明原初涨落受到显著阻尼。
• 由于反作用效应，ζ 的三谱被抑制，非高斯性水平低于当前来自宇宙微波背景（CMB）的观测约束。
• ζ 的高阶统计特性揭示出累积量级联增长（κ₅ > κ₄ > κ₃ > 1），表明标准n点函数无法捕捉 ζ 的完整非高斯形状。
• 模拟结果表明，ζ 的完整概率密度函数（PDF）是非高斯且不对称的，具有长尾特征，这对PBH形成至关重要，但传统微扰方法无法获取。
• 在整个模拟过程中能量守恒保持在1%以内，验证了RK4积分器与离散化方案的数值精度。
• 晶格框架成功捕捉了规范场产生与反作用的非微扰动力学，为早期宇宙提供了实空间、非微扰的描述。