Uncertainty quantification of holographic transport and energy loss for the hot and baryon-dense QGP

This paper presents an open-source C++ implementation of a holographic Einstein-Maxwell-Dilaton model with improved numerical stability to quantify uncertainties in various transport coefficients and energy loss parameters across the QCD phase diagram by propagating lattice QCD constraints via Bayesian sampling.

The Gist

Imagine you are trying to understand the behavior of a substance so strange and hot that it doesn't exist anywhere else in the universe today. This substance is called Quark-Gluon Plasma (QGP). It's the "soup" that existed just a fraction of a second after the Big Bang, before the universe cooled down enough to form protons and neutrons.

Scientists create tiny droplets of this soup by smashing heavy atoms together at nearly the speed of light in giant particle accelerators like the Large Hadron Collider (LHC). But here's the problem: this soup is incredibly difficult to study. It's too hot, too dense, and it disappears in a blink of an eye.

This paper is like a super-advanced weather forecast for that cosmic soup, but with a twist: it predicts what happens when you squeeze the soup with extra pressure (adding "baryon density," which is like adding more matter to the mix).

Here is a simple breakdown of what the authors did, using everyday analogies:

1. The Problem: We Can't See the Whole Picture

Scientists have two main ways to study this soup:

• The "Real World" Lab (Lattice QCD): They run supercomputer simulations based on the laws of physics. It's like taking a high-resolution photo of the soup, but the camera only works well when the soup is "thin" (low pressure). When they try to take a picture of the soup when it's "thick" and dense, the camera gets blurry and the math breaks down.
• The "Theoretical Map" (Holography): They use a mathematical trick called "Holography." Imagine you have a 3D object (the soup), but you can only see its 2D shadow on a wall. By studying the shadow (which is easier to calculate), you can figure out what the 3D object is doing. This works great for dense soup, but it's a bit of a guesswork map.

2. The Solution: Merging the Photo and the Map

The authors of this paper decided to combine the best of both worlds. They took their "Holographic Map" and forced it to match the "Real World Photos" where the photos were clear (low pressure).

• The Analogy: Imagine you are trying to draw a map of a mountain range. You have a very accurate satellite photo of the base of the mountain, but the top is covered in fog. You also have a rough sketch of the whole mountain.
  • They took their rough sketch (the Holographic model) and adjusted it until the base matched the satellite photo perfectly.
  • Once the base was locked in, they used the logic of the sketch to predict what the foggy top (the dense, high-pressure region) looks like.

3. The New Tool: A Better Ruler

One of the biggest headaches in this field is that the math gets "noisy" and unstable when you try to calculate things at the edges of the universe (the "boundary" in their math). It's like trying to measure the temperature of a fire with a ruler that melts when it gets too close.

• The Innovation: The team invented a new numerical method (a new way of doing the math) that acts like a heat-resistant, self-calibrating ruler. Instead of guessing the temperature at the edge, they used a "relaxation" technique that smooths out the noise. This made their calculations much faster and much more stable, allowing them to run thousands of simulations to get a reliable average.

4. The "Bayesian" Safety Net

Because the math involves many variables, there's always a chance of being slightly wrong. To handle this, they used a statistical method called Bayesian Inference.

• The Analogy: Imagine you are trying to guess the weight of a mystery box. You don't just guess one number; you ask 1,000 experts. Some say 5kg, some say 6kg. Instead of picking one, you create a "confidence band" (a range) that says, "We are 95% sure the weight is between 5.2 and 5.8 kg."
• The authors did this with physics. They didn't just give one prediction for how the soup flows; they gave a range of possibilities with error bars, showing exactly how confident they are in their results.

5. What Did They Find? (The Weather Report)

They calculated several "transport coefficients," which are fancy words for "how the soup moves and reacts." Here is what they found:

• Viscosity (Stickiness): They found that as you squeeze the soup (add more baryon density), it becomes less sticky and flows more like a perfect liquid. It's like honey that suddenly turns into water when you press it hard. This suggests the soup becomes even more "perfect" under pressure.
• Energy Loss (The Jet Quenching): When a fast particle (like a jet of energy) flies through the soup, it slows down. They found that in the dense soup, the soup is thicker and stickier to these fast particles, causing them to lose energy much faster. It's like running through a crowd; the denser the crowd, the harder it is to get through.
• The Critical Point: They mapped out where the soup might undergo a sudden phase change (like water turning to ice, but for plasma). They found a "Critical Point" where the soup behaves strangely, and they predicted exactly where this happens in terms of temperature and pressure.

6. Why Does This Matter?

This research helps scientists interpret data from experiments like the Beam Energy Scan at RHIC (a particle collider in the US) and the CBM experiment planned in Germany.

• The Big Picture: By knowing exactly how this "perfect fluid" behaves under different pressures, scientists can better understand the history of our universe (how it cooled down after the Big Bang) and the physics inside neutron stars (which are essentially giant balls of this dense matter).

Summary

In short, this paper built a super-accurate, 3D simulation of the universe's hottest, densest soup. They fixed the math to make it stable, used real-world data to calibrate it, and provided a detailed "weather forecast" showing how this soup flows, how it resists movement, and how it loses energy under extreme pressure. It's a major step forward in understanding the fundamental building blocks of our reality.

Technical Summary

1. Problem Statement

Understanding the properties of the Quark-Gluon Plasma (QGP) across the full QCD phase diagram, particularly in the regime of high baryon density ($\mu_B$), remains a central challenge in nuclear physics.

• Limitations of Lattice QCD: While Lattice QCD (LQCD) provides precise results at zero baryon chemical potential, it suffers from the "sign problem" at finite $\mu_B$, making first-principles calculations of the equation of state (EoS) and transport coefficients unreliable in the baryon-rich regime.
• Theoretical Gaps: Effective models often lack quantitative constraints, while standard holographic (AdS/CFT) models, though capable of describing strongly coupled systems, often rely on analytical approximations or lack rigorous statistical uncertainty quantification when extrapolating to finite density.
• Specific Need: There is a need for a framework that can systematically propagate LQCD uncertainties (at $\mu_B=0$) to predictions of transport coefficients (viscosities, conductivities, energy loss) across the crossover, the predicted Critical End Point (CEP), and the first-order phase transition line at finite $\mu_B$.

2. Methodology

The authors employ a Holographic Einstein-Maxwell-Dilaton (EMD) model combined with Bayesian inference and novel numerical techniques.

A. Theoretical Framework

• EMD Model: A bottom-up holographic model based on gauge-gravity duality. The bulk action includes gravity, a Maxwell field (dual to baryon number), and a dilaton field (breaking conformal symmetry).
• Functional Forms: Two distinct Ansätze for the dilaton potential $V(\phi)$ and the Maxwell-dilaton coupling $f(\phi)$ are used:
  1. Poly-Hyperbolic Ansatz (PHA): A polynomial-hyperbolic form.
  2. Parametric Ansatz (PA): A form utilizing hyperbolic tangents.
• Bayesian Inference: The model parameters are constrained using LQCD data for entropy density and baryon susceptibility at $\mu_B = 0$ (from the $N_f = 2+1$ flavor continuum limit). This generates a posterior distribution of model parameters.

B. Novel Numerical Schemes

To enable efficient Bayesian sampling and high-precision transport calculations, the authors introduced two key algorithmic improvements in their C++ implementation (within the MUSES Framework):

1. Relaxational UV Extraction: A new method to extract ultraviolet (UV) asymptotic coefficients (specifically the dilaton coefficient $\phi_A$) from near-boundary data. Instead of noisy least-squares fitting, they introduce an auxiliary field $C(r)$ governed by a relaxational differential equation. This averages out numerical noise, significantly improving stability and computational speed.
2. CEP Location Algorithm: A fast iterative algorithm to locate the Critical End Point (CEP) by identifying intersections of constant $\phi_0$ lines (where $\phi_0$ is the horizon value of the dilaton) in the $(T, \mu_B)$ plane. This allows for the rapid mapping of the phase diagram for every sample in the Bayesian posterior.

C. Transport Coefficients Calculation

Using the posterior samples, the authors computed a wide array of observables:

• Thermodynamics: Baryon susceptibility, entropy density, pressure.
• Transport: Baryon conductivity ($\sigma_B$), thermal conductivity ($\sigma_T$), baryon diffusion ($D_B$), shear viscosity ($\eta$), and bulk viscosity ($\zeta$).
• Energy Loss: Heavy quark drag force, Langevin diffusion coefficients ($\kappa_\parallel, \kappa_\perp$), and the jet-quenching parameter ($\hat{q}$).

3. Key Contributions

1. First Bayesian Uncertainty Bands at Finite $\mu_B$: This is the first study to provide fully quantified Bayesian uncertainty bands for a broad set of transport coefficients across the entire phase diagram (crossover, CEP, and first-order line) constrained by LQCD data.
2. Numerical Innovation: The introduction of the relaxational method for UV asymptotics solves long-standing numerical stability issues in holographic transport calculations, enabling the large-scale sampling required for Bayesian analysis.
3. Comprehensive Observable Set: The work extends beyond previous studies (which often focused only on thermodynamics or zero density) to include thermal conductivity, baryon diffusion, and heavy-quark energy loss at finite density.
4. Robustness Check: By comparing results from two different functional Ansätze (PHA and PA), the authors demonstrate that the qualitative features of the transport coefficients are robust against the specific choice of the holographic potential.

4. Key Results

A. Phase Diagram and Critical Point

• The Bayesian analysis confirms the existence of a CEP. The posterior distribution tightly clusters the CEP location within $T \approx 101\text{--}108$ MeV and $\mu_B \approx 560\text{--}625$ MeV (95% confidence level).
• The phase diagram exhibits a crossover at low $\mu_B$, terminating at the CEP, followed by a first-order phase transition line at higher $\mu_B$.

B. Transport Coefficients

• Shear Viscosity ($\eta/s$): Consistent with the universal holographic bound $\eta/s = 1/4\pi$. However, the dimensionless ratio $\eta T / (\epsilon + P)$ decreases with increasing $\mu_B$, indicating the medium becomes "more perfect" (closer to a perfect fluid) in the baryon-dense regime.
• Bulk Viscosity ($\zeta$): Shows a peak near the crossover region. The peak height is suppressed as $\mu_B$ increases, and the peak shifts toward the CEP. The results show good agreement with independent Bayesian analyses of heavy-ion data (JETSCAPE Collaboration).
• Baryon Transport:
  • Conductivity ($\sigma_B$): Increases with temperature. At fixed $T$, it increases with $\mu_B$ at low temperatures but decreases at high temperatures.
  • Diffusion ($D_B$): Generally suppressed by increasing baryon density. At the CEP, $D_B \to 0$ because the baryon susceptibility diverges while conductivity remains finite.
  • Thermal Conductivity: Exhibits a minimum near the crossover, which evolves into a cusp at the CEP.
• Energy Loss (Heavy Quarks & Jets):
  • Drag Force & Diffusion: Both the heavy quark drag force and Langevin diffusion coefficients ($\kappa_\parallel, \kappa_\perp$) increase significantly with baryon density. This suggests that partons lose energy more rapidly in baryon-rich matter.
  • Jet Quenching ($\hat{q}$): The jet-quenching parameter also increases with $\mu_B$ and temperature, indicating stronger jet suppression in the baryon-dense regime.
  • Critical Behavior: All energy loss observables exhibit a peak or cusp near the CEP and a discontinuity gap across the first-order phase transition line.

C. Comparison with Data

• At $\mu_B = 0$, the holographic predictions for bulk viscosity and jet-quenching parameters align well with the Bayesian constraints derived from heavy-ion collision data by the JETSCAPE Collaboration and the Duke Group.
• The agreement validates the phenomenological applicability of the EMD model for describing the QGP.

5. Significance

• Phenomenological Impact: The results provide essential input for hydrodynamic simulations of heavy-ion collisions at lower beam energies (e.g., RHIC Beam Energy Scan, FAIR/CBM), where baryon density is significant. Specifically, the suppression of diffusion and enhancement of energy loss in dense matter will affect observables like proton/antiproton flow and jet suppression patterns.
• Methodological Benchmark: The paper sets a new standard for holographic modeling by integrating rigorous statistical uncertainty quantification. It demonstrates that holographic models can be systematically constrained by LQCD to make reliable predictions in regimes where LQCD fails.
• Critical Point Physics: The identification of sharp features (peaks, divergences, discontinuities) in transport coefficients near the CEP offers potential experimental signatures for the search of the QCD critical point.

In summary, this work bridges the gap between first-principles lattice data and phenomenological models of dense QCD matter, providing a robust, uncertainty-quantified map of how the QGP's transport properties evolve from the crossover region into the baryon-rich phase transition.