Constraining hot and cold nuclear matter properties from heavy-ion collisions and deep-inelastic scattering

This paper presents a global analysis combining deep-inelastic scattering and heavy-ion collision data within a saturation-based QCD framework to constrain the early-time shear viscosity to entropy density ratio ($\eta/s$) of the quark-gluon plasma.

The Gist

Imagine the universe as a giant kitchen. Usually, the ingredients (atoms) are like solid, frozen blocks of ice. But if you turn the heat up to an unimaginable degree—like the temperature inside a star or the moment right after the Big Bang—those blocks melt into a super-hot, super-dense soup. Physicists call this soup Quark-Gluon Plasma (QGP). It's a state of matter where the tiny building blocks of protons and neutrons (quarks and gluons) are free to swim around instead of being stuck together.

This paper is like a team of detectives trying to figure out how "thick" or "runny" this cosmic soup is. In physics, this "thickness" is called viscosity. If the soup is very runny (low viscosity), it flows easily. If it's thick (high viscosity), it resists flowing. Knowing this helps scientists understand how the universe behaved in its very first moments.

Here is how the authors solved the mystery, using a step-by-step detective story:

1. The Three Clues (The Data)

To figure out the properties of this soup, the team didn't just look at one thing. They combined three different types of clues, like a detective cross-referencing a fingerprint, a witness statement, and a security camera:

• Clue A: The "Cold" Snapshot (HERA): They looked at data from smashing electrons into protons (deep-inelastic scattering). Think of this as taking a high-speed photo of a single, cold proton to understand its internal structure before it gets smashed. This tells them how the "ingredients" are packed when things are calm.
• Clue B: The "Small" Smashes (p+p and p+Pb): They looked at collisions where a proton hits another proton or a light lead nucleus. These are like small-scale experiments that help them calibrate their measuring tools without the soup getting too messy.
• Clue C: The "Big" Smashes (Pb+Pb): Finally, they looked at heavy lead nuclei smashing into each other at the Large Hadron Collider (LHC). This is where the real "soup" is made. They measured how many particles came out of the crash.

2. The Recipe (The Model)

The team used a theoretical "recipe" based on a concept called Color Glass Condensate (CGC).

• The Analogy: Imagine the proton isn't a solid ball, but a fuzzy cloud of tiny, fast-moving gluons (like a swarm of bees). When you smash two of these clouds together, the bees get squished and the energy explodes.
• The authors built a computer model that simulates this explosion. They started with the "cold" snapshot (Clue A) to set the initial conditions, then used the "small" smashes (Clue B) to adjust the scale of the explosion (a factor they call K).

3. The Shortcut (The Estimator)

Simulating the entire explosion of a heavy-ion collision is incredibly hard and slow, like trying to simulate every single water molecule in a tsunami.

• The Trick: The team realized that the number of particles produced (the "multiplicity") is directly linked to how much energy was dumped into the soup at the start.
• They created a shortcut formula. Instead of running a full, slow simulation every time, they used this formula to estimate the final result based on the initial energy. They "calibrated" this shortcut by running a few full simulations first to make sure the math worked.

4. The Big Reveal (The Results)

By combining all these clues and running their model against the real data from the ALICE experiment at the LHC, they found the answer to the "thickness" question.

• The Viscosity: They determined the ratio of viscosity to entropy (a measure of disorder) for this early-stage soup. Their result is 0.31.
  • What does this mean? It suggests the quark-gluon plasma is a very "perfect" fluid—extremely runny, almost like a super-fluid. It flows with very little resistance.
• The Temperature: They also estimated the temperature of this soup during the very early, chaotic phase. It's incredibly hot, around 500 MeV (which is roughly 5.8 trillion degrees Celsius).

Why This Matters

The authors emphasize that this is a "proof-of-principle" study. They showed that you can figure out the properties of this extreme, hot matter by carefully connecting the dots between cold proton data, small collisions, and big collisions.

They found that their result (0.31) matches up well with other theoretical predictions from super-computers (Lattice QCD) and high-energy math (perturbative QCD). This gives them confidence that their model of the early universe is on the right track.

In short: The team built a bridge between the cold, quiet world of single protons and the hot, chaotic world of heavy-ion collisions. By crossing that bridge, they measured the "runniness" of the universe's first soup, finding it to be an incredibly fluid substance.

Technical Summary

Technical Summary: Constraining Hot and Cold Nuclear Matter Properties from Heavy-Ion Collisions and Deep-Inelastic Scattering

Problem Statement
Understanding the transport properties of the quark-gluon plasma (QGP), specifically the shear viscosity to entropy density ratio ($\eta/s$), is central to characterizing strongly interacting matter under extreme conditions. While hydrodynamic models successfully describe the late-stage evolution of heavy-ion collisions, extracting $\eta/s$ during the early-time, non-equilibrium pre-equilibrium phase remains challenging. This phase is sensitive to the initial state of the collision, which is governed by high-energy QCD dynamics, including gluon saturation. A significant gap exists in simultaneously constraining the properties of "cold" nuclear matter (probed via Deep-Inelastic Scattering, DIS) and "hot" nuclear matter (probed via heavy-ion collisions) within a unified framework. Previous approaches often treat the initial state energy density as a free parameter or rely on complex, computationally expensive dynamical simulations for every data point in a global fit.

Methodology
The authors perform a global analysis combining data from three distinct regimes:

1. Deep-Inelastic Scattering (DIS): Inclusive $e+p$ scattering data from HERA.
2. Small Systems: Transverse energy distributions in $p+p$ and $p+Pb$ collisions.
3. Large Systems: Charged hadron multiplicities in central $Pb+Pb$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV from ALICE.

The theoretical framework is based on the Color Glass Condensate (CGC) effective field theory. The initial state is modeled using the Golec-Biernat and Wüsthoff (GBW) saturation model, extended to include transverse position dependence and nucleon-level fluctuations.

The analysis proceeds in a multi-step fitting procedure:

• Step 1 (Cold Matter): A Bayesian analysis of HERA DIS data constrains the parameters of the proton saturation scale ($Q_{s,0}$ and $\lambda$). A Gaussian Process (GP) emulator is trained to efficiently map these parameters to observables.
• Step 2 (Initial State Calibration): The model is applied to $p+p$ and $p+Pb$ collisions. By comparing the model's predicted transverse energy per unit rapidity ($dE_\perp/dy$) with ALICE measurements, a global $K$-factor is determined to account for higher-order corrections in the initial energy deposition.
• Step 3 (Hot Matter & Pre-equilibrium): For $Pb+Pb$ collisions, the authors utilize a multiplicity estimator derived from the initial transverse energy density. This estimator relates the initial energy to the final charged particle multiplicity ($dN_{ch}/d\eta$) via the pre-equilibrium entropy production. To avoid running full dynamical simulations for every fit iteration, the authors calibrate this estimator using a hybrid simulation chain (KøMPøST pre-equilibrium + MUSIC hydrodynamics + iSS particlization + SMASH hadronic afterburner). This calibration yields a proportionality factor $C(\eta/s)$ that corrects the analytical estimator.
• Step 4 (Extraction): The calibrated estimator is used to extract $\eta/s$ by fitting the model to the most central ($0-2.5\%$) $Pb+Pb$ multiplicity data, varying $\eta/s$ while propagating uncertainties from the DIS and small-system fits.

Key Results

• Parameter Constraints: The DIS fit yields $Q_{s,0} = 0.393 \pm 0.010$ GeV and $\lambda = 0.219 \pm 0.007$. The initial state $K$-factor is determined to be $K = 1.93 \pm 0.04$.
• Shear Viscosity: The extracted value for the non-equilibrium shear viscosity to entropy density ratio is $\langle \eta/s \rangle_{\tau_{hydro}} = 0.31 \pm 0.05$. This value is found to be insensitive to the switching time ($\tau_{hydro}$) between the pre-equilibrium and hydrodynamic stages (tested in the range $0.4$ to $1.0$ fm).
• Effective Temperature: By analyzing the entropy production during the pre-equilibrium stage, the authors estimate an effective temperature scale of $\langle T_{eff} \rangle \approx 470 - 560$ MeV, depending on the switching time.
• Model Performance: The model provides a statistically consistent description of the HERA data ($\chi^2/\text{d.o.f.} = 1.02$) and reproduces the charged hadron multiplicity in $Pb+Pb$ collisions well, particularly in the most central bin.

Significance and Claims
The paper presents a "proof-of-principle" study demonstrating that a saturation-based initial state model, constrained by cold nuclear matter data (DIS) and small-system data, can be used to extract transport properties of hot nuclear matter without treating the initial energy density as a free parameter.

The authors claim that their approach offers a computationally efficient alternative to full Bayesian analyses of heavy-ion collisions, which typically require thousands of event-by-event dynamical simulations. By decoupling the initial state constraints (from DIS and $p+p/p+Pb$) from the final state extraction (via a calibrated estimator), they achieve tight constraints on $\eta/s$ at high temperatures.

The extracted value of $\eta/s \approx 0.31$ at an effective temperature of $\sim 0.5$ GeV is consistent with values from Lattice QCD and perturbative QCD (NLO). The authors emphasize that because their extraction occurs in the pre-equilibrium phase (where the system is not yet in thermodynamic equilibrium), this value represents an upper bound for $\eta/s$ at these temperatures. They note that their result differs qualitatively from some other Bayesian extractions (e.g., Ref. [59]) which show a decreasing $\eta/s$ with temperature, highlighting the sensitivity of the extraction to the modeling of the pre-equilibrium phase.

The work establishes a pathway for future comprehensive analyses that could incorporate flow observables and temperature-dependent viscosities while retaining the rigorous constraints provided by DIS and small-system data.