SimpleBounce: A Simple Package for False Vacuum Decay

Basic Information

• Paper ID: 1908.10868
• Title: SimpleBounce: A Simple Package for False Vacuum Decay
• Author: Ryosuke Sato (Deutsches Elektronen-Synchrotron (DESY))
• Classification: hep-ph astro-ph.CO hep-th
• Publication Date: DESY 19-148 (2019)
• Paper Link: https://arxiv.org/abs/1908.10868

Abstract

This paper introduces SimpleBounce, a C++ software package for computing bounce solutions in false vacuum decay. The package is based on gradient flow equations proposed by the author and solves the Coleman-Glaser-Martin minimization problem: minimizing kinetic energy while fixing the potential energy. Bounce configurations are obtained through scaling transformations of solutions to this problem. For models containing 1-8 scalar fields, the package can compute the bounce action with O(0.1)% precision in approximately 0.1 seconds.

Research Background and Motivation

Importance of the Problem

False vacuum decay is an important phenomenon in particle physics and cosmology. The lifetime of a false vacuum can be calculated through the Euclidean path integral formalism, where the bounce solution serves as a saddle point of the action, providing the leading contribution to the decay width.

Limitations of Existing Methods

Several numerical software packages for computing bounce solutions already exist, including:

• CosmoTransitions
• AnyBubble
• BubbleProfiler

However, these methods still have room for improvement in computational efficiency and numerical stability.

Research Motivation

The author proposes a new framework based on gradient flow equations, recently introduced by Chigusa, Moroi, and Shoji. This paper implements the specific flow equations proposed by the author in reference 1, aiming to provide a fast, stable, and user-friendly tool for computing bounce solutions.

Core Contributions

1. Development of SimpleBounce Package: A C++ implementation based on gradient flow equations for computing bounce solutions in false vacuum decay
2. Implementation of Numerical Solution for Coleman-Glaser-Martin Minimization Problem: Obtaining bounce solutions by minimizing kinetic energy while fixing potential energy
3. Provision of Efficient Computational Performance: For models with 1-8 scalar fields, achieving O(0.1)% precision in O(0.1) seconds
4. Verification of Method Accuracy: Validating results through comparison with existing tools such as CosmoTransitions

Methodology Details

Problem Definition

Given an Euclidean action containing nφ scalar fields: $S_E[\phi] = T[\phi] + V[\phi]$

where: $T[\phi] = \sum_{i=1}^{n_\phi} \int d^d x \frac{1}{2}(\nabla\phi_i)^2$$V[\phi] = \int d^d x [V(\phi) - V(\phi_{FV})]$

The objective is to find bounce solutions satisfying the equations of motion and boundary conditions.

Core Methodological Framework

1. Coleman-Glaser-Martin Minimization Problem

Instead of directly solving the equations of motion: $-\nabla^2\phi_i + \frac{\partial V}{\partial \phi_i} = 0, \quad \lim_{|x|\to\infty}\phi_i(x) = \phi_{FV,i}$

the method solves the minimization problem for kinetic energy Tφ while fixing negative potential energy Vφ < 0.

2. Gradient Flow Equations

Introducing the function φᵢ(r,τ), whose flow is described by: $\frac{\partial}{\partial \tau}\phi_i(r,\tau) = \nabla^2\phi_i - \lambda[\phi]\frac{\partial V(\phi)}{\partial \phi_i}$

where: $\lambda[\phi] = \frac{\sum_{i=1}^{n_\phi}\int_0^\infty dr r^{d-1}\frac{\partial V(\phi)}{\partial \phi_i}\nabla^2\phi_i}{\sum_{i=1}^{n_\phi}\int_0^\infty dr r^{d-1}\left(\frac{\partial V(\phi)}{\partial \phi_i}\right)^2}$

3. Scaling Transformation

Bounce solutions are obtained through scaling transformations of fixed points: $\phi_B(r) = \lim_{\tau\to\infty}\phi(\lambda^{-1/2}r, \tau)$

Technical Innovations

1. Numerical Stability: Since the potential energy V is fixed, the fixed point is not a saddle point, ensuring numerical stability
2. Adaptive Algorithm: Adjusting initial condition parameters r₀ and σ to ensure boundary conditions are satisfied
3. Efficient Discretization: Using specialized Laplacian operator discretization schemes suitable for spherically symmetric problems

Experimental Setup

Test Models

The paper employs multiple benchmark models for testing:

• Single scalar field model: V(φ) = ½φ² - ⅓φ³
• Eight benchmark models from reference 10
• Four benchmark models from reference 24

Evaluation Metrics

• Euclidean Action SE: The Euclidean action of the bounce solution
• Computation Time: Algorithm runtime
• Precision: Relative error compared to reference results

Comparison Methods

Primary comparison with CosmoTransitions, the widely-used standard tool in this field.

Implementation Details

• Spatial dimension: d = 3 or 4
• Grid points: n = 100-800
• Time step: δτ ~ O(δr²) to ensure numerical stability
• Boundary conditions: φᵢ,ₙ = φFV,i

Experimental Results

Main Results

Table 1: Comparison with CosmoTransitions

Results demonstrate that SimpleBounce's computational results are highly consistent with CosmoTransitions, and in most cases exhibit faster computation speeds.

Convergence Analysis

Figure 2 demonstrates the effects of different grid numbers n and flow time τ₁ on results:

• SE values converge as n and τ₁ increase
• Using n=100 achieves O(0.1)% precision in O(0.1) seconds
• Runtime scales as n³

Algorithm Characteristics

• Stability with respect to initial configuration selection
• Thin-wall bounces require more computation time than thick-wall bounces
• Computation time scales approximately linearly with the number of scalar fields

Related Work

Existing Numerical Methods

1. CosmoTransitions: Based on finite difference methods for solving boundary value problems
2. AnyBubble: Efficient numerical solver specialized for multi-field cases
3. BubbleProfiler: Bubble wall calculation tool combining semi-analytical techniques

Theoretical Foundations

• Coleman-Glaser-Martin Theory: Establishes connection between bounce solutions and minimization problems
• Gradient Flow Methods: New framework proposed by Chigusa et al.
• O(d) Symmetry: Spherical symmetry characteristics of bounce solutions

Conclusions and Discussion

Main Conclusions

1. SimpleBounce provides a fast and accurate method for computing bounce solutions
2. The gradient flow-based method exhibits good numerical stability
3. The software package demonstrates superior computational efficiency compared to existing tools

Limitations

1. Thin-Wall Limit: Computation time increases significantly for thin-wall bounces
2. Grid Dependence: Higher precision requirements necessitate finer grids, increasing computational cost
3. Scalar Field Limitation: Current testing primarily focuses on 1-8 scalar fields

Future Directions

1. Optimize computational efficiency for thin-wall bounces
2. Extend to cases with more scalar fields
3. Integration into larger physics computation frameworks

In-Depth Evaluation

Strengths

1. Methodological Innovation: The gradient flow approach provides new insights for numerical solution methods
2. High Practical Value: Open-source software package facilitates community use and verification
3. Excellent Performance: Significantly improves computational efficiency while maintaining precision
4. Solid Theoretical Foundation: Built upon Coleman-Glaser-Martin theory

Weaknesses

1. Limited Testing Scope: Testing primarily focuses on relatively simple models
2. Thin-Wall Performance: Efficiency for thin-wall bounce handling requires further improvement
3. Documentation Completeness: As a software package, user documentation could be more comprehensive

Impact

1. Academic Contribution: Provides new tools for false vacuum decay calculations
2. Practical Value: Likely to become a standard computational tool in this field
3. Reproducibility: Open-source code ensures result reproducibility

Applicable Scenarios

• Vacuum stability analysis in particle physics
• Cosmological phase transition research
• Numerical calculations in scalar field theory
• Theoretical research requiring rapid bounce solution computation

References

Primary references include classical theoretical works by Coleman et al., papers on existing numerical methods, and the author's previous theoretical research on gradient flow equations. The complete software package implementation is available on GitHub.