Hard thermal contributions to phase transition observables at NNLO

Technical Summary: Hard Thermal Contributions to Phase Transition Observables at NNLO

Problem Statement
The potential observation of a stochastic gravitational wave (GW) background from cosmological first-order phase transitions (PTs) requires precise theoretical predictions of thermodynamic parameters. While effective field theory (EFT) methods, specifically high-temperature dimensional reduction (DR), are standard tools for studying these transitions, existing literature often truncates the expansion at leading order (LO) or next-to-leading order (NLO) in the 3d EFT. This truncation omits higher-dimensional operators (beyond super-renormalizable terms) and higher-loop corrections, leading to significant uncertainties in nucleation computations, particularly regarding gauge and scale dependence. For moderately strong to strong transitions, the inclusion of marginal operators (e.g., dimension-6) and higher-order loop corrections is necessary to reduce these uncertainties and ensure the validity of the perturbative expansion.

Methodology
The authors perform a systematic calculation of hard thermal contributions to the Abelian Higgs model (U(1) gauge theory coupled to a complex fundamental scalar) up to order $O(g^6)$ in the gauge coupling. The methodology involves:

1. Dimensional Reduction (DR): Integrating out hard modes (Matsubara frequencies of order $\pi T$) to construct a static 3d EFT for the soft bosonic modes. The authors utilize the "DR-B" approach, where the effective potential is computed at the soft scale, avoiding the breakdown of the operator expansion seen in the "DR-A" approach for strong transitions.
2. High-Order Matching:
   • Masses: Computing the scalar and Debye masses at the three-loop level (NNLO). This involves matching the poles of the static propagators in the 4d theory to the 3d EFT.
   • Couplings: Computing the quartic couplings ($\lambda_3, h_3, \kappa_3$) at the two-loop level.
   • Dimension-6 Operators: Including one-loop contributions to dimension-6 operators in the 3d EFT.
3. Renormalization and Gauge Independence: The calculation employs the $\overline{\text{MS}}$ scheme. The authors explicitly demonstrate the cancellation of gauge dependence order-by-order in loops, a crucial consistency check for physical observables.
4. Master Integrals: The computation relies on a basis of three-loop bosonic sum-integrals. The authors utilize integration-by-parts (IBP) relations to reduce diagrams to a set of master integrals. They identify a potential mismatch in the $\mathcal{O}(\epsilon^{-1})$ and $\mathcal{O}(\epsilon^0)$ parts of a specific master integral ($I_{31111-2}$) required for the scalar mass renormalization.
5. Thermodynamic Analysis: Using the derived effective potential (including LO, NLO, NNLO, and dimension-6 terms), the authors calculate critical parameters such as the critical temperature ($T_c$), latent heat, and transition strength ($\alpha_c$). They compare the magnitude of higher-dimensional operator effects against higher-loop corrections.

Key Contributions and Results

• Three-Loop Masses and Two-Loop Couplings: The paper provides the first complete derivation of the three-loop scalar and Debye masses for U(1) and SU(N) gauge-Higgs models, as well as the two-loop quartic couplings for the Abelian case, achieving $O(g^6)$ accuracy.
• Gauge Independence: The authors prove that the gauge dependence cancels completely in the physical 3d parameters (masses and couplings) once the three-loop matching contributions are included. This resolves a residual gauge dependence found in previous $O(g^6)$ calculations.
• Master Integral Mismatch: A significant technical finding is the identification of a missing contribution in the master integral $I_{31111-2}$. The standard expression for this integral leads to a residual $1/\epsilon$ pole in the renormalized scalar mass. The authors show that setting a specific parameter $X = -2/5$ (modifying the integral's finite parts) restores finiteness and consistency in both U(1) and SU(N) theories. This suggests a previously unaccounted error in the literature regarding this specific integral.
• Dominance of Dimension-6 Operators: In the regime of strong phase transitions (characterized by small $x = \lambda/g^2$), the authors quantify that one-loop dimension-6 effects typically dominate over two- and three-loop corrections to super-renormalizable parameters. Specifically, for $x_{LO} \lesssim 1$, the corrections from dimension-6 operators exceed those from three-loop hard-scale diagrams.
• Thermodynamic Parameters: The study provides expressions for critical parameters ($T_c$, $\Delta \langle \phi^\dagger \phi \rangle$, $\Delta \langle (\phi^\dagger \phi)^2 \rangle$) up to N2LO + dimension-6. These results are implemented in the software DRalgo.

Significance and Claims
The paper claims to provide the necessary theoretical infrastructure for accurate predictions of gravitational wave signals from strong first-order phase transitions. By establishing the $O(g^6)$ EFT, the authors argue that:

1. Precision: The inclusion of dimension-6 operators is not merely a higher-order correction but is often the dominant effect in the phenomenologically interesting regime of strong transitions, necessitating their inclusion to reduce theoretical uncertainties in GW predictions.
2. Consistency: The work demonstrates that the 3d EFT is fully gauge-independent up to $O(g^6)$, validating the perturbative approach for these calculations.
3. Technical Correction: The identification of the mismatch in $I_{31111-2}$ corrects a subtle error in the renormalization of thermal masses, ensuring the consistency of 4d and 3d renormalizability.
4. Future Work: The results pave the way for computing the bubble nucleation rate at $O(g^6)$ accuracy and studying non-equilibrium dynamics. The authors note that while dimension-6 operators dominate near the critical temperature, their impact may be less pronounced at the nucleation temperature ($T_n$) for classically conformal models, a topic reserved for future investigation.

The authors maintain a modest stance, noting that their analysis is specific to the Abelian Higgs model as a prototype, though the methods and findings regarding the master integrals generalize to SU(N) theories. They do not claim to have solved the non-perturbative regime but rather to have bridged the gap between weak and strong transitions within the perturbative framework.