An optimisation based algorithm for finding the nucleation temperature of cosmological phase transitions

DOI10.1016/J.CPC.2021.108252arXiv2103.01985OpenAlexW3217034304MaRDI QIDQ6159604FDOQ6159604

Authors: Michael Bardsley 

Publication date: 20 June 2023

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: We present the OptiBounce algorithm, a new and fast method for finding the bounce action for cosmological phase transitions. This is done by direct solution of the "reduced" minimisation problem proposed by Coleman, Glaser, and Martin. By using a new formula for the action, our method avoids the rescaling step used in other algorithms based on this formulation. The bounce path is represented using a pseudo-spectral Gauss-Legendre collocation scheme leading to a non-linear optimisation problem over the collocation coefficients. Efficient solution of this problem is enabled by recent advances in automatic differentiation, sparse matrix representation and large scale non-linear programming. The algorithm is optimised for finding nucleation temperatures by sharing model initialisation work between instances of the calculation when operating at different temperatures. We present numerical results on a range of potentials with up to 20 scalar fields, demonstrating O(1%) agreement with existing codes and highly favourable performance characteristics.

Full work available at URL: https://arxiv.org/abs/2103.01985

zbMATH Keywords

phase transitionsgravitational wavesEuclidean actionvacuum stabilitybounce solutionbaryogensis

Mathematics Subject Classification ID

Numerical analysis (65-XX) Operations research, mathematical programming (90-XX)

Cites Work

• CasADi: a software framework for nonlinear optimization and optimal control
• On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming
• BubbleProfiler: finding the field profile and action for cosmological phase transitions
• Adaptive mesh refinement method for optimal control using nonsmoothness detection and mesh size reduction
• Direct trajectory optimization and costate estimation of finite-horizon and infinite-horizon optimal control problems using a Radau pseudospectral method
• Tunneling in theories with many fields
• Efficient numerical solution to vacuum decay with many fields
• A fresh look at the calculation of tunneling actions in multi-field potentials
• \texttt{FindBounce}: package for multi-field bounce actions
• \texttt{SimpleBounce}: a simple package for the false vacuum decay
• Uniqueness of a negative mode about a bounce solution

Cited In (2)

• Computing the gauge-invariant bubble nucleation rate in finite temperature effective field theory
• \texttt{BubbleDet}: a Python package to compute functional determinants for bubble nucleation

This page was built for publication: An optimisation based algorithm for finding the nucleation temperature of cosmological phase transitions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6159604)