Matching higher-dimensional operators at finite temperature for general models

This paper presents an automated framework, implemented as an extension of the DRalgo Mathematica package, for matching generic dimension-five and -six operators in three-dimensional effective field theories derived from arbitrary models containing scalars, fermions, and gauge fields at finite temperature.

The Gist

The Big Picture: Simplifying a Hot Mess

Imagine the universe just after the Big Bang as a giant, boiling pot of soup. This soup is filled with all kinds of particles (like ingredients) moving around at incredibly high speeds. Physicists want to understand how this soup behaves, especially during "phase transitions"—moments where the soup suddenly changes its state, like water turning into ice or steam.

To study this, scientists use a technique called Dimensional Reduction. Think of this like taking a complex, 3D movie and compressing it into a 2D cartoon. The 3D movie (the real, high-temperature universe) is too complicated to calculate directly. So, physicists create a simpler "effective" version (the 2D cartoon) that captures the most important behaviors but ignores the tiny, fast-moving details that don't matter for the big picture.

The Problem: Missing the "Spices"

For a long time, scientists had a good recipe for this 2D cartoon. They knew how to handle the main ingredients (the basic particles and forces). However, they were missing the "spices"—the subtle, high-level interactions that only show up when the soup is boiling very hard.

In physics terms, these are called higher-dimensional operators.

• The Old Way: They could only calculate the main flavors (super-renormalizable operators).
• The New Problem: When the phase transition is very strong (like a violent explosion rather than a gentle freeze), those missing "spices" become crucial. If you ignore them, your prediction of the explosion is wrong.
• The Challenge: Calculating these spices by hand is like trying to solve a Sudoku puzzle while juggling chainsaws. It is incredibly tedious, prone to human error, and takes forever.

The Solution: A New "Auto-Chef"

The authors of this paper have built a new tool inside a software package called DRalgo (Dimensional Reduction algorithm). Think of this software as an automated chef.

Previously, the chef could only chop the main vegetables. Now, with this new update (version 1.5.0), the chef can:

1. Identify the missing spices: It automatically figures out exactly which complex interactions (dimension-5 and dimension-6 operators) need to be added to the 2D cartoon.
2. Calculate the amounts: It does the heavy math to determine exactly how much of each "spice" is needed based on the original 3D recipe.
3. Do it for any model: Whether you are cooking a simple soup (a scalar-Yukawa model), a spicy curry (Hot QCD), or a massive banquet (the full Standard Model), this tool can handle it.

How It Works (The Analogy)

Imagine you have a complex blueprint for a skyscraper (the 4D theory). You want to build a model of it on a table (the 3D effective theory).

• The Old Method: You had to manually measure every window, door, and beam, then write down the instructions for the model. If you missed a tiny detail, the model would collapse.
• The New Method: You feed the blueprint into a 3D printer (the DRalgo software). The printer automatically scans the blueprint, realizes, "Oh, this skyscraper has these specific, weirdly shaped balconies that only show up when the building is hot," and automatically prints the instructions for those balconies into the model.

What They Actually Did

The paper doesn't just talk about the tool; they tested it on three specific "recipes":

1. Scalar-Yukawa Model: A simple theoretical soup. They checked their tool against known results and confirmed it worked perfectly.
2. Hot QCD (Quantum Chromodynamics): This is the physics of the strong nuclear force (what holds atoms together). They calculated the "spices" for this hot environment, including how the "temporal" parts of the force fields behave.
3. The Standard Model: This is the complete recipe of our known universe (electrons, quarks, Higgs boson, etc.). They successfully calculated the complex interactions that mix the strong force with the weak and electromagnetic forces, and even found interactions that violate "parity" (a type of symmetry, like how your left hand is a mirror image of your right).

Key Takeaways for the Reader

• Automation is Key: The math required to find these high-level interactions is too hard for humans to do reliably. This software automates the process.
• Accuracy for Strong Events: If a phase transition in the early universe was violent, these new calculations are necessary to get the physics right.
• Gauge Dependence: The authors noted that some of these calculations look different depending on how you "view" the math (gauge dependence), but when you put the final pieces together, the result is consistent and correct.
• Availability: They didn't keep the tool secret. They made the code and the example "recipes" available on GitHub for anyone to use.

What They Did Not Do

The paper is strictly about building the tool and testing it.

• They did not use this tool to predict a specific new particle that we will find tomorrow.
• They did not claim this solves the mystery of why the universe exists (though it helps us understand the conditions that could lead to it).
• They did not apply this to medical or clinical scenarios.

In short, they built a better calculator for theoretical physicists so that when they study the hot, early universe, they don't have to worry about missing the subtle details that could change the whole story.

Technical Summary

Technical Summary: Matching Higher-Dimensional Operators at Finite Temperature for General Models

Problem Statement
High-temperature dimensional reduction serves as a systematic effective field theory (EFT) framework for analyzing finite-temperature thermodynamics and cosmological phase transitions. While the matching of super-renormalizable operators in the resulting three-dimensional (3d) effective theories is well-established, the matching of higher-dimensional operators (specifically dimension-five and dimension-six) has historically been a manual and complex process. These operators, arising from subleading terms in the high-temperature expansion, are phenomenologically critical for describing strong first-order phase transitions, where they quantify the convergence of the expansion and impact non-perturbative phenomena. The complexity of including these operators necessitates automation, particularly for models with large field contents like the Standard Model (SM) and Beyond the Standard Model (BSM) scenarios.

Methodology
The authors present an extension to the Mathematica package DRalgo (version 1.5.0) that automates the matching of generic 3d dimension-five and dimension-six operators for arbitrary models containing scalars, fermions, and gauge fields. The methodology involves:

1. Operator Basis Definition: The authors define a redundant operator basis for the 3d effective theory obtained by integrating out the hard scale $\pi T$. This basis includes:
   • Dimension-5: Operators involving scalar fields ($R_s$), temporal gauge fields ($A_0^a$), spatial field strengths ($F_{ij}^a$), and covariant derivatives ($D_i$). This includes parity-violating terms and operators mixing scalar and gauge sectors.
   • Dimension-6: A comprehensive set of 20 independent operators for the Standard Model, including sextic scalar terms, quartic gauge terms, and mixed derivative interactions.
2. Tensor Notation and Matching: The matching procedure utilizes a compact tensor notation for couplings (e.g., $g_{abc}$, $Y_{sIJ}$, $\lambda_{s_1...s_4}$) implemented in DRalgo. The Wilson coefficients of the 3d operators are determined by matching $n$-point correlation functions.
3. Thermal Integrals: The calculation relies on hard thermal one-loop master integrals ($I_{b,f}$) evaluated in the background-field $R_\xi$ gauge. The results are expressed in terms of these integrals and the renormalization scale $\Lambda$.
4. Implementation: New functions Dimension5Matching and Dimension6Matching allow users to input group tensor structures and automatically solve for the corresponding Wilson coefficients. Helper functions (ODIM5, ODIM6, HardThermal1LoopInt, Contract, SymmetrizeTensor) facilitate the construction of operator tensors and the evaluation of thermal integrals.
5. Handling Redundancy and Gauge Dependence: The paper acknowledges that the chosen operator basis is redundant (some operators are proportional to equations of motion). Consequently, Wilson coefficients may exhibit explicit dependence on the gauge-fixing parameter $\xi$. The authors retain this dependence to verify that it cancels upon appropriate field redefinitions to reach a physical (on-shell) basis.

Key Contributions and Results
The paper provides the following specific technical outputs:

• Automated Matching Code: The extension of DRalgo allows for the systematic derivation of dimension-five and dimension-six operators for any model implementable within the package.
• Explicit Examples:
  • Scalar-Yukawa Model: The authors reproduce known results for dimension-five and dimension-six operators, validating the code against previous literature.
  • Hot QCD: The dimension-six Lagrangian for hot QCD (EQCD) is derived, including fermionic contributions. The authors explicitly relate the redundant basis coefficients ($\alpha_i$) to the physical basis coefficients ($\beta_i$), demonstrating the cancellation of gauge-parameter dependence in the physical basis. They also compute the fermionic contribution to the magnetostatic wavefunction renormalization at two loops, identifying a UV divergence relevant for $O(g^8)$ calculations of the electrostatic gauge coupling.
  • Standard Model (SM): The full dimensional reduction of the SM up to dimension six is presented. This includes operators mixing the strong and electroweak sectors and parity-violating operators. The authors note that while these parity-violating operators are generated in the general computation, their Wilson coefficients vanish for the physical hypercharge assignments of the SM due to anomaly cancellation conditions.
• Technical Guidelines: The paper details technical aspects regarding the handling of spacetime dimensions (specifically the limitations of fermionic matching in $d \neq 3$ due to Chisholm identities), the necessity of correct symmetry properties for input group tensors, and power-counting schemes.
• Performance: Benchmarks indicate that the matching of dimension-six operators for the full Standard Model takes approximately 2045 seconds on an Apple M3 MacBook Air, with the majority of time spent on solving for Wilson coefficients rather than tensor contractions.

Significance and Claims
The authors claim that this work completes the set of matching relations required to determine all operators up to and including marginal ones in 3d EFTs. The primary significance lies in the ability to:

1. Systematically Test Validity: With the full operator basis available, researchers can now systematically test the convergence of the high-temperature expansion in large-scale phase-transition studies, particularly for strong transitions where the high-temperature assumption might break down.
2. Investigate CP Violation: The automated matching of parity-violating operators enables a systematic investigation of CP-odd effects at high temperatures in generic models. This is highlighted as a new handle for probing BSM baryogenesis scenarios in a model-independent way, given that CP violation is a necessary ingredient for baryogenesis.
3. Facilitate Future Work: The code and example files are publicly available, lowering the barrier for future studies involving higher-dimensional operators, two-loop matching, and automated field redefinitions.

The paper remains modest regarding immediate experimental proposals, focusing instead on providing the theoretical infrastructure and computational tools necessary for precise EFT descriptions of finite-temperature phase transitions.