Impact of muons on the bulk viscosity of neutron star matter metamodels

This paper investigates how the inclusion of muons in neutron star matter metamodels significantly alters the bulk viscosity, revealing that increasing the symmetry energy slope $L$ can change the viscosity by orders of magnitude and induce a double-peak structure in its frequency dependence, with potential implications for neutron star merger dynamics.

The Gist

The Big Picture: The Cosmic "Sticky" Problem

Imagine the universe as a giant kitchen. Inside this kitchen, you have the most extreme objects imaginable: Neutron Stars. These are the leftover cores of massive stars that have exploded. They are so dense that a teaspoon of their material would weigh a billion tons on Earth.

When two of these stars crash into each other (a "merger"), it creates a massive explosion and sends ripples through space-time called Gravitational Waves. Scientists use these waves to try to understand what's inside the stars.

However, there's a problem. When these stars crash, the matter inside them doesn't just flow smoothly like water. It acts like a very strange, thick fluid. This "thickness" is called viscosity. Just like honey resists being stirred, this star-matter resists being squished and stretched. This resistance is called Bulk Viscosity.

The authors of this paper are trying to figure out exactly how "sticky" this cosmic honey is, and they discovered that a tiny, invisible particle called the muon changes the recipe completely.

---

The Ingredients: The "Metamodel" Kitchen

To understand the star, you need a recipe. In physics, this recipe is called the Equation of State (EoS). It tells you how the pressure, density, and temperature of the star relate to each other.

For a long time, scientists used specific, rigid recipes. But this paper uses a "Metamodel."

• The Analogy: Imagine you are baking a cake. Instead of following one specific recipe (e.g., "2 cups of flour, 1 egg"), you use a flexible formula that says, "The cake depends on how much sugar (Parameter A) and how much butter (Parameter B) you use."
• The Goal: The authors vary these "ingredients" (specifically a parameter called the Symmetry Energy Slope, or $L$) to see how the final cake (the star's behavior) changes. They want to know: If we tweak the recipe slightly, does the cake become a sponge or a brick?

The New Ingredient: The Muon

Most people know about electrons (the tiny particles orbiting atoms). But in the super-dense core of a neutron star, there is enough pressure to create muons.

• The Analogy: Think of a crowded elevator. Usually, it's full of people (protons and electrons). But if you push the elevator down hard enough (high density), a new type of person (the muon) suddenly squeezes in.
• The Discovery: Previous studies looked at the elevator with just the original people. This paper asks: What happens to the crowd's movement when the muons show up?

The Main Finding: The "Double-Boom" Effect

The authors found that adding muons changes the "stickiness" (viscosity) of the star in two surprising ways:

1. The "Volume Knob" Effect

They found that the "stickiness" of the star is incredibly sensitive to the recipe parameter $L$.

• The Analogy: Imagine a volume knob on a stereo. Usually, turning it a little bit makes the music a little louder. But in this star, turning the knob a tiny bit makes the music go from a whisper to a deafening roar (or vice versa).
• The Result: Changing the symmetry slope ($L$) by a small amount can change the viscosity by thousands of times. This means the star could be incredibly fluid or incredibly thick depending on the exact physics we don't fully understand yet.

2. The "Double-Hump" Resonance

This is the most exciting part. When they calculated how the star reacts to being squished (oscillations), they found a new pattern.

• The Analogy: Imagine pushing a child on a swing. Usually, there is one perfect rhythm where the swing goes highest (one peak).
• The Discovery: With muons present, the swing has two perfect rhythms. If you push at the right speed, the swing goes high. If you push at a different speed, it goes high again.
• Why it matters: This "double peak" happens because the muons and electrons react at different speeds. It creates a "sweet spot" where the star absorbs energy (dissipates) much more efficiently. This only happens in specific density windows that exist inside real neutron stars.

Why Should We Care? (The Merger Connection)

Why does any of this matter to us?

1. Gravitational Waves: When two neutron stars merge, they vibrate like a bell. The "stickiness" (viscosity) of the material inside damps these vibrations.
2. The Signal: If the star is very sticky, the vibrations die out fast. If it's slippery, they last longer. This changes the "sound" of the gravitational wave detected by Earth.
3. The Muon Impact: Because muons create this "double peak" and massive changes in stickiness, they could leave a distinct fingerprint on the gravitational waves we detect.

The Bottom Line

This paper is like a chef realizing that a secret ingredient (the muon) completely changes how a cake bakes.

• Before: We thought the star's "stickiness" was predictable based on a few main ingredients.
• Now: We know that if you include muons, the stickiness can jump by orders of magnitude, and the star can vibrate in a complex "double-beat" rhythm.

The Takeaway: To accurately predict what happens when neutron stars crash (and to decode the gravitational waves they send us), we must include muons in our models. Ignoring them is like trying to bake a cake without knowing about sugar; the result will be wrong, and we might miss the true flavor of the universe.

Technical Summary

1. Problem Statement

Neutron star (NS) mergers are key sources of gravitational waves (GWs), and the dynamics of these events are heavily influenced by the transport properties of dense nuclear matter, particularly bulk viscosity. While recent studies have established that bulk viscosity is highly sensitive to the nuclear symmetry energy (specifically its slope $L$), previous calculations often neglected the presence of muons or relied on specific microscopic Equation of State (EoS) models.

The core problem addressed is the lack of a unified understanding of how muons, which appear at densities near nuclear saturation, affect the bulk viscosity in the neutrino-transparent regime. Specifically, the authors investigate:

• How the presence of muons alters the chemical equilibration rates compared to a neutron-proton-electron ($npe$) system.
• The quantitative and qualitative impact of the symmetry energy slope ($L$) on the frequency-dependent bulk viscosity when muons are included.
• Whether muons introduce new resonant structures in the viscosity profile that could significantly alter the damping of density oscillations during NS mergers.

2. Methodology

The authors employ a metamodeling approach to describe the Equation of State (EoS) of neutron star matter, allowing for a systematic exploration of nuclear parameters rather than relying on a single microscopic model.

• Equation of State (EoS): The EoS is constructed via a Taylor expansion of the energy density around nuclear saturation density ($n_0 \approx 0.15 \text{ fm}^{-3}$). It includes contributions from nucleons (neutrons and protons) and leptons (electrons and muons) treated as ideal Fermi gases. The EoS is parametrized by nuclear parameters: binding energy ($B_{sat}$), incompressibility ($K$), symmetry energy ($J$), and the symmetry energy slope ($L$).
• Chemical Equilibrium & Imbalance: The system is analyzed in the neutrino-transparent regime (low temperature, $T \lesssim 10$ MeV). The authors define two independent chemical imbalances ($\mu_1, \mu_2$) corresponding to deviations from $\beta$-equilibrium for electron and muon channels, respectively.
• Electroweak Processes: The bulk viscosity is driven by the competition between the expansion/compression of the star and the weak interaction rates that restore equilibrium. The authors consider:
  • Direct Urca (dUrca): $n \leftrightarrow p + l + \bar{\nu}_l$ (allowed only above specific density thresholds).
  • Modified Urca (mUrca): $n + N \leftrightarrow p + N + l + \bar{\nu}_l$ (always allowed but slower).
  • Muon Decay: $\mu \to e + \bar{\nu}_e + \nu_\mu$ (found to be subdominant and neglected in numerical results).
• Hydrodynamic Formalism:
  • The evolution of bulk viscous pressure is described by a Burgers equation, derived from second-order hydrodynamics (Israel-Stewart formalism).
  • This formalism yields second-order transport coefficients (relaxation times $\tau_\pm$ and viscosity components $\zeta_\pm$).
  • The frequency-dependent bulk viscosity $\zeta(\omega)$ is derived as the real part of the response function, expressed as a sum of two Lorentzian-like terms corresponding to the two relaxation modes.

3. Key Contributions

• Unified Metamodel Analysis: The paper provides a systematic analysis of bulk viscosity as a function of the symmetry energy slope $L$ (ranging from 50 to 110 MeV) using a metamodel, bridging the gap between specific microscopic models (like DDME2 or NL3) and general nuclear physics constraints.
• Inclusion of Muon Channels: Unlike previous $npe$-only studies, this work explicitly incorporates muons. It demonstrates that muons introduce a second independent relaxation channel ($\lambda_2$) distinct from the electron channel ($\lambda_1$).
• Discovery of Double-Peak Structure: The authors identify a novel double-peak resonance in the frequency-dependent bulk viscosity as a function of temperature. This structure arises specifically when the density thresholds for electron dUrca and muon dUrca processes are separated, creating two distinct relaxation times ($\tau_+$ and $\tau_-$) that are sufficiently different.
• Quantitative Impact of $L$: The study quantifies how increasing $L$ lowers the density thresholds for dUrca processes, drastically changing the viscosity by orders of magnitude.

4. Key Results

A. Density Thresholds and Reaction Rates

• The density threshold for dUrca processes decreases as the symmetry energy slope $L$ increases.
• Crucially, the threshold for electron dUrca is always lower than that for muon dUrca.
• For high values of $L$ (e.g., $L=110$ MeV), there exists a density window where electron dUrca is active, but muon dUrca is still forbidden (or just opening). In this window, the system is dominated by electron dUrca and mUrca for muons, leading to a large separation between the two relaxation times.

B. Relaxation Times and Viscosity Components

• Relaxation Times ($\tau_\pm$): When only one dUrca channel is active (e.g., electrons only), the two relaxation times become widely separated (differing by $\sim 4$ orders of magnitude). When both channels are active, they tend to converge.
• Viscosity Components ($\zeta_\pm$): The total bulk viscosity $\zeta = \zeta_+ + \zeta_-$. The components show piecewise behavior at the density thresholds of dUrca processes.
  • The emergence of electron dUrca causes a sharp drop in $\zeta_+$ but a rise in $\zeta_-$.
  • The emergence of muon dUrca causes a massive drop (up to 7 orders of magnitude) in $\zeta_-$.

C. Frequency-Dependent Bulk Viscosity $\zeta(\omega)$

• Single vs. Double Peaks:
  • In standard $npe$ matter or when muon channels are frozen, $\zeta(\omega)$ exhibits a single peak.
  • In the presence of muons, specifically in density windows where electron dUrca is active but muon dUrca is not (or just opening), $\zeta(\omega)$ exhibits a double-peak structure.
• Magnitude: The presence of muons and the specific value of $L$ can change the bulk viscosity by several orders of magnitude. For instance, increasing $L$ from 50 to 110 MeV can shift the viscosity profile significantly.
• Temperature Dependence: The double-peak structure is most prominent at specific temperatures where the inverse relaxation times match the oscillation frequency ($\omega \approx 1/\tau_\pm$).

D. Damping Times

• The damping time of density oscillations ($\tau_\zeta$) is calculated.
• Results show that for high $L$ and specific densities, damping times can be as short as milliseconds (at $\omega/2\pi = 10$ kHz), which is comparable to merger timescales.
• The lowest damping times do not necessarily correlate with the highest viscosity values but are strongly influenced by the incompressibility of the matter.

5. Significance and Implications

• Gravitational Wave Astronomy: The findings suggest that bulk viscosity could leave a detectable imprint on the gravitational wave signal from binary neutron star mergers, particularly in the post-merger phase. The double-peak structure and the sensitivity to $L$ imply that GW observations could potentially constrain the symmetry energy slope and the presence of muons in the NS core.
• Merger Simulations: Current numerical simulations of NS mergers often neglect transport coefficients. This work provides the necessary formalism and parameterized results to include bulk viscosity (specifically the muon-induced effects) in future simulations, potentially altering the predicted GW waveform and the stability of the remnant.
• Nuclear Physics Constraints: The strong dependence of viscosity on $L$ reinforces the idea that astrophysical observations (GWs) can serve as a complementary probe to terrestrial experiments for constraining the nuclear symmetry energy.

In conclusion, the paper establishes that muons are not a negligible correction but a fundamental component that qualitatively changes the transport properties of neutron star matter, introducing complex resonant behaviors in bulk viscosity that are critical for modeling the dynamics of neutron star mergers.