Higher-dimensional operators and Polyakov loop in hot Scalar QED from the heat kernel

Abstract: Using the finite-temperature heat kernel method, we compute the gauge-invariant effective Lagrangian up to dimension-six for massive hot scalar QED. We propose two complementary methods: integrating out heavy modes at finite temperature, and deriving the finite-temperature heat kernel coefficients from the zero-temperature ones. We show that in the static limit, both lead to the same three-dimensional effective operators. We also compute the gauge-invariant Coleman-Weinberg effective potential for a constant background at finite temperature. We further examine how the Polyakov loop modifies the matching coefficients and assess its impact together with the higher-dimensional operators on the thermodynamics of cosmological first-order phase transitions, which in turn can affect an associated gravitational-wave spectrum.

• The paper constructs the finite-temperature one-loop effective action for massive Scalar QED using the heat kernel formalism to systematically include higher-dimensional operators up to dimension six.
• It demonstrates the equivalence of direct finite-temperature integration and zero-temperature matching methods, ensuring gauge invariance and precise operator extraction.
• It reveals that resumming Polyakov loop effects is crucial for accurately predicting phase transition observables, impacting gravitational wave signals from cosmological first-order phase transitions.

Higher-Dimensional Operators and Polyakov Loop in Finite-Temperature Scalar QED via the Heat Kernel

Introduction and Motivation

This work presents a systematic construction of the finite-temperature one-loop effective action for massive scalar quantum electrodynamics (QED) utilizing the heat kernel formalism, with explicit focus on operators up to dimension six. The motivation is rooted in the need for high-precision calculations of thermodynamic parameters relevant to cosmological first-order phase transitions (FOPTs), which are anticipated targets of future gravitational-wave observatories such as LISA. Specifically, precision in quantities like transition strength $\alpha$, transition rate $\beta/H$, and bubble wall velocities is essential for robust predictions—necessitating the inclusion of higher-dimensional operators suppressed by the inverse temperature scale.

The standard approach employing dimensional reduction typically retains only renormalizable three-dimensional operators, sufficient for leading-order thermodynamics. However, subleading contributions arise from an infinite tower of higher-dimensional terms, which become non-negligible for next-generation gravitational wave predictions. While diagrammatic techniques for extracting such operator bases exist, the heat kernel formalism provides an algorithmic, gauge-covariant procedure to generate the effective action, inherently capturing the infinite operator tower.

Formalism and Heat Kernel Approaches

Two complementary methodologies are expounded for deriving the finite-temperature effective action:

1. Direct finite-temperature integration: The non-zero Matsubara modes are integrated out directly in the compactified Euclidean spacetime, with the heat kernel constructed to respect the temporal $S^1_\beta$ topology and operator mixing.
2. Matching from zero-temperature coefficients: Finite-temperature heat kernel coefficients are constructed from their zero-temperature analogs using systematic matching relations, effectively mapping UV contributions from the full theory to the three-dimensional EFT.

The equivalence of these methods is explicitly demonstrated in the static limit, ensuring consistency of the resultant operator bases and Wilson coefficients.

The operator basis is reduced to a minimal, gauge-invariant set employing field redefinitions and equations of motion (EOMs), paralleling the reduction protocols used for the Abelian Higgs model in dimensional reduction.

Dimension-Six Operator Basis and Wilson Coefficient Extraction

The complete one-loop effective action is presented up to and including dimension-six operators, for both scalar and gauge sectors, with explicit expressions for the Wilson coefficients. These are derived for both mass-degenerate and non-degenerate cases, highlighting the smoothness of the degenerate mass limit and the absence of singularities. The contributions of both local (potential, kinetic, and higher-derivative) and non-local Polyakov loop-dressed terms are detailed, with the master integrals evaluated via thermal Matsubara sums.

Crucially, the heat kernel formalism automatically resums Polyakov loop effects, incorporating the temporal background holonomy as a non-perturbative input parameter, and yielding genuinely gauge-invariant matching coefficients that depend on the gauge charge.

The Coleman-Weinberg potential is also computed within this framework, including finite-temperature corrections in terms of background-dependent masses. Kinetic term corrections impacting the tunneling/bounce action relevant for phase transition dynamics are retained.

Polyakov Loop Effects and Phase Transition Thermodynamics

A significant result of the work is the identification and quantification of the effect of nontrivial Polyakov loop backgrounds (holonomy) on thermodynamic observables. By treating the Polyakov loop as an external parameter, the analysis demonstrates its impact on the phase structure:

• Critical Temperature and Transition Strength: Increasing the Polyakov holonomy monotonically increases the critical temperature $T_c$ and suppresses the transition strength $\alpha$. At a maximal holonomy, the transition becomes second order and $\alpha$ vanishes.
• Implications for Gravitational Wave Phenomenology: The explicit dependence of $T_c$ and $\alpha$ on the Polyakov loop translates directly into the predicted gravitational-wave spectra from cosmological FOPTs. Neglecting the Polyakov loop (i.e., assuming trivial holonomy) leads to underestimation of $T_c$ and overestimation of $\alpha$, rendering the accurate resummation of holonomy indispensable for precision cosmological predictions.

Comparison with Diagrammatic and Automated Dimensional Reduction

A thorough comparison is performed between the operator bases and matching coefficients obtained via the heat kernel and established diagrammatic (e.g., DRalgo) dimensional reduction methods. Agreement is found in the static and expanded local limit, with differences arising from:

• Treatment of Holonomy: Diagrammatic approaches yield explicit towers of $\beta/H$0 operators, while the heat kernel resums all Polyakov loop effects into master sums, capturing nonlocal holonomy physics automatically.
• Redundant Operator Elimination: Heat kernel outputs include EOM-redundant terms, requiring systematic reduction to the minimal operator basis.
• Generality and Automation: The heat kernel approach is algorithmic and generalizable to broad classes of gauge theories (non-Abelian, SMEFT), providing a robust path to further automation.

Implications and Outlook

The analysis substantiates the critical role of higher-dimensional operators generated at finite temperature for high-precision studies of thermal phase transitions and associated gravitational wave signals. The robust capture of Polyakov loop physics by the heat kernel is shown to be essential for correct thermodynamics in gauge theories at finite temperature, signaling the necessity for future automated tools to incorporate full holonomy dependence in their matching procedures.

The techniques developed are immediately extensible to non-Abelian settings and Standard Model Effective Field Theory (SMEFT) at finite temperature. The explicit demonstration that matching from zero-temperature coefficients suffices to produce the correct thermal operators streamlines the construction of effective actions for a wide array of physical contexts.

Future research directions include:

• Integration of heat kernel-based functional matching into automated EFT toolchains.
• Extension to cases with dynamical fermionic content in the loops.
• Detailed study of the impact of holonomy-dressed higher-dimensional operator corrections on bubble nucleation dynamics and gravitational-wave spectra across parameter space, particularly for models probed by LISA.

Conclusion

This paper establishes a precise, gauge-invariant framework for constructing higher-dimensional operators in finite-temperature EFTs using the heat kernel formalism, with explicit treatment of Polyakov loop holonomy. The results clarify the interplay of operator matching, holonomy-resummed thermal effects, and renormalization in Abelian gauge theories and set a new standard for the theoretical control of effective actions governing cosmological phase transitions and gravitational wave predictions. The outlined methodology paves the way for further algorithmic advances in the computation of finite-temperature effective actions in gauge theories relevant for new physics searches in cosmology and beyond.