Consistent Treatment of Muons in Binary Neutron Star Mergers

This study presents numerical-relativity simulations of binary neutron star mergers incorporating muons and muonic reactions, revealing that their inclusion has a relatively minor impact on remnant evolution and ejecta properties compared to previous literature, suggesting less severe consequences for nucleosynthetic yields and electromagnetic counterparts.

The Gist

Imagine two neutron stars as the universe's most extreme weightlifters. They are so dense that a single teaspoon of their material would weigh a billion tons on Earth. When these two giants collide, they create a cosmic explosion so violent it ripples through space-time itself (gravitational waves) and lights up the sky with a flash of light called a "kilonova."

For decades, scientists have tried to simulate these crashes on supercomputers to understand what happens inside. But there's a problem: their simulations were missing a key ingredient. They treated the matter inside the stars as if it only contained electrons (the tiny particles that orbit atoms). They ignored muons.

What is a Muon?
Think of a muon as a "heavy electron." It's like a cousin to the electron, but it's about 200 times heavier. In the crushing pressure of a neutron star collision, the environment gets so intense that these heavy cousins are forced into existence.

This paper, titled "Consistent Treatment of Muons in Binary Neutron Star Mergers," is the first to run a full, high-definition simulation of these collisions while properly accounting for these heavy muons and the neutrinos (ghostly particles) they interact with.

Here is the breakdown of what they found, using some everyday analogies:

1. The "Heavy Cousin" Effect

In previous simulations, scientists assumed the matter inside the collision was like a crowd of light people (electrons). This paper realized that in the densest parts of the crash, the crowd actually includes heavy bodybuilders (muons).

• The Analogy: Imagine a crowded elevator. If you only count the children (electrons), you think the elevator is light. But if you realize there are also sumo wrestlers (muons) in there, the weight distribution changes.
• The Result: The presence of these "sumo wrestlers" makes the equation of state (the rules of how the matter behaves) "softer." It's like the matter becomes slightly more squishy. This allows the core of the merged star to get a tiny bit denser, but not by much.

2. The "Ghostly Neutrino" Traffic Jam

Neutrinos are like ghosts; they pass through almost everything. However, in the super-hot, super-dense center of a merger, they get stuck, like cars in a massive traffic jam.

• The Old Way: Previous simulations treated these traffic jams roughly, often ignoring the specific interactions involving muons.
• The New Way: The authors built a sophisticated "traffic control system." They tracked not just electron-neutrinos, but also muon-neutrinos. They used a new mathematical trick called the "Two Timescales Approach."
• The Analogy: Think of the matter and the neutrinos as two groups of dancers. Sometimes they dance in perfect sync (equilibrium), and sometimes they are out of step. The old models struggled to keep the rhythm when the music changed fast. The new "two timescales" method is like a conductor who knows exactly when to speed up the music and when to slow it down, ensuring the dancers (matter and radiation) eventually find their rhythm without the simulation crashing.

3. The Big Surprise: It Doesn't Change Much!

This is the most important part of the paper. Before this study, some researchers thought that adding muons would completely rewrite the story of neutron star mergers. They feared it would drastically change:

• How much stuff gets blasted out into space (ejecta).
• What elements are created (like gold and platinum).
• How bright the resulting explosion (kilonova) would look.

The Finding: The authors found that while muons are physically present and do their job, they don't change the big picture.

• The Analogy: Imagine you are baking a cake. You thought that adding a secret ingredient (muons) would turn the cake into a completely different dessert, like a soufflé. Instead, you find out that the cake still tastes 95% the same. The secret ingredient changes the texture slightly, but the cake is still a cake.
• The Numbers: The amount of material blasted out changed by less than 17%. The speed and temperature of the debris changed by less than 6%. The "flavor" of the debris (the electron fraction) remained almost identical.

4. Why This Matters

You might ask, "If it doesn't change much, why did they spend so much effort?"

1. Validation: It confirms that the simpler models scientists have been using for years are actually quite good. We don't need to panic and throw away all our previous data on neutron star collisions.
2. Accuracy: It fixes the "bugs" in the code. Even if the final result is similar, getting there the right way ensures that our predictions for future telescope observations (like the James Webb Space Telescope) are rock solid.
3. Methodology: They developed a new way to handle the "pair processes" (when particles and anti-particles annihilate each other). Think of this as inventing a new, better recipe for handling the most chaotic parts of the kitchen, ensuring the simulation doesn't "burn the cake" (crash) when things get too hot.

The Bottom Line

This paper is a "stress test" for our understanding of the universe. The scientists asked: "What happens if we finally get the physics of heavy particles right?"

The answer is reassuring: The universe is robust. Even with the inclusion of heavy muons and complex neutrino interactions, the story of how neutron stars merge, how they create heavy elements, and how they light up the sky remains largely the same. The "red" and "blue" colors of the kilonova, and the gold we find on Earth, are safe from being completely rewritten by this new discovery.

Technical Summary

1. Problem Statement

Binary Neutron Star (BNS) mergers are the primary astrophysical sites for $r$-process nucleosynthesis, producing heavy elements and powering kilonova electromagnetic counterparts. The properties of the ejected material (mass, velocity, composition, and electron fraction $Y_e$) are critical for interpreting these signals.

Current numerical-relativity simulations of BNS mergers typically rely on a 3-neutrino species (3-$\nu$) approximation, considering only electron neutrinos ($\nu_e$), electron antineutrinos ($\bar{\nu}_e$), and a single heavy-lepton species ($\nu_x$) representing muon and tau neutrinos. This approximation assumes matter contains only electrons ($e^-$) and positrons ($e^+$) as charged leptons.

However, in the dense, hot environments of BNS merger remnants, the electron chemical potential ($\mu_e$) often exceeds the muon rest mass ($m_\mu \approx 105.7$ MeV). This triggers the formation of muons ($\mu^-$) and antimuons ($\mu^+$), necessitating a 5-neutrino species (5-$\nu$) treatment that explicitly tracks muon neutrinos ($\nu_\mu, \bar{\nu}_\mu$) and the net muon fraction ($Y_\mu$). Previous studies using simplified transport schemes (e.g., leakage schemes) or approximate rates suggested that muons could drastically alter merger outcomes, potentially reducing ejecta masses by up to 50% and significantly cooling the remnant. The authors aim to resolve discrepancies in the literature by performing a consistent treatment of muons and muonic weak interactions using advanced neutrino transport and realistic reaction rates.

2. Methodology and Setups

The authors performed four numerical-relativity simulations of equal-mass BNS mergers ($M_{tot} = 2.5 M_\odot$) using the BAM code.

• Equations of State (EOS): Two baryonic baselines were used: SFHo and DD2. These were extended to include muons and antimuons for the 5-$\nu$ cases, while the 3-$\nu$ cases included only electrons/positrons.
• Neutrino Transport:
  • Implemented a gray (energy-integrated) M1 moment scheme.
  • Utilized an Implicit-Explicit (IMEX) time integrator (RK(4,4) explicit + I42L implicit) to handle the stiff coupling between matter and radiation fields.
  • Extended the transport to track 5 species: $\nu_e, \bar{\nu}_e, \nu_\mu, \bar{\nu}_\mu, \nu_x$.
• Reaction Rates and Opacities:
  • Employed full kinematics for $\beta$-reactions (e.g., $\nu_e + n \leftrightarrow p + e^-$), including mean-field potentials and weak magnetism.
  • Pair Processes: A novel approach was used for pair processes (electron-positron annihilation, nucleon-nucleon bremsstrahlung, and inverse muon decay). Instead of standard Kirchhoff's law which can lead to divergent opacities for muon neutrinos at low densities, the authors constructed finite, well-behaved gray opacities using reaction kernels and assuming Local Thermodynamic Equilibrium (LTE) for pairing neutrinos.
  • Included inverse muon decay ($\nu_\mu + e^- \to \mu^- + \nu_e$) for the first time in BNS merger simulations.
• Equilibration Scheme:
  • Introduced a "two timescales approach" to handle the equilibration of matter and radiation. This method separately tracks the equilibrium timescales for electron ($\tau_e$) and muon ($\tau_\mu$) neutrinos, allowing for a consistent treatment of partial trapping states (Free-streaming, Partially-trapped, Trapped) for both lepton families.
• Simulation Parameters:
  • Total mass: $2.5 M_\odot$ (individual $1.25 M_\odot$).
  • Grid resolution: $\Delta x \approx 185$ m in the finest level.
  • Magnetic fields were neglected to isolate neutrino physics.

3. Key Contributions

1. First Consistent 5-$\nu$ BNS Simulations: This is the first study to combine a gray M1 transport scheme, realistic full-kinematics rates, and a consistent treatment of muon decay and pair processes in a full BNS merger simulation.
2. Novel Opacity Treatment: Developed a kernel-based method to compute pair-process opacities that remain finite even when muon neutrino degeneracy is large and negative, avoiding the numerical instabilities and divergences seen in previous works.
3. Two-Timescales Equilibration: Extended the Perego et al. (2019) equilibration prescription to the 5-$\nu$ case, ensuring accurate capture of the transition between trapped and free-streaming regimes for both electron and muon neutrinos.
4. Inverse Muon Decay: Explicitly included inverse muon decay, which acts as a heating mechanism for the fluid, counteracting some of the cooling effects of muon production.

4. Key Results

The authors compared 3-$\nu$ and 5-$\nu$ simulations for both SFHo and DD2 EOSs.

• Remnant Evolution:
  • Density: The presence of muons softens the EOS, leading to slightly higher maximum densities ($\sim 1-3\%$ increase) in the immediate post-merger phase.
  • Temperature: Muonic reactions provide additional cooling, but the effect is moderated by inverse muon decay. Maximum temperatures differ by at most $\sim 2$ MeV compared to 3-$\nu$ cases.
  • Luminosity: Combined luminosities of heavy-lepton neutrinos ($\nu_\mu, \bar{\nu}_\mu, \nu_x$) in 5-$\nu$ cases converge to the 3-$\nu$ $\nu_x$ luminosity values by the end of the simulation.
• Muon Distribution:
  • Muons are primarily confined to high-density regions ($\rho \gtrsim 10^{14}$ g cm$^{-3}$) and hot annuli.
  • The "two timescales" approach successfully captures the equilibration of matter and radiation, showing that deviations from equilibrium are small ($\sim 1-3\%$) in the core but larger in low-density outflows.
• Ejecta Properties (Crucial Finding):
  • Ejecta Mass: The inclusion of muons reduces the total ejecta mass by at most $\sim 17\%$. This is significantly less severe than the $\sim 50\%$ reduction reported in previous leakage-scheme studies (Gieg et al. 2025a) and the $\sim 20-30\%$ reduction in low-$Y_e$ ejecta reported by Ng et al. (2025).
  • Composition ($Y_e$): Average electron fractions in the ejecta differ by less than $\sim 6\%$ between 3-$\nu$ and 5-$\nu$ cases. The muonic runs do not produce the extreme neutron-rich ($Y_e \lesssim 0.05$) ejecta fractions seen in Ng et al. (2025).
  • Velocity and Temperature: Asymptotic velocities and temperatures of the outflows are consistent within $\sim 5\%$.
• Discrepancy Resolution: The authors attribute the differences with Ng et al. (2025) to the treatment of neutrino degeneracies and pair-process opacities. Ng et al. rescaled degeneracies to avoid divergence, which inadvertently suppressed electronization and enhanced muonization in low-density regions. The current kernel-based approach avoids this artificial suppression.

5. Significance and Conclusion

The study concludes that neglecting muons in BNS merger simulations has limited impact on the macroscopic outcomes (remnant stability, ejecta mass, and composition) for small-mass, symmetric systems.

• Astrophysical Implications: The presence of muons and muonic weak reactions does not significantly alter nucleosynthetic yields or the resulting kilonova electromagnetic signals compared to standard 3-$\nu$ simulations. The "red" vs. "blue" kilonova components are likely not as drastically affected by muons as previously feared by some literature.
• Methodological Impact: The paper establishes a robust framework for future 5-$\nu$ simulations, demonstrating that consistent treatment of pair processes and inverse muon decay is essential to avoid unphysical results (divergent opacities or artificial cooling).
• Future Work: The authors suggest extending these methods to include magnetic fields, asymmetric mass ratios, neutrino oscillations, and more complete microphysics (e.g., pions).

In summary, while muons are physically present and dynamically relevant in the dense cores of merger remnants, their inclusion in state-of-the-art simulations with consistent transport does not lead to the dramatic suppression of ejecta or extreme compositional changes predicted by earlier, less rigorous models.