[论文解读] How fast does the WallGo? A package for computing wall velocities in first-order phase transitions

Andreas Ekstedt, Oliver Gould+5 more|arXiv (Cornell University)|Jan 1, 2024

Advanced Thermodynamics and Statistical Mechanics被引用 2

一句话总结

WallGo 是一个开源的 Python 软件包，通过求解耦合的标量场方程、流体流体动力学方程和玻尔兹曼方程，以高精度计算一阶宇宙相变中的气泡壁速度。它采用谱方法在基多项式中实现指数收敛，并包含专门的模块用于矩阵元（WallGoMatrix）和碰撞积分（WallGoCollision），从而能够准确建模非平衡态粒子对阻尼力和引力波产生的贡献。

ABSTRACT

WallGo is an open source software for the computation of the bubble wall velocity in first-order cosmological phase transitions. It also computes the energy budget available for the generation of gravitational waves. The main part of WallGo, built in Python, determines the wall velocity by solving the scalar-field(s) equation of motion, the Boltzmann equations and energy-momentum conservation for the fluid velocity and temperature. WallGo also includes two auxiliary modules: WallGoMatrix, which computes matrix elements for out-of-equilibrium particles, and WallGoCollision, which performs higher-dimensional integrals for Boltzmann collision terms. Users can implement custom models by defining an effective potential and specifying a list of out-of-equilibrium particles and their interactions. As the first public software to compute the wall velocity including out-of-equilibrium contributions, WallGo improves the precision of the computation compared to common assumptions in earlier computations. It utilises a spectral method for the deviation from equilibrium and collision terms that provides exponential convergence in basis polynomials, and supports multiple out-of-equilibrium particles, allowing for Boltzmann mixing terms. WallGo is tailored for non-runaway wall scenarios where leading-order coupling effects dominate friction. While this work introduces the software and the underlying theory, a more detailed documentation can be found in https://wallgo.readthedocs.io.

研究动机与目标

• 开发一种精确、开源的工具，用于计算一阶相变中的终端壁速度，解决先前近似方法的局限性。
• 包含非平衡态粒子对阻尼力的贡献，提高文献中标准假设的准确性。
• 提供模块化、可扩展的框架，使用户可通过有效势和粒子相互作用实现自定义模型。
• 实现早期宇宙相变中可用于引力波产生的能量预算的精确计算。
• 支持以主导性耦合效应为特征的非发散壁情景，适用于现实模型构建。

提出的方法

• 使用变分法和谱基中的试探函数求解标量场运动方程。
• 在能量-动量守恒条件下求解流体流体动力学方程，以描述速度和温度的演化。
• 通过基多项式对玻尔兹曼方程进行谱分解，实现指数收敛。
• 使用 WallGoCollision 对高维碰撞项进行自适应蒙特卡洛积分。
• 利用 WallGoMatrix 计算非平衡过程中的矩阵元，包括费米子和标量相互作用。
• 支持多个非平衡态粒子，通过自旋旋量-手征性形式处理自旋依赖的振幅，并引入混合项。

实验结果

研究问题

• RQ1 当包含非平衡态粒子贡献时，一阶相变中壁速度的计算精度如何？
• RQ2 谱方法在求解玻尔兹曼方程和标量场方程时，对收敛速度和精度有何影响？
• RQ3 多粒子过程的碰撞积分在数值上如何缩放？何种积分策略能确保鲁棒性？
• RQ4 与标准平衡近似相比，非平衡效应在多大程度上改变了终端壁速度？
• RQ5 该框架能否推广至任意模型，支持用户自定义的有效势和粒子内容？

主要发现

• 由于采用谱方法，WallGo 在基多项式数量上实现指数收敛，显著提升了计算效率和精度。
• 通过谱玻尔兹曼解包含非平衡贡献，其壁速度预测精度高于常见的平衡近似。
• WallGoCollision 通过动量守恒狄拉克函数将九维碰撞积分简化为五维积分，实现鲁棒的数值评估。
• 代码库支持多个非平衡态粒子并包含混合项，可实现多组分离子等离子体中摩擦力的现实建模。
• 收敛性测试证实了在不同网格参数、平均自由程尺度以及场/温度变化尺度下的稳定性与精度。
• 该框架已具备生产就绪性且可扩展，通过用户自定义的有效势和相互作用拉格朗日量，提供完整文档支持。