⚛️ phenomenology

Thermal masses in the Standard Model at three loops

This paper computes the Higgs thermal mass and electroweak and colour Debye mass parameters to full $\mathcal{O}(g^6)$ accuracy in dimensional reduction by employing new methods to evaluate previously unknown mixed-signature sum-integrals.

Original authors: Mikael Chala

Published 2026-09-29

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Mikael Chala

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the earliest moments of the universe, long before stars ignited or atoms formed, the cosmos was a seething, superheated soup of fundamental particles. In this extreme environment, particles do not behave as they do in the cold vacuum of space today. Instead, they interact with a dense thermal bath, a collective state of matter that fundamentally alters their properties. One of the most significant changes is that particles, which might be massless or light in a vacuum, acquire an effective mass simply by moving through this hot plasma. This phenomenon is not unique to the early universe; it occurs whenever charged particles move through a hot gas, causing the gas to rearrange itself and partially shield the particle's influence, much like how a crowd might shift to cushion a person walking through it. Physicists call this shielding effect Debye screening, and the resulting "thermal mass" is a crucial number for understanding how the universe evolved, how energy moves through hot matter, and how the fundamental forces of nature behaved at high temperatures.

For decades, scientists have worked to calculate these thermal masses with increasing precision, but the mathematics involved is notoriously difficult. The challenge lies in the fact that the interactions between particles happen at different energy scales simultaneously. Some interactions are very strong and fast, while others are weaker and slower. To get an accurate picture, researchers must account for all these layers of interaction at once. Until now, the most complete calculations for the Standard Model—the theory describing all known fundamental particles and forces—had only reached a certain level of accuracy, leaving out subtle but important effects that only appear when the math is pushed to its absolute limit.

A new study by Mikael Chala at the University of Granada has finally pushed these calculations to their full three-loop limit, a technical term meaning the researchers have accounted for the most complex web of particle interactions possible within the current framework of the theory. In this work, Chala computed the thermal masses for the Higgs boson and the gauge bosons that carry the electroweak and strong nuclear forces. The goal was to determine these masses with a precision that includes terms up to the sixth power of the coupling constants, which represent the strength of the forces. This level of detail is essential because thermal masses are not just abstract numbers; they are the keys to predicting how the universe behaved during phase transitions, such as the moment when the electromagnetic and weak forces separated, and they influence calculations for the pressure of the early universe and how particles annihilate.

To achieve this, the researcher had to solve a specific mathematical hurdle that had previously blocked progress: the evaluation of "mixed-signature sum-integrals." In the language of the theory, these are complex calculations that involve summing up the contributions of both bosons and fermions, two distinct families of particles, as they move through the thermal bath. While simpler versions of these calculations had been solved before, the specific combinations required for the full Standard Model at this high level of precision had remained unknown. Chala developed new computational methods to tackle these previously unsolved integrals, combining several advanced mathematical techniques to break the problem down into manageable pieces. By doing so, he was able to evaluate these integrals numerically with high precision, effectively clearing the path to the final result.

The findings reveal the precise values for the thermal masses of the Higgs field and the color-charged gluons that mediate the strong force. The study provides numerical expressions for these masses that include contributions from the top quark's interaction and the various gauge forces. One of the most significant outcomes is the dramatic improvement in stability when predicting how these masses change with energy scales. When the results are plotted against temperature, the range of uncertainty narrows significantly compared to previous, less precise calculations. For instance, at a temperature of one trillion electronvolts, the uncertainty band for the color Debye mass shrinks by approximately seventy percent compared to earlier estimates. This reduction in uncertainty is vital for making reliable predictions about the thermal history of the universe.

The paper also confirms that the new, more complex calculations are consistent with the fundamental rules of the theory. The researchers verified that the final results are free of mathematical infinities that often plague such high-order calculations, a sign that the method is sound and the physics is correctly captured. While the new, highly precise values for the thermal masses do not drastically change the qualitative picture of the early universe, they provide a much firmer foundation for future work. This includes the study of phase transitions in the early universe and the development of theories that go beyond the Standard Model. By establishing these precise benchmarks, the work ensures that any future discoveries in high-energy physics or cosmology will be built upon a solid and accurate understanding of how particles behave in the hottest environments imaginable.

Technical Summary: Thermal Masses in the Standard Model at Three Loops

Problem Statement
Elementary particles acquire thermal effective masses ( $m_{\text{eff}}$ ) when interacting with a thermal bath, a phenomenon known as Debye screening. In Quantum Field Theory (QFT), these masses encode quantum corrections to two-point particle propagators. Accurate determination of these masses is critical for predicting observables such as pressure, spatial string tension, static heavy-quark screening, jet broadening, and particle annihilation rates. Furthermore, recent studies indicate that $O(g^6)$ hard matching is essential for precision studies of thermal phase transitions within dimensionally reduced Effective Field Theories (EFTs).

While $O(g^4)$ Higgs and Standard Model (SM) Debye masses were previously computed, and $O(g^6)$ results existed for pure Yang-Mills and partial QCD cases, a full $O(g^6)$ computation for the Higgs and electroweak/colour Debye masses in the complete Standard Model remained incomplete. Specifically, certain mixed-signature (bosonic-fermionic) three-loop sum-integrals required for the full SM calculation had not been evaluated.

Methodology
The author employs the dimensional reduction program, which exploits the separation of energy scales (hard $E \sim \pi T$ , soft $E \sim gT$ , and ultrasoft $E \sim g^2 T$ ) to integrate out the hard scale. This results in a 3-dimensional EFT containing only bosonic fields, as fermionic fields possess no zero modes due to anti-periodic boundary conditions.

The calculation proceeds by matching the 3D EFT coefficients to the hard-region expansion of static two-point correlators in the full 4D SM. The process involves:

Lagrangian Setup: Using the Minkowski SM Lagrangian in the background-field gauge, splitting gauge bosons into background and quantum fluctuations.
Perturbative Expansion: Computing off-shell masses up to three loops, kinetic terms up to two loops, and four-derivative operators at one loop. The power counting assumes $\mu^2, \lambda, y_t^2, g_1^2, g_2^2, g_s^2 \sim g^2$ .
Sum-Integral Evaluation: The core technical challenge involves evaluating three-loop sum-integrals. While two-loop integrals factorize into products of one-loop integrals, three-loop integrals generally do not. The author utilizes integration-by-parts (IBP) relations and the software SIRENA to reduce these to master integrals.
New Numerical Techniques: For previously unevaluated mixed-statistics (bosonic-fermionic) three-loop sum-integrals, the author develops a numerical evaluation method combining:
- Schwinger parametrization of propagators.
- Sector decomposition to split integration ranges.
- Partial Poisson summation to swap sums over Matsubara frequencies for winding modes (applied only to "small lines").
- Analytic endpoint subtraction to isolate and remove $1/\epsilon$ poles, leaving finite integrals for the Laurent coefficients.

Key Contributions

Full $O(g^6)$ Computation: The paper provides the first complete $O(g^6)$ calculation for the Higgs thermal mass and the electroweak and colour Debye mass parameters in the full Standard Model.
Evaluation of Mixed Sum-Integrals: The author presents the first numerical determination of Laurent coefficients for mixed bosonic-fermionic three-loop thermal sum-integrals. These results are tabulated (Table I) and made available in an ancillary file (matching.txt).
Renormalization and Consistency: The calculation demonstrates that the renormalized masses are completely pole-free. The author verifies this by showing that explicit poles from factorizable integrals cancel against those from non-factorizable three-loop integrals.
Scale Dependence Analysis: The paper analyzes the renormalization scale dependence of the masses, showing a significant reduction in the uncertainty band at higher orders.

Results
The author provides numerical expressions for the renormalized Higgs mass ( $m_\phi^2$ ) and the colour Debye mass ( $m_{G0}^2$ ) through $O(g^6)$ in the limit where $g_1, \lambda, \mu^2 \to 0$ and $g_2^2 \sim g_3$ .

Higgs Mass: The result includes terms proportional to $g_2^2$ , $y_t^2$ , and various mixed terms up to $y_t^6$ and $g_s^4 y_t^2$ , with explicit dependence on the logarithm of the renormalization scale $\ell = \log(\mu/T)$ .
Colour Debye Mass: The result includes terms up to $g_s^6$ and mixed terms involving $y_t$ . Notably, the finite pieces of the new mixed sum-integrals do not affect the colour Debye mass, as they enter the matching proportional to $(d-3)$ , a finding consistent with pure QCD results.
Scale Stability: As shown in Figures 2 and 3, the renormalization-scale band for the colour Debye mass narrows by approximately 70% at $T = 1$ TeV when moving from $O(g^4)$ to $O(g^6)$ accuracy.

Significance
The paper claims that achieving full $O(g^6)$ accuracy substantially extends previous results and is instrumental for precision studies of thermal phase transitions within dimensionally reduced EFTs. The methodology for evaluating mixed-signature sum-integrals provides a systematic route for higher-order calculations and is applicable to new-physics scenarios, such as the determination of thermal masses in the Standard Model Effective Field Theory (SMEFT), which the author notes will be addressed in forthcoming work. The work serves as a foundational step for improving the theoretical precision of thermal observables in the early universe and heavy-ion collision contexts.

More like this