Transport coefficients of chiral fluid dynamics using low-energy effective models

This paper calculates the first-order transport coefficients, specifically bulk and shear viscosities, for a chiral fluid of quasiparticles with temperature-dependent masses by applying a Chapman-Enskog expansion within a relativistic Boltzmann framework using the relaxation time approximation and thermal masses derived from the linear sigma and NJL models.

The Gist

Imagine you are trying to understand how a pot of water boils, but instead of water, it's a super-hot, super-dense soup of subatomic particles called quarks and gluons. This soup, known as the Quark-Gluon Plasma (QGP), is what existed just microseconds after the Big Bang and is recreated today in massive particle colliders like the Large Hadron Collider.

Physicists want to know how this soup flows. Does it flow like honey (thick and sticky)? Or like water (thin and slippery)? To answer this, they need to calculate "transport coefficients"—mathematical numbers that describe how easily the fluid moves, how it resists being squished, and how it conducts heat.

This paper is like a new recipe book for calculating these numbers, but with a special twist: it accounts for how the "ingredients" (the particles) change their weight as the soup gets hotter.

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

1. The Problem: The Particles Change Their Weight

In most physics models, particles are treated like billiard balls with a fixed weight. But in this hot soup, the particles are quasiparticles. Think of them like swimmers in a pool.

• When the water is cold, the swimmers wear heavy winter coats (they have a large "mass").
• As the water heats up, the coats melt away, and the swimmers become light and fast (their "mass" drops).

The authors used two different "rulebooks" (called the Linear Sigma Model and the NJL Model) to predict exactly how fast these coats melt as the temperature rises. They found that one rulebook predicts the coats melt very suddenly, while the other predicts a slower, smoother melt.

2. The Method: A Traffic Jam Analogy

To figure out how the fluid flows, the authors used a method called Kinetic Theory. Imagine a crowded highway:

• The Boltzmann Equation: This is the rulebook that tracks every single car (particle) on the highway, where they are going, and how they bump into each other.
• The Collision Term: This describes the traffic jams. When cars bump, they slow down or change direction.
• The Relaxation Time: This is the "cooling off" period. It's the time it takes for a car that just got bumped to get back to its normal speed and lane.

The Innovation:
Previous studies used a shortcut (the Anderson-Witting approximation) to calculate these traffic jams. The authors say this shortcut was like assuming all cars take the exact same amount of time to recover, regardless of how fast they were going. This broke the laws of physics (specifically, the conservation of energy).

In this paper, the authors used a new, smarter shortcut. They made sure their math respected the laws of physics even when the "recovery time" changed depending on how fast the particle was moving. It's like a traffic system that knows a speeding car takes longer to stabilize than a slow one.

3. The Results: What Happens When the Coats Melt?

The authors calculated the "stickiness" (viscosity) of the fluid using their new method and the two different rulebooks. Here is what they found:

• Shear Viscosity (The "Slipperiness"): This measures how easily the fluid slides past itself.

  • The Finding: As the temperature rises and the particles lose their "coats" (mass), the fluid becomes incredibly slippery. Once the particles are massless, the slipperiness levels off and stays constant. Both rulebooks agreed on this.

• Bulk Viscosity (The "Sponginess"): This measures how much the fluid resists being squeezed or expanded.

  • The Finding: This is where the magic happens. As the particles shed their coats (a process called Chiral Symmetry Restoration), the fluid stops resisting being squeezed. The "sponginess" drops to almost zero.
  • The Difference: The model where the coats melt suddenly (Linear Sigma Model) showed a sharp, dramatic drop in sponginess right before the critical temperature. The model where the coats melt slowly (NJL Model) showed a gentler, smoother decline.

• The Speed of Sound:

  • The speed at which a "ripple" travels through the fluid also changed. In the model with the sudden coat-melting, the speed of sound dipped sharply right before the transition, whereas the other model was smoother.

4. Why Does This Matter?

Think of the Quark-Gluon Plasma as a new state of matter that behaves like a "perfect fluid." By understanding exactly how the particles change their mass and how that affects the fluid's flow, scientists can better interpret the data from particle colliders.

If the fluid flows exactly as these new calculations predict, it tells us that the "melting coat" theory is correct. If the data doesn't match, we might need to look for other forces at play (like the authors suggest, perhaps involving "Polyakov loops," which are like invisible strings tying the particles together).

Summary

In short, this paper built a more accurate traffic simulator for the hottest soup in the universe. They showed that as the universe cools down (or heats up in a collider), the particles change their "weight," and this change dramatically alters how the fluid flows. They proved that the way this weight changes (suddenly vs. slowly) leaves a distinct fingerprint on the fluid's behavior, specifically on how "spongy" it is when you try to squeeze it.

Technical Summary

1. Problem Statement

The paper addresses the challenge of describing the transport properties (specifically viscosity) of the Quark-Gluon Plasma (QGP) and hot dense matter produced in heavy-ion collisions. While relativistic hydrodynamics has been successful in modeling the QGP, a rigorous microscopic derivation of its transport coefficients from quantum field theory remains difficult, particularly near the chiral phase transition.

Key issues addressed include:

• Incorporating Chiral Symmetry: Standard kinetic theories often neglect phase transitions. The authors aim to incorporate the effects of chiral symmetry breaking and restoration, which significantly alter the effective mass of quasiparticles (quarks) in the medium.
• Consistency of Approximations: Previous studies often used the Anderson-Witting Relaxation Time Approximation (RTA) to simplify collision terms. However, this standard RTA violates energy-momentum conservation when the relaxation time depends on particle energy or when specific matching conditions are used.
• Model Dependence: There is a need to understand how different low-energy effective models (specifically the Linear Sigma Model with constituent quarks and the Nambu-Jona-Lasinio model) influence transport coefficients via their temperature-dependent thermal masses.

2. Methodology

The authors employ a multi-step theoretical framework combining finite-temperature field theory with relativistic kinetic theory:

• Effective Models for Thermal Mass:

  • They utilize two chiral effective models: the Linear Sigma Model coupled with constituent quarks (LSMq) and the Nambu-Jona-Lasinio (NJL) model.
  • These models are solved at the mean-field level to extract the temperature-dependent effective quark mass, $M(T)$, which serves as the input for the kinetic theory.
  • The LSMq and NJL models predict different behaviors for chiral restoration: LSMq exhibits a faster restoration (sharper transition) compared to the smoother transition in NJL.

• Relativistic Boltzmann Equation with Background Field:

  • The system is described by a relativistic Boltzmann equation where the particle mass $M(T)$ varies with temperature.
  • To preserve energy-momentum conservation in the presence of a temperature-dependent mass, a temperature-dependent background field $B(T)$ is introduced into the energy-momentum tensor ($T^{\mu\nu} = \langle p^\mu p^\nu \rangle + g^{\mu\nu}B$).
  • A specific matching condition is chosen to define the local temperature and four-velocity, ensuring the particle number density matches the equilibrium distribution.

• Improved Relaxation Time Approximation (RTA):

  • The collision term is simplified using a novel RTA (based on Ref. [27]) that is consistent with fundamental conservation laws (energy and momentum).
  • Unlike the standard Anderson-Witting approximation, this improved RTA projects out the zero-eigenvalue modes of the collision operator, ensuring conservation even when the relaxation time $\tau_R$ has an energy dependence ($\tau_R \propto (E_p/T)^\gamma$).
  • The parameter $\gamma$ is treated as a phenomenological parameter representing the momentum dependence of the interaction.

• Chapman-Enskog Expansion:

  • The authors perform a first-order Chapman-Enskog expansion of the Boltzmann equation to derive the non-equilibrium distribution function.
  • This allows for the analytical calculation of first-order transport coefficients: shear viscosity ($\eta$), bulk viscosity ($\zeta$), and the energy deviation coefficient ($\chi$).

3. Key Contributions

• Consistent Derivation: The paper provides a self-consistent framework for calculating transport coefficients in chiral fluid dynamics without violating energy-momentum conservation, a common pitfall in previous RTA applications.
• Background Field Formalism: It explicitly demonstrates how to incorporate the temperature-dependent background field $B(T)$ required to maintain conservation laws when using quasiparticles with thermal masses derived from chiral models.
• Model Comparison: It offers a comparative analysis of transport coefficients derived from two distinct chiral models (LSMq and NJL), highlighting how the speed of chiral restoration affects physical observables.
• Analytical Expressions: The authors derive closed-form expressions for transport coefficients in terms of the thermal mass $M(T)$ and the relaxation time parameter $\gamma$.

4. Key Results

• Shear Viscosity ($\eta$):
  • Both models yield similar results for shear viscosity.
  • Above the critical temperature ($T_c$), as chiral symmetry is restored and $M(T) \to 0$, the shear viscosity approaches a constant value determined by $\gamma$. It becomes mass-independent in the high-temperature limit.
• Bulk Viscosity ($\zeta$) and Energy Deviation ($\chi$):
  • These coefficients are highly sensitive to the chiral phase transition.
  • As chiral symmetry is restored ($M(T) \to 0$), both $\zeta$ and $\chi$ drop sharply toward zero.
  • LSMq vs. NJL: Because the LSMq model exhibits a faster chiral restoration (narrower temperature range), the drop in bulk viscosity is more abrupt and pronounced compared to the smoother decrease seen in the NJL model.
  • Just before $T_c$, the bulk viscosity exhibits an abrupt increase, a feature linked to the rapid change in the thermal mass.
• Speed of Sound ($c_s^2$):
  • The speed of sound decreases sharply before $T_c$ in the LSMq model, whereas it varies more smoothly in the NJL model. This correlates with the rate of chiral restoration.
• Ratio $\zeta/\eta$:
  • The ratio of bulk to shear viscosity is small at high temperatures but shows significant structure near the phase transition, reflecting the dynamics of the chiral condensate.
• Dependence on $\gamma$:
  • The magnitude of the transport coefficients depends on the phenomenological parameter $\gamma$. Higher values of $\gamma$ generally lead to smaller bulk-to-shear ratios at high temperatures.

5. Significance

• Phenomenological Relevance: The results provide crucial inputs for hydrodynamic simulations of heavy-ion collisions. Understanding how bulk viscosity behaves near the chiral critical point is essential for interpreting experimental data (e.g., flow harmonics) from RHIC and the LHC.
• Theoretical Robustness: By resolving the conservation law violations in the RTA, this work establishes a more reliable baseline for future studies of transport in strongly coupled systems.
• Insight into Chiral Dynamics: The study demonstrates that the dynamics of the phase transition (how fast the mass drops) are as important as the transition itself in determining transport properties. The difference between LSMq and NJL results suggests that the specific nature of the chiral transition (crossover vs. first-order, width of the transition) leaves a distinct imprint on the fluid's viscosity.
• Future Directions: The authors propose extending this work to second-order transport coefficients and including the dynamical evolution of the chiral field (non-equilibrium corrections to the mass), which would further refine the description of the QGP.