The polarization of thermal dileptons emitted in high-energy heavy-ion collisions

This paper presents a comprehensive framework combining next-to-leading-order spectral functions and iEBE-MUSIC hydrodynamic simulations to demonstrate that thermal dilepton polarization in LHC Pb+Pb collisions is a sensitive probe of quark-gluon plasma properties, revealing its dependence on collision frames, pre-equilibrium effects, and establishing a direct mapping between dielectron and dimuon polarization.

The Gist

The Big Picture: A Cosmic Soup and Its "Glow"

Imagine smashing two heavy atoms (like lead) together at nearly the speed of light. This collision creates a tiny, incredibly hot, and dense drop of "soup" called a Quark-Gluon Plasma (QGP). This is the state of matter that existed just microseconds after the Big Bang.

Usually, scientists study this soup by looking at the "debris" (particles like protons and neutrons) that fly out when the soup cools down. But this paper focuses on something different: light that escapes the soup while it is still hot.

Specifically, the authors are looking at dileptons. Think of a dilepton as a pair of particles (like an electron and a positron, or a muon and an antimuon) that are born together from a "ghostly" flash of light (a virtual photon) inside the soup. Because these particles don't interact much with the soup, they fly straight out, carrying a perfect snapshot of what the soup was like at the exact moment they were born.

The Main Discovery: The "Polarization" of the Light

The paper isn't just about how many of these pairs are made; it's about how they are oriented.

The Analogy: The Spinning Top
Imagine the virtual photon (the parent of the dilepton pair) is like a spinning top.

• Polarization is the direction the top is spinning or leaning.
• In a calm, still room, the tops might spin in random directions.
• But in this "soup," the fluid is rushing and swirling. The authors found that the direction the tops lean (their polarization) is heavily influenced by the flow and movement of the soup itself.

The paper calculates exactly how this "leaning" changes based on the speed of the soup and the energy of the particles. They found that the polarization acts like a sensitive compass, pointing out the internal properties of the plasma that other measurements miss.

The Tools: A High-Definition Simulation

To figure this out, the authors built a massive computer simulation.

1. The Engine (Hydrodynamics): They used a model called iEBE-MUSIC to simulate the explosion. Think of this as a high-end video game engine that tracks every tiny drop of the soup as it expands, cools, and swirls.
2. The Physics (NLO): They didn't just use the basic rules of physics. They used "Next-to-Leading Order" (NLO) calculations.
   • Analogy: If a basic calculation is like a sketch of a car, the NLO calculation is like a 3D blueprint that includes the engine, the tires, and the air resistance. It accounts for complex interactions, like when a "gluon" (a particle that holds the soup together) bumps into a quark and changes the outcome.

Key Findings in Plain English

1. The "Frame of Reference" Matters
The authors looked at the polarization from different "camera angles" (called frames).

• The Helicity Frame (HX): Imagine looking at the spinning top from the side.
• The Collins-Soper Frame (CS): Imagine looking at it from a different angle, perhaps from the direction of the colliding beams.
• The Result: The polarization looks very different depending on which angle you choose. However, the authors discovered a special mathematical combination of these angles that stays the same no matter how you look at it. This is a "universal truth" about the soup that doesn't depend on your viewpoint.

2. The "Early Morning" vs. "Late Night" Soup
The soup changes over time.

• Pre-equilibrium (The "Early Morning"): Right after the collision, before the soup settles into a smooth flow, it's chaotic. The authors modeled this chaotic phase and found that the dileptons born here have a very strong polarization signal.
• Hydrodynamic Stage (The "Late Night"): As the soup flows smoothly, the signal changes.
• The Takeaway: By measuring the polarization of the particles, scientists might be able to tell if they are seeing the "chaotic early morning" or the "smooth late night" of the collision.

3. Electrons vs. Muons: The Same Story
The paper looked at two types of particle pairs: electrons (lightweight) and muons (heavier).

• The Result: Even though muons are heavier, the "leaning" (polarization) of the muon pairs is mathematically locked to the electron pairs. If you know how the electrons are leaning, you can perfectly predict how the muons are leaning. This is a strict "one-to-one" rule.

4. The "Background Noise"
At very high energies, there is another source of these particle pairs called the Drell-Yan process (created by hard collisions at the very start). The authors showed that this background noise has a different polarization signature than the thermal soup. This helps scientists separate the "signal" (the soup) from the "noise" (the initial crash).

Summary

This paper is a theoretical guidebook for future experiments. It tells scientists:

• "If you measure the direction (polarization) of these particle pairs, you can learn about the flow and temperature of the Quark-Gluon Plasma."
• "Don't just count the particles; look at how they are oriented."
• "We have calculated exactly how this works using the most advanced physics tools available, so when you look at data from the Large Hadron Collider (LHC), you will know what to expect."

In short, they turned the "spin" of these escaping particles into a new way to take a temperature and flow measurement of the universe's hottest, densest matter.

Technical Summary

Technical Summary: The Polarization of Thermal Dileptons Emitted in High-Energy Heavy-Ion Collisions

Problem Statement
Relativistic heavy-ion collisions (HICs) generate a quark-gluon plasma (QGP), a state of matter where quarks and gluons are deconfined. While hadronic observables primarily reflect the late-stage conditions of the collision, electromagnetic (EM) radiation, such as photons and dileptons, escapes the medium almost unscathed, carrying direct information about the early-stage temperature and medium properties. Although the spectra of thermal dileptons have been extensively studied, their polarization observables remain less explored despite offering unique insights into the QGP's in-medium properties, such as momentum anisotropy and external electromagnetic fields. This work addresses the need for a comprehensive theoretical framework to calculate thermal dilepton polarization at next-to-leading-order (NLO) in the strong coupling, specifically within the dynamic environment of heavy-ion collisions at Large Hadron Collider (LHC) energies.

Methodology
The authors develop a framework combining perturbative QCD calculations with state-of-the-art hydrodynamic simulations:

1. Theoretical Framework:

   • Spectral Functions: The calculation utilizes virtual photon spectral functions complete at NLO in the strong coupling ($\alpha_s$). This includes contributions from quark-antiquark annihilation, Compton scattering, and modified annihilation processes, as well as 1-loop corrections and Landau-Pomeranchuck-Migdal (LPM) resummation.
   • Polarization Formalism: The authors derive a Lorentz-covariant description of dilepton production rates and polarization. They define polarization coefficients ($\lambda_\theta, \lambda_\phi, \lambda_{\theta\phi}, \dots$) based on the angular distribution of the lepton pair in the virtual photon rest frame. The formalism accounts for the fluid velocity of the medium ($u^\mu$) and distinguishes between different reference frames, specifically the Helicity (HX) and Collins-Soper (CS) frames.
   • Invariant Observables: A rotationally invariant combination of polarization coefficients, $\tilde{\lambda}$, is defined and proven to be frame-independent even after averaging over the entire space-time evolution of the fireball.

2. Hydrodynamic Simulation:

   • Model: The space-time evolution of Pb+Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV is simulated using the iEBE-MUSIC framework.
   • Multi-Stage Approach:
     • Initial State: The IP-Glasma model (based on the Color Glass Condensate effective field theory) simulates the primordial evolution from $\tau = 0^+$ to $0.1$ fm.
     • Pre-equilibrium: The KøMPøST model simulates the pre-equilibrium stage ($0.1$ to $0.8$ fm) using pure-gluonic QCD effective kinetic theory.
     • Hydrodynamic Evolution: At $\tau = 0.8$ fm, the system is matched to relativistic viscous hydrodynamics (MUSIC) using Landau matching. The evolution includes temperature-dependent bulk viscosity and a constant shear viscosity to entropy density ratio ($\eta/s = 0.12$).
     • Hadronic Afterburner: The Cooper-Frye prescription with viscous corrections converts the fluid to hadrons, which are further processed by the UrQMD model.

3. Pre-equilibrium Dilepton Estimation:

   • Since a first-principles description of pre-equilibrium dileptons is not yet available, the authors employ a phenomenological approach. They estimate the pre-equilibrium dilepton production rate by applying a suppression factor ($SF$) to the equilibrium rates to account for the reduced quark abundance. Different suppression indices ($\alpha$) are tested for Leading Order (LO) and NLO channels to model the varying sensitivity of different production mechanisms to the quark density.

Key Contributions and Results

• NLO Effects on Polarization: The inclusion of NLO corrections significantly alters the predicted polarization coefficients compared to LO calculations. Specifically, the NLO spectral functions enhance the thermal dilepton yields in the low-mass region (LMR) and intermediate-mass region (IMR). In the HX frame, the inclusion of NLO processes leads to a visibly different $\lambda_\theta$ spectrum compared to LO.
• Frame Dependence and Flow Sensitivity:
  • The polarization coefficients $\lambda_\theta$ and $\lambda_\phi$ exhibit distinct behaviors in the HX and CS frames. For a static medium, these are directly correlated; however, in a realistic, expanding QGP, the CS frame polarization is found to be more sensitive to the longitudinal flow velocity of the medium.
  • The study confirms that the invariant combination $\tilde{\lambda}$ remains frame-independent even after integrating over the fireball's space-time history and transverse momentum ($p_T$), validating its utility as a robust observable.
• Pre-equilibrium Dynamics: The polarization in the IMR ($1 < M < 3$ GeV) is shown to be sensitive to the pre-equilibrium stage. The momentum dependence ($p_T$) of the polarization coefficient $\lambda_\theta$ in this mass window is particularly sensitive to the suppression of quark production during the pre-equilibrium phase. The "gluon effect" (helicity flipping via gluon emission in NLO channels) and the "temperature effect" (higher effective temperature in early stages) are identified as competing mechanisms influencing the net polarization.
• Drell-Yan vs. Thermal Contributions: In the high-mass region (HMR), Drell-Yan (DY) processes dominate. The paper demonstrates that DY and thermal contributions have opposite signs for the polarization quadrupole moment. While DY dominates the yield at $M > 4.5$ GeV, the thermal polarization signature remains considerable in the IMR.
• Dielectron-Dimuon Mapping: A strict one-to-one mapping is derived between dielectron ($e^+e^-$) and dimuon ($\mu^+\mu^-$) polarization coefficients. This relationship holds regardless of the local flow velocity or the choice of reference frame, provided the approximation $m_e \approx 0$ is used. This allows experimental measurements of one lepton species to be directly inferred for the other.

Significance
The paper establishes a comprehensive theoretical baseline for interpreting future dilepton polarization measurements in heavy-ion collisions. By integrating NLO spectral functions with multi-stage hydrodynamic simulations, the authors demonstrate that polarization observables are sensitive probes of:

1. In-medium properties: Specifically, the spectral function difference $\rho_\Delta = \rho_T - \rho_L$, which vanishes in vacuum but is non-zero in the medium.
2. Early-stage dynamics: The sensitivity of the polarization spectrum to pre-equilibrium quark abundance and the transition from a gluon-dominated to a chemically equilibrated system.
3. Flow anisotropy: The differential sensitivity of HX and CS frames to the longitudinal expansion of the QGP.

The work provides a necessary theoretical tool for upcoming experimental analyses at the LHC, offering a way to disentangle thermal emission from hard scattering processes (Drell-Yan) and to extract information about the QGP's early evolution that is inaccessible through invariant mass spectra alone. The derived dielectron-dimuon correspondence further aids in cross-validating experimental results across different detector capabilities.