Critical slowing down and bulk viscosity in binary neutron star mergers

This paper demonstrates that if neutron star merger matter passes near a QCD critical point, critical slowing down can significantly enhance bulk viscosity to levels comparable to electroweak contributions, potentially leaving observable imprints on the merger's hydrodynamic evolution.

The Gist

Imagine two giant, super-dense stars (neutron stars) dancing a violent waltz, spiraling closer and closer until they crash into each other. This cosmic collision is one of the most energetic events in the universe. To predict what happens next—how they merge, what sound waves (gravitational waves) they emit, and what light they flash—scientists use a set of rules called hydrodynamics.

Think of hydrodynamics like the rules for how water flows. Usually, when we model water, we assume it flows smoothly. But inside these crashing stars, the "fluid" is made of subatomic particles (quarks and gluons) behaving in extreme ways.

This paper asks a fascinating question: What if the fluid inside these stars hits a "traffic jam" caused by a hidden phase transition, making it stickier than we thought?

Here is the breakdown of the paper's ideas using everyday analogies:

1. The Usual Rules: The "Slow Clock"

Normally, scientists assume there is a clear hierarchy of time in these stars:

• Strong interactions (particles bumping into each other) happen super fast.
• Weak interactions (particles changing flavor, like a neutron turning into a proton) happen slower.
• Hydrodynamics (the big flow of the star) happens the slowest.

Because the "weak" interactions are the slowest, they act like the bottleneck. They determine how much bulk viscosity (internal friction) the star has. Think of viscosity like the thickness of honey. If the honey is thick, it resists moving. In neutron stars, this "thickness" is usually determined by how fast particles can change their identity.

2. The Twist: The "Critical Point" Traffic Jam

The authors propose a scenario where the matter in the star passes near a Critical Point (CP).

• The Analogy: Imagine a highway that is usually clear. Suddenly, you hit a specific spot where the road conditions change drastically. Instead of cars just slowing down, the traffic starts to "stall" completely. The cars (particles) start to fluctuate wildly, trying to decide whether to be in one lane or another, but they can't make up their minds.
• Critical Slowing Down: Near this Critical Point, the particles get "confused" and take a very long time to settle down. This is called Critical Slowing Down.
• The Result: Because the particles are stuck in this confused state, the "friction" (viscosity) of the star skyrockets. It's like the honey suddenly turning into solid glue.

3. The Limits: Why isn't it infinite?

You might ask, "If the friction gets so high, does the star stop moving?"
The paper explains that nature puts a few "brakes" on this effect:

• Time Limit: The star merger happens very fast (milliseconds). The particles don't have infinite time to get confused; they only get confused for a short while.
• Size Limit: The "confused zone" (correlation length) can't get bigger than the star itself.
• The "Smoothness" Limit: Hydrodynamics assumes the fluid is smooth. If the "confused zone" gets too big, the fluid looks bumpy and lumpy, and the smooth-flow rules break down.

The authors calculate that even with these limits, the "confused zone" can grow to about 10 to 100 nanometers. That's tiny (smaller than a virus), but in the world of subatomic physics, it's huge!

4. The Big Reveal: Glue vs. Honey

The team did the math to see how much this "critical friction" adds to the star's viscosity.

• The Old View: The friction was mostly due to weak nuclear forces (the "honey").
• The New View: If the star passes near this Critical Point, the "critical friction" can become thousands or even millions of times stronger than the usual friction.
• The Catch: This only happens if the star's path takes it right through the "confused zone." However, the authors show that this zone could be large enough (hundreds of meters) to matter in computer simulations of the merger.

5. Why Should We Care?

If this "super-sticky" effect happens, it changes the story of the merger:

• Gravitational Waves: The way the stars crash and ring like a bell would sound different. The "glue" would dampen the vibrations, changing the signal we detect on Earth.
• Finding the Critical Point: Neutron star mergers might be the only place in the universe hot and dense enough to find this elusive "Critical Point" in the laws of physics. By listening to the "sound" of the crash, we might be able to map the phase diagram of the universe.

Summary

This paper suggests that inside colliding neutron stars, there might be a hidden "traffic jam" of particles. If the stars pass through this jam, the internal friction of the star explodes, potentially overpowering all other forces. This could leave a unique fingerprint on the gravitational waves we detect, offering us a new way to test the fundamental laws of matter under extreme pressure.

In short: The universe might be using neutron star collisions to show us a "glitch" in the matrix of physics, where matter gets temporarily stuck, creating a super-viscous fluid that changes the sound of the crash.

Technical Summary

1. Problem Statement

Binary neutron star (BNS) mergers are typically modeled using relativistic hydrodynamics, which relies on a strict hierarchy of timescales: $\tau_{\text{strong}} \ll \tau_{\text{weak}} \ll \tau_{\text{hydro}}$. Under this assumption, the system is treated as an ideal fluid where dissipative effects (viscosity) are negligible or dominated by weak interaction processes (specifically Urca processes).

However, this hierarchy is not guaranteed. If the thermodynamic trajectory of the merging matter passes near a QCD Critical Point (CP), Critical Slowing Down (CSD) occurs. This phenomenon introduces a new, slow relaxation timescale ($\tau_{\text{CSD}}$) associated with the order parameter fluctuations. If $\tau_{\text{CSD}} \sim \tau_{\text{hydro}}$, critical fluctuations can dominate the bulk viscosity ($\zeta$), potentially altering the hydrodynamic evolution and gravitational wave emission of the merger. The paper investigates whether such a CP, if it exists within the phase space of a BNS merger, can generate a bulk viscosity large enough to rival or exceed standard electroweak contributions.

2. Methodology

The authors employ a multi-step theoretical framework combining critical phenomena theory, hydrodynamics, and nuclear equation of state (EoS) modeling:

• Universality Class Mapping: They assume the CP belongs to the static universality class of the 3D Ising model and the dynamic universality class of Model H. They map QCD variables (baryon chemical potential $\mu_B$, isospin chemical potential $\mu_I$, temperature $T$) to Ising variables (reduced temperature $r$, magnetic field $h$) using a linear transformation.
• Bounding the Correlation Length ($\xi$): In equilibrium, $\xi$ diverges at the CP. However, in a dynamic merger, $\xi$ is limited by:
  1. Finite Time: The time the system spends near the CP limits relaxation. Using $\tau \sim \xi^z$, they estimate a maximum $\xi$ based on the merger duration ($\sim 1$ ms).
  2. Finite Size/Gradient: $\xi$ cannot exceed the scale over which thermodynamic fields (like temperature) vary significantly ($l \sim \Delta T / |\nabla T|$).
  3. Hydrodynamic Validity: If $\xi$ becomes comparable to the coarse-graining scale of the simulation, hydrodynamics breaks down.
• Bulk Viscosity Calculation: They calculate the critical contribution to bulk viscosity ($\zeta_{\text{CSD}}$) using the Kawasaki approximation for Model H dynamics. The formula involves the correlation length, susceptibility, and a driving frequency ($\omega \sim 1$ kHz).
• Equation of State (EoS): To model the non-critical background, they use an Excluded Volume Hadron Resonance Gas (EV-HRG) model. This accounts for repulsive interactions between hadrons and includes the PDG21+ hadron list.
• Parameter Scanning: They scan across various critical temperatures ($T_c$) and chemical potentials ($\mu_{B,c}$), varying the non-universal mapping parameter $\rho w$ to account for uncertainties in the microscopic details of the CP.

3. Key Contributions

• Quantification of CSD Limits: The paper provides rigorous order-of-magnitude estimates for the maximum achievable correlation length in a BNS merger, constrained by finite-time and finite-size effects. They find $\xi_{\text{max}}$ is limited to the nanometer scale ($\sim 10$ nm to $100$ nm), rather than diverging to infinity.
• Hydrodynamic Validity Check: A crucial contribution is determining the regime where hydrodynamics remains valid despite large critical fluctuations. They demonstrate that even with significant viscosity enhancements, hydrodynamics holds for coarse-graining scales typical of merger simulations (tens to hundreds of meters), provided the non-universal parameter $\rho w$ is not excessively large.
• Comparison with Electroweak Baseline: The study explicitly compares the critical bulk viscosity against the standard electroweak bulk viscosity ($\zeta_{\text{EW}}$) in both neutrino-trapped and neutrino-transparent regimes.

4. Results

• Maximum Correlation Length:
  • Finite-time constraints limit $\xi$ to $\sim 10$ nm.
  • Finite-temperature gradient constraints allow $\xi$ up to $\sim 100$ nm.
  • The effective limit is the minimum of these two, generally $\sim 10$ nm.
• Bulk Viscosity Enhancement:
  • The critical bulk viscosity $\zeta_{\text{CSD}}$ scales strongly with temperature and correlation length.
  • For a broad range of parameters (specifically when the non-universal mapping parameter $\rho w \gtrsim 100$), the critical contribution can reach $10^{25} - 10^{33}$ g/cm·s.
  • This value significantly exceeds the electroweak baseline ($\zeta_{\text{EW}}$), particularly in the neutrino-trapped regime ($T_c > 10$ MeV).
• Spatial Extent:
  • The region where $\zeta_{\text{CSD}} \gtrsim \zeta_{\text{EW}}$ is not infinitesimal. For $\rho w \gtrsim 1$, this enhancement occurs over macroscopic scales comparable to or larger than typical fluid cells in BNS simulations.
  • Hydrodynamics breaks down only if the coarse-graining scale is reduced to $\ll 10^{-4}$ m (for $\rho w \sim 100$) or $\ll 10$ m (for $\rho w \sim 10^6$), which is far smaller than the resolution of current merger simulations.
• Neutrino-Transparent Regime: If the CP exists at very low temperatures ($T_c \sim 2$ MeV, neutrino-transparent), the enhancement is less distinct because the background electroweak viscosity is higher, but a meaningful enhancement still exists for $\rho w \gtrsim 1$.

5. Significance

• Observational Imprints: The findings suggest that if a QCD CP exists within the phase space of a BNS merger, it could leave observable imprints on the hydrodynamic evolution. This could manifest as modified gravitational wave signals or changes in the post-merger remnant's lifetime and stability due to enhanced bulk viscous dissipation.
• New Search Strategy: The paper proposes a novel avenue for searching for the QCD Critical Point using astrophysical observations of compact objects, complementing terrestrial heavy-ion collision experiments (like those at RHIC or LHC).
• Theoretical Robustness: By demonstrating that critical effects can dominate viscosity without invalidating the hydrodynamic framework (for realistic simulation resolutions), the authors establish that critical slowing down is a viable and potentially dominant physical mechanism in neutron star mergers, necessitating its inclusion in future high-precision simulations.

In conclusion, the paper establishes that critical slowing down near a QCD critical point can generate bulk viscosity orders of magnitude larger than standard weak-interaction processes, potentially dominating the dissipative dynamics of binary neutron star mergers and offering a new probe for the QCD phase diagram.