Non-perturbative Renormalization of the EMT in Full QCD

✨

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

The Big Picture: Measuring the "Sticky" Universe

Imagine the universe, specifically the Quark-Gluon Plasma (QGP), which is the super-hot, super-dense soup of particles that existed just after the Big Bang. Scientists want to know how "sticky" or "fluid" this soup is. In physics, this stickiness is called viscosity.

To measure this, physicists need to calculate a specific mathematical object called the Energy-Momentum Tensor (EMT). Think of the EMT as a giant, complex dashboard that tells you exactly how much energy is moving, where it's going, and how much pressure it's exerting.

The Problem:
Calculating this dashboard on a computer is incredibly hard. The computer uses a grid (like a chessboard) to simulate space and time. However, nature doesn't actually use a grid; it's smooth and continuous. When you force nature onto a grid, you break some of the fundamental rules of symmetry (like how a circle looks the same no matter how you rotate it).

Because the grid breaks these rules, the numbers on our "dashboard" come out wrong. They are "renormalized" (a fancy word for "fixed") by adding invisible correction factors. In a simple world, you could figure out these corrections easily. But in the complex world of Full QCD (which includes both gluons and quarks), there are too many unknowns. It's like trying to solve a math problem with two missing numbers but only one equation.

The Solution: The "Imaginary Chemical Potential" Trick

The authors of this paper found a clever way to solve this "two unknowns, one equation" problem. They used a trick involving Imaginary Chemical Potential.

Here is the analogy:
Imagine you are trying to weigh a heavy backpack (the Gluons) and a light feather (the Quarks) together on a scale, but the scale is broken and only gives you the total weight. You can't tell how much is the backpack and how much is the feather.

To fix this, the scientists created a special "magic room" (a specific theoretical setup) where the feather suddenly becomes invisible or weighs almost nothing, while the backpack stays exactly the same.

The Setup: They used a concept called the Roberge-Weiss transition. In this specific setup, by tweaking a dial (the imaginary chemical potential), the quarks (the feather) are suppressed by a massive factor (about 100 times less effective).
The Result: Now, when they look at the scale, the reading is almost entirely due to the backpack (gluons).
The Comparison: They then compare this "feather-suppressed" reading with a normal reading where both the backpack and feather are present.
The Math: By looking at the difference between the two readings, they can mathematically isolate exactly how much the "feather" (quarks) was contributing. This allows them to finally solve for the missing correction factors (the renormalization coefficients).

The Journey: From Pure Glue to Full QCD

The paper describes a step-by-step journey:

Step 1 (Pure Glue): They first mastered the math for a world with only gluons (no quarks). This was like learning to ride a bike with training wheels.
Step 2 (Adding Quarks): They added quarks back in. Suddenly, the bike got wobbly because quarks introduce new, messy variables.
Step 3 (The Fix): They applied the "magic room" trick (imaginary chemical potential) to stabilize the bike.

The Findings: What They Discovered

After running massive simulations on supercomputers, they found some surprising things:

The Gluons are "Heavy": They expected the pressure from the gluons to be close to what a simple "free gas" theory predicts. Instead, they found the gluons are contributing much less pressure than expected. It's as if the backpack is heavier and more sluggish than anyone thought.
The Quarks are "Light": The quarks, on the other hand, are behaving very much like a free gas (the feather is behaving normally).
Success: They successfully calculated the correction factors needed to fix the EMT dashboard.

Why This Matters

This paper is a crucial stepping stone. Before this, scientists could calculate the viscosity (stickiness) of the "glue-only" universe, but not the "real" universe with quarks.

Now that they have figured out how to fix the dashboard (renormalize the EMT) in the full theory, they are ready to do the final calculation: determining the shear viscosity of the Quark-Gluon Plasma.

In short: They built a new, more accurate ruler to measure the "stickiness" of the universe's first moments. By using a clever trick to separate the ingredients, they finally have the tools to understand how the early universe flowed.

1. Problem Statement

The Energy-Momentum Tensor (EMT) is a fundamental conserved current associated with space-time translation symmetry, essential for calculating transport coefficients like shear viscosity ( $\eta$ ) and bulk viscosity ( $\zeta$ ) in the quark-gluon plasma (QGP). While the EMT is well-defined in the continuum limit of Quantum Chromodynamics (QCD), its formulation on the lattice is non-trivial due to the breaking of continuous $SO(4)$ rotational and translational symmetries by the discrete lattice regulator.

The Core Challenge: On the lattice, the EMT operators are not uniquely defined. Unlike the continuum where Ward identities fix coefficients, the lattice requires distinct renormalization coefficients for each irreducible representation (irrep) of the symmetry group.
Specific Difficulty in Full QCD: While renormalization in pure-gauge (gluonic) theory has been established, extending this to $(2+1)$ -flavor QCD (including dynamical fermions) introduces significant complexity. The inclusion of fermions adds new operator structures, creating an underdetermined system where the number of unknown renormalization coefficients exceeds the number of available physical constraints (equations).

2. Methodology

The authors employ Gradient Flow to construct the EMT operators and a novel Imaginary Chemical Potential strategy to resolve the underdetermined system.

A. Gradient Flow Formalism

The EMT is constructed using the gradient flow method, which smears operators over a radius $\sqrt{8\tau_F}$ (where $\tau_F$ is the flow time). This effectively suppresses UV divergences, rendering the operators finite at finite flow time.

Pure-Gauge Sector: The EMT decomposes into two irreps ( $T^1_{\mu\nu}$ and $T^2_{\mu\nu}$ ) with renormalization coefficients $c_1$ and $c_2$ .
Full QCD Sector: The inclusion of dynamical fermions introduces three additional irreps ( $T^3_{\mu\nu}, T^4_{\mu\nu}, T^5_{\mu\nu}$ $T_{μν}^{3}, T_{μν}^{4}, T_{μν}^{5}$ ), resulting in five operators with coefficients $Z_1$ $Z_{1}$ through $Z_5$ $Z_{5}$ .
- $Z_1$ and $Z_3$ are the critical coefficients for the traceless components relevant to shear viscosity.
- $Z_4$ and $Z_5$ are linked via equations of motion and can be fixed via scale dependence.

B. The Underdetermined System

To determine the renormalization coefficients, the authors utilize the enthalpy density ( $\epsilon + p$ ) as a physical constraint.

In pure-gauge theory, enthalpy density provides a single equation to solve for the single unknown $c_1$ (since the trace term cancels out).
In Full QCD, the enthalpy density equation (Eq. 11) contains two unknowns ( $Z_1$ and $Z_3$ ) but only one equation, making the system underdetermined.

C. The Solution: Imaginary Isospin Chemical Potential

To generate a second independent equation, the authors utilize Imaginary Chemical Potential ( $\mu_I$ ) to suppress the fermionic contribution to the pressure while keeping the gluonic contribution relatively stable.

Mechanism: They exploit the Roberge-Weiss (RW) transition and $Z_3$ center symmetry. By setting $\mu_u/T = -\mu_d/T = i\theta$ and $\mu_s/T = 0$ , they induce a phase transition at $\theta = 2\pi/3$ .
Effect: At this transition point, the fermionic contribution to the pressure is suppressed by a factor of $\approx 81$ (in the massless limit) or $>100$ (at physical masses), effectively isolating the gluonic sector.
Strategy: The authors perform measurements on two distinct ensembles:
1. Standard Ensemble: $\mu = 0$ (Full fermionic contribution).
2. RW Ensemble: $\mu_u/T = -\mu_d/T = i2\pi/3$ (Suppressed fermionic contribution).
  By comparing the enthalpy density between these two states, they can disentangle and solve for both $Z_1$ and $Z_3$ simultaneously.

3. Key Contributions

Extension to Full QCD: The paper provides the first detailed renormalization procedure for the EMT in $(2+1)$ -flavor QCD using gradient flow, moving beyond the previously solved pure-gauge case.
Novel Constraint Technique: The introduction of the imaginary isospin chemical potential to create a "fermion-suppressed" ensemble is a key methodological innovation for solving the underdetermined system of renormalization coefficients.
Multi-Beta Analysis: The study performs calculations at multiple lattice couplings ( $\beta = 7.373$ and $\beta = 7.596$ ) at approximately the same temperature ( $T \approx 400$ MeV) to facilitate a reliable continuum extrapolation.

4. Results and Analysis

The authors present results for number density, pressure, energy density, and enthalpy density at $T \approx 400$ MeV (significantly above the pseudo-critical temperature $T_{pc} \approx 155$ MeV).

Consistency: The thermodynamic quantities (number density, pressure, etc.) show high consistency between the two lattice spacings ( $\beta$ values), validating the setup for continuum extrapolation.
Pressure Suppression: As $\theta \to 2\pi/3$ $θ \to 2 π /3$ , the pressure drops significantly. The observed suppression factor is $\approx 10$ $\approx 10$ , which is stronger than the theoretical free-field prediction of $\approx 3$ $\approx 3$ .
- Analysis of Suppression: The authors find that the fermionic pressure approaches the Stefan-Boltzmann limit, but the gluonic pressure is significantly suppressed (far below the free-theory value). This suggests the gap between free and interacting pressure in QCD at $T=400$ MeV arises primarily from gluonic suppression, not quark interactions.
Enthalpy Density: The enthalpy density was measured without gradient flow (relying on its plateau region), while EMT operators were measured at finite flow time. The results are consistent across the different ensembles.
Current Status: The renormalization coefficients are being finalized. The primary source of uncertainty remains the zero-chemical-potential pressure contribution, which is being refined.

5. Significance and Outlook

Shear Viscosity Determination: The successful determination of the renormalization coefficients ( $Z_1, Z_3$ ) is the critical prerequisite for calculating the Euclidean two-point correlation functions of the EMT. This will enable the first-ever non-perturbative determination of shear viscosity in $(2+1)$ -flavor QCD on the lattice.
Theoretical Insight: The findings regarding the suppression of gluonic pressure provide new insights into the nature of the QGP at temperatures around 400 MeV, highlighting the dominance of non-perturbative gluonic effects over fermionic ones in this regime.
Future Work: The team is incorporating an additional lattice spacing ( $\beta=7.825$ ) to finalize the continuum limit and refine the error budget to complete the shear viscosity calculation.

In summary, this work bridges a critical gap in lattice QCD by solving the renormalization problem for the EMT in full QCD, paving the way for precise calculations of transport properties in the quark-gluon plasma.