Non-perturbative determination of the QCD Equation of State up to the electroweak scale

This paper presents a non-perturbative determination of the QCD Equation of State for three massless quark flavors across a temperature range from 3 GeV to the electroweak scale, revealing that even at these extreme temperatures, perturbative predictions remain insufficient without incorporating non-perturbative contributions.

The Gist

Imagine the universe as a giant, cosmic soup. In the very beginning, just after the Big Bang, this soup was incredibly hot and dense. It wasn't made of atoms or molecules like the soup in your kitchen; it was a "primordial broth" made of the most fundamental building blocks of matter: quarks and gluons.

Physicists call this state of matter Quark-Gluon Plasma (QGP). To understand how the universe expanded and cooled, scientists need to know the "recipe" for this soup: specifically, how much pressure it exerts and how much energy it holds at different temperatures. This recipe is called the Equation of State (EoS).

For decades, scientists could only calculate this recipe for "cool" soup (temperatures around 150 million degrees) or for "very hot" soup using math approximations that often broke down. But there was a huge gap in the middle: the scorching temperatures between 3 billion degrees and 165 billion degrees (the scale of the Electroweak era).

This paper, by Michele Pepe and his team, is like a master chef finally tasting that missing, super-hot part of the soup and writing down the exact recipe using a method that doesn't rely on guesswork.

Here is how they did it, explained through simple analogies:

1. The Problem: The "Zoom" Dilemma

Imagine trying to take a photo of a tiny ant (a quark) while it's running on a giant treadmill (the hot universe).

• To see the ant clearly, you need a camera with a very high zoom (a tiny grid size).
• To see the whole treadmill, you need a wide-angle lens (a huge space).
• The Catch: If you zoom in too much to see the ant, the treadmill looks tiny and you can't see the motion. If you zoom out to see the treadmill, the ant disappears.
• The Old Way: Previous calculations tried to use "math approximations" (perturbation theory) to guess what the ant was doing at high speeds. But at these extreme temperatures, the math gets messy, like trying to predict the weather by only looking at the wind speed and ignoring the humidity. The predictions were off.

2. The Solution: The "Moving Room" Trick

The team used a clever trick called Shifted Boundary Conditions.

• The Analogy: Imagine you are in a room with a clock on the wall. Usually, you stand still and watch the clock tick. But what if you start walking sideways while looking at the clock?
• In physics, this "walking" (shifting the boundary) changes how time and space look to you. It's like a Lorentz boost (a fancy term for moving fast).
• By "walking" in their computer simulation, the team could measure the "entropy" (a measure of disorder or the number of ways the soup can be arranged) directly.
• Why it helps: Normally, to measure this, you have to subtract a huge amount of "background noise" (vacuum energy) which is like trying to hear a whisper in a rock concert. The "moving room" trick cancels out the noise automatically, letting them hear the whisper clearly without needing to subtract the roar.

3. The Map: Finding the "Constant Path"

To simulate the universe cooling down, they needed to know exactly how to tune their computer settings as they changed the temperature.

• The Analogy: Imagine you are hiking up a mountain (increasing energy). You need a map that tells you exactly where you are. Usually, hikers use landmarks like "the big oak tree" or "the red barn." But at the top of the mountain (high energy), there are no trees or barns; it's just snow and rock.
• The New Map: Instead of using landmarks, they used a "compass" based on how a specific force (the coupling constant) changes as you climb. They followed a "Line of Constant Physics," ensuring that no matter how high they went, they were always measuring the same physical reality, just at different scales. This allowed them to bridge the gap from the "forest floor" (low energy) to the "snowy peak" (electroweak scale).

4. The Results: The Recipe is Complex

Once they ran the simulations, they found something surprising:

• The Expectation: They thought that at such high temperatures, the soup would behave like a simple, ideal gas (like air in a balloon). The math suggested the "recipe" should be simple.
• The Reality: Even at 165 billion degrees, the soup is not simple. The team found that the "ideal gas" math was missing crucial ingredients. There are hidden, complex interactions (non-perturbative contributions) that only show up when you look closely.
• The Takeaway: Even when things are super hot and fast, the universe still has "secret sauce" that simple math can't predict. You need the full, complex recipe to get it right.

Why Does This Matter?

This isn't just about cooking soup.

1. Cosmic History: It helps us understand exactly how the universe expanded in its first fractions of a second.
2. Gravitational Waves: The way this soup behaved affects the "ripples" in space-time (gravitational waves) that might still be echoing today. If we get the recipe wrong, we might miss the signal of the Big Bang.
3. Future Tech: This method is so robust that it can be applied to other types of matter (with 4 or 5 types of quarks), helping us understand the fundamental laws of nature better than ever before.

In a nutshell: The team built a super-precise computer simulation that "walks" through time to measure the properties of the universe's hottest moments. They discovered that even at the hottest temperatures imaginable, the universe is more complex and interesting than our simple math models predicted. They finally wrote down the true recipe for the early universe.

Technical Summary

Here is a detailed technical summary of the paper "Non-perturbative determination of the QCD Equation of State up to the electroweak scale" by Michele Pepe.

1. Problem Statement

The Equation of State (EoS) of Quantum Chromodynamics (QCD) describes the thermodynamic properties of strongly interacting matter. While well-understood near the QCD crossover ($T \sim 150$ MeV) and in the low-temperature hadronic regime, the EoS at high temperatures ($T \gg 1$ GeV) remains a challenge.

• Limitations of Perturbation Theory: At high temperatures, standard perturbative expansions in the strong coupling constant $g$ converge poorly. Even at temperatures approaching the electroweak scale ($T \sim 100$ GeV), the expansion becomes sensitive to ultrasoft modes that are inherently non-perturbative. Previous studies suggest that perturbative predictions deviate significantly from the true result even near the electroweak scale.
• Limitations of Lattice QCD: Traditional lattice determinations of the EoS have been limited to $T \lesssim 1$ GeV. At higher temperatures, the "window problem" arises: the lattice must be fine enough to resolve the thermal scale ($T^{-1}$) while large enough to contain hadronic scales ($M_{had}^{-1}$), requiring computationally prohibitive volumes. Furthermore, standard methods (like the integral method) require subtracting large ultraviolet divergences, which becomes numerically unstable when comparing widely separated temperature scales.

2. Methodology

The authors present a fully non-perturbative lattice QCD computation covering the temperature range 3 GeV to 165 GeV for $N_f=3$ massless quark flavors. The strategy relies on two key innovations:

A. Lines of Constant Physics (LCPs) via Finite-Volume Schemes

To address the scale-setting problem at high energies without relying on hadronic observables (which are irrelevant at $T \gg M_{had}$), the authors define LCPs using a renormalized coupling computed in a finite volume.

• They utilize the Schrödinger Functional (SF) scheme.
• The renormalization scale is set by the inverse box size ($\mu = 1/L_0$).
• By following the non-perturbative running of the coupling $\bar{g}^2_{SF}$, they define a sequence of bare parameters that connect the hadronic scale to the electroweak scale, allowing for precise control over the continuum limit.

B. Shifted Boundary Conditions (Moving Reference Frame)

To avoid the numerical difficulties of the integral method (trace anomaly subtraction), the authors formulate QCD in a moving reference frame using shifted boundary conditions in the Euclidean time direction.

• Setup: Periodic boundary conditions are applied in spatial directions. In the temporal direction, fields satisfy $U_\mu(x_0 + L_0, \mathbf{x}) = U_\mu(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$, where $\boldsymbol{\xi}$ is a constant spatial shift vector.
• Physical Interpretation: This setup is equivalent to a Lorentz boost in Euclidean space. The effective temperature is $T = \gamma / L_0$, where $\gamma = 1/\sqrt{1+\xi^2}$.
• Advantage: This geometric setup allows for the direct determination of the entropy density ($s$) via the derivative of the free energy with respect to the shift parameter $\xi$, eliminating the need for temperature-dependent vacuum subtractions.
  $$\frac{s}{T^3} = \frac{1+\xi^2}{\xi_k T^4} \frac{\partial f_\xi}{\partial \xi_k}$$

C. Numerical Implementation

The entropy density calculation is decomposed into two parts:

1. Pure Gauge Contribution ($SU(3)$): Evaluated by integrating the expectation value of the gauge action density over the bare coupling $g_0$.
2. Fermionic Contribution ($N_f=3$): Evaluated by integrating the scalar density $\langle \bar{\psi}\psi \rangle$ over the quark mass from the chiral limit to the static limit.

• Simulations were performed on lattices with spatial size $L/a = 144$ and temporal resolutions $L_0/a = 4, 6, 8, 10$.
• Continuum extrapolation was performed using a global fit to data at $L_0/a = 6, 8, 10$, removing $O(a)$ and $O(g^2)$ artifacts.

3. Key Contributions

• Unprecedented Temperature Range: The first non-perturbative determination of the QCD EoS extending from 3 GeV up to 165 GeV (the electroweak scale).
• Methodological Framework: A robust, generalizable framework combining shifted boundary conditions with finite-volume renormalization (SF scheme) to tackle high-temperature thermodynamics without hadronic inputs.
• Precision: Results are achieved with a relative error of 0.5% – 1.0% across the entire temperature range.

4. Results

• Entropy Density ($s/T^3$): The continuum-extrapolated entropy density shows a smooth approach to the Stefan-Boltzmann limit but reveals significant deviations from perturbative predictions.
• Breakdown of Perturbation Theory:
  • Comparison with standard perturbation theory (up to $O(g^6 \ln g^2)$) and Hard Thermal Loop (HTL) resummation shows that neither captures the non-perturbative data accurately, even at $T \sim 100$ GeV.
  • The quadratic coefficient in the perturbative expansion ($s_2$) deviates significantly from the lattice data ($-5.1(9)$ vs. $-8.438$).
• Non-Perturbative Contributions: The data can only be accurately described by including terms beyond the known perturbative series (specifically $O(g^6)$ and higher), which contain intrinsic non-perturbative contributions.
• Phenomenological Fit: The authors provide a parametrization valid for $T \ge 500$ MeV that incorporates known perturbative coefficients up to $O(g^6)$ and fits the unknown higher-order terms using the lattice data. This fit smoothly connects to low-temperature lattice results.

5. Significance

• Cosmological Implications: The QCD EoS governs the expansion rate of the early Universe and the evolution of the effective number of relativistic degrees of freedom ($N_{eff}$) between the electroweak scale and the QCD crossover. Accurate knowledge of this EoS is crucial for predicting the spectrum of primordial gravitational waves and understanding Big Bang Nucleosynthesis (BBN).
• Validation of Non-Perturbative Dynamics: The study confirms that complex non-perturbative dynamics persist up to the electroweak scale, challenging the assumption that the system becomes purely perturbative at such high energies.
• Future Applications: The methodology is general and can be extended to QCD with $N_f=4$ and $N_f=5$ massive quark flavors, paving the way for a complete non-perturbative description of the Standard Model's thermodynamic evolution in the early Universe.
• Transport Coefficients: The framework sets the stage for calculating transport coefficients (e.g., shear viscosity) by determining the non-perturbative renormalization constants of the energy-momentum tensor.