Helium as an Indicator of the Neutron-Star Merger Remnant Lifetime and its Potential for Equation of State Constraints

By analyzing helium abundance constraints in the kilonova AT2017gfo, this study infers that the binary neutron star merger remnant collapsed into a black hole within 20–30 ms, thereby ruling out numerous nuclear equation of state models and establishing tight upper limits on neutron star radii and maximum mass.

The Gist

Imagine two neutron stars as a pair of incredibly dense, city-sized marbles made of pure atomic nuclei, spinning around each other in the dark. When they finally crash into each other, it's the most violent event in the universe since the Big Bang. But what happens immediately after the crash? Do they instantly smash together into a black hole, or do they bounce around for a while as a super-hot, spinning "super-star" before collapsing?

This paper, written by a team of astrophysicists, proposes a clever new way to answer that question by looking at the "smoke" left behind: the helium.

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

1. The Cosmic "Smoke" (The Kilonova)

When these stars collide, they fling out a cloud of debris. This debris glows brightly for a few days, creating what we call a kilonova (like a supernova, but smaller and made of heavy elements). Think of this debris cloud as a giant, expanding fog.

Usually, scientists look at the light from this fog to figure out what's inside. But this team realized that the composition of the fog tells a secret story about how long the "super-star" survived after the crash.

2. The Helium "Timer"

The authors argue that if the super-star survives for even a little bit longer (say, a few seconds), it acts like a giant, glowing blowtorch. It shoots out a powerful wind of particles (neutrinos) that hits the debris cloud.

This wind has a specific job: it turns the heavy elements in the cloud into helium.

• Short Life: If the super-star collapses into a black hole almost instantly (within a blink of an eye, or ~20-30 milliseconds), it doesn't have time to blow much helium into the cloud.
• Long Life: If it survives for a second or more, it pumps out a massive amount of helium, turning the cloud into a helium-rich soup.

3. The Detective Work (Looking at the Light)

The team looked at the light from the famous crash event GW170817 (which happened in 2017). They specifically looked at a specific "fingerprint" in the light spectrum—a dark line that appears when helium is present.

They found something surprising: There was almost no helium.

The light showed that the helium content was extremely low (less than 5%). If the super-star had survived for even a second, the helium would have been much higher, and that dark line would have been huge. The fact that the line was missing (or very faint) meant the "blowtorch" was turned off almost immediately.

The Verdict: The super-star didn't just survive for a second; it collapsed into a black hole in about 20 to 30 milliseconds. That is faster than a human eye can blink, and even faster than a computer can process a single frame of a video.

4. What This Tells Us About the "Rules of the Universe"

Why does this matter? Because the speed at which the star collapses depends on the Equation of State (EoS).

Think of the Equation of State as the "recipe" for neutron star matter. It tells us how "squishy" or "stiff" the matter is.

• If the matter is very stiff (like a hard rock), the star can hold up more weight before collapsing.
• If the matter is soft (like a marshmallow), it collapses under its own weight much easier.

By proving the star collapsed so quickly, the authors put a strict limit on how "stiff" the recipe can be. They calculated that:

• Neutron stars cannot be too big (their radius is likely around 11–12 km, not 14 km).
• They cannot be too heavy (the maximum weight a neutron star can hold is likely around 2.3 times the mass of our Sun).

This rules out many existing theories that suggested neutron stars could be much larger or heavier. It's like finding out that a specific type of building material can only support a 3-story building, not a 10-story one, which forces architects to redesign their blueprints.

5. The "Engine" of the Explosion

The crash also produced a short gamma-ray burst (a flash of high-energy light). Scientists have debated what powered this flash:

• Scenario A: A super-strong magnetic star (a magnetar) spinning wildly.
• Scenario B: A black hole with a swirling disk of gas (a torus).

Because the star collapsed so fast (in 20ms), it didn't have time to become a magnetar. Therefore, the "engine" that fired the gamma-ray burst must have been the black hole and its swirling gas disk. It's like realizing a car didn't have time to warm up its engine before the crash, so the explosion must have come from the fuel tank, not the engine.

Summary

In short, this paper is a cosmic detective story. By noticing that the "smoke" (debris) from a stellar crash didn't smell like helium, the scientists deduced that the "fire" (the super-star) went out almost instantly. This tiny detail allows them to rewrite the rules of how matter behaves at the highest densities in the universe, telling us that neutron stars are smaller, lighter, and more fragile than many scientists previously thought.

It's a powerful reminder that even the smallest chemical clues (like a tiny bit of helium) can unlock the biggest secrets of the cosmos.

Technical Summary

1. Problem Statement

The lifetime of the hypermassive neutron star (HMNS) remnant formed after a binary neutron star (NS) merger is a critical but poorly constrained parameter. This lifetime ($\tau_{BH}$) determines the threshold binary mass ($M_{thres}$) for prompt black hole (BH) formation. Knowing $M_{thres}$ allows physicists to constrain the nuclear Equation of State (EoS) of dense matter, specifically the maximum mass ($M_{max}$) and radius ($R$) of non-rotating neutron stars.

Previous constraints relied on:

• Kilonova brightness: Arguing that a bright kilonova (like AT2017gfo) implies a delayed collapse (ruling out prompt collapse), setting a lower limit on $M_{thres}$.
• Gamma-ray burst (GRB) timing: Arguing that the short delay of GRB170817A implies a short lifetime if powered by a BH-torus, but this is degenerate with the "magnetar" scenario.
• Nucleosynthesis patterns: Previous studies suggested strontium overproduction in long-lived remnants, but uncertainties in strontium modeling and its production conditions made this inconclusive.

The Gap: There was no robust, direct observational method to place an upper limit on the remnant lifetime or the threshold mass, leaving the EoS constraints incomplete and allowing for models with simultaneously large radii and large maximum masses.

2. Methodology

The authors propose a novel method linking the helium abundance in the merger ejecta to the remnant lifetime, using the spectral features of the kilonova AT2017gfo.

A. Observational Constraint (Spectroscopy)

• Target: The kilonova AT2017gfo (from GW170817) at 4.4 days post-merger.
• Feature: The spectral region around 800–1200 nm, specifically the He I $\lambda$1083.3 nm line.
• Modeling:
  • Used Non-Local Thermodynamic Equilibrium (NLTE) radiative transfer modeling to simulate the helium level populations (specifically the $1s2s \ ^3S$ state).
  • Calculated P Cygni profiles to compare synthetic spectra with the observed VLT/X-shooter spectrum.
  • Result: The absence of a strong absorption feature in the 800–1000 nm range implies that the helium mass fraction in the polar ejecta (velocity range $0.19c \lesssim v \lesssim 0.3c$) must be low: $X_{He} \lesssim 0.05$.

B. Theoretical Constraint (Hydrodynamics & Nucleosynthesis)

• Simulations: Utilized a suite of "end-to-end" neutrino-hydrodynamic simulations (from Just et al. [133]) covering dynamical ejecta, the HMNS phase, and the post-collapse BH-torus phase.
• Mechanism:
  • Neutrino-driven winds from a long-lived HMNS are rich in neutrinos, leading to high electron fractions ($Y_e \approx 0.5$) and high entropies.
  • These conditions are optimal for helium production via the r-process (specifically the $\alpha$-process).
  • In contrast, dynamical ejecta and BH-torus ejecta typically have lower $Y_e$ and produce negligible helium.
• Correlation: The simulations show that the helium mass fraction ($X_{He}$) in the polar ejecta increases steeply with the HMNS lifetime ($\tau_{BH}$). Long-lived remnants ($\tau_{BH} > 30$ ms) produce $X_{He} \gg 0.05$, while short-lived remnants produce very little helium.

C. Linking to EoS

• Lifetime to Mass: The remnant lifetime is a steep function of the total binary mass ($M_{tot}$). A short lifetime implies $M_{tot}$ is very close to the threshold mass for prompt collapse ($M_{thres}$).
• EoS Constraints: Using empirical fit formulae (Bauswein et al. [40]) that relate $M_{thres}$ to $M_{max}$ and $R$, the authors derived upper limits on these stellar parameters.

3. Key Contributions

1. First Spectral EoS Constraint: This is the first study to use specific kilonova spectral features (helium lines) to constrain the nuclear EoS, moving beyond lightcurve modeling or gravitational wave data alone.
2. Upper Limit on Remnant Lifetime: Establishes that the HMNS in GW170817 collapsed within $\tau_{BH} \lesssim 20\text{--}30$ ms.
3. Upper Limit on Threshold Mass: Infers that the threshold mass for prompt collapse is $M_{thres} \lesssim 2.93 M_\odot$ (assuming $\Delta M \approx 0.2 M_\odot$ above the GW170817 total mass).
4. Combined EoS Exclusion: Demonstrates that the combination of a short lifetime (upper limit on $R$ and $M_{max}$) and the "no prompt collapse" argument (lower limit on $R$) creates a narrow "allowed window" for NS parameters, ruling out a significant number of current EoS models.

4. Key Results

• Helium Abundance Limit: The helium mass fraction in the polar ejecta of AT2017gfo is constrained to $X_{He} \lesssim 0.05$ (90% confidence).
• Remnant Lifetime: To satisfy the helium limit, the HMNS must have collapsed within $\tau_{BH} \lesssim 20\text{--}30$ ms.
• Threshold Mass: The threshold mass for prompt collapse is estimated as $M_{thres} \lesssim 2.93 M_\odot$.
• Neutron Star Parameters:
  • Maximum Mass ($M_{max}$): The upper limit is constrained to $M_{max} \lesssim 2.3 M_\odot$ (depending on the assumed $\Delta M$ and causality arguments).
  • Radius ($R_{1.6}$): For a standard $M_{max} = 2.0 M_\odot$, the radius of a $1.6 M_\odot$ NS is constrained to $12 \pm 1$ km. For $M_{max} = 2.15 M_\odot$, the range tightens to $11.5 \pm 1$ km.
  • Exclusion: Models with simultaneously large radii ($R > 13$ km) and large maximum masses ($M_{max} > 2.3 M_\odot$) are ruled out.
• Central Engine: The short lifetime strongly favors a Black Hole-Torus central engine for GRB170817A, ruling out the long-lived magnetar scenario.
• GW190814 Interpretation: The derived upper limit on $M_{max}$ ($\lesssim 2.3 M_\odot$) suggests the secondary object in GW190814 (2.5–2.67 $M_\odot$) is likely a black hole, unless it is rapidly spinning.

5. Significance and Implications

• Tightening EoS Constraints: By combining an upper limit (from helium) and a lower limit (from kilonova brightness), the paper significantly narrows the allowed parameter space for the high-density EoS, effectively ruling out "stiff" EoS models that predict both large radii and high maximum masses.
• Central Engine of sGRBs: Provides strong evidence that short gamma-ray bursts are powered by black hole accretion disks (Blandford-Znajek mechanism) rather than magnetars, at least for events similar to GW170817.
• Future Observations:
  • Prompt Collapse Candidates: Future events with $M_{tot} \gtrsim M_{GW170817} + \Delta M$ should show prompt collapse (dim kilonovae, no post-merger GW signal, and no helium features).
  • Long-Lived Candidates: Events with lower total mass should exhibit strong helium spectral features due to longer-lived HMNSs.
  • Orphan Kilonovae: The presence or absence of helium features in kilonovae without GW counterparts could indicate the total mass relative to GW170817.

6. Limitations and Future Work

The authors acknowledge that this is an exploratory study. Key uncertainties include:

• Radiative Transfer: The NLTE modeling of helium and the contribution of other elements (like Strontium) to the 1 $\mu$m feature.
• Hydrodynamics: The dependence of helium production on the specific EoS, mass ratio, and the treatment of neutrino transport and magnetic fields in simulations.
• Binary Mass Ratio: The constraints are sensitive to the mass ratio ($q$) of GW170817; a highly asymmetric system would weaken the constraints.

Future work aims to refine the NLTE models, improve hydrodynamic simulations with full MHD and Boltzmann neutrino transport, and apply this method to future multi-messenger events.