Burgers equation for the bulk viscous pressure of quark matter

This paper derives a Burgers-type evolution equation for the bulk viscous pressure in unpaired quark matter, identifying four key transport coefficients based on electroweak decay rates and evaluating them across different equations of state to provide a new framework for simulating dissipation in compact star mergers.

The Gist

Imagine a cosmic dance floor where two massive stars, made of incredibly dense matter, crash into each other. This is a neutron star merger, an event so violent it ripples through the universe as gravitational waves.

Now, imagine that inside these stars, the matter isn't just solid or liquid like water or iron. Under such extreme pressure, the atoms break apart, and their building blocks (quarks) float freely in a "soup" called quark matter.

This paper is about understanding how this "quark soup" behaves when it gets squeezed and stretched during the crash. Specifically, the authors are looking at bulk viscosity.

The Analogy: The Sticky Honey vs. The Bouncy Rubber Band

To understand viscosity, think of two fluids:

1. Water: If you stir it, it moves easily and stops quickly. It has low viscosity.
2. Honey: If you stir it, it resists. It's thick and "sticky." It takes time to settle. This is viscosity.

In the context of these stars, bulk viscosity is the resistance the matter offers when it is rapidly compressed and expanded (like a bell being rung). If the matter is very "sticky" (high viscosity), it acts like a shock absorber, damping out the vibrations and stopping the star from ringing like a bell. This affects the gravitational waves we detect on Earth.

The Problem: Old Maps Don't Fit New Territory

For a long time, scientists used a standard rulebook (called the Israel-Stewart equation) to describe how this "sticky" matter behaves. Think of this rulebook as a simple map for a flat road. It works okay for normal situations.

However, the authors discovered that unpaired quark matter (the specific type of soup inside these stars) is more complex. It's not just a flat road; it's a terrain with hills, valleys, and two different types of traffic.

The old map (Israel-Stewart) assumes the matter reacts in one simple way. But the authors found that quark matter actually behaves like a two-component fluid. It has two different "modes" of reacting to being squeezed, each with its own speed and stickiness.

The Solution: The "Burgers" Equation

The authors derived a new mathematical rule, which they call the Burgers equation.

• The Old Way (Israel-Stewart): Imagine a single spring that stretches and snaps back. It has one speed.
• The New Way (Burgers): Imagine a system with two springs connected in a tricky way. One spring is stiff and snaps back fast; the other is loose and takes a long time to settle.

The paper shows that to accurately predict how the star vibrates and how the gravitational waves look, you must use this two-spring (Burgers) model, not the single-spring model.

The Four Key Ingredients

To make this new map work, the authors calculated four specific numbers (transport coefficients) that act like the settings on a thermostat:

1. Two "Relaxation Times": How long does it take for the matter to stop vibrating? One is fast, one is slow.
2. Two "Viscosity Coefficients": How "thick" or "sticky" is the matter for each of those two speeds?

They calculated these numbers for two different scenarios:

• The "Bag Model": Like a heavy, dense bag of marbles (good for lower densities).
• Perturbative QCD: A high-tech, theoretical model for extremely high densities (like the very center of the star).

The "Switch" in Behavior

One of the most interesting findings is a "switch" that happens based on temperature.

• Cold & Slow: At lower temperatures, the "stickiness" is dominated by one type of particle interaction (non-leptonic). It's like driving on a highway where only one lane is open.
• Hot & Fast: As the temperature rises, a second type of interaction (semi-leptonic) kicks in and becomes the dominant factor. It's like a second lane opening up, completely changing the traffic flow.

The authors found a specific "crossing temperature" where this switch happens. If you are simulating a star merger, you need to know which "lane" is open to get the physics right.

Why Does This Matter?

When two neutron stars merge, they create a chaotic, vibrating mess.

• If the matter is too sticky, the vibrations die out quickly, and the gravitational wave signal is short and quiet.
• If the matter is less sticky, the vibrations last longer, creating a longer, louder signal.

By using this new Burgers equation, scientists can build better computer simulations. Instead of guessing how the star behaves, they can use this precise "two-spring" model. This helps astronomers interpret the signals from detectors like LIGO and Virgo, potentially revealing if the cores of these stars are made of normal matter or this exotic quark soup.

In a Nutshell

The authors took a complex problem (how exotic matter resists being squeezed) and realized the old, simple math wasn't enough. They found that this matter acts like a two-speed shock absorber. By writing down the new math (the Burgers equation) and calculating the specific settings for different star densities, they gave astronomers a better tool to decode the sounds of the universe.

Technical Summary

1. Problem Statement

The paper addresses the challenge of modeling bulk viscous dissipation in unpaired quark matter (matter composed of deconfined quarks without color superconductivity) under astrophysical conditions, specifically during compact star mergers and stellar oscillations.

• Context: In dense matter, compression and rarefaction drive the system out of chemical equilibrium. The return to equilibrium via weak interactions (electroweak processes) generates bulk viscosity, which damps oscillations and influences gravitational wave signals.
• Limitation of Current Models: Standard first-order hydrodynamics (Navier-Stokes) suffers from acausality and instability. Second-order hydrodynamics (Israel-Stewart or IS) introduces a relaxation time but assumes a single relaxation mode. Previous studies often calculated bulk viscosity as a frequency-dependent coefficient, which is difficult to implement in time-domain numerical simulations of relativistic hydrodynamics.
• Gap: While Gavassino [36] showed that hadronic matter with specific conditions follows a Burgers-type equation, it was unclear if unpaired quark matter (which has different reaction topologies and chemical potentials) could be described similarly.

2. Methodology

The authors derive a closed evolution equation for the bulk viscous pressure ($\Pi$) starting from first principles of kinetic theory and hydrodynamics.

• System Setup:
  • Particles: Massless $u, d$ quarks, massive $s$ quarks, massless electrons ($e$), and neutrinos ($\nu$).
  • Processes: The system is driven out of equilibrium by five electroweak processes: one nonleptonic ($u+d \leftrightarrow u+s$) and four semileptonic (involving neutrinos).
  • Chemical Potentials: The deviation from equilibrium is characterized by two independent chemical potential differences: $\mu_1 = \mu_s - \mu_d$ and $\mu_2 = \mu_u - \mu_d + \mu_e$.
• Derivation Steps:
  1. Linearization: The evolution equations for particle fractions are linearized around chemical equilibrium.
  2. Matrix Formulation: The system of differential equations for the chemical potentials is expressed in matrix form involving the inverse susceptibility matrix ($\chi^{-1}$) and the reaction rate matrix ($\lambda$).
  3. Reduction: By taking the second convective derivative and utilizing the properties of the $2\times2$ transport matrix $T$, the authors derive a second-order differential equation for the bulk scalar $\Pi$.
  4. Resulting Equation: The derived equation is identified as a Burgers-type equation:
     $$\tau_+ \tau_- D_u^2 \Pi + (\tau_+ + \tau_-) D_u \Pi + \Pi = -\zeta \vartheta - \xi D_u \vartheta$$
     Where $D_u$ is the convective derivative, $\vartheta = \partial_\mu u^\mu$ is the expansion scalar, and $\tau_\pm, \zeta, \xi$ are transport coefficients.
• Equations of State (EOS): The transport coefficients are explicitly evaluated for two distinct regimes:
  1. Modified MIT Bag Model: Valid for intermediate densities ($\mu_d \approx 400\text{--}600$ MeV).
  2. Perturbative QCD (pQCD): Valid for high densities ($\mu_d \gtrsim 850$ MeV).
• Reaction Rates: The authors use leading-order electroweak rates ($\lambda_1, \lambda_2, \lambda_3$) for nonleptonic and semileptonic processes, accounting for finite temperature and strange quark mass effects.

3. Key Contributions

• Burgers Equation for Quark Matter: The primary contribution is proving that the bulk viscous pressure in unpaired quark matter is governed by a two-component Burgers equation, rather than the standard single-relaxation-time Israel-Stewart equation.
• Identification of Transport Coefficients: The authors define four key transport coefficients expressed in terms of equilibrium parameters and reaction rates:
  • Two relaxation times: $\tau_+$ and $\tau_-$.
  • Two bulk viscosity components: $\zeta_+$ and $\zeta_-$ (derived from $\zeta$ and $\xi$).
• Green's Function Analysis: By analyzing the Green's function of the Burgers equation, the authors identify two distinct relaxation modes. The system behaves as a single-component fluid (IS limit) only when one relaxation time is effectively infinite (e.g., at low temperatures where semileptonic processes are suppressed).
• Regime Mapping: The paper delineates the temperature and density regions where nonleptonic processes dominate dissipation versus regions where semileptonic processes become crucial.

4. Key Results

• Relaxation Times ($\tau_\pm$):
  • $\tau_+$ scales linearly with temperature ($T$) in log-log plots for both EOSs.
  • $\tau_-$ exhibits complex behavior: it is piecewise linear with an inflection point. In the pQCD regime, $\tau_-$ saturates to a density-independent value at high $T$, whereas in the Bag model, it remains strongly density-dependent and can become comparable to $\tau_+$.
• Bulk Viscosity Components ($\zeta_\pm$):
  • $\zeta_-$ (associated with the slower mode) saturates to a constant at low $T$.
  • $\zeta_+$ (associated with the faster mode) overtakes $\zeta_-$ at approximately $T \approx 0.1$ MeV.
• Crossing Temperature ($T_{cross}$):
  • The authors define a crossing temperature where the contributions of the two Green's function modes ($G_+$ and $G_-$) are equal.
  • Below $T_{cross}$: Dissipation is dominated by nonleptonic processes ($G_+$), and the system can be approximated by a single-component IS fluid.
  • Above $T_{cross}$: Both channels contribute significantly, requiring the full Burgers equation description.
  • For merger-relevant timescales ($10^{-4}$ to $10^{-2}$ s), this crossing occurs at temperatures relevant to neutron star mergers (typically $< 1$ MeV for high densities).
• Parameter Sensitivity:
  • The bulk viscosity is highly sensitive to the strange quark mass ($m_s$) in the Bag model, scaling roughly as $m_s^5$ for $\zeta_-$ at low temperatures.
  • Uncertainties in pQCD (renormalization scale) significantly affect $\zeta_-$ but less so other coefficients.

5. Significance

• Numerical Simulations: The derived Burgers equation provides a causal, stable, and implementable framework for incorporating bulk viscous effects into numerical relativity simulations of neutron star mergers. Unlike frequency-dependent formulations, this time-domain equation can be directly coupled to hydrodynamic solvers.
• Astrophysical Observables: By accurately modeling the damping of stellar oscillations and post-merger dynamics, this work helps refine predictions for gravitational wave signals. The transition between nonleptonic and semileptonic dominance could leave distinct signatures in the damping rates of oscillation modes.
• Theoretical Unification: The work unifies the description of bulk viscosity in quark matter, showing that the Israel-Stewart limit is merely a special case (when one relaxation time diverges) of a more general Burgers-type behavior inherent to multi-reaction systems.
• Future Extensions: The framework is designed to be adaptable, allowing for future inclusion of more accurate EOSs, updated reaction rates, and potentially superconducting phases of quark matter.

In summary, this paper establishes that unpaired quark matter requires a two-relaxation-time description (Burgers equation) to accurately capture bulk viscous dissipation, offering a critical tool for interpreting future gravitational wave data from compact star mergers.