Higher-dimensional operators at finite-temperature affect gravitational-wave predictions

We investigate the effect of higher-dimensional marginal operators on the thermodynamics of cosmological phase transitions. Focusing on the Abelian Higgs model, we systematically match these operators, which arise at higher orders in the underlying high-temperature expansion of thermal effective field theory, and use field redefinitions to construct a complete, minimal, and gauge-invariant operator basis. We argue that for strong transitions, temporal gauge modes, which enhance the transition strength, should be treated on equal footing with spatial gauge modes in perturbation theory. Marginal operators are found to weaken the transition strength and induce significant uncertainties for strong transitions. For even stronger transitions that could potentially produce a gravitational wave background detectable by LISA, the validity of the high temperature expansion is uncertain, which may impact the applicability of effective theory techniques, including their use in non-perturbative lattice studies.