⚛️ general relativity

Numerical simulations of black hole-neutron star mergers with equal and near-equal mass ratios

This paper presents numerical simulations of black hole-neutron star mergers with near-equal mass ratios to address gaps in parameter space, revealing limitations in current gravitational waveform models while confirming the accuracy of remnant mass predictions and demonstrating that these systems produce detectable kilonovae.

Original authors: Ivan Markin, Mattia Bulla, Tim Dietrich

Published 2026-01-28

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Ivan Markin, Mattia Bulla, Tim Dietrich

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant dance floor where massive objects spin around each other until they crash together. Usually, when a Black Hole (the heavy, invisible dancer) meets a Neutron Star (the dense, heavy dancer), they are very different sizes. The Black Hole is usually much bigger, like a sumo wrestler dancing with a toddler.

But recently, scientists detected a crash (called GW230529) where the two dancers were much closer in size—more like a sumo wrestler dancing with a heavyweight boxer. This paper asks: What happens when these two "heavyweights" crash?

Here is a simple breakdown of what the researchers did and found:

1. The Simulation: A Cosmic Crash Test

Since we can't actually go to space and watch these crashes in real-time, the scientists built a super-accurate computer simulation. They created 12 different scenarios where a Black Hole and a Neutron Star had equal or nearly equal mass.

Think of this like a video game developer testing a new physics engine. They wanted to see if the current "rules" (mathematical models) scientists use to predict these crashes were actually correct for these specific, equal-sized partners.

2. The Soundtrack: The "Music" Was Off

When these objects crash, they send out ripples in space-time called gravitational waves. It's like the sound of the crash.

The Finding: The scientists compared their new, high-definition simulation "soundtrack" against the existing "rules" (the models).
The Result: The existing rules were wrong. The models predicted the crash would happen slightly earlier or later than it actually did in the simulation. It's like trying to predict the beat of a song, but your prediction is off by a whole drumbeat.
Why it matters: If we use the wrong rules, we might misidentify what kind of objects are crashing in the real universe.

3. The Debris: The "Splatter" and the "Ring"

When the Neutron Star hits the Black Hole, it doesn't just vanish. It gets shredded.

The Ejecta (The Splatter): Some material gets thrown out into space. The scientists found that the existing math formulas for how much stuff gets thrown out were actually pretty good.
The Disk (The Ring): Most of the Neutron Star gets swallowed but forms a swirling ring of hot, glowing gas around the Black Hole, like water going down a drain.
- Equal Mass (q=1): If the dancers are the same size, the ring forms almost instantly and becomes a perfect circle.
- Unequal Mass (q=1/2): If one is slightly smaller, the ring is messy at first, with spiral waves crashing into each other, before it finally settles down.

4. The Aftermath: The "Heartbeat" of the Ring

The scientists looked closely at how this ring of gas behaved.

The Pulse: They found that the ring doesn't just sit there; it "breathes." It has global oscillations (vibrations) that act like a heartbeat.
The Effect: These vibrations actually control how fast the gas falls into the Black Hole. It's like a faucet that opens and closes rhythmically because the water in the pipe is sloshing back and forth.
The Connection: This rhythmic "sloshing" might create a specific signal in the light (gamma rays) we see from these crashes, similar to a heartbeat monitor.

5. The Light Show: Will We See It?

When the Neutron Star is shredded, it creates a "kilonova"—a bright flash of light caused by the radioactive material flying out.

The Prediction: The scientists modeled how bright this flash would be.
The Result: If these crashes happen within about 200 million light-years of us, the flash would be bright enough for our biggest telescopes (like the Vera C. Rubin Observatory) to see it within a few days.
The Difference: The brightness depends on how "stiff" the Neutron Star is. A "stiffer" star creates a bigger, brighter explosion. A "softer" one creates a dimmer one.

Summary

This paper is essentially a "quality control" check for our understanding of the universe.

The Good News: We can predict how much debris is thrown out and how heavy the final Black Hole will be.
The Bad News: Our current models for the "sound" of the crash (gravitational waves) are inaccurate for these equal-sized partners. We need to update our math.
The New Discovery: These crashes create a unique, rhythmic "heartbeat" in the gas ring, and they produce light flashes bright enough for us to spot with our next generation of telescopes.

The authors conclude that to truly understand these cosmic collisions, we need better "maps" (waveform models) and more simulations to fill in the gaps where the old rules don't work.

Technical Summary: Numerical Simulations of Black Hole-Neutron Star Mergers with Equal and Near-Equal Mass Ratios

Problem and Motivation
The detection of the gravitational wave event GW230529, interpreted as a black hole-neutron star (BHNS) merger with a mass ratio $q \approx 1/2.6$ , suggests the existence of more symmetric BHNS systems than previously anticipated. While standard stellar evolution models predict mass ratios $q \lesssim 1/3.5$ , dynamical formation channels in dense clusters or second-generation mergers may produce systems with $q \sim 1/2$ or even $q \approx 1$ . These symmetric systems are theoretically expected to produce massive accretion disks, potentially driving bright electromagnetic transients like kilonovae. However, the parameter space for equal and near-equal mass ratios is sparsely populated by numerical-relativity (NR) simulations. This lack of data limits the accuracy of gravitational waveform models and analytical fitting formulas for remnant properties, potentially leading to misclassification of real events and missed opportunities to identify electromagnetic counterparts.

Methodology
The authors perform twelve NR simulations of non-spinning BHNS mergers using the BAM code. The simulations cover mass ratios $q \in \{1, 1/2, 1/3\}$ and neutron star (NS) masses $m_{NS} \in \{0.8, 1.2, 1.4, 2.2\} M_\odot$ . Two equations of state (EOS) are employed: the softer SLy and the stiffer DD2. The simulations utilize piecewise-polytropic representations for the EOS and evolve the systems for approximately nine orbits using the conformal Z4 formulation for spacetime and the Valencia formulation for hydrodynamics. Initial data are constructed with the FUKA code, with eccentricity reduced to below $10^{-3}$ .

The study involves:

Gravitational Waveform Analysis: Extracting waveforms at multiple radii, performing second-order extrapolation, and correcting for center-of-mass drift. These are compared against the IMRPhenomNSBH and SEOBNRv4_ROM_NRTidalv2_NSBH waveform models.
Remnant Property Verification: Calculating the masses of dynamical ejecta and accretion disks, as well as the mass and spin of the remnant black hole (BH) at 10 ms post-merger. These are compared against analytical fitting formulas from Refs. [28, 31, 40, 41].
Disk Dynamics: Analyzing the morphology, rotation laws, and global density oscillations of the remnant disks for configurations with significant disk mass.
Accretion Modeling: Categorizing accretion behaviors and identifying modulation by disk oscillations.
Kilonova Modeling: Using the POSSIS radiative transfer code to model kilonova lightcurves based on dynamical ejecta profiles and an added spherical wind component (30% of disk mass).

Key Results

Gravitational Waveforms: The current BHNS waveform models (IMRPhenomNSBH and SEOBNRv4) show poor agreement with the NR waveforms in the symmetric mass ratio regime. Most cases exhibit dephasing of approximately 1 radian at the time of merger. The models often predict earlier merger times and extended ringdown stages compared to NR data. Notably, for the subsolar-mass NS case ( $m_{NS}=0.8 M_\odot$ ), the tidal deformability exceeds the valid parameter space of existing models; when extrapolated, the models show significant phase and amplitude discrepancies.
Remnant Properties:
- Ejecta and Disk Mass: Analytical fitting formulas for dynamical ejecta mass and disk mass show good agreement with the simulation data.
- Remnant BH Mass: The analytical formulas predict the remnant BH mass with excellent accuracy (errors generally $< 1\%$ ).
- Remnant BH Spin: The updated fitting formulas are less accurate for predicting the remnant spin, particularly for systems with large tidal deformability (DD2 EOS). For the subsolar NS case, the spin prediction deviates by nearly 50%. The authors suggest this may be due to overfitting on systems with lower tidal deformabilities and more asymmetric mass ratios.
Disk Structure and Dynamics:
- Morphology: Equal-mass ( $q=1$ ) disks circularize almost immediately ( $\sim 3$ ms post-merger), whereas $q=1/2$ disks remain perturbed by spiral density waves and fallback matter for longer, settling into axisymmetry only after $\sim 8$ ms.
- Rotation Law: The specific angular momentum follows a power law $j \propto \Omega^{-n}$ . The exponent $n$ varies radially, transitioning from constant angular momentum ( $n=0$ ) inside the ISCO to sub-Keplerian values in the outer disk. This contradicts common assumptions of fixed rotation laws in long-term axisymmetric simulations.
- Oscillations: Global trapped $g$ -mode density oscillations are identified in disks with mass $> 10^{-2} M_\odot$ . Their frequencies align with the maximum prograde epicyclic frequency of the remnant BH.
- Accretion: Three accretion types are distinguished. Type B ( $q=1/2$ ) shows accretion rate variability directly modulated by the global $g$ -mode oscillations. Type A ( $q=1/2$ ) and Type C (direct plunge) show different behaviors. The $g$ -mode frequencies for some systems are compatible with high-frequency quasi-periodic oscillations (QPOs) observed in gamma-ray bursts.
Kilonova Emission: The modeled kilonovae are primarily dominated by the disk wind ejecta. Systems with the stiffer DD2 EOS produce brighter kilonovae than those with the softer SLy EOS due to more massive disks. For a distance of 200 Mpc, most systems would be detectable by the Vera C. Rubin Observatory within four days and by DECam within two days in the $r$ and $i$ bands.

Significance and Claims
The paper claims to fill a critical gap in the NR simulation parameter space for symmetric BHNS mergers. The primary significance lies in demonstrating that current gravitational waveform models are insufficiently accurate for equal and near-equal mass ratios, potentially leading to biases in parameter estimation for future detections like GW230529. The study validates analytical fitting formulas for remnant masses but highlights their limitations in predicting remnant spins for high-tidal-deformability systems. Furthermore, the detailed analysis of disk dynamics reveals complex rotation laws and oscillation-driven accretion modulation, suggesting that $g$ -modes could leave observable imprints on gamma-ray burst lightcurves. Finally, the kilonova modeling indicates that symmetric BHNS mergers are promising targets for electromagnetic follow-up with current wide-field surveys, offering a higher probability of detection compared to the asymmetric systems that have dominated previous observations.