TransitionListener v2.0 -- Robust gravitational wave predictions for cosmological phase transitions

This paper introduces TransitionListener v2.0, a robust Python framework that provides an end-to-end pipeline for precision gravitational wave predictions from cosmological phase transitions by incorporating self-consistent transition dynamics, improved physical modeling for challenging regimes, and comprehensive tools for detector sensitivity analysis and Bayesian inference.

The Gist

Imagine the early universe as a giant, super-hot pot of soup. As this soup cools down, it doesn't just get colder; it changes its state, much like water turning into ice. In physics, this is called a phase transition.

The paper introduces a new, upgraded software tool called TransitionListener v2.0. Think of this tool as a high-tech, super-accurate weather forecast for the universe's past. Its job is to predict what happened during these cosmic "freezing" events and, more importantly, to calculate the "sound" they would have made.

Here is a breakdown of what the paper says, using simple analogies:

1. The Cosmic "Pop" (Gravitational Waves)

When the universe cooled, it didn't always freeze smoothly. Sometimes, it got stuck in a "false" state (like supercooled water that hasn't frozen yet). Suddenly, bubbles of the "true" state would form and expand, crashing into each other like bubbles in a boiling pot.

When these bubbles collide and the surrounding fluid sloshes around, it creates ripples in space and time called Gravitational Waves. The paper explains that TransitionListener v2.0 is designed to predict exactly how loud these ripples would be and what frequency they would have, so scientists can know which telescopes (like LISA or the Einstein Telescope) might hear them.

2. The Problem with Old Maps

Before this new version, scientists used older tools to predict these events. The authors say these old tools were like using a paper map to navigate a stormy ocean. They worked okay for calm waters (weak transitions), but they broke down when things got wild (strong, super-cooled transitions).

The old tools made simplifying assumptions, like:

• "The universe expands at a steady, predictable rate."
• "The temperature drops in a simple, straight line."
• "The bubbles form and merge instantly."

The paper argues that in the most interesting, energetic scenarios, these assumptions are wrong. If you use the old map, you might predict a loud roar when the event was actually a whisper, or vice versa.

3. What's New in Version 2.0?

TransitionListener v2.0 is like upgrading from a paper map to a real-time GPS with a live camera feed. It fixes the old problems by:

• Self-Correction: It realizes that as bubbles form, they release energy that actually changes how the universe expands. The new tool calculates this "backreaction" on the fly, rather than ignoring it.
• The "Bubble Count": Instead of guessing how big the bubbles are, it counts them and measures the average distance between them directly from the simulation. This gives a much more accurate size for the "ripples" they create.
• Reheating the Soup: When bubbles collide, they release heat. The new tool calculates exactly how much the universe "reheats" during this process, which changes the final sound of the waves.
• Handling the Slow-Motion: Some transitions happen so slowly that the universe gets stuck. The new tool can handle these "ultraslow" scenarios where older tools would just crash or give up.

4. The "Recipe Book"

The software is built as a flexible pipeline. You can feed it a "recipe" (a specific theory of physics with new particles), and it will:

1. Trace the path: Figure out how the universe cooled and where the bubbles formed.
2. Simulate the crash: Calculate how the bubbles expanded and collided.
3. Predict the sound: Generate the gravitational wave signal.
4. Check the forecast: Compare that signal against the sensitivity of real-world detectors to see if we could actually hear it today.

5. Why It Matters

The authors tested their new tool against an existing popular tool (BSMPT). They found that for standard, mild transitions, both tools agreed. However, for the strongest, most energetic transitions (the ones most likely to be detected by future telescopes), the old tool gave very different, and likely inaccurate, results.

The paper concludes that TransitionListener v2.0 provides a much more reliable way to connect theoretical physics models with the real data that upcoming gravitational wave observatories will collect. It ensures that when we finally "hear" the universe's past, we know exactly what story that sound is telling us.

In short: This paper presents a smarter, more robust calculator that helps physicists predict the "soundtrack" of the early universe, ensuring that when we build telescopes to listen to the cosmos, we know exactly what to listen for.

Technical Summary

Technical Summary: TransitionListener v2.0 – Robust gravitational wave predictions for cosmological phase transitions

Problem Statement
The detection of gravitational wave (GW) backgrounds from strong first-order cosmological phase transitions (FOPTs) is a primary target for current and future observatories, including LISA, the Einstein Telescope, and Pulsar Timing Arrays (PTAs). However, translating specific Beyond-the-Standard-Model (BSM) theories into precise, testable GW predictions remains a significant challenge. Existing numerical tools often rely on approximations that break down in the regime of strongest GW signals: strongly supercooled and ultraslow transitions. Key limitations in previous codes include inconsistent treatments of thermodynamics during supercooling, simplified assumptions about the time-temperature relation (e.g., radiation domination), and the use of proxy scaling for reheating temperatures. Furthermore, handling multi-field potentials and seamlessly integrating physical calculations with observational likelihoods for large-scale parameter scans has been difficult.

Methodology
The paper introduces TransitionListener v2.0, a modular Python framework designed to provide an end-to-end pipeline from a user-defined scalar potential to GW spectra and signal-to-noise ratios (SNR). The methodology is built on a self-consistent treatment of the transition dynamics:

1. Effective Potential and Phase Tracing: The code constructs the 1-loop finite-temperature effective potential, including Coleman-Weinberg corrections, thermal functions, and daisy resummation. It employs a phase tracing algorithm to identify vacuum minima and possible transitions across a temperature range.
2. Nucleation and Percolation: Instead of relying on fixed nucleation criteria (e.g., $S_3/T \approx 140$), the code solves the percolation integral iteratively. It computes the evolution of the true-vacuum fraction ($P_t$) and its backreaction on the Hubble expansion rate ($H$). This involves solving a system of coupled ordinary differential equations (ODEs) for the percolation integral and the energy density of the true vacuum, accounting for local energy conservation at the bubble wall.
3. Thermodynamics and Hydrodynamics: The framework calculates thermodynamic quantities (pressure, energy density, entropy) directly from the full effective potential, including field-independent radiation terms. It determines the bubble wall velocity ($v_w$) using the Local Thermal Equilibrium (LTE) approximation, which serves as an upper bound, and computes efficiency factors for bubble collisions, sound waves, and turbulence based on the transition strength ($\alpha$) and wall velocity.
4. GW Spectrum and Observability: The code maps the nucleation history to GW spectral templates (collisions, sound waves, turbulence) using the mean bubble separation ($R_{sep}$) as the primary length scale, rather than the inverse duration parameter $(\beta/H)_{S3}$, which can be ill-defined for slow transitions. It evaluates observability by computing SNRs against built-in sensitivity curves for various detectors and interfaces with PTArcade for PTA likelihoods.
5. Scanning and Inference: The framework supports single-point studies, line scans, grid scans, random scans, and Bayesian inference via UltraNest, allowing for high-dimensional parameter space exploration.

Key Contributions

• Self-Consistent Percolation: The code implements a self-consistent solution for the true-vacuum fraction evolution, including its backreaction on the Hubble rate and the reheating process during percolation. This avoids the unphysical assumption that the universe remains radiation-dominated during the transition.
• Mean Bubble Separation ($R_{sep}$): The framework adopts the direct computation of mean bubble separation as the default length scale for GW templates. This provides a more faithful mapping to simulation-calibrated spectra, particularly in regimes where the conventional transition speed parameter $(\beta/H)_{S3}$ becomes negative or physically meaningless (ultraslow transitions).
• Improved Thermodynamics: Reheating temperatures are computed via local energy conservation rather than the approximation $T_{reh} \approx T_{perc}(1+\alpha)^{1/4}$. The code also correctly handles field-independent radiation terms in the energy density, which is crucial for strongly supercooled transitions.
• Multi-Field and Model Flexibility: The framework supports multi-field effective potentials with automated one-loop corrections and includes built-in models such as the Dark Abelian Higgs model, Conformal dark U(1)′, Dark flipflop, and the Two-Higgs Doublet Model (2HDM).
• Robust Error Handling: A detailed error-code system identifies non-physical scenarios (e.g., vacuum trapping, eternal inflation) and numerical instabilities, distinguishing them from physical failures.

Results and Validation
The authors validate TransitionListener v2.0 through benchmark points across different energy scales and by comparing it with the established code BSMPT v3.

• Benchmark Points: The code successfully reproduces GW spectra for four benchmark points, ranging from sub-GeV dark sector transitions (potentially detectable by $\mu$Ares) to electroweak-scale transitions (LISA) and high-scale dark Higgs transitions (Einstein Telescope).
• Comparison with BSMPT: For intermediate transition strengths, TransitionListener shows excellent agreement with BSMPT. However, in the regime of strong supercooling ($\alpha \gg 1$, $\beta/H \lesssim 100$), significant discrepancies arise. TransitionListener predicts different reheating temperatures (up to 20% difference) and transition strengths due to its more rigorous treatment of thermodynamics and the absence of certain bugs in BSMPT's daisy potential derivative.
• Numerical Stability: The code demonstrates improved numerical stability in the ultraslow regime where conventional approximations fail. It correctly identifies cases of vacuum trapping where the transition does not complete, whereas other codes may fail to converge or produce unphysical results.
• GW Signal Predictions: The study highlights that approximations used in previous tools can lead to $O(1)$ uncertainties in predicted SNRs for strong transitions, directly affecting the assessment of observability.

Significance
TransitionListener v2.0 addresses critical gaps in the physical consistency and numerical stability of GW predictions for cosmological phase transitions. By moving beyond standard approximations that break down in the most promising GW signal regimes (strong supercooling), the framework provides a robust tool for connecting fundamental BSM theories to observational data. The authors claim that this tool is essential for the upcoming era of GW astronomy, enabling both precision phenomenological studies and large-scale Bayesian inference required to interpret signals from LISA, PTAs, and future detectors. The code is released as open-source software (GPLv3) to facilitate reproducibility and further development in the field.