Matching higher-dimensional operators at finite temperature for general models

Abstract: High-temperature dimensional reduction provides a systematic effective field theory framework for studying finite-temperature thermodynamics and cosmological phase transitions. While the matching of super-renormalizable operators in the resulting three-dimensional effective theories is well established, the matching of higher-dimensional operators has recently been reinvigorated. These operators become phenomenologically relevant in strong first-order phase transitions where they quantify the convergence of the high-temperature expansion. This work automates the matching of generic three-dimensional dimension-five and -six operators for arbitrary models containing scalars, fermions, and gauge fields, implemented as an extension of the Mathematica package DRalgo. We present the operator basis, the matching procedure, and explicit examples including a scalar-Yukawa model, hot QCD, and the full Standard Model up to dimension six, covering operators mixing the strong and electroweak sectors as well as parity-violating contributions. Redundant operators, gauge dependence, and the corresponding field redefinitions are discussed in detail. The code and example model files are publicly available at https://github.com/DR-algo/DRalgo.

• The paper proposes a systematic framework for automated matching of dimension-5 and 6 operators in finite-temperature effective field theories.
• It implements new Mathematica functions to compute Wilson coefficients across various models, ensuring accurate gauge dependence tracking and field redefinitions.
• Validation against scalar-Yukawa, hot QCD, and Standard Model cases underscores its impact on analyzing strong phase transitions and BSM physics.

Automated Matching of Higher-Dimensional Operators in Thermal Effective Field Theory

Context and Motivation

High-temperature dimensional reduction is a standard method in thermal quantum field theory for simplifying the analysis of finite-temperature dynamics and cosmological phase transitions. The process entails integrating out the Matsubara frequencies associated with the hard scale ($\sim \pi T$), yielding a three-dimensional Euclidean effective field theory (EFT). While established techniques exist for matching super-renormalizable operators between the parent theory and its $3d$ EFT, the proper treatment and automation of subleading, higher-dimensional operators (e.g., dimension five and six) has emerged as a necessity for quantitative studies involving strong first-order phase transitions and BSM physics, especially where these operators impact convergence and reliability of the high-temperature expansion.

This paper provides a systematic framework for automated matching of generic dimension-five and dimension-six operators arising from the hard scale in arbitrary models containing scalars, fermions, and gauge fields. The approach is realized as a significant extension of the DRalgo Mathematica package (2605.15176) and supports generic model files, explicit operator bases, and gauge/field redefinition handling.

General Formalism and Operator Basis

The theoretical foundation involves expressing the four-dimensional Lagrangian in a compact tensor notation, covering real scalars, gauge fields, Weyl fermions, and ghosts. The interaction structure makes extensive use of coupling tensors to generalize group-theoretic input.

Dimensional reduction is performed by integrating out hard modes ($\sim \pi T$), restricting attention to the corresponding contributions since soft-mode integrations (including $A_0$ temporal scalars of gauge fields) can be unreliable for strong transitions. The resulting $3d$ operator bases are enumerated for dimension-five and dimension-six, accounting for contributions mixing gauge and electroweak sectors and parity-violating terms. The paper rigorously discusses operator redundancies, gauge dependence, and the necessity of proper field redefinitions, but presents the results in a redundant basis by default, with planned extension to a strictly non-redundant basis.

Algorithmic Automation and Implementation

All operator matching is automated via new Mathematica functions Dimension5Matching and Dimension6Matching, which compute the Wilson coefficients of higher-dimensional operators for arbitrary input group tensors, operator indices, and dimension specifications. Auxiliary routines for inspecting operator shapes (ODIM5, ODIM6), evaluating thermal loop integrals, and constructing group tensors (including symmetrization and contraction routines) are provided. The matching utilizes tensor-notation couplings and thermal one-loop master integrals, and explicit gauge-dependence is tracked.

Validation is carried out across a range of models:

• Scalar-Yukawa Model: Matches known results for singlet fermionic dark matter, reproducing strong numerical results for Wilson coefficients.
• Hot QCD: Fully automates bosonic and fermionic dimension-six operator matching, including explicit gauge-fixing parameter ($\xi$) dependence and field redefinition effects. Fermionic contributions to EQCD are detailed and are relevant for magnetostatic coupling computations at higher orders.
• Standard Model: Extended matching to operators mixing strong and electroweak sectors, and systematic evaluation of parity-violating terms, demonstrating that anomaly-free hypercharge combinations result in vanishing Wilson coefficients, in accordance with previous literature.

Theoretical and Computational Significance

The operator basis construction and matching capabilities enable systematic quantitative assessment of the convergence of the high-temperature operator expansion for strong phase transitions, including analyses where the expansion may break down. The automation of matching for arbitrary models removes practical bottlenecks encountered in analytical treatments and facilitates rapid model exploration, as evidenced by the provided timing benchmarks for both simple and complex models.

Importantly, the inclusion of parity-violating and CP-odd operators allows for model-independent, systematic investigation of baryogenesis scenarios in the presence of finite-temperature effects, which are otherwise highly challenging.

The treatment of gauge dependence in redundant bases and its elimination via field redefinitions is central to ensuring physical results, and the paper lays the groundwork for future implementation of automated non-redundant physical bases.

Implications and Future Directions

The practical implications are considerable: the ability to reliably automate higher-dimensional operator matching at finite temperature is instrumental for Monte-Carlo studies, non-perturbative EFT analyses, and model-building for gravitational wave cosmology. The package's modular structure allows easy extension to new operator classes and two-loop matching, as well as direct integration with lattice simulation pipelines and advanced renormalization-group calculations.

For theoretical developments, the framework enables:

• Exploration of breakdown regimes for the high-temperature expansion by systematically quantifying higher-dimensional contributions to thermodynamic observables.
• Comparison between $3d$ and $4d$ lattice simulations for strong phase transitions in BSM models.
• Extension to soft-to-softer scale matching with full operator bases, and inclusion of automated field redefinitions for non-redundant basis extraction.
• Investigation of CP-odd effects in QCD and the Standard Model at high temperature, relevant for baryogenesis.

Conclusion

This work provides a robust formal and computational framework for the automated matching of higher-dimensional operators (dimension five and six) at finite temperature for general models, realized as an extension to DRalgo. The toolset and formalism facilitate systematic studies of strong cosmological phase transitions, sharpen model-building capabilities for thermal EFTs, and open new avenues for rigorous analysis of BSM baryogenesis scenarios and IR behavior in hot gauge theories. The approach is validated against existing literature and ready for high-impact applications in thermal field theory, cosmological phenomenology, and computational physics (2605.15176).