Diffusion of multiple conserved charges from entropy production

Using the Chapman-Enskog method within kinetic theory, this paper derives first- and second-order dissipative relativistic hydrodynamic equations for a multi-component quark-gluon plasma with baryon, electric, and strangeness charges, explicitly calculating the temperature and chemical potential dependence of the resulting diffusion matrix elements and their ratio to shear viscosity.

The Gist

Imagine a massive, chaotic party where thousands of guests are dancing, bumping into each other, and moving around a room. In the world of physics, this "party" is a Quark-Gluon Plasma (QGP)—a super-hot, super-dense soup of particles created when heavy atomic nuclei smash together at nearly the speed of light.

This paper is like a detailed instruction manual for predicting how this chaotic party evolves over time. Specifically, the authors are trying to figure out how different "types" of guests move and mix when the party isn't perfectly balanced (which is always the case in real life).

Here is a breakdown of their work using simple analogies:

1. The Three Types of "Guests" (Conserved Charges)

In this particle party, every guest carries three specific "ID tags" that cannot be lost or created out of thin air:

• Baryon Number (B): Think of this as a "Guest Count" tag. It keeps track of how many matter particles are there versus anti-matter particles.
• Electric Charge (Q): This is the "Positive/Negative" tag.
• Strangeness (S): This is a special "Exotic Flavor" tag carried only by certain particles (strange quarks).

In previous studies, scientists often only tracked the "Guest Count" (Baryon number). However, the authors of this paper realized that to truly understand the party, you have to track all three tags simultaneously because they influence each other.

2. The Problem: The "Traffic Jam" of Diffusion

When the party is out of balance (for example, if there are too many "Guests" in one corner of the room), they naturally try to spread out to even things out. This spreading process is called diffusion.

The authors discovered something tricky: The tags are connected.
Imagine you are trying to move a crowd of people who are holding red, blue, and green balloons. If you push the red balloons to the left, the blue and green balloons might accidentally get pushed to the right or left as well, depending on how the crowd is tangled.

• In physics terms, the movement of "Baryon Number" can cause "Electric Charge" to move, and vice versa.
• The paper calculates a "Diffusion Matrix." Think of this as a complex map or a traffic control chart that tells you exactly how much one type of charge will move when you try to move another type.

3. The Method: The "Relaxation Time" Guess

To solve the math of how these particles move, the authors used a method called the Chapman-Enskog expansion.

• The Analogy: Imagine trying to predict how a crowd moves after a sudden push. Instead of tracking every single person's footstep (which is impossible), you assume the crowd has a "relaxation time." This is like saying, "If the crowd is pushed, it will take this amount of time to settle back into a calm, organized flow."
• They used this "relaxation" idea to write down equations that describe how the "traffic" of charges flows, first in a simple, immediate way (like a car braking instantly) and then in a more complex, delayed way (like a car taking a moment to react before braking).

4. The Key Findings: The "Heat" of the Matter

The authors ran simulations to see how these diffusion rules change based on two main factors: Temperature (how hot the party is) and Chemical Potential (how crowded the room is with specific types of guests).

• The "Cross-Talk": They found that the "cross-diffusion" (how one charge drags another along) is significant. It's not just a straight line; the movement of one charge creates ripples that affect the others.

• The Competition: They found that diffusion is a tug-of-war between two forces:

  1. The Kinetic Term: How fast the particles are zipping around due to heat.
  2. The Thermodynamic Term: How the density and pressure of the crowd push back.
  • Result: At very high temperatures, the heat wins, and the particles move freely. But as the crowd gets denser (higher chemical potential), the "push back" from the crowd becomes so strong that the diffusion slows down significantly.

• Viscosity vs. Diffusion: They compared the "stickiness" of the fluid (viscosity) to the "spreading" ability (diffusion). They found that as the crowd gets denser, the fluid becomes "stickier" (viscosity dominates), making it harder for the charges to diffuse through the medium.

5. Why This Matters (According to the Paper)

The paper doesn't claim to cure diseases or build new engines. Instead, it provides the mathematical foundation for understanding the early moments of heavy-ion collisions (like those at the Large Hadron Collider).

By creating these detailed equations for how Baryon, Electric, and Strangeness charges move together, the authors provide a better "rulebook" for physicists to simulate what happens in these high-energy collisions. This is crucial for understanding the QCD Critical Point—a theoretical "phase transition" in the universe where matter changes state, which scientists are actively hunting for in experiments.

In summary: The authors built a sophisticated traffic model for a super-hot particle soup, showing that the movement of different particle "tags" is deeply interconnected, and that the density of the crowd plays a massive role in how fast or slow these tags can spread out.

Technical Summary

Technical Summary: Diffusion of Multiple Conserved Charges from Entropy Production

Problem Statement
Relativistic heavy-ion collision (HIC) experiments at facilities like RHIC and the LHC have successfully modeled the evolution of Quark-Gluon Plasma (QGP) in the high-temperature, low-baryon-chemical-potential regime. However, there is a critical need to understand the QCD phase structure at finite baryon density, particularly in the context of the search for the QCD critical endpoint. In this regime, the produced medium is baryon-rich, and the dynamics of conserved charges—specifically baryon number ($B$), electric charge ($Q$), and strangeness ($S$)—become paramount.

While previous hydrodynamic frameworks often treated net-baryon number as the sole conserved charge, a more complete description requires accounting for the thermalization of all light quark flavors ($u, d, s$). In such a multi-component system, the diffusion of one conserved charge is non-trivially coupled to the gradients of others. These cross-diffusion effects are encoded in the diffusion matrix elements ($\kappa_{qq'}$). A rigorous derivation of the dissipative hydrodynamic equations, including second-order terms necessary for causality and stability, and a quantitative estimation of these multi-component diffusion coefficients within a kinetic theory framework, remains an open challenge.

Methodology
The authors derive dissipative relativistic hydrodynamic equations for a mixture of massive fermions (quarks $u, d, s$) and massless bosons (gluons $g$) carrying multiple conserved charges ($B, Q, S$). The derivation proceeds within the kinetic theory approach using the following steps:

1. Relaxation Time Approximation (RTA): The relativistic Boltzmann equation is solved using the RTA with a momentum-independent relaxation time ($\tau_R$) for the collision term.
2. Chapman-Enskog (CE) Expansion: The out-of-equilibrium distribution function is expanded around the local equilibrium distribution in powers of space-time gradients. The authors derive both first-order (Navier-Stokes limit) and second-order evolution equations for the dissipative currents.
3. Entropy Production: Utilizing Boltzmann's H-theorem, the authors calculate the entropy production rate by taking the divergence of the entropy current. They demand that entropy production be non-negative for a dissipative system.
4. Transport Coefficient Extraction: By matching the entropy production expression to the standard form involving dissipative currents (bulk pressure $\Pi$, shear stress $\pi^{\mu\nu}$, and charge diffusion currents $n_q^\mu$), the authors identify the transport coefficients. This includes the shear viscosity ($\eta$), bulk viscosity ($\zeta$), and the full diffusion matrix elements ($\kappa_{qq'}$).
5. Thermodynamic Integrals: The coefficients are expressed in terms of thermodynamic integrals ($I_{nq}^{\pm}$ and $J_{nq}^{\pm}$) derived from the equilibrium distribution functions, which depend on temperature ($T$) and chemical potentials ($\mu_q$).

Key Contributions

• Derivation of Second-Order Equations: The paper provides a systematic derivation of second-order hydrodynamic evolution equations for a multi-component system with three conserved charges. These equations include relaxation terms and coupling between different dissipative currents, ensuring a causal and stable formulation.
• Diffusion Matrix Elements: Beyond the standard shear and bulk viscosities, the authors explicitly derive the diffusion matrix elements ($\kappa_{qq'}$) for the $B, Q, S$ system. They demonstrate that the diagonal components are manifestly positive and that the off-diagonal (cross-diffusion) terms arise naturally from the coupling of charge gradients.
• Kinetic vs. Thermodynamic Competition: The authors decompose the diffusion coefficients into a "kinetic term" (dependent on the distribution function and particle masses) and a "thermodynamic term" (dependent on macroscopic variables like energy density and pressure). They analyze how the competition between these terms dictates the behavior of the diffusion coefficients.
• Numerical Estimates: The paper provides numerical estimates for the temperature and chemical potential dependence of the diffusion matrix elements for a $(2+1)$ flavor QGP.

Results
The authors present numerical results for the scaled diffusion coefficients ($\kappa_{qq'}/(\tau_R T^3)$) and the ratio of charge conductivity to shear viscosity ($\kappa_{qq'}T/\eta$) across a range of temperatures ($150-500$ MeV) and chemical potentials:

• Diagonal Components ($\kappa_{BB}, \kappa_{QQ}, \kappa_{SS}$):
  • $\kappa_{BB}$ increases with temperature but is suppressed by finite baryon chemical potential ($\mu_B$) due to the negative contribution of the thermodynamic term.
  • $\kappa_{QQ}$ and $\kappa_{SS}$ are primarily controlled by the kinetic term. $\kappa_{QQ}$ shows a slow rise with temperature, while $\kappa_{SS}$ increases with $\mu_B$ because the kinetic term (driven by strange quark population) dominates the suppression from the enthalpy term in the denominator.
• Off-Diagonal Components ($\kappa_{QB}, \kappa_{BS}, \kappa_{QS}$):
  • Cross-diffusion coefficients exhibit complex behavior due to the sign differences in quark charges. For instance, $\kappa_{BS}$ is negative because the baryonic and strangeness charges of the strange quark have opposite signs.
  • The behavior of cross-coefficients changes significantly when varying $\mu_S$ versus $\mu_B$, reflecting the specific charge composition of the dominant quark species.
• Scaling with Shear Viscosity:
  • The ratio $\kappa_{qq'}T/\eta$ saturates to constant values at low $\mu_q/T$ but approaches zero at large $\mu_q/T$. This suppression at high chemical potential is attributed to the rapid growth of enthalpy ($\varepsilon + P$) in the denominator and Pauli blocking effects.
  • The flavor dependence ($N_f = 1, 2, 2+1$) shows that adding flavors generally increases $\kappa_{BB}T/\eta$ (due to additive baryonic charges) but suppresses $\kappa_{QQ}T/\eta$ (due to opposing electric charges).
• Thermal Conductivity: The ratio of thermal conductivity to shear viscosity ($R_T$) is found to match AdS/CFT predictions ($C_f$) in the limits of small and large $\mu_B/T$, with deviations in the intermediate range attributed to mass effects and multiple conserved charges.

Significance and Claims
The paper claims to provide a necessary theoretical framework for modeling multi-component diffusion dynamics in the initial state of heavy-ion collisions, particularly those relevant to the Beam Energy Scan program where baryon stopping is significant. By deriving the second-order evolution equations and estimating the diffusion matrix elements, the work addresses the need for stable and causal hydrodynamic simulations of the QGP at finite baryon density.

The authors emphasize that the presence of non-trivial cross-diffusion coefficients ($\kappa_{qq'}$ where $q \neq q'$) introduces novelty to the hydrodynamic framework, as the diffusion current of a specific charge is no longer solely determined by its own gradient. This is crucial for understanding correlations and fluctuations of conserved charges near the QCD critical point. The study serves as a baseline for future phenomenological investigations, though the authors explicitly note that their current work is limited to systems with momentum-independent relaxation times and does not yet include hadronic medium effects or magnetic field influences, which are identified as directions for future research.