Fact-checked by Grok 2 weeks ago

Neutron star merger

A neutron star merger is the violent collision of two neutron stars orbiting each other in a binary system, which gradually spiral inward over millions of years due to energy loss via gravitational radiation.^[1] These stars, the ultra-dense remnants of massive stellar explosions with diameters of about 10-20 kilometers and masses typically 1-2 times that of the Sun, merge in milliseconds, unleashing gravitational waves detectable across vast cosmic distances.^[2]^[3] Such events are multimessenger phenomena, producing not only gravitational waves but also electromagnetic signals including short gamma-ray bursts—triggered within seconds of the merger—and a subsequent kilonova from the radioactive decay of ejected neutron-rich material, which synthesizes heavy elements like gold and platinum through rapid neutron capture.^[3]^[2] The merger's outcome depends on the total mass: if exceeding roughly 2.5-3 solar masses, it forms a black hole, as evidenced by the low post-merger emissions in observed cases; otherwise, a hypermassive neutron star may briefly persist before collapsing.^[4] The landmark detection of GW170817 on August 17, 2017, by LIGO, Virgo, and numerous telescopes, originating from a system 130 million light-years away in galaxy NGC 4993, confirmed neutron stars as sources of short gamma-ray bursts and advanced multi-messenger astronomy.^[3]^[2] These mergers illuminate extreme regimes of general relativity, nuclear physics, and cosmic element formation, with ongoing detections by gravitational-wave observatories revealing their frequency and diversity.^[1]

Fundamentals

Neutron stars

Neutron stars are the collapsed remnants of massive stars that have undergone supernova explosions, primarily composed of neutrons with densities approaching that of atomic nuclei. These objects form when the core of a progenitor star, initially more massive than about 8 solar masses, collapses under gravity after exhausting its nuclear fuel. The resulting neutron star typically has a mass between 1.4 and 2 solar masses but is compressed into a sphere with a radius of roughly 10–15 km, making it one of the densest forms of matter in the universe excluding black holes.^[5]^[6]^[7] Key properties of neutron stars include their extreme density, on the order of nuclear saturation density at approximately $10^{17} kg/m³, where matter behaves as a degenerate Fermi gas of neutrons. Many neutron stars possess strong magnetic fields, reaching up to $10^{15} Gauss in the case of magnetars, which are a subclass powered by magnetic energy decay. Additionally, rapidly rotating neutron stars known as pulsars can spin at rates up to 716 Hz, emitting beamed radiation that appears as pulses when viewed from Earth. These properties arise from the intense gravitational compression and conservation of angular momentum during formation.^[8]^[9]^[10] The internal structure and stability of neutron stars depend on the equation of state (EOS) of ultra-dense matter, which remains uncertain due to the extreme conditions unattainable in laboratories. This EOS governs possible compositions, ranging from conventional neutron-rich nuclear matter to exotic states like hyperon matter or strange quark matter. A fundamental limit on neutron star masses is provided by the Tolman-Oppenheimer-Volkoff (TOV) equation, which describes hydrostatic equilibrium in general relativity and predicts a maximum stable mass of approximately 2–3 solar masses, beyond which collapse to a black hole occurs. The TOV equation is given by

\frac{dP}{dr} = -\frac{G M(r) \rho(r)}{r^2} \left(1 + \frac{P}{\rho c^2}\right) \left(1 + \frac{4\pi r^3 P}{M c^2}\right) \left(1 - \frac{2 G M}{r c^2}\right)^{-1},

where P is pressure, \rho is energy density, M(r) is the enclosed mass, G is the gravitational constant, and c is the speed of light.^[6]^[11]

Binary systems

Neutron star binaries form primarily through two channels: the evolution of isolated massive binary star systems and dynamical interactions in dense stellar environments. In the isolated channel, a pair of massive stars (typically 8–20 solar masses each) evolves such that both undergo core-collapse supernovae, producing neutron stars; this often involves a common envelope phase where the expanding envelope of the first supernova progenitor engulfs the companion, leading to orbital tightening via drag and envelope ejection. Dynamical formation predominates in globular clusters, where neutron stars capture companions through three-body encounters or exchange interactions, resulting in tighter orbits prone to further evolution.^[12] After formation, double neutron star binaries undergo evolutionary stages driven by angular momentum loss, which causes orbital decay over timescales of 10^7 to 10^9 years. In systems where one neutron star has been recycled via accretion—often observed as low-mass X-ray binaries earlier in evolution—mass transfer from a low-mass companion can spin up the pulsar, but post-supernova double neutron star systems primarily shrink due to gravitational wave emission, with supplementary losses from magnetic braking in the pulsar wind. The iconic binary pulsar PSR B1913+16 exemplifies this, with its observed orbital period decrease matching general relativity predictions for gravitational wave-driven inspiral.^[13] Stability criteria for these binaries hinge on their total mass, analogous to a Chandrasekhar-like limit; systems with combined masses exceeding approximately 2.9 solar masses reach a threshold during late inspiral where the merged remnant transitions from a stable hypermassive neutron star to prompt collapse into a black hole, depending on the equation of state. Population statistics derived from binary pulsar surveys, including PSR B1913+16 and others, yield an estimated merger rate for neutron star binaries of roughly 10–1000 Gpc⁻³ yr⁻¹, calibrated against observed double neutron star systems and informed by population synthesis models.^[14]

Merger dynamics

Inspiral phase

The inspiral phase of a neutron star merger is dominated by the gradual orbital decay of the binary system due to energy and angular momentum loss through gravitational wave emission. As the two neutron stars orbit each other, the quadrupole moment of the system varies, producing gravitational waves that carry away energy, causing the orbital separation to decrease and the orbital frequency to increase progressively. This process begins at low frequencies, typically in the millihertz range for wide binaries detectable by space-based observatories, and accelerates to kilohertz frequencies in the final stages observable by ground-based detectors like LIGO and Virgo. The overall duration of this phase spans approximately $10^4 to $10^8 years for the transition from millihertz to merger, depending on the binary's initial separation and masses, though the final detectable seconds in the sensitive band of current instruments last only tens to hundreds of seconds. The rate of energy loss is described by the Peters-Mathews formula, derived in the quadrupole approximation for circular orbits:

\left\langle \frac{dE}{dt} \right\rangle = -\frac{32}{5} \frac{G^4}{c^5} \frac{\mu^2 M^3}{a^5},

where G is the gravitational constant, c is the speed of light, \mu = M_1 M_2 / (M_1 + M_2) is the reduced mass, M = M_1 + M_2 is the total mass, and a is the orbital separation. This can be equivalently expressed in terms of the chirp mass M_\mathrm{chirp} = (M_1 M_2)^{3/5} / (M_1 + M_2)^{1/5} and orbital frequency f_\mathrm{orb}:

\left\langle \frac{dE}{dt} \right\rangle \propto \frac{(G M_\mathrm{chirp})^{5/3} (\pi f_\mathrm{orb})^{10/3}}{c^5},

highlighting how the energy loss rate scales with the binary's intrinsic properties and accelerates as the frequency rises. For typical neutron star binaries with masses around 1.4 M_\odot each, this leads to a coalescence timescale \tau = \frac{5}{256} \frac{c^5}{(G M_\mathrm{chirp})^{5/3} (\pi f_\mathrm{gw})^{8/3}} (where f_\mathrm{gw} = 2 f_\mathrm{orb}), consistent with observed systems like PSR J0737-3039 having \tau \sim 85 Myr from current separation.^[15] As the separation shrinks to the Roche limit, typically around 10–20 km for neutron stars with radii of about 10–12 km, tidal forces become significant, deforming the stars and potentially leading to mass transfer or partial disruption before full merger. The Roche lobe radius for each star scales roughly as R_\mathrm{L} \approx 0.49 q^{2/3} a / (0.6 q^{2/3} + \ln(1 + q^{1/3})) (where q = M_2 / M_1), and when R_\mathrm{L} \approx R_\mathrm{NS}, tidal bulging induces deviations from point-mass behavior, imprinting finite-size effects on the gravitational waveform through tidal deformability parameters. These effects cause an earlier phase advance in the waveform compared to point-particle predictions, providing probes of the neutron star equation of state. Modeling the inspiral requires a transition from post-Newtonian (PN) approximations, valid for early weak-field stages up to 3.5PN order, to full numerical relativity (NR) simulations near merger where strong-field general relativity dominates. PN expansions describe the orbital evolution analytically, incorporating corrections for tidal interactions via effective-one-body (EOB) frameworks that match to NR waveforms for hybrid models. Full GR NR simulations, using formulations like BSSN (Baumgarte-Shapiro-Shibata-Nakamura), resolve the late inspiral over tens of orbits, capturing nonlinear tidal deformations and the approach to contact with high accuracy, as demonstrated in early works evolving irrotational binaries to merger. These simulations confirm that tidal disruption begins just prior to the dynamical merger, with gravitational wave frequencies reaching ~1–2 kHz.

Merger and ringdown

The merger phase of a binary neutron star system occurs rapidly once the stars come into contact, typically over a timescale of approximately 0.1–1 ms, during which tidal deformation causes significant distortion of the stellar structures.^[16] This violent collision releases a substantial fraction of the system's binding energy in the form of gravitational waves, with total radiated energy on the order of $10^{53} erg, corresponding to roughly 3–5% of the total rest mass energy for typical systems with component masses around 1.4 M_\odot.^[16] Numerical relativity simulations indicate that the peak gravitational wave luminosity during this phase can reach up to $10^{56} erg/s, dominating the signal as the stars coalesce into a compact remnant.^[16] The outcome of the merger is either a hypermassive neutron star (HMNS) supported temporarily by differential rotation or a prompt black hole formation, depending primarily on the total baryonic mass of the system.^[16] Following the merger, the ringdown phase ensues as the remnant settles into a quasi-equilibrium state, characterized by oscillations that emit gravitational waves in the form of quasi-normal modes (QNMs). These modes damp exponentially over a timescale of 10–100 ms, with dominant frequencies in the range of 1–4 kHz for HMNS remnants, reflecting the excitation of fundamental oscillation modes.^[16] If the remnant collapses to a black hole, the ringdown is described by perturbations in the Kerr metric, governed by the Teukolsky equation, which models the damped sinusoidal waves from the spinning black hole's horizon. During this phase, additional gravitational wave energy of approximately $10^{50}–$10^{51} erg is radiated, providing a clean probe of the remnant's final spin and mass.^[16] A key dynamical feature of the merger is the expulsion of neutron-rich ejecta, driven by tidal torques and shock heating, with typical masses of 0.01–0.1 M_\odot and velocities ranging from 0.1c to 0.3c.^[17] This material, characterized by low electron fractions (Y_e \lesssim 0.2), undergoes rapid neutron capture (r-process) nucleosynthesis, synthesizing heavy elements beyond the iron peak, such as lanthanides and gold.^[17] The ejecta dynamics are highly sensitive to the binary mass ratio and equation of state, with more asymmetric systems producing greater amounts of fast-moving, neutron-rich outflow.^[17] The final outcome hinges on the post-merger remnant's stability: prompt collapse to a black hole occurs if the total mass exceeds approximately 2.5–3 M_\odot, leading to immediate horizon formation without a prolonged HMNS phase.^[18] For lower masses (around 2.2–2.7 M_\odot), a metastable HMNS forms, which may undergo delayed collapse after 10–100 ms due to angular momentum redistribution and neutrino cooling, or remain stable if supramassive and supported by rotation.^[16] These scenarios are delineated through numerical relativity, with the threshold mass varying by 10–20% depending on the neutron star equation of state.^[18]

Detection methods

Gravitational waves

Gravitational waves from neutron star mergers are produced as the binary system loses orbital energy through quadrupole radiation, manifesting as spacetime ripples detectable by specialized interferometers. The signal is characterized by three distinct phases: an inspiral phase featuring a chirp-like waveform where the frequency and amplitude increase as the neutron stars spiral closer, typically sweeping through the sensitive band of ground-based detectors over durations of about 100 seconds; a brief merger phase producing a high-amplitude burst lasting milliseconds; and a ringdown phase with an exponentially decaying waveform as the remnant settles into a stable configuration.^[19]^[20] These signals are primarily detected using ground-based laser interferometers such as the Laser Interferometer Gravitational-Wave Observatory (LIGO), Virgo, and Kamioka Gravitational Wave Detector (KAGRA), which operate in the frequency range of approximately 10–1000 Hz, aligning with the late inspiral, merger, and early ringdown of neutron star binaries.^[21] Future space-based detectors like the Laser Interferometer Space Antenna (LISA), sensitive to millihertz frequencies, will probe the earlier inspiral phases of such systems at lower frequencies, enabling detection of galactic and nearby extragalactic events.^[22]^[23] Parameter estimation from the detected waveform relies on Bayesian inference and matched filtering techniques, which compare the observed signal to theoretical templates to extract properties such as the chirp mass (a combination of the individual masses), component spins, and luminosity distance to the source. For instance, the event GW170817 yielded a chirp mass of approximately 1.188 solar masses, near-zero spins for both components, and a distance of about 40 megaparsecs, with a signal-to-noise ratio of around 32 achieved through optimal matched filtering that maximizes detection in noisy data.^[24] Observations up to the fourth Gravitational-Wave Transient Catalog (GWTC-4) in 2025 constrain the local merger rate density of binary neutron stars to 7.6–250 Gpc^{-3} yr^{-1}, based on hierarchical Bayesian analyses of detected and subthreshold candidates.^[21] This rate informs expectations for future observing runs, such as O5, where improved sensitivity is projected to yield several detections annually.^[25]

Electromagnetic signals

Neutron star mergers are expected to produce a variety of electromagnetic (EM) signals across multiple wavelengths, arising from the violent dynamics of the collision and its aftermath. These signals include short gamma-ray bursts (sGRBs) from relativistic jets, kilonovae powered by the radioactive decay of merger ejecta, and longer-lasting afterglows from interactions with the surrounding interstellar medium (ISM). While not strictly electromagnetic, potential neutrino bursts from the hot merger remnant complement these observations in multi-messenger studies, though they remain undetected. The detection of these EM counterparts relies on rapid follow-up to gravitational wave (GW) triggers using observatories such as Fermi for gamma rays, Hubble for optical imaging, Chandra for X-rays, and radio telescopes like the Very Large Array.^[26] Short gamma-ray bursts represent one of the most energetic EM counterparts, occurring in at least 40% of mergers due to the launch of ultra-relativistic jets powered by accretion onto a central black hole or magnetar remnant. These jets, collimated within opening angles of ~5–20 degrees, produce prompt gamma-ray emission through internal dissipation mechanisms such as shocks or magnetic reconnection in the jet outflow. The emission peaks in the gamma-ray band (1 keV to 10 MeV) with isotropic-equivalent energies around 10^{46–52} erg and durations of 0.01–2 seconds, followed by a softer tail. Beaming effects limit detectability to favorable viewing angles, with the local merger rate implying a beamed sGRB rate of ~1000 Gpc^{-3} yr^{-1}.^[26] Kilonovae provide a more isotropic EM signal, emerging from the radioactive decay of heavy r-process elements synthesized in the sub-relativistic ejecta expelled during the merger. This ejecta, with masses of ~0.01–0.1 M_\sun and velocities ~0.1–0.3c, becomes optically thick due to high opacities from lanthanides, leading to a "red" kilonova in the infrared. The luminosity peaks at ~10^{40–41} erg s^{-1} in the optical to infrared bands days after the merger, with peak times of 1–7 days and total radiated energy ~10^{47} erg, fading over weeks as the ejecta expands and cools to temperatures ~5000 K. Lanthanide-free "blue" components may appear earlier in the ultraviolet/optical if present.^[26] Afterglow emission follows the prompt signals, originating from synchrotron radiation as the relativistic jet or its cocoon interacts with the ISM, producing structured outflows observable over wider angles. This non-thermal emission spans X-ray, optical, and radio wavelengths, with X-ray fluxes peaking at ~100 μJy around 100 days post-merger and decaying as t^{-2}, optical fading rapidly but with potential infrared excesses at ~1 week, and radio rising to peaks of ~20 μJy at 100–400 days in dense environments (n \geq 0.01 cm^{-3}). The afterglow persists for weeks to years, providing constraints on jet energetics and ISM density.^[26] Neutrino emissions, while not electromagnetic, are predicted from the hot (~10 MeV) merger remnant and neutrino-driven winds, with total energies up to ~10^{51} erg in all flavors, primarily in the GeV–EeV range from beta-decays and pair processes. These bursts could last milliseconds to seconds but have evaded detection, with upper limits <10^{50} erg from joint observatory searches.^[26]^[27] In multi-messenger astronomy, EM signals are correlated with GW detections to localize mergers and probe physics inaccessible to GWs alone, such as ejecta composition and jet launching efficiency. Follow-up campaigns using telescopes like Fermi, Chandra, Hubble, and ground-based optical/radio arrays enable rapid sky surveys within minutes to hours of a GW alert, enhancing event characterization.^[26]^[28]

Observed events

GW170817

GW170817 was the first binary neutron star merger detected through gravitational waves, observed on August 17, 2017, by the Advanced LIGO and Advanced Virgo observatories.^[29] The event originated in the elliptical galaxy NGC 4993 at a luminosity distance of approximately 40 megaparsecs.^[30] Analysis of the gravitational wave signal indicated component masses of about 1.17 and 1.60 solar masses, with a chirp mass of 1.188 solar masses.^[31]^[29] The multimessenger observations of GW170817 included a gravitational wave inspiral signal followed by an electromagnetic counterpart. Approximately 1.7 seconds after the merger, the short gamma-ray burst GRB 170817A was detected by the Fermi Gamma-ray Burst Monitor and the INTEGRAL satellite, marking the first direct association between gravitational waves and a short gamma-ray burst.^[29] An optical transient, designated SSS17a, was identified 10.9 hours post-merger, evolving from a blue-dominated spectrum in the early phase—attributed to lanthanide-poor ejecta—to a redder component days later, consistent with a kilonova powered by r-process nucleosynthesis in neutron-rich material. A structured jet produced a broad afterglow detectable across the electromagnetic spectrum, from gamma rays to radio waves, with observations spanning X-ray, ultraviolet, optical, infrared, and radio bands over months.^[32] The merger remnant was inferred to be a short-lived hypermassive neutron star that promptly collapsed into a black hole, based on the absence of prolonged gravitational wave emission or persistent electromagnetic signatures indicative of a stable remnant.^[33] No evidence supported a long-lived magnetar, as the post-merger signal lacked the expected high-frequency ringdown modes.^[31] GW170817 provided transformative insights into astrophysical phenomena. It confirmed the gravitational wave origin of short gamma-ray bursts through the observed coincidence and off-axis jet structure. As a standard siren, the event enabled a direct measurement of the Hubble constant at approximately 70 km/s/Mpc, independent of the cosmic distance ladder, by combining the luminosity distance from gravitational waves with the redshift of the host galaxy.^[34]

Recent candidates

Following the landmark detection of GW170817, subsequent gravitational-wave observations have identified several candidate neutron star mergers, though most lack confirmed electromagnetic counterparts. One of the earliest post-2017 events was GW190425, detected on April 25, 2019, by the LIGO Hanford and Livingston observatories, with Virgo contributing marginally. This signal is consistent with a binary neutron star merger at a luminosity distance of approximately 159 Mpc, featuring a total mass of about 3.4 M_⊙ consistent with two neutron stars each of mass approximately 1.4 M_⊙ (90% credible intervals 1.12–2.52 M_⊙). No electromagnetic counterpart was identified despite extensive follow-up searches across optical, infrared, and radio wavelengths, likely due to the event's distance and sky localization challenges.^[35]^[36]^[37] In May 2023, the LIGO Livingston detector alone captured GW230529, a signal from the merger of two compact objects with masses of roughly 1.5–1.9 M_⊙ and 2.5–4.5 M_⊙, placing it in the mass gap between the heaviest neutron stars and lightest black holes. While initially interpreted as a possible neutron star-black hole merger, analyses suggest it could represent a binary neutron star system if the heavier object is an unusually massive neutron star, challenging models of the neutron star equation of state and the mass gap's boundaries. The event's low signal-to-noise ratio and single-detector origin limited precise sky localization, preventing electromagnetic follow-up confirmation.^[38] By 2025, ongoing observations during the LIGO-Virgo-KAGRA O4 run yielded additional candidates, including the low-significance event S250818k on August 18, 2025, reported as a potential subsolar-mass binary neutron star merger with a false alarm rate indicating sub-threshold confidence. This candidate, with inferred component masses below typical neutron star values, was associated with the transient AT2025ulz, discovered by the Zwicky Transient Facility, which exhibited early-time optical and infrared emission suggestive of a kilonova. Follow-up with the James Webb Space Telescope in the infrared band revealed spectral features consistent with r-process nucleosynthesis in neutron-rich ejecta, though the association remains tentative due to the event's marginal detection.^[39]^[40] The fourth Gravitational-Wave Transient Catalog (GWTC-4), released in August 2025, incorporated data from the initial phase of O4 (up to early 2024), adding 128 new compact binary coalescence candidates to the previous tally, none of which are confident binary neutron star mergers, highlighting the ongoing scarcity of such detections based on mass distributions and waveform morphology. These include events with total masses under 3 M_⊙ and no evidence of post-merger remnants exceeding neutron star limits, though none have confirmed multimessenger signals.^[21]^[41] Observing these recent candidates has highlighted persistent challenges, such as the scarcity of electromagnetic counterparts, often attributable to poor sky localization from single- or dual-detector triggers and the faintness of distant kilonovae. Merger rates inferred from GWTC-4 suggest binary neutron star coalescences occur at 13–170 Gpc⁻³ yr⁻¹ (90% credible interval), higher than pre-2017 estimates from binary pulsar surveys, indicating that such events may be more common but harder to localize without full network sensitivity.^[21]^[42]

Implications

Heavy element production

Neutron star mergers serve as primary sites for the rapid neutron-capture process (r-process) nucleosynthesis, where neutron-rich material ejected during the event captures neutrons at high rates to form heavy elements beyond iron.^[43] In the ejected material, characterized by a low electron fraction (Y_e < 0.25), free neutrons enable rapid successive captures, bypassing the slow beta-decay bottlenecks of lighter element formation and producing isotopes of lanthanides, gold, platinum, uranium, and other actinides.^[44] Typical yields from these events range from 0.01 to 0.05 solar masses (M_⊙) of r-process material, with dynamical ejecta contributing around 10^{-3} to 10^{-2} M_⊙ and neutrino-driven winds from the post-merger accretion disk adding up to 0.05 M_⊙.^[44] The electromagnetic counterpart to these mergers, known as kilonovae, reveals the stratification of ejecta velocities and compositions through distinct spectral components. The blue kilonova, peaking early at bluer wavelengths, arises from faster outflows (∼0.1–0.2c) with higher Y_e (∼0.3–0.4), producing lighter r-process elements like strontium with lower opacity.^[45] In contrast, the red kilonova emerges later from slower, equatorial ejecta (∼0.05–0.1c) rich in heavy r-process nuclei (Y_e ∼0.1–0.2), featuring high opacity from lanthanides that reddens and dims the emission.^[45] This velocity-dependent layering, with outflows reaching up to 0.3c, powers the kilonova luminosity via radioactive decay of the freshly synthesized isotopes.^[44] Observations of the GW170817 merger provided direct confirmation of r-process production, with total ejecta mass estimated at ∼0.04 M_⊙, including lanthanide-rich components that matched spectral models for heavy element synthesis.^[45] These findings align with simulations distinguishing dynamical ejecta—launched tidally and by shocks during the inspiral and merger—from neutrino-irradiated winds off the remnant disk, where varying Y_e influences the isotopic yields.^[43] Opacity models incorporating lanthanide absorption lines are essential for interpreting kilonova light curves, as they account for the suppressed early emission from heavy-element regions.^[44] In terms of galactic chemical evolution, neutron star mergers are the dominant source of r-process elements observed in the Milky Way.^[46] This dominance arises because merger ejecta are neutron-richer (Y_e ≲ 0.25) than other potential sources, enabling efficient production of elements heavier than barium.^[44]

Remnant formation

The outcome of a neutron star merger is determined primarily by the total mass of the binary system and the nuclear equation of state, which governs the maximum stable mass of a neutron star. For systems with a total gravitational mass exceeding approximately 2.8 M_\odot, where M_\odot denotes the solar mass, the merger leads to the prompt formation of a black hole immediately following the collision, as the combined object surpasses the threshold for gravitational collapse. This threshold, often expressed as roughly 1.41 times the maximum mass of a non-rotating neutron star, varies slightly with the equation of state stiffness but typically falls in the range of 2.7–3.0 M_\odot for realistic models.^[47]^[48] In mergers with total masses between about 2.0 and 2.8 M_\odot, a hypermassive neutron star forms as the initial remnant. This object, exceeding the maximum mass of a stable neutron star but temporarily supported by differential rotation from the merger dynamics, survives for milliseconds to seconds before collapsing into a black hole due to angular momentum redistribution and viscous dissipation. The exact lifetime depends on the remnant's rotation profile and equation of state, with softer equations of state promoting quicker collapse. Formation of a long-lived stable neutron star is rare, occurring only in low-mass binaries where the total mass is below approximately twice the Tolman–Oppenheimer–Volkoff limit (around 2.0–2.3 M_\odot for typical neutron star masses of 1.0–1.4 M_\odot each), which represents a small fraction of expected merger events based on observed binary populations.^[47]^[49]^[50] If the merger remnant is a rapidly rotating neutron star with an initial spin period of approximately 1 ms—achievable through conservation of orbital angular momentum—it can develop into a magnetar. In this scenario, magnetohydrodynamic instabilities amplify the magnetic field to strengths of about $10^{16} G within the first few milliseconds post-merger, powered by the remnant's differential rotation and winding of pre-existing seed fields. Such magnetars are potential progenitors of soft gamma repeaters, which emit recurrent bursts of soft gamma rays due to magnetic field decay and crustal instabilities over longer timescales.^[51]^[52]^[53] For prompt black hole formation, the remnant typically has a mass of 2–3 M_\odot, reflecting the binary's total mass minus a small fraction lost to gravitational waves and ejecta (about 3–5%). The dimensionless spin parameter a (where $0 \leq a \leq 1) ranges from 0.7 to 0.9, arising predominantly from the orbital angular momentum transferred during the merger, with contributions from the neutron stars' intrinsic spins being minor unless they are highly aligned. Per the no-hair theorem of general relativity, these black holes are fully described by their mass and spin alone, as electric charge is negligible in astrophysical contexts, leading to axisymmetric Kerr metrics that dictate their future evolution and surrounding spacetime.^[54]^[55]^[56] Numerical simulations have provided critical insights into remnant formation, particularly in 2025 advancements from the Max Planck Institute for Gravitational Physics. These general relativistic magnetohydrodynamic models, incorporating neutrino transport and realistic equations of state, track the post-merger evolution for up to 1.5 seconds using massive computational resources (e.g., 130 million CPU hours). For high-mass binaries (e.g., 1.25 + 1.65 M_\odot), they demonstrate prompt black hole formation accompanied by rapid magnetic field amplification in the surrounding accretion disk, which launches relativistic jets along the black hole's spin axis—key to interpreting gamma-ray burst associations. These simulations highlight how remnant properties influence multi-messenger signals, with fast-spinning black holes enhancing jet collimation and energy extraction efficiency.^[57]^[58]

Broader effects

Neutron star mergers pose no significant threat to Earth due to their immense distances and the minuscule amplitude of the gravitational waves they produce. For instance, the strain from the detected event GW170817, located approximately 40 megaparsecs away, was on the order of h \sim 10^{-21}, representing a fractional change in spacetime metric far too weak to cause any measurable physical effects on planetary bodies or life. Similarly, the associated electromagnetic radiation and potential cosmic ray contributions from such distant mergers are heavily diluted, rendering any spikes in cosmic rays negligible at Earth's location. Within galaxies like the Milky Way, neutron star mergers serve as tracers of underlying stellar populations, primarily occurring in the disks and bulges where binary systems form from massive stars. Their locations reflect the distribution of these stellar components, with simulations indicating that merger sites align with the galactic stellar mass density in these regions.^[59] Over the lifetime of a galaxy, approximately 10^5 to 10^6 such mergers are estimated to have contributed to the observed abundances of heavy elements, providing a key mechanism for galactic chemical enrichment without requiring alternative dominant sources.^[60] On cosmological scales, these mergers enable precise measurements of the Hubble constant H_0 through the standard siren method, where the gravitational wave signal directly infers luminosity distance, potentially resolving the longstanding tension between early- and late-universe estimates of cosmic expansion. As of 2025, analyses of multiple detected mergers have yielded Hubble constant estimates with uncertainties around 5-10%, with ongoing observations expected to resolve the tension between early- and late-universe measurements to better than 2% precision.^[61] Observed merger rates further constrain models of stellar evolution by testing binary formation channels and mass-loss prescriptions, while also placing limits on dark matter interactions within neutron stars, as deviations in merger signals could indicate dark matter accumulation altering stellar remnants.^[62] Although rare, the gamma-ray bursts occasionally produced by neutron star mergers could pose a sterilization risk to planetary atmospheres if occurring nearby, within about 10 kiloparsecs, by depleting ozone layers and exposing surfaces to lethal ultraviolet radiation. However, the probability of such an event in a typical galaxy remains low, estimated at roughly one per gigayear, making it an infrequent hazard over planetary timescales.^[63]

References

[1]
Sources and Types of Gravitational Waves | LIGO Lab | Caltech
This means that black-hole merger signals are much shorter in duration than neutron star merger signals, and the differences are quite striking. For example, ...
[2]
When Dead Stars Collide! - NASA Science
Oct 16, 2017 · Neutron stars are the crushed leftover cores of massive stars (stars more than 8 times the mass of our sun) that long ago exploded as supernovae ...Missing: definition | Show results with:definition
[3]
GW170817 Press Release | LIGO Lab | Caltech
At the moment of collision, the bulk of the two neutron stars merged into one ultradense object, emitting a “fireball” of gamma rays.
[4]
Gravitational Wave Event Likely Signaled Creation of a Black Hole
May 31, 2018 · These two neutron stars were so close that gravitational wave radiation pulled them together until they merged and collapsed into a black hole.Missing: definition | Show results with:definition
[5]
Introduction to neutron stars - UMD Astronomy
What is its structure like? First, the typical mass of a neutron star is about 1.4 solar masses, and the radius is probably about 10 km. By the way, the " ...
[6]
[PDF] The physics of neutron stars - arXiv
Jun 22, 2011 · Neutron stars, the most compact stars in the Universe, were given this name because their interior is largely composed of neutrons.
[7]
Neutron Stars and White Dwarfs | Center for Astrophysics | Harvard ...
In the range between 8 and 20 times the Sun's mass, the remaining collapsed core gives birth to a neutron star. These objects are more massive than the Sun, but ...
[8]
[PDF] Neutron Star Proper-es and the Role of the Nuclear Symmetry Energy
Jan 16, 2014 · 10 – 15 km. Magne-c Field: 1010 – 1015 Gauss. Pressure: ~1029 atm. Density: ~1017 kg/m3. Temperature: ~106 K. Period: ~(ms to sec). Surface ...
[9]
Magnetars: Special Stars With That Attractive Charm
Estimates from the March 5, 1979 event indicated that SGRs had magnetic fields ranging from 100 trillion to a quadrillion (1015) gauss. Theorists began to ...
[10]
Astronomers Discover Fastest-Spinning Pulsar
Jan 12, 2006 · The newly-discovered pulsar is spinning 716 times per second, or at 716 Hertz (Hz), readily beating the previous record of 642 Hz from a pulsar discovered in ...
[11]
Higher mass limits of neutron stars from the equation of states in ...
Sep 20, 2021 · In order to solve the Tolman-Oppenheimer-Volkoff equations for neutron stars, one routinely uses the equation of states which are computed in ...
[12]
Formation and evolution of compact binaries in globular clusters
Jun 27, 2007 · We study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing neutron stars (NSs) in dense globular ...
[13]
An increased estimate of the merger rate of double neutron stars ...
Dec 2, 2003 · A reliable estimate of the double-neutron-star merger rate in the Galaxy is crucial in order to predict whether current gravity wave detectors ...Missing: 10-1000 Gpc^ 3} yr^ 1}
[14]
Double neutron stars: merger rates revisited - Oxford Academic
Within our study of the parameter space, we find solutions that allow for DNS merger rates as high as R l o c a l ≈ 600 − 300 + 600 Gpc−3 yr−1 which are thus ...
[15]
https://wwwusers.ts.infn.it/~milotti/Didattica/GravitationalWaves/handouts/20-CBC.pdf
[16]
Dynamical mass ejection from binary neutron star mergers
Mar 26, 2015 · We find that the properties of the dynamical ejecta of the merger such as total mass, average electron fraction, and thermal energy depend ...
[17]
[1801.05680] An explosion is triggered by the late collapse of ... - arXiv
Jan 17, 2018 · Abstract:It is known that a binary neutron star merger produces a hypermassive neutron star. The lifetime of this compact remnant depends on ...
[18]
Multimessenger Binary Mergers Containing Neutron Stars ... - arXiv
Feb 5, 2021 · In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general ...
[19]
Gravitational waves from binary neutron stars | Arabian Journal of ...
Jan 16, 2022 · I review the current global status of research on gravitational waves emitted from mergers of binary neutron star systems, ...Missing: Peters- formula
[20]
GWTC-4.0: Population Properties of Merging Compact Binaries - arXiv
Aug 25, 2025 · ... mergers: binary neutron star, neutron star--black hole binary, and binary black hole mergers. We resolve multiple over- and under-densities ...
[21]
Identifying galactic binary systems of neutron stars and black holes ...
Jul 15, 2025 · The Laser Interferometer Space Antenna (LISA) will detect ~ 100 galactic binary systems comprised of black holes (BHs) and neutron stars (NSs).
[22]
How can LISA probe a population of GW190425-like binary neutron ...
If GW190425-like binaries exist in the Milky Way, LISA will detect them within the volume of our Galaxy and will measure their chirp masses to better than 10 ...
[23]
[PDF] Gravitational Wave Analysis : Matched Filtering
This document briefly presents the way to compute the signal-to-noise ratio (SNR) through the Optimal Matched Filtering and to find gravitational waves ...
[24]
Detecting Electromagnetic Counterparts to LIGO/Virgo/KAGRA ...
Oct 22, 2025 · In this work, we make end-to-end simulations of GW events during the current O4 run and the upcoming O5 run (expected to start in 2028) using ...
[25]
https://iopscience.iop.org/article/10.3847/1538-4357/ae0336
[26]
Linking the rates of neutron star binaries and short gamma-ray bursts
Jan 21, 2022 · We use current estimates for the BNS and NSBH merger rates to calculate the fraction of observable short gamma-ray bursts produced through each channel.
[27]
Search for High-energy Neutrinos from Binary Neutron Star Merger ...
We searched for high-energy neutrinos from the merger in the GeV-EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories.
[28]
Multimessenger Binary Mergers Containing Neutron Stars - Frontiers
GW170817: observation of gravitational waves from a binary neutron star inspiral. ... Binary neutron-star mergers: a review of Einstein's richest laboratory.
[29]
Observation of Gravitational Waves from a Binary Neutron Star Inspiral
Oct 16, 2017 · Binary neutron star mergers may result in a short- or long-lived neutron star remnant that could emit gravitational waves following the merger
[30]
A Precise Distance to the Host Galaxy of the Binary Neutron Star ...
This is the most precise distance to NGC 4993 available to date. Combining our distance measurement with the corrected recession velocity of NGC 4993 implies a ...
[31]
Properties of the Binary Neutron Star Merger GW170817
Jan 2, 2019 · The gravitational-wave signal provides information on the masses and spins of the binary's components, the distance and location of the source, ...
[32]
The Panchromatic Afterglow of GW170817: The Full Uniform Data ...
Nov 26, 2021 · The afterglow is optically thin throughout its evolution, consistent with a single spectral index (−0.584 ± 0.002) across all epochs. This gives ...
[33]
When Did the Remnant of GW170817 Collapse to a Black Hole?
Jan 14, 2019 · We determine that the remnant of GW170817 must have collapsed to a black hole after t_{\rm coll}=0.98_{-0.26}^{+0.31}\,{\rm s}.
[34]
A gravitational-wave standard siren measurement of the Hubble ...
Oct 16, 2017 · We determine the Hubble constant to be about 70 kilometres per second per megaparsec. This value is consistent with existing measurements ...
[35]
LIGO-Virgo Network Catches Another Neutron Star Collision
Jan 6, 2020 · When two neutron stars spiral together, they undergo a violent merger that sends gravitational shudders through the fabric of space and time.
[36]
Is GW190425 Consistent with Being a Neutron Star-Black Hole ...
GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total ...
[37]
LIGO-Virgo-KAGRA (LVK) Collaboration Detected a Remarkable ...
Apr 5, 2024 · Astronomers receiving the alert were informed that a neutron star and a black hole most likely merged about 650 million light-years from Earth.
[38]
https://www.ligo.caltech.edu/news/ligo20240405
[39]
https://arxiv.org/abs/2510.24620
[40]
GWTC-4.0: Updated Gravitational-Wave Catalog Released | LIGO Lab
Aug 26, 2025 · August 26th 2025 marks an important date for the international LIGO, Virgo, and KAGRA (LVK) gravitational-wave network. The collaborations are ...Missing: constraints | Show results with:constraints
[41]
GWTC-4.0: Population Properties of Merging Compact Binaries
Aug 26, 2025 · We infer the local merger rates (i.e., at redshift z=0 z = 0 ) to be 7.6 7.6 --250Gpc−3yr−1 250 G p c − 3 y r − 1 for binary neutron stars, 9.1 ...Missing: constraints | Show results with:constraints
[42]
https://dcc-llo.ligo.org/LIGO-P2400004/public
[43]
https://arxiv.org/pdf/1801.01141.pdf
[44]
https://arxiv.org/pdf/2109.09162.pdf
[45]
Prompt Black Hole Formation in Binary Neutron Star Mergers - arXiv
Feb 16, 2024 · The paper analyzes prompt collapse in binary neutron star mergers, finding a critical mass that separates prompt and non-prompt collapse, ...
[46]
Systematics of prompt black-hole formation in neutron star mergers
Oct 9, 2020 · Abstract:This study addresses the collapse behavior of neutron star (NS) mergers expressed through the binary threshold mass M_thr for ...
[47]
Multimessenger Binary Mergers Containing Neutron Stars
In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general ...
[48]
General relativistic simulations of high-mass binary neutron star ...
The possible outcomes of a BNS merger can be either a prompt collapse to a black hole (BH), or the formation of a short-lived hypermassive neutron star ...
[49]
High-energy Neutrinos from Millisecond Magnetars Formed from the ...
The merger of a neutron star (NS) binary may result in the formation of a long-lived, or indefinitely stable, millisecond magnetar remnant surrounded by a low- ...
[50]
Tayler-Spruit dynamo in binary neutron star merger remnants
The dynamo's growth timescale depends on the initial large-scale magnetic field right after the merger. In the case of a long-lived hypermassive neutron star, ...
[51]
Soft Gamma Repeaters - an overview | ScienceDirect Topics
neutron stars with an anomalously strong magnetic field of ∼ 1014 − −1015 Gs at the ejection stage (E).
[52]
Black-hole formation in binary neutron star mergers - ResearchGate
The dimensionless spin of the remnant BH formed as a result of BNS mergers are typically 0.7-0.8 [106][107] [108] unless the NS spins are extremely large [109].<|control11|><|separator|>
[53]
Coalescence of black hole–neutron star binaries | Living Reviews in ...
Dec 6, 2021 · We review the current status of general relativistic studies for coalescences of black hole–neutron star binaries.
[54]
[1208.5869] The Black Hole Remnant of Black Hole-Neutron Star ...
Aug 29, 2012 · We present a model for determining the dimensionless spin parameter and mass of the black hole remnant of black hole-neutron star mergers with ...
[55]
Breakthrough in simulating how neutron stars collide
May 30, 2025 · ... binary neutron star merger, lasting 1.5 seconds of real time. The ... During this inspiral phase, they fall towards each other as they ...
[56]
Supercomputer simulation reveals how merging neutron stars form ...
Jun 2, 2025 · Our new simulation follows the binary all the way through its evolution: inspiral, merger, and the post-merger phase, including the jet formation.
[57]
Finding the Remnants of the Milky Way's Last Neutron Star Mergers
Jul 19, 2019 · We assume that the “birth” places of the binary NS systems trace the stellar mass distribution of the MW bulge and disk. We model the latter ...
[58]
The role of neutron star mergers in the chemical evolution of the ...
Observations of heavy element abundances in Galactic halo stars provide important constraints on the astrophysical site(s) of r-process nucleosynthesis.
[59]
MERGER RATES OF DOUBLE NEUTRON STARS AND STELLAR ...
We investigate the impact on the predicted merger rates of neutron stars and black holes. Despite the changes with previous assumptions, we only find an ...
[60]
The Resilience of Life to Astrophysical Events | Scientific Reports
Jul 14, 2017 · We consider the relative likelihood of complete global sterilisation events from three astrophysical sources – supernovae, gamma-ray bursts, large asteroid ...