We study first-order electroweak phase transitions in the real-singlet-extended Standard Model, for which nonzero mixing between the Higgs field and the singlet can efficiently strengthen the transitions. We perform large-scale parameter-space scans of the model using two-loop effective potential at next-to-next-to-leading order in the high-temperature expansion, greatly improving description of phase transition thermodynamics over existing one-loop studies. We find that (1) two-loop corrections to the effective potential lead to narrower regions of strong first-order transitions and significantly smaller critical temperatures, (2) transitions involving a discontinuity in the singlet expectation value are significantly stronger at two-loop order, and (3) high-temperature expansion is accurate for a wide range of parameter space that allows strong transitions, although it is less reliable for the very strongest transitions. These findings suggest revisiting past studies that connect the possibility of a first-order electroweak phase transition with future collider phenomenology.

We complete the perturbative program for equilibrium thermodynamics of cosmological first-order phase transitions of gauge-Higgs theories that map into the (three-dimensional) superrenormalizable SU(2)+doublet effective theory at high temperatures. To this end, we determine their finite-temperature effective potential at next-to-next-to-next-to-next-to-leading order (N4LO). The computation of the three-loop effective potential required to reach this order is also presented for U(1) gauge theories and is readily extendable to generic models in dimensionally reduced effective theories. Our N4LO result is the last perturbative order before confinement renders electroweak gauge-Higgs theories nonperturbative at four loops. By contrasting our analysis with nonperturbative lattice results, we find a remarkable agreement. As a direct application for predictions of gravitational waves produced by a first-order transition, our computation provides the final fully perturbative results for the phase transition strength and speed of sound.

Models with radiative symmetry breaking typically feature strongly supercooled first-order phase transitions, which result in an observable stochastic gravitational wave background. In this work, we analyse the role of higher-order thermal corrections for these transitions, applying high-temperature dimensional reduction to a theory with dimensional transmutation. In particular, we study to what extent high-temperature effective field theories (3D EFT) can be used. We find that despite significant supercooling down from the critical temperature, the high-temperature expansion for the bubble nucleation rate can be applied using the 3D EFT framework, and we point out challenges in the EFT description. We compare our findings to previous studies and find that the next-to-leading order corrections obtained in this work have a significant effect on the predictions for GW observables, motivating a further exploration of higher-order thermal effects.

Guided by previous non-perturbative lattice simulations of a two-step electroweak phase transition, we reformulate the perturbative analysis of equilibrium thermodynamics for generic cosmological phase transitions in terms of effective field theory (EFT) expansions. Based on thermal scale hierarchies, we argue that the scale of many interesting phase transitions is in-between the soft and ultrasoft energy scales, which have been the focus of studies utilising high-temperature dimensional reduction. The corresponding EFT expansions provide a handle to control the perturbative expansion, and allow us to avoid spurious infrared divergences, imaginary parts, gauge dependence and renormalisation scale dependence that have plagued previous studies. As a direct application, we present a novel approach to two-step electroweak phase transitions, by constructing separate effective descriptions for two consecutive transitions. Our approach provides simple expressions for effective potentials separately in different phases, a numerically inexpensive method to determine thermodynamics, and significantly improves agreement with the non-perturbative lattice simulations.

DRalgo is an algorithmic implementation that constructs an effective, dimensionally reduced, high-temperature field theory for generic models. The corresponding Mathematica package automatically performs the matching to next-to-leading order. This includes two-loop thermal corrections to scalar and Debye masses as well as one-loop thermal corrections to couplings. DRalgo also allows for integrating out additional heavy scalars. Along the way, the package provides leading-order beta functions for general gauge-charges and fermion-families; both in the fundamental and in the effective theory. Finally, the package computes the finite-temperature effective potential within the effective theory. The article explains the theory of the underlying algorithm while introducing the software on a pedagogical level.

We present a gauge-invariant framework for bubble nucleation in theories with radiative symmetry breaking at high temperature. As a procedure, this perturbative framework establishes a practical, gauge-invariant computation of the leading order nucleation rate, based on a consistent power counting in the high-temperature expansion. In model building and particle phenomenology, this framework has applications such as the computation of the bubble nucleation temperature and the rate for electroweak baryogenesis and gravitational wave signals from cosmic phase transitions.

For computing thermodynamics of the electroweak phase transition, we discuss a minimal approach that reconciles both gauge invariance and thermal resummation. Such a minimal setup consists of a two-loop dimensional reduction to three-dimensional effective theory, a one-loop computation of the effective potential and its expansion around the leading-order minima within the effective theory. This approach is tractable and provides formulae for resummation that are arguably no more complicated than those that appear in standard techniques ubiquitous in the literature. In particular, we implement renormalisation group improvement related to the hard thermal scale. Despite its generic nature, we present this approach for the complex singlet extension of the Standard Model which has interesting prospects for high energy collider phenomenology and dark matter predictions. The presented expressions can be used in future studies of phase transition thermodynamics and gravitational wave production in this model.

A gauge-invariant framework for computing bubble nucleation rates at finite temperature in the presence of radiative barriers was presented and advocated for model-building and phenomenological studies in an accompanying article [1]. Here, we detail this computation using the Abelian Higgs Model as an illustrative example. Subsequently, we recast this approach in the dimensionally-reduced high-temperature effective field theory for nucleation. This allows for including several higher order thermal resummations and furthermore delineate clearly the approach’s limits of validity. This approach provides for robust perturbative treatments of bubble nucleation during possible first-order cosmic phase transitions, with implications for electroweak baryogenesis and production of a stochastic gravitational wave background. Furthermore, it yields a sound comparison between results of perturbative and non-perturbative computations. 

The energy budget for gravitational waves of a cosmological first order phase transitions depends on the speed of sound in the thermal plasma in both phases around the bubble wall. Working in the real-singlet augmented Standard Model, which admits a strong two-step electroweak phase transition, we compute higher order corrections to the pressure and sound speed. We compare our result to lower-order approximations to the sound speed and the energy budget and investigate the impact on the gravitational wave signal. We find that deviations in the speed of sound from c = 1/3 are enhanced up to O(5%) in our higher-order computation. This results in a suppression in the energy budget of up to O(50%) compared to approximations assuming c = 1/3. The effect is most significant for hybrid and detonation solutions. We generalise our discussion to the case of multiple inert scalars and the case of a reduced number of fermion families in order to mimic hypothetical dark sector phase transitions. In this sector with modified field content, the sound speed can receive significant suppression, with potential order-of-magnitude impact on the gravitational wave amplitude.

Beyond the Standard Model physics is required to explain both dark matter and the baryon asymmetry of the universe, the latter possibly generated during a strong first-order electroweak phase transition. While many proposed models tackle these problems independently, it is interesting to inquire whether the same model can explain both. In this context, we link state-of-the-art perturbative assessments of the phase transition thermodynamics with the extraction of the dark matter energy density. These techniques are applied to a next-to-minimal dark matter model containing an inert Majorana fermion that is coupled to Standard Model leptons via a scalar mediator, where the mediator interacts directly with the Higgs boson. For dark matter masses 180 GeV < M_χ ≲ 300 GeV, we discern regions of the model parameter space that reproduce the observed dark matter energy density and allow for a first-order phase transition, while evading the most stringent collider constraints.