Jet-environment interaction after delayed collapse in binary neutron star mergers

This study presents general relativistic magnetohydrodynamic simulations demonstrating that the lifetime of a metastable massive neutron star remnant critically shapes the polar outflow environment, thereby significantly influencing the propagation and final properties of the subsequent black hole-driven jet in binary neutron star mergers.

The Gist

Imagine two neutron stars—city-sized spheres of matter so dense that a teaspoon of them would weigh a billion tons—spiraling toward each other and colliding. This cosmic crash is one of the most violent events in the universe. When they smash together, they don't just make a mess; they can create a "short gamma-ray burst" (GRB), which is essentially a laser beam of pure energy shooting out into space, powerful enough to be seen across the entire universe.

But how does this happen? And why do some crashes create a powerful laser while others just fizzle out?

This paper is like a high-speed, supercomputer movie that simulates exactly what happens in the seconds after these stars collide. The researchers, led by Jay Kalinani, wanted to solve a mystery: How does the "engine" of the explosion change depending on how long the merged star survives before turning into a black hole?

Here is the story of their discovery, broken down into simple concepts:

1. The Two-Act Play

When the two stars merge, they don't immediately become a black hole. Instead, they form a giant, spinning, super-hot "monster star" (called a Massive Neutron Star or MNS). This monster star has two possible fates:

• The Quick Collapse: It survives for a short time (about 25 milliseconds) and then implodes into a black hole.
• The Long Haul: It survives longer (about 50 milliseconds) before collapsing.
• The Survivor: It never collapses at all (in one of their test cases).

The researchers found that how long the monster star survives changes everything about the explosion that follows.

2. The "Pre-Game" Warm-up (The Polar Outflow)

Before the black hole even forms, the spinning monster star acts like a cosmic sprinkler. It shoots out massive streams of hot gas and magnetic fields from its North and South poles.

• Analogy: Imagine a spinning ice skater throwing off water droplets from their hands. The faster they spin, the further the water flies.
• The Problem: If the monster star lives longer, it throws off more water. By the time the black hole forms, the space around it is clogged with this "pre-game" debris.

3. The "Jet" vs. The "Traffic Jam"

Once the monster star collapses into a black hole, a new engine kicks in. The black hole, spinning with a disk of material around it, tries to launch a super-fast jet of energy (the GRB).

• The Short Life Scenario (25ms): The monster star didn't have time to throw off much debris. The path is relatively clear. The black hole's jet blasts through easily, accelerating to near the speed of light. It's like a Formula 1 car on an empty track.
• The Long Life Scenario (50ms): The monster star lived longer and filled the track with thick, heavy gas. When the black hole tries to launch its jet, it hits a wall of "traffic." The jet has to fight its way through this dense fog. It slows down, gets clogged with heavy particles (baryons), and loses energy. It's like that same Formula 1 car trying to drive through a swamp.

4. The "Shock" and the "Flash"

When the jet crashes into the heavy gas left behind by the monster star, it creates a massive shockwave.

• Analogy: Think of a supersonic jet breaking the sound barrier, creating a sonic boom. Here, the jet is breaking through a wall of gas, creating a "shock heating" effect.
• The Result: This collision might create a "precursor" signal—a flash of light that happens before the main gamma-ray burst. It's like seeing the flash of a gun before hearing the bang. This could explain strange signals astronomers have seen before the main explosion.

5. The Technical Secret Sauce

The researchers had to be incredibly clever with their computer code to get this right.

• The "Floor" Problem: In computer simulations, you can't have "zero" density (empty space) because the math breaks. Usually, scientists set a "floor"—a minimum amount of fake gas everywhere to keep the math working.
• The Innovation: Previous simulations used a floor that was too "thick" (too much fake gas). This would have accidentally slowed down the jet in the computer, making the results wrong.
• The Fix: These researchers used a floor that gets thinner and thinner the further you go from the center (like a fog that clears up as you drive away). This allowed them to simulate the jet traveling thousands of kilometers without the "fake gas" getting in the way. It was like clearing the air so perfectly they could see the jet all the way to the edge of the simulation.

The Big Takeaway

The main lesson of this paper is that timing is everything.

• If the merged star collapses quickly, the path is clear, and we get a bright, fast, successful gamma-ray burst.
• If the merged star lingers, it clogs the path with its own debris. The jet might get choked off, or it might struggle so hard that it produces a different kind of signal (like a precursor flash) or a dimmer burst.

This helps astronomers understand why some neutron star collisions create brilliant cosmic fireworks, while others might be "choked" and fail to produce a visible burst. It turns out, the life story of the merged star dictates the finale of the show.

Technical Summary

1. Problem Statement

Binary neutron star (BNS) mergers are the leading progenitors of short gamma-ray bursts (GRBs). While General Relativistic Magnetohydrodynamic (GRMHD) simulations have successfully demonstrated that black hole (BH)–disk systems can launch jets with terminal Lorentz factors ($\Gamma_\infty \gtrsim 100$) compatible with GRBs, significant uncertainties remain regarding the formation conditions and propagation environment of these jets.

Specifically, realistic BNS remnants often form a metastable massive neutron star (MNS) that survives for tens of milliseconds before collapsing into a BH. During this delay, the MNS can launch massive polar outflows. The paper addresses the lack of self-consistent simulations investigating how a subsequent BH-driven jet interacts with these pre-collapse MNS outflows. Key questions include:

• How does the MNS lifetime (delayed collapse) affect the density and structure of the environment the jet must traverse?
• Can a jet successfully break out of this dense, baryon-loaded cocoon to power a GRB?
• What are the electromagnetic (EM) signatures of the interaction between the jet and the pre-collapse wind?

2. Methodology

The authors performed high-resolution GRMHD simulations using the Spritz code, evolving the Einstein equations (BSSN formulation) and fluid equations.

• Physical Setup:
  • System: Equal-mass BNS ($1.3624 M_\odot$ each, chirp mass $1.186 M_\odot$, matching GW170817).
  • Equation of State (EOS): A hybrid piecewise-polytropic EOS (APR4) with a thermal component. Crucially, the adiabatic index $\gamma_{th}$ was set to $1.8$ for $\rho > 10^{10}$ g/cm$^3$ and $4/3$ for $\rho < 10^{10}$ g/cm$^3$ to mimic realistic microphysical behavior at low densities.
  • Magnetic Fields: Dipolar fields initialized with $B_{pole} \approx 10^{15}$ G (unrealistically strong to compensate for resolution limits, but treated as a perturbation).
• Simulation Cases:
  • Case A: Collapse induced at $\approx 25$ ms post-merger.
  • Case B: No collapse (MNS survives indefinitely).
  • Case C: Same as A but lower resolution (used for long-term evolution up to 160 ms).
  • Case D: Collapse induced at $\approx 50$ ms post-merger.
• Technical Innovation (Density Floor):
  • A major technical hurdle in jet simulations is the "numerical density floor" ($\rho_{atm}$), which can artificially impede jet acceleration.
  • The authors employed an unprecedentedly low floor scaling as $\rho_{atm} \propto r^{-6}$ (down to $\sim 6.17 \times 10^4$ g/cm$^3$ at 74 km). This ensures the floor does not affect outflows even at distances of $\sim 10^4$ km, unlike previous works using shallower profiles (e.g., $r^{-3}$).

3. Key Results

A. Jet Launching and Magnetization

• Pre-collapse Phase: All cases produced a helical magnetic tower extending $>500$ km, driving polar outflows with velocities up to $\sim 0.2c$.
• Post-collapse Phase: Upon BH formation, a new, highly collimated, magnetically dominated jet emerged.
  • Magnetization ($b^2/2\rho$): Near the BH, magnetization reached $>10^3$ (Cases A/C) and a few hundred (Case D), implying a potential terminal Lorentz factor $\Gamma_\infty \gtrsim 100$.
  • Poynting Flux: The jet luminosity ($L_{EM}$) stabilized at $\sim 3 \times 10^{50}$ erg/s (Cases A/C) and $\sim 10^{51}$ erg/s (Case D), consistent with the Blandford-Znajek mechanism.

B. Impact of Collapse Time (MNS Lifetime)

The delay between merger and collapse fundamentally altered the jet's propagation environment:

• Early Collapse (25 ms, Cases A/C): The MNS-driven outflow was less massive and extended. The BH-driven jet encountered less resistance, accelerating rapidly to $v_r \approx 0.77c$ and reaching $\sim 3000$ km by 75 ms.
• Delayed Collapse (50 ms, Case D): The longer-lived MNS launched a significantly more massive and extended polar outflow. This created a "heavier cocoon" of baryon-loaded material.
  • The jet struggled to penetrate this environment, retaining higher densities and lower terminal velocities.
  • The jet front stalled repeatedly, showing an intermittent advancement pattern rather than continuous acceleration.
  • Conclusion: The collapse time is a critical parameter determining whether a jet can successfully breakout or remains choked.

C. Jet-Environment Interaction

• Shock Heating: As the BH-driven jet pushed through the dense MNS outflow, it generated strong shocks.
  • These shocks caused localized heating (increased internal energy density $\rho\epsilon$) and sharp jumps in entropy ($s_{eff}$).
  • The interaction was confined to small opening angles ($\lesssim 15^\circ$) around the jet axis.
• Electromagnetic Signatures: The shock heating and energy dissipation during the "stalling" phase could produce distinct EM signals, potentially observable as GRB precursors or contributing to the early emission of the burst.

D. Non-Collapsing Case (Case B)

• In the absence of a BH, the MNS drives a steady, collimated outflow.
• However, this outflow remains much denser and slower ($v_r \lesssim 0.3c$) compared to BH-driven jets, suggesting that powering a typical short GRB without a central BH engine is unlikely under these conditions.

4. Significance and Conclusions

1. Jet Diversity in BNS Mergers: The study provides a coherent framework for explaining the diversity of observed BNS outcomes (from failed jets/kilonovae to successful GRBs). The MNS lifetime is identified as the key variable: shorter lifetimes favor successful jet breakout, while longer lifetimes create dense environments that may choke the jet or delay breakout.
2. GRB Precursors: The self-consistent capture of the jet interacting with the pre-collapse wind suggests that shock heating during the breakout phase could generate observable precursor signals before the main GRB emission.
3. Technical Advancement: The use of the $r^{-6}$ density floor is a critical methodological breakthrough. It allows for the reliable simulation of jet propagation to large scales ($\sim 10^4$ km) without artificial suppression, a prerequisite for connecting numerical models to observed GRB phenomenology.
4. Future Implications: The results imply that multi-wavelength diagnostics (e.g., photospheric emission delays) could potentially be used to infer the collapse time of the MNS remnant in future GW/EM events.

In summary, the paper demonstrates that while BH-driven jets are capable of reaching GRB-compatible Lorentz factors at their base, their ability to escape and power an observable burst is heavily dependent on the density of the environment left behind by the metastable neutron star, which is directly controlled by the collapse timescale.