Fermi-liquid view of viscosity in cold and dense nucleon matter

This paper develops a Fermi-liquid framework to derive leading-order expressions for shear and bulk viscosities in cold, dense nucleon matter, demonstrating that the bulk-to-shear viscosity ratio scales as $(T/\mu^*)^4$ in the degenerate regime and remains robust against quasiparticle mass corrections when coupled with a Walecka-type equation of state.

The Gist

Imagine a crowded dance floor where everyone is moving in perfect, synchronized rhythm. This is what physicists call "cold, dense nucleon matter"—a state of matter found inside neutron stars or created briefly in particle accelerators, where protons and neutrons are packed tightly together and moving very slowly relative to their energy.

In this paper, the authors act like engineers trying to understand how this "dance floor" resists being pushed, squished, or twisted. They are calculating two specific types of "stickiness" or resistance, known as viscosity:

1. Shear Viscosity (The "Twist" Resistance): Imagine trying to slide one layer of the dance floor past another, like shuffling a deck of cards. The resistance you feel is shear viscosity.
2. Bulk Viscosity (The "Squeeze" Resistance): Imagine trying to compress the entire dance floor into a smaller ball or expand it like a balloon. The resistance to this volume change is bulk viscosity.

The Problem They Solved

In previous studies, scientists had a tool (a mathematical framework based on "Fermi-liquid theory") to calculate these resistances, but it had a glitch. When they tried to calculate the "squeeze" resistance (bulk viscosity), the math sometimes gave a negative number.

In the real world, resistance cannot be negative (you can't have a fluid that helps you squeeze it while you're trying to squeeze it; that would violate the laws of physics). The authors realized this happened because they hadn't set the "rules of the game" correctly for how the particles interact with their environment.

The Fix: They introduced a set of "Landau matching conditions." Think of this as calibrating a scale. Before weighing an object, you must ensure the scale reads zero when empty. Similarly, the authors ensured that their mathematical model correctly accounted for the fact that the particles' mass and energy change depending on how crowded the room is. Once they fixed this calibration, they proved mathematically that the "squeeze" resistance is always positive (or zero), fixing the glitch.

The Big Discovery: The "Silent" Squeeze

Once the math was fixed, they looked at what happens when the temperature is extremely low (which is the case for the dense matter they are studying).

They found a massive difference between the two types of resistance:

• Shear Viscosity (Twisting): Even at very low temperatures, the fluid still resists being twisted. It's like trying to stir honey; it's thick and slow.
• Bulk Viscosity (Squeezing): This resistance essentially disappears. It becomes so tiny it's almost zero.

The Analogy:
Imagine the dance floor is made of perfectly round, hard marbles packed tight.

• If you try to twist the floor (Shear), the marbles have to roll over each other. Because they are packed so tightly, they can't move easily, creating a lot of friction (high viscosity).
• If you try to squeeze the floor (Bulk), the marbles just shift slightly to fit the new shape. Because they are already in a perfect, efficient arrangement (the "Fermi surface"), they can rearrange themselves without losing any energy. It's like a perfectly organized bookshelf; you can slide the books slightly to make room, but you don't need to force them or generate heat.

The authors found that as the system gets colder, the "squeeze" resistance drops off incredibly fast—much faster than the "twist" resistance. In fact, the ratio of squeeze-resistance to twist-resistance shrinks by the fourth power of the temperature. This means in the cold, dense world of neutron stars or heavy-ion collisions, squeezing the matter is almost frictionless, but twisting it is very difficult.

Why Does This Matter?

The authors applied their new, corrected math to a specific model of nuclear matter (the Walecka model) to predict how real neutron matter behaves.

They conclude that for experiments trying to study this dense matter (like those at the Electron-Ion Collider or in heavy-ion collisions), scientists should focus on the "twist" (shear) effects. The "squeeze" (bulk) effects are so small in this cold, dense regime that they are likely too weak to be noticed or to affect the outcome of the experiment.

In short: The authors built a better ruler to measure how "sticky" dense nuclear matter is. They proved that while this matter is very hard to twist, it is almost perfectly easy to squeeze when it is cold and dense, correcting a previous mathematical error that made the "squeeze" resistance look weird or impossible.

Technical Summary

Technical Summary: Fermi-liquid view of viscosity in cold and dense nucleon matter

Problem Statement
The transport properties of cold, dense nucleonic matter are critical for understanding non-equilibrium dynamics in intermediate-energy heavy-ion collisions and the cooling mechanisms of neutron stars. While transport theory for hot QCD matter is well-developed, a controlled description for cold, dense systems near the Fermi surface remains challenging. A specific theoretical gap exists in the quasiparticle picture of Fermi liquids: without rigorous constraints, the sign of the bulk viscosity ($\zeta$) is not automatically guaranteed to be positive, and the relationship between transport theory and hydrodynamics requires precise matching conditions. Furthermore, the impact of medium-dependent quasiparticle masses on the ratio of bulk to shear viscosity ($\zeta/\eta$) in the degenerate regime needed clarification.

Methodology
The authors develop a kinetic framework based on the Fermi-liquid theory at the mean-field level. The core methodology involves:

1. Linearized Relativistic Boltzmann Equation: The study starts with the linearized relativistic Boltzmann equation tailored for quasiparticles possessing a medium-dependent dispersion relation ($E^*_p = \sqrt{m^{*2} + p^2}$) and effective chemical potential ($\mu^*$).
2. Relaxation-Time Approximation (RTA): The collision kernel is treated within the RTA, where the collision term is simplified to $-\delta f / \tau_{rel}$.
3. Landau Matching Conditions: To resolve the non-uniqueness of the solution caused by the zero modes of the collision operator (associated with particle number and energy-momentum conservation), the authors implement Landau matching conditions. These conditions enforce that deviations in local thermodynamic quantities do not alter the local energy-momentum density or conserved charge density.
4. Low-Temperature Expansion: A low-temperature expansion (Sommerfeld expansion) is applied to derive leading-order expressions for transport coefficients in the degenerate regime ($T \ll \mu^*$).
5. Equation of State (EoS): The kinetic framework is coupled to a Walecka-type mean-field EoS, incorporating scalar ($\sigma$) and vector ($\omega$) meson exchanges to model nucleon interactions.

Key Contributions

• Proof of Non-Negative Bulk Viscosity: The paper establishes a self-consistent scheme proving that the bulk viscosity $\zeta$ is manifestly non-negative. This is achieved by deriving the general solution for the non-equilibrium distribution function and applying Landau matching conditions to fix the coefficients of the zero modes.
• Robustness of the $\zeta/\eta$ Scaling: The authors demonstrate that the leading-order behavior $\zeta/\eta \propto (T/\mu^*)^4$ in the degenerate regime is robust against corrections arising from the medium-dependent quasiparticle mass ($m^*$).
• Resolution of Zero-Mode Ambiguity: The work clarifies that the arbitrariness in the solution of the Boltzmann equation is associated solely with the bulk viscosity term, while the shear viscosity remains uniquely determined.
• Numerical Implementation: The framework is applied to cold, dense nucleon matter using the Walecka model, providing numerical results for shear and bulk viscosities across relevant densities.

Results

• Viscosity Expressions: The leading-order (LO) expressions derived are $\eta_{LO} \propto p_F^{*5} \tau_{rel} / \mu^*$ and $\zeta_{LO} \propto m^{*4} p_F^* T^4 / (\mu^{*5} \tau_{rel})$.
• Suppression of Bulk Dissipation: In the degenerate regime, bulk viscosity is parametrically negligible compared to shear viscosity ($\zeta \ll \eta$). As $T \to 0$, $\eta$ diverges (due to $\tau_{rel} \propto T^{-2}$) while $\zeta$ vanishes.
• Physical Interpretation: The vanishing of bulk viscosity at zero temperature is attributed to the fact that isotropic compression or expansion merely rescales the Fermi sphere without generating entropy or energy dissipation, as the system remains on the equilibrium manifold. Shear flow, by contrast, distorts the Fermi sphere, requiring momentum-randomizing collisions that are suppressed at $T=0$.
• Numerical Findings: Calculations using the Walecka model at $T=8$ MeV confirm that dissipation in the bulk channel is negligible compared to the shear channel. The study also notes that apparent pathologies (such as sign changes or divergences in the response coefficient $\kappa^*$) observed in mean-field solutions are artifacts of the saddle-point nature of the effective potential and do not indicate genuine instabilities in the properly matched hydrodynamics.

Significance
The paper provides a dedicated, theoretically consistent framework for nucleonic transport at low temperatures, ensuring that transport coefficients are fundamentally consistent with the symmetries, thermodynamics, and microscopic interactions of nuclear matter. By proving the non-negativity of bulk viscosity and establishing the robustness of the $\zeta/\eta$ scaling, the work addresses a critical gap in the understanding of quasiparticle systems. The results are directly relevant for interpreting signals in intermediate-energy heavy-ion experiments (where density gradients and compressional flows are sensitive to transport coefficients) and for modeling the hydrodynamic behavior of cold, dense matter in neutron stars and merger events. The authors suggest that while the current mean-field approach is a significant step, future work should incorporate microscopic collision integrals, mesonic fluctuations beyond the mean field, and isospin asymmetry to further refine these transport inputs.