Thermal static Potential at Finite Density in (2+1)-flavor QCD

This paper investigates the thermal static potential in (2+1)-flavor QCD at finite density using a Taylor expansion, revealing an enhancement of in-medium screening in both real and imaginary parts that offers new constraints on heavy-quark interactions relevant to the RHIC Beam Energy Scan and future FAIR experiments.

The Gist

Imagine you are trying to understand how a heavy object, like a bowling ball, behaves when it's dropped into a thick, hot soup. In the world of particle physics, this "soup" is called Quark-Gluon Plasma (QGP). It's the state of matter that existed just microseconds after the Big Bang and is recreated today in massive particle colliders like the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC).

Inside this soup, heavy particles called quarks (specifically "heavy quarks" like charm and bottom) try to stick together in pairs, forming "heavy quarkonia" (like a tiny, heavy atom). Scientists want to know: Does the hot soup melt these pairs apart? And if the soup gets denser (more crowded), does that change how they melt?

This paper by Jishnu Goswami and colleagues is a sophisticated attempt to answer that question using a technique called Lattice QCD (Quantum Chromodynamics). Here is the breakdown in everyday terms:

1. The Problem: The "Noisy Soup"

To study these particles, scientists use supercomputers to simulate the universe on a grid (a lattice). However, there's a catch.

• The Noise: The data coming from these simulations is incredibly "noisy," like trying to hear a whisper in a hurricane.
• The Density Issue: Most previous studies looked at the soup at a specific temperature but assumed it had zero "density" (no extra particles crowding it). But in real heavy-ion collisions (like those at RHIC), the soup is dense. We need to know how adding "crowd" (density) changes the rules.

2. The Method: The "Taylor Expansion" Recipe

Since the computer can't easily handle a super-dense soup directly, the scientists used a mathematical trick called a Taylor expansion.

• The Analogy: Imagine you are baking a cake. You know exactly how the cake tastes with 0 cups of sugar. You want to know how it changes if you add a little sugar. Instead of baking a whole new cake for every possible amount of sugar, you calculate the rate of change. You figure out how the taste changes for the first pinch of sugar, the second pinch, etc.
• In the Paper: They started with the "zero density" soup and calculated how the physics changes when you add a tiny bit of "chemical potential" (which is just a fancy word for density). They focused on the second-order effect (the "sugar squared" term), which tells them how the interaction changes as the soup gets slightly denser.

3. The Discovery: The Soup Gets "Thinner" and "Messier"

The team looked at the Static Potential, which is essentially the "glue" holding the heavy quark pair together. This glue has two parts:

1. The Real Part (The Glue): How strong the attraction is.
2. The Imaginary Part (The Friction): How much the soup bumps into the pair, causing them to wobble or break apart.

What they found:

• Enhanced Screening (The Glue Weakens): As the density of the soup increases, the "glue" holding the quarks together gets weaker, especially when the quarks are a bit far apart.
  • Analogy: Imagine two magnets trying to stick together. If you put a thick layer of foam between them, they can't feel each other as well. The denser the soup, the thicker the foam, and the easier it is for the magnets to drift apart.
• Enhanced Broadening (The Friction Increases): The "imaginary" part of the potential also grew. This means the dense soup is bumping into the quark pairs more aggressively, making them unstable and more likely to dissolve.

4. Why This Matters

This research is a crucial step for the Beam Energy Scan program at RHIC and future experiments at FAIR (Germany).

• The Goal: These experiments are trying to map out the "phase diagram" of the universe—basically, a map of how matter behaves under extreme heat and pressure.
• The Impact: By understanding how density changes the "glue" between heavy quarks, scientists can better interpret what happens in real collisions. It helps them figure out exactly when and how heavy particles melt in the early universe or in neutron stars.

Summary

Think of this paper as a detailed weather report for a microscopic storm. The scientists didn't just measure the wind (temperature); they figured out how the wind changes when the air gets crowded (density). They discovered that crowding makes the "glue" between heavy particles weaker and the environment more chaotic, which helps us understand the life and death of matter in the most extreme conditions in the universe.

They used a "mathematical recipe" to predict these changes without having to simulate the impossible complexity of a fully dense soup, providing the first solid lattice-based evidence that density makes it harder for heavy quarks to stay together.

Technical Summary

1. Problem Statement

The central goal of this research is to understand the behavior of heavy quark bound states (quarkonia, such as charmonium and bottomonium) within the Quark-Gluon Plasma (QGP) under conditions of finite baryon density.

• Context: Heavy quarkonia serve as probes for the QGP created in heavy-ion collisions (e.g., at RHIC and LHC). Their suppression is governed by color screening and collisional effects, described by a complex thermal potential ($V = V_{Re} + iV_{Im}$).
• Gap: While the thermal potential at zero chemical potential ($\mu_B = 0$) is well-studied, its behavior at nonzero density (relevant for the Beam Energy Scan at RHIC and future FAIR experiments) remains poorly constrained.
• Challenge: Extracting the thermal potential from Euclidean lattice QCD data is an ill-posed problem requiring analytic continuation. Extending this to finite density introduces the "sign problem," making direct simulations difficult.

2. Methodology

The authors employ Lattice QCD with a Taylor expansion approach to circumvent the sign problem at small chemical potentials.

• Simulation Setup:

  • Action: (2+1)-flavor Highly Improved Staggered Quark (HISQ) action.
  • Parameters: Temperature fixed at $T = 151.92$ MeV (near the pseudo-critical temperature $T_{pc} \approx 156.5$ MeV). Lattice size $64^3 \times 16$.
  • Quark Masses: Degenerate light quarks ($m_u=m_d$) and physical strange quark mass, with a ratio $m_l/m_s = 1/27$.
  • Statistics: Approximately 7,014 gauge configurations.

• Taylor Expansion Strategy:

  • Instead of simulating at finite $\mu$, the Wilson-line correlator $W(r, \tau)$ is expanded around vanishing chemical potentials ($\hat{\mu} = \mu/T = 0$) up to second order ($\hat{\mu}^2$).
  • Channels: The expansion is performed in two bases:
    1. Flavor Basis: Light ($u, d$) and Strange ($s$) quark chemical potentials.
    2. Conserved Charge Basis: Baryon number ($B$), Electric charge ($Q$), and Strangeness ($S$).
  • Curvature Coefficients: The second-order coefficients ($Y_2$) of the logarithm of the correlator are calculated. These are derived from second-derivative estimators involving traces of inverse Dirac matrices.
  • Noise Reduction: Stochastic estimation using Gaussian random noise vectors (2000 vectors for the noisiest term $(D_u^1)^2$, 500 for others).

• Potential Extraction & Renormalization:

  • Ansatz: A specific parametrization (Eq. 7) is used to model the $\tau$-dependence of the Wilson line, allowing the extraction of $V_{Re}$ and $V_{Im}$.
  • Renormalization: The bare potential contains a divergent self-energy term dependent on lattice spacing. The authors use gradient flow to regulate this and subtract the zero-temperature self-energy contribution ($\delta V_0$) to obtain a renormalized potential.

3. Key Contributions

1. First Lattice Constraints on Finite-Density Potential: This work provides the first extraction of the $\hat{\mu}^2$ corrections to both the real and imaginary parts of the static thermal potential in (2+1)-flavor QCD.
2. Dual Channel Analysis: The study simultaneously analyzes flavor-specific channels ($u, s$) and conserved charge channels ($B, Q$), offering a comprehensive view of how different density sources affect the potential.
3. Imaginary Part at Finite Density: A significant contribution is the determination of the density dependence of the imaginary part of the potential, which relates to the thermal width and dissociation rates of quarkonia.
4. Gradient Flow Subtraction: The implementation of a flow-based subtraction scheme to handle the additive self-energy divergence, ensuring the potential is properly renormalized for comparison across different flow times.

4. Results

• Enhanced Screening: The $\hat{\mu}^2$ coefficient (the curvature) shows a clear enhancement of in-medium screening effects at intermediate and large separations ($r \gtrsim 0.5$ fm).
  • The magnitude of the correction increases with separation distance.
  • The effect is visible in both the real part (screening) and the imaginary part (broadening).
• Flavor Dependence: The light-quark contribution ($Y_u$) is larger than the strange-quark contribution ($Y_s$), consistent with the lighter mass of up/down quarks.
• Real Part Behavior: The real part of the potential decreases (becomes more screened) as the chemical potential increases. This effect is small at short distances but becomes pronounced at larger separations ($r \gtrsim 0.5$ fm) and higher $\mu/T$.
• Imaginary Part Behavior: The $\hat{\mu}^2$ contribution to the imaginary part suggests that finite density enhances in-medium broadening, which is critical for understanding quarkonium dissociation.
• Flow Time Dependence: The results show weak dependence on the gradient flow time ($\tau_f/a^2 = 0.2$ vs $0.4$) at intermediate and large separations, though short-distance physics still requires extrapolation to $\tau_f \to 0$.

5. Significance and Outlook

• Phenomenological Impact: These results provide essential input for Schrödinger equation models used to predict quarkonium yields in low-energy heavy-ion collisions (RHIC Beam Energy Scan, FAIR). The enhanced screening and broadening at finite density imply that quarkonia may dissociate more readily in dense matter than previously estimated from zero-density studies.
• Theoretical Advancement: By moving beyond static free energies (which only give the real part) to the full complex potential, this work offers a more rigorous framework for describing the dynamical evolution of heavy quarks in the QGP.
• Future Directions: The authors plan to:
  • Increase statistics to reduce errors at large separations.
  • Perform a controlled extrapolation to the zero-flow-time limit.
  • Estimate the continuum limit.
  • Extend the analysis to a wider range of temperatures.

In summary, this paper establishes a foundational lattice QCD framework for quantifying how finite baryon density modifies the heavy-quark interaction, revealing that density significantly enhances both the screening and the thermal width of quarkonia in the QGP.