Binary neutron star mergers with a subsolar mass star

✨

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, cosmic dance floor. Usually, when two heavy dancers (neutron stars) spin toward each other to merge, they do so in a predictable rhythm. Scientists have built a "rulebook" (mathematical models) to predict exactly how they will move, what sound they will make (gravitational waves), and what happens when they finally collide.

But what if one of the dancers is much lighter than usual? What if one is a "sub-solar" neutron star—one that weighs less than our Sun?

This paper asks: If we find a pair like this, will our current rulebook still work, or will the dance look completely different?

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

1. The "Super-Sticky" Light Dancer

Neutron stars are incredibly dense. Usually, they are hard and stiff, like a rock. But the lighter a neutron star is, the "fluffier" and more stretchy it becomes.

The authors simulated a heavy star (1.7 times the Sun's mass) dancing with a very light one (0.8 times the Sun's mass). Because the light star is so stretchy, the heavy star's gravity pulls on it like taffy.

The Analogy: Imagine a heavy bowling ball spinning next to a giant, soft marshmallow. As they get close, the marshmallow doesn't just stay round; it stretches out into a long, thin tail long before they actually touch.
The Result: The light star starts "leaking" its mass onto the heavy star way earlier than scientists expected. This happens at a lower frequency (a lower pitch) than our current models predicted.

2. The "Explosion" of Debris

When two heavy stars merge, they usually fling out a small amount of debris (like dust from a car crash).

The Analogy: Think of two cars crashing. Usually, you get a few crumpled parts flying off. But in this scenario, because the light star is so stretchy and gets ripped apart early, it's like the crash happens while the marshmallow is still being pulled into a long stream.
The Result: The simulation showed that this "sub-solar" merger throws out 30 times more debris than a standard merger of two heavy stars. This is a huge amount of cosmic material, which could create very bright flashes of light (electromagnetic signals) that telescopes could see.

3. Did We Miss the Music? (The Waveform Problem)

Scientists listen for these mergers using detectors like LIGO and Virgo. They use "templates"—pre-recorded songs of what a merger should sound like—to find the real signals in the noise.

The Fear: Since the light star stretches out early, the "song" it sings is different from the standard rulebook. Scientists worried that if a real sub-solar merger happened, our detectors might think it's just background noise because the "song" didn't match the template.
The Finding: The authors checked this carefully. They found that while the end of the song (the actual crash) sounds different, the part of the song that our detectors can actually hear (the lower frequencies) is still very similar to the standard models.
The Verdict: We probably wouldn't miss these events. Our current detectors are sensitive enough to hear them, even if our "rulebook" isn't perfect for the very last second of the dance.

4. The "Black Hole" Confusion

One of the biggest mysteries in astronomy is distinguishing between a tiny black hole and a tiny neutron star. They both look the same if you only measure their weight.

The Clue: Black holes are perfectly round and rigid; they don't stretch. Neutron stars stretch.
The Finding: If we see a merger where the lighter object stretches significantly (like our marshmallow), we know for sure it's a neutron star, not a black hole. The paper shows that even with current technology, if the signal is loud enough, we can tell the difference. However, if the signal is very faint, we might get confused and think the light object is a black hole, which would lead us to calculate the wrong mass for the heavy star.

Summary: What Does This Mean for Us?

New Physics: If we find these light neutron stars, it proves that nature can make stars much lighter than we thought possible through normal star death.
Better Models: The paper tells us that our current computer models need a little update to account for these "stretchy" light stars, especially for the very loud, high-energy events we might see in the future.
No Panic: Don't worry that we've been missing these events all along. Our current detectors are good enough to catch them, and we can likely tell they are special even with our current tools.

In short, the universe might be full of these "lightweight" neutron stars, and when they dance, they splash a lot more cosmic water than we thought, but our ears are still tuned in well enough to hear the music.

1. Problem Statement

The paper addresses the potential existence of subsolar-mass neutron stars (NSs with mass $M < 1 M_\odot$ ). While theoretical models (e.g., fragmentation in collapsar disks) and tentative observational evidence suggest such objects could exist, they have not been conclusively detected. A major challenge in identifying them via gravitational waves (GWs) is distinguishing them from subsolar-mass black holes (BHs).

Current GW data analysis relies on waveform models calibrated primarily on equal-mass or near-equal-mass systems with masses $>1 M_\odot$ . Subsolar-mass NSs are predicted to have extremely large tidal deformabilities ( $\Lambda \sim O(10^4)$ to $O(10^6)$ ), significantly larger than standard NSs. The authors investigate:

How these extreme tidal effects alter merger dynamics and GW signals.
Whether current waveform models (used by LIGO/Virgo/KAGRA) can accurately detect and characterize these systems without significant bias or loss of sensitivity.
Whether current detectors can distinguish a subsolar NS from a subsolar BH.

2. Methodology

The authors performed fully general-relativistic hydrodynamic simulations of binary neutron star (BNS) mergers.

Configuration:
- Masses: A primary star of $1.8 M_\odot$ merging with a subsolar secondary of $0.7 M_\odot$ (Mass ratio $q \approx 0.39$ ).
- Equations of State (EOS): Two cold piecewise polytropic EOSs were used: BSk21 (stiff, $R_{1.4} \approx 12.6$ km) and SLy2 (softer, $R_{1.4} \approx 11.8$ km).
- Tidal Deformability: The $0.7 M_\odot$ star exhibited tidal deformabilities of $\Lambda \approx 1.5 \times 10^4$ (SLy2) and $1.88 \times 10^4$ (BSk21).
- Comparison: An equal-mass ( $q=1$ ) simulation with the same total mass ( $2.5 M_\odot$ ) was also performed for comparison.
Numerical Setup:
- Code: Einstein equations coupled to hydrodynamics using the methods described in Ref. [66].
- Resolution: Adaptive mesh refinement (7 levels), with finest grid spacing $\sim 0.033 M$ .
- Initial Data: Quasi-circular orbits with eccentricity reduced to $<10^{-3}$ .
Waveform Analysis:
- Hybridization: Numerical Relativity (NR) waveforms (covering late inspiral to post-merger) were stitched with semi-analytical inspiral waveforms (starting at 20 Hz) to create complete "hybrid" waveforms.
- Model Comparison: These hybrids were compared against state-of-the-art waveform models: IMRPhenomXAS_NRTidalv3, SEOBNRv5_ROM_NRTidalv3, SEOBNRv5THM, and the baseline BBH model SEOBNRv5_ROM.
- Parameter Estimation (PE): Injection studies were conducted using the Bilby Bayesian inference framework with LIGO/Virgo design sensitivity. Signals were injected with a network Signal-to-Noise Ratio (SNR) of $\sim 35$ .

3. Key Contributions & Results

A. Merger Dynamics and Ejecta

Mass Transfer: Due to the extreme tidal deformability of the subsolar star, significant mass transfer occurs during the late inspiral (starting $\sim 30$ ms before merger). The smaller star loses $\sim 4\%$ of its mass to the primary before the actual merger.
Disruption Frequency: The merger/disruption frequency is significantly lower than in equal-mass systems, occurring around 700–1050 Hz (within the sensitive band of current detectors).
Dynamical Ejecta: The amount of unbound matter (ejecta) is drastically enhanced compared to equal-mass systems.
- Unequal Mass ( $q=0.39$ ): $\sim 3.5 \times 10^{-2} M_\odot$ to $4 \times 10^{-2} M_\odot$ .
- Equal Mass ( $q=1$ ): $\sim 1.3 \times 10^{-3} M_\odot$ .
- Finding: The ejecta mass is $\sim 30$ times larger in the subsolar case. This contradicts existing phenomenological fitting formulas derived from super-solar mass simulations, which underestimate ejecta for high mass ratios.

B. Waveform Model Performance

Mismatch: The authors calculated the mismatch ( $\bar{F}$ $\overset{ˉ}{F}$ ) between their NR hybrids and current waveform models.
- All tidal models (IMRPhenomXAS, SEOBNRv5 variants) showed mismatches below $4 \times 10^{-3}$ across the frequency band.
- Even the BBH-only model (SEOBNRv5_ROM) showed low mismatches at lower frequencies because the signal power at merger frequencies (where tidal effects dominate) is low for current detectors.
Sensitivity Loss: The study concludes that current waveform models do not suffer significant sensitivity loss for detecting these events. The missing physics (early mass transfer) occurs at frequencies where current detectors have low sensitivity.

C. Parameter Estimation and Biases

Intrinsic Parameters: For SNRs $\lesssim 100$ (and up to $\sim 10^5$ for tidal deformability), current models recover the chirp mass, mass ratio, and effective tidal deformability ( $\tilde{\Lambda}$ ) without significant bias (within $1\sigma$ ).
Distinguishing NS vs. BH:
- If the subsolar object is modeled as a Black Hole (tidal deformability = 0), the analysis recovers the mass ratio and chirp mass with small biases, but the data slightly favors a BH-NS (BHNS) model over a BNS model (Log Bayes Factor $\approx 0.4$ ).
- Crucial Constraint: However, if one assumes the BHNS interpretation, the inferred tidal deformability of the primary (the $1.8 M_\odot$ star) becomes unphysically high ( $\Lambda \approx 990$ ) to compensate for the missing tidal signal of the secondary. This value is inconsistent with constraints from other observations (e.g., NICER, GW170817).
- Conclusion: Once realistic EOS constraints are applied, the BHNS interpretation is disfavored, and the BNS interpretation is strongly preferred (Log Bayes Factor $\approx 8.2$ ).

4. Significance

Detection Viability: The study confirms that current GW detectors (LIGO/Virgo) are capable of detecting subsolar-mass BNS mergers using existing waveform templates, despite the extreme tidal effects.
New Physics in Ejecta: The finding that unequal-mass subsolar mergers produce $\sim 30\times$ more ejecta than equal-mass systems suggests that such events could be significant sources of electromagnetic counterparts (kilonovae), potentially brighter than previously predicted for low-mass mergers.
Model Validation: The work validates that current waveform models are robust enough for initial searches but highlights the need for improved models in the late-inspiral/merger regime for future high-SNR events (e.g., in the O5 observing run) to avoid systematic biases in parameter estimation.
Identification: It provides a robust method to distinguish subsolar NSs from subsolar BHs: a large measured tidal deformability is incompatible with a BH interpretation, and the resulting tension in the primary star's inferred properties rules out the BHNS scenario.

5. Conclusion

The paper demonstrates that binary neutron star mergers involving a subsolar-mass component are dynamically distinct, characterized by early mass transfer and massive ejecta. While current waveform models do not significantly hinder detection or parameter estimation for typical SNRs, the unique signature of large tidal deformability offers a powerful tool to confirm the existence of subsolar neutron stars and rule out primordial black hole interpretations. Future observations with higher sensitivity will require refined waveform models to fully exploit the information content of these extreme events.