PT2GWFinder: A Package for Cosmological First-Order Phase Transitions and Gravitational Waves

This paper introduces PT2GWFinder, a Mathematica package that automates the computation of phase transition parameters and gravitational wave spectra for arbitrary scalar theories exhibiting first-order phase transitions, featuring integration with DRalgo for high-temperature effective potentials and validated through analytical and numerical case studies.

The Gist

Imagine the universe as a giant, cosmic pot of water. When you heat water, it's liquid. When you cool it down, it freezes into ice. But what if, instead of freezing smoothly, the water suddenly decided to snap into ice in a chaotic, explosive way? In the early universe, something similar happened, but with the fundamental building blocks of reality instead of water.

This paper introduces a new digital tool called PT2GWFinder (think of it as a "Phase Transition Treasure Finder") designed to help scientists predict the echoes of these cosmic explosions.

Here is the story of the paper, broken down into simple concepts:

1. The Big Bang's "Snap" (First-Order Phase Transitions)

In the very early universe, things were incredibly hot. As the universe expanded and cooled, it went through "phase transitions," much like water turning to ice.

• The Smooth Way: Usually, these changes happen gradually (like water slowly getting colder).
• The "Snap" Way: Sometimes, the universe gets stuck in a "false vacuum" (a metastable state, like supercooled water that hasn't frozen yet). Suddenly, it snaps into a "true vacuum" (the stable ice state).

This snap doesn't happen everywhere at once. Instead, tiny bubbles of the new "true" reality start popping into existence, like bubbles in boiling water. These bubbles expand rapidly, smash into each other, and create ripples in the fabric of space and time. These ripples are Gravitational Waves.

2. The Problem: It's Hard to Predict the Ripples

Scientists know these bubbles should exist if certain theories about the universe are true. But calculating exactly how loud the "crash" of these bubbles would be is incredibly difficult. It involves complex math, simulating how bubbles grow, how they collide, and how they shake the cosmic fluid.

It's like trying to predict the exact sound of a thousand balloons popping in a hurricane. You need a super-computer and a very specific set of instructions to get the answer right.

3. The Solution: PT2GWFinder

This is where the new package comes in. PT2GWFinder is a software tool (written for a program called Mathematica) that acts as an automated cosmic weather forecaster.

• The Input: You give the software a mathematical recipe for a specific universe model (the "potential").
• The Process:
  1. Phase Tracing: It looks at the recipe and figures out when the universe would get stuck in the "false vacuum" and when it would snap.
  2. Bubble Simulation: It calculates how fast the bubbles would grow and how hard they would hit each other.
  3. The Echo: It predicts the specific "song" (the gravitational wave spectrum) that these collisions would create today.
• The Output: It tells you: "If this universe model is real, we should hear a gravitational wave signal at this specific frequency with this specific loudness."

4. Why Do We Care? (The Detective Work)

We have detectors on Earth (like LIGO) and in space (like the future LISA mission) listening for these gravitational waves.

• If we hear a signal that matches the "song" predicted by PT2GWFinder, it's proof that new physics exists beyond our current understanding of the universe.
• It's like finding a fingerprint at a crime scene. If the fingerprint matches the suspect (the new physics model), we know who did it.

5. How the Tool Works (The "Magic" Inside)

The paper explains that the tool uses some clever tricks to do the heavy lifting:

• The "Bounce" Solver: It calculates the path a bubble takes when it tunnels through the energy barrier to form. It's like calculating the exact trajectory of a ball rolling over a hill.
• The "Drain" Helper: It can connect with other advanced tools (like DRalgo) to handle complex math that simplifies the universe's behavior at high temperatures.
• The "Efficiency" Calculator: It figures out how much of the explosion's energy turns into sound waves (ripples) versus how much is wasted.

6. The Test Drive

The authors didn't just build the tool; they tested it.

• They used it on two known models: one simple "toy" model and a more complex "Dark Abelian Higgs" model (which involves invisible, dark matter particles).
• They compared their results with other existing tools and found that PT2GWFinder was just as accurate, but often faster and easier to use.
• They even proved that for simple cases, their computer calculations matched the "perfect" math formulas exactly.

Summary

PT2GWFinder is a new, user-friendly calculator for cosmologists. It takes a theoretical model of the early universe, simulates the violent "bubble collisions" that happened billions of years ago, and predicts the gravitational wave signal we might hear today.

It's a bridge between abstract math (what the universe could be) and real-world data (what our telescopes might hear), helping us solve the mystery of what happened in the first split-second of the universe's life.

Technical Summary

1. Problem Statement

The detection of gravitational waves (GWs) has opened a new window into the early universe. One of the most promising sources of a stochastic gravitational wave background (SGWB) is a cosmological first-order phase transition (FOPT), which occurs when the universe decays from a metastable false vacuum to a stable true vacuum via bubble nucleation.

However, predicting the GW signal from a specific Beyond-the-Standard-Model (BSM) scenario is computationally complex. It requires:

• Constructing the finite-temperature effective potential (handling infrared divergences via daisy resummation or dimensional reduction).
• Tracing the evolution of vacuum phases as a function of temperature.
• Solving the Euclidean equations of motion to find the "bounce" solution and calculating the Euclidean action ($S_3$).
• Determining critical temperatures ($T_c$), nucleation temperatures ($T_n$), and percolation temperatures ($T_p$).
• Calculating GW efficiency factors and the resulting power spectrum from bubble collisions, sound waves, and turbulence.

Existing tools (e.g., BSMPTv3, PhaseTracer2, ELENA) often have limitations in flexibility, speed, or the integration of advanced methods like dimensional reduction. There is a need for a robust, user-friendly, and comprehensive pipeline that connects a particle physics model directly to the observable GW spectrum.

2. Methodology

The authors introduce PT2GWFinder, a Mathematica package designed to automate the entire pipeline from a scalar field theory to the GW spectrum. The methodology is structured as follows:

• Input: The package requires the user to supply the effective thermal potential ($V(\phi, T)$) and the bubble wall velocity ($v_w$). It is agnostic to how the potential was constructed (e.g., via daisy resummation or dimensional reduction).
• Phase Tracing: The package identifies the symmetric and broken phases of the potential across a temperature range. It uses numerical minimization (FindMinimum) or semi-analytical methods (NSolve) to locate minima and determine the critical temperature ($T_c$) where phases are degenerate.
• Bounce Solution & Action:
  • Utilizes the FindBounce package, which employs the polygonal bounce method. This approximates the potential as linear segments between vacua to analytically solve the bounce profile on each segment.
  • Implements numerical stability strategies, such as exit point shifting (optimizing the segmentation range when the exit point is far from the true vacuum) and action filtering/refining to handle numerical instabilities near inflection points.
  • Fits the Euclidean action $S_3(T)$ using a piecewise function: a cubic spline at lower temperatures and a Laurent polynomial at high temperatures (near $T_c$).
• Temperature Determination:
  • Nucleation Temperature ($T_n$): Calculated using the integral criterion $\int \frac{\Gamma}{H^4} dT \approx 1$ or the approximate condition $S_3/T \approx 140$.
  • Percolation Temperature ($T_p$): Determined by solving the JMAK equation where the false vacuum fraction drops to $\approx 0.71$ (34% of volume converted).
• Gravitational Wave Spectrum:
  • Computes the GW power spectrum from three sources: bubble collisions, sound waves, and MHD turbulence.
  • Calculates efficiency factors ($\kappa_{col}, \kappa_{sw}$) based on the energy released and the bubble wall velocity.
  • Uses templates from the LISA Cosmology Working Group and includes sensitivity curves for current and future detectors (LISA, DECIGO, BBO, etc.).
• Dimensional Reduction Integration: Includes a helper module, DRTools, to interface with the DRalgo package. This allows users to import dimensionally reduced effective potentials (3D EFT) directly into PT2GWFinder, facilitating high-temperature calculations.

3. Key Contributions

• Comprehensive Pipeline: PT2GWFinder is the first Mathematica package to provide a unified, end-to-end workflow for single-field FOPTs, covering phase tracing, bounce computation, temperature determination, and GW spectrum generation.
• Advanced Numerical Techniques:
  • Implementation of exit point shifting to improve the accuracy of the bounce solution in regions with low potential barriers.
  • Robust action fitting combining splines and Laurent polynomials to accurately capture the divergence of $S_3/T$ near $T_c$.
  • Analytical control via thin-wall approximation for polynomial potentials, allowing for benchmarking and faster estimates.
• DRTools Helper: A dedicated module to seamlessly integrate the output of DRalgo (dimensional reduction) into the main package, supporting high-precision thermal field theory calculations.
• User-Friendly Interface: The package features an object-oriented design (Transition objects) with autocompletion, easy plotting functions, and the ability to store metadata. It supports both pre-loaded benchmark models and custom user-defined potentials.
• Validation: The package was rigorously tested against literature results and other codes (CosmoTransitions, BSMPTv3) using two distinct models.

4. Results and Validation

The authors validated PT2GWFinder using two phenomenologically relevant models:

1. Coupled Fluid-Scalar Field (CFF) Model:

   • A simple quartic polynomial potential widely used for numerical simulations.
   • Result: The package reproduced results from Ref. [59] with high precision.
   • Analytical Agreement: For this model, the authors derived the bounce action and nucleation temperature analytically using the thin-wall approximation. The numerical results from PT2GWFinder showed excellent agreement with these analytical solutions across the entire temperature range.
   • Performance: Analyzed benchmark points in ~45 seconds.

2. Dark Abelian Higgs Model:

   • A complex scalar charged under a dark $U(1)'$ gauge group.
   • Comparison 1 (Daisy Resummation): Results were compared with CosmoTransitions. The phase transition strength ($\alpha$) and inverse duration ($\beta/H$) matched well, though PT2GWFinder showed slightly better stability in fitting the action near the inflection point.
   • Comparison 2 (Dimensional Reduction): The authors used DRalgo to generate the effective potential and DRTools to process it. They found that the dimensional reduction approach predicted a stronger and slower transition compared to the daisy resummation approach for the same parameters, highlighting the importance of the method used to construct the potential.
   • Performance: The full scan of the Dark Abelian Higgs parameter space took approximately 6 minutes and 40 seconds using parallelization.

5. Significance

• Accessibility: By providing a high-level, modular Mathematica package, PT2GWFinder lowers the barrier to entry for researchers studying cosmological phase transitions, allowing them to focus on model building rather than numerical implementation details.
• Reliability: The inclusion of advanced numerical fixes (exit point shifting, action filtering) ensures that the package can handle "difficult" potentials where other codes might fail, particularly near inflection points or in the thin-wall regime.
• Future-Proofing: The modular architecture allows for easy integration of new GW templates, efficiency factor calculations, or multi-field extensions (planned for future versions).
• Phenomenological Impact: As GW detectors (LISA, DECIGO, BBO) come online, PT2GWFinder provides a critical tool for interpreting potential signals, helping to constrain BSM physics and understand the thermal history of the universe.

In summary, PT2GWFinder represents a significant advancement in computational cosmology, offering a robust, flexible, and validated tool for connecting theoretical particle physics models to observable gravitational wave signatures.