⚛️ phenomenology

Strong First-Order Electroweak Phase Transitions and Gravitational Waves in the Normal Two-Higgs-Doublet Model: A Comparative Study of the Four Yukawa Types and Thermal Resummation Schemes

This paper presents a global analysis of strong first-order electroweak phase transitions and associated gravitational waves in the Normal Two-Higgs-Doublet Model, revealing that while both the Parwani and Arnold-Espinosa thermal resummation schemes favor the Higgs alignment limit, the Arnold-Espinosa prescription severely restricts the viable parameter space and heavy-scalar mass scales compared to the more stable Parwani approach, thereby highlighting the critical role of theoretical scheme choices in predicting observable signals for future experiments like LISA and colliders.

Original authors: Jin-Hwan Cho, Dongjoo Kim, Jinheung Kim, Soojin Lee, Jeonghyeon Song

Published 2026-07-14

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jin-Hwan Cho, Dongjoo Kim, Jinheung Kim, Soojin Lee, Jeonghyeon Song

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the early universe as a giant, super-hot pot of soup. As it cools down, it's supposed to undergo a dramatic "phase change," like water turning into ice. In our current understanding of physics, this change happens smoothly, like butter softening. But for the universe to have the matter we see today, this change needed to be a violent, explosive "first-order" transition—like water suddenly boiling into steam with a massive pop. This violent pop would create ripples in spacetime called gravitational waves, which we hope to catch with future space telescopes.

This paper is a massive, computer-simulated detective story about a specific theory called the "Two-Higgs-Doublet Model" (2HDM). Think of the Standard Model of physics as a recipe with one type of flour (one Higgs boson). The 2HDM suggests there's actually a second, hidden type of flour. The authors asked: "If we add this second flour, can we get that violent 'pop' we need, and will it make a sound loud enough for our future telescopes to hear?"

The Great Recipe Debate: Two Ways to Cook
The biggest surprise in this study isn't just about the ingredients, but about how you calculate the heat. The authors compared two different mathematical "recipes" for handling the heat in the early universe soup: the Parwani method and the Arnold–Espinosa (AE) method.

It turns out, the choice of recipe changes everything.

The Parwani Method: This approach suggests a much more generous playground. It says the "heavy" extra particles (the new Higgs bosons) could be as heavy as 1.6 TeV (that's 1,600 billion electron volts) and still cause the violent pop. The results are smooth and continuous, like a well-organized map.
The Arnold–Espinosa (AE) Method: This is the stricter, more cautious recipe. It suggests the heavy particles can't be nearly as heavy; they must be below about 800 GeV. Worse, when the authors used this method, the "good" spots on the map didn't look like a smooth region. Instead, they looked like scattered islands in a sea of empty space, with weird holes and gaps. The authors found this method so sensitive that tiny changes in the particle masses caused the results to jump around wildly, creating a fragmented mess.

Because the Parwani method seems more stable and matches other advanced theories better, the authors decided to focus their main story on the results from the Parwani "recipe."

The Four Flavors of the Model
The 2HDM comes in four different "flavors" (Type-I, Type-II, Type-X, and Type-Y), depending on how these new particles interact with other matter. The authors ran a massive simulation, checking 1 million possible settings for each of these four flavors.

Here is what they found:

One Step vs. Two Steps: In the "Normal Scenario" (where the lightest new particle is the one we already found at 125 GeV), the universe almost always takes a one-step path to the violent transition. It's rare for it to take a two-step path. This is the opposite of what happens in a different version of the theory (the "Inverted Scenario"), where two-step transitions are common.
The Weight Limit: For the violent transition to happen, the new heavy particles cannot be too heavy. The simulation shows a hard ceiling: they must be lighter than 1.6 TeV. If they are heavier, the transition becomes too smooth to be useful.
The "Alignment" Rule: The new particles must behave very much like the Higgs boson we already know. They have to be "aligned" with it. The authors found that the new particles can't deviate too much from the behavior of the known Higgs, or the violent transition won't happen.

The Sound of the Universe
Even if the universe has a violent transition, will it make a sound loud enough for the LISA telescope (a future space-based detector) to hear? The authors defined a "loud" signal as having a signal-to-noise ratio (SNR) of 10 or higher after four years of listening.

The Good News: Some settings do produce loud signals. In the "Type-I" flavor, the loudest signal could reach an SNR of 95. In "Type-X," it could reach 60.
The Bad News: The "Type-II" and "Type-Y" flavors are much quieter. In Type-II, the maximum signal was only 19, and in Type-Y, it was also 19. In fact, when using the stricter AE recipe, Type-II and Type-Y produced almost zero loud signals.
The Heavy Mass Cut: To get a loud signal, the heavy particles must be even lighter than the general limit. While the transition can happen with particles up to 1.6 TeV, to get a detectable sound, they usually need to be below 700 GeV.

The "Short-Lived" Sound
One of the most interesting findings is why the signal is sometimes weak. The authors discovered that the "sound" (acoustic waves in the early universe soup) doesn't last very long.
Usually, scientists assumed these waves would ring out for a long time, like a bell. But in these simulations, the waves crash into each other and turn into turbulence very quickly. The sound source dies out in just a tiny fraction of a second (less than 1% of the time it takes for the universe to expand significantly). This "short lifetime" acts like a volume knob turned way down, suppressing the signal.

The Verdict
The paper concludes that while the Two-Higgs-Doublet Model can explain the violent transition we need, it's a very picky theory.

It suggests that the new heavy particles must be relatively light (under 1.6 TeV, and under 700 GeV for a loud signal).
It suggests that the "Type-I" and "Type-X" flavors are the most likely to produce a signal we can hear, while "Type-II" and "Type-Y" are much harder to detect.
It warns us that our mathematical tools matter: using the wrong "recipe" (like the AE method) might make us think the theory is broken or fragmented when it might actually be fine.

Ultimately, the authors suggest that if we build the LISA telescope and start listening, and if we build bigger particle colliders to hunt for these heavy particles, we might finally catch a glimpse of this violent, ancient pop in the history of our universe. But we have to be careful about how we calculate the heat, or we might miss the sound entirely.

Technical Summary: Strong First-Order Electroweak Phase Transitions and Gravitational Waves in the Normal Two-Higgs-Doublet Model

Problem Statement
The origin of the baryon asymmetry of the Universe requires a departure from thermal equilibrium, typically achieved via a strong first-order electroweak phase transition (SFOEWPT). In the Standard Model (SM), the observed 125 GeV Higgs boson mass implies a smooth crossover, failing to satisfy this condition. Extended Higgs sectors, such as the Two-Higgs-Doublet Model (2HDM), offer a mechanism for SFOEWPTs and the generation of stochastic gravitational-wave (GW) backgrounds detectable by future space-based interferometers like LISA.

While previous studies have explored the "Inverted Scenario" of the 2HDM (where the heavier CP-even scalar is the 125 GeV Higgs), the "Normal Scenario" (where the lighter CP-even scalar $h$ is the 125 GeV Higgs) remains less comprehensively analyzed, particularly regarding the interplay between different Yukawa structures and thermal resummation schemes. Furthermore, the theoretical uncertainty associated with the treatment of thermal daisy diagrams in the finite-temperature effective potential has not been systematically quantified across the full parameter space of the Normal 2HDM.

Methodology
The authors perform a high-statistics global scan of the CP-conserving 2HDM with softly broken $Z_2$ symmetry under the Normal Scenario ( $m_h = 125$ GeV). The study covers all four Yukawa types (Type-I, II, X, and Y).

Parameter Space: The scan utilizes $10^6$ physically viable parameter points per Yukawa type, satisfying theoretical constraints (boundedness from below, vacuum stability, unitarity, perturbativity) and experimental limits (electroweak precision observables, flavor constraints like $B \to X_s\gamma$ , and LHC Higgs signal strengths).
Thermal Resummation: A comparative analysis is conducted between two thermal resummation prescriptions for daisy diagrams: the Parwani scheme (global mass replacement) and the Arnold–Espinosa (AE) scheme (selective resummation of Matsubara zero modes). The analysis is performed using BSMPT v3.1.8.
Classification: Parameter points are categorized into:
- SFOEWPT points: Configurations satisfying the strong transition criterion $\xi_p \equiv \Delta v(T_p)/T_p > 1$ .
- GW points: SFOEWPT configurations yielding a four-year LISA signal-to-noise ratio (SNR) $\ge 10$ .
Hydrodynamics: The GW signal is modeled including acoustic sound waves, bubble collisions, and turbulence. Special attention is paid to the sound-wave lifetime suppression factor $\Upsilon$ , which accounts for the premature termination of the acoustic source due to turbulence.

Key Contributions

Comprehensive Comparative Study: This is the first global study of SFOEWPTs and GW signatures in the Normal Scenario across all four Yukawa types, contrasting with prior work limited to the Inverted Scenario or single Yukawa types.
Thermal Resummation Uncertainty Quantification: The paper provides a systematic comparison of the Parwani and AE schemes, revealing that the choice of resummation prescription induces significant theoretical uncertainty in the viable parameter space and macroscopic phase-transition parameters.
Hydrodynamic Analysis: The authors investigate the plasma hydrodynamics underlying GW signals, specifically identifying the prevalence of the "short-lifetime" acoustic regime and its impact on signal amplitude.
High-Statistics Baseline: The analysis handles a total of 8 million distinct thermodynamic histories, establishing a statistically definitive baseline for the 2HDM phase-transition landscape.

Results

Resummation Scheme Dependence:
- Viable Parameter Space: The Parwani scheme yields a significantly larger number of SFOEWPT points ( $O(10^4)$ – $O(10^5)$ per type) compared to the AE scheme. The AE scheme drastically suppresses the viable space, particularly for Type-II and Type-Y, where the yield drops to $O(0.1)$ %– $O(1)$ % of the Parwani yield.
- Mass Bounds: Under the Parwani scheme, heavy scalar masses ( $M_H, M_A, M_{H^\pm}$ ) are bounded by $\sim 1.6$ TeV. The AE scheme imposes much stricter bounds, restricting masses to $\lesssim 800$ GeV.
- Structural Differences: The AE prescription produces fragmented, irregular distributions in the heavy-scalar mass planes due to extreme parametric sensitivity to the $\xi_p > 1$ threshold. In contrast, the Parwani scheme yields continuous, stable distributions.
- Macroscopic Parameters: The Parwani scheme predicts larger transition strengths ( $\alpha_{GW} \sim O(1)$ ) and faster transitions ( $\beta_{GW}/H_* \gtrsim 300$ ), whereas the AE scheme restricts $\alpha_{GW} < 0.3$ and $\beta_{GW}/H_* < 300$ .
Yukawa Type Dependence (Parwani Scheme):
- Transition Dynamics: Single-step transitions overwhelmingly dominate the Normal Scenario across all types, contrasting with the Inverted Scenario where multi-step transitions are frequent.
- Mass Spectra: A universal upper bound of $\sim 1.6$ TeV applies to heavy scalars. Lower bounds are driven by experimental constraints (e.g., $M_{H^\pm} \gtrsim 580$ GeV for Type-II/Y) or thermodynamics (e.g., $M_A \gtrsim 200$ GeV for Type-I/X).
- Alignment: All types favor the Higgs alignment limit ( $|s_{\beta-\alpha}| \to 1$ ). Type-I allows the largest deviations ( $\sim 10^{-2}$ ), while Type-II and Type-Y are restricted to $\sim 10^{-3}$ .
- Multi-step Transitions: Two-step transitions are exceedingly rare in the Normal Scenario, appearing only in highly localized regions of Type-I.
Gravitational Wave Signatures:
- Detectability: The GW parameter space is a strict subset of the SFOEWPT space. Type-I yields the largest number of detectable GW points ( $N_{GW} = 1789$ ), followed by Type-X. Type-II and Type-Y are heavily suppressed by FCNC constraints and the SNR requirement.
- Mass Constraints: The requirement of a detectable GW signal ( $SNR \ge 10$ ) imposes tighter mass limits than the SFOEWPT condition alone, restricting heavy scalars to $\lesssim 700$ GeV.
- Hydrodynamic Suppression: All viable GW points operate in a short-lifetime acoustic regime ( $H_*\tau_{sw} \ll 1$ ), leading to a severe suppression factor $\Upsilon \sim O(10^{-2})$ . The signal amplitude is maximized not by the strongest transitions, but by an intermediate range of the order parameter ( $2.8 \lesssim \xi_p \lesssim 3.85$ ) where the trade-off between transition duration and acoustic kinetic energy is optimal.

Significance
The paper establishes that the cosmological viability of the Normal 2HDM and the predictability of its GW signatures are highly sensitive to the theoretical treatment of thermal resummation. The authors argue that the Parwani scheme, which aligns better with partial-dressing formalisms, provides a more stable and physically consistent description of the phase transition landscape compared to the AE scheme.

The results highlight a strong complementarity between future space-based GW observations (LISA) and high-energy collider searches. While LISA can probe the hydrodynamic regime of the early Universe, collider experiments are essential for testing the specific mass ranges and Yukawa structures identified as viable. The study concludes that the detectable GW population in the 2HDM is confined to a specific hydrodynamic regime characterized by short-lived acoustic sources, and that the viable parameter space is significantly more restricted than previously assumed when accounting for thermal resummation uncertainties.

More like this